Search

The high quality spectra required for radial velocity planet searches are well-suited to providing abundances for a wide array of elements in large samples of stars. Abundance ratios of the most common elements relative to Fe are observed to vary by more than a factor of two in planet host candidates. This level of variation has a substantial impact on the evolution of the host star and the extent of its habitable zone. We present stellar models of 1 solar mass stars with custom compositions representing the full range of these non-solar abundance ratios. We find that the effects derived from variation over the observed range of [O/Fe] are particularly dramatic. Habitability lifetimes for some classes of orbits can vary by gigayears for the observed range in [O/Fe]., Comment: 3 figures, accepted for publication in the Astrophysical Journal Letters

It is suggested in observations of supernova remnants that a number of large- and small-scale structures form at various points in the explosion. Multidimensional modeling of core-collapse supernovae has been undertaken since SN1987A, and both simulations and observations suggest/show that Rayleigh-Taylor instabilities during the explosion is a main driver for the formation of structure in the remnants. We present a case study of structure formation in 3D in a \msol{15} supernova for different parameters. We investigate the effect of moderate asymmetries and different resolutions of the formation and morphology of the RT unstable region, and take first steps at determining typical physical quantities (size, composition) of arising clumps. We find that in this progenitor the major RT unstable region develops at the He/OC interface for all cases considered. The RT instabilities result in clumps that are overdense by 1-2 orders of magnitude with respect to the ambient gas, have size scales on the level of a few % of the remnant diameter, and are not diffused after the first $\sim30$ yrs of the remnant evolution, in the absence of a surrounding medium., Comment: 59 pages, 34 figures

An algorithm for the direct inversion of the linear systems arising from Nystrom discretization of integral equations on one-dimensional domains is described. The method typically has O(N) complexity when applied to boundary integral equations (BIEs) in the plane with non-oscillatory kernels such as those associated with the Laplace and Stokes' equations. The scaling coefficient suppressed by the "big-O" notation depends logarithmically on the requested accuracy. The method can also be applied to BIEs with oscillatory kernels such as those associated with the Helmholtz and Maxwell equations; it is efficient at long and intermediate wave-lengths, but will eventually become prohibitively slow as the wave-length decreases. To achieve linear complexity, rank deficiencies in the off-diagonal blocks of the coefficient matrix are exploited. The technique is conceptually related to the H and H^2 matrix arithmetic of Hackbusch and co-workers, and is closely related to previous work on Hierarchically Semi-Separable matrices.

We compare the yields of Ti44 and Ni56 produced from post-processing the thermodynamic trajectories from three different core-collapse models -- a Cassiopeia A progenitor, a double shock hypernova progenitor, and a rotating 2D explosion -- with the yields from exponential and power-law trajectories. The peak temperatures and densities achieved in these core-collapse models span several of the distinct nucleosynthesis regions we identify, resulting in different trends in the Ti44 and Ni56 yields for different mass elements. The Ti44 and Ni56 mass fraction profiles from the exponential and power-law profiles generally explain the tendencies of the post-processed yields, depending on which regions are traversed by the model. We find integrated yields of Ti44 and Ni56 from the exponential and power-law trajectories are generally within a factor 2 or less of the post-process yields. We also analyze the influence of specific nuclear reactions on the Ti44 and Ni56 abundance evolution. Reactions that affect all yields globally are the 3a, p(e-,nu)n and n(e+,nubar)p. The rest of the reactions are ranked according to their degree of impact on the synthesis of Ti44. The primary ones include Ti44(a,p)V47, Ca40(a,g)Ti44, V45(p,g)Cr46, Ca40(a,p)Sc43, F17(a,p)Ne20, Na21(a,p)Mg24, Sc41(p,g)Ti42, Sc43(p,g)Ti44, Ti44(p,g)V45, and Ni57(p,g)Cu58, along with numerous weak reactions. Our analysis suggests that not all Ti44 need be produced in an a-rich freeze-out in core-collapse events, and that reaction rate equilibria in combination with timescale effects for the expansion profile may account for the paucity of Ti44 observed in supernovae remnants., Comment: Accepted for publication in the Astrophysical Journal Supplement

