The onset of cool massive winds in evolved giants is correlated with an evolutionary feature on the red giant branch known as the bump. Also at the bump, shear instability in the star leads to magnetic fields that occur preferentially on small length scales. Pneuman (1983) has suggested that the emergence of small scale flux tubes in the Sun can give rise to enhanced acceleration of the solar wind as a result of plasmoid acceleration (the melon seed mechanism). In this paper, we examine the Pneuman formalism to determine if it may shed some light on the process that drives mass loss from stars above the bump. Because we do not currently have detailed information for some of the relevant physical parameters, we are not yet able to derive a detailed model. Instead, our goal in this paper is to explore a proof of concept. Using parameters that are known to be plausible in cool giants, we find that the total mass loss rate from such stars can be replicated. Moreover, we find that the radial profile of the wind speed in such stars can be steep or shallow depending on the fraction of the mass loss which is contained in the plasmoids. This is consistent with empirical data which indicate that the velocity profiles of winds from cool giants range from shallow to steep.

It has been suggested that the absence of macroturbulence in the atmosphere of NGC 1624 - 2 is due its strong magnetic field (the strongest known for a massive O star) suppressing convection in its outer layers, removing the mechanism thought responsible for the observed macroturbulence in stars with lower field strengths. Here, we develop and apply a criterion for a uniform magnetic field to suppress convection in stellar envelopes in which radiation pressure is a significant contributor to hydrostatic balance. We find upper mass limits of ~55 Msun and ~30 Msun for magnetic suppression to be possible in Zero Age Main Sequence and Terminal Age Main Sequence stars, respectively. For evolved stars, magnetic suppression of convection can significantly alter the stars' evolution. For NGC 1624 - 2, we find a polar dipole strength of 16.5 +/- 5.9 kG is required to suppress convection, in good agreement with the value ~20 kG measured by spectropolarimetry.

Jurcak et al (2018) have reported that, in a sample of more than 100 umbral cores in sunspots, the umbral-penumbral boundary (UPB) is characterized by a remarkably narrowly-defined numerical value (1867 G) of the vertical component of the magnetic field. Gough and Tayler (1966), in their study of magneto-convection, showed that the onset of convection in the presence of a magnetic field is controlled by a parameter {\delta} which also depends on the vertical component of the field. Combining the Jurcak et al result with various empirical models of sunspots leads us to propose the following hypothesis: the UPB occurs where the vertical field is strong enough to increase the effective adiabatic temperature gradient by at least 100% above its non-magnetic value.

We construct evolutionary models of Trappist-1 in which magnetic fields impede the onset of convection according to a physics-based criterion. In the models that best fit all observational constraints, the photospheric fields in Tr-1 are found to be in the range 1450-1700 G. These are weaker by a factor of about 2 than the fields we obtained in previous magnetic models of two other cool dwarfs (GJ65A/B). Our results suggest that Tr-1 possesses a global poloidal field which is some one hundred times stronger than in the Sun. In the context of exoplanets in orbit around Tr-1, the strong poloidal fields on the star may help to protect the planets from the potentially destructive effects of coronal mass ejections. This, in combination with previous arguments about beneficial effects of flare photons in ultraviolet and visible portions of the spectrum, suggests that conditions on Tr-1 are not necessarily harmful to life on a planet in the habitable zone of Tr-1.

We show that the recently discovered narrow gap in the Gaia Hertzsprung - Russell Diagram near MG =10 can be explained by standard stellar evolution models and results from a dip in the luminosity function associated with mixing of 3He during merger of envelope and core convection zones that occurs for a narrow range of masses.

GJ 65 is a M dwarf binary system consisting of the two flare stars BL Cet (GJ 65A) and UV Cet (GJ 65B). Two teams of investigators have recently reported total magnetic fluxes corresponding to fields of 4.5 and 5.2 kG for GJ65A, and 5.8 and 6.7 kG for GJ65B: for each component, the magnetic results obtained by the two teams agree with each other within 1{\sigma}. For the first time, we can directly compare the predictions of our magneto-convective models, based on fitting observed stellar parameters, with measured field strengths. We find that our models agree with the observed field strengths provided the effects of finite conductivity are accounted for. Thus, GJ65 provides us an opportunity to use observations of field strengths to distinguish between the predictions of our models that assume perfect electrical conductivity and those that allow for finite conductivity.

