Showing 1–50 of 72 results for author: Forgan, D

The young protostellar disk in IRAS16293-2422 B is hot and shows signatures of gravitational instability

Authors: Joaquin Zamponi, María José Maureira, Bo Zhao, Hauyu Baobab Liu, John D. Ilee, Duncan Forgan, Paola Caselli

Abstract: Deeply embedded protostars are actively fed from their surrounding envelopes through their protostellar disk. The physical structure of such early disks might be different from that of more evolved sources due to the active accretion. We present 1.3 and 3\,mm ALMA continuum observations at resolutions of 6.5\,au and 12\,au respectively, towards the Class 0 source IRAS 16293-2422 B. The resolved br… ▽ More Deeply embedded protostars are actively fed from their surrounding envelopes through their protostellar disk. The physical structure of such early disks might be different from that of more evolved sources due to the active accretion. We present 1.3 and 3\,mm ALMA continuum observations at resolutions of 6.5\,au and 12\,au respectively, towards the Class 0 source IRAS 16293-2422 B. The resolved brightness temperatures appear remarkably high, with $T_{\rm b} >$ 100\,K within $\sim$30\,au and $T_{\rm b}$ peak over 400\,K at 3\,mm. Both wavelengths show a lopsided emission with a spectral index reaching values less than 2 in the central $\sim$ 20\,au region. We compare these observations with a series of radiative transfer calculations and synthetic observations of magnetohydrodynamic and radiation hydrodynamic protostellar disk models formed after the collapse of a dense core. Based on our results, we argue that the gas kinematics within the disk may play a more significant role in heating the disk than the protostellar radiation. In particular, our radiation hydrodynamic simulation of disk formation, including heating sources associated with gravitational instabilities, is able to generate the temperatures necessary to explain the high fluxes observed in IRAS 16293B. Besides, the low spectral index values are naturally reproduced by the high optical depth and high inner temperatures of the protostellar disk models. The high temperatures in IRAS 16293B imply that volatile species are mostly in the gas phase, suggesting that a self-gravitating disk could be at the origin of a hot corino. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 16 pages, 13 figures. Accepted for publication in MNRAS

The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300 au with SPHERE

Authors: A. Vigan, C. Fontanive, M. Meyer, B. Biller, M. Bonavita, M. Feldt, S. Desidera, G. -D. Marleau, A. Emsenhuber, R. Galicher, K. Rice, D. Forgan, C. Mordasini, R. Gratton, H. Le Coroller, A. -L. Maire, F. Cantalloube, G. Chauvin, A. Cheetham, J. Hagelberg, A. -M. Lagrange, M. Langlois, M. Bonnefoy, J. -L. Beuzit, A. Boccaletti , et al. (86 additional authors not shown)

Abstract: The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses betwee… ▽ More The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 24 pages, 14 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 651, A72 (2021)

Massive discs around low-mass stars

Authors: Thomas J. Haworth, James Cadman, Farzana Meru, Cassandra Hall, Emma Albertini, Duncan Forgan, Ken Rice, James E. Owen

Abstract: 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. T… ▽ More 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. △ Less

Submitted 27 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, Accepted for publication MNRAS

Pebble accretion in self-gravitating protostellar discs

Authors: Duncan H Forgan

Abstract: 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… ▽ More 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. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: 10 pages, 7 figures, accepted for publication in MNRAS

G11.92-0.61 MM 1: A fragmented Keplerian disk surrounding a proto-O star

Authors: J. D. Ilee, C. J. Cyganowski, C. L. Brogan, T. R. Hunter, D. H. Forgan, T. J. Haworth, C. J. Clarke, T. J. Harries

Abstract: 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 M… ▽ More 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. △ Less

Submitted 13 November, 2018; originally announced November 2018.

Comments: 8 pages, 3 figures, accepted for publication in ApJ Letters

Climate modelling of hypothetical moon-moons in the Kepler-1625b system

Authors: Duncan Forgan

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2018; originally announced October 2018.

Comments: Published as an AAS Research Note: 3 pages, 1 figure

Observing substructure in circumstellar discs around massive young stellar objects