A scheme for rapidly and accurately computing solutions to boundary integral equations (BIEs) on rotationally symmetric surfaces in three dimensions is presented. The scheme uses the Fourier transform to reduce the original BIE defined on a surface to a sequence of BIEs defined on a generating curve for the surface. It can handle loads that are not necessarily rotationally symmetric. Nystrom discretization is used to discretize the BIEs on the generating curve. The quadrature used is a high-order Gaussian rule that is modified near the diagonal to retain high-order accuracy for singular kernels. The reduction in dimensionality, along with the use of high-order accurate quadratures, leads to small linear systems that can be inverted directly via, e.g., Gaussian elimination. This makes the scheme particularly fast in environments involving multiple right hand sides. It is demonstrated that for BIEs associated with Laplace's equation, the kernel in the reduced equations can be evaluated very rapidly by exploiting recursion relations for Legendre functions. Numerical examples illustrate the performance of the scheme; in particular, it is demonstrated that for a BIE associated with Laplace's equation on a surface discretized using 320 000 points, the set-up phase of the algorithm takes 2 minutes on a standard desktop, and then solves can be executed in 0.5 seconds., Comment: 28 pages, 4 figures

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy., Comment: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the stellar structure which results is similar to that obtained by YREC and Garching models, the theoretical picture differs. A new quantitative connection is made between macro-turbulence, micro-turbulence, and the convective velocity scale at the photosphere, which has finite values. The "geometric parameter" in MLT is found to correspond more reasonably with the size of the strong downward plumes which drive convection (Stein and Nordlund 1998), and thus has a physical interpretation even in MLT. Use of 3D simulations of both adiabatic convection and stellar atmospheres will allow the determination of the dissipation length and the geometric parameter (i.e., the entropy jump), with no astronomical calibration. A physically realistic treatment of convection in stellar evolution will require additional modifications beyond MLT, including effects of kinetic energy flux, entrainment (the most dramatic difference from MLT found by Meakin and Arnett 2007), rotation, and magnetic fields (Balbus 2009}., Comment: 16 pages, 5 figures, submitted to ApJ

Astronomers have proposed a number of mechanisms to produce supernova explosions. Although many of these mechanisms are now not considered primary engines behind supernovae, they do produce transients that will be observed by upcoming ground-based surveys and NASA satellites. Here we present the first radiation-hydrodynamics calculations of the spectra and light curves from three of these "failed" supernovae: supernovae with considerable fallback, accretion induced collapse of white dwarfs, and energetic helium flashes (also known as type .Ia supernovae)., Comment: 33 pages, 14 figures

We report extensive radial-velocity measurements of the two giant components of the detached, 104-day period binary system of Capella. Our highly accurate three-dimensional orbital solution based on all existing spectroscopic and astrometric observations including our own yields much improved masses of 2.466 +/- 0.018 M_Sun and 2.443 +/- 0.013 M_Sun for the primary and secondary (relative errors of 0.7% and 0.5%). Improved values are derived also for the radii (11.87 +/- 0.56 R_Sun and 8.75 +/- 0.32 R_Sun), effective temperatures (4920 +/- 70 K and 5680 +/- 70 K), and luminosities (79.5 +/- 4.8 L_Sun and 72.1 +/- 3.6 L_Sun). The distance is determined to be 13.042 +/- 0.028 pc. Capella is unique among evolved stars in that, in addition to all of the above, the chemical composition is known, including the overall metallicity [m/H], the carbon isotope ratio 12C/13C for the primary, and the lithium abundance and C/N ratios for both components. The latter three quantities are sensitive diagnostics of evolution, and change drastically for giants as a result of the deepening of the convective envelope during the first dredge-up. The secondary is crossing the Hertzprung gap, while the primary is believed to be in the longer-lived core-helium burning phase. However, we find that current stellar evolution models are unable to match all of the observations for both components at the same time, and at a single age, for any evolutionary state of the primary. Similar problems are found when testing the rotational synchronization, spin axis alignment, and orbital circularization of the system against tidal theory. We conclude that our understanding of the advanced stages of stellar evolution is still very incomplete. [Abridged], Comment: 42 pages in emulateapj format, including figures and tables. To appear in the 2009 August 10 issue of The Astrophysical Journal

We present new reddening measurements towards the young supernova remnant Cassiopeia A, using two techniques not previously applied to this object. Our observations of the near-infrared [Fe II] 1.257 micron and 1.644 micron lines show the extinction to be highly variable across the remnant, increasing towards the west and the south, consistent with previous radio and X-ray observations. While the absolute value of AV as determined by the [Fe II] lines is uncertain due to conflicting calculations and observations of their intrinsic flux ratio, parts of the remnant without previous optical measurements show comparatively higher reddening. We find AV = 6.2 +/- 0.6 from the broadband shape of the infrared synchrotron emission of a knot within 13" of the expansion center. Given this reddening, the apparent faintness of the nascent supernova, and iron mass constraints from X-ray observations, we estimate an ejected mass of 56Ni of 0.058 - 0.16 MSun. Taken with gamma-ray observations of the 44Ti decay chain, this nickel mass is broadly consistent with the solar 44Ca/56Fe ratio., Comment: 8 pages, 6 figures, accepted for publication in The Astrophysical Journal

