Protoplanetary Discs are frisbees of matter formed around stars at the same time that solar systems are born.  They are the spawning grounds of the planets - we must unlock their secrets if we are to understand the solar systems we see in our Galaxy.  I use SPH simulations to study the evolution of these discs, in particular studying the evolution of gravitational instability.  Can discs form planets by fragmentation of the disc? We must explore how the instability varies with mass, radius, and thermodynamic regimes to answer this question fully.

I have augmented SPH with a radiative transfer algorithm to better model the disc thermodynamics.  We couple this with a basic equation of state to model the physical chemistry of the gas.  This is done at a minimal computational cost. 

I have applied this code to the issue of discs perturbed by stellar encounters.  If the disc cannot fragment in isolation, can it fragment under perturbation? We have confirmed that fragmentation is difficult both in isolation and under perturbation.  The primary factor in both cases is the thermodynamics - the ability of the disc to radiate and cool is most crucial.  We have also shown that encounters produce observational signatures cosmetically similar to an outburst (though too infrequent to explain all the outbursts we see).

Monte Carlo Radiative Transfer tracks the passage of photons through a medium using stochastic methods.  It is especially useful for hydrodynamic simulations as it allows the effects of a non-trivial density distribution to be traced.  We can then use the photons’ absorption and scattering in the medium to both reconstruct the system’s temperature distribution and create an image of the system (much like real telescopes image the heavens).

I have been using this method in protostellar disc scenarios.  In particular, I have been adapting this method for use in SPH density fields (typically, the method requires a density grid to operate - I have devised a method by which such gridding is not required).  This allows us to create observational signatures from our simulations, connecting us directly to observers.

The Search for Extraterrestrial Intelligence (SETI) has been largely influenced by the original numerical analyses of Frank Drake and Enrico Fermi (viz. the Drake Equation and Fermi’s Paradox respectively).  The current influx of data regarding extrasolar planets is providing new information to these arguments, but specific data is difficult to incorporate into such simplistic models.

I have developed a means of quantifying the number of intelligent civilisations in the Galaxy over its lifetime, using Monte Carlo Realisation Techniques.  Its advantage over traditional methods is that new exoplanet data can be incorporated immediately into its framework, and different hypotheses for the origin of life can hence be compared.