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.

In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9-100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models.

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $\chi_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.

The past four years have seen a scientific revolution through the birth of a new field: gravitational-wave astronomy. The first detection of gravitational waves---recognised by the 2017 Nobel Prize in Physics---provided unprecedented tests of general relativity while unveiling a previously unknown class of massive black holes, thirty times more massive than the Sun. The subsequent detection of gravitational waves from a merging binary neutron star confirmed the hypothesised connection between binary neutron stars and short gamma-ray bursts while providing an independent measurement of the expansion of the Universe. The discovery enabled precision measurement of the speed of gravity while shedding light on the origin of heavy elements. At the time of writing, the Laser Interferometer Gravitational-wave Observatory (LIGO) and its European partner, Virgo, have published the detection of eleven gravitational-wave events. New, not-yet-published detections are announced on a nearly weekly basis. This fast-growing catalogue of gravitational-wave transients is expected to yield insights into a number of topics, from the equation of state of matter at supra-nuclear densities to the fate of massive stars. The science potential of 3G observatories is enormous, enabling measurements of gravitational waves from the edge of the Universe and precise determination of the neutron star equation of state. Australia is well-positioned to help develop the required technology. The Mid-term Review for the Decadal plan for Australian astronomy 2016-2025 should consider investment in a scoping study for an Australian Gravitational-Wave Pathfinder that develops and validates core technologies required for the global 3G detector network.

We compare the results of a large grid of N-body simulations with the surface brightness and velocity dispersion profiles of the globular clusters $\omega$ Cen and NGC 6624. Our models include clusters with varying stellar-mass black hole retention fractions and varying masses of a central intermediate-mass black hole (IMBH). We find that an $\sim 45,000$ M$_\odot$ IMBH, whose presence has been suggested based on the measured velocity dispersion profile of $\omega$ Cen, predicts the existence of about 20 fast-moving, $m>0.5$ M$_\odot$ main-sequence stars with a (1D) velocity $v>60$ km/sec in the central 20 arcsec of $\omega$ Cen. However no such star is present in the HST/ACS proper motion catalogue of Bellini et al. (2017), strongly ruling out the presence of a massive IMBH in the core of $\omega$ Cen. Instead, we find that all available data can be fitted by a model that contains 4.6% of the mass of $\omega$ Cen in a centrally concentrated cluster of stellar-mass black holes. We show that this mass fraction in stellar-mass BHs is compatible with the predictions of stellar evolution models of massive stars. We also compare our grid of $N$-body simulations with NGC 6624, a cluster recently claimed to harbor a 20,000 M$_\odot$ black hole based on timing observations of millisecond pulsars. However, we find that models with $M_{IMBH}>1,000$ M$_\odot$ IMBHs are incompatible with the observed velocity dispersion and surface brightness profile of NGC 6624,ruling out the presence of a massive IMBH in this cluster. Models without an IMBH provide again an excellent fit to NGC 6624.

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.

The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $\rho \lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated version of the Binary Stellar Evolution (BSE) code for single and binary star evolution. We find that the recently suggested empirical prescription of consequential angular momentum loss (CAML) brings into agreement predicted and observed space densities of CVs and period bouncers. To progress with our understanding of CV evolution it is crucial to understand the physical mechanism behind empirical CAML. Our changes to the BSE code are also provided in details, which will allow the community to accurately model mass transfer in interacting binaries in which degenerate objects accrete from low-mass main-sequence donor stars.

The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number found in earlier analytical work and commonly used in the literature. Moreover, the evaporation rate derived with these simulations is higher than values previously published. These models also show that the value of the mass evaporation rate depends on the strength of the tidal field. We give an analytical estimate of the mass evaporation rate as a function of time and galactocentric distance xi(R_{GC},t). Upon extrapolating this formula to smaller R_{GC} values, our results provide tentative evidence for a very high xi value at small R_{GC}. Our results suggest that the corresponding mass loss in the inner Galactic potential could be high and it should be accounted for when star clusters pass within it. This has direct relevance to nuclear cluster formation/growth via the infall of globular clusters through dynamical friction. As an illustrative example, we estimate how the evaporation rate increases for a ~10^5 solar masses globular cluster that decays through dynamical friction into the galactic centre. We discuss the findings of this work in relation to the formation of nuclear star clusters by inspiraling globular clusters.

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.

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

Explaining the origin and evolution of exoplanetary "hot Jupiters" remains a significant challenge. One possible mechanism for their production is planet-planet interactions, which produces hot Jupiters from planets born far from their host stars but near their dynamical stability limits. In the much more likely case of planets born far from their dynamical stability limits, can hot Jupiters can be formed in star clusters? Our N-body simulations of planetary systems inside star clusters answer this question in the affirmative, and show that hot Jupiter formation is not a rare event. We detail three case studies of the dynamics-induced births of hot Jupiters on highly eccentric orbits that can only occur inside star clusters. The hot Jupiters' orbits bear remarkable similarities to those of some of the most extreme exoplanets known: HAT-P-32 b, HAT-P-2 b, HD 80606 b and GJ 876 d. If stellar perturbations formed these hot Jupiters then our simulations predict that these very hot, inner planets are often accompanied by much more distant gas giants in highly eccentric orbits.

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.