Abridged: The golden age of astrophysics is upon us with both grand discoveries (extra-solar planets, dark matter, dark energy) and precision cosmology. Fundamental understanding of the working of stars and galaxies is within reach, thanks to newly available precision measurements. We highlight the importance of distances and model independent distances and masses. Distances are fundamental in astrophysics and their knowledge can change our perception of phenomena dramatically: e.g., in antiquity, the Heliocentric model was rejected because the predicted stellar parallaxes were not observed. Distance measurements are directly related to the history & fate of the uni- verse as they provide 2 of 3 methods available to date the universe. The 1st method is based on the ages of stars, which can be ascertained if their lumi- nosities (distances) are accurately known. The 2nd method relies on cosmolo- gical methods. To 1st order, the age of the universe is the inverse H_0. As stressed by the previous decadal report, "the fundamental goal of ... astrophysics is to understand how the universe ... galaxies [and] stars ... formed, how they evolved, and what their destiny will be." These questions can be answered partly by micro-arcsecond astrometry: 1) Galactic archeology: a detailed reconstruction of the formation history of the Milky Way and other Local Group galaxies, 2) the very oldest stars in the Milky Way and the age of the Universe, and 3) H_0 and concordance cosmology. These goals are achievable by combining muas-arcsecond astrometry from the proposed SIM-Lite mission supplemented with ground-based spectroscopy. The results of our proposed project will force the biggest reassessment of stellar astrophysics in 50 years, which will affect most branches of astrophysics., Comment: 8 pages; submission for the US 2010 Decadal Review

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low$-$mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole. Here we present the results of the first three-dimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion. The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields. We compare the nucleosynthetic yields to those of hypernovae. Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates. Especially in studies of integrated stellar yields, simplifications reduce these uncertainties. But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis. Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them., Comment: To appear in the Conference Proceedings for the "10th Symposium on Nuclei in the Cosmos (NIC X)", July 27 - August 1 2008, Mackinack Island, Michigan, USA

The currently-favored model for long-duration gamma-ray bursts (GRBs) invokes explosions from the collapse of a massive star down to a black hole: either directly or through fallback. Those GRBs forming via fallback will produce much less radioactive nickel, and hence it has been argued (without any real calculation) that these systems produce dim supernovae. These fallback black-hole GRBs have been recently been argued as possible progenitors of a newly discovered set of GRBs lacking any associated supernovae. Here we present the first ever radiation-hydrodynamics calculations of the light-curves produced in the hypernova explosion by a delayed-fallback gamma-ray burst. We find that the bolometric light-curve is dominated by shock-deposited energy, not the decay of radioactive elements. As such, observations of such bursts actually probe the density in the progenitor wind more than it does the production of radioactive nickel., Comment: 11 pages (including 3 figures), submitted to ApJ, comments welcome

We use population synthesis to construct distributions of gamma-ray bursts (GRBs) for different proposed progenitor models. We use a description of star formation that takes into account the evolution of metallicity with redshift and galaxy mass, the evolution of galaxy mass with redshift, and the star formation rate with galaxy mass and redshift. We compare predicted distributions with redshift and metallicity to observations of GRB host galaxies and find that the the simple models cannot produce the observed distributions, but that current theoretical models can reproduce the observations within some constraints on the fraction of fallback black holes that produce GRBs., Comment: 17 pages, 9 figures, submitted to the Astrophysical Journal

Accurate stellar ages remain one of the most poorly constrained, but most desired, astronomical quantities. Here we briefly summarize some recent efforts to improve the stellar age scale from a subset of talks from the ``Stellar Ages'' splinter session at the "14th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun". The topics discussed include both the apparent successes and alarming discrepancies in using Li depletion to age-date clusters, sources of uncertainty in ages due to input physics in evolutionary models, and recent results from asteroseismology and gyrochronology., Comment: 8 pages, 2 figures, Summary of splinter session on "Stellar Ages" for the 14th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, 6-10 November 2006

Theoretical nucleosynthetic yields from supernovae are sensitive to both the details of the progenitor star and the explosion calculation. We attempt to comprehensively identify the sources of uncertainties in these yields. In this paper we concentrate on the variations in yields from a single progenitor arising from common 1-dimensional methods of approximating a supernova explosion. Subsequent papers will examine 3-dimensional effects in the explosion and the progenitor, and trends in mass and composition. For the 1-dimensional explosions we find that both elemental and isotopic yields for Si and heavier elements are a sensitive function of explosion energy. Also, piston-driven and thermal bomb type explosions have different yields for the same explosion energy. Yields derived from 1-dimensional explosions are non-unique., Comment: 87 pages, 22 figures, subitted to the Astrophysical Journal

