We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.

Pebble accretion has become a popular component to core accretion models of planet formation, and is especially relevant to the formation of compact, resonant terrestrial planetary systems. Pebbles initially form in the inner protoplanetary disc, sweeping outwards in a radially expanding front, potentially forming planetesimals and planetary cores via migration and the streaming instability. This pebble front appears at early times, in what is typically assumed to be a low mass disc. We argue this picture is in conflict with the reality of young circumstellar discs, which are massive and self-gravitating. We apply standard pebble accretion and streaming instability formulae to self-gravitating protostellar disc models. Fragments will open a gap in the pebble disc, but they will likely fail to open a gap in the gas, and continue rapid inward migration. If this does not strongly perturb the pebble disc, our results show that disc fragments will accrete pebbles efficiently. We find that in general the pebble-to-gas-density ratio fails to exceed 0.01, suggesting that the streaming instability will struggle to operate. It may be possible to activate the instability if 10 cm grains are available, and spiral structures can effectively concentrate them in regions of low gravito-turbulence. If this occurs, lunar mass cores might be assembled on timescales of a few thousand years, **but this is likely to be rare, and is far from proven**. In any case, this work highlights the need for study of how self-gravitating protostellar discs define the distribution and properties of solid bodies, for future planet formation by core accretion.

The recent announcement of a Neptune-sized exomoon candidate orbiting the Jupiter-sized object Kepler-1625b has forced us to rethink our assumptions regarding both exomoons and their host exoplanets. In this paper I describe calculations of the habitable zone for Earthlike exomoons in orbit of Kepler-1625b under a variety of assumptions. I find that the candidate exomoon, Kepler-1625b-i, does not currently reside within the exomoon habitable zone, but may have done so when Kepler-1625 occupied the main sequence. If it were to possess its own moon (a "moon-moon") that was Earthlike, this could potentially have been a habitable world. If other exomoons orbit Kepler-1625b, then there are a range of possible semimajor axes/eccentricities that would permit a habitable surface during the main sequence phase, while remaining dynamically stable under the perturbations of Kepler-1625b-i. This is however contingent on effective atmospheric CO$_2$ regulation.

We present high resolution ($\sim$300 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment of MM 1 in 1.3 mm continuum emission and CH$_{3}$CN emission for the first time. The object divides into two main sources - MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0.57'' (1920 au) to the South-East. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of $0.141'' \times 0.050''$ ($480\times170$ au). The gas kinematics toward MM 1a probed via CH$_{3}$CN trace a variety of scales. The lower energy $J=12-11$ $K=3$ line traces extended, rotating gas within the outflow cavity, while the $v$8=1 line shows a clearly-resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is $40\pm5$ M$_{\odot}$ (where between 2.2-5.8 M$_{\odot}$ is attributed to the disk), while MM 1b is $<0.6$ M$_{\odot}$. The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star.

Simulations of massive star formation predict the formation of discs with significant substructure, such as spiral arms and clumps due to fragmentation. Here we present a semi-analytic framework for producing synthetic observations of discs with substructure, in order to determine their observability in interferometric observations. Unlike post-processing of hydrodynamical models, the speed inherent to our approach permits us to explore a large parameter space of star and disc parameters, and thus constrain properties for real observations. We compute synthetic dust continuum and molecular line observations probing different disc masses, distances, inclinations, thermal structures, dust distributions, and number and orientation of spirals and fragments. With appropriate spatial and kinematic filtering applied, our models predict that ALMA observations of massive YSOs at <5 kpc distances should detect spirals in both gas and dust in strongly self-gravitating discs (i.e. discs with up to two spiral arms and strong kinematic perturbations). Detecting spirals will be possible in discs of arbitrary inclination, either by directly spatially resolving them for more face-on discs (inclinations up to ~50 degrees), or through a kinematic signature otherwise. Clumps resulting from disc fragmentation should be detectable in the continuum, if the clump is sufficiently hotter than the surrounding disc material.

If the exomoon candidate orbiting Kepler-1625b truly exists, it is much more massive than the moons observed in the Solar system (Teachey et al. 2017; Teachey & Kipping 2018). This exomoon would be sufficiently large to stably host its own satellite. This has sparked discussion of a new category of celestial object - a moon-moon (Forgan 2018) or submoon (Kollmeier & Raymond 2018). In this Note, I describe initial results of climate modelling of a hypothetical moon-moon in the Kepler-1625b system, calculated using the OBERON code, which jointly computes 1D latitudinal energy balance models for individual worlds alongside the dynamical evolution of the system they inhabit (Forgan 2016a, DOI:10.5281/ZENODO.61236). Both the code and the parameter files used in these runs are available at github.com/dh4gan/oberon.

Recent high-sensitivity observations carried out with ALMA have revealed the presence of complex organic molecules (COMs) such as methyl cyanide (CH$_{\rm 3}$CN) and methanol (CH$_{\rm 3}$OH) in relatively evolved protoplanetary discs. The behaviour and abundance of COMs in earlier phases of disc evolution remains unclear. Here we combine a smoothed particle hydrodynamics simulation of a fragmenting, gravitationally unstable disc with a gas-grain chemical code. We use this to investigate the evolution of formamide (NH$_{\rm 2}$CHO), a pre-biotic species, in both the disc and in the fragments that form within it. Our results show that formamide remains frozen onto grains in the majority of the disc where the temperatures are $<$100 K, with a predicted solid-phase abundance that matches those observed in comets. Formamide is present in the gas-phase in three fragments as a result of the high temperatures ($\geq$200\,K), but remains in the solid-phase in one colder ($\leq$150 K) fragment. The timescale over which this occurs is comparable to the dust sedimentation timescales, suggesting that any rocky core which is formed would inherit their formamide content directly from the protosolar nebula.

We revisit an alternate explanation for the turbulent nature of molecular clouds - namely, that velocity dispersions matching classical predictions of driven turbulence can be generated by the passage of clumpy material through a shock. While previous work suggested this mechanism can reproduce the observed Larson relation between velocity dispersion and size scale ($\sigma \propto L^{\Gamma}$ with $\Gamma \approx 0.5$), the effects of self-gravity and magnetic fields were not considered. We run a series of smoothed particle magnetohydrodynamics experiments, passing clumpy gas through a shock in the presence of a combination of self-gravity and magnetic fields. We find powerlaw relations between $\sigma$ and $L$ throughout, with indices ranging from $\Gamma=0.3-1.2$. These results are relatively insensitive to the strength and geometry of magnetic fields, provided that the shock is relatively strong. $\Gamma$ is strongly sensitive to the angle between the gas' bulk velocity and the shock front, and the shock strength (compared to the gravitational boundness of the pre-shock gas). If the origin of the $\sigma-L$ relation is in clumpy shocks, deviations from the standard Larson relation constrain the strength and behaviour of shocks in spiral galaxies.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

The spiral waves detected in the protostellar disc surrounding Elias 2-27 have been suggested as evidence of the disc being gravitationally unstable. However, previous work has shown that a massive, stable disc undergoing an encounter with a massive companion are also consistent with the observations. We compare the spiral morphology of smoothed particle hydrodynamic simulations modelling both cases. The gravitationally unstable disc produces symmetric, tightly wound spiral arms with constant pitch angle, as predicted by the literature. The companion disc's arms are asymmetric, with pitch angles that increase with radius. However, these arms are not well-fitted by standard analytic expressions, due to the high disc mass and relatively low companion mass. We note that differences (or indeed similarities) in morphology between pairs of spirals is a crucial discriminant between scenarios for Elias 2-27, and hence future studies must fit spiral arms individually. If Elias 2-27 continues to show symmetric tightly wound spiral arms in future observations, then we posit that it is the first observed example of a gravitationally unstable protostellar disc.

A recent ALMA observation of the Elias 2-27 system revealed a two-armed structure extending out to ~300 au in radius. The protostellar disc surrounding the central star is unusually massive, raising the possibility that the system is gravitationally unstable. Recent work has shown that the observed morphology of the system can be explained by disc self-gravity, so we examine the physical properties of the disc necessary to detect self-gravitating spiral waves. Using three-dimensional Smoothed Particle Hydrodynamics, coupled with radiative transfer and synthetic ALMA imaging, we find that observable spiral structure can only be explained by self-gravity if the disc has a low opacity (and therefore efficient cooling), and is minimally supported by external irradiation. This corresponds to a very narrow region of parameter space, suggesting that, although it is possible for the spiral structure to be due to disc self-gravity, other explanations, such as an external perturbation, may be preferred.

We demonstrate numerical techniques for automatic identification of individual spiral arms in hydrodynamic simulations of astrophysical discs. Building on our earlier work, which used tensor classification to identify regions that were "spiral-like", we can now obtain fits to spirals for individual arm elements. We show this process can even detect spirals in relatively flocculent spiral patterns, but the resulting fits to logarithmic "grand-design" spirals are less robust. Our methods not only permit the estimation of pitch angles, but also direct measurements of the spiral arm width and pattern speed. In principle, our techniques will allow the tracking of material as it passes through an arm. Our demonstration uses smoothed particle hydrodynamics simulations, but we stress that the method is suitable for any finite-element hydrodynamics system. We anticipate our techniques will be essential to studies of star formation in disc galaxies, and attempts to find the origin of recently observed spiral structure in protostellar discs.

Human influence on the biosphere has been evident at least since the development of widespread agriculture, and some stratigraphers have suggested that the activities of modern civilization indicate a geological epoch transition. The study of the anthropocene as a geological epoch, and its implication for the future of energy-intensive civilizations, is an emerging transdisciplinary field in which astrobiology can play a leading role. Habitability research of Earth, Mars, and exoplanets examines extreme cases relevant for understanding climate change as a planetary process. Energy-intensive civilizations will also face thermodynamic limits to growth, which provides an important constraint for estimating the longevity of human civilization and guiding the search for extraterrestrial intelligence. We recommend that missions concepts such as LUVOIR, HabEx, and OST be pursued in order to make significant progress toward understanding the future evolution of life on our planet and the possible evolution of technological, energy-intensive life elsewhere in the universe.

We have previously argued that targeted interstellar communication has a physical optimum at narrowband X-ray wavelengths $\lambda\approx1\,$nm, limited by the surface roughness of focusing devices at the atomic level (arXiv:1711.05761). We search 24,247 archival X-ray spectra (of 6,454 unique objects) for such features and present 19 sources with monochromatic signals. Close examination reveals that these are most likely of natural origin. The ratio of artificial to natural sources must be <0.01%. This first limit can be improved in future X-ray surveys.

The optimal frequency for interstellar communication, using "Earth 2017" technology, was derived in papers I and II of this series (arXiv:1706.03795, arXiv:1706.05570). The framework included models for the loss of photons from diffraction (free space), interstellar extinction, and atmospheric transmission. A major limit of current technology is the focusing of wavelengths $\lambda<300\,$nm (UV). When this technological constraint is dropped, a physical bound is found at $\lambda\approx1\,$nm ($E\approx\,$keV) for distances out to kpc. While shorter wavelengths may produce tighter beams and thus higher data rates, the physical limit comes from surface roughness of focusing devices at the atomic level. This limit can be surpassed by beam-forming with electromagnetic fields, e.g. using a free electron laser, but such methods are not energetically competitive. Current lasers are not yet cost efficient at nm wavelengths, with a gap of two orders of magnitude, but future technological progress may converge on the physical optimum. We recommend expanding SETI efforts towards targeted (at us) monochromatic (or narrow band) X-ray emission at 0.5-2 keV energies.

It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragment's interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.

Upon entering a star system, light sails are subject to both gravitational forces and radiation pressure, and can use both in concert to modify their trajectory. Moreover, stars possess significant magnetic fields, and if the sail is in any way charged, it will feel the Lorentz force also. We investigate the dynamics of so-called "photogravimagnetic assists" of sailcraft around $\alpha$ Centauri A, a potential first destination en route to Proxima Centauri (the goal of the Breakthrough Starshot program). We find that a 10m$^2$ sail with a charge-to-mass-ratio of around 10 $\mu$C/g or higher will need to take account of magnetic field effects during orbital maneouvres. The magnetic field can provide an extra source of deceleration and deflection, and allow capture onto closer orbits around a target star. However, flipping the sign of the sailcraft's charge can radically change resulting trajectories, resulting in complex loop-de-loops around magnetic field lines and essentially random ejection from the star system. Even on well-behaved trajectories, the field can generate off-axis deflections at $\alpha$ Centauri that, while minor, can result in very poor targetting of the final destination (Proxima) post-assist. Fortunately for Breakthrough Starshot, nanosails are less prone to charging en route than their heavier counterparts, but can still accrue relatively high charge at both the origin and destination, when travelling at low speeds. Photogravimagnetic assists are highly non-trivial, and require careful course correction to mitigate against unwanted changes in trajectory.

The process of planet formation offers a rich source of dust production via grain growth in protostellar discs, and via grinding of larger bodies in debris disc systems. Chemical evolution models, designed to follow the build up of metals and dust in galaxies, do not currently account for planet formation. We consider the possibility that the apparent under-prediction of dust mass in high redshift galaxies by chemical evolution models could be in part, due to these models neglecting this process, specifically due to their assumption that a large fraction of the dust mass is removed from the interstellar medium during star formation (so-called astration). By adding a planet formation phase into galaxy chemical evolution, we demonstrate that the dust budget crisis can be partially ameliorated by a factor of 1.3-1.5 only if a) circumstellar discs prevent a large fraction of the dust mass entering the star during its birth, and b) that dust mass is preferentially liberated via jets, winds and outflows rather than accreted into planetary-mass bodies.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

