In this thesis I have studied how discs around young stars evolve and disperse. In particular, I build models which combine viscous evolution with photoevaporation, as previous work suggests that photoevaporation can reproduce the observed disc evolution and dispersal time-scales. The main question this thesis attempts to address is: Can photoevaporation provide a dominant dispersal mechanism for the observed population of young stars? Photoevaporation arises from the heating that high energy (UV and X-ray) photons provide to the surface layers of a disc. Before I started this work, only photoevaporation from a pure EUV radiation field was described within a hydrodynamic framework. Therefore, I start by building a hydrodynamic solution to the pure X-ray photoevaporation problem, and then extend this solution to the entire high energy spectrum. This hydrodynamic model leads me to conclude that it is the X-ray radiation field that sets the mass-loss rates. These mass-loss rates scale linearly with X-ray luminosity, are independent of the underlying disc structure and explicitly independent of stellar mass. I build a radiation-hydrodynamic algorithm, based on previous work, to describe the process of X-ray heating in discs. I then use this algorithm to span the full range of observed parameter space, to fully solve the X-ray photoevaporation problem. I further extend the algorithm to roughly approximate the heating an FUV radiation field would have on the photoevaporative flow, as well as separately testing the effect an EUV radiation field will have. These numerical tests are in agreement with the hydrodynamic model derived. Specifically, it is the X-rays that are driving the photoevaporative flow from the inner disc. Armed with an accurate description of the photoevaporative mass-loss rates from young stars, I consider the evolution of a population of disc-bearing, young (

Observations of gravitationally microlensed quasars offer a unique opportunity to probe quasar structure on extremely small scales. In this thesis, we conduct extensive microlensing simulations and compare with observational data to constrain quasar accretion discs, and conduct preliminary probes of broad emission line region structure. This analysis is done using a new single-epoch imaging technique that requires very little telescope time, and yet produces results that are comparable to those obtained from long-term monitoring campaigns. / We begin by exploring the impact of variable smooth matter percentage and source size on microlensing simulations. Adding a smooth matter component affects minimum and saddle point images differently, broadening the magnification distribution for the saddle point image significantly. However, increasing the radius of the background source washes out this difference. The observation of suppressed saddle point images can therefore only be explained by microlensing with a smooth matter component if the background source is sufficiently small. We use these simulations, in combination with I-band imaging of the lensed quasar MG 0414+0534 to constrain the radius of the quasar source. This demonstrates the viability of a single-epoch imaging method for constraining quasar structure. / This technique is then expanded to single-epoch multi-band observations, in order to constrain the radial profile of quasar accretion discs as a function of observed wavelength. We present new Magellan observations of two gravitationally lensed quasars: MG 0414+0534 and SDSS J0924+0219. We also analyse two epochs of Q2237+0305 data obtained from the literature. Our results are compared with four fidicial accretion disc models. At the 95 per cent level, only SDSS J0924+0219 is inconsistent with any of the accretion disc models. When we combine the results from all three quasars -- a first step towards assembling a statistical sample -- we find that the two steepest accretion disc models are ruled out with 68 per cent confidence. / In addition, we are also able to use our microlensing simulations to constrain the smooth matter percentages in the lenses at the image positions. In both MG 0414+0534 and SDSS J0924+0219 we find smooth matter percentages that are inconsistent with zero. A smooth matter percentage of approximately 50 per cent is preferred in MG 0414+0534, and approximately 80 per cent in SDSS J0924+0219. Q2237+0305 is usually assumed to have a smooth matter percentage of zero at the image positions, as they lie in the bulge of the lensing galaxy. Though our measurement is consistent with a zero smooth matter percentage, there is a peak in the probability distribution at a value 20 per cent. This is perhaps a hint of additional intervening structures along the line of sight to the background quasar. / We test the sensitivity of our accretion disc constraints to a range of modelling parameters. These include errors in lens modelling, Bayesian prior probability selection, errors in observational data, and the microlens mass function. Constraints on the power-law index relating source radius to observed wavelength are found to be relatively unaffected by changes in the modelling parameters. Constraints on source radii are found to be more strongly affected. / Finally, the broad emission line region of Q2237+0305 is examined. Gemini IFU observations are presented clearly showing differential microlensing across the velocity profile of the CIII] emission line. To analyse this signature, we present three simple broad emission line region models: a biconical outflow, a Keplerian disc, and spherical infall. A method is developed to compare the shapes of simulated flux ratio spectra with the observed spectrum. We are unable to discriminate between the biconical outflow and Keplerian disc models when we average over all viewing angles and orientations. The spherical infall model, however, does not fit the observed data. We also find that for the non-spherically symmetric geometries, low inclination angles are essentially incompatible with the observations. This analysis offers hope that with sufficiently high signal-to-noise observations, differential microlensing signatures may allow us to constrain the geometry and kinematics of this poorly understood portion of quasar structure.