The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our attention to the evolution of star clusters. We have carried out our simulations in two steps: (1) we create a simulation of dissolution and mixing processes which yields a close fit to the present-day Milky Way dynamics and (2) we have evolved three sets of stellar clusters with masses of 400, 1000 and 15000 M$_{\odot}$ to dissolution. The birth location of these sets was 4, 6, 8 and 10 kpc for the 400 and 1000 M$_{\odot}$ clusters and 4, 6, 8, 10 and 12 kpc for the 15000 M$_{\odot}$. We have focused our efforts on studying the state of the escapers from these clusters after 4.5 Gyr of evolution with particular attention to stars that reach the Solar annulus, i.e. 7.5 $\le$ R$_{\rm{gc}}$ $\le$ 8.5 kpc. We give results for Solar twins and siblings over a wide range of radii and cluster masses for two dissolution mechanisms. From kinematics alone, we conclude that the Sun was $\sim$50 per cent more likely to have been born near its current Galactocentric radius, rather than have migrated (radially) $\sim$2 kpc since birth. We conclude our analysis by calculating magnitudes and colours of our single stars for comparison with the samples that the $Gaia$, $Gaia$-ESO and GALAH-AAO surveys will obtain. In terms of reconstructing dissolved star clusters we find that on short time-scales we cannot rely on kinematic evolution alone and thus it will be necessary to extend our study to include information on chemical space.

The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3 - 2.2 MSun) in NGC 1818 have a binary frequency that decreases towards the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in the core). We show here, through detailed N-body modeling, that both clusters could have formed with the same primordial binary frequency and with binary orbital elements and masses drawn from the same distributions (defined from observations of open clusters and the field of our Galaxy). The observed radial trends in binary frequency for both clusters are best matched with models that have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency that rises only towards the core, as is observed in old Milky Way star clusters). This recovery rate increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar dynamical sequence in the radial distribution of the binary frequency as a function of binary primary mass.

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]

A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently incorporated into a degenerate core. We synthesize a population of binary systems to investigate the stellar merging hypothesis for observed single HFMWDs. Our calculations provide mass distribution and the fractions of white dwarfs that merge during a common envelope phase or as double degenerate systems in a post common envelope phase. We vary the common envelope efficiency parameter alpha and compare with observations. We find that this hypothesis can explain both the observed incidence of magnetism and the mass distribution of HFMWDs for a wide range of alpha. In this model, the majority of the HFMWDs are of the Carbon Oxygen type and merge within a common envelope. Less than about a quarter of a per cent of HFMWDs originate from double degenerate stars that merge after common envelope evolution and these populate the high-mass tail of the HFMWD mass distribution.

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.

We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

Direct N-body simulations of globular clusters in a realistic Milky Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relationship between disk mass and star cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5x10^10 solar masses to 10x10^10 solar masses halves the dissolution time of a satellite star cluster orbiting the host galaxy at 6 kpc from the galactic center. Different geometries in a disk of identical mass can determine either the survival or dissolution of a star cluster orbiting within the inner 6 kpc of the galactic center. Furthermore, disk geometry has measurable effects on the mass loss of star clusters up to 15 kpc from the galactic center. N-body simulations performed with a fine output time step show that at each disk crossing the outer layers of star clusters experience an increase in velocity dispersion of ~5% of the average velocity dispersion in the outer section of star clusters. This leads to an enhancement of mass-loss -- a clearly discernable effect of disk shocking. By running models with different inclinations we determine that star clusters with an orbit perpendicular to the Galactic plane have larger mass loss rates than both clusters evolving in the Galactic plane or in an inclined orbit.

We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries.

The massive (13,000-26,000 Msolar), young (15-30 Myr) Large Magellanic Cloud star cluster NGC 1818 reveals an unexpected increasing binary frequency with radius for F-type stars (1.3-2.2 Msolar). This is in contrast to many older star clusters that show a decreasing binary frequency with radius. We study this phenomenon with sophisticated N-body modeling, exploring a range of initial conditions, from smooth virialized density distributions to highly substructured and collapsing configurations. We find that many of these models can reproduce the cluster's observed properties, although with a modest preference for substructured initial conditions. Our models produce the observed radial trend in binary frequency through disruption of soft binaries (with semi-major axes, a > 3000 AU), on approximately a crossing time (~ 5.4 Myr), preferentially in the cluster core. Mass segregation subsequently causes the binaries to sink towards the core. After roughly one initial half-mass relaxation time (t_rh(0) ~ 340 Myr) the radial binary frequency distribution becomes bimodal, the innermost binaries having already segregated towards the core, leaving a minimum in the radial binary frequency distribution that marches outwards with time. After 4-6 t_rh(0), the rising distribution in the halo disappears, leaving a radial distribution that rises only towards the core. Thus, both a radial binary frequency distribution that falls towards the core (as observed for NGC 1818) and one that rises towards the core (as for older star clusters) can arise naturally from the same evolutionary sequence owing to binary disruption and mass segregation in rich star clusters.

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

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.

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.

We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.

