Engineering topics

Demarque, Pierre

Publications and source records attributed to Demarque, Pierre.

At least 19 records

DOE OSTI2017

The Sirius System and Its Astrophysical Puzzles: Hubble Space Telescope and Ground-based Astrometry

Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, the 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 M{sub ⊙} and 1.018±0.011 M{sub ⊙} 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 M{sub Jup}. The location of Sirius B in the Hertzsprung–Russell 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 on 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 subsolar metallicity, of about 0.85 Z{sub ⊙}, to fit its location on 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

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 Letter, we show that the entropy of the deeper, adiabatic layers in these simulations can be expressed as a simple function of logg and logT{sub eff}, which holds potential for calibrating stellar models in a simple and more general manner.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

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-year period by the faint DQZ white dwarf (WD) 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 M{sub ⊙} and 0.592 ± 0.006 M{sub ⊙} 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 WDs 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 M{sub ⊙}, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

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-year period by the faint DQZ white dwarf (WD) 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 plus or minus 0.012M and 0.592 plus or minus 0.006M 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 approximately 2.7 Gyr. Procyon B's location in the H-R diagram is in excellent agreement with theoretical cooling tracks for WDs 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.2M, depending on its amount of core overshoot. Several astrophysical puzzles remain. In the progenitor system, the stars at periastron were separated by only approximately 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 (approximately 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.

PROCYON↗

DOE OSTI2015

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 six 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. The 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