Two fundamental problems for extraterrestrial intelligences (ETIs) attempting to establish interstellar communication are timing and energy consumption. Humanity's study of exoplanets via their transit across the host star highlights a means of solving both problems. An ETI 'A' can communicate with ETI 'B' if B is observing transiting planets in A's star system, either by building structures to produce artificial transits observable by B, or by emitting signals at B during transit, at significantly lower energy consumption than typical electromagnetic transmission schemes. This can produce a network of interconnected civilisations, establishing contact via observing each other's transits. Assuming that civilisations reside in a Galactic Habitable Zone (GHZ), I conduct Monte Carlo Realisation simulations of the establishment and growth of this network, and analyse its properties in the context of graph theory. I find that at any instant, only a few civilisations are correctly aligned to communicate via transits. However, we should expect the true network to be cumulative, where a "handshake" connection at any time guarantees connection in the future via e.g. electromagnetic signals. In all our simulations, the cumulative network connects all civilisations together in a complete network. If civilisations share knowledge of their network connections, the network can be fully complete on timescales of order a hundred thousand years. Once established, this network can connect any two civilisations either directly, or via intermediate civilisations, with a path much less than the dimensions of the GHZ.

We present a new method of identifying protostellar disc fragments in a simulation based on density derivatives, and analyse our data using this and the existing CLUMPFIND method, which is based on an ordered search over all particles in gravitational potential energy. Using smoothed particle hydrodynamics, we carry out 9 simulations of a $0.25$ M$_{\odot}$ disc around a 1 M$_{\odot}$ star, all of which fragment to form at least 2 bound objects. We find that when using all particles ordered in gravitational potential space, only fragments that survive the duration of the simulation are detected. When we use the density derivative method, all fragments are detected, so the two methods are complementary, as using the two methods together allows us to identify all fragments, and to then determine those that are likely to be destroyed. We find a tentative empirical relationship between the dominant azimuthal wavenumber in the disc $m$ and the maximum semi-major axis a fragment may achieve in a simulation, such that $a_{\mathrm{max}}\propto\frac{1}{m}$. We find the fragment destruction rate to be around half that predicted from population synthesis models. This is due to fragment-fragment interactions in the early gas phase of the disc, which can cause scattering and eccentricity pumping on short timescales, and affects the fragment's internal structure. We therefore caution that measurements of eccentricity as a function of semi-major axis may not necessarily constrain the formation mechanism of giant planets and brown dwarfs.

An exoplanet-exomoon system presents a superposition of phase curves to observers - the dominant component varies according to the planetary period, and the lesser varies according to both the planetary and the lunar period. If the spectra of the two bodies differs significantly, then it is likely there are wavelength regimes where the contrast between the moon and planet is significantly larger. In principle, this effect could be used to isolate periodic oscillations in the combined phase curve. Being able to detect the exomoon component would allow a characterisation of the exomoon radius, and potentially some crude atmospheric data. We run a parameter survey of combined exoplanet-exomoon phase curves, which show that for most sets of planet-moon parameters, the lunar component of the phase curve is undetectable to current state-of-the-art transit observations. Even with future transit survey missions, measuring the exomoon signal will most likely require photometric precision of 10 parts per million or better. The only exception to this is if the moon is strongly tidally heated or in some way self-luminous. In this case, measurements of the phase curve at wavelengths greater than a few microns can be dominated by the lunar contribution. Instruments like the James Webb Space Telescope and its successors are needed to make this method feasible.