The collapsar engine behind long-duration gamma-ray bursts extracts the energy released from the rapid accretion of a collapsing star onto a stellar-massed black hole. In a collapsing star, this black hole can form in two ways: the direct collapse of the stellar core into a black hole and the delayed collapse of a black hole caused by fallback in a weak supernova explosion. In the case of a delayed-collapse black hole, the strong collapsar-driven explosion overtakes the weak supernova explosion before shock breakout, and it is very difficult to distinguish this black hole formation scenario from the direct collapse scenario. However, the delayed-collapse mechanism, with its double explosion, produces explosive nucleosynthetic yields that are very different from the direct collapse scenario. We present 1-dimensional studies of the nucleosynthetic yields from both black hole formation scenarios, deriving differences and trends in their nucleosynthetic yields., Comment: 47 pages, submitted to ApJ

Gamma-ray burst (GRB) afterglow observations have allowed us to significantly constrain the engines producing these energetic explosions. Te redshift and position information provided by these afterglows have already allowed us to limit the progenitors of GRBs to only a few models. The afterglows may also provide another observation that can place further constraints on the GRB progenitor: measurements telling us about the environments surrounding GRBs. Current analyses of GRB afterglows suggest that roughly half of long-duration gamma-ray bursts occur in surroundings with density profiles that are uniform. We study the constraints placed by this observation on both the classic ``collapsar'' massive star progenitor and its relative, the ``helium-merger'' progenitor. We study several aspects of wind mass-loss and find that our modifications to the standard Wolf-Rayet mass-loss paradigm are not sufficient to produce constant density profiles. Although this does not rule out the standard ``collapsar'' progenitor, it does suggest a deficiency with this model. We then focus on the He-merger models and find that such progenitors can fit this particular constraint well. We show how detailed observations can not only determine the correct progenitor for GRBs, but also allow us to study binary evolution physics., Comment: 44 pages including 11 figures

We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints., Comment: 31 pages, 4 figures, accepted for publication to the Astrophysical Journal

The Sirius AB binary system has masses that are well determined from many decades of astrometric measurements. Because of the well-measured radius and luminosity of Sirius A, we employed the TYCHO stellar evolution code to determine the age of the Sirius A,B binary system accurately, at 225-250 Myr. Note that this fit requires the assumption of solar abundance, and the use of the new Asplund et al. primordial solar metallicity. No fit to Sirius A's position is possible using the old Grevesse & Sauval scale. Because the Sirius B white dwarf parameters have also been determined accurately from space observations, the cooling age could be determined from recent calculations by Fontaine et al. or Wood to be 124 +/- 10 Myr. The difference of the two ages yields the nuclear lifetime and mass of the original primary star, 5.056 +0.374/-0.276 solar masses. This result yields in principle the most accurate data point at relatively high masses for the initial-final mass relation. However, the analysis relies on the assumption that the primordial abundance of the Sirius stars was solar, based on membership in the Sirius supercluster. A recent study suggests that its membership in the group is by no means certain., Comment: 12 pages, one figure (reduced resolution JPG posted), accepted for publication in the ApJ Letters. Replaced with revised version correcting a typo in the text (the derived age of Sirius A on page 5); conclusions unchanged

Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially., Comment: 12 pages, 2 figures, accepted for publication in the Astrophysical Journal Letters

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test., Comment: 33 pages, 8 figures, accepted for publication in the Astrophysical Journal

The detection of NaCl at large radii in the Egg Nebula, CRL 2688, requires densities of 10^7 - 10^8 cm^-3 in a thick shell of r ~ a few X 10^17 cm. To explain these results, a mechanism is needed for producing high densities at a considerable distance from the central star. In two dimensional simulations of the interaction of the fast wind with an ambient medium, the material becomes thermally unstable. The resulting clumps can achieve the requisite conditions for NaCl excitation. We present 2D models with simple physics as proof-of-principle calculations to show that the clumping behavior is robust. Clumping is a natural outcome of cooling in the colliding wind model and comparable to that inferred from observations., Comment: 12 pages, 3 figures. Accepted for publication in the Astrophysical Journal Letters

We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of ``overshooting'' gives extra mixing relative to the treatment traditionally used in stellar evolutionary codes. The interface between convectively stable and unstable regions is dynamic and nonspherical, so that the nonturbulent material is driven into motion, even in the absence of ``penetrative overshoot.'' These motions may be described by the theory of nonspherical stellar pulsations, and are related to motion measured by helioseismology. Multi-dimensional numerical simulations of convective flow show puzzling features which we explain by this simplified physical model. Gravity waves generated at the interface are dissipated, resulting in slow circulation and mixing seen outside the formal convection zone. The approach may be extended to deal with rotation and composition gradients. Tests of this description in the stellar evolution code TYCHO produce carbon stars on the asymptotic giant branch (AGB), an isochrone age for the Hyades and three young clusters with lithium depletion ages from brown dwarfs, and lithium and beryllium depletion consistent with observations of the Hyades and Pleiades, all without tuning parameters. The insight into the different contributions of rotational and hydrodynamic mixing processes could have important implications for realistic simulation of supernovae and other questions in stellar evolution., Comment: 27 pages, 5 figures, accepted to the Astrophysical Journal