In this thesis, an assessment is made of the value of optical CMDs as a useful diagnostic of the accretion properties of young stars. An analysis has been made of the phenomena that we observe and their effect on the position of stars in the CMD. Limitations and potential biases have been identified and evaluated. Variability causes some luminosity spread at a given colour in optical CMDs. A detailed characterisation of variability has been performed which places strong constraints on the magnitudes and the timescales on which the variability is seen. On timescales 15 minutes, almost no variability is detected (at levels greater than ≈ 0.2%) in the i band for a sample of ≈ 700 disc-bearing young stellar objects (YSOs). This suggests that the variability predicted by some accretion shock models is either very weak or not present. On hours to days timescales the optical variability in most stars is well described by a simple power law. The amplitude of the variability, a ∝ f−k, where f is the frequency of the variability in days. Disc-bearing and discless YSOs exhibit median values of k of 0.85 ± 0.02 and 0.95 ± 0.03 respectively, the uncertainity being the error on the median. The power law is valid up to a certain timescale (tmax) at which point the variability amplitude does not increase any further. tmax is found to be 1.50 ± 0.07 days and 1.41 ± 0.10 days for disc-bearing and discless stars respectively. Disc-bearing stars show greater variability amplitudes than the discless stars. However, it is notable that the variability timescale and power spectrum exponent are remarkably similar. This implies that the amplitude of the variability is driven by the physics of the underlying process, but that the timescales are instead driven by geometric effects. For disc-bearing stars, the highest amplitude variables are the accreting stars, which often appear to vary in the CMD along lines that correspond to changes in accretion luminosity. Four disc-bearing stars (approximately 0.5% of the disc-bearing sample) in Cep OB3b show extreme variability on timescales of years. Three (possible EXor candidates), show long-timescale changes that have a dramatic effect on their CMD position. However their small numbers mean that the overall impact on the CMDs of young associations is small. Variability on timescales of the rotational period and shorter adds uncertainty to age estimates of individual stars that are calculated by comparison with PMS models. Having provided a detailed description of variability and its impact on the CMD, it is clear that there are further significant mechanisms that affect the positions of YSOs in the CMD. I show that the spread in luminosity seen in the Orion Nebula Cluster and NGC 2264 could not be explained by accretion at rates of M ̇ ≥ 5 × 10−4 M⊙ yr−1 occurring within the protostellar phase of YSO evolution. Thus it appears that CMDs are not a useful diagnostic for study of the accretion histories of YSOs. The wavelength dependence of the extinction by dust within the inner regions of YSO discs is shown to differ from that seen in the ISM. Typically the wavelength dependence of the extinction is given by RV ≈5-8, compared with the value of RV ≈3.1 typical of the ISM. The interpretation is that grain growth has occurred. The location of this material within the ‘snow line’ implies that grains have coalesced rather than simply gaining an ice mantle. This is evidence for the beginning of planet formation. The effect of the high value of RV on the CMD is to add additional uncertainty of 0.1 mag to photometric measurements that have been corrected for the effects of extinction. Accretion luminosity is shown to be the dominant signal in the luminosity spread seen in CMDs of young associations. Stars which exhibit excess flux in the U band or Hα are displaced in CMD space. The accretion vector is shown to be a significant blueward shift in colour accompa- nied by a modest brightening in the g, g − i CMD. Accretion results in a luminosity spread as stars are displaced blueward below the PMS locus. This effect is not seen in non-accreting disc-bearing stars. Examination of the underlying excess luminosity spectrum for 15 accreting stars shows that the colour of the emission excess is not consistent across the sample. Thus to quantify the effect of accretion luminosity on CMD positions for individual stars, moderate resolution spectra are required with a large range in wavelength. This accretion luminosity may systematically bias estimates of PMS ages. A simple mitigation is to exclude accreting stars from age analysis. U band and Hα flux excesses are shown to vary independently by ≈ 1 dex on timescales shorter than the rotation period of the star. The relation between U band flux excess and veiling at 7000Å also appears to be variable. This implies that single epoch measurements of these parameters will add an uncertainty of ≈ 1 dex on accretion rates derived from them. Accretion rates derived from either U or Hα excess should be calculated from a mean of several photometric measurements, separated by significant fractions of the rotation period of the star. In most stars, the veiling at 7000Å is shown not to be a good measure for the calculation of the accretion rate. Despite providing a detailed characterisation of phenomena that influence the positions of YSOs in the CMD, there exists some residual luminosity spread at a given Teff that cannot be explained by variability on any timescale, extinction uncertainties or accretion luminosity. This residual spread should provide an opportunity to study an as-yet uncharacterised aspect of young stars.

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and secondly, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here.