Understanding the formation and evolution of giant planets ($\ge$1 $M_{Jup}$) at wide orbital separation ($\ge$5 AU) is one of the goals of direct imaging. Over the past 15 years, many surveys have placed strong constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including 3 stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75-5.7% at the 68% confidence level (CL) within 20-300 AU and 0.5-75 $M_{Jup}$, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario by Forgan & Rice (2013), with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1-75 $M_{Jup}$ range (95% CL). With the 3 sub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0-8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA, http://cesam.lam.fr/diva/), a public database that aims at gathering the results of past, present, and future direct imaging surveys.

We investigate how star formation efficiency can be significantly decreased by the removal of a molecular cloud's envelope by feedback from an external source. Feedback from star formation has difficulties halting the process in dense gas but can easily remove the less dense and warmer envelopes where star formation does not occur. However, the envelopes can play an important role keeping their host clouds bound by deepening the gravitational potential and providing a constraining pressure boundary. We use numerical simulations to show that removal of the cloud envelopes results in all cases in a fall in the star formation efficiency (SFE). At 1.38 free-fall times our 4 pc cloud simulation experienced a drop in the SFE from 16 to six percent, while our 5 pc cloud fell from 27 to 16 per cent. At the same time, our 3 pc cloud (the least bound) fell from an SFE of 5.67 per cent to zero when the envelope was lost. The star formation efficiency per free-fall time varied from zero to $\approx$ 0.25 according to $\alpha$, defined to be the ratio of the kinetic plus thermal to gravitational energy, and irrespective of the absolute star forming mass available. Furthermore the fall in SFE associated with the loss of the envelope is found to even occur at later times. We conclude that the SFE will always fall should a star forming cloud lose its envelope due to stellar feedback, with less bound clouds suffering the greatest decrease.

We investigate the role of magnetic fields in the fragmentation of self-gravitating discs using 3D global ideal magnetohydrodynamic simulations performed with the "phantom" smoothed particle hydrodynamics code. For initially toroidal fields, we find two regimes. In the first, where the cooling time is greater than five times the dynamical time, magnetic fields reduce spiral density wave amplitudes, which in turn suppresses fragmentation. This is the case even if the magnetic pressure is only a tenth of the thermal pressure. The second regime occurs when the cooling time is sufficiently short that magnetic fields cannot halt fragmentation. We find that magnetised discs produce more massive fragments, due to both the additional pressure exerted by the magnetic field, and the additional angular momentum transport induced by Maxwell stresses. The fragments are confined to a narrower range of initial semimajor axes than those in unmagnetised discs. The orbital eccentricity and inclination distributions of unmagnetised and magnetised disc fragments are similar. Our results suggest the fragmentation boundary could be at cooling times a factor of two lower than predicted by purely hydrodynamical models.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

The habitability of planets in binary star systems depends not only on the radiation environment created by the two stars, but also on the perturbations to planetary orbits and rotation produced by the gravitational field of the binary and neighbouring planets. Habitable planets in binaries may therefore experience significant perturbations in orbit and spin. The direct effects of orbital resonances and secular evolution on the climate of binary planets remain largely unconsidered. We present latitudinal energy balance modelling of exoplanet climates with direct coupling to an N Body integrator and an obliquity evolution model. This allows us to simultaneously investigate the thermal and dynamical evolution of planets orbiting binary stars, and discover gravito-climatic oscillations on dynamical and secular timescales. We investigate the Kepler-47 and Alpha Centauri systems as archetypes of P and S type binary systems respectively. In the first case, Earthlike planets would experience rapid Milankovitch cycles (of order 1000 years) in eccentricity, obliquity and precession, inducing temperature oscillations of similar periods (modulated by other planets in the system). These secular temperature variations have amplitudes similar to those induced on the much shorter timescale of the binary period. In the Alpha Centauri system, the influence of the secondary produces eccentricity variations on 15,000 year timescales. This produces climate oscillations of similar strength to the variation on the orbital timescale of the binary. Phase drifts between eccentricity and obliquity oscillations creates further cycles that are of order 100,000 years in duration, which are further modulated by neighbouring planets.

There have been several recent detections of candidate Keplerian discs around massive young protostars. Given the relatively large disc-to-star mass ratios in these systems, and their young ages, it is worth investigating their propensity to becoming self-gravitating. To this end, we compute self-consistent, semi-analytic models of putative self-gravitating discs for five candidate disc systems. Our aim is not to fit exactly the observations, but to demonstrate that the expected dust continuum emission from marginally unstable self-gravitating discs can be quite weak, due to high optical depth at the midplane even at millimetre wavelengths. In the best cases, the models produce "observable" disc masses within a factor of <1.5 of those observed, with midplane dust temperatures comparable to measured temperatures from molecular line emission. We find in two cases that a self-gravitating disc model compares well with observations. If these discs are self-gravitating, they satisfy the conditions for disc fragmentation in their outer regions. These systems may hence have as-yet-unresolved low mass stellar companions, and are thus promising targets for future high angular resolution observations.

The formation process of massive stars is not well understood, and advancement in our understanding benefits from high resolution observations and modelling of the gas and dust surrounding individual high-mass (proto)stars. Here we report sub-arcsecond (<1550 au) resolution observations of the young massive star G11.92-0.61 MM1 with the SMA and VLA. Our 1.3 mm SMA observations reveal consistent velocity gradients in compact molecular line emission from species such as CH$_3$CN, CH$_3$OH, OCS, HNCO, H$_2$CO, DCN and CH$_3$CH$_2$CN, oriented perpendicular to the previously reported bipolar molecular outflow from MM1. Modelling of the compact gas kinematics suggests a structure undergoing rotation around the peak of the dust continuum emission. The rotational profile can be well fit by a model of a Keplerian disc, including infall, surrounding an enclosed mass of 30-60M$_{\odot}$, of which 2-3M$_{\odot}$ is attributed to the disc. From modelling the CH$_3$CN emission, we determine that two temperature components, of 150 K and 230 K, are required to adequately reproduce the spectra. Our 0.9 and 3.0cm VLA continuum data exhibit an excess above the level expected from dust emission; the full centimetre-submillimetre wavelength spectral energy distribution of MM1 is well reproduced by a model including dust emission, an unresolved hypercompact H{\i}{\i} region, and a compact ionised jet. In combination, our results suggest that MM1 is an example of a massive proto-O star forming via disc accretion, in a similar way to that of lower mass stars.

Circumstellar discs are thought to be self-gravitating at very early times. If the disc is relatively cool, extended and accreting sufficiently rapidly, it can fragment into bound objects of order a few Jupiter masses and upwards. Given that the fragment's initial angular momentum is non-zero, and it will continue to accrete angular momentum from the surrounding circumstellar disc, we should expect that the fragment will also possess a relatively massive disc at early times. Therefore, we can ask: is disc fragmentation a hierarchical process? Or, can a disc fragment go on to produce its own self-gravitating circumfragmentary disc that produces daughter fragments? We investigate this using a set of nested 1D self-gravitating disc models. We calculate the radial structure of a marginally stable, self-gravitating circumstellar disc, and compute its propensity to fragmentation. We use this data to construct the local fragment properties at this radius. For each circumstellar disc model that results in fragmentation, we then compute a marginally stable self-gravitating circumfragmentary disc model. In general, the circumfragmentary discs are geometrically thick and truncated inside the Hill radius, and are hence stable against daughter fragmentation. The typical steady-state accretion rate is between 0.1 and 10 percent of the local circumstellar disc accretion rate. The lifetime of the circumfragmentary discs' self-gravitating phase is somewhat less than 0.1 Myr, quite comparable with that of the circumstellar disc. We should therefore expect that disc fragments will not produce satellites via gravitational instability, but equally the self-gravitating phase of circumfragmentary discs is likely to affect the properties of satellites subsequently formed via core accretion.

We use a simple, self-consistent, self-gravitating semi-analytic disc model to conduct an examination of the parameter space in which self-gravitating discs may exist. We then use Monte-Carlo radiative transfer to generate synthetic ALMA images of these self-gravitating discs to determine the subset of this parameter space in which they generate non-axisymmetric structure that is potentially detectable by ALMA. Recently, several transition discs have been observed to have non-axisymmetric structure that extends out to large radii. It has been suggested that one possible origin of these asymmetries could be spiral density waves induced by disc self-gravity. We use our simple model to see if these discs exist in the region of parameter space where self-gravity could feasibly explain these spiral features. We find that for self-gravity to play a role in these systems typically requires a disc mass around an order of magnitude higher than the observed disc masses for the systems. The spiral amplitudes produced by self-gravity in the local approximation are relatively weak when compared to amplitudes produced by tidal interactions, or spirals launched at Lindblad resonances due to embedded planets in the disc. As such, we ultimately caution against diagnosing spiral features as being due to self-gravity, unless the disc exists in the very narrow region of parameter space where the spiral wave amplitudes are large enough to produce detectable features, but not so large as to cause the disc to fragment.