Authors: Marija R. Jankovic, Thomas J. Haworth, John D. Ilee, Duncan H. Forgan, Claudia J. Cyganowski, Catherine Walsh, Crystal L. Brogan, Todd R. Hunter, Subhanjoy Mohanty

Abstract: 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 inhere… ▽ More 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. △ Less

Submitted 6 November, 2018; v1 submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS, models available online (DOI 10.5281/zenodo.1408072)

The habitable zone for Earthlike exomoons orbiting Kepler-1625b

Authors: Duncan Forgan

Abstract: 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… ▽ More 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. △ Less

Submitted 28 November, 2018; v1 submitted 4 October, 2018; originally announced October 2018.

Comments: 17 pages, 5 figures, accepted to IJA. arXiv admin note: text overlap with arXiv:1601.00789, arXiv:1404.4493

The fate of formamide in a fragmenting protoplanetary disc

Authors: David Quénard, John D. Ilee, Izaskun Jiménez-Serra, Duncan H. Forgan, Cassandra Hall, Ken Rice

Abstract: 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 fragmenti… ▽ More 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. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: 10 pages, 5 figures, accepted for publication in ApJ

Clumpy shocks as the driver of velocity dispersion in molecular clouds: the effects of self-gravity and magnetic fields

Authors: Duncan H. Forgan, Ian A. Bonnell

