Source author record

Duncan Forgan

Duncan Forgan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2020arXiv

Massive discs around low-mass stars

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.

preprint2018arXiv

The habitable zone for Earthlike exomoons orbiting Kepler-1625b

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.

preprint2016arXiv

#FoundThem - 21st Century Pre-Search and Post-Detection SETI Protocols for Social and Digital Media

The transmission of news stories in global culture has changed fundamentally in the last three decades. The general public are alerted to breaking stories on increasingly rapid timescales, and the discussion/distortion of facts by writers, bloggers, commenters and Internet users can also be extremely fast. The narrative of a news item no longer belongs to a small cadre of conventional media outlets, but is instead synthesised to some level by the public as they select where and how they consume news. The IAA Search for Extraterrestrial Intelligence (SETI) post-detection protocols, initially drafted in 1989 and updated in 2010, were written to guide SETI scientists in the event of detecting evidence of extraterrestrial intelligence, but do not give guidance as to how scientists should prepare to navigate this media maelstrom. The protocols assume communications channels between scientists and the public still resemble those of 1989, which were specifically one-way with a narrative controlled by a select few media outlets. Modern SETI researchers must consider this modern paradigm for consumption of news by the public, using social media and other non-traditional outlets, when planning and executing searches for extraterrestrial intelligence. We propose additions to the post-detection protocols as they pertain to the use of internet and social media, as well as pre-search protocols. It is our belief that such protocols are necessary if there is to be a well-informed, sane global conversation amongst the world's citizens following the discovery of intelligent life beyond the Earth.

preprint2016arXiv

Daughter Fragmentation is Unlikely To Occur in Self-Gravitating Circumstellar Discs

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.

preprint2016arXiv

Directly observing continuum emission from self-gravitating spiral waves

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.

preprint2016arXiv

Exomoon Climate Models with the Carbonate-Silicate Cycle and Viscoelastic Tidal Heating

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.

preprint2016arXiv

Tensor classification of structure in smoothed particle hydrodynamics density fields

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.

preprint2015arXiv

Disc fragmentation rarely forms planetary-mass objects

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.

preprint2014arXiv

The Dynamical Fate of Self-Gravitating Disc Fragments After Tidal Downsizing

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.

preprint2014arXiv

The Effect of Planetary Illumination on Climate Modelling of Earthlike Exomoons

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.

preprint2013arXiv

A Self-Gravitating Disc Around L1527 IRS?

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.

preprint2013arXiv

Assessing Circumbinary Habitable Zones using Latitudinal Energy Balance Modelling

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.

preprint2013arXiv

Dynamical Effects on the Habitable Zone for Earth-like Exomoons

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.

preprint2013arXiv

Slingshot Dynamics for Self Replicating Probes and the Effect on Exploration Timescales

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.

preprint2013arXiv

The Effect of Irradiation on the Jeans Mass in Fragmenting Self-Gravitating Protostellar Discs

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.

preprint2013arXiv

The Hunt for Exomoons with Kepler (HEK): III. The First Search for an Exomoon around a Habitable-Zone Planet

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.

preprint2013arXiv

Towards a Population Synthesis Model of Objects formed by Self-Gravitating Disc Fragmentation and Tidal Downsizing

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)

preprint2013arXiv

Virulence as a Model for Interplanetary and Interstellar Colonisation - Parasitism or Mutualism

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.

preprint2012arXiv

Oscillations in the Habitable Zone around Alpha Centauri B

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.

preprint2011arXiv

A Lower Angular Momentum Limit for Self-Gravitating Protostellar Disc Fragmentation

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.

preprint2011arXiv

Extrasolar Asteroid Mining as Forensic Evidence for Extraterrestrial Intelligence

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.

preprint2011arXiv

The Jeans Mass as a Fundamental Measure of Self-Gravitating Disc Fragmentation and Initial Fragment Mass

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.

preprint2010arXiv

Native Synthetic Imaging of Smoothed Particle Hydrodynamics density fields using gridless Monte Carlo Radiative Transfer

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.

preprint2010arXiv

The Nature of Angular Momentum Transport in Radiative Self-Gravitating Protostellar Discs

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.

preprint2009arXiv

Stellar Encounters in the Context of Outburst Phenomena

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.

preprint2009arXiv

Stellar Encounters: A Stimulus for Disc Fragmentation?

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.

preprint2008arXiv

A Numerical Testbed for Hypotheses of Extraterrestrial Life and Intelligence

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.