As hydrodynamic simulations increase in scale and resolution, identifying structures with non-trivial geometries or regions of general interest becomes increasingly challenging. There is a growing need for algorithms that identify a variety of different features in a simulation without requiring a "by-eye" search. We present tensor classification as such a technique for smoothed particle hydrodynamics (SPH). These methods have already been used to great effect in N-Body cosmological simulations, which require smoothing defined as an input free parameter. We show that tensor classification successfully identifies a wide range of structures in SPH density fields using its native smoothing, removing a free parameter from the analysis and preventing the need for tesselation of the density field, as required by some classification algorithms. As examples, we show that tensor classification using the tidal tensor and the velocity shear tensor successfully identifies filaments, shells and sheet structures in giant molecular cloud simulations, as well as spiral arms in discs. The relationship between structures identified using different tensors illustrates how different forces compete and co-operate to produce the observed density field. We therefore advocate the use of multiple tensors to classify structure in SPH simulations, to shed light on the interplay of multiple physical processes.

The habitable zone for exomoons with Earth-like properties is a non-trivial manifold, compared to that of Earth-like exoplanets. The presence of tidal heating, eclipses and planetary illumination in the exomoon energy budget combine to produce both circumstellar and circumplanetary habitable regions. Analytical calculations suggest that the circumplanetary habitable region is defined only by an inner edge (with its outer limits determined by orbital stability). Subsequent calculations using 1D latitudinal climate models indicated that the combined effect of eclipses and ice-albedo feedback can produce an outer edge to the circumplanetary habitable zone. But is this outer edge real, or an artefact of the climate model's relative simplicity? We present an upgraded 1D climate model of Earth-like exomoon climates, containing the carbonate-silicate cycle and viscoelastic tidal heating. We conduct parameter surveys of both the circumstellar and circumplanetary habitable zones, and we find that the outer circumplanetary habitable edge remains provided the moon's orbit is not inclined relative to that of the planet. Adding the carbonate-silicate cycle pushes the circumplanetary habitable zone outward, by allowing increases in atmospheric partial pressure of carbon dioxide to boost the greenhouse effect. Viscoelastic tidal heating widens the habitable zone compared to standard, fixed-Q models. Weakening the tidal heating effect due to melting allows moons to be habitable at higher eccentricity, and pushes the inner circumstellar and circumplanetary habitable zone boundary inward.

We present the first model that couples high-resolution simulations of the formation of Local Group galaxies with calculations of the galactic habitable zone (GHZ), a region of space which has sufficient metallicity to form terrestrial planets without being subject to hazardous radiation. These simulations allow us to make substantial progress in mapping out the asymmetric three-dimensional GHZ and its time evolution for the Milky Way (MW) and Triangulum (M33) galaxies, as opposed to works that generally assume an azimuthally symmetric GHZ. Applying typical habitability metrics to MW and M33, we find that while a large number of habitable planets exist as close as a few kiloparsecs from the galactic centre, the probability of individual planetary systems being habitable rises as one approaches the edge of the stellar disc. Tidal streams and satellite galaxies also appear to be fertile grounds for habitable planet formation. In short, we find that both galaxies arrive at similar GHZs by different evolutionary paths, as measured by the first and third quartiles of surviving biospheres. For the Milky Way, this interquartile range begins as a narrow band at large radii, expanding to encompass much of the galaxy at intermediate times before settling at a range of 2-13kpc. In the case of M33, the opposite behaviour occurs - the initial and final interquartile ranges are quite similar, showing gradual evolution. This suggests that galaxy assembly history strongly influences the time evolution of the GHZ, which will affect the relative time lag between biospheres in different galactic locations. We end by noting the caveats involved in such studies and demonstrate that high resolution cosmological simulations will play a vital role in understanding habitability on galactic scales, provided that these simulations accurately resolve chemical evolution.

It is now reasonably clear that disc fragmentation can only operate in the outer parts of protostellar discs ($r > 50$ au). It is also expected that any object that forms via disc fragmentation will have an initial mass greater than that of Jupiter. However, whether or not such a process actually operates, or can play a significant role in the formation of planetary-mass objects, is still unclear. We do have a few examples of directly imaged objects that may have formed in this way, but we have yet to constrain how often disc fragmentation may actually form such objects. What we want to consider here is whether or not we can constrain the likely population of planetary-mass objects formed via disc fragmentation by considering how a population of objects at large radii ($a > 50$) au - if they do exist - would evolve under perturbations from more distant stellar companions. We find that there is a specific region of parameter space to which such objects would be scattered and show that the known exoplanets in that region have properties more consistent with that of the bulk exoplanet population, than with having been formed via disc fragmentation at large radii. Along with the scarcity of directly-imaged objects at large radii, our results provide a similar, but independent, constraint on the frequency of objects formed via disc fragmentation.

We address the possibility that intelligent civilisations that destroy themselves could present signatures observable by humanity. Placing limits on the number of self-destroyed civilisations in the Milky Way has strong implications for the final three terms in Drake's Equation, and would allow us to identify which classes of solution to Fermi's Paradox fit with the evidence (or lack thereof). Using the Earth as an example, we consider a variety of scenarios in which humans could extinguish their own technological civilisation. Each scenario presents some form of observable signature that could be probed by astronomical campaigns to detect and characterise extrasolar planetary systems. Some observables are unlikely to be detected at interstellar distances, but some scenarios are likely to produce significant changes in atmospheric composition that could be detected serendipitously with next-generation telescopes. In some cases, the timing of the observation would prove crucial to detection, as the decay of signatures is rapid compared to humanity's communication lifetime. In others, the signatures persist on far longer timescales.

The gravitational instability model of planet/brown dwarf formation proposes that protostellar discs can fragment into objects with masses above a few Jupiter masses at large semimajor axis. Tidal downsizing may reduce both the object mass and semimajor axis. However, most studies of tidal downsizing end when the protostellar disc disperses, while the system is embedded in its parent star-forming region. To compare disc fragment descendants with exoplanet and brown dwarf observations, the subsequent dynamical evolution must be explored. We carry out N-Body integrations of fragment-fragment scattering in multi-object star systems, and star systems embedded in substructured clusters. In both cases, we use initial conditions generated by population synthesis models of tidal downsizing. The scattering simulations produce a wide range of eccentricities. The ejection rate is around 25%. The ejecta mass distribution is similar to that for all objects, with a velocity dispersion consistent with those produced by full hydrodynamic simulations. The semimajor axis distribution after scattering extends to parsec scales. In the cluster simulations, 13% of objects are ejected from their planetary system, and around 10% experience significant orbit modification. A small number of objects are recaptured on high eccentricity, high inclination orbits. The velocity distribution of ejecta is similar to that produced by fragment-fragment scattering. If fragment-fragment scattering and cluster stripping act together, then disc fragmentation should be efficient at producing free-floating substellar objects, and hence characterising the free-floating planet population will provide strong constraints on the frequency of disc fragmentation.

The Optical Search for Extraterrestrial Intelligence (OSETI) attempts to detect collimated, narrowband pulses of electromagnetic radiation. These pulses may either consist of signals intentionally directed at the Earth, or signals between two star systems with a vector that unintentionally intersects the Solar System, allowing Earth to intercept the communication. But should we expect to be able to intercept these unintentional signals? And what constraints can we place upon the frequency of intelligent civilisations if we do? We carry out Monte Carlo Realisation simulations of interstellar communications between civilisations in the Galactic Habitable Zone (GHZ) using collimated beams. We measure the frequency with which beams between two stars are intercepted by a third. The interception rate increases linearly with the fraction of communicating civilisations, and as the cube of the beam opening angle, which is somewhat stronger than theoretical expectations, which we argue is due to the geometry of the GHZ. We find that for an annular GHZ containing 10,000 civilisations, intersections are unlikely unless the beams are relatively uncollimated. These results indicate that optical SETI is more likely to find signals deliberately directed at the Earth than accidentally intercepting collimated communications. Equally, civilisations wishing to establish a network of communicating species may use weakly collimated beams to build up the network through interception, if they are willing to pay a cost penalty that is lower than that meted by fully isotropic beacons. Future SETI searches should consider the possibility that communicating civilisations will attempt to strike a balance between optimising costs and encouraging contact between civilisations, and look for weakly collimated pulses as well as narrow-beam pulses directed deliberately at the Earth.

