Search

We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes., Comment: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica., Comment: Accepted for publication in PASA

Saturn's polar stratosphere exhibits the seasonal growth and dissipation of broad, warm, vortices poleward of $\sim75^\circ$ latitude, which are strongest in the summer and absent in winter. The longevity of the exploration of the Saturn system by Cassini allows the use of infrared spectroscopy to trace the formation of the North Polar Stratospheric Vortex (NPSV), a region of enhanced temperatures and elevated hydrocarbon abundances at millibar pressures. We constrain the timescales of stratospheric vortex formation and dissipation in both hemispheres. Although the NPSV formed during late northern spring, by the end of Cassini's reconnaissance (shortly after northern summer solstice), it still did not display the contrasts in temperature and composition that were evident at the south pole during southern summer. The newly-formed NPSV was bounded by a strengthening stratospheric thermal gradient near $78^\circ$N. The emergent boundary was hexagonal, suggesting that the Rossby wave responsible for Saturn's long-lived polar hexagon - which was previously expected to be trapped in the troposphere - can influence the stratospheric temperatures some 300 km above Saturn's clouds., Comment: 51 pages, 12 figures, published in Nature Communications

The evolution of star clusters is largely affected by the tidal field generated by the host galaxy. It is thus in principle expected that under the assumption of an "universal" initial cluster mass function the properties of the evolved present-day mass function of star cluster systems should show a dependency on the properties of the galactic environment in which they evolve. To explore this expectation a sophisticated model of the tidal field is required in order to study the evolution of star cluster systems in realistic galaxies. Along these lines, in the present work we first describe a method developed for coupling $N$-body simulations of galaxies and star clusters. We then generate a database of galaxy models along the Hubble sequence and calibrate evolutionary equations to the results of direct $N$-body simulations of star clusters in order to predict the clusters' mass evolution as function of the galactic environment. We finally apply our methods to explore the properties of evolved "universal" initial cluster mass functions and any dependence on the host galaxy morphology and mass distribution. The preliminary results show that an initial power-law distribution of the masses "universally" evolves into a log-normal distribution, with the properties correlated with the stellar mass and stellar mass density density of the host galaxy., Comment: 19 Pages, 16 Figures. Accepted for publication in MNRAS

We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour diagrams and through comparison to simple stellar population models infer that the GCs have old ages consistent with being 9 Gyr or older, while their metallicities are in the range between -1.6 < [Fe/H] < -0.4. We conduct a kinematic analysis of the GC population and find tentative evidence for weak net rotation of the GC system, in the same sense as that exhibited by the underlying spheroid. The most likely amplitude of rotation is ~10 km/s, approximately half the magnitude of the observed velocity dispersion. Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within 11 kpc and we formally find it to be in the range between (3-4)10^9 Msun. This implies an overall mass-to-light ratio in the range of ~ 30-40 and indicates that NGC 6822 is highly dark matter dominated. The mass and the corresponding mass-to-light ratio estimates are affected by various additional systematic effects due to limitations of the data and the model that are not necessary reflected in the formal uncertainties., Comment: 15 pages, 5 figures, 8 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

The seasonal evolution of Saturn's polar atmospheric temperatures and hydrocarbon composition is derived from a decade of Cassini Composite Infrared Spectrometer (CIRS) 7-16 $\mu$m thermal infrared spectroscopy. We construct a near-continuous record of atmospheric variability poleward of 60$^\circ$ from northern winter/southern summer (2004, $L_s=293^\circ$) through the equinox (2009, $L_s=0^\circ$) to northern spring/southern autumn (2014, $L_s=56^\circ$). The hot tropospheric polar cyclones and the hexagonal shape of the north polar belt are both persistent features throughout the decade of observations. The hexagon vertices rotated westward by $\approx30^\circ$ longitude between March 2007 and April 2013, confirming that they are not stationary in the Voyager-defined System III longitude system as previously thought. The extended region of south polar stratospheric emission has cooled dramatically poleward of the sharp temperature gradient near 75$^\circ$S, coinciding with a depletion in the abundances of acetylene and ethane, and suggestive of stratospheric upwelling with vertical wind speeds of $w\approx+0.1$ mm/s. This is mirrored by a general warming of the northern polar stratosphere and an enhancement in acetylene and ethane abundances that appears to be most intense poleward of 75$^\circ$N, suggesting subsidence at $w\approx-0.15$ mm/s. However, the sharp gradient in stratospheric emission expected to form near 75$^\circ$N by northern summer solstice (2017, $L_s=90^\circ$) has not yet been observed, so we continue to await the development of a northern summer stratospheric vortex. North polar minima in tropospheric and stratospheric temperatures were detected in 2008-2010 (lagging one season, or 6-8 years, behind winter solstice); south polar maxima appear to have occurred before the start of the Cassini observations (1-2 years after summer solstice). [Abridged], Comment: Preprint of article accepted for publication in Icarus

