Oscillations in a maturation model of blood cell production.

We present a mathematical model of blood cell production which describes both the development of cells through the cell cycle, and the maturation of these cells as they differentiate to form the various mature blood cell types. The model differs from earlier similar ones by considering primitive stem cells as a separate population from the differentiating cells, and this formulation removes an apparent inconsistency in these earlier models. Three different controls are included in the model: proliferative control of stem cells, proliferative control of differentiating cells, and peripheral control of stem cell committal rate. It is shown that an increase in sensitivity of these controls can cause oscillations to occur through their interaction with time delays associated with proliferation and differentiation, respectively. We show that the characters of these oscillations are quite distinct and suggest that the model may explain an apparent superposition of fast and slow oscillations which can occur in cyclical neutropenia. © 2006 Society for Industrial and Applied Mathematics

Asymmetric supernova remnants generated by Galactic, massive runaway stars

After the death of a runaway massive star, its supernova shock wave interacts with the bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its spherical symmetry before expanding into the local interstellar medium (ISM). We investigate whether the initial mass and space velocity of these progenitors can be associated with asymmetric supernova remnants. We run hydrodynamical models of supernovae exploding in the pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks that accumulate more than about 1.5 Mo generate asymmetric remnants. The shock wave first collides with these bow shocks 160-750 yr after the supernova, and the collision lasts until 830-4900 yr. The shock wave is then located 1.35-5 pc from the center of the explosion, and it expands freely into the ISM, whereas in the opposite direction it is channelled into the region of undisturbed wind material. This applies to an initially 20 Mo progenitor moving with velocity 20 km/s and to our initially 40 Mo progenitor. These remnants generate mixing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from the center of the explosion. Their lightcurves are dominated by emission from optically-thin cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely to be observed in the [OIII] lambda 5007 spectral line emission or in the soft energy-band of X-rays. Finally, we discuss our results in the context of observed Galactic supernova remnants such as 3C391 and the Cygnus Loop.Comment: 21 pages, 16 figure

On the Exponentials of Some Structured Matrices

In this note explicit algorithms for calculating the exponentials of important structured 4 x 4 matrices are provided. These lead to closed form formulae for these exponentials. The techniques rely on one particular Clifford Algebra isomorphism and basic Lie theory. When used in conjunction with structure preserving similarities, such as Givens rotations, these techniques extend to dimensions bigger than four.Comment: 19 page

Models of the circumstellar medium of evolving, massive runaway stars moving through the Galactic plane

At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass loss and space velocity of massive runaway stars affect the morphology of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [OIII]. The Ha emission of the bow shocks around hot stars originates from near their contact discontinuity. The H$\alpha$ emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically-thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.Comment: 22 pages, 24 figure

On the stability of bow shocks generated by red supergiants: the case of IRC-10414

In this Letter, we explore the hypothesis that the smooth appearance of bow shocks around some red supergiants (RSGs) might be caused by the ionization of their winds by external sources of radiation. Our numerical simulations of the bow shock generated by IRC-10414 (the first-ever RSG with an optically detected bow shock) show that the ionization of the wind results in its acceleration by a factor of two, which reduces the difference between the wind and space velocities of the star and makes the contact discontinuity of the bow shock stable for a range of stellar space velocities and mass-loss rates. Our best fit model reproduces the overall shape and surface brightness of the observed bow shock and suggests that the space velocity and mass-loss rate of IRC-10414 are $\approx$50 ${\rm km} \, {\rm s}^{-1}$ and $\approx$$10^{-6}$ $M_\odot \, {\rm yr}^{-1}$, respectively, and that the number density of the local ISM is $\approx$3 ${\rm cm}^{-3}$. It also shows that the bow shock emission comes mainly from the shocked stellar wind. This naturally explains the enhanced nitrogen abundance in the line-emitting material, derived from the spectroscopy of the bow shock. We found that photoionized bow shocks are $\approx$15$-$50 times brighter in optical line emission than their neutral counterparts, from which we conclude that the bow shock of IRC-10414 must be photoionized.Comment: 5 pages, 5 figures. Accepted for publication in MNRAS Letter

Young accreted globular clusters in the outer halo of M31