We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits. The stellar main sequence of the SMC is measured in the F606W, F814W color-magnitude diagram (CMD) down to ~30th magnitude, and is cleanly separated from the foreground star cluster population using proper motions. We simulate the SMC population by extracting stellar masses (single and unresolved binaries) from specific IMFs, and converting those masses to luminosities in our bandpasses. The corresponding photometry for these simulated stars is drawn directly from a rich cloud of 4 million artificial stars, thereby accounting for the real photometric scatter and completeness of the data. Over a continuous and well populated mass range of M = 0.37 - 0.93 Msun (i.e., down to a ~75% completeness limit at F606W = 28.7), we demonstrate that the IMF is well represented by a single power-law form with slope \alpha = -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^{\alpha}). This is shallower than the Salpeter slope of \alpha = -2.35, which agrees with the observed stellar luminosity function at higher masses. Our results indicate that the IMF does {\it not} turn over to a more shallow power-law form within this mass range. We discuss implications of this result for the theory of star formation, the inferred masses of galaxies, and the (lack of a) variation of the IMF with metallicity.

While tens or hundreds of stellar-remnant black holes are expected to form in globular star clusters, it is still unclear how many of those will be retained upon formation, and how many will be ejected through subsequent dynamical interactions. No such black holes have been found in any Milky Way globular cluster until the recent discovery of stellar-mass black holes in the globular cluster M22 (NGC 6656) with now an estimated population of 5-100 black holes. We present a direct N-body model of a star cluster of the same absolute and dynamical age as M22. Imposing an initial retention fraction of approx. 10% for black holes, 16 stellar-remnant black holes are retained at a cluster age of 12 Gyr, in agreement with the estimate for M22. Of those 16 BHs, two are in a binary system with a main sequence star each while also one pure black hole binary is present. We argue that multiple black holes can be present in any Milky Way cluster with an extended core radius, such as M22 or the model presented here.

(Abridged) Following on from a recently completed radial-velocity survey of the old (7 Gyr) open cluster NGC 188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we investigate the dynamical evolution of NGC 188 through a sophisticated N-body model. We employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, to guide our choices for parameters of the initial binary population. At 7 Gyr the main-sequence solar-type hard-binary population in the model matches that of NGC 188 in both binary frequency and distributions of orbital parameters. This agreement between the model and observations is in a large part due to the similarities between the NGC 188 and M35 solar-type binaries. Indeed, among the 7 Gyr main-sequence binaries in the model, only those with P>1000 days show potentially observable evidence for modifications by dynamical encounters. This emphasizes the importance of defining accurate initial conditions for star cluster models, which we propose is best accomplished through comparisons with observations of young open clusters like M35. Furthermore, this suggests that observations of the present-day binaries in even old open clusters can provide valuable information on their primordial binary populations. However, despite the model's successes at matching the true cluster, the model underproduces blue stragglers and produces an overabundance of long-period circular main-sequence--white-dwarf binaries as compared to NGC 188. We conclude that improvements in the physics of mass transfer and common envelope may in fact solve both discrepancies with the observations. This project highlights the unique accessibility of open clusters to both comprehensive observational surveys and full-scale N-body simulations, and underscores the importance of open clusters to the study of star cluster dynamics.

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.

Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolution. The models start with N = 100000 stars and include primordial binaries. We also take metallicity dependent stellar evolution and an external tidal field into account. We find no significant difference for the half-mass radii of those models, indicating that the clusters are structurally similar. However, utilizing observational tools to fit half-light (or effective) radii confirms that metallicity effects related to stellar evolution combined with dynamical effects such as mass segregation produce an apparent size difference of 17% on average. The metallicity effect on the overall cluster luminosity also leads to higher mass-to-light ratios for metal-rich clusters.

We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80% of its stars in the process. It reaches the end of the main core-collapse phase at 10.5 Gyr and experiences core oscillations from that point onwards -- direct numerical confirmation of this phenomenon. However, we find that after a further 1 Gyr the core oscillations are halted by the ejection of a massive binary comprised of two black holes from the core, producing a core that shows no signature of the prior core-collapse. We also show that the results of previous studies with N ranging from 500 to 100 000 scale well to this new model with larger N. In particular, the timescale to core-collapse (in units of the relaxation timescale), mass segregation, velocity dispersion, and the energies of the binary population all show similar behaviour at different N.

Direct N-body simulations of star clusters in a realistic Milky Way-like potential are carried out using the code NBODY6. Based on these simulations a new relationship between scale size and galactocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent of the galactocentric distance. The half-mass radius of star clusters increases systematically with galactocentric distance but levels off when star clusters orbit the galaxy beyond ~40 kpc. These simulations show that the half-mass radius of individual star clusters varies significantly as they evolve over a Hubble time, more so for clusters with shorter relaxation times, and remains constant through several relaxation times only in certain situations when expansion driven by the internal dynamics of the star cluster and the influence of the host galaxy tidal field balance each other. Indeed, the radius of a star cluster evolving within the inner 20 kpc of a realistic galactic gravitational potential is severely truncated by tidal interactions and does not remain constant over a Hubble time. Furthermore, the half-mass radius of star clusters measured with present day observations bears no memory of the original cluster size. Stellar evolution and tidal stripping are the two competing physical mechanisms that determine the present day size of globular clusters. These simulations also show that extended star clusters can form at large galactocentric distances while remaining fully bound to the host galaxy. There is thus no need to invoke accretion from an external galaxy to explain the presence of extended clusters at large galactocentric distances in a Milky Way-type galaxy.