Doctor of Philosophy / We investigate the linear growth and vertical structure of the magnetorotational instability (MRI) in weakly ionised, stratified accretion discs. The magnetic field is initially vertical and perturbations have vertical wavevectors only. Solutions are obtained at representative radial locations from the central protostar for different choices of the initial magnetic field strength, sources of ionisation, disc structure and configuration of the conductivity tensor. The MRI is active over a wide range of magnetic field strengths and fluid conditions in low conductivity discs. For the minimum-mass solar nebula model, incorporating cosmic ray and x-ray ionisation and assuming that charges are carried by ions and electrons only, perturbations grow at 1 AU for B < 8G. For a significant subset of these strengths (200mG < B < 5 G), the growth rate is of order the ideal MHD rate (0.75 Omega). Hall conductivity modifies the structure and growth rate of global unstable modes at 1 AU for all magnetic field strengths that support MRI. As a result, at this radius, modes obtained with a full conductivity tensor grow faster and are active over a more extended cross-section of the disc, than perturbations in the ambipolar diffusion limit. For relatively strong fields (e.g. B > 200 mG), ambipolar diffusion alters the envelope shapes of the unstable modes, which peak at an intermediate height, instead of being mostly flat as modes in the Hall limit are in this region of parameter space. Similarly, when cosmic rays are assumed to be excluded from the disc by the winds emitted by the magnetically active protostar, unstable modes grow at this radius for B < 2 G. For strong fields, perturbations exhibit a kink at the height where x-ray ionisation becomes active. Finally, for R = 5 AU (10 AU), unstable modes exist for B < 800 mG (B < 250 mG) and the maximum growth rate is close to the ideal-MHD rate for 20 mG < B < 500 mG (2 mG < B < 50 mG). Similarly, perturbations incorporating Hall conductivity have a higher wavenumber and grow faster than solutions in the ambipolar diffusion limit for B < 100 mG (B < 10 mG). Unstable modes grow even at the midplane for B > 100 mG (B ~ 1 mG), but for weaker fields, a small dead region exists. When a population of 0.1 um grains is assumed to be present, perturbations grow at 10 AU for B < 10 mG. We estimate that the figure for R = 1 AU would be of order 400 mG. We conclude that, despite the low magnetic coupling, the magnetic field is dynamically important for a large range of fluid conditions and field strengths in protostellar discs. An example of such magnetic activity is the generation of MRI unstable modes, which are supported at 1 AU for field strengths up to a few gauss. Hall diffusion largely determines the structure and growth rate of these perturbations for all studied radii. At radii of order 1 AU, in particular, it is crucial to incorporate the full conductivity tensor in the analysis of this instability, and more generally, in studies of the dynamics of astrophysical discs.