We have previously reported on models of low-mass stars in which the presence of inflated radii is ascribed to magnetic fields which impede the onset of convection (e.g. MacDonald & Mullan [2017a] and citations therein). Some of our magneto-convection models have been criticized because, when they were first reported by Mullan & MacDonald (2001), the deep interior fields were found to be very large (50-100 MG). Such large fields are now known to be untenable. For example, Browning et al. (2016) used stability arguments to suggest that interior fields in low-mass stars cannot be larger than ~1 MG. Moreover, 3D models of turbulent stellar dynamos suggest that fields generated in low-mass interiors may be not much stronger than 10-20 kG (Browning 2008). In the present paper, we present magneto-convective models of inflated low-mass stars in which the interior fields are not permitted to be stronger than 10 kG. These models are used to fit empirical data for 15 low-mass stars for which precise masses and radii have been measured. We show that our 10 kG magneto-convective models can replicate the empirical radii and effective temperatures for 14 of the stars. In the case of the remaining star (in the Praesepe cluster), two different solutions have been reported in the literature. We find that one of these solutions (by Gillen et al. 2017) can be fitted well with our model using the nominal age of Praesepe (800 Myr). However, the second solution (by Kraus et al. 2017) cannot be fitted unless the star's age is assumed to be much younger (~ 150 Myr).