Abstract: 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 ($σ\propto L^Γ$ with… ▽ More 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 ($σ\propto L^Γ$ with $Γ\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 $σ$ and $L$ throughout, with indices ranging from $Γ=0.3-1.2$. These results are relatively insensitive to the strength and geometry of magnetic fields, provided that the shock is relatively strong. $Γ$ 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 $σ-L$ relation is in clumpy shocks, deviations from the standard Larson relation constrain the strength and behaviour of shocks in spiral galaxies. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Comments: 12 pages, 22 figures, accepted to MNRAS

Are Elias 2-27's spiral arms driven by self-gravity, or by a companion? A comparative spiral morphology study

Authors: Duncan H. Forgan, John D. Ilee, Farzana Meru

Abstract: 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 cas… ▽ More 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. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: 9 pages, 3 figures, accepted for publication in ApJL

Is the spiral morphology of the Elias 2-27 circumstellar disc due to gravitational instability?

Authors: Cassandra Hall, Ken Rice, Giovanni Dipierro, Duncan Forgan, Tim Harries, Richard Alexander

Abstract: 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 propertie… ▽ More 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. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: 12 pages, 5 figures

Classifying and modelling spiral structures in hydrodynamic simulations of astrophysical discs

Authors: D. H. Forgan, F. G. Ramón-Fox, I. A. Bonnell

Abstract: 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,… ▽ More 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. △ Less

Submitted 5 February, 2018; originally announced February 2018.

Comments: 13 pages, 19 figures, accepted to MNRAS

The Astrobiology of the Anthropocene

Authors: Jacob Haqq-Misra, Sanjoy Som, Brendan Mullan, Rafael Loureiro, Edward Schwieterman, Lauren Seyler, Haritina Mogosanu, Adam Frank, Eric Wolf, Duncan Forgan, Charles Cockell, Woodruff Sullivan

Abstract: 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… ▽ More 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. △ Less

Submitted 6 January, 2018; v1 submitted 29 December, 2017; originally announced January 2018.

Comments: A white paper on "Astrobiology Science Strategy" submitted to the NAS

Interstellar communication. VI. Searching X-ray spectra for narrowband communication

Authors: Michael Hippke, Duncan H. Forgan

Abstract: We have previously argued that targeted interstellar communication has a physical optimum at narrowband X-ray wavelengths $λ\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 th… ▽ More We have previously argued that targeted interstellar communication has a physical optimum at narrowband X-ray wavelengths $λ\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. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 3 pages, 2 figures, 1 table

Photogravimagnetic assists of light sails: a mixed blessing for Breakthrough Starshot?

Authors: Duncan H. Forgan, René Heller, Michael Hippke

Abstract: 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 $α$ Centauri A, a… ▽ More 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 $α$ 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 $μ$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 $α$ 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. △ Less

Submitted 30 October, 2017; originally announced November 2017.

Comments: 10 pages, 8 figures, accepted for publication in MNRAS

Interstellar communication. III. Optimal frequency to maximize data rate

Authors: Michael Hippke, Duncan H. Forgan

Abstract: 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 $λ<300\,$nm (UV). When t… ▽ More 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 $λ<300\,$nm (UV). When this technological constraint is dropped, a physical bound is found at $λ\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. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Towards a population synthesis model of self-gravitating disc fragmentation and tidal downsizing II: The effect of fragment-fragment interactions

Authors: D. H. Forgan, C. Hall, F. Meru, W. K. M. Rice

Abstract: 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 f… ▽ More 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. △ Less

Submitted 7 November, 2017; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 14 pages, 17 figures, accepted for publication in MNRAS

Can planet formation resolve the dust budget crisis in high redshift galaxies?

Authors: D. H. Forgan, K. Rowlands, H. L. Gomez, E. L. Gomez, S. P. Schofield, L. Dunne, S. Maddox

Abstract: 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… ▽ More 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. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Comments: 9 pages, 11 figures, accepted for publication in MNRAS

The chemistry of protoplanetary fragments formed via gravitational instabilities

Authors: J. D. Ilee, D. H. Forgan, M. G. Evans, C. Hall, R. Booth, C. J. Clarke, W. K. M. Rice, A. C. Boley, P. Caselli, T. W. Hartquist, J. M. C. Rawlings

Abstract: 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… ▽ More 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. △ Less

Submitted 5 August, 2017; originally announced August 2017.

Comments: 17 pages, 7 figures, accepted for publication in MNRAS

Exoplanet Transits as the Foundation of an Interstellar Communications Network

Authors: Duncan H Forgan

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: 17 pages, 21 figures, accepted for publication in the International Journal of Astrobiology

Identifying and Analysing Protostellar Disc Fragments in Smoothed Particle Hydrodynamics Simulations

Authors: Cassandra Hall, Duncan Forgan, Ken Rice

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: 22 pages, 22 figures

On The Feasibility of Exomoon Detection Via Exoplanet Phase Curve Spectral Contrast

Authors: Duncan H. Forgan

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 May, 2017; originally announced May 2017.

Comments: 12 pages, 26 figures, accepted for publication by MNRAS

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. IV. Gravitational instability rarely forms wide, giant planets

Authors: A. Vigan, M. Bonavita, B. Biller, D. Forgan, K. Rice, G. Chauvin, S. Desidera, J. -C. Meunier, P. Delorme, J. E. Schlieder, M. Bonnefoy, J. Carson, E. Covino, J. Hagelberg, T. Henning, M. Janson, A. -M. Lagrange, S. P. Quanz, A. Zurlo, J. -L. Beuzit, A. Boccaletti, E. Buenzli, M. Feldt, J. H. V. Girard, R. Gratton , et al. (12 additional authors not shown)

Abstract: 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,… ▽ More 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. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 20 pages, 11 figures. Accepted for publication in A&A

Journal ref: A&A 603, A3 (2017)

Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics

Authors: Daniel J. Price, James Wurster, Terrence S. Tricco, Chris Nixon, Stéven Toupin, Alex Pettitt, Conrad Chan, Daniel Mentiplay, Guillaume Laibe, Simon Glover, Clare Dobbs, Rebecca Nealon, David Liptai, Hauke Worpel, Clément Bonnerot, Giovanni Dipierro, Giulia Ballabio, Enrico Ragusa, Christoph Federrath, Roberto Iaconi, Thomas Reichardt, Duncan Forgan, Mark Hutchison, Thomas Constantino, Ben Ayliffe , et al. (2 additional authors not shown)

Abstract: 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 bir… ▽ More 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. △ Less

Submitted 15 June, 2018; v1 submitted 13 February, 2017; originally announced February 2017.

Comments: 88 pages, 60 figures, accepted to PASA. Code available from https://phantomsph.bitbucket.io/

Can the removal of molecular cloud envelopes by external feedback affect the efficiency of star formation?

Authors: William Lucas, Ian Bonnell, Duncan Forgan

Abstract: 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 ho… ▽ More 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 $α$, 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. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

On the fragmentation boundary in magnetised self-gravitating discs

Authors: Duncan Forgan, Daniel J. Price, Ian Bonnell

Abstract: 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 amplitud… ▽ More 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. △ Less

Submitted 19 December, 2016; originally announced December 2016.

Comments: 12 pages, 22 figures, accepted for publication in MNRAS

Milankovitch Cycles of Terrestrial Planets in Binary Star Systems

Authors: Duncan H Forgan

Abstract: 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 resonanc… ▽ More 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. △ Less

Submitted 31 August, 2016; v1 submitted 19 August, 2016; originally announced August 2016.

Comments: 14 pages, 32 figures, accepted for publication in MNRAS (revised to add OBERON to acknowledgements)

G11.92-0.61 MM1: A Keplerian disc around a massive young proto-O star

Authors: J. D. Ilee, C. J. Cyganowski, P. Nazari, T. R. Hunter, C. L. Brogan, D. H. Forgan, Q. Zhang

Abstract: 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 ve… ▽ More 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ıı 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. △ Less

Submitted 19 August, 2016; originally announced August 2016.

Comments: 18 pages, 9 figures. Accepted to MNRAS

Self-gravitating disc candidates around massive young stars

Authors: D. H. Forgan, J. D. Ilee, C. J. Cyganowski, C. L. Brogan, T. R. Hunter

Abstract: 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. O… ▽ More 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. △ Less

Submitted 19 August, 2016; originally announced August 2016.

Comments: 9 pages, 20 figures, accepted for publication in MNRAS

Grand challenges in protoplanetary disc modelling

Authors: Thomas J. Haworth, John D. Ilee, Duncan H. Forgan, Stefano Facchini, Daniel J. Price, Community authors, Dominika M. Boneberg, Richard A. Booth, Cathie J. Clarke, Jean-François Gonzalez, Mark A. Hutchison, Inga Kamp, Guillaume Laibe, Wladimir Lyra, Farzana Meru, Subhanjoy Mohanty, Olja Panić, Ken Rice, Takeru Suzuki, Catherine Walsh, Peter Woitke

Abstract: 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 t… ▽ More 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. △ Less

Submitted 20 September, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: Accepted for publication in PASA. 25 pages

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

Authors: Duncan Forgan

Abstract: 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 expec… ▽ More 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. △ Less

Submitted 22 February, 2016; originally announced February 2016.

Comments: 7 pages, 4 figures, submitted to the Open Journal of Astrophysics

Directly observing continuum emission from self-gravitating spiral waves

Authors: Cassandra Hall, Duncan Forgan, Ken Rice, Tim J. Harries, Pamela D. Klaassen, Beth Biller

Abstract: 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 detec… ▽ More 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. △ Less

Submitted 3 February, 2016; originally announced February 2016.

Comments: 13 pages, 9 figures

Tensor classification of structure in smoothed particle hydrodynamics density fields

Authors: Duncan Forgan, Ian Bonnell, William Lucas, Ken Rice

Abstract: 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 hydrodyna… ▽ More 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. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 14 pages, 28 figures, accepted for publication in MNRAS

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

Authors: Duncan Forgan, Vera Dobos

Abstract: 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 in… ▽ More 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. △ Less

Submitted 5 January, 2016; originally announced January 2016.

Comments: 10 pages, 15 figures, accepted for publication in MNRAS

Evaluating Galactic Habitability Using High Resolution Cosmological Simulations of Galaxy Formation

Authors: Duncan Forgan, Pratika Dayal, Charles Cockell, Noam Libeskind

Abstract: 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 an… ▽ More 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. △ Less

Submitted 5 November, 2015; originally announced November 2015.

Comments: 24 pages, 33 figures, accepted for publication in the International Journal of Astrobiology

Disc fragmentation rarely forms planetary-mass objects

Authors: Ken Rice, Eric Lopez, Duncan Forgan, Beth Biller

Abstract: 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.… ▽ More 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. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 8 pages, 5 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Observational Signatures of Self-Destructive Civilisations

Authors: Adam Stevens, Duncan Forgan, Jack O'Malley-James

Abstract: 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). U… ▽ More 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. △ Less

Submitted 30 July, 2015; originally announced July 2015.

Comments: 35 pages, 2 figures, accepted for publication in the International Journal of Astrobiology

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

Authors: Duncan Forgan, Richard Parker, Ken Rice

Abstract: 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… ▽ More 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. △ Less

Submitted 26 November, 2014; originally announced November 2014.

Comments: 11 pages, 16 figures, accepted for publication in MNRAS

Can Collimated Extraterrestrial Signals be Intercepted?

Authors: Duncan H. Forgan

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 October, 2014; originally announced October 2014.

Comments: 12 pages, 7 figures, accepted for publication in JBIS

Surface Flux Patterns on Planets in Circumbinary Systems, and Potential for Photosynthesis

Authors: Duncan H. Forgan, Alexander Mead, Charles S. Cockell, John A. Raven

Abstract: 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… ▽ More 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. △ Less

Submitted 22 August, 2014; originally announced August 2014.

Comments: 26 pages, 43 figures, accepted for publication in the International Journal of Astrobiology

The Effect of Planetary Illumination on Climate Modelling of Earthlike Exomoons

Authors: Duncan Forgan, Vergil Yotov

Abstract: 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 a… ▽ More 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. △ Less

Submitted 17 April, 2014; originally announced April 2014.

Comments: 11 pages, 16 figures, accepted for publication in MNRAS

Photosynthetic Potential of Planets in 3:2 Spin Orbit Resonances

Authors: S. P. Brown, A. J. Mead, D. H. Forgan, J. A. Raven, C. S. Cockell

Abstract: 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 un… ▽ More 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. △ Less

Submitted 20 February, 2014; originally announced February 2014.

Comments: 18 pages, 17 figures, accepted for publication in the International Journal of Astrobiology

Triple trouble for XZ Tau: deep imaging with the Jansky Very Large Array

Authors: D. Forgan, R. J. Ivison, B. Sibthorpe, J. S. Greaves, E. Ibar

Abstract: 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 orbita… ▽ More 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. △ Less

Submitted 10 February, 2014; originally announced February 2014.

Comments: 4 pages, 4 figures, accepted for publication in MNRAS

Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity

Authors: W. K. M. Rice, S. -J. Paardekooper, D. H. Forgan, P. J. Armitage

Abstract: 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 long… ▽ More 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. △ Less

Submitted 13 January, 2014; v1 submitted 28 November, 2013; originally announced November 2013.

Comments: 11 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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

Authors: Jonathan Starling, Duncan Forgan

Abstract: 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… ▽ More 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. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 18 pages, 7 figures, published in the International Journal of Astrobiology

Assessing Circumbinary Habitable Zones using Latitudinal Energy Balance Modelling

Authors: Duncan Forgan

Abstract: 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 (f… ▽ More 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. △ Less

Submitted 14 October, 2013; originally announced October 2013.

Comments: 10 pages, 17 figures, accepted for publication in MNRAS

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

Authors: Arwen Nicholson, Duncan Forgan

Abstract: 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 whe… ▽ More 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. △ Less

Submitted 5 July, 2013; originally announced July 2013.

Comments: Accepted for publication in the International Journal of Astrobiology, 13 pages, 7 figures

On the Possibility of Detecting Class A Stellar Engines Using Exoplanet Transit Curves

Authors: Duncan H. Forgan

Abstract: 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 m… ▽ More 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. △ Less

Submitted 7 June, 2013; originally announced June 2013.

Comments: 20 pages, 10 figures, accepted for publication in the Journal of the British Interplanetary Society

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

Authors: David M. Kipping, Duncan Forgan, Joel Hartman, David Nesvorny, Gáspár Á. Bakos, Allan R. Schmitt, Lars A. Buchhave

Abstract: 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 avai… ▽ More 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. △ Less

Submitted 9 September, 2013; v1 submitted 6 June, 2013; originally announced June 2013.

Comments: 19 pages, 11 figures, 7 tables. Accepted in ApJ. Planet-moon transit animations available at https://www.cfa.harvard.edu/~dkipping/kepler22.html

Journal ref: Astrophys.J.777:134-150,2013