In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.

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.

We use direct N-body simulations to investigate the evolution of star clusters with large size-scales with the particular goal of understanding the so-called extended clusters observed in various Local Group galaxies, including M31 and NGC6822. The N-body models incorporate a stellar mass function, stellar evolution and the tidal field of a host galaxy. We find that extended clusters can arise naturally within a weak tidal field provided that the tidal radius is filled at the start of the evolution. Differences in the initial tidal filling-factor can produce marked differences in the subsequent evolution of clusters and the size-scales that would be observed. These differences are more marked than any produced by internal evolution processes linked to the properties of cluster binary stars or the action of an intermediate-mass black hole, based on models performed in this work and previous work to date. Models evolved in a stronger tidal field show that extended clusters cannot form and evolve within the inner regions of a galaxy such as M31. Instead our results support the suggestion many extended clusters found in large galaxies were accreted as members of dwarf galaxies that were subsequently disrupted. Our results also enhance the recent suggestion that star clusters evolve to a common sequence in terms of their size and mass.

The work presented here examines populations of double compact binary systems and tidally enhanced collapsars. We make use of BINPOP and BINKIN, two components of a recently developed population synthesis package. Results focus on correlations of both binary and spatial evolutionary population characteristics. Pulsar and long duration gamma-ray burst observations are used in concert with our models to draw the conclusions that: double neutron star binaries can merge rapidly on timescales of a few million years (much less than that found for the observed double neutron star population), common envelope evolution within these models is a very important phase in double neutron star formation, and observations of long gamma-ray burst projected distances are more centrally concentrated than our simulated coalescing double neutron star and collapsar Galactic populations. Better agreement is found with dwarf galaxy models although the outcome is strongly linked to the assumed birth radial distribution. The birth rate of the double neutron star population in our models range from 4-160 Myr^-1 and the merger rate ranges from 3-150 Myr^-1. The upper and lower limits of the rates results from including electron capture supernova kicks to neutron stars and decreasing the common envelope efficiency respectively. Our double black hole merger rates suggest that black holes should receive an asymmetric kick at birth.

We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.

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.

We study the dynamical evolution of the young star cluster Arches and its dependence on the assumed initial stellar mass function (IMF). We perform many direct $N$-body simulations with various initial conditions and two different choices of IMFs. One is a standard Kroupa IMF without any mass segregation. The other is a radially dependent IMF, as presently observed in the Arches. We find that it is unlikely for the Arches to have attained the observed degree of mass segregation at its current age starting from a standard non-segregated Kroupa IMF. We also study the possibility of a collisional runaway developing in the first $\sim 2-3 \rm{Myr}$ of dynamical evolution. We find that the evolution of this cluster is dramatically different depending on the choice of IMF: if a primordially mass segregated IMF is chosen, a collisional runaway should always occur between $2-3 \rm{Myr}$ for a broad range of initial concentrations. In contrast, for a standard Kroupa IMF no collisional runaway is predicted. We argue that if Arches was created with a mass segregated IMF similar to what is observed today then at the current cluster age a very unusual, high-mass star should be created. However, whether a collisional runaway leads to the formation of an intermediate-mass black hole (IMBH) depends strongly on the mass loss rate via winds from massive stars. Growth of stellar mass through collisions can be quenched by strong wind mass loss. In that case, the inter-cluster as well as intra-cluster medium are expected to have a significant Helium enrichment which may be observed via Helium recombination lines. The excess amount of gas lost in winds may also be observed via X-ray observations as diffused X-ray sources.

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.

We produce synthetic populations of pulsars within our Galaxy and calculate the resulting scale heights as well as the radial and space velocity distributions of the pulsars. Results are presented for isolated pulsars, binary pulsars and millisecond pulsars. We also test the robustness of the outcomes to variations in the assumed form of the Galactic potential, the birth distribution of binary positions, and the strength of the velocity kick given to neutron stars at birth. We find that isolated pulsars have a greater scale height than binary pulsars. This is also true when restricted to millisecond pulsars unless we allow for low-mass stars to be ablated by radiation from their pulsar companion in which case the isolated and binary scale heights are comparable. Double neutron stars are found to have a large variety of space velocities, in particular, some systems have speeds similar to the Sun. We look in detail at the predicted Galactic population of millisecond pulsars with black hole companions, including their formation pathways, and show where the short-period systems reside in the Galaxy. Some of our population predictions are compared in a limited way to observations but the full potential of this aspect will be realised in the near future when we complete our population synthesis code with the selection effects component.

Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 suggest that WDs receive a kick of a few km/s shortly before they are born. Using our Monte Carlo cluster evolution code, which includes accurate treatments of all relevant physical processes operating in globular clusters, we study the effects of the kicks on their host cluster and on the WD population itself. We find that in clusters whose velocity dispersion is comparable to the kick speed, WD kicks are a significant energy source for the cluster, prolonging the initial cluster core contraction phase significantly so that at late times the cluster core to half-mass radius ratio is a factor of up to ~ 10 larger than in the no-kick case. WD kicks thus represent a possible resolution of the large discrepancy between observed and theoretically predicted values of this key structural parameter. Our modeling also reproduces the observed trend for younger WDs to be more extended in their radial distribution in the cluster than older WDs.