We propose an evolutionary model to describe the dynamical evolution of star cluster systems in tidal fields, in which we calibrated the parametric equations defining the model by running direct N-body simulations of star clusters with a wide range of initial masses and set of orbital parameters, living within the external tidal field generated by a disc-like galaxy. We derived a new method to solve numerically the evolutionary equations, allowing us to infer constraints on the mass of a star cluster from its age, present-day mass, orbital parameters and external gravitational potential. The result has been applied to the metal-rich subsample of Galactic globular clusters, being a good representation of a disc-bulge population. We reconstructed the initial mass function of these objects from the present-day mass function, finding that a lognormal distribution is well preserved during the evolution of the globular cluster system. The evolution of a power-law initial mass function has been evaluated, confirming that it transforms into a lognormal distribution of the cluster masses within an Hubble time. Our results are consistent with a formation scenario in which metal-rich Galactic globular clusters formed from giant molecular clouds in high-pressure regions during the early phases of the evolution of the Galactic disc and bulge., Comment: 16 pages, 22 figures, accepted for publication in MNRAS

Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397., Comment: Main Article: 10 pages, 4 figures; Supplementary Info 15 pages, 5 figures. Nature, Aug 1, 2013

A significant proportion of Milky Way stars are born in stellar clusters, which dissolve over time so that the members become part of the disc and halo populations of the Galaxy. In the present work we will assume that these young stellar clusters live mainly within the disc of the Galaxy and that they can have primordial binary percentages ranging from 0% to as high as 70%. We have evolved models of such clusters to an age of 4 Gyr through N-body simulations, paying attention to the stars and binaries that escape in the process. We have quantified the contribution of these escaping stars to the Galaxy population by analysing their escape velocity and evolutionary stage at the moment of escape. In this way we could analyse the mechanisms that produced these escapers, whether evaporation through weak two- body encounters, energetic close encounters or stellar evolution events, e.g. supernovae. In our models we found that the percentage of primordial binaries in a star cluster does not produce significant variations in the velocities of the stars that escape in the velocity range of 0-20 km/s. However, in the high-velocity 20-100 km/s range the number of escapers increased markedly as the primordial binary percentage increased. We could also infer that dissolving stellar clusters such as those that we have modelled can populate the Galactic halo with giant stars for which the progenitors were stars of up to 2.4 Msun. Furthermore, choices made for the velocity kicks of remnants do influence the production of hyper-velocity stars - and to a lesser extent stars in the high-velocity range - but once again the difference for the 99% of stars in the 0-20 km/s range is not significant., Comment: 24 pages, 20 figures, Accepted for publication in MNRAS

Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr., Comment: 25 pages, 20 figures, accepted for publication in Astrophysical Journal

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample., Comment: Accepted to The Astronomical Journal, 13 Figures, 2 Tables. Figures 3 and 6 are figure sets, each composed of 59 subfigures (to appear in the electronic journal). This is a Companion paper to the article ID: submit/0375614

We present a direct N-body simulation modeling the evolution of the old (7 Gyr) open cluster NGC 188. This is the first N-body open cluster simulation whose initial binary population is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare the simulated color-magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue stragglers produced in the simulation. We compare the solar-type main sequence binary period and eccentricity distributions of the simulation to detailed observations of similar binaries in NGC 188. These results demonstrate the importance of detailed observations in guiding N-body open cluster simulations. Finally, we discuss the implications of a few discrepancies between the NGC 188 model and observations and suggest a few methods for bringing N-body open cluster simulations into better agreement with observations., Comment: Contributed talk at IAUS 266 'Star clusters: basic galactic building blocks', R. de Grijs and J.R.D. Lepine, eds. 6 pages, 4 figures

We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters using direct N-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which the star cluster loses mass and we find that increasing halo mass and concentration drives enhanced mass loss rates and in principle shorter cluster disruption timescales. In addition, we examine disruption timescales in a three-component model of a galaxy (bulge, disk and dark matter halo) and find good agreement with results based on an empirical model of the Galactic potential if we assume a halo mass of ~1e12 solar masses. In general, dark matter halos are expected to contribute significantly to the masses of galaxies and should not be ignored when modelling the evolution of star clusters. We extend our results to discuss how this can have a potentially profound effect on the disruption timescales of globular clusters, suggesting that we may underestimate the rate at which primordial globular clusters are disrupted., Comment: 8 pages, 4 figures, accepted for publication in New Astronomy

Stellar collisions are an important formation channel for blue straggler stars in globular and old open clusters. Hydrodynamical simulations have shown that the remnants of such collisions are out of thermal equilibrium, are not strongly mixed and can rotate very rapidly. Detailed evolution models of collision products are needed to interpret observed blue straggler populations and to use them to probe the dynamical history of a star cluster. We expand on previous studies by presenting an efficient procedure to import the results of detailed collision simulations into a fully implicit stellar evolution code. Our code is able to evolve stellar collision products in a fairly robust manner and allows for a systematic study of their evolution. Using our code we have constructed detailed models of the collisional blue stragglers produced in the $N$-body simulation of M67 performed by Hurley \emph{et al.} in 2005. We assume the collisions are head-on and thus ignore the effects of rotation in this paper. Our detailed models are more luminous than normal stars of the same mass and in the same stage of evolution, but cooler than homogeneously mixed versions of the collision products. The increased luminosity and inefficient mixing decrease the remaining main-sequence lifetimes of the collision products, which are much shorter than predicted by the simple prescription commonly used in $N$-body simulations., Comment: To be published in A&A