Recently, the Kepler Space Telescope has detected several planets in orbit around a close binary star system. These so-called circumbinary planets will experience non-trivial spatial and temporal distributions of radiative flux on their surfaces, with features not seen in their single-star orbiting counterparts. Earthlike circumbinary planets inhabited by photosynthetic organisms will be forced to adapt to these unusual flux patterns. We map the flux received by putative Earthlike planets (as a function of surface latitude/longitude and time) orbiting the binary star systems Kepler-16 and Kepler-47, two star systems which already boast circumbinary exoplanet detections. The longitudinal and latitudinal distribution of flux is sensitive to the centre of mass motion of the binary, and the relative orbital phases of the binary and planet. Total eclipses of the secondary by the primary, as well as partial eclipses of the primary by the secondary add an extra forcing term to the system. We also find that the patterns of darkness on the surface are equally unique. Beyond the planet's polar circles, the surface spends a significantly longer time in darkness than latitudes around the equator, due to the stars' motions delaying the first sunrise of spring (or hastening the last sunset of autumn). In the case of Kepler-47, we also find a weak longitudinal dependence for darkness, but this effect tends to average out if considered over many orbits. In the light of these flux and darkness patterns, we consider and discuss the prospects and challenges for photosynthetic organisms, using terrestrial analogues as a guide.

From analytical studies of tidal heating, eclipses and planetary illumination, it is clear that the exomoon habitable zone (EHZ) - the set of moon and host planet orbits that permit liquid water on an Earthlike moon's surface - is a manifold of higher dimension than the planetary HZ. This paper outlines the first attempt to produce climate models of exomoons which possess all the above sources and sinks of energy. We expand on our previous 1D latitudinal energy balance models (LEBMs), which follow the evolution of the temperature on an Earthlike moon orbiting a Jupiterlike planet, by adding planetary illumination. We investigate the EHZ in four dimensions, running two separate suites of simulations. The first investigates the EHZ by varying the planet's orbit, keeping the moon's orbit fixed, to compare the EHZ with planetary habitable zones. In general, planetary illumination pushes EHZs slightly further away from the star. Secondly, we fix the planet's orbit and vary the moon's orbit, to investigate the circumplanetary inner habitable edge. We demonstrate that an outer edge can exist due to eclipses (rather than merely orbital stability), but this edge may be pushed outwards when the effect of the carbonate-silicate cycle is taken into account.

Photosynthetic life requires sufficient photosynthetically active radiation (PAR) to metabolise. On Earth, plant behaviour, physiology and metabolism are sculpted around the night-day cycle by an endogenous biological circadian clock. The evolution of life was influenced by the Earth-Sun orbital dynamic, which generates the photo-environment incident on the planetary surface. In this work the unusual photo-environment of an Earth-like planet (ELP) in 3:2 spin orbit resonance is explored. Photo-environments on the ELP are longitudinally differentiated, in addition to differentiations relating to latitude and depth (for aquatic organisms) which are familiar on Earth. The light environment on such a planet could be compatible with Earth's photosynthetic life although the threat of atmospheric freeze-out and prolonged periods of darkness would present significant challenges. We emphasise the relationship between the evolution of life on a planetary body with its orbital dynamics.

We present new observations of the XZ Tau system made at high angular resolution (55 milliarcsec) with the Karl G. Jansky Very Large Array (VLA) at a wavelength of 7 mm. Observations of XZ Tau made with the VLA in 2004 appeared to show a triple system, with XZ Tau A resolved into two sources, XZ Tau A and XZ Tau C. The angular separation of XZ Tau A and C (0.09 arcsec) suggested a projected orbital separation of around 13 AU with a possible orbital period of around 40 yr. Our observations were obtained approximately 8 yr later, a fifth of this putative orbital period, and should therefore allow us to constrain the orbit of XZ Tau C, and evaluate the possibility that a recent periastron passage of C coincided with the launch of extended optical outflows from XZ Tau A. Despite improved sensitivity and resolution compared with previous observations, we find no evidence of XZ Tau C in our data. Components A and B are detected with a signal-to-noise ratio greater than ten; their orbits are consistent with previous studies, although the emission from XZ Tau A appears to be weaker. Three possible interpretations are offered: either XZ Tau C is transiting XZ Tau A, which is broadly consistent with the periastron passage hypothesis, or the emission seen in 2004 was that of a transient, or XZ Tau C does not exist. A fourth interpretation, that XZ Tau C was ejected from the system, is dismissed due to the lack of angular momentum redistribution in the orbits of XZ Tau A and XZ Tau B that would result from such an event. Our observations are insufficient to distinguish between the remaining possibilities, at least not until we obtain further VLA observations at a sufficiently later time. A further non-detection would allow us to reject the transit hypothesis, and the periastron passage of XZ Tau C as agent of XZ Tau A's outflows.

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

In the light of current scientific assessments of human-induced climate change, we investigate an experimental model to inform how resource-use strategies may influence interplanetary and interstellar colonisation by intelligent civilisations. In doing so, we seek to provide an additional aspect for refining the famed Fermi Paradox. The model described is necessarily simplistic, and the intent is to simply obtain some general insights to inform and inspire additional models. We model the relationship between an intelligent civilisation and its host planet as symbiotic, where the the relationship between the symbiont and the host species (the civilisation and the planets ecology, respectively) determines the fitness and ultimate survival of both organisms. We perform a series of Monte Carlo Realisation simulations, where civilisations pursue a variety of different relationships/strategies with their host planet, from mutualism to parasitism, and can consequently 'infect' other planets/hosts. We find that parasitic civilisations are generally less effective at survival than mutualist civilisations, provided that interstellar colonisation is inefficient (the maximum velocity of colonisation/infection is low). However, as the colonisation velocity is increased, the strategy of parasitism becomes more successful, until they dominate the 'population'. This is in accordance with predictions based on island biogeography and r/K selection theory. While heavily assumption dependent, we contend that this provides a fertile approach for further application of insights from theoretical ecology for extraterrestrial colonisation - while also potentially offering insights for understanding the human-Earth relationship and the potential for extraterrestrial human colonisation.

Previous attempts to describe circumbinary habitable zones have been concerned with the spatial extent of the zone, calculated analytically according to the combined radiation field of both stars. By contrast to these "spatial HZs", we present a numerical analysis of the "orbital HZ", a habitable zone defined as a function of planet orbital elements. This orbital HZ is better equipped to handle (for example) eccentric planet orbits, and is more directly connected to the data returned by exoplanet observations. Producing an orbital HZ requires a large number of climate simulations to be run to investigate the parameter space - we achieve this using Latitudinal Energy Balance Models (LEBMs), which handle the insolation of the planet by both stars (including mutual eclipses), as well as the planetary atmosphere's ability to absorb, transfer and lose heat. We present orbital HZs for several known circumbinary planetary systems: Kepler-16, Kepler-34, Kepler-35, Kepler-47 and PH-1. Generally, the orbital HZs at zero eccentricity are consistent with spatial HZs derived by other authors, although we detect some signatures of variability that coincide with resonances between the binary and planet orbital periods. We confirm that Earthlike planets around Kepler-47 with Kepler-47c's orbital parameters could possess liquid water, despite current uncertainties regarding its eccentricity. Kepler-16b is found to be outside the habitable zone, as well as the other circumbinary planets investigated.

Kepler-22b is the first transiting planet to have been detected in the habitable-zone of its host star. At 2.4 Earth radii, Kepler-22b is too large to be considered an Earth-analog, but should the planet host a moon large enough to maintain an atmosphere, then the Kepler-22 system may yet possess a telluric world. Aside from being within the habitable-zone, the target is attractive due to the availability of previously measured precise radial velocities and low intrinsic photometric noise, which has also enabled asteroseismology studies of the star. For these reasons, Kepler-22b was selected as a target-of-opportunity by the 'Hunt for Exomoons with Kepler' (HEK) project. In this work, we conduct a photodynamical search for an exomoon around Kepler-22b leveraging the transits, radial velocities and asteroseismology plus several new tools developed by the HEK project to improve exomoon searches. We find no evidence for an exomoon around the planet and exclude moons of mass >0.5 Earth masses to 95% confidence. By signal injection and blind retrieval, we demonstrate that an Earth-like moon is easily detected for this planet even when the time-correlated noise of the data set is taken into account. We provide updated parameters for the planet Kepler-22b including a revised mass of <53 Earth masses to 95% confidence and an eccentricity of 0.13(-0.13)(+0.36) by exploiting Single-body Asterodensity Profiling (SAP). Finally, we show that Kepler-22b has a >95% probability of being within the empirical habitable-zone but a <5% probability of being within the conservative habitable-zone.