An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic stellar evolution code SSE. We find that, providing the physics package within STARS is set up correctly to match the parameters of the models used to construct SSE, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in SSE would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.

We present the first results of a project aimed at following the formation and long-term dynamical evolution of star clusters within the potential of a host galaxy. Here we focus on a model evolved within a simplified potential representing the Large Magellanic Cloud. This demonstrates for the first time the self-consistent formation of a bound star cluster from a giant molecular cloud. The model cluster reproduces the density profiles and structural characteristics of observed star clusters.

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.

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad-hoc manner. Thus no complete and direct comparison with observations of the pulsar population within the Galactic disk has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed \textit{typical} MSP formation process. Results also focus on reproducing the observed $P\dot{P}$ diagram for Galactic pulsars and how this precludes short timescales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.

We describe here the reduction methods that we developed to study the faintest red dwarfs and white dwarfs in an outer field of NGC6397, which was observed by \hst for 126 orbits in 2005. The particular challenge of this data set is that the faintest stars are not readily visible in individual exposures, so special care must be taken to combine the information in all the exposures in order to identify and measure them. Unfortunately, it is hard to find the faintest stars without also finding a large number of faint galaxies, so we developed specialized tools to distinguish between the point-like stars and the barely resolved galaxies. We found that artificial-star tests, while obviously necessary for completeness determination, can also play an important role in helping us optimize our finding and measuring algorithms. Although this paper focuses on this data set specifically, many of the techniques are new and might find application in other work, particularly when a large number of images is available for a single field.

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

We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.

We outline the steps needed in order to incorporate the evolution of single and binary stars into a particular Monte Carlo code for the dynamical evolution of a star cluster. We calibrate the results against N-body simulations, and present models for the evolution of the old open cluster M67 (which has been studied thoroughly in the literature with N-body techniques). The calibration is done by choosing appropriate free code parameters. We describe in particular the evolution of the binary, white dwarf and blue straggler populations, though not all channels for blue straggler formation are represented yet in our simulations. Calibrated Monte Carlo runs show good agreement with results of N-body simulations not only for global cluster parameters, but also for e.g. binary fraction, luminosity function and surface brightness. Comparison of Monte Carlo simulations with observational data for M67 shows that is possible to get reasonably good agreement between them. Unfortunately, because of the large statistical fluctuations of the numerical data and uncertainties in the observational data the inferred conclusions about the cluster initial conditions are not firm.

The majority of pulsar population synthesis studies performed to date have focused on isolated pulsar evolution. Those that have incorporated pulsar evolution within binary systems have tended to either treat binary evolution poorly of evolve the pulsar population in an ad-hoc manner. Here we present the first model of the Galactic field pulsar population that includes a comprehensive treatment of both binary and pulsar evolution. Synthetic observational surveys mimicking a variety of radio telescopes are then performed on this population. As such, a complete and direct comparison of model data with observations of the pulsar population within the Galactic disk is now possible. The tool used for completing this work is a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Here we give a brief overview of the method and assumptions involved with each component. Some preliminary results are also presented as well as plans for future applications of the code.

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.

We present the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.

There is currently much interest in the possible presence of intermediate-mass black holes in the cores of globular clusters. Based on theoretical arguments and simulation results it has previously been suggested that a large core radius -- or particularly a large ratio of the core radius to half-mass radius -- is a promising indicator for finding such a black hole in a star cluster. In this study N-body models of 100000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data is analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that is directly comparable to the Galactic globular cluster sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an intermediate-mass black hole. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a black-hole binary (comprised of stellar mass black holes) rather than a single massive black hole. This work does not rule out the existence of intermediate-mass black holes in the cores of at least some star clusters.

N-body simulations give us a rough idea of how the shape of a simulated object appears in three-dimensional space. From an observational point of view this may give us a misleading picture. The faint stars may be the most common stars in the system but the morphological information obtained by observations of an object may be dominated by the color properties of the bright stars. Due to dynamical effects, such as energy equipartition, different masses of stars may populate different regions in the object. Since stars are evolving in mass the stellar evolution may also influence the dynamics of the system. Hence, if one is interested in simulating what the morphology will look like through a telescope, one needs to simulate in addition evolving stars and weight them by their luminosity. Therefore we need to combine simulations of the dynamical evolution and a stellar population synthesis at the same time. For the dynamical evolution part we use a parallel version of a direct N-body code, NBODY6++. This code also includes stellar evolution. We link the parameters from this stellar evolution routine to the BaSeL 2.0 stellar library. For isolated clusters we found results assimilable to standard stellar population synthesis codes such as the PEGASE code. For clusters in a tidal field we found that the integrated colors become relatively bluer due to energy equipartition effects. In the time shortly before dissolution of the cluster the stellar M/L ratio becomes lower compared to isolated clusters. We compared the results of our simulations to integrated spectra of galactic globular clusters. For the cluster NGC 1851 we found good agreement between simulation and observation. For extragalactic clusters in M81 and M31 we looked at medium band spectral energy distributions and found for some clusters also a good agreement.