Doctor of Philosophy / We investigate the linear growth and vertical structure of the magnetorotational instability (MRI) in weakly ionised, stratified accretion discs. The magnetic field is initially vertical and perturbations have vertical wavevectors only. Solutions are obtained at representative radial locations from the central protostar for different choices of the initial magnetic field strength, sources of ionisation, disc structure and configuration of the conductivity tensor. The MRI is active over a wide range of magnetic field strengths and fluid conditions in low conductivity discs. For the minimum-mass solar nebula model, incorporating cosmic ray and x-ray ionisation and assuming that charges are carried by ions and electrons only, perturbations grow at 1 AU for B < 8G. For a significant subset of these strengths (200mG < B < 5 G), the growth rate is of order the ideal MHD rate (0.75 Omega). Hall conductivity modifies the structure and growth rate of global unstable modes at 1 AU for all magnetic field strengths that support MRI. As a result, at this radius, modes obtained with a full conductivity tensor grow faster and are active over a more extended cross-section of the disc, than perturbations in the ambipolar diffusion limit. For relatively strong fields (e.g. B > 200 mG), ambipolar diffusion alters the envelope shapes of the unstable modes, which peak at an intermediate height, instead of being mostly flat as modes in the Hall limit are in this region of parameter space. Similarly, when cosmic rays are assumed to be excluded from the disc by the winds emitted by the magnetically active protostar, unstable modes grow at this radius for B < 2 G. For strong fields, perturbations exhibit a kink at the height where x-ray ionisation becomes active. Finally, for R = 5 AU (10 AU), unstable modes exist for B < 800 mG (B < 250 mG) and the maximum growth rate is close to the ideal-MHD rate for 20 mG < B < 500 mG (2 mG < B < 50 mG). Similarly, perturbations incorporating Hall conductivity have a higher wavenumber and grow faster than solutions in the ambipolar diffusion limit for B < 100 mG (B < 10 mG). Unstable modes grow even at the midplane for B > 100 mG (B ~ 1 mG), but for weaker fields, a small dead region exists. When a population of 0.1 um grains is assumed to be present, perturbations grow at 10 AU for B < 10 mG. We estimate that the figure for R = 1 AU would be of order 400 mG. We conclude that, despite the low magnetic coupling, the magnetic field is dynamically important for a large range of fluid conditions and field strengths in protostellar discs. An example of such magnetic activity is the generation of MRI unstable modes, which are supported at 1 AU for field strengths up to a few gauss. Hall diffusion largely determines the structure and growth rate of these perturbations for all studied radii. At radii of order 1 AU, in particular, it is crucial to incorporate the full conductivity tensor in the analysis of this instability, and more generally, in studies of the dynamics of astrophysical discs.

Emission from X-ray binaries is variable on a wide range of timescales. On long timescales, changes in mass accretion rate drive changes in spectral state. There is also rapid variability, the power spectrum of which consists of a low frequency quasi-periodic oscillation (QPO) superimposed on a broad band noise continuum. Here I investigate a model intended to quantitatively explain the observed spectral and variability properties. I consider a truncated disc geometry whereby the inner regions of an optically thick, geometrically thin accretion disc evaporate to form an optically thin, large scale height accretion flow. The QPO is driven by Lense-Thirring precession of the entire hot flow and the broad band noise is due to fluctuations in mass accretion rate which propagate towards the central object. Mass conservation ties these two processes together, enabling me to define a model for the QPO and broad band noise which uses only one set of parameters. I am thus able fit the model to data. The accretion rate fluctuations drive fluctuations in the precession frequency, giving rise to a quasi-periodic oscillation rather than a pure periodicity. The model thus predicts recent observations which show the QPO frequency to correlate with flux on short timescales. I then investigate a more unique model prediction. As the flow precesses, the patch of the disc preferentially illuminated by the flow rotates such that a non face on observer sees a quasi-periodic shift between blue and red shift in the iron K alpha line. An observation of such an effect would constitute excellent evidence for the model.

We investigate the origin of the parsec-scale radio emission from the changing-look active galactic nucleus (AGN) of Mrk 590, and examine whether the radio power has faded concurrently with the dramatic decrease in accretion rates observed between the 1990s and the present. We detect a compact core at 1.6 and 8.4 GHz using new Very Long Baseline Array observations, finding no significant extended, jet-like features down to similar to 1 pc scales. The flat spectral index (alpha(8.4)(1.6) = 0.03) and high brightness temperature (T-b similar to 10(8) K) indicate self-absorbed synchrotron emission from the AGN. The radio to X-ray luminosity ratio of log(L-R/L-X) similar to -5, similar to that in coronally active stars, suggests emission from magnetized coronal winds, although unresolved radio jets are also consistent with the data. Comparing new Karl G. Jansky Very Large Array measurements with archival and published radio flux densities, we find 46 per cent, 34 per cent, and (insignificantly) 13 per cent flux density decreases between the 1990s and the year 2015 at 1.4 GHz, 5 GHz and 8.4 GHz, respectively. This trend, possibly due to the expansion and fading of internal shocks within the radio-emitting outflow after a recent outburst, is consistent with the decline of the optical-UV and X-ray luminosities over the same period. Such correlated variability demonstrates the AGN accretion-outflow connection, confirming that the changing-look behaviour in Mrk 590 originates from variable accretion rates rather than dust obscuration. The present radio and X-ray luminosity correlation, consistent with low/hard state accretion, suggests that the black hole may now be accreting in a radiatively inefficient mode.