Stellar oscillations appear all across the Hertzsprung-Russell diagram. Recent theoretical studies support their existence also in the atmospheres of M dwarfs. These studies predict for them short periodicities ranging from 20~min to 3~h. Our Cool Tiny Beats (CTB) programme aims at finding these oscillations for the very first time. With this goal, CTB explores the short time domain of M dwarfs using radial velocity data from the HARPS-ESO and HARPS-N high-precision spectrographs. Here we present the results for the two most long-term stable targets observed to date with CTB, GJ~588 and GJ~699 (i.e. Barnard's star). In the first part of this work we detail the correction of several instrumental effects. These corrections are specially relevant when searching for sub-night signals. Results show no significant signals in the range where M dwarfs pulsations were predicted. However, we estimate that stellar pulsations with amplitudes larger than $\sim0.5\,\mathrm{m\,s}^{-1}$ can be detected with a 90% completeness with our observations. This result, along with the excess of power regions detected in the periodograms, open the possibility of non-resolved very low amplitude pulsation signals. Next generation more precise instrumentation would be required to detect such oscillations. However, the possibility of detecting pulsating M-dwarf stars with larger amplitudes is feasible due to the short size of the analysed sample. This motivates the need for completeness of the CTB survey.

LSPM J1314+1320 (=NLTT 33370) is a binary star system consisting of two nearly identical pre-main sequence stars of spectral type M7. The system is remarkable among ultracool dwarfs for being the most luminous radio emitter over the widest frequency range. Masses and luminosities are at first sight consistent with the system being coeval at age ~80 Myr according to standard (non-magnetic) evolutionary models. However, these models predict an average effective temperature, 2950 +/- 5 K, which is 180 K hotter than the empirical value. Thus, the empirical radii are oversized relative to the standard models by ~13%. We demonstrate that magnetic stellar models can account quantitatively for the oversizing. As a check on our models, we note that the radio emission limits the surface magnetic field strengths: the limits depend on identifying the radio emission mechanism. We find that the field strengths required by our magnetic models are too strong to be consistent with gyrosynchrotron emission, but are consistent with electron cyclotron maser emission.

Two eclipsing binaries in the USco association have recently yielded precise values of masses and radii for 4 low-mass members of the association. Standard evolution models would require these dM4.5 - dM5 stars to have ages which are younger than the ages of more massive stars in the association by factors which appear (in extreme cases) to be as large as ~3. Are the stars in the association therefore non-coeval? We suggest that the answer is No: by incorporating the effects of magnetic inhibition of convective onset, we show that the stars in USco can be restored to coevality provided that the 4 low-mass member stars have vertical surface fields in the range 200 - 700 G. Fields of such magnitude have already been measured on the surface of certain solar-type stars in other young clusters.

KIC 7177553 is a quadruple system containing two binaries of orbital periods 16.5 and 18 d. All components have comparable masses and are slowly rotating, non-evolved stars of spectral type near G2V. The longer period binary is eclipsing with component masses and radii, M1 = 1.043 +/- 0.014 M_sun, R1 = 0.940 +/- 0.005 R_sun, and M2 = 0.986 +/- 0.015 M_sun, R2 = 0.941 +/- 0.005 R_sun. The essentially equal radii measurements are inconsistent with the two stars being on the main sequence at the same age using standard stellar evolution models. Instead a consistent scenario is found if the stars are in their pre-main sequence phase of evolution and have age 33 - 36 Myr. Such a young age has important implications for the detectability of the massive planet indicated by eclipse time variations.

Empirically, the X-ray luminosity LX from M dwarfs has been found to have an upper limit of about 0.2% of the bolometric flux Lbol. In the limit where magnetic fields in M dwarfs are generated in equipartition with convective motions, we use stellar models to calculate the energy flux of Alfven waves FA as a function of depth in the sub-surface convection zone. Since Alfven waves have the optimal opportunity for wave modes to reach the corona, we suggest that FA sets an upper limit on the mechanical flux Fmech which causes coronal heating. This suggestion accounts quantitatively for the saturated values of LX/Lbol which have been reported empirically for M dwarfs.

The first stars continue to elude modern telescopes, but much has been accomplished in observing the glow of the first galaxies. As detection capabilities improve we will eventually resolve these galaxies, but hopes of observing an individual star remains dim for the foreseeable future. However, our first view of an individual first star may be possible due to its explosion. In this work, we present evolution calculations for Population III (Pop III) stars and their subsequent supernovae explosions. Our evolution models include a mass range of 15 - 100 M$_\odot$, each with initial heavy element abundance Z = 10-14. Our models are evolved from pre-main sequence through formation of an iron core, and thus near to core collapse. We find that modelling the evolution of these stars is very sensitive to the choice of convection criterion; here we provide evolution results using both the Schwarzschild and Ledoux criteria. We also use the final structure from our models for numerical simulation of their supernovae light curves using a radiation hydrodynamics code. In doing so, we estimate a lower bound of initial model mass that may be possible to observe in the near future. We also find that our 40 M$_\odot$ Schwarzschild evolution model produces the brightest supernova peak and statistically should be the most frequently observed. At our highest redshift z= 15, only the 60 M$_\odot$ Schwarzschild model at peak magnitude starts to rival the 40 M$_\odot$ model in brightness.

Magnetic fields impede the onset of convection, thereby altering the thermal structure of a convective envelope in a low mass star: this has an effect on the amount of lithium depletion in a magnetized star. In order to quantify this effect, we have applied a magneto-convective model to two low mass stars for which lithium abundances and precise structural parameters are known: YY Gem and CU Cnc. For both stars, we have obtained models which satisfy empirical constraints on the following parameters: R, L, surface magnetic field strength, and Li abundance. In the case of YY Gem, we have obtained a model which satisfies the empirical constraints with an internal magnetic field of several megagauss: such a field strength is within the range of a dynamo where the field energy is in equipartition with rotational energy deep inside the convection zone. However, in the case of CU Cnc, the Li requires an internal magnetic field which is probably too strong for a dynamo origin: we suggest possible alternatives which might account for the reported Li abundance in CU Cnc.

We present the analysis of four M dwarf stars -plus one M giant that seeped past our selection criteria- observed in Cycle 3 of Kepler Guest Observer program (GO3) in a search for intrinsic pulsations. Stellar oscillations in M dwarfs were theoretically predicted by Rodr\'iguez-L\'opez et al. (2012) to be in the range ~20-40 min and ~4-8 h, depending on the age and the excitation mechanism. We requested Kepler short cadence observations to have an adequate sampling of the oscillations. The targets were chosen on the basis of detectable rotation in the initial Kepler results, biasing towards youth.The analysis reveals no oscillations attributable to pulsations at a detection limit of several parts per million, showing that either the driving mechanisms are not efficient in developing the oscillations to observable amplitudes, or that if pulsations are driven, the amplitudes are very low. The size of the sample, and the possibility that the instability strip is not pure, allowing the coexistence of pulsators and non-pulsators, prevent us from deriving definite conclusions. Inmediate plans include more M dwarfs photometric observations of similar precision with Kepler K2 mission and spectroscopic searches already underway within the Cool Tiny Beats Project (Anglada-Escud\'e et al. 2014, Berdi\~nas et al. 2014) with the high-resolution spectrographs HARPS and HARPS-N.

The overstability of the fundamental radial mode in M dwarf models was theoretically predicted by Rodr\'iguez-L\'opez et al. (2012). The periods were found to be in the ranges ~25-40 min and ~4-8 h, depending on stellar age and excitation mechanism. We have extended our initial M dwarf model grid in mass, metallicity, and mixing length parameter. We have also considered models with boundary conditions from PHOENIX NextGen atmospheres to test their influence on the pulsation spectra. We find instability of non-radial modes with radial orders up to k=3, degree l=0-3, including p and g modes, with the period range extending from 20 min up to 11 h. Furthermore, we find theoretical evidence of the potential of M dwarfs as solar-like oscillators.

Precision modeling of M dwarfs has become worthwhile in recent years due to the increasingly precise values of masses and radii which can be obtained from eclipsing binary studies. Torres (2013) has identified 4 prime M dwarf pairs with the most precise empirical determinations of masses and radii. The measured radii are consistently larger than standard stellar models predict. We have previously modeled M dwarfs in the context of a criterion due to Gough & Tayler in which magnetic fields inhibit the onset of convection according to a physics-based prescription. New constraints on the models of M dwarfs are now provided by measurements of lithium abundances. The key aspect of Li in terms of setting constraints on magnetic modeling is that Li burning starts at T = 2.5 MK, and temperatures of just such magnitude are associated with the base of the convection zone: magnetic inhibition of convective onset can shift this base slightly closer to the surface, i.e. to slightly lower temperatures, thereby reducing the amount of Li depletion compared to a non-magnetic model. In the present paper, we consider how our magneto-convection models handle the new test of stellar structure provided by Li measurements. Among the prime systems listed by Torres, we find that plausible magnetic models work well for CM Dra and YY Gem but not for CU Cnc. (The fourth system in Torres's list does not yet have enough information to warrant magnetic modeling.) For CU Cnc, we suggest that the observed lithium may have been accreted from a circumstellar disk. We find that our magneto-convection models of CM Dra, YY Gem and CU Cnc yield results which are consistent with the observed correlation between magnetic flux and X-ray luminosity.

By combining rotation periods with spectroscopic determinations of projected rotation velocity, Jackson, Jeffries & Maxted (2009) have found that the mean radii for low-mass M-dwarfs in the young, open cluster NGC 2516 are larger than model predictions at a given absolute I magnitude or I - K color and also larger than measured radii of magnetically inactive M-dwarfs. The relative radius difference is correlated with magnitude, increasing from a few per cent at MI = 7 to greater than 50 per cent for the lowest luminosity stars in their sample at MI about 9.5. Jackson et al (2009) have suggested that a two-temperature star spot model is capable of explaining the observations, but their model requires spot coverage fractions of at least 50 per cent in rapidly rotating M-dwarfs. Here we examine these results in terms of stellar models that include the inhibiting effects of magnetic fields on convective energy transport, with and without the effects of star spots. We find that a pure spot model is inconsistent with the color - magnitude diagram. The observations of radii versus color and radii versus absolute magnitude in NGC 2516 are consistent with models which include only magnetic inhibition or a combination of magnetic inhibition and spots. At a given mass we find a large dispersion in the strength of the vertical component of the magnetic field in the stellar photosphere but the general trend is that the vertical field increases with decreasing mass from a few hundred Gauss at 0.65 Msun to 600 - 900 Gauss, depending on spot coverage, in the lowest mass stars in the sample at 0.25 Msun.

Spectroscopic analysis of the Galactic halo star SDSS J102915+172927 has shown it to have a very low heavy element abundance, Z < 7.4 10-7, with [Fe/H] = -4.89 plus/minus 0.10 and an upper limit on the C abundance of [C/H] < -4.5. The low C/Fe ratio distinguishes this object from most other extremely metal poor stars. The effective temperature and surface gravity have been determined to be Teff = 5811 plus/minus 150 K and log g = 4.0 plus/minus 0.5. The surface gravity estimate is problematical in that it places the star between the main sequence and the subgiants in the Hertzsprung-Russell diagram. If it is assumed that the star is on the main sequence, its mass and are estimated to be M = 0.72 plus/minus 0.06 Msun and L = 0.45 plus/minus 0.10 Lsun, placing it at a distance of 1.35 plus/minus 0.16 kpc. The upper limit on the lithium abundance, A(Li) < 0.9, is inconsistent with the star being a dwarf, assuming that mixing is due only to convection. In this paper, we propose that SJ102915 is a sub-giant that formed with significantly higher Z than currently observed, in agreement with theoretical predictions for the minimum C and/or O abundances needed for low mass star formation. In this scenario, extremely low Z and low Li abundance result from gravitational settling on the main sequence followed by incomplete convective dredge-up during subgiant evolution. The observed Fe abundance requires the initial Fe abundance to be enhanced compared to C and O, which we interpret as formation of SJ102915 occurring in the vicinity of a type Ia supernova.

We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models are then used to determine the rate at which objects are leaving the main sequence, which must be the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result is an empirically derived relation between temperature ($T_{eff}$) and time ($t$) on the white dwarf cooling sequence. Comparing this result to theoretical cooling models, we find general agreement with the expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition to the Mestel cooling rate of $T_{eff} \propto t^{-0.4}$ is found to occur at hotter temperatures, and more abruptly than is predicted by any of these models.

We present the results of the first theoretical non-radial non-adiabatic pulsational study of M dwarf stellar models with masses in the range 0.1 to 0.5M_solar. We find the fundamental radial mode to be unstable due to an \epsilon mechanism caused by deuterium (D-) burning for the young 0.1 and 0.2M_solar models, by non-equilibrium He^3 burning for the 0.2 and 0.25M_solar models of 10^4Myr, and by a flux blocking mechanism for the partially convective 0.4 and 0.5M_solar models once they reach the age of 500Myr. The periods of the overstable modes excited by the D-burning are in the range 4.2 to 5.2h for the 0.1M_solar models and is of order 8.4h for the 0.2M_solar models. The periods of the modes excited by He^3 burning and flux blocking are in the range 23 to 40min. The more massive and oldest models are more promising for the observational detection of pulsations, as their ratio of instability e-folding time to age is more favourable.

The eclipsing binary CM Dra contains nearly identical red dwarfs of spectral class dM4.5. Their masses and radii have been reported with unprecedentedly small statistical errors. When compared with standard stellar models of appropriate age (\approx4 Gyr), the empirical results indicate that both components are larger in R, and lower in luminosity L, by several standard deviations. Here, we attempt at first to model the components in the context of standard stellar models using differ

Although the means of the ages of stars in young groups determined from Li depletion often agree with mean ages determined from Hertzsprung - Russell diagram isochrones, there are often statistically significant differences in the ages of individual stars determined by the two methods. We find that inclusion of the effects of inhibition of convection due to the presence of magnetic fields leads to consistent ages for the individual stars. We illustrate how age consistency arises by applying our results to the \beta Pictoris moving group. We find that, although magnetic inhibition of convection leads to increased ages from the Hertzsprung - Russell diagram isochrones for all stars, Li ages are decreased for fully convective M stars and increased for stars with radiative cores. Our consistent age determination for the \beta Pictoris moving group of 40 Myr is larger than previous determinations by a factor of about two.

We present the first description of mode trapping for sdO models. Mode trapping of gravity modes caused by the He/H chemical transition is found for a particular model, providing a selection effect for high radial order trapped modes. Low- and intermediate-radial order {\em p}-modes (mixed modes with a majority of nodes in the P-mode region) are found to be trapped by the C-O/He transition, but with no significant effects on the driving. This region seems to have also a subtle effect on the trapping of low radial order {\em g}-modes (mixed modes with a majority of nodes in the G-mode region), but again with no effect on the driving. We found that for mode trapping to have an influence on the driving of sdO modes (1) the mode should be trapped in a way that the amplitude of the eigenfunctions is lower in a damping region and (2) in this damping region significant energy interchange has to be produced.

We have explored the possibility of driving pulsation modes in models of sdO stars in which the effects of element diffusion, gravitational settling and radiative levitation have been neglected so that the distribution of iron-peak elements remains uniform throughout the evolution. The stability of these models was determined using a non-adiabatic oscillations code. We analysed 27 sdO models from 16 different evolutionary sequences and discovered the first ever sdO models capable of driving high-radial order g-modes. In one model, the driving is by a classical kappa-mechanism due to the opacity bump from iron-peak elements at temperature ~200,000 K. In a second model, the driving result from the combined action of kappa-mechanisms operating in three distinct regions of the star: (i) a carbon-oxygen partial ionization zone at temperature ~2 10^6 K, (ii) a deeper region at temperature ~2 10^7 K, which we attribute to ionization of argon, and (iii) at the transition from radiative to conductive opacity in the core of the star.

A spectroscopic analysis of SDSS J160043.6+074802.9, a binary system containing a pulsating subdwarf-O (sdO) star with a late-type companion, yields Teff = 70 000 +/- 5000 K and log g = 5.25 +/- 0.30, together with a most likely type of K3V for the secondary star. We compare our results with atmospheric parameters derived by Fontaine et al. (2008) and in the context of existing evolution models for sdO stars. New and more extensive photometry is also presented which recovers most, but not all, frequencies found in an earlier paper. It therefore seems probable that some pulsation modes have variable amplitudes. A non-adiabatic pulsation analysis of uniform metallicity sdO models show those having log g > 5.3 to be more likely to be unstable and capable of driving pulsation in the observed frequency range.

Contrary to a common argument that a small increase in the strength of the strong force would lead to destruction of all hydrogen in the big bang due to binding of the diproton and the dineutron with a catastrophic impact on life as we know it, we show that provided the increase in strong force coupling constant is less than about 50% substantial amounts of hydrogen remain. The reason is that an increase in strong force strength leads to tighter binding of the deuteron, permitting nucleosynthesis to occur earlier in the big bang at higher temperature than in the standard big bang. Photodestruction of the less tightly bound diproton and dineutron delays their production to after the bulk of nucleosynthesis is complete. The decay of the diproton can, however, lead to relatively large abundances of deuterium.

We present the results of a multisite photometric campaign on the pulsating sdB star Balloon 090100001. The star is one of the two known hybrid hot subdwarfs with both long- and short-period oscillations. The campaign involved eight telescopes with three obtaining UBVR data, four B-band data, and one Stromgren uvby photometry. The campaign covered 48 nights, providing a temporal resolution of 0.36microHz with a detection threshold of about 0.2mmag in B-filter data. Balloon 090100001 has the richest pulsation spectrum of any known pulsating subdwarf B star and our analysis detected 114 frequencies including 97 independent and 17 combination ones. The strongest mode (f_1) in the 2.8mHz region is most likely radial while the remaining ones in this region form two nearly symmetric multiplets: a triplet and quintuplet, attributed to rotationally split \ell=1 and 2 modes, respectively. We find clear increases of splitting in both multiplets between the 2004 and 2005 observing campaigns, amounting to 15% on average. The observed splittings imply that the rotational rate in Bal09 depends on stellar latitude and is the fastest on the equator. We use a small grid of models to constrain the main mode (f_1), which most likely represents the radial fundamental pulsation. The groups of p-mode frequencies appear to lie in the vicinity of consecutive radial overtones, up to the third one. Despite the large number of g-mode frequencies observed, we failed to identify them, most likely because of the disruption of asymptotic behaviour by mode trapping. The observed frequencies were not, however, fully exploited in terms of seismic analysis which should be done in the future with a larger grid of reliable evolutionary models of hot subdwarfs.

We investigate whether angular momentum transport due to unstable pulsation modes can play a significant role in the rotational evolution of massive stars. We find that these modes can redistribute appreciable angular momentum, and moreover trigger shear-instability mixing in the molecular weight gradient zone adjacent to stellar cores, with significant evolutionary impact.

We present stellar evolution calculations for Population III stars for both single and binary star evolution. Our models include 10 Msun and 16.5 Msun single stars and a 10 Msun model star that undergoes an episode of accretion resulting in a final mass of 16.1 Msun. For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near solar heavy metal abundance stars and not including both the evolution and supernovae explosions in one work.

We examine the stability of l=1 and l=2 g modes in a pair of nitrogen-rich Wolf-Rayet stellar models characterized by differing hydrogen abundances. We find that modes with intermediate radial orders are destabilized by a kappa mechanism operating on an opacity bump at an envelope temperature log T ~ 6.25. This `deep opacity bump' is due primarily to L-shell bound-free transitions of iron. Periods of the unstable modes span ~ 11-21 hr in the model containing some hydrogen, and ~ 3-12 hr in the hydrogen-depleted model. Based on the latter finding, we suggest that self-excited g modes may be the source of the 9.8 hr-periodic variation of WR 123 recently reported by Lefevre et al. (2005).

We present a new set of stellar interior and synthesis models for predicting the integrated emission from stellar populations in star clusters and galaxies of arbitrary age and metallicity. This work differs from existing spectral synthesis codes in a number of important ways, namely (1) the incorporation of new stellar evolutionary tracks, with sufficient resolution in mass to sample rapid stages of stellar evolution; (2) a physically consistent treatment of evolution in the HR diagram, including the approach to the main sequence and the effects of mass loss on the giant and horizontal-branch phases. Unlike several existing models, ours yield consistent ages when used to date a coeval stellar population from a wide range of spectral features and colour indexes. We rigorously discuss degeneracies in the age-metallicity plane and show that inclusion of spectral features blueward of 4500 AA, suffices to break any remaining degeneracy and that with moderate S/N spectra (10 per 20AA, resolution element) age and metallicity are not degenerate. We also study sources of systematic errors in deriving the age of a single stellar population and conclude that they are not larger than 10-15%. We illustrate the use of single stellar populations by predicting the colors of primordial proto-galaxies and show that one can first find them and then deduce the form of the IMF for the early generation of stars in the universe. Finally, we provide accurate analytic fitting formulas for ultra fast computation of colors of single stellar populations. The models can be found at http://www.physics.upenn.edu/~raulj

We estimate the cosmic production rate of helium relative to metals ($\Delta Y/\Delta Z$) using K dwarf stars in the Hipparcos catalog with accurate spectroscopic metallicities. The best fitting value is $\Delta Y/\Delta Z=2.1 \pm 0.4$ at the 68% confidence level. Our derived value agrees with determinations from HII regions and with theoretical predictions from stellar yields with standard assumptions for the initial mass function. The amount of helium in stars determines how long they live and therefore how fast they will enrich the insterstellar medium with fresh material.

BVRI CCD-based and near-IR (J) imaging, together with unfiltered photometry of the hot subdwarf B star PG 0856+121 are reported. Two close, faint, red, point-like sources are resolved. They account for the previously reported IR excess observed in this hot subdwarf. In addition, the new unfiltered differential photometry of PG 0856+121 confirms its previously reported pulsational nature. A comparison with the oscillation modes of stellar models suggests the possible presence of g modes.

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed for following the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available ten years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al (1986). We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between $10^{-3}$ and $10^{-2} M_{\odot}$, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of $0.6 M_{\odot}$. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate.

We have applied a new method to analyze the horizontal branch (HB) morphology in relation to the distribution of stars near the red giant branch (RGB) tip for the globular clusters M22, M5, M68, M107, M72, M92, M3 and 47 Tuc. This new method permits determination of the cluster ages to greater accuracy than conventional isochrone fitting. Using the method in conjunction with our new high-quality photometric data for RGB and HB stars in the first five of these clusters, we discuss the origins of the spread in color on the HB and its relation to the `second parameter' problem. The oldest clusters in our sample are found to have relatively low ages ($13.5 \pm 2$ Gyr). A 1$\sigma$ uncertainty in each of the parameters of mass and helium content combined with the effects of helium diffusion gives a lower limit for the age of the oldest clusters of 9.7 Gyr.

We present here results of stellar evolution calculations that include the latest advances in radiative opacities and neutrino cooling, and discuss on the basis of these models how the internal stellar structure responds to mass-loss from the stellar surface. This problem has particular importance for the development of semi-analytical algorithms for effi cient calculation of synthetic stellar populations with realistic (and hence complex) mass-loss scenarios. We therefore compare our numerical results with test calculations based on a semi-analytical stellar evolution method developed by us. Although small, but important, differences between results from the two methods are revealed, the evolutionary tracks in the HR-diagram predicted by the two approaches are almost identical.