Using Hubble Space Telescope (HST) observations of the globular cluster NGC6397, we constrain the cluster's binary fraction. The observations consist of two fields: the primary science field, a single ACS pointing centered approximately 5' from the cluster center; and the parallel field, a single WFPC2 field centered on the cluster center. Using the exquisite photometric precision of these observations, we determine the binary fraction in these regions of the cluster by examining stars lying off the main sequence. The binary fraction is constrained to be 0.012 +/- 0.004 in the ACS field, and to be 0.051 +/- 0.010 in the WFPC field. N-body simulations by Hurley et al. (2007) suggest that the binary fraction remains nearly constant beyond the half-mass radius for the lifetime of the cluster. In the context of these simulations, our results suggest that NGC 6397 had a primordial binary fraction of only ~1%., Comment: 19 pages, 7 figures, accepted in AJ

We have examined the radial distribution of white dwarfs over a single HST/ACS field in the nearby globular cluster NGC 6397. In relaxed populations, such as in a globular cluster, stellar velocity dispersion, and hence radial distribution, is directly dependent on stellar masses. The progenitors of very young cluster white dwarfs had a mass of ~0.8 solar masses, while the white dwarfs themselves have a mass of ~0.5 solar masses. We thus expect young white dwarfs to have a concentrated radial distribution (like that of their progenitors) that becomes more extended over several relaxation times to mimic that of ~0.5 solar mass main-sequence stars. However, we observe young white dwarfs to have a significantly extended radial distribution compared to both the most massive main sequence stars in the cluster and also to old white dwarfs., Comment: 13 pages including 1 table and 3 figures. Accepted for publication in the MNRAS Letters

Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems., Comment: 3 pages, no figs 1 table. Published in proceedings of Active OB stars: Laboratories for Stellar and Circumstellar Physics (Aug 29th to Sep 2, 2005; Hokkaido Uni, Sapporo, Japan)

This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified., Comment: Accepted by New Astronomy

We present colour-magnitude diagrams for two rich (~10^4 Msun) Large Magellanic Cloud star clusters with ages ~10^7 years, constructed from optical and near-infrared data obtained with the Hubble Space Telescope. These data are part of an HST project to study LMC clusters with a range of ages. In this paper we investigate the massive star content of the young clusters, and determine the cluster ages and metallicities, paying particular attention to Be star and blue straggler populations and evidence of age spreads. We compare our data to detailed stellar population simulations to investigate the turn-off structure of ~25 Myr stellar systems, highlighting the complexity of the blue straggler phenomenon., Comment: 15 pages, accepted for publication in MNRAS, postscript version with higher quality figures available from http://www.ast.cam.ac.uk/~raj/biblio.html

(shortened) Direct N-body calculations are presented of the formation of Galactic clusters using GasEx, which is a variant of the code Nbody6. The calculations focus on the possible evolution of the Orion Nebula Cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10^4 stars, and resembles the Pleiades Cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions (DIVA and GAIA). These missions should also detect the original Pleiades members as an associated expanding young Galactic-field sub-population. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations., Comment: MNRAS, in press, 36 pages, 15 figures; repl.vers. contains adjustments for consistency with published version

We present analytic formulae that approximate the evolution of stars for a wide range of mass and metallicity. Stellar luminosity, radius and core mass are given as a function of age, M and Z, for all phases from the zero-age main-sequence up to, and including, the remnant stages. For the most part we find continuous formulae accurate to within 5% of detailed models. These formulae are useful for purposes such as population synthesis that require very rapid but accurate evaluation of stellar properties, and in particular for use in combination with N-body codes. We describe a mass loss prescription that can be used with these formulae and investigate the resulting stellar remnant distribution., Comment: 29 pages, 20 figures, submitted for publication in MNRAS

This report examines the rural principal's role by investigating the attitudes of rural school teachers who have "principal potential." Interviews were carried out with 25 teachers from 5 rural school districts who were identified by their peers as having school leadership potential. Respondents were asked to reflect on advantages and disadvantages associated with the principal's role. Only five of the teachers interviewed expressed interest in becoming principals. Others either were not interested or stated that the principal's role would have to change significantly before they would consider it. They viewed the role of the principal as: (1) too distant from the instructional core; (2) involving too many non-instructional duties; (3) too large to perform duties effectively; (4) requiring too much time away from family; (5) providing too little direct contact with students; (6) requiring behaviors that do not match personality; (7) affecting existing relationships with teachers; (8) lacking autonomy and control; (9) seldom offered to a woman; (10) disciplinarian, affecting relationships with students; (11) too political; and (12) needing redefinition. To attract teacher-leaders into the rural school principalship, the principal's role should focus more on instructional leadership and activities that directly affect students. School restructuring suggestions should be carried out so that more educational decisions are made at the school level, school cultures become more collegial, and more female teacher-leaders are encouraged to become principals. (LP)