We investigate the evolution of binary fractions in star clusters using N-body models of up to 100000 stars. Primordial binary frequencies in these models range from 5% to 50%. Simulations are performed with the NBODY4 code and include a full mass spectrum of stars, stellar evolution, binary evolution and the tidal field of the Galaxy. We find that the overall binary fraction of a cluster almost always remains close to the primordial value, except at late times when a cluster is near dissolution. A critical exception occurs in the central regions where we observe a marked increase in binary fraction with time -- a simulation starting with 100000 stars and 5% binaries reached a core binary frequency as high as 40% at the end of the core-collapse phase (occurring at 16 Gyr with ~20000 stars remaining). Binaries are destroyed in the core by a variety of processes as a cluster evolves, but the combination of mass-segregation and creation of new binaries in exchange interactions produces the observed increase in relative number. We also find that binaries are cycled into and out of cluster cores in a manner that is analogous to convection in stars. For models of 100000 stars we show that the evolution of the core-radius up to the end of the initial phase of core-collapse is not affected by the exact value of the primordial binary frequency (for frequencies of 10% or less). We discuss the ramifications of our results for the likely primordial binary content of globular clusters.

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.

We present the results of a deep Hubble Space Telescope (HST) exposure of the nearby globular cluster NGC6397, focussing attention on the cluster's white dwarf cooling sequence. This sequence is shown to extend over 5 magnitudes in depth, with an apparent cutoff at magnitude F814W=27.6. We demonstrate, using both artificial star tests and the detectability of background galaxies at fainter magnitudes, that the cutoff is real and represents the truncation of the white dwarf luminosity function in this cluster. We perform a detailed comparison between cooling models and the observed distribution of white dwarfs in colour and magnitude, taking into account uncertainties in distance, extinction, white dwarf mass, progenitor lifetimes, binarity and cooling model uncertainties. After marginalising over these variables, we obtain values for the cluster distance modulus and age of \mu_0 = 12.02 \pm 0.06 and T_c = 11.47 \pm 0.47Gyr (95% confidence limits). Our inferred distance and white dwarf initial-final mass relations are in good agreement with other independent determinations, and the cluster age is consistent with, but more precise than, prior determinations made using the main sequence turnoff method. In particular, within the context of the currently accepted \Lambda CDM cosmological model, this age places the formation of NGC6397 at a redshift z=3, at a time when the cosmological star formation rate was approaching its peak.

NGC 6397 is the second closest globular star cluster to the Sun. Using 5 days of time on the Hubble Space Telescope, we have constructed the deepest ever color-magnitude diagram for this cluster. We see a clear truncation in each of its two major stellar sequences. Faint red main sequence stars run out well above our observational limit and near to the theoretical prediction for the lowest mass stars capable of stable hydrogen-burning in their cores. We also see a truncation in the number counts of faint blue stars, namely white dwarfs. This reflects the limit to which the bulk of the white dwarfs can cool over the lifetime of the cluster. There is also a turn towards bluer colors in the least luminous of these objects. This was predicted for the very coolest white dwarfs with hydrogen-rich atmospheres as the formation of H2 causes their atmospheres to become largely opaque to infrared radiation due to collision-induced absorption.

As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.

We perform binary population synthesis calculations to investigate the incidence of low-mass X-ray binaries and their birth rate in the Galaxy. We use a binary evolution algorithm that models all the relevant processes including tidal circularization and synchronization. Parameters in the evolution algorithm that are uncertain and may affect X-ray binary formation are allowed to vary during the investigation. We agree with previous studies that under standard assumptions of binary evolution the formation rate and number of black-hole low-mass X-ray binaries predicted by the model are more than an order of magnitude less than what is indicated by observations. We find that the common-envelope process cannot be manipulated to produce significant numbers of black-hole low-mass X-ray binaries. However, by simply reducing the mass-loss rate from helium stars adopted in the standard model, to a rate that agrees with the latest data, we produce a good match to the observations. Including low-mass X-ray binaries that evolve from intermediate-mass systems also leads to favourable results. We stress that constraints on the X-ray binary population provided by observations are used here merely as a guide as surveys suffer from incompleteness and much uncertainty is involved in the interpretation of results.

We describe a simple computer package which illustrates a method of combining stellar dynamics with stellar evolution. Though the method is intended for elaborate applications (especially the dynamical evolution of rich star clusters) it is illustrated here in the context of three-body scattering, i.e. interactions between a binary star and a field of single stars. We describe the interface between the dynamics and the two independent packages which describe the internal evolution of single stars and binaries. We also give an example application, and introduce a stand alone utility for the visual presentation of simulation results.