Interstellar probes can carry out slingshot manoeuvres around the stars they visit, gaining a boost in velocity by extracting energy from the star's motion around the Galactic Centre. These maneouvres carry little to no extra energy cost, and in previous work it has been shown that a single Voyager-like probe exploring the galaxy does so 100 times faster when carrying out these slingshots than when navigating purely by powered flight (Forgan et al. 2012). We expand on these results by repeating the experiment with self-replicating probes. The probes explore a box of stars representative of the local Solar neighbourhood, to investigate how self-replication affects exploration timescales when compared with a single non-replicating probe. We explore three different scenarios of probe behaviour: i) standard powered flight to the nearest unvisited star (no slingshot techniques used), ii) flight to the nearest unvisited star using slingshot techniques, and iii) flight to the next unvisited star that will give the maximum velocity boost under a slingshot trajectory. In all three scenarios we find that as expected, using self-replicating probes greatly reduces the exploration time, by up to three orders of magnitude for scenario i) and iii) and two orders of magnitude for ii). The second case (i.e. nearest-star slingshots) remains the most time effective way to explore a population of stars. As the decision-making algorithms for the fleet are simple, unanticipated "race conditions" amongst probes are set up, causing the exploration time of the final stars to become much longer than necessary. From the scaling of the probes' performance with star number, we conclude that a fleet of self-replicating probes can indeed explore the Galaxy in a sufficiently short time to warrant the existence of the Fermi Paradox.

The Class A stellar engine (also known as a Shkadov thruster) is a spherical arc mirror, designed to use the impulse from a star's radiation pressure to generate a thrust force, perturbing the star's motion. If this mirror obstructs part of the stellar disc during the transit of an exoplanet, then this may be detected by studying the shape of the transit light curve, presenting another potential means by which the action of extraterrestrial intelligence (ETI) can be discerned. We model the light curves produced by exoplanets transiting a star which possesses a Shkadov thruster, and show how the parameters of the planet and the properties of the thruster can be disentangled. provided that radial velocity follow-up measurements are possible, and that other obscuring phenomena typical to exoplanet transit curves (such as the presence of starspots or intrinsic stellar noise) do not dominate. These difficulties aside, we estimate the a priori probability of detecting a Shkadov thruster during an exoplanet transit, which even given optimistic assumptions remains stubbornly low. Despite this, many exoplanet transit surveys designed for radial velocity follow-up are on the horizon, so we argue that this remains a useful serendipitous SETI technique. At worst, this technique will place an upper limit on the number of Class A stellar engines in the Solar neighbourhood; at best, this could help identify unusual transiting exoplanet systems as candidates for further investigation with other SETI methods.

Recent observations of the Class 0 protostar L1527 IRS have revealed a rotationally supported disc with an outer radius of at least 100 au. Measurements of the integrated flux at 870 microns suggest a disc mass that is too low for gravitational instability to govern angular momentum transport. However, if parts of the disc are optically thick at sub-mm wavelengths, the sub-mm fluxes will underestimate the disc mass, and the disc's actual mass may be substantially larger, potentially sufficient to be self-gravitating. We investigate this possibility using simple self-gravitating disc models. To match the observed mass accretion rates requires a disc-to-star mass ratio of at least ~0.5, which produces sub-mm fluxes that are similar to those observed for L1527 IRS in the absence of irradiation from the envelope or central star. If irradiation is significant, then the predicted fluxes exceed the observed fluxes by around an order of magnitude. Our model also indicates that the stresses produced by the gravitational instability are low enough to prevent disc fragmentation. As such, we conclude that observations do not rule out the possibility that the disc around L1527 IRS is self-gravitating, but it is more likely that despite being a very young system, this disc may already have left the self-gravitating phase.

Recently, the gravitational instability (GI) model of giant planet and brown dwarf formation has been revisited and recast into what is often referred to as the "tidal downsizing" hypothesis. The fragmentation of self-gravitating protostellar discs into gravitationally bound embryos - with masses of a few to tens of Jupiter masses, at semi major axes above 30 - 40 au - is followed by grain sedimentation inside the embryo, radial migration towards the central star and tidal disruption of the embryo's upper layers. The properties of the resultant object depends sensitively on the timescales upon which each process occurs. Therefore, GI followed by tidal downsizing can theoretically produce objects spanning a large mass range, from terrestrial planets to giant planets and brown dwarfs. Whether such objects can be formed in practice, and what proportions of the observed population they would represent, requires a more involved statistical analysis. We present a simple population synthesis model of star and planet formation via GI and tidal downsizing. We couple a semi-analytic model of protostellar disc evolution to analytic calculations of fragmentation, initial embryo mass, grain growth and sedimentation, embryo migration and tidal disruption. While there are key pieces of physics yet to be incorporated, it represents a first step towards a mature statistical model of GI and tidal downsizing as a mode of star and planet formation. We show results from four runs of the population synthesis model, varying the opacity law and the strength of migration, as well as investigating the effect of disc truncation during the fragmentation process. Our early results suggest that GI plus tidal downsizing is not the principal mode of planet formation, but remains an excellent means of forming gas giant planets, brown dwarfs and low mass stars at large semimajor axes. (Abridged)

With the detection of extrasolar moons (exomoons) on the horizon, it is important to consider their potential for habitability. If we consider the circumstellar Habitable Zone (HZ, often described in terms of planet semi-major axis and orbital eccentricity), we can ask, "How does the HZ for an Earth-like exomoon differ from the HZ for an Earth-like exoplanet?" For the first time, we use 1D latitudinal energy balance modelling to address this question. The model places an Earth-like exomoon in orbit around a Jupiter mass planet, which in turn orbits a Sun-like star. The exomoon's surface temperature is evolved under the action of stellar insolation, atmospheric cooling, heat diffusion, eclipses and tidal heating. We use this model to carry out two separate investigations. In the first, four test cases are run to investigate in detail the dependence of the exomoon climate on orbital direction, orbital inclination, and on the frequency of stellar eclipse by the host planet. We find that lunar orbits which are retrograde to the planetary orbit exhibit greater climate variations than prograde orbits, with global mean temperatures around 0.1 K warmer due to the geometry of eclipses. If eclipses become frequent relative to the atmospheric thermal inertia timescale, climate oscillations become extremely small. In the second investigation, we carry out an extensive parameter study, running the model many times to study the habitability of the exomoon in the 4-dimensional space composed of the planet semi-major axis and eccentricity, and the moon semi-major axis and eccentricity. We find that for zero moon eccentricity, frequent eclipses allow the moon to remain habitable in regions of high planet eccentricity, but tidal heating severely constrains habitability in the limit of high moon eccentricity, making the habitable zone a sensitive function of moon semi-major axis.

When a self-gravitating disc is subject to irradiation, its propensity to fragmentation will be affected. The strength of self-gravitating disc stresses is expected to dictate disc fragmentation: as the strength of these torques typically decrease with increasing sound speed, it is reasonable to assume, to first-order, that disc fragmentation is suppressed when compared to the non-irradiated case, although previous work has shown that the details are complicated by the source of the irradiation. We expand on previous analysis of the Jeans mass inside spiral structures in self-gravitating discs, incorporating the effects of stellar irradiation and background irradiation. If irradiation is present, fragmentation is suppressed for marginally unstable discs at low accretion rates (compared to the no-irradiation case), but these lower accretion rates correspond to higher mass discs. Fragmentation can still occur for high accretion rates, but is consequently suppressed at lower disc surface densities, and the subsequent Jeans mass is boosted. These results further bolster the consensus that, without subsequent fragment disruption or mass loss, the gravitational instability is more likely to form brown dwarfs and low-mass stars than gas giant planets.

The travel time required for one civilisation to explore the Milky Way using probes is a crucial component of Fermi's Paradox. Previous attempts to estimate this travel time have assumed that the probe's motion is simple, moving at a constant maximum velocity, with powered flight producing the necessary change in velocity required at each star to complete its chosen trajectory. This approach ignores lessons learned from interplanetary exploration, where orbital slingshot maneouvres can provide significant velocity boosts at little to no fuel cost. It is plausible that any attempt to explore the Galaxy would utilise such economising techniques, despite there being an upper limit to these velocity boosts, related to the escape velocity of the object being used to provide the slingshot. In order to investigate the effects of these techniques, we present multiple realisations of single probes exploring a small patch of the Milky Way. We investigate 3 separate scenarios, studying the slingshot effect on trajectories defined by simple heuristics. These scenarios are: i) standard powered flight to the nearest unvisited star without using slingshot techniques; ii) flight to the nearest unvisited star using slingshot techniques, and iii) flight to the next unvisited star which provides the maximal velocity boost under a slingshot trajectory. We find that adding slingshot velocity boosts can decrease the travel time by up to two orders of magnitude over simple powered flight. In the third case, selecting a route which maximises velocity boosts also reduces the travel time relative to powered flight, but by a much reduced factor. From these simulations, we suggest that adding realistic probe trajectories tends to strengthen Fermi's Paradox.