We report on Gemini/GMOS observations of two newly discovered globular clusters in the outskirts of M31. These objects, PAndAS-7 and PAndAS-8, lie at a galactocentric radius of ~87 kpc and are projected, with separation ~19 kpc, onto a field halo substructure known as the South-West Cloud. We measure radial velocities for the two clusters which confirm that they are almost certainly physically associated with this feature. Colour-magnitude diagrams reveal strikingly short, exclusively red horizontal branches in both PA-7 and PA-8; both also have photometric [Fe/H] = -1.35 +/- 0.15. At this metallicity, the morphology of the horizontal branch is maximally sensitive to age, and we use the distinctive configurations seen in PA-7 and PA-8 to demonstrate that both objects are very likely to be at least 2 Gyr younger than the oldest Milky Way globular clusters. Our observations provide strong evidence for young globular clusters being accreted into the remote outer regions of M31 in a manner entirely consistent with the established picture for the Milky Way, and add credence to the idea that similar processes play a central role in determining the composition of globular cluster systems in large spiral galaxies in general.Comment: 14 pages, 8 figures, accepted for publication in MNRA

Deterministic Brownian motion generated from differential delay equations

This paper addresses the question of how Brownian-like motion can arise from the solution of a deterministic differential delay equation. To study this we analytically study the bifurcation properties of an apparently simple differential delay equation and then numerically investigate the probabilistic properties of chaotic solutions of the same equation. Our results show that solutions of the deterministic equation with randomly selected initial conditions display a Gaussian-like density for long time, but the densities are supported on an interval of finite measure. Using these chaotic solutions as velocities, we are able to produce Brownian-like motions, which show statistical properties akin to those of a classical Brownian motion over both short and long time scales. Several conjectures are formulated for the probabilistic properties of the solution of the differential delay equation. Numerical studies suggest that these conjectures could be "universal" for similar types of "chaotic" dynamics, but we have been unable to prove this.Comment: 15 pages, 13 figure

Chemical evolution of star clusters

I discuss the chemical evolution of star clusters, with emphasis on old globular clusters, in relation to their formation histories. Globular clusters clearly formed in a complex fashion, under markedly different conditions from any younger clusters presently known. Those special conditions must be linked to the early formation epoch of the Galaxy and must not have occurred since. While a link to the formation of globular clusters in dwarf galaxies has been suggested, present-day dwarf galaxies are not representative of the gravitational potential wells within which the globular clusters formed. Instead, a formation deep within the proto-Galaxy or within dark-matter minihaloes might be favoured. Not all globular clusters may have formed and evolved similarly. In particular, we may need to distinguish Galactic halo from Galactic bulge clusters.Comment: 27 pages, 2 figures. To appear as invited review article in a special issue of the Phil. Trans. Royal Soc. A: Ch. 6 "Star clusters as tracers of galactic star-formation histories" (ed. R. de Grijs). Fully peer reviewed. LaTeX, requires rspublic.cls style fil

The adrenocortical response to synthetic ACTH following a trek to high altitude

Background: Gradual ascent to high altitude is typically associated with reduced resting aldosterone and unchanged cortisol, features that may facilitate acclimatization but are poorly understood. Aims: To investigate the cortisol and aldosterone response to adrenocorticotrophic hormone at altitude. Methods: Eleven subjects underwent a 250 µg short synacthen test at sea-level and again after trekking to 3600 m in Nepal. Cortisol and aldosterone were measured by conventional assay from blood samples taken immediately prior to the administration of synacthen (T0) and then 30 (T30) and 60 (T60) minutes later. Results: At 3600 m resting basal cortisol and aldosterone levels were both significantly lower than they were at sea-level (p=0.004, p=0.003 respectively). Cortisol values at T30 and T60 were no different between sea-level and 3600 m but the increment after synacthen was significantly (p=0.041) greater at 3600 m due to a lower basal value. Aldosterone at T30 and T60 was significantly lower (p=0.003 for both) at 3600 m than at sea-level and the increment following synacthen was also significantly (p=0.003) less at 3600 m. Conclusions: At 3600 m there appears to be a divergent adrenal response to synthetic adrenocorticotrophic hormone with an intact cortisol response but a reduced aldosterone response, relative to sea-level. This may reflect a specific effect of hypoxia on aldosterone synthesis and may be beneficial to acclimatization