Researcher profile

J. MacDonald

J. MacDonald contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
15works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2019arXiv

Mass loss on the red giant branch: plasmoid-driven winds above the RGB bump

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.

preprint2016arXiv

Mechanical energy fluxes associated with saturated coronal heating in M dwarfs: comparison with predictions of a turbulent dynamo

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.

preprint2015arXiv

Magneto-convective models of red dwarfs: constraints imposed by the lithium abundance

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.

preprint2015arXiv

The effects of convection criteria on the evolution of population III stars and the detectability of their supernovae

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.

preprint2014arXiv

M dwarf search for pulsations within Kepler GO program

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íguez-López 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é et al. 2014, Berdiñas et al. 2014) with the high-resolution spectrographs HARPS and HARPS-N.

preprint2013arXiv

The extremely low metallicity star SDSS J102915+172927: a subgiant scenario

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.

preprint2013arXiv

The theoretical instability strip of M dwarf stars

The overstability of the fundamental radial mode in M dwarf models was theoretically predicted by Rodríguez-López 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.

preprint2012arXiv

An Empirical Measure of the Rate of White Dwarf Cooling in 47 Tucanae

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.

preprint2011arXiv

Pulsations in M dwarf stars

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 ε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.

preprint2010arXiv

Magneto-Convection and Lithium Age Estimates of the βPictoris Moving Group

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 β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 βPictoris moving group of 40 Myr is larger than previous determinations by a factor of about two.

preprint2009arXiv

Study of sdO models: mode trapping

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.

preprint2009arXiv

Study of sdO models. Pulsation Analysis

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.

preprint2009arXiv

The rapidly pulsating sdO star, SDSS J160043.6+074802.9

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.

preprint1996arXiv

Ages of globular clusters: a new approach

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&#39; problem. The oldest clusters in our sample are found to have relatively low ages ($13.5 \pm 2$ Gyr). A 1$σ$ 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.

preprint1995arXiv

Stellar evolution with mass loss - comparison of numerical and semi-analytical computations.

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.