The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and binary evolution. Our preferred model starts with 12000 single stars and 12000 binaries placed in a Galactic tidal field at 8.0 kpc from the Galactic Centre. Our choices for the initial conditions and for the primordial binary population are explained in detail. At 4 Gyr, the age of M67, the total mass has reduced by 90% as a result of mass loss and stellar escapes. The mass and half-mass radius of luminous stars in the cluster are a good match to observations although the model is more centrally concentrated than observations indicate. The stellar mass and luminosity functions are significantly flattened by preferential escape of low-mass stars. We find that M67 is dynamically old enough that information about the initial mass function is lost, both from the current luminosity function and from the current mass fraction in white dwarfs. The model contains 20 blue stragglers at 4 Gyr which is slightly less than the 28 observed in M67. Nine are in binaries. The blue stragglers were formed by a variety of means and we find formation paths for the whole variety observed in M67. Both the primordial binary population and the dynamical cluster environment play an essential role in shaping the population. A substantial population of short-period primordial binaries (with periods less than a few days) is needed to explain the observed number of blue stragglers in M67.

Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst binary stellar populations with and without binary interactions. The ages of the populations are in the range $1-15 $Gyr and the metallicities are in the range $0.0001-0.03$. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the integrated $\rm U-B$, $\rm B-V$, $\rm V-R$ and $\rm R-I$ colours and all Lick/IDS spectral absorption indices (except for $\rm H_\beta$) substantially smaller. In other words binary evolution makes a population appear bluer. This effect raises the derived age and metallicity of the population. We calculate several sets of additional solar-metallicity binary stellar populations to explore the influence of input parameters to the binary evolution algorithm (the common-envelope ejection efficiency and the stellar wind mass-loss rate) on the resulting integrated colours. We also look at the dependence on the choice of distribution functions used to generate the initial binary population. The results show that variations in the choice of input model parameters and distributions can significantly affect the results. However, comparing the discrepancies that exist between the colours of various models, we find that the differences are less than those produced between the models with and those without binary interactions. Therefore it is very necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.

Using evolutionary population synthesis, we present integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst single stellar populations (SSPs). The ages of the SSPs are in the range 1-19 Gyr and the metallicities [Fe/H] are in the range -2.3 - 0.2. Our models use the rapid single stellar evolution algorithm of Hurley, Pols and Tout for the stellar evolutionary tracks, the empirical and semi-empirical calibrated BaSeL-2.0 model of Lejeune, Cuisinier and Buser for the library of stellar spectra and the empirical fitting functions of Worthey, Faber, Gonzalez and Burstein for the Lick/IDS spectral absorption-feature indices. Applying our synthetic Lick/IDS absorption-line indices to the merit function, we obtain the age and the metallicity of the central region of M32, which can be well explained by an instantaneous SSP with an age of 6.5 Gyr and a metallicity similar to solar. Applying the derived age and the metallicity from the merit function to a number of index-index diagrams, we find that the plots of Hbeta-Fe5015 and Hbeta-Fe5782 are the best index-index diagrams from which we can directly obtain reasonable age and metallicity.

We present integrated colours, integrated spectral energy distributions, and absorption-line indices, for instantaneous burst solar-metallicity binary stellar populations with ages in the range 1 - 15Gyr. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the appearance of the population substantially bluer -- this is the case for each of the quantities we have considered. This effect raises the derived age and metallicity of the population. Therefore it is necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.

We use N-body simulations of star clusters to investigate the possible dynamical origins of the observed spread in core radius among intermediate-age and old star clusters in the Large Magellanic Cloud (LMC). Two effects are considered, a time-varying external tidal field and variations in primordial hard binary fraction. Simulations of clusters orbiting a point-mass galaxy show similar core radius evolution for clusters on both circular and elliptical orbits and we therefore conclude that the tidal field of the LMC has not yet significantly influenced the evolution of the intermediate-age clusters. The presence of large numbers of hard primordial binaries in a cluster leads to core radius expansion; however, the magnitude of the effect is insufficient to explain the observations. Further, the range of binary fractions required to produce significant core radius growth is inconsistent with the observational evidence that all the LMC clusters have similar stellar luminosity functions.

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.

We use the results of realistic N-body simulations to investigate the appearance of the white dwarf population in dense star clusters. We show that the presence of a substantial binary population in a star cluster, and the interaction of this population with the cluster environment, has serious consequences for the morphology of the observed white dwarf sequence and the derived white dwarf cooling age of the cluster. We find that over time the dynamical evolution of the cluster -- mass-segregation, stellar interactions and tidal stripping -- hampers the use of white dwarfs as tracers of the initial mass function, and also leads to a significant enhancement of the white dwarf mass fraction. Future observations of star clusters should be conducted slightly interior to the half-mass radius of the cluster in order to best obtain information about the cluster age and initial mass function from the white dwarf luminosity function. The evolution of binary stars and the cluster environment must necessarily be accounted for when studying the white dwarf populations of dynamically evolved star clusters.

