We present theoretical models for stellar black hole (BH) properties in young, massive star clusters. Using a Monte Carlo code for stellar dynamics, we model realistic star clusters with $N\simeq 5\times10^5$ stars and significant binary fractions (up to 50%) with self-consistent treatments of stellar dynamics and stellar evolution. We compute the formation rates and characteristic properties of single and binary BHs for various representative ages, cluster parameters, and metallicities. Because of dynamical interactions and supernova (SN) kicks, more single BHs end up retained in clusters compared to BHs in binaries. We also find that the ejection of BHs from a cluster is a strong function of initial density. In low-density clusters (where dynamical effects are negligible), it is mainly SN kicks that eject BHs from the cluster, whereas in high-density clusters (initial central density $\rho_c(0) \sim 10^5 \, M_\odot\, {\rm pc}^{-3} $ in our models) the BH ejection rate is enhanced significantly by dynamics. Dynamical interactions of binary systems in dense clusters also modify the orbital period and eccentricity distributions while also increasing the probability of a BH having a more massive companion.

A model for the formation and evolution of binary millisecond radio pulsars in systems with low mass companions (< 0.1 Msun) is investigated using a binary population synthesis technique. Taking into account the non conservative evolution of the system due to mass loss from an accretion disk as a result of propeller action and from the companion via ablation by the pulsar, the transition from the accretion powered to rotation powered phase is investigated. It is shown that the operation of the propeller and ablation mechanisms can be responsible for the formation and evolution of black widow millisecond pulsar systems from the low mass X-ray binary phase at an orbital period of ~0.1 day. For a range of population synthesis input parameters, the results reveal that a population of black widow millisecond pulsars characterized by orbital periods as long as ~0.4 days and companion masses as low as ~0.005 Msun can be produced. The orbital periods and minimum companion mass of this radio millisecond pulsar population critically depend on the thermal bloating of the semi-degenerate hydrogen mass losing component, with longer orbital periods for a greater degree of bloating. Provided that the radius of the companion is increased by about a factor of 2 relative to a fully degenerate, zero temperature configuration, an approximate agreement between observed long orbital periods and theoretical modeling of hydrogen rich donors can be achieved. We find no discrepancy between the estimated birth rates for LMXBs and black widow systems, which on average are $\sim~1.3\times10^{-5} {\rm yr}^{-1}$ and $1.3\times10^{-7} {\rm yr}^{-1}$ respectively.

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

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.

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.

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