Protoplanetary discs, composed of gas and dust, usually surround young stellar objects and serve two main purposes: they determine the accretion of matter onto the central object and also represent sites of planet formation. The accretion proceeds through the transport of angular momentum outwards allowing the disc matter to fall towards the centre. A mechanism responsible for the transport can be turbulence, waves or other coherent structures originating from various instabilities in discs that could, in addition, play a role in the planet formation process. For an understanding of these instabilities, it is necessary to study perturbation dynamics in differentially rotating, or sheared media. Thus, this thesis focuses on new aspects in the perturbation dynamics in non-magnetised protoplanetary discs that arise due to their self-gravity and velocity shear associated with the disc’s differential rotation. The analysis is carried out in the framework of the widely employed local shearing box approximation. We start with 2D discs and then move on to 3D ones. In 2D discs, there are two basic perturbation types/modes – spiral density waves and vortices – that are responsible for angular momentum transport and that can also contribute to accelerating planet formation. First, in the linear regime, we demonstrate that the vortical mode undergoes large growth due to self-gravity and in this process generates density waves via shear-induced linear mode coupling phenomenon. This is noteworthy, because commonly only density waves are considered in self-gravitating discs. Then we investigate vortex dynamics in the non-linear regime under the influence of self-gravity by means of numerical simulations. It is shown that vortices are no longer well-organised and long-lived structures, unlike those occurring in non-self-gravitating discs. They undergo recurring phases (lasting for a few disc rotation periods) of formation, growth and eventual destruction. We also discuss the dust trapping capability of such transient vortices. Perturbation dynamics in 3D vertically stratified discs is richer, as there are more mode types. We first consider non-axisymmetric modes in non-self-gravitating discs and then only axisymmetric modes in the more complicated case when self-gravity is present. Specifically, in non-self-gravitating discs with superadiabatic vertical stratification, motivated by the recent results on the transport properties of incompressible convection, we show that when compressibility is taken into account, the non-axisymmetric convective mode excites density waves via the same shear-induced linear mode coupling mechanism mentioned above. These generated density waves transport angular momentum outwards in the trailing phase, and we suggest that they may aid and enhance the transport due solely to convection in the non-linear regime, where the latter becomes outward. In the final part of the thesis, we carry out a linear analysis of axisymmetric vertical normal modes in stratified self-gravitating discs. Although axisymmetric modes do not display shear-induced couplings, their analysis provides insight into how gravitational instabilities develop in the 3D case and their onset criterion. We examine how the structure of dispersion curves and eigenfunctions of 3D modes are influenced by self-gravity, which mode first becomes gravitationally unstable and thus determines the onset criterion and nature of the gravitational instability in stratified discs. We also contrast the more exact instability criterion obtained with our 3D model with that of density waves in 2D discs. Based on these findings, we discuss the origin of 3D behaviour of perturbations involving noticeable disc surface distortions, as seen in some numerical simulations of self-gravitating discs.

In this thesis we analyse X-ray data of accretion powered low mass and high mass X-ray binaries to understand the nature of their accretion mechanisms by searching for some clues of viscous time-scales of their accretion discs, if they have, in their low frequency power density spectra created from their long-term X-ray observations, or by doing pulse timing analysis with much shorter X-ray data to detect the effects of torque fluctuations caused by the accreting material on the pulsar. The low mass and high mass X-ray binaries we analysed have breaks in their power density spectra, which are attributed to the role of viscosity in the formation of accretion discs. Although, the time-scales corresponding to these break frequencies are smaller than the predictions of the Standard theory of accretion discs, the sources give consistent results among themselves by displaying the expected correlation between their break and orbital frequencies. The correlation curve of LMXBs implies thicker appearing accretion discs than those assumed by the theory. The dichotomy of the HMXBs on this curve points out the different origins of accretion that these sources may have, and offers a way to distinguish the stellar-wind fed systems from the Roche-lobe overflow systems. The timing and spectral analysis of Swift J1626.6-5156 reveal a correlation between the spin-up rate and the luminosity of the source implying that the pulsar is accretion-powered. This correlation together with the characteristics of the X-ray spectra enables us to estimate the magnetic field and the distance of the source. The AXP 1E 2259+586 does not display any signs of viscous time-scale in its low frequency power density spectra, and its pulse timing analysis gives a much smaller torque noise value than that expected from accretion powered pulsars. In addition, the analysis results presented in this thesis reveal magnetar-like glitches which differ than those of radio pulsars, due to the presence of the strong magnetic field of the pulsar. These results eliminate the possibility that the AXP is an accretion-powered pulsar.