THE EFFECT OF METALLICITY-DEPENDENT T-τ 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, α, 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 α as well, and a recent study of stars in the Kepler field of view found α to vary with metallicity. Interpreting α obtained from calibrated stellar models, however, is complicated by the fact that the value for α depends on the surface boundary condition of the stellar model, or T-τ relation. Calibrated models that use typical T-τ relations, which are static and insensitive to chemical composition, do not include the complete effect of metallicity on α. We use three-dimensional radiation-hydrodynamic simulations to extract metallicity-dependent T-τ relations and use them in calibrated stellar models. We find the previously reported α-metallicity trend to be robust, and not significantly affected by the surface boundary condition of the stellar models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

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 using a grid of three-dimensional 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. 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 depend 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 superadiabatic layer compared to that of helium abundance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

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 three-dimensional radiation hydrodynamic simulations with four metallicities (Z = 0.040, 0.020, 0.010, 0.001), and spanning a range in effective temperature (4950 < T{sub 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 shows that the amount of superadiabaticity, which tracks the transition from efficient to inefficient convection, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

COMPARING THE EFFECT OF RADIATIVE TRANSFER SCHEMES ON CONVECTION SIMULATIONS

We examine the effect of different radiative transfer schemes on the properties of three-dimensional (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.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

The seismology of eta Bootes

Some p-mode frequencies and other observations were used to determine the mass, the age and the helium abundance of eta Bootes. It is shown how, by direct application, the p-mode frequencies and stellar seismological tools help in constraining the physical parameters of eta Boo. The existence of mode bumping is confirmed and it is discussed how it may be used to refine the estimate of the eta Boo's age. The effect of the OPAL equation of state on the p-mode frequencies is described.

Demarque, Pierre↗

NASA NTRS1995

Sensitivity of solar g-modes to varying G cosmologies

The sensitivity of the solar g-mode oscillation spectrum to variability in the universal gravitational constant G is described. Solar models in varying G cosmologies were constructed by evolving a zero-age main-sequence stellar model to the Sun's current age, while allowing the value of G to change according to the power law G(t) proportional to t(exp -beta), where Beta approximately equals delta G/GH and H is the Hubble constant. All solar models were constrained to the observed luminosity and radius at the current age of the Sun by adjusting the helium abundance and the mixing-length parameter of the models in the usual way for standard stellar models. Low-l g-mode oscillation periods were calculated for each of the models and compared to the claimed observation of the solar g-mode oscillation spectrum by Hill & Gu (1990). If one accepts Hill & Gu's claims, then within the uncertainties of the physics of the solar model calculation, our models rule out all but (delta G/GH) less than approximately 0.05. In other words, we conclude that G could not have varied by more than 2% over the past 4.5 Gyr, the lifetime of the present-day Sun. This result lends independent support to the validity of the standard solar model.

Guenther, D. B.↗

NASA NTRS1995

Modeling of shallow and inefficient convection in the outer layers of the Sun using realistic physics

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

Kim, Yong-Cheol↗

NASA NTRS1994

The gap in the color-magnitude diagram of NGC 2420: A test of convective overshoot and cluster age

Theoretical isochrones have been constructed using the OPAL opacities specifically to study the color-magnitude diagram of the open star cluster NGC 2420. This cluster provides a rare test of core convection in intermediate-mass stars. At the same time, its age is of interest because of its low metallicity and relatively high Galactic latitude for an open cluster. The excellent color-magnitude diagram constructed by Anthony-Twarog et al. (1990) allows a detailed fit of the isochrones to the photometric data. We discuss the importance of convective overshoot at the convective core edge in determining the morphology of the gap located near the main-sequence turnoff. We find that given the assumptions made in the models, a modest amount of overshoot (0.23 H(sub p)) is required for the best fit. Good agreement is achieved with all features of the turnoff gap for a cluster age of 2.4 +/- 0.2 Gyr. We note that a photometrically complete luminosity function near the main-sequence turnoff and subgiant branch would also provide an important test of the overshoot models.

Demarque, Pierre↗

NASA NTRS1994

The horizontal-branch stars in globular clusters. 2: The second parameter phenomenon

Using synthetic horizontal-branch models, we have investigated the origin of the systematic variation in horizontal-branch (HB) morphology with galactocentric distance (R(sub G)) among globular clusters. The variations in He abundance, CNO abundance, and core mass required separately to explain this effect are inconsistent with either the observed properties of the RR Lyrae variables or the observed main-sequence turnoffs in the clusters. There is also no clear evidence that the trend with R(sub G) is related to the central concentrations, central densities, or absolute magnitudes of the clusters. The variations in cluster age required to explain this effect are not in conflict with any observations. A detailed comparison of our synthetic HB calculations with pairs of clusters of very different HB morphology but similar (Fe/H) reveals reasonably good agreement between the age differences inferred from HB morphology and the main-sequence turnoff. The major source of uncertainty is the need for ad hoc hypotheses in the modeling of the HB morphologies of a few peculiar clusters (e.g., NGC 6752). Nonetheless, there is firm evidence for age variations of several gigayears (as much as approximately 5 Gyr) among the halo globular clusters. Our results support the hypothesis of Searle & Zinn that the inner halo is more uniform in age and is older in the mean than the outer halo, and we estimate this difference to be approximately 2 Gyr.

Lee, Young-Wook↗

NASA NTRS1994

Mass loss during the RR Lyrae phase of the horizontal branch: Mass dispersion on the horizontal branch and RR Lyrae period changes

Mass loss on the horizontal branch has been invoked in the literature to explain such phenomena as the color (mass) dispersion of the horizontal branch and the observed distribution of period changes in RR Lyrae stars. To test these claims, the Yale stellar evolution code was used to evolve horizontal branch models of masses 0.64, 0.66, 0.68, 0.70, and 0.72 solar mass with Z of 0.001, core mass of 0.4893, main-sequence helium abundance of 0.23, and constant mass loss rates of 0, 10(exp -10), 5 x 10(exp -10), and 10(exp -9) solar mass/yr. Mass loss was assumed to occur only in the instability strip, where a mechanism is most likely to exist. Synthetic horizontal branches, constructed from the models, show that mass loss on the horizontal branch cannot produce the observed color dispersion even for the highest mass-loss rate of 10(exp -9) solar mass/yr. Mass loss is unlikely to occur at a higher rate without significant effects on the horizontal branch morphology, which would destroy the good agreement between standard synthetic models without mass loss and observed horizontal branches. Periods and period changes were calculated for all models. The period changes are not significantly larger for models with mass loss. The effect of mass loss in clusters of other metallicities is discussed.

Koopmann, Rebecca A.↗

NASA NTRS1993

The effect of Livermore OPAL opacities on the evolutionary masses of RR Lyrae stars

We have investigated the effect of the new Livermore OPAL opacities on the evolution of horizontal-branch (HB) stars. This work was motivated by the recent stellar pulsation calculations using the new Livermore opacities, which suggest that the masses of double-mode RR Lyrae stars are 0.1-0.2 solar mass larger than those based on earlier opacities. Unlike the pulsation calculations, we find that the effect of opacity change on the evolution of HB stars is not significant. In particular, the effect of the mean masses of RR Lyrae stars is very small, showing a decrease of only 0.01-0.02 solar mass compared to the models based on old Cox-Stewart opacities. Consequently, with the new Livermore OPAL opacities, both the stellar pulsation and evolution models now predict approximately the same masses for the RR Lyrae stars. Our evolutionary models suggest that the mean masses of the RR Lyrae stars are about 0.76 and about 0.71 solar mass for M15 (Oosterhoff group II) and M3 (group I), respectively. If (alpha/Fe) = 0.4, these values are decreased by about 0.03 solar mass. Variations of the mean masses of RR Lyrae stars with HB morphology and metallicity are also presented.

Yi, Sukyoung↗

NASA NTRS1992

Ages of globular clusters and helium diffusion

Evolutionary tracks have been calculated with alpha-enhanced compositions which cover the entire globular cluster metallicity range and have constructed isochrones which include the effects of microscopic diffusion of helium. The turnoff magnitudes from the isochrones were combined with the theoretical RR Lyrae magnitudes from Lee to determine the ages of 32 Galactic globular clusters using the magnitude difference between the turnoff and horizontal branch. It is found that including the effects of helium diffusion has a negligible effect on the derived ages of globular clusters. Regardless of the inclusion of helium diffusion, a significant age spread of 5 Gyr among the globular clusters is found. The oldest globular clusters studied here are 17 +/- 2 Gyr old.

Chaboyer, Brian↗

NASA NTRS1992

Evolutionary models and the p-mode oscillation spectrum of Alpha Centauri A and B

Spectroscopically observed abundances of Furenlid and Meylan (1990) and Los Alamos Opacity Library opacities specifically constructed for this mixture are used to construct the most detailed models to date for Alpha Cen A and B. Models including the effects of the diffusion of helium in Alpha Cen A's interior were constructed. Although the effect of helium diffusion is small in Alpha Cen A, it is slightly larger than in the sun because of the larger mass and therefore shallower convection zone of Alpha Cen A. Keeping the heavy-element mixture unchanged but varying Z within the uncertainties yields limits to the helium abundance of Y = 0.295 and 0.305, respectively. The derived age of Alpha Cen is 4.6 +/-0.4 Gyr, very similar to the solar age of 4.5 Gyr.

Edmonds, Peter↗