The progenitors of Type Ia supernovae (SNe Ia) have not been identified. Though they are no longer fashionable we investigate the consequences if a significant number of SNe Ia were edge-lit detonations (ELDs) of carbon/oxygen white dwarfs that have accreted a critical mass of helium. Our best understanding of the Phillips relation between light curve speed and peak luminosity assigns both these phenomena to the amount of $^{56}$Ni produced. In ELDs there are two sites of $^{56}$Ni synthesis. If the peak luminosity is determined primarily by the C/O ratio in the core it is primarily a function of its progenitor's initial mass. If the light curve decay speed is determined by the total mass of iron group elements ejected this is a function of the total mass of the ELD at the time of explosion. In general these two masses are correlated and an empirical relation between peak luminosity and light curve shape can be expected. However when we perform population synthesis for progenitors of different metallicities we find a systematic shift in this relation that would make distant SNe Ia fainter than those nearby. The abundances of alpha-rich isotopes, such as $^{44}$Ca, in the solar system indicate that only about 40 per cent of SNe Ia are edge-lit so any systematic effect that could be present would be correspondingly diluted. If we examine only the small subset of ELDs that accrete from a naked helium star, rather than a He white dwarf, the systematic effect disappears.

We find a remarkably enhanced production rate in star clusters (relative to the field) of very short period, massive double-white-dwarf stars and of giant-white dwarf binaries. These results are based on N-body simulations performed with the new GRAPE-6 special purpose hardware and are important in identifying and characterizing the progenitors of type Ia supernovae. The high incidence of very close double-white-dwarf systems is the result of dynamical encounters between (mostly) primordial binaries and other cluster stars. Orbital hardening rapidly drives these degenerate binaries to periods under ~10 hours. Gravitational radiation emission and mergers producing supra-Chandrasekhar objects follow in less than a Hubble time. If most stars are born in clusters then estimates of the double white dwarf merger rates in galaxies (due to cluster dynamical interaction) must be increased more than tenfold. A majority of the Roche lobe overflow giant-white dwarf binaries are not primordial; they are produced in exchange reactions. Most cases resulted in a common-envelope and formation of a double-white-dwarf binary rather than Supersoft X-ray sources leading possibly to a type Ia supernova.

We present a rapid binary evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but because orbits generally circularize before Roche-lobe overflow the outcome of the interactions of systems with the same semi-latus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries, however, the initial separations should be distributed according to the observed distribution of semi-latera recta rather than periods or semi-major axes.

The recent availability of special purpose computers designed for calculating gravitational interactions of N-bodies at extremely high speed has provided the means to model globular clusters on a star-by-star basis for the first time. By endeavouring to make the N-body codes that operate on these machines as realistic as possible, the addition of stellar evolution being one example, much is being learnt about the interaction between the star cluster itself and the stars it contains. A fascinating aspect of this research is the ability to follow the orbits of individual stars in detail and to document the formation of observed exotic systems. This has revealed that many stars within a star cluster lead wildly promiscuous lives, interacting, often intimately and in rapid succession, with a variety of neighbours.

We find that free-floating planets can remain bound to a star cluster for much longer than was previously assumed: of the order of the cluster half-mass relaxation timescale as opposed to the crossing-time. This result is based on N-body simulations performed with the new GRAPE-6 special purpose hardware and is important in the context of the preliminary detection of a population of free-floating sub-stellar objects in the globular cluster M22. The planets in our N-body study are of Jupiter mass and are initially placed in circular orbits of between 0.05 and 50 AU about a parent star whose mass is chosen from a realistic initial mass function. The presence of the free-floating planets is the result of dynamical encounters between planetary systems and the cluster stars. Most planets are liberated from their parent star in, or near, the cluster core, and then drift outwards on a timescale of ~ 10^{8-9} year. This still implies the existence of many (~100) planets per star if the M22 result is confirmed.

We derive stellar ages, from evolutionary tracks, and metallicities, from Stromgren photometry, for a sample of 5828 dwarf and sub-dwarf stars from the Hipparcos Catalogue. This stellar disk sample is used to investigate the age-metallicity diagram in the solar neighbourhood. Such diagrams are often used to derive a so called age-metallicity relation. Because of the size of our sample, we are able to quantify the impact on such diagrams, and derived relations, due to different selection effects. Some of these effects are of a more subtle sort, giving rise to erroneous conclusions. In particular we show that [1] the age-metallicity diagram is well populated at all ages and especially that old, metal-rich stars do exist, [2] the scatter in metallicity at any given age is larger than the observational errors, [3] the exclusion of cooler dwarf stars from an age-metallicity sample preferentially excludes old, metal-rich stars, depleting the upper right-hand corner of the age-metallicity diagram, [4] the distance dependence found in the Edvardsson et al. sample by Garnett & Kobulnicky is an expected artifact due to the construction of the original sample. We conclude that, although some of it can be attributed to stellar migration in the galactic disk, a large part of the observed scatter is intrinsic to the formation processes of stars.

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

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.

We present a state-of-the-art N-body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N-body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N-body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N = 10000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.

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.

We show how the second sequence seen lying above the main sequence in cluster colour magnitude diagrams results from binaries with a large range of mass ratios and not just from those with equal masses. We conclude that the presence of a densely populated second sequence, with only sparse filling in between it and the single star main sequence, does not necessarily imply that binary mass ratios are close to unity.

We have identified a candidate 1-2 x 10^5 year old luminous white dwarf in NGC 1818, a young star cluster in the Large Magellanic Cloud. This discovery strongly constrains the boundary mass M_c at which stars stop forming neutron stars and start forming white dwarfs, to M_c > 7.6 Msun.