We use smoothed particle hydrodynamics simulations of massive protostellar discs to investigate the predicted broadening of molecular lines from discs in which self-gravity is the dominant source of angular momentum transport. The simulations include radiative transfer, and span a range of disc-to-star mass ratios between 0.25 and 1.5. Subtracting off the mean azimuthal flow velocity, we compute the distribution of the in-plane and perpendicular peculiar velocity due to large scale structure and turbulence induced by self-gravity. For the lower mass discs, we show that the characteristic peculiar velocities scale with the square root of the effective turbulent viscosity parameter, as expected from local turbulent-disc theory. The derived velocities are anisotropic, with substantially larger in-plane than perpendicular values. As the disc mass is increased, the validity of the locally determined turbulence approximation breaks down, and this is accompanied by anomalously large in-plane broadening. There is also a high variance due to the importance of low-m spiral modes. For low-mass discs, the magnitude of in-plane broadening is, to leading order, equal to the predictions from local disc theory and cannot constrain the source of turbulence. However, combining our results with prior evaluations of turbulent broadening expected in discs where the magnetorotational instability (MRI) is active, we argue that self-gravity may be distinguishable from the MRI in these systems if it is possible to measure the anisotropy of the peculiar velocity field with disc inclination. Furthermore, for large mass discs, the dominant contribution of large-scale modes is a distinguishing characteristic of self-gravitating turbulence versus MRI driven turbulence.

The Alpha Centauri AB system is an attractive one for radial velocity observations to detect potential exoplanets. The high metallicity of both Alpha Centauri A and B suggest that they could have possessed circumstellar discs capable of forming planets. As the closest star system to the Sun, with well over a century of accurate astrometric measurements (and Alpha Centauri B exhibiting low chromospheric activity) high precision surveys of Alpha Centauri B's potential exoplanetary system are possible with relatively cheap instrumentation. Authors studying habitability in this system typically adopt habitable zones (HZs) based on global radiative balance models that neglect the radiative perturbations of Alpha Centauri A. We investigate the habitability of planets around Alpha Centauri B using 1D latitudinal energy balance models (LEBMs), which fully incorporate the presence of Alpha Centauri A as a means of astronomically forcing terrestrial planet climates. We find that the extent of the HZ is relatively unchanged by the presence of Alpha Centauri A, but there are variations in fractional habitability for planets orbiting at the boundaries of the zone due to Alpha Centauri A, even in the case of zero eccentricity. Temperature oscillations of a few K can be observed at all planetary orbits, the strength of which varies with the planet's ocean fraction and obliquity.

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

We attempt to verify recent claims (made using semi-analytic models) that for the collapse of spherical homogeneous molecular clouds, fragmentation of the self-gravitating disc that subsequently forms can be predicted using the cloud's initial angular momentum alone. In effect, this condition is equivalent to requiring the resulting disc be sufficiently extended in order to fragment, in line with studies of isolated discs. We use smoothed particle hydrodynamics with hybrid radiative transfer to investigate this claim, confirming that in general, homogeneous spherical molecular clouds will produce fragmenting self-gravitating discs if the ratio of rotational kinetic energy to gravitational potential energy is greater than ~ 5e-3, where this result is relatively insensitive to the initial thermal energy. This condition begins to fail at higher cloud masses, suggesting that sufficient mass at large radii governs fragmentation. While these results are based on highly idealised initial conditions, and may not hold if the disc's accretion from the surrounding envelope is sufficiently asymmetric, or if the density structure is perturbed, they provide a sensible lower limit for the minimum angular momentum required to fragment a disc in the absence of significant external turbulence.

As a formation route for objects such as giant planets and low-mass stars in protostellar discs (as well as stars in AGN discs), theories of self-gravitating disc fragmentation need to be able to predict the initial masses of fragments. We describe a means by which the local Jeans mass inside the spiral structure of a self-gravitating disc can be estimated. If such a self-gravitating disc satisfies the criteria for disc fragmentation, this estimate provides a lower limit for the initial mass of any fragments formed. We apply this approach to a series of self-gravitating protostellar disc models, to map out the typical masses of fragments produced by this formation mode. We find a minimum fragment mass of around 3 Jupiter masses, which is insensitive to the stellar mass, and that - within the parameter space surveyed - fragments with masses between 10 and 20 Jupiter masses are the most common. We also describe how the Jeans mass allows us to derive a more general criterion for disc fragmentation, which accounts for the processes of viscous heating, radiative cooling, accretion and the disc's thermal history. We demonstrate how such a criterion can be determined, and show that in limiting cases it recovers several fragmentation criteria that have been posited in the past, including the minimum cooling time/maximum stress criterion.

The Zoo Hypothesis posits that we have not detected extraterrestrial intelligences (ETIs) because they deliberately prevent us from detecting them. While a valid solution to Fermi's Paradox, it is not particularly amenable to rigorous scientific analysis, as it implicitly assumes a great deal about the sociological structure of a plurality of civilisations. Any attempt to assess its worth must begin with its most basic assumption - that ETIs share a uniformity of motive in shielding Earth from extraterrestrial contact. This motive is often presumed to be generated by the influence of the first civilisation to arrive in the Galaxy. I show that recent work on inter-arrival time analysis, while necessary, is insufficient to assess the validity of the Zoo Hypothesis (and its related variants). The finite speed of light prevents an early civilisation from exerting immediate cultural influence over a later civilisation if they are sufficiently distant. I show that if civilisation arrival times and spatial locations are completely uncorrelated, this strictly prevents the establishment of total hegemony throughout the Galaxy. I finish by presenting similar results derived from more realistic Monte Carlo Realisation simulations (where arrival time and spatial locations are partially correlated). These also show that total hegemony is typically broken, even when the total population of civilisations remains low. In the terminology of previous studies of solutions to Fermi's Paradox, this confirms the Zoo Hypothesis as a "soft" solution. However, an important question to be resolved by future work is the extent to which many separate hegemonies are established, and to what extent this affects the Zoo Hypothesis.

The development of civilisations like ours into spacefaring, multi-planet entities requires significant raw materials to construct vehicles and habitats. Interplanetary debris, including asteroids and comets, may provide such a source of raw materials. In this article we present the hypothesis that extraterrestrial intelligences (ETIs) engaged in asteroid mining may be detectable from Earth. Considering the detected disc of debris around Vega as a template, we explore the observational signatures of targeted asteroid mining (TAM), such as unexplained deficits in chemical species, changes in the size distribution of debris and other thermal signatures which may be detectable in the spectral energy distribution (SED) of a debris disc. We find that individual observational signatures of asteroid mining can be explained by natural phenomena, and as such they cannot provide conclusive detections of ETIs. But, it may be the case that several signatures appearing in the same system will prove harder to model without extraterrestrial involvement. Therefore signatures of TAM are not detections of ETI in their own right, but as part of "piggy-back" studies carried out in tandem with conventional debris disc research, they could provide a means of identifying unusual candidate systems for further study using other SETI techniques.

Semi-analytic models of self-gravitating discs often approximate the angular momentum transport generated by the gravitational instability using the phenomenology of viscosity. This allows the employment of the standard viscous evolution equations, and gives promising results. It is, however, still not clear when such an approximation is appropriate. This paper tests this approximation using high resolution 3D smoothed particle hydrodynamics (SPH) simulations of self-gravitating protostellar discs with radiative transfer. The nature of angular momentum transport associated with the gravitational instability is characterised as a function of both the stellar mass and the disc-to-star mass ratio. The effective viscosity is calculated from the Reynolds and gravitational stresses in the disc. This is then compared to what would be expected if the effective viscosity were determined by assuming local thermodynamic equilibrium or, equivalently, that the local dissipation rate matches the local cooling rate. In general, all the discs considered here settle into a self-regulated state where the heating generated by the gravitational instability is modulated by the local radiative cooling. It is found that low-mass discs can indeed be represented by a local "alpha-parametrisation", provided that the disc aspect ratio is small (H/R < 0.1) which is generally the case when the disc-to-star mass ratio q <0.5. However, this result does not extend to discs with masses approaching that of the central object. These are subject to transient burst events and global wave transport, and the effective viscosity is not well modelled by assuming local thermodynamic equilibrium. In spite of these effects, it is shown that massive (compact) discs can remain stable and not fragment, evolving rapidly to reduce their disc-to-star mass ratios through stellar accretion and radial spreading.

The Square Kilometre Array (SKA) will operate in frequency ranges often used by military radar and other communications technology. It has been shown that if Extraterrestrial Intelligences (ETIs) communicate using similar technology, then the SKA should be able to detect such transmissions up to distances of ~100 pc (~300 light years) from Earth. However, Mankind has greatly improved its communications technology over the last century, dramatically reducing signal leakage and making the Earth "radio quiet". If ETIs follow the same pattern as the human race, will we be able to detect their signal leakage before they become radio quiet? We investigate this question using Monte Carlo Realisation techniques to simulate the growth and evolution of intelligent life in the Galaxy. We show that if civilisations are "human" in nature (i.e. they are only "radio loud" for ~100 years, and can only detect each other with an SKA-like instrument out to 100 pc, within a maximum communication time of 100 years), then the probability for such civilisations accidentally detecting each other is low (~10^{-7}), much lower than if other, dedicated communication techniques are permissible (e.g. optical SETI or neutrino communication).

An algorithm for creating synthetic telescope images of Smoothed Particle Hydrodynamics (SPH) density fields is presented, which utilises the adaptive nature of the SPH formalism in full. The imaging process uses Monte Carlo Radiative Transfer (MCRT) methods to model the scattering and absorption of photon packets in the density field, which then exit the system and are captured on a pixelated image plane, creating a 2D image (or a 3D datacube, if the photons are also binned by their wavelength). The algorithm is implemented on the density field directly: no gridding of the field is required, allowing the density field to be described to an identical level of accuracy as the simulations that generated it. Some applications of the method to star and planet formation simulations are presented to illustrate the advantages of this new technique, and suggestions as to how this framework could support a Radiative Equilibrium algorithm are also given as an indication for future development.

We begin with the premise that the law of entropy could prove to be fundamental for the evolution of intelligent life and the advent of technological civilization. Building on recent theoretical results, we combine a modern approach to evolutionary theory with Monte Carlo Realization Techniques. A numerical test for a proposed significance of the law of entropy within the evolution of intelligent species is performed and results are compared with a neutral test hypothesis. Some clarifying aspects on the emergence of intelligent species arise and are discussed in the framework of contemporary astrobiology.

The Search for Extraterrestrial Intelligence (SETI) has thus far failed to provide a convincing detection of intelligent life. In the wake of this null signal, many "contact pessimistic" hypotheses have been formulated, the most famous of which is the Rare Earth Hypothesis. It postulates that although terrestrial planets may be common, the exact environmental conditions that Earth enjoys are rare, perhaps unique. As a result, simple microbial life may be common, but complex metazoans (and hence intelligence) will be rare. This paper uses Monte Carlo Realisation Techniques to investigate the Rare Earth Hypothesis, in particular the environmental criteria considered imperative to the existence of intelligence on Earth. By comparing with a less restrictive, more optimistic hypothesis, the data indicates that if the Rare Earth hypothesis is correct, intelligent civilisation will indeed be relatively rare. Studying the separations of pairs of civilisations shows that most intelligent civilisation pairs (ICPs) are unconnected: that is, they will not be able to exchange signals at lightspeed in the limited time that both are extant. However, the few ICPs that are connected are strongly connected, being able to participate in numerous exchanges of signals. This may provide encouragement for SETI researchers: although the Rare Earth Hypothesis is in general a contact-pessimistic hypothesis, it may be a "soft" or "exclusive" hypothesis, i.e. it may contain facets that are latently contact-optimistic.

Young stellar systems are known to undergo outbursts, where the star experiences an increased accretion rate, and the system's luminosity increases accordingly. The archetype is the FU Orionis (FU Ori) outburst, where the accretion rate can increase by three orders of magnitude (and the brightness of the system by five magnitudes). The cause appears to be instability in the circumstellar disc, but there is currently some debate as to the nature of this instability (e.g. thermal, gravitational, magneto-rotational). This paper details high resolution Smoothed Particle Hydrodynamics (SPH) simulations that were carried out to investigate the influence of stellar encounters on disc dynamics. Star-star encounters (where the primary has a self-gravitating, marginally stable protostellar disc) were simulated with various orbital parameters to investigate the resulting disc structure and dynamics. Crucially, the simulations include the effects of radiative transfer to realistically model the resulting thermodynamics. Our results show that the accretion history and luminosity of the system during the encounter displays many of the features of outburst phenomena. In particular, the magnitudes and decay times seen are comparable to those of FU Ori. There are two caveats to this assertion: the first is that these events are not expected to occur frequently enough to explain all FU Ori or EX Lupi; the second is that the inner discs of these simulations are subject to numerical viscosity, which will act to reduce the accretion rate (although it has less of an effect on the total mass accreted). In short, these results cannot rule out binary interactions as a potential source of some FU Ori-esque outbursts.

An interaction between a star-disc system and another star will perturb the disc, possibly resulting in a significant modification of the disc structure and its properties. It is still unclear if such an encounter can trigger fragmentation of the disc to form brown dwarfs or gas giant planets. This paper details high resolution Smoothed Particle Hydrodynamics (SPH) simulations investigating the influence of stellar encounters on disc dynamics. Star-star encounters (where the primary has a self-gravitating, marginally stable protostellar disc, and the secondary has no disc) were simulated with various orbital parameters to investigate the resulting disc structure and dynamics. This work is the first of its kind to incorporate realistic radiative transfer techniques to realistically model the resulting thermodynamics. The results suggest that the effect of stellar encounters is to prohibit fragmentation - compressive and shock heating stabilises the disc, and the radiative cooling is insufficient to trigger gravitational instability. The encounter strips the outer regions of the disc (either through tidal tails or by capture of matter to form a disc around the secondary), which triggers a readjustment of the primary disc to a steeper surface density profile (and a flatter Toomre Q profile). The disc around the secondary plays a role in the potential capture of the secondary to form a binary. However, this applies only to orbits that are parabolic - hyperbolic encounters do not form a secondary disc, and are not captured.

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods - the polytropic cooling approximation as devised by Stamatellos et al (2007), and flux limited diffusion (e.g. Mayer et al 2007). This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of molecular hydrogen dissociation, atomic hydrogen ionisation, He0 and He+ ionisation, ice evaporation, dust sublimation, molecular absorption, bound-free and free-free transitions and electron scattering. The method is tested in several scenarios, including: (1) the evolution of a 0.07 solar mass protoplanetary disc surrounding a 0.5 solar mass star; (2) the collapse of a 1 solar mass protostellar cloud, and (3) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.

The search for extraterrestrial intelligence (SETI) has been heavily influenced by solutions to the Drake Equation, which returns an integer value for the number of communicating civilisations resident in the Milky Way, and by the Fermi Paradox, glibly stated as: "If they are there, where are they?". Both rely on using average values of key parameters, such as the mean signal lifetime of a communicating civilisation. A more accurate answer must take into account the distribution of stellar, planetary and biological attributes in the galaxy, as well as the stochastic nature of evolution itself. This paper outlines a method of Monte Carlo realisation which does this, and hence allows an estimation of the distribution of key parameters in SETI, as well as allowing a quantification of their errors (and the level of ignorance therein). Furthermore, it provides a means for competing theories of life and intelligence to be compared quantitatively.

We present a new method of incorporating radiative transfer into Smoothed Particle Hydrodynamics (SPH). There have been many recent attempts at radiative transfer in SPH (Stamatellos et al 2005, 2005, Mayer et al 2007, Whitehouse and Bate 2006), however these are becoming increasingly complex, with some methods requiring the photosphere to be mapped (which is often of non-trivial geometric shape), and extra conditions to be applied there (matching atmospheres as in Cai et al (2008), or specifying cooling at the photosphere as in Mayer et al (2007)). The method of identifying the photosphere is usually a significant addition to the total simulation runtime, and often requires extra free parameters, the changing of which will affect the final results. Our method is not affected by such concerns, as the photosphere is constructed implicitly by the algorithm without the need for extra free parameters. The algorithm used is a synergy of two current formalisms for radiative effects: a) the polytropic cooling formalism proposed by Stamatellos et al (2007), and b) flux-limited diffusion, used by many authors to simulate radiation transport in the optically thick regime (e.g. Mayer et al 2007). We present several tests of this method: (1) The evolution of a 0.07 solar mass protoplanetary disc around a 0.5 solarmass star (Pickett et al 2003, Mejia et al 2005, Boley et al 2006, Cai et al 2008); (2) The collapse of a non-rotating 1 solar mass molecular cloud (Masunaga & Inutsuka 2000, Stamatellos et al 2007); (3) The thermal relaxation of temperature fluctuations in an static homogeneous sphere (Masunaga et al 1998, Spiegel 1957, Stamatellos et al 2007)