Showing 1–50 of 90 results for author: Parker, R J

No signature of the birth environment of exoplanets from their host stars' Mahalanobis phase space

Authors: George A. Blaylock-Squibbs, Richard J. Parker, Emma C. Daffern-Powell

Abstract: The architectures of extrasolar planetary systems often deviate considerably from the ``standard" model for planet formation, which is largely based on our own Solar System. In particular, gas giants on close orbits are not predicted by planet formation theory and so some process(es) are thought to move the planets closer to their host stars. Recent research has suggested that Hot Jupiter host sta… ▽ More The architectures of extrasolar planetary systems often deviate considerably from the ``standard" model for planet formation, which is largely based on our own Solar System. In particular, gas giants on close orbits are not predicted by planet formation theory and so some process(es) are thought to move the planets closer to their host stars. Recent research has suggested that Hot Jupiter host stars display a different phase space compared to stars that do not host Hot Jupiters. This has been attributed to these stars forming in star-forming regions of high stellar density, where dynamical interactions with passing stars have perturbed the planets. We test this hypothesis by quantifying the phase space of planet-hosting stars in dynamical N-body simulations of star-forming regions. We find that stars that retain their planets have a higher phase space than non-hosts, regardless of their initial physical density. This is because an imprint of the kinematic substructure from the regions birth is retained, as these stars have experienced fewer and less disruptive encounters than stars whose planets have been liberated and become free-floating. However, host stars whose planets remain bound but have had their orbits significantly altered by dynamical encounters are also primarily found in high phase space regimes. We therefore corroborate other research in this area which has suggested the high phase space of Hot Jupiter host stars is not caused by dynamical encounters or stellar clustering, but rather reflects an age bias in that these stars are (kinematically) younger than other exoplanet host stars. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 11 pages, 4 figures. Accepted for publication in ApJ

The dynamical evolution of star-forming regions measured with INDICATE

Authors: George A. Blaylock-Squibbs, Richard J. Parker

Abstract: Observations of star-forming regions provide snapshots in time of the star formation process, and can be compared with simulation data to constrain the initial conditions of star formation. In order to make robust inferences, different metrics must be used to quantify the spatial and kinematic distributions of stars. In this paper, we assess the suitability of the INDICATE (INdex to Define Inheren… ▽ More Observations of star-forming regions provide snapshots in time of the star formation process, and can be compared with simulation data to constrain the initial conditions of star formation. In order to make robust inferences, different metrics must be used to quantify the spatial and kinematic distributions of stars. In this paper, we assess the suitability of the INDICATE (INdex to Define Inherent Clustering And TEndencies) method as a diagnostic to infer the initial conditions of star-forming regions that subsequently undergo dynamical evolution. We use INDICATE to measure the degree of clustering in N-body simulations of the evolution of star-forming regions with different initial conditions. We find that the clustering of individual stars, as measured by INDICATE, becomes significantly higher in simulations with higher initial stellar densities, and is higher in subvirial star-forming regions where significant amounts of dynamical mixing has occurred. We then combine INDICATE with other methods that measure the mass segregation, relative stellar surface density ratio and the morphology (Q-parameter) of star-forming regions, and show that the diagnostic capability of INDICATE increases when combined with these other metrics. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 17 pages (including figures), 12 Figures. Accepted for publication in MNRAS

Devolatilization of extrasolar planetesimals by 60Fe and 26Al heating

Authors: Joseph W. Eatson, Tim Lichtenberg, Richard J. Parker, Taras V. Gerya

Abstract: Whilst the formation of Solar system planets is constrained by meteoritic evidence, the geophysical history of low-mass exoplanets is much less clear. The bulk composition and climate states of rocky exoplanets may vary significantly based on the composition and properties of the planetesimals they form from. An important factor influenced by planetesimal composition is water content, where the de… ▽ More Whilst the formation of Solar system planets is constrained by meteoritic evidence, the geophysical history of low-mass exoplanets is much less clear. The bulk composition and climate states of rocky exoplanets may vary significantly based on the composition and properties of the planetesimals they form from. An important factor influenced by planetesimal composition is water content, where the desiccation of accreting planetesimals impacts the final water content of the resultant planets. While the inner planets of the Solar system are comparatively water-poor, recent observational evidence from exoplanet bulk densities and planetary formation models suggest that rocky exoplanets engulfed by substantial layers of high-pressure ices or massive steam atmospheres could be widespread. Here we quantify variations in planetesimal desiccation due to potential fractionation of the two short-lived radioisotopes 26Al and 60Fe relevant for internal heating on planetary formation timescales. We focus on how order of magnitude variations in 60Fe can affect the water content of planetesimals, and how this may alter the formation of extrasolar ocean worlds. We find that heating by 26Al is the dominant cause of planetesimal heating in any Solar system analogue scenario, thus validating previous works focussing only on this radioisotope. However, 60Fe can become the primary heating source in the case of high levels of supernova enrichment in massive star-forming regions. These diverging scenarios can affect the formation pathways, bulk volatile budget, and climate diversity of low-mass exoplanets. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 12 pages, 13 figures

JWST/NIRSpec Observations of the Coldest Known Brown Dwarf

Authors: K. L. Luhman, P. Tremblin, C. Alves de Oliveira, S. M. Birkmann, I. Baraffe, G. Chabrier, E. Manjavacas, R. J. Parker, J. Valenti

Abstract: We present 1-5um spectroscopy of the coldest known brown dwarf, WISE J085510.83-071442.5 (WISE 0855), performed with the Near-Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST). NIRSpec has dramatically improved the measurement of spectral energy distribution of WISE 0855 in terms of wavelength coverage, signal-to-noise ratios, and spectral resolution. We have performed… ▽ More We present 1-5um spectroscopy of the coldest known brown dwarf, WISE J085510.83-071442.5 (WISE 0855), performed with the Near-Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST). NIRSpec has dramatically improved the measurement of spectral energy distribution of WISE 0855 in terms of wavelength coverage, signal-to-noise ratios, and spectral resolution. We have performed preliminary modeling of the NIRSpec data using the ATMO 2020 models of cloudless atmospheres, arriving at a best fitting model that has T_eff=285 K. That temperature is ~20 K higher than the value derived by combining our luminosity estimate with evolutionary models (i.e., the radius in the model fit to the SED is somewhat smaller than expected from evolutionary models). Through comparisons to the model spectra, we detect absorption in the fundamental band of CO, which is consistent with an earlier detection in a ground-based spectrum and indicates the presence of vertical mixing. Although PH_3 is expected in Y dwarfs that experience vertical mixing, it is not detected in WISE 0855. Previous ground-based M-band spectroscopy of WISE 0855 has been cited for evidence of H_2O ice clouds, but we find that the NIRSpec data in that wavelength range are matched well by our cloudless model. Thus, clear evidence of H_2O ice clouds in WISE 0855 has not been identified yet, but it may still be present in the NIRSpec data. The physical properties of WISE 0855, including the presence of H_2O clouds, can be better constrained by more detailed fitting with both cloudless and cloudy models and the incorporation of unpublished 5-28um data from the Mid-infrared Instrument on JWST. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Astronomical Journal, in press

A dependence of binary and planetary system destruction on subtle variations in the substructure in young star-forming regions

Authors: Richard J. Parker

Abstract: Simulations of the effects of stellar fly-bys on planetary systems in star-forming regions show a strong dependence on subtle variations in the initial spatial and kinematic substructure of the regions. For similar stellar densities, the more substructured star-forming regions disrupt up to a factor of two more planetary systems. We extend this work to look at the effects of substructure on stella… ▽ More Simulations of the effects of stellar fly-bys on planetary systems in star-forming regions show a strong dependence on subtle variations in the initial spatial and kinematic substructure of the regions. For similar stellar densities, the more substructured star-forming regions disrupt up to a factor of two more planetary systems. We extend this work to look at the effects of substructure on stellar binary populations. We present $N$-body simulations of substructured, and non-substructured (smooth) star-forming regions in which we place different populations of stellar binaries. We find that for binary populations that are dominated by close ($<$100au) systems, a higher proportion are destroyed in substructured regions. However, for wider systems ($>$100au), a higher proportion are destroyed in smooth regions. The difference is likely due to the hard-soft, or fast-slow boundary for binary destruction. Hard (fast/close) binaries are more likely to be destroyed in environments with a small velocity dispersion (kinematically substructured regions), whereas soft (slow/wide) binaries are more likely to be destroyed in environments with higher velocity dispersions (non-kinematically substructured regions). Due to the vast range of stellar binary semimajor axes in star-forming regions ($10^{-2} - 10^4$au) these differences are small and hence unlikely to be observable. However, planetary systems have a much smaller initial semimajor axis range (likely $\sim$1 -- 100au for gas giants) and here the difference in the fraction of companions due to substructure could be observed if the star-forming regions that disrupt planetary systems formed with similar stellar densities. △ Less

Submitted 10 August, 2023; originally announced August 2023.

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

On the origin of planetary-mass objects in NGC1333

Authors: Richard J. Parker, Catarina Alves de Oliveira

Abstract: The dominant formation mechanism of brown dwarfs and planetary mass objects in star-forming regions is presently uncertain. Do they form like stars, via the collapse and fragmentation of cores in Giant Molecular clouds, or do they form like planets in the discs around stars and are ejected via dynamical interactions? In this paper, we quantify the spatial distribution of substellar objects in NGC1… ▽ More The dominant formation mechanism of brown dwarfs and planetary mass objects in star-forming regions is presently uncertain. Do they form like stars, via the collapse and fragmentation of cores in Giant Molecular clouds, or do they form like planets in the discs around stars and are ejected via dynamical interactions? In this paper, we quantify the spatial distribution of substellar objects in NGC1333, in particular focusing on planetary-mass objects that have been the target of recent deep imaging observations. We find that these objects have a spatial distribution that is indistinguishable from the stars, and more massive brown dwarfs. We also analyse N-body simulations and find that a population of ejected planets would have a significantly different spatial and kinematic distribution to stars, and brown dwarfs that also formed through gravitational collapse and fragmentation. We therefore conclude that the low-mass substellar objects in NGC1333 formed more like stars than planets, although we predict that a population of hitherto undetected ejected planetary mass objects may be lurking in this, and other star-forming regions. △ Less

Submitted 2 August, 2023; originally announced August 2023.

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

Isotopic enrichment of planetary systems from Asymptotic Giant Branch stars

Authors: Richard J. Parker, Christina Schoettler

Abstract: Short-lived radioisotopes, in particular 26-Al and 60-Fe, are thought to contribute to the internal heating of the Earth, but are significantly more abundant in the Solar System compared to the Interstellar Medium. The presence of their decay products in the oldest Solar System objects argues for their inclusion in the Sun's protoplanetary disc almost immediately after the star formation event tha… ▽ More Short-lived radioisotopes, in particular 26-Al and 60-Fe, are thought to contribute to the internal heating of the Earth, but are significantly more abundant in the Solar System compared to the Interstellar Medium. The presence of their decay products in the oldest Solar System objects argues for their inclusion in the Sun's protoplanetary disc almost immediately after the star formation event that formed the Sun. Various scenarios have been proposed for their delivery to the Solar System, usually involving one or more core-collapse supernovae of massive stars. An alternative scenario involves the young Sun encountering an evolved Asymptotic Giant Branch (AGB) star. AGBs were previously discounted as a viable enrichment scenario for the Solar System due to the presumed low probability of an encounter between an old, evolved star and a young pre-main sequence star. We report the discovery in Gaia data of an interloping AGB star in the star-forming region NGC2264, demonstrating that old, evolved stars can encounter young forming planetary systems. We use simulations to calculate the yields of 26-Al and 60-Fe from AGBs and their contribution to the long-term geophysical heating of a planet, and find that these are comfortably within the range previously calculated for the Solar System. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 8 pages, 4 figures, accepted for publication in ApJL

Mother of Dragons: A Massive, quiescent core in the dragon cloud (IRDC G028.37+00.07)

Authors: A. T. Barnes, J. Liu, Q. Zhang, J. C. Tan, F. Bigiel, P. Caselli, G. Cosentino, F. Fontani, J. D. Henshaw, I. Jiménez-Serra, D-S. Kalb, C. Y. Law, S. N. Longmore, R. J. Parker, J. E. Pineda, A. Sánchez-Monge, W. Lim, K. Wang

Abstract: Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c… ▽ More Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature. △ Less

Submitted 31 March, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

Short-lived radioisotope enrichment in star-forming regions from stellar winds and supernovae

Authors: Richard J. Parker, Tim Lichtenberg, Miti Patel, Cheyenne K. M. Polius, Matthew Ridsdill-Smith

Abstract: The abundance of the short-lived radioisotopes 26-Al and 60-Fe in the early Solar system is usually explained by the Sun either forming from pre-enriched material, or the Sun's protosolar disc being polluted by a nearby supernova explosion from a massive star. Both hypotheses suffer from significant drawbacks: the former does not account for the dynamical evolution of star-forming regions, while i… ▽ More The abundance of the short-lived radioisotopes 26-Al and 60-Fe in the early Solar system is usually explained by the Sun either forming from pre-enriched material, or the Sun's protosolar disc being polluted by a nearby supernova explosion from a massive star. Both hypotheses suffer from significant drawbacks: the former does not account for the dynamical evolution of star-forming regions, while in the latter the time for massive stars to explode as supernovae can be similar to, or even longer than, the lifetime of protoplanetary discs. In this paper, we extend the disc enrichment scenario to include the contribution of 26-Al from the winds of massive stars before they explode as supernovae. We use N-body simulations and a post-processing analysis to calculate the amount of enrichment in each disc, and we vary the stellar density of the star-forming regions. We find that stellar winds contribute to disc enrichment to such an extent that the Solar system's 26-Al/60-Fe ratio is reproduced in up to 50 per cent of discs in dense (rho = 1000Msun pc^-3) star-forming regions. When winds are a significant contributor to the SLR enrichment, we find that Solar system levels of enrichment can occur much earlier (before 2.5 Myr) than when enrichment occurs from supernovae, which start to explode at later ages (>4 Myr). We find that Solar system levels of enrichment all but disappear in low-density star-forming regions (rho < 10Msun pc^-3), implying that the Solar system must have formed in a dense, populous star-forming region if 26-Al and 60-Fe were delivered directly to the protosolar disc from massive-star winds and supernovae. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 16 pages, 10 figures, 1 appendix with 3 further figures, accepted for publication in MNRAS

The evolution of phase space densities in star-forming regions

Authors: George A. Blaylock-Squibbs, Richard J. Parker

Abstract: The multi-dimensional phase space density (both position and velocity) of star-forming regions may encode information on the initial conditions of star and planet formation. Recently, a new metric based on the Mahalanobis distance has been used to show that hot Jupiters are more likely to be found around exoplanet host-stars in high 6D phase space density, suggesting a more dynamic formation envir… ▽ More The multi-dimensional phase space density (both position and velocity) of star-forming regions may encode information on the initial conditions of star and planet formation. Recently, a new metric based on the Mahalanobis distance has been used to show that hot Jupiters are more likely to be found around exoplanet host-stars in high 6D phase space density, suggesting a more dynamic formation environment for these planets. However, later work showed that this initial result may be due to a bias in the age of hot Jupiters and the kinematics of their host stars. We test the ability of the Mahalanobis distance and density to differentiate more generally between star-forming regions with different morphologies by applying it to static regions that are either substructured or smooth and centrally concentrated. We find that the Mahalanobis distance is unable to distinguish between different morphologies, and that the initial conditions of the N-body simulations cannot be constrained using only the Mahalanobis distance or density. Furthermore, we find that the more dimensions in the phase space the less effective the Mahalanobis density is at distinguishing between different initial conditions. We show that a combination of the mean three-dimensional (x, y, z) Mahalanobis density and the Q-parameter for a region can constrain its initial virial state. However this is due to the discriminatory power of the Q-parameter and not from any extra information imprinted in the Mahalanobis density. We therefore recommend continued use of multiple diagnostics for determining the initial conditions of star-forming regions, rather than relying on a single multi-dimensional metric. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 13 pages + appendices, 10 figures, accepted for publication in MNRAS

Evaporation before disruption: comparing timescales for Jovian planets in star-forming regions

Authors: Emma C. Daffern-Powell, Richard J. Parker

Abstract: Simulations show that the orbits of planets are readily disrupted in dense star-forming regions; planets can be exchanged between stars, become free-floating and then be captured by other stars. However, dense star-forming regions also tend to be populous, containing massive stars that emit photoionising radiation, which can evaporate the gas in protoplanetary discs. We analyse N-body simulations… ▽ More Simulations show that the orbits of planets are readily disrupted in dense star-forming regions; planets can be exchanged between stars, become free-floating and then be captured by other stars. However, dense star-forming regions also tend to be populous, containing massive stars that emit photoionising radiation, which can evaporate the gas in protoplanetary discs. We analyse N-body simulations of star-forming regions containing Jovian-mass planets and determine the times when their orbits are altered, when they become free-floating, and when they are stolen or captured. Simultaneously, we perform calculations of the evolution of protoplanetary discs when exposed to FUV radiation fields from massive stars in the same star-forming regions. In almost half (44 per cent) of the planetary systems that are disrupted - either altered, captured, stolen or become free-floating, we find that the radius of the protoplanetary disc evolves inwards, or the gas in the disc is completely evaporated, before the planets' orbits are disrupted. This implies that planets that are disrupted in dense, populous star-forming regions are more likely to be super Earths or mini Neptunes, as Jovian mass planets would not be able to form due to mass loss from photoevaporation. Furthermore, the recent discoveries of distant Jovian mass planets around tightly-packed terrestrial planets argue against their formation in populous star-forming regions, as photoevaporation would preclude gas giant planet formation at distances of more than a few au. △ Less

Submitted 27 September, 2022; originally announced September 2022.

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

Dynamics of young stellar clusters as planet forming environments

Authors: Megan Reiter, Richrd J. Parker

Abstract: Most stars and thus most planetary systems do not form in isolation. The larger star-forming environment affects protoplanetary disks in multiple ways: gravitational interactions with other stars truncate disks and alter the architectures of exoplanet systems; external irradiation from nearby high-mass stars truncates disks and shortens their lifetimes; and remaining gas and dust in the environmen… ▽ More Most stars and thus most planetary systems do not form in isolation. The larger star-forming environment affects protoplanetary disks in multiple ways: gravitational interactions with other stars truncate disks and alter the architectures of exoplanet systems; external irradiation from nearby high-mass stars truncates disks and shortens their lifetimes; and remaining gas and dust in the environment affects dynamical evolution (if removed by feedback processes) and provides some shielding for disks from external irradiation. The dynamical evolution of the region regulates when and how long various feedback mechanisms impact protoplanetary disks. Density is a key parameter that regulates the intensity and duration of UV irradiation and the frequency of dynamical encounters. The evolution of larger star-forming complexes may also play an important role by mixing populations. Observations suggest that clusters are not a single-age population but multiple populations with small age differences which may be key to resolving several timescale issues (i.e., proplyd lifetimes, enrichment). In this review, we consider stellar clusters as the ecosystems in which most stars and therefore most planets form. We review recent observational and theoretical results and highlight upcoming contributions from facilities expected to begin observations in the next five years. Looking further ahead, we argue that the next frontier is large-scale surveys of low-mass stars in more distant high-mass star-forming regions. The future of ecosystem studies is bright as faint low-mass stars in more distant high-mass star-forming regions will be routinely observable in the era of Extremely Large Telescopes (ELTs). △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: invited review accepted for publication in the EpJ Plus Focus Point on Environmental and Multiplicity Effects on Planet Formation. Guest editors: G. Lodato, C.F. Manara

Making BEASTies: dynamical formation of planetary systems around massive stars

Authors: Richard J. Parker, Emma C. Daffern-Powell

Abstract: Exoplanets display incredible diversity, from planetary system architectures around Sun-like stars that are very different to our Solar System, to planets orbiting post-main sequence stars or stellar remnants. Recently the B-star Exoplanet Abundance STudy (BEAST) reported the discovery of at least two super-Jovian planets orbiting massive stars in the Sco Cen OB association. Whilst such massive st… ▽ More Exoplanets display incredible diversity, from planetary system architectures around Sun-like stars that are very different to our Solar System, to planets orbiting post-main sequence stars or stellar remnants. Recently the B-star Exoplanet Abundance STudy (BEAST) reported the discovery of at least two super-Jovian planets orbiting massive stars in the Sco Cen OB association. Whilst such massive stars do have Keplerian discs, it is hard to envisage gas giant planets being able to form in such hostile environments. We use N-body simulations of star-forming regions to show that these systems can instead form from the capture of a free-floating planet, or the direct theft of a planet from one star to another, more massive star. We find that this occurs on average once in the first 10Myr of an association's evolution, and that the semimajor axes of the hitherto confirmed BEAST planets (290 and 556au) are more consistent with capture than theft. Our results lend further credence to the notion that planets on more distant (>100au) orbits may not be orbiting their parent star. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 5 pages, 4 figures, published in MNRAS Letters

Journal ref: MNRAS, 2022, 516, L91

The evolution of protoplanetary disc radii and disc masses in star-forming regions

Authors: Bridget Marchington, Richard J. Parker

Abstract: Protoplanetary discs are crucial to understanding how planets form and evolve, but these objects are subject to the vagaries of the birth environments of their host stars. In particular, photoionising radiation from massive stars has been shown to be an effective agent in disrupting protoplanetary discs. External photoevaporation leads to the inward evolution of the radii of discs, whereas the int… ▽ More Protoplanetary discs are crucial to understanding how planets form and evolve, but these objects are subject to the vagaries of the birth environments of their host stars. In particular, photoionising radiation from massive stars has been shown to be an effective agent in disrupting protoplanetary discs. External photoevaporation leads to the inward evolution of the radii of discs, whereas the internal viscous evolution of the disc causes the radii to evolve outwards. We couple N-body simulations of star-forming regions with a post-processing analysis of disc evolution to determine how the radius and mass distributions of protoplanetary discs evolve in young star-forming regions. To be consistent with observations, we find that the initial disc radii must be of order 100au, even though these discs are readily destroyed by photoevaporation from massive stars. Furthermore, the observed disc radii distribution in the Orion Nebula Cluster is more consistent with moderate initial stellar densities (100M$_\odot$ pc$^{-3}$), in tension with dynamical models that posit much higher inital densities for the ONC. Furthermore, we cannot reproduce the observed disc radius distribution in the Lupus star-forming region if its discs are subject to external photoevaporation. A more detailed comparison is not possible due to the well-documented uncertainties in determining the ages of pre-main sequence (disc-hosting) stars. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 17 pages + 2 pages of Appendices, 10 figures, accpeted for publication in MNRAS

Quantifying kinematic substructure in star-forming regions with statistical tests of spatial autocorrelation

Authors: Becky Arnold, Nicholas J. Wright, Richard J. Parker

Abstract: We investigate whether spatial-kinematic substructure in young star-forming regions can be quantified using Moran's $I$ statistic. Its presence in young star clusters would provide an indication that the system formed from initially substructured conditions, as expected by the hierarchical model of star cluster formation, even if the cluster were spatially smooth and centrally concentrated. Its ab… ▽ More We investigate whether spatial-kinematic substructure in young star-forming regions can be quantified using Moran's $I$ statistic. Its presence in young star clusters would provide an indication that the system formed from initially substructured conditions, as expected by the hierarchical model of star cluster formation, even if the cluster were spatially smooth and centrally concentrated. Its absence, on the other hand, would be evidence that star clusters form monolithically. The Moran's $I$ statistic is applied to $N$-body simulations of star clusters with different primordial spatial-velocity structures, and its evolution over time is studied. It is found that this statistic can be used to reliably quantify spatial-kinematic substructure, and can be used to provide evidence as to whether the spatial-kinematic structure of regions with ages $\lesssim$ 6 Myr is best reproduced by the hierarchical or monolithic models of star formation. Moran's $I$ statistic is also able to conclusively say whether the data are $not$ consistent with initial conditions that lack kinematic substructure, such as the monolithic model, in regions with ages up to, and potentially beyond, 10 Myrs. This can therefore provide a kinematic signature of the star cluster formation process that is observable for many Myr after any initial spatial structure has been erased. △ Less

Submitted 21 July, 2022; originally announced July 2022.

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

The Great Planetary Heist: Theft and capture in star-forming regions

Authors: Emma C. Daffern-Powell, Richard J. Parker, Sascha P. Quanz

Abstract: Gravitational interactions in star-forming regions are capable of disrupting and destroying planetary systems, as well as creating new ones. In particular, a planet can be stolen, where it is directly exchanged between passing stars during an interaction; or captured, where a planet is first ejected from its birth system and is free-floating for a period of time, before being captured by a passing… ▽ More Gravitational interactions in star-forming regions are capable of disrupting and destroying planetary systems, as well as creating new ones. In particular, a planet can be stolen, where it is directly exchanged between passing stars during an interaction; or captured, where a planet is first ejected from its birth system and is free-floating for a period of time, before being captured by a passing star. We perform sets of direct N-body simulations of young, substructured star-forming regions, and follow their evolution for 10 Myr in order to determine how many planets are stolen and captured, and their respective orbital properties. We show that in high density star-forming regions, stolen and captured planets have distinct properties. The semimajor axis distribution of captured planets is significantly skewed to wider orbits compared to the semimajor axis distribution of stolen planets and planets that are still orbiting their parent star (preserved planets). However, the eccentricity and inclination distributions of captured and stolen planets are similar, but in turn very different to the inclination and eccentricity distributions of preserved planets. In low-density star-forming regions these differences are not as distinct but could still, in principle, be used to determine whether observed exoplanets have likely formed in situ or have been stolen or captured. We find that the initial degree of spatial and kinematic substructure in a star-forming region is as important a factor as the stellar density in determining whether a planetary system will be altered, disrupted, captured or stolen. △ Less

Submitted 16 May, 2022; originally announced May 2022.

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

Constraints on star formation in NGC2264

Authors: Richard J. Parker, Christina Schoettler

Abstract: We quantify the spatial distribution of stars for two subclusters centred around the massive/intermediate mass stars S Mon and IRS1/2 in the NGC2264 star-forming region. We find that both subclusters are have neither a substructured, nor a centrally concentrated distribution according to the Q-parameter. Neither subcluster displays mass segregation according to the $Λ_{\rm MSR}$ ratio, but the mos… ▽ More We quantify the spatial distribution of stars for two subclusters centred around the massive/intermediate mass stars S Mon and IRS1/2 in the NGC2264 star-forming region. We find that both subclusters are have neither a substructured, nor a centrally concentrated distribution according to the Q-parameter. Neither subcluster displays mass segregation according to the $Λ_{\rm MSR}$ ratio, but the most massive stars in IRS1/2 have higher relative surface densities according to the $Σ_{\rm LDR}$ ratio. We then compare these quantities to the results of N-body simulations to constrain the initial conditions of NGC2264, which are consistent with having been dense ($\tildeρ \sim 10^4$M$_\odot$pc$^{-3}$), highly substructured and subvirial. These initial conditions were also derived from a separate analysis of the runaway and walkaway stars in the region, and indicate that star-forming regions within 1kpc of the Sun likely have a broad range of initial stellar densities. In the case of NGC2264, its initial stellar density could have been high enough to cause the destruction or truncation of protoplanetary discs and fledgling planetary systems due to dynamical encounters between stars in the early stages of its evolution. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 13 pages, 5 figures, 4 tables, 1 short appendix. Accepted for publication in MNRAS. This paper complements the results from Schoettler et al 2021 (see arXiv:2111.14892)

Investigating the structure of star-forming regions using INDICATE

Authors: George A. Blaylock-Squibbs, Richard J. Parker, Anne S. M. Buckner, Manuel Guedel

Abstract: The ability to make meaningful comparisons between theoretical and observational data of star-forming regions is key to understanding the star formation process. In this paper we test the performance of INDICATE, a new method to quantify the clustering tendencies of individual stars in a region, on synthetic star-forming regions with sub-structured, and smooth, centrally concentrated distributions… ▽ More The ability to make meaningful comparisons between theoretical and observational data of star-forming regions is key to understanding the star formation process. In this paper we test the performance of INDICATE, a new method to quantify the clustering tendencies of individual stars in a region, on synthetic star-forming regions with sub-structured, and smooth, centrally concentrated distributions. INDICATE quantifies the amount of stellar affiliation of each individual star, and also determines whether this affiliation is above random expectation for the star-forming region in question. We show that INDICATE cannot be used to quantify the overall structure of a region due to a degeneracy when applied to regions with different geometries. We test the ability of INDICATE to detect differences in the local stellar surface density and its ability to detect and quantify mass segregation. We then compare it to other methods such as the mass segregation ratio $Λ_{\rm{MSR}}$, the local stellar surface density ratio $Σ_{\rm{LDR}}$ and the cumulative distribution of stellar positions. INDICATE detects significant differences in the clustering tendencies of the most massive stars when they are at the centre of a smooth, centrally concentrated distribution, corresponding to areas of greater stellar surface density. When applied to a subset of the 50 most massive stars we show INDICATE can detect signals of mass segregation. We apply INDICATE to the following nearby star-forming regions: Taurus, ONC, NGC 1333, IC 348 and $ρ$ Ophiuchi and find a diverse range of clustering tendencies in these regions. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 16 pages + appendices, 10 figures, accepted for publication in MNRAS

Constraining the initial conditions of NGC 2264 using ejected stars found in Gaia DR2

Authors: Christina Schoettler, Richard J. Parker, Jos de Bruijne

Abstract: Fast, ejected stars have been found around several young star-forming regions, such as the Orion Nebula Cluster (ONC). These ejected stars can be used to constrain the initial density, spatial and kinematic substructure when compared to predictions from $N$-body simulations. We search for runaway and slower walkaway stars using $Gaia$ DR2 within 100 pc of NGC 2264, which contains subclustered regi… ▽ More Fast, ejected stars have been found around several young star-forming regions, such as the Orion Nebula Cluster (ONC). These ejected stars can be used to constrain the initial density, spatial and kinematic substructure when compared to predictions from $N$-body simulations. We search for runaway and slower walkaway stars using $Gaia$ DR2 within 100 pc of NGC 2264, which contains subclustered regions around higher-mass OB-stars (S Mon, IRS 1 and IRS 2). We find five runaways and nine walkaways that trace back to S Mon and six runaways and five walkaways that trace back to IRS 1/2 based on their 3D-kinematics. We compare these numbers to a range of $N$-body simulations with different initial conditions. The number of runaways/walkaways is consistent with initial conditions with a high initial stellar density ($\sim$10 000 M$_{\odot}$ pc$^{-3}$), a high initial amount of spatial substructure and either a subvirial or virialised ratio for all subclusters. We also confirm the trajectories of our ejected stars using the data from $Gaia$ EDR3, which reduces the number of runaways from IRS 1/2 from six to four but leaves the number of runaways from S Mon unchanged. The reduction in runaways is due to smaller uncertainties in the proper motion and changes in the parallax/distance estimate for these stars in $Gaia$ EDR3. We find further runaway/walkaway candidates based on proper motion alone in $Gaia$ DR2, which could increase these numbers once radial velocities are available. We also expect further changes in the candidate list with upcoming $Gaia$ data releases. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 21 pages, 6 figures, accepted for publication in MNRAS

Long-term stability of planets in and around binary stars

Authors: Harry A. Ballantyne, Tore Espaas, Bethan Z. Norgrove, Bethany A. Wootton, Benjamin R. Harris, Isaac L. Pepper, Richard D. Smith, Rosie E. Dommett, Richard J. Parker

Abstract: Planets are observed to orbit the component star(s) of stellar binary systems on so-called circumprimary or circumsecondary orbits, as well as around the entire binary system on so-called circumbinary orbits. Depending on the orbital parameters of the binary system a planet will be dynamically stable if it orbits within some critical separation of the semimajor axis in the circumprimary case, or b… ▽ More Planets are observed to orbit the component star(s) of stellar binary systems on so-called circumprimary or circumsecondary orbits, as well as around the entire binary system on so-called circumbinary orbits. Depending on the orbital parameters of the binary system a planet will be dynamically stable if it orbits within some critical separation of the semimajor axis in the circumprimary case, or beyond some critical separation for the circumbinary case. We present N-body simulations of star-forming regions that contain populations of primordial binaries to determine the fraction of binary systems that can host stable planets at various semimajor axes, and how this fraction of stable systems evolves over time. Dynamical encounters in star-forming regions can alter the orbits of some binary systems, which can induce long-term dynamical instabilities in the planetary system and can even change the size of the habitable zone(s) of the component stars. However, the overall fraction of binaries that can host stable planetary systems is not greatly affected by either the assumed binary population, or the density of the star-forming region. Instead, the critical factor in determining how many stable planetary systems exist in the Galaxy is the stellar binary fraction - the more stars that are born as singles in stellar nurseries, the higher the fraction of stable planetary systems. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 13 pages, 7 figures, 2 short appendices, accepted for publication in MNRAS. A video abstract created by Helena Gibbon is available here: https://youtu.be/76kANnTK9-s

External photoevaporation of protoplanetary discs: does location matter?

Authors: Richard J. Parker, Hayley L. Alcock, Rhana B. Nicholson, Olja Panić, Simon P. Goodwin

Abstract: Many theoretical studies have shown that external photoevaporation from massive stars can severely truncate, or destroy altogether, the gaseous protoplanetary discs around young stars. In tandem, several observational studies report a correlation between the mass of a protoplanetary disc and its distance to massive ionising stars in star-forming regions, and cite external photoevaporation by the m… ▽ More Many theoretical studies have shown that external photoevaporation from massive stars can severely truncate, or destroy altogether, the gaseous protoplanetary discs around young stars. In tandem, several observational studies report a correlation between the mass of a protoplanetary disc and its distance to massive ionising stars in star-forming regions, and cite external photoevaporation by the massive stars as the origin of this correlation. We present N-body simulations of the dynamical evolution of star-forming regions and determine the mass-loss in protoplanetary discs from external photoevaporation due to far ultraviolet (FUV) and extreme ultraviolet (EUV) radiation from massive stars. We find that projection effects can be significant, in that low-mass disc-hosting stars that appear close to the ionising sources may be fore- or background stars in the star-forming region. We find very little evidence in our simulations for a trend in increasing disc mass with increasing distance from the massive star(s), even when projection effects are ignored. Furthermore, the dynamical evolution of these young star-forming regions moves stars whose discs have been photoevaporated to far-flung locations, away from the ionising stars, and we suggest that any correlation between disc mass and distance the ionising star is either coincidental, or due to some process other than external photoevaporation. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 21 pages, 12 figures, 1 short appendix, accepted for publication in ApJ

ALMA-IRDC: Dense gas mass distribution from cloud to core scales

Authors: A. T. Barnes, J. D. Henshaw, F. Fontani, J. E. Pineda, G. Cosentino, J. C. Tan, P. Caselli, I. Jiménez-Serra, C. Y. Law, A. Avison, F. Bigiel, S. Feng, S. Kong, S. N. Longmore, L. Moser, R. J. Parker, Á. Sánchez-Monge, K. Wang

Abstract: Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm… ▽ More Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm). These observations have a high angular resolution (~3arcsec or ~0.05pc), and high continuum and spectral line sensitivity (~0.15mJy/beam and ~0.2K per 0.1km/s channel at the N2H+(1-0) transition). From the dust continuum emission, we identify 96 cores ranging from low- to high-mass (M = 3.4 to 50.9Msun) that are gravitationally bound (alpha_vir = 0.3 to 1.3) and which would require magnetic field strengths of B = 0.3 to 1.0mG to be in virial equilibrium. We combine these results with a homogenised catalogue of literature cores to recover the hierarchical structure within these clouds over four orders of magnitude in spatial scale (0.01pc to 10pc). Using supplementary observations at an even higher angular resolution, we find that the smallest fragments (<0.02pc) within this hierarchy do not currently have the mass and/or the density required to form high-mass stars. Nonetheless, the new ALMA observations presented in this paper have facilitated the identification of 19 (6 quiescent and 13 star-forming) cores that retain >16Msun without further fragmentation. These high-mass cores contain trans-sonic non-thermal motions, are kinematically sub-virial, and require moderate magnetic field strengths for support against collapse. The identification of these potential sites of high-mass star formation represents a key step in allowing us to test the predictions from high-mass star and cluster formation theories. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 21(+6 in appendix) pages, 12(+3) figures, 3(+3) tables. Machine-readable versions of Table A1, A2 and A3 are available online. Accepted for publication in MNRAS

Far and extreme ultraviolet radiation fields and consequent disc destruction in star-forming regions

Authors: Richard J. Parker, Rhana B. Nicholson, Hayley L. Alcock

Abstract: The first stages of planet formation usually occur when the host star is still in a (relatively) dense star-forming region, where the effects of the external environment may be important for understanding the outcome of the planet formation process. In particular, star-forming regions that contain massive stars have strong far ultraviolet (FUV) and extreme ultraviolet (EUV) radiation fields, which… ▽ More The first stages of planet formation usually occur when the host star is still in a (relatively) dense star-forming region, where the effects of the external environment may be important for understanding the outcome of the planet formation process. In particular, star-forming regions that contain massive stars have strong far ultraviolet (FUV) and extreme ultraviolet (EUV) radiation fields, which can induce mass-loss from protoplanetary discs due to photoevaporation. In this paper we present a parameter-space study of the expected FUV and EUV fields in N-body simulations of star-forming regions with a range of initial conditions. We then use recently published models to determine the mass-loss due to photoevaporation from protoplanetary discs. In particular, we focus on the effects of changing the initial degree of spatial structure and initial virial ratio in the star-forming regions, as well as the initial stellar density. We find that the FUV fields in star-forming regions are much higher than in the interstellar medium, even when the regions have stellar densities as low as in the Galactic field, due to the presence of intermediate-mass, and massive, stars (>5Msun). These strong radiation fields lead to the destruction of the gas component in protoplanetary discs within 1 Myr, implying that gas giant planets must either form extremely rapidly (<1 Myr), or that they exclusively form in star-forming regions like Taurus, which contain no intermediate-mass or massive stars. The latter scenario is in direct tension with meteoritic evidence from the Solar system that suggests the Sun and its protoplanetary disc was born in close proximity to massive stars. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 18 pages (including two short appendices), 6 Fgiures in main paper, accepted for publication in MNRAS

Double trouble: Gaia reveals (proto)-planetary systems that may experience more than one dense star-forming environment

Authors: Christina Schoettler, Richard J. Parker

Abstract: Planetary systems appear to form contemporaneously around young stars within young star-forming regions. Within these environments, the chances of survival, as well as the long-term evolution of these systems, are influenced by factors such as dynamical interactions with other stars and photoevaporation from massive stars. These interactions can also cause young stars to be ejected from their birt… ▽ More Planetary systems appear to form contemporaneously around young stars within young star-forming regions. Within these environments, the chances of survival, as well as the long-term evolution of these systems, are influenced by factors such as dynamical interactions with other stars and photoevaporation from massive stars. These interactions can also cause young stars to be ejected from their birth regions and become runaways. We present examples of such runaway stars in the vicinity of the Orion Nebula Cluster (ONC) found in Gaia DR2 data that have retained their discs during the ejection process. Once set on their path, these runaways usually do not encounter any other dense regions that could endanger the survival of their discs or young planetary systems. However, we show that it is possible for star-disc systems, presumably ejected from one dense star-forming region, to encounter a second dense region, in our case the ONC. While the interactions of the ejected star-disc systems in the second region are unlikely to be the same as in their birth region, a second encounter will increase the risk to the disc or planetary system from malign external effects. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 6 pages, 1 figure, 2 tables, published in MNRAS Letters

The birth environment of planetary systems

Authors: Richard J. Parker

Abstract: Star and planet formation are inextricably linked. In the earliest phases of the collapse of a protostar a disc forms around the young star and such discs are observed for the first several million years of a star's life. It is within these circumstellar, or protoplanetary, discs that the first stages of planet formation occur. Recent observations from ALMA suggest that planet formation may alread… ▽ More Star and planet formation are inextricably linked. In the earliest phases of the collapse of a protostar a disc forms around the young star and such discs are observed for the first several million years of a star's life. It is within these circumstellar, or protoplanetary, discs that the first stages of planet formation occur. Recent observations from ALMA suggest that planet formation may already be well under way after only 1 Myr of a star's life. However, stars do not form in isolation; they form from the collapse and fragmentation of giant molecular clouds several parsecs in size. This results in young stars forming in groups - often referred to as 'clusters'. In these star-forming regions the stellar density is much higher than the location of the Sun, and other stars in the Galactic disc that host exoplanets. As such, the environment where stars form has the potential to influence the planet formation process. In star-forming regions, protoplanetary discs can be truncated or destroyed by interactions with passing stars, as well as photoevaporation from the radiation fields of very massive stars. Once formed, the planets themselves can have their orbits altered by dynamical encounters - either directly from passing stars or through secondary effects such as the Kozai-Lidov mechanism. In this contribution, I review the different processes that can affect planet formation and stability in star-forming regions. I discuss each process in light of the typical range of stellar densities observed for star-forming regions. I finish by discussing these effects in the context of theories for the birth environment of the Solar System. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 18 pages (excluding references), 2 figures. Submitted for publication to a Royal Society collection of review articles. Comments, suggestions and clarifications welcome

Exoplanet detection and its dependence on stochastic sampling of the stellar Initial Mass Function

Authors: Amy L. Bottrill, Molly E. Haigh, Madeleine R. A. Hole, Sarah C. M. Theakston, Rosa B. Allen, Liam P. Grimmett, Richard J. Parker

Abstract: Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising l… ▽ More Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising locations to search for young forming exoplanets. However, due to their low numbers of stars, stochastic sampling of the IMF means their stellar populations could vary significantly. We determine the range of planet-hosting stars (spectral types A, G and M) possible from sampling the IMF multiple times, and find that some YMGs appear deficient in M-dwarfs. We then use these data to show that the expected probability of detecting terrestrial magma ocean planets is highly dependent on the exact numbers of stars produced through stochastic sampling of the IMF. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 11 pages (7 pages of main text), 2 figures (+ 3 in appendix), 3 tables. ApJ in press

Unlocking Galactic Wolf-Rayet stars with $\textit{Gaia}$ DR2 II: Cluster and Association membership

Authors: Gemma Rate, Paul A. Crowther, Richard J. Parker

Abstract: Galactic Wolf-Rayet (WR) star membership of star forming regions can be used to constrain the formation environments of massive stars. Here, we utilise $\textit{Gaia}$ DR2 parallaxes and proper motions to reconsider WR star membership of clusters and associations in the Galactic disk, supplemented by recent near-IR studies of young massive clusters. We find that only 18$-$36% of 553 WR stars exter… ▽ More Galactic Wolf-Rayet (WR) star membership of star forming regions can be used to constrain the formation environments of massive stars. Here, we utilise $\textit{Gaia}$ DR2 parallaxes and proper motions to reconsider WR star membership of clusters and associations in the Galactic disk, supplemented by recent near-IR studies of young massive clusters. We find that only 18$-$36% of 553 WR stars external to the Galactic Centre region are located in clusters, OB associations or obscured star-forming regions, such that at least 64% of the known disk WR population are isolated, in contrast with only 13% of O stars from the Galactic O star Catalogue. The fraction located in clusters, OB associations or star-forming regions rises to 25$-$41% from a global census of 663 WR stars including the Galactic Centre region. We use simulations to explore the formation processes of isolated WR stars. Neither runaways, nor low mass clusters, are numerous enough to account for the low cluster membership fraction. Rapid cluster dissolution is excluded as mass segregation ensures WR stars remain in dense, well populated environments. Only low density environments consistently produce WR stars that appeared to be isolated during the WR phase. We therefore conclude that a significant fraction of WR progenitors originate in low density association-like surroundings which expand over time. We provide distance estimates to clusters and associations host to WR stars, and estimate cluster ages from isochrone fitting. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 19 pages, 7 figures. Accepted for publication in MNRAS

Runaway and walkaway stars from the ONC with Gaia DR2

Authors: Christina Schoettler, Jos de Bruijne, Eero Vaher, Richard J. Parker

Abstract: Theory predicts that we should find fast, ejected (runaway) stars of all masses around dense, young star-forming regions. $N$-body simulations show that the number and distribution of these ejected stars could be used to constrain the initial spatial and kinematic substructure of the regions. We search for runaway and slower walkaway stars within 100 pc of the Orion Nebula Cluster (ONC) using… ▽ More Theory predicts that we should find fast, ejected (runaway) stars of all masses around dense, young star-forming regions. $N$-body simulations show that the number and distribution of these ejected stars could be used to constrain the initial spatial and kinematic substructure of the regions. We search for runaway and slower walkaway stars within 100 pc of the Orion Nebula Cluster (ONC) using $Gaia$ DR2 astrometry and photometry. We compare our findings to predictions for the number and velocity distributions of runaway stars from simulations that we run for 4 Myr with initial conditions tailored to the ONC. In $Gaia$ DR2, we find 31 runaway and 54 walkaway candidates based on proper motion, but not all of these are viable candidates in three dimensions. About 40 per cent are missing radial velocities, but we can trace back 9 3D-runaways and 24 3D-walkaways to the ONC, all of which are low/intermediate-mass (<8 M$_{\odot}$). Our simulations show that the number of runaways within 100 pc decreases the older a region is (as they quickly travel beyond this boundary), whereas the number of walkaways increases up to 3 Myr. We find fewer walkaways in $Gaia$ DR2 than the maximum suggested from our simulations, which may be due to observational incompleteness. However, the number of $Gaia$ DR2 runaways agrees with the number from our simulations during an age of $\sim$1.3-2.4 Myr, allowing us to confirm existing age estimates for the ONC (and potentially other star-forming regions) using runaway stars. △ Less

Submitted 28 April, 2020; originally announced April 2020.

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

Dynamical evolution of fractal structures in star-forming regions

Authors: Emma C. Daffern-Powell, Richard J. Parker

Abstract: The Q-parameter is used extensively to quantify the spatial distributions of stars and gas in star-forming regions as well as older clusters and associations. It quantifies the amount of structure using the ratio of the average length of a minimum spanning tree, mbar, to the average length within the complete graph, sbar. The interpretation of the Q-parameter often relies on comparing observed val… ▽ More The Q-parameter is used extensively to quantify the spatial distributions of stars and gas in star-forming regions as well as older clusters and associations. It quantifies the amount of structure using the ratio of the average length of a minimum spanning tree, mbar, to the average length within the complete graph, sbar. The interpretation of the Q-parameter often relies on comparing observed values of Q, mbar and sbar to idealised synthetic geometries, where there is little or no match between the observed star-forming regions and the synthetic regions. We measure Q, mbar, and sbar over 10 Myr in N-body simulations which are compared to IC 348, NGC 1333, and the ONC. For each star-forming region we set up simulations that approximate their initial conditions for a combination of different virial rations and fractal dimensions. We find that dynamical evolution of idealised fractal geometries can account for the observed Q, mbar, and sbar values in nearby star-forming regions. In general, an initially fractal star-forming region will tend to evolve to become more smooth and centrally concentrated. However, we show that initial conditions, as well as where the edge of the region is defined, can cause significant differences in the path that a star-forming region takes across the mbar-sbar plot as it evolves. We caution that the observed Q-parameter should not be directly compared to idealised geometries. Instead, it should be used to determine the degree to which a star-forming region is either spatially substructured or smooth and centrally concentrated. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 12 pages, 10 figures, acccepted for publication in MNRAS

On the mass segregation of cores and stars

Authors: Hayley L. Alcock, Richard J. Parker

Abstract: Observations of pre-/proto-stellar cores in young star-forming regions show them to be mass segregated, i.e. the most massive cores are centrally concentrated, whereas pre-main sequence stars in the same star-forming regions (and older regions) are not. We test whether this apparent contradiction can be explained by the massive cores fragmenting into stars of much lower mass, thereby washing out a… ▽ More Observations of pre-/proto-stellar cores in young star-forming regions show them to be mass segregated, i.e. the most massive cores are centrally concentrated, whereas pre-main sequence stars in the same star-forming regions (and older regions) are not. We test whether this apparent contradiction can be explained by the massive cores fragmenting into stars of much lower mass, thereby washing out any signature of mass segregation in pre-main sequence stars. Whilst our fragmentation model can reproduce the stellar initial mass function, we find that the resultant distribution of pre-main sequence stars is mass segregated to an even higher degree than that of the cores, because massive cores still produce massive stars if the number of fragments is reasonably low (between one and five). We therefore suggest that the reason cores are observed to be mass segregated and stars are not is likely due to dynamical evolution of the stars, which can move significant distances in star-forming regions after their formation. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 10 pages (including a short appendix), 11 figures, accepted for publication in MNRAS

Kinematic signatures of cluster formation from cool collapse in the Lagoon Nebula cluster NGC 6530

Authors: Nicholas J. Wright, Richard J. Parker

Abstract: We examine the mass-dependence of the velocity dispersion of stars in the young cluster NGC 6530 to better understand how it formed. Using a large sample of members we find that the proper motion velocity dispersion increases with stellar mass. While this trend is the opposite to that predicted if the cluster were developing energy equipartition, it is in agreement with recent N-body simulations t… ▽ More We examine the mass-dependence of the velocity dispersion of stars in the young cluster NGC 6530 to better understand how it formed. Using a large sample of members we find that the proper motion velocity dispersion increases with stellar mass. While this trend is the opposite to that predicted if the cluster were developing energy equipartition, it is in agreement with recent N-body simulations that find such a trend develops because of the Spitzer instability. In these simulations the massive stars sink to the centre of the cluster and form a self-gravitating system with a higher velocity dispersion. If the cluster has formed by the cool collapse of an initially substructured distribution then this occurs within 1-2 Myr, in agreement with our observations of NGC 6530. We therefore conclude that NGC 6530 formed from much more extended initial conditions and has since collapsed to form the cluster we see now. This cluster formation model is inconsistent with the idea that all stars form in dense, compact clusters and provides the first dynamical evidence that star clusters can form by hierarchical mergers between subclusters. △ Less

Submitted 29 August, 2019; originally announced August 2019.

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

Dynamical evolution of star-forming regions: III. Unbound stars and predictions for Gaia

Authors: Christina Schoettler, Richard J. Parker, Becky Arnold, Liam P. Grimmett, Jos de Bruijne, Nicholas J. Wright

Abstract: We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We… ▽ More We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia. △ Less

Submitted 24 May, 2019; originally announced May 2019.

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

A tale of two clusters: dynamical history determines disc survival in Tr14 and Tr16 in the Carina Nebula

Authors: Megan Reiter, Richard J. Parker

Abstract: Understanding how the birthplace of stars affects planet-forming discs is important for a comprehensive theory of planet formation. Most stars are born in dense star-forming regions where the external influence of other stars, particularly the most massive stars, will affect the survival and enrichment of their planet-forming discs. Simulations suggest that stellar dynamics play a central role in… ▽ More Understanding how the birthplace of stars affects planet-forming discs is important for a comprehensive theory of planet formation. Most stars are born in dense star-forming regions where the external influence of other stars, particularly the most massive stars, will affect the survival and enrichment of their planet-forming discs. Simulations suggest that stellar dynamics play a central role in regulating how external feedback affects discs, but comparing models to observations requires an estimate of the initial stellar density in star-forming regions. Structural analyses constrain the amount of dynamical evolution a star-forming region has experienced; regions that maintain substructure and do not show mass segregation are likely dynamically young, and therefore close to their birth density. In this paper, we present a structural analysis of two clusters in the Carina Nebula, Tr14 and Tr16. We show that neither cluster shows evidence for mass segregation or a centrally concentrated morphology, suggesting that both regions are dynamically young. This allows us to compare to simulations from Nicholson et al. (2019) who predict disc survival rates in star-forming regions of different initial densities. The surviving disc fractions in Tr14 and Tr16 are consistent with their predictions (both are $\sim 10$%), supporting a growing body of evidence that the star-forming environment plays an important role in the survival and enrichment of protoplanetary discs. △ Less

Submitted 22 April, 2019; originally announced April 2019.

Comments: MNRAS, accepted

The Gaia-ESO Survey: Asymmetric expansion of the Lagoon Nebula cluster NGC 6530 from GES and Gaia DR2

Authors: Nicholas J. Wright, R. D. Jeffries, R. J. Jackson, A. Bayo, R. Bonito, F. Damiani, V. Kalari, A. C. Lanzafame, E. Pancino, R. J. Parker, L. Prisinzano, S. Randich, J. S. Vink, E. J. Alfaro, M. Bergemann, E. Franciosini, G. Gilmore, A. Gonneau, A. Hourihane, P. Jofré, S. E. Koposov, J. Lewis, L. Magrini, G. Micela, L. Morbidelli , et al. (3 additional authors not shown)

Abstract: The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D… ▽ More The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories. △ Less

Submitted 28 March, 2019; originally announced March 2019.

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

Enlarging habitable zones around binary stars in hostile environments

Authors: Bethany A. Wootton, Richard J. Parker

Abstract: Habitable zones are regions around stars where large bodies of liquid water can be sustained on a planet or satellite. As many stars form in binary systems with non-zero eccentricity, the habitable zones around the component stars of the binary can overlap and be enlarged when the two stars are at periastron (and less often when the stars are at apastron). We perform N-body simulations of the evol… ▽ More Habitable zones are regions around stars where large bodies of liquid water can be sustained on a planet or satellite. As many stars form in binary systems with non-zero eccentricity, the habitable zones around the component stars of the binary can overlap and be enlarged when the two stars are at periastron (and less often when the stars are at apastron). We perform N-body simulations of the evolution of dense star-forming regions and show that binary systems where the component stars originally have distinct habitable zones can undergo interactions that push the stars closer together, causing the habitable zones to merge and become enlarged. Occasionally, overlapping habitable zones can occur if the component stars move further apart, but the binary becomes more eccentric. Enlargement of habitable zones happens to 1-2 binaries from an average initial total of 352 in each simulated star-forming region, and demonstrates that dense star-forming regions are not always hostile environments for planet formation and evolution. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: 6 pages, 2 figures, accepted for publication in MNRAS Letters. Press release including artist's impression at: https://ras.ac.uk/news-and-press/research-highlights/chances-life-expand-when-passing-stars-push-binaries-together

Rapid destruction of protoplanetary discs due to externalphotoevaporation in star-forming regions

Authors: Rhana B. Nicholson, Richard J. Parker, Ross P. Church, Melvyn B. Davies, Niamh M. Fearon, Sam R. J. Walton

Abstract: We analyse N-body simulations of star-forming regions to investigate the effects of external far and extreme ultra-violet photoevaporation from massive stars on protoplanetary discs. By varying the initial conditions of simulated star-forming regions, such as the spatial distribution, net bulk motion (virial ratio), and density, we investigate which parameters most affect the rate at which discs a… ▽ More We analyse N-body simulations of star-forming regions to investigate the effects of external far and extreme ultra-violet photoevaporation from massive stars on protoplanetary discs. By varying the initial conditions of simulated star-forming regions, such as the spatial distribution, net bulk motion (virial ratio), and density, we investigate which parameters most affect the rate at which discs are dispersed due to external photoevaporation. We find that disc dispersal due to external photoevaporation is faster in highly substructured star-forming regions than in smooth and centrally concentrated regions. Sub-virial star-forming regions undergoing collapse also show higher rates of disc dispersal than regions that are in virial equilibrium or are expanding. In moderately dense ($\sim$100 M$_{\odot}$ pc$^{-3}$) regions, half of all protoplanetary discs with radii $\geq$ 100 AU are photoevaporated within 1 Myr, three times faster than is currently suggested by observational studies. Discs in lower-density star-forming regions ($\sim$10 M$_{\odot}$ pc$^{-3}$) survive for longer, but half are still dispersed on short timescales ($\sim$2 Myr). This demonstrates that the initial conditions of the star forming regions will greatly impact the evolution and lifetime of protoplanetary discs. These results also imply that either gas giant planet formation is extremely rapid and occurs before the gas component of discs is evaporated, or gas giants only form in low-density star-forming regions where no massive stars are present to photoevaporate gas from protoplanetary discs. △ Less

Submitted 28 February, 2019; originally announced February 2019.

Comments: 14 pages, 5 figures

Weighing Melnick 34: the most massive binary system known

Authors: Katie A. Tehrani, Paul A. Crowther, Joachim M. Bestenlehner, Stuart P. Littlefair, A. M. T. Pollock, Richard J. Parker, Olivier Schnurr

Abstract: Here we confirm Melnick 34, an X-ray bright star in the 30 Doradus region of the Large Magellanic Cloud, as an SB2 binary comprising WN5h+WN5h components. We present orbital solutions using 26 epochs of VLT/UVES spectra and 22 epochs of archival Gemini/GMOS spectra. Radial-velocity monitoring and automated template fitting methods both reveal a similar high eccentricity system with a mass ratio cl… ▽ More Here we confirm Melnick 34, an X-ray bright star in the 30 Doradus region of the Large Magellanic Cloud, as an SB2 binary comprising WN5h+WN5h components. We present orbital solutions using 26 epochs of VLT/UVES spectra and 22 epochs of archival Gemini/GMOS spectra. Radial-velocity monitoring and automated template fitting methods both reveal a similar high eccentricity system with a mass ratio close to unity, and an orbital period in agreement with the 155.1 +/- 1 day X-ray light curve period previously derived by Pollock et al. Our favoured solution derived an eccentricity of 0.68 +/- 0.02 and mass ratio of 0.92 +/- 0.07, giving minimum masses of Ma_sin^{3}(i) = 65 +/- 7 Msun and Mb_sin^{3}(i) = 60 +/- 7 Msun. Spectral modelling using WN5h templates with CMFGEN reveals temperatures of T ~53 kK for each component and luminosities of log(La/Lsun) = 6.43 +/- 0.08 and log(Lb/Lsun) = 6.37 +/- 0.08, from which BONNSAI evolutionary modelling gives masses of Ma = 139 (+21,-18) Msun and Mb = 127 (+17,-17) Msun and ages of ~0.6 Myrs. Spectroscopic and dynamic masses would agree if Mk34 has an inclination of i ~50°, making Mk34 the most massive binary known and an excellent candidate for investigating the properties of colliding wind binaries. Within 2-3 Myrs, both components of Mk34 are expected to evolve to stellar mass black holes which, assuming the binary system survives, would make Mk34 a potential binary black hole merger progenitor and gravitational wave source. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Evolution of spatio-kinematic structures in star-forming regions: are Friends of Friends worth knowing?

Authors: Richard J. Parker, Nicholas J. Wright

Abstract: The Friends of Friends algorithm identifies groups of objects with similar spatial and kinematic properties, and has recently been used extensively to quantify the distributions of gas and stars in young star-forming regions. We apply the Friends of Friends algorithm to $N$-body simulations of the dynamical evolution of subvirial (collapsing) and supervirial (expanding) star-forming regions. We fi… ▽ More The Friends of Friends algorithm identifies groups of objects with similar spatial and kinematic properties, and has recently been used extensively to quantify the distributions of gas and stars in young star-forming regions. We apply the Friends of Friends algorithm to $N$-body simulations of the dynamical evolution of subvirial (collapsing) and supervirial (expanding) star-forming regions. We find that the algorithm picks out a wide range of groups (1 -- 25) for statistically identical initial conditions, and cannot distinguish between subvirial and supervirial regions in that we obtain similar mode and median values for the number of groups it identifies. We find no correlation between the number of groups identified initially and either the initial or subsequent spatial and kinematic tracers of the regions' evolution, such as the amount of spatial substructure, dynamical mass segregation, or velocity dispersion. We therefore urge caution in using the Friends of Friends algorithm to quantify the initial conditions of star formation. △ Less

Submitted 3 September, 2018; originally announced September 2018.

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

On the spatial distributions of dense cores in Orion B

Authors: Richard J. Parker

Abstract: We quantify the spatial distributions of dense cores in three spatially distinct areas of the Orion B star-forming region. For L1622, NGC2068/NGC2071 and NGC2023/NGC2024 we measure the amount of spatial substructure using the $\mathcal{Q}$-parameter and find all three regions to be spatially substructured ($\mathcal{Q} < 0.8$). We quantify the amount of mass segregation using $Λ_{\rm MSR}$ and fin… ▽ More We quantify the spatial distributions of dense cores in three spatially distinct areas of the Orion B star-forming region. For L1622, NGC2068/NGC2071 and NGC2023/NGC2024 we measure the amount of spatial substructure using the $\mathcal{Q}$-parameter and find all three regions to be spatially substructured ($\mathcal{Q} < 0.8$). We quantify the amount of mass segregation using $Λ_{\rm MSR}$ and find that the most massive cores are mildly mass segregated in NGC2068/NGC2071 ($Λ_{\rm MSR} \sim 2$), and very mass segregated in NGC2023/NGC2024 ($Λ_{\rm MSR} = 28^{+13}_{-10}$ for the four most massive cores). Whereas the most massive cores in L1622 are not in areas of relatively high surface density, or deeper gravitational potentials, the massive cores in NGC2068/NGC2071 and NGC2023/NGC2024 are significantly so. Given the low density (10 cores pc$^{-2}$) and spatial substructure of cores in Orion B, the mass segregation cannot be dynamical. Our results are also inconsistent with simulations in which the most massive stars form via competitive accretion, and instead hint that magnetic fields may be important in influencing the primordial spatial distributions of gas and stars in star-forming regions. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 14 pages, 10 figures (including 1 appendix), accepted for publication in MNRAS

Was Planet 9 captured in the Sun's natal star-forming region?

Authors: Richard J. Parker, Tim Lichtenberg, Sascha P. Quanz

Abstract: The presence of an unseen `Planet 9' on the outskirts of the Solar system has been invoked to explain the unexpected clustering of the orbits of several Edgeworth--Kuiper Belt Objects. We use $N$-body simulations to investigate the probability that Planet 9 was a free-floating planet (FFLOP) that was captured by the Sun in its birth star-formation environment. We find that only 1 - 6 per cent of F… ▽ More The presence of an unseen `Planet 9' on the outskirts of the Solar system has been invoked to explain the unexpected clustering of the orbits of several Edgeworth--Kuiper Belt Objects. We use $N$-body simulations to investigate the probability that Planet 9 was a free-floating planet (FFLOP) that was captured by the Sun in its birth star-formation environment. We find that only 1 - 6 per cent of FFLOPs are ensnared by stars, even with the most optimal initial conditions for capture in star-forming regions (one FFLOP per star, and highly correlated stellar velocities to facilitate capture). Depending on the initial conditions of the star-forming regions, only 5 - 10 of 10000 planets are captured onto orbits that lie within the constraints for Planet 9. When we apply an additional environmental constraint for Solar system formation - namely the injection of short-lived radioisotopes into the Sun's protoplanetary disc from supernovae - we find that the probability for the capture of Planet 9 to be almost zero. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: 5 pages, 3 figures, 1 table, accepted for publication in MNRAS Letters (31.08.2017)

Hierarchical formation of Westerlund 1: a collapsing cluster with no primordial mass segregation?

Authors: Mario Gennaro, Simon P. Goodwin, Richard J. Parker, Richard J. Allison, Wolfgang Brandner

Abstract: We examine the level of substructure and mass segregation in the massive, young cluster Westerlund 1. We find that it is relatively smooth, with little or no mass segregation, but with the massive stars in regions of significantly higher than average surface density. While an expanding or bouncing-back scenario for the evolution of Westerlund 1 cannot be ruled out, we argue that the most natural m… ▽ More We examine the level of substructure and mass segregation in the massive, young cluster Westerlund 1. We find that it is relatively smooth, with little or no mass segregation, but with the massive stars in regions of significantly higher than average surface density. While an expanding or bouncing-back scenario for the evolution of Westerlund 1 cannot be ruled out, we argue that the most natural model to explain these observations is one in which Westerlund 1 formed with no primordial mass segregation and at a similar or larger size than we now observe. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: Accepted for publication in MNRAS (11 pages, 4 figures)

How do binary clusters form?

Authors: Becky Arnold, Simon P. Goodwin, Daniel W. Griffiths, Richard J. Parker

Abstract: Approximately 10 per cent of star clusters are found in pairs, known as binary clusters. We propose a mechanism for binary cluster formation; we use N-body simulations to show that velocity substructure in a single (even fairly smooth) region can cause binary clusters to form. This process is highly stochastic and it is not obvious from a region's initial conditions whether a binary will form and,… ▽ More Approximately 10 per cent of star clusters are found in pairs, known as binary clusters. We propose a mechanism for binary cluster formation; we use N-body simulations to show that velocity substructure in a single (even fairly smooth) region can cause binary clusters to form. This process is highly stochastic and it is not obvious from a region's initial conditions whether a binary will form and, if it does, which stars will end up in which cluster. We find the probability that a region will divide is mainly determined by its virial ratio, and a virial ratio above 'equilibrium' is generally necessary for binary formation. We also find that the mass ratio of the two clusters is strongly influenced by the initial degree of spatial substructure in the region. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 10 pages, 9 figures, 1 table. Accepted for publication in MNRAS (accepted 07/07/17)

Structure and mass segregation in Galactic stellar clusters

Authors: Sami Dib, Stefan Schmeja, Richard J. Parker

Abstract: We quantify the structure of a very large number of Galactic open clusters and look for evidence of mass segregation for the most massive stars in the clusters. We characterise the structure and mass segregation ratios of 1276 clusters in the Milky Way Stellar Cluster (MWSC) catalogue containing each at least 40 stars and that are located at a distance of up to $\approx 2$ kpc from the Sun. We use… ▽ More We quantify the structure of a very large number of Galactic open clusters and look for evidence of mass segregation for the most massive stars in the clusters. We characterise the structure and mass segregation ratios of 1276 clusters in the Milky Way Stellar Cluster (MWSC) catalogue containing each at least 40 stars and that are located at a distance of up to $\approx 2$ kpc from the Sun. We use an approach based on the calculation of the minimum spanning tree of the clusters, and for each one of them, we calculate the structure parameter \Q\ and the mass segregation ratio $Λ_{\rm MSR}$. Our findings indicate that most clusters possess a \Q\ parameter that falls in the range 0.7-0.8 and are thus neither strongly concentrated nor do they show significant substructure. Only 27\% can be considered centrally concentrated with \Q\ values $> 0.8$. Of the 1276 clusters, only 14\% show indication of significant mass segregation ($Λ_{\rm MSR} > 1.5$). Furthermore, no correlation is found between the structure of the clusters or the degree of mass segregation with their position in the Galaxy. A comparison of the measured \Q\ values for the young open clusters in the MWSC to N-body numerical simulations that follow the evolution of the \Q\ parameter over the first 10 Myrs of the clusters life suggests that the young clusters found in the MWSC catalogue initially possessed local mean volume densities of $ρ_{*} \approx 10-100$ M$_{\odot}$ pc$^{-3}$. △ Less

Submitted 15 September, 2017; v1 submitted 3 July, 2017; originally announced July 2017.

Comments: Accepted to MNRAS

Journal ref: MNRAS 473, 849 (2018)

No preferential spatial distribution for massive stars expected from their formation

Authors: Richard J. Parker, James E. Dale

Abstract: We analyse N-body and Smoothed Particle Hydrodynamic (SPH) simulations of young star-forming regions to search for differences in the spatial distributions of massive stars compared to lower-mass stars. The competitive accretion theory of massive star formation posits that the most massive stars should sit in deeper potential wells than lower-mass stars. This may be observable in the relative surf… ▽ More We analyse N-body and Smoothed Particle Hydrodynamic (SPH) simulations of young star-forming regions to search for differences in the spatial distributions of massive stars compared to lower-mass stars. The competitive accretion theory of massive star formation posits that the most massive stars should sit in deeper potential wells than lower-mass stars. This may be observable in the relative surface density or spatial concentration of the most massive stars compared to other, lower-mass stars. Massive stars in cool--collapse N-body models do end up in significantly deeper potentials, and are mass segregated. However, in models of warm (expanding) star-forming regions, whilst the massive stars do come to be in deeper potentials than average stars, they are not mass segregated. In the purely hydrodynamical SPH simulations, the massive stars do come to reside in deeper potentials, which is due to their runaway growth. However, when photoionisation and stellar winds are implemented in the simulations, these feedback mechanisms regulate the mass of the stars and disrupt the inflow of gas into the clouds' potential wells. This generally makes the potential wells shallower than in the control runs, and prevents the massive stars from occupying deeper potentials. This in turn results in the most massive stars having a very similar spatial concentration and surface density distribution to lower-mass stars. Whilst massive stars do form via competitive accretion in our simulations, this rarely translates to a different spatial distribution and so any lack of primordial mass segregation in an observed star-forming region does not preclude competitive accretion as a viable formation mechanism for massive stars. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS (12.05.2017)

Dynamical histories of the IC348 and NGC1333 star-forming regions in Perseus

Authors: Richard J. Parker, Catarina Alves de Oliveira

Abstract: We present analyses of the spatial distributions of stars in the young (1 - 3 Myr) star-forming regions IC348 and NGC1333 in the Perseus Giant Molecular Cloud. We quantify the spatial structure using the $\mathcal{Q}$-parameter and find that both IC348 and NGC1333 are smooth and centrally concentrated with $\mathcal{Q}$-parameters of 0.98 and 0.89 respectively. Neither region exhibits mass segrega… ▽ More We present analyses of the spatial distributions of stars in the young (1 - 3 Myr) star-forming regions IC348 and NGC1333 in the Perseus Giant Molecular Cloud. We quantify the spatial structure using the $\mathcal{Q}$-parameter and find that both IC348 and NGC1333 are smooth and centrally concentrated with $\mathcal{Q}$-parameters of 0.98 and 0.89 respectively. Neither region exhibits mass segregation ($Λ_{\rm MSR} = 1.1^{+0.2}_{-0.3}$ for IC348 and $Λ_{\rm MSR} = 1.2^{+0.4}_{-0.3}$ for NGC1333, where $Λ_{\rm MSR} \sim 1$ corresponds to no mass segregation), nor do the most massive stars reside in areas of enhanced stellar surface density compared to the average surface density, according to the $Σ_{\rm LDR}$ method. We then constrain the dynamical histories and hence initial conditions of both regions by comparing the observed values to $N$-body simulations at appropriate ages. Stars in both regions likely formed with sub-virial velocities which contributed to merging of substructure and the formation of smooth clusters. The initial stellar densities were no higher than $ρ\sim 100 - 500$M$_\odot$pc$^{-3}$ for IC348 and $ρ\sim 500 - 2000$M$_\odot$pc$^{-3}$ for NGC1333. These initial densities, in particular that of NGC1333, are high enough to facilitate dynamical interactions which would likely affect $\sim$10 per cent of protoplanetary discs and binary stars. △ Less

Submitted 24 March, 2017; originally announced March 2017.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS (23.03.2017)

The Gaia-ESO Survey: Structural and dynamical properties of the young cluster Chamaeleon I

Authors: G. G. Sacco, L. Spina, S. Randich, F. Palla, R. J. Parker, R. D. Jeffries, R. Jackson, M. R. Meyer, M. Mapelli, A. C. Lanzafame, R. Bonito, F. Damiani, E. Franciosini, A. Frasca, A. Klutsch, L. Prisinzano, E. Tognelli, S. Degl'Innocenti, P. G. Prada Moroni, E. J. Alfaro, G. Micela, T. Prusti, D. Barrado, K. Biazzo, H. Bouy , et al. (22 additional authors not shown)

Abstract: The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersi… ▽ More The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $γ$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density Σand the level of mass segregation $Λ_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: Accepted for publication on Astronomy and Astrophysics

Journal ref: A&A 601, A97 (2017)

Supernova enrichment of planetary systems in low-mass star clusters

Authors: Rhana B. Nicholson, Richard J. Parker

Abstract: The presence and abundance of short lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe in chondritic meteorites implies that the Sun formed in the vicinity of one or more massive stars that exploded as supernovae (SNe). Massive stars are more likely to form in massive star clusters ($>$1000 M$_{\odot}$) than lower mass clusters. However, photoevaporation of protoplanetary discs from massive stars… ▽ More The presence and abundance of short lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe in chondritic meteorites implies that the Sun formed in the vicinity of one or more massive stars that exploded as supernovae (SNe). Massive stars are more likely to form in massive star clusters ($>$1000 M$_{\odot}$) than lower mass clusters. However, photoevaporation of protoplanetary discs from massive stars and dynamical interactions with passing stars can inhibit planet formation in clusters with radii of $\sim$1 pc. We investigate whether low-mass (50 - 200 M$_{\odot}$) star clusters containing one or two massive stars are a more likely avenue for early Solar system enrichment as they are more dynamically quiescent. We analyse $N$-body simulations of the evolution of these low-mass clusters and find that a similar fraction of stars experience supernova enrichment than in high mass clusters, despite their lower densities. This is due to two-body relaxation, which causes a significant expansion before the first supernova even in clusters with relatively low (100 stars pc$^{-3}$) initial densities. However, because of the high number of low mass clusters containing one or two massive stars, the absolute number of enriched stars is the same, if not higher than for more populous clusters. Our results show that direct enrichment of protoplanetary discs from supernovae occurs as frequently in low mass clusters containing one or two massive stars (>20 M$_{\odot}$) as in more populous star clusters (1000 M$_\odot$). This relaxes the constraints on the direct enrichment scenario and therefore the birth environment of the Solar System. △ Less

Submitted 17 October, 2016; originally announced October 2016.

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

The effects of supernovae on the dynamical evolution of binary stars and star clusters

Authors: Richard J. Parker

Abstract: In this chapter I review the effects of supernovae explosions on the dynamical evolution of (1) binary stars and (2) star clusters. (1) Supernovae in binaries can drastically alter the orbit of the system, sometimes disrupting it entirely, and are thought to be partially responsible for `runaway' massive stars - stars in the Galaxy with large peculiar velocities. The ejection of the lower-mass s… ▽ More In this chapter I review the effects of supernovae explosions on the dynamical evolution of (1) binary stars and (2) star clusters. (1) Supernovae in binaries can drastically alter the orbit of the system, sometimes disrupting it entirely, and are thought to be partially responsible for `runaway' massive stars - stars in the Galaxy with large peculiar velocities. The ejection of the lower-mass secondary component of a binary occurs often in the event of the more massive primary star exploding as a supernova. The orbital properties of binaries that contain massive stars mean that the observed velocities of runaway stars (10s - 100s km s$^{-1}$) are consistent with this scenario. (2) Star formation is an inherently inefficient process, and much of the potential in young star clusters remains in the form of gas. Supernovae can in principle expel this gas, which would drastically alter the dynamics of the cluster by unbinding the stars from the potential. However, recent numerical simulations, and observational evidence that gas-free clusters are observed to be bound, suggest that the effects of supernova explosions on the dynamics of star clusters are likely to be minimal. △ Less

Submitted 29 September, 2016; v1 submitted 19 September, 2016; originally announced September 2016.

Comments: 16 pages, to appear in the 'Handbook of Supernovae', eds. Paul Murdin and Athem Alsabti. This version replaces an earlier version that contained several typos

Isotopic enrichment of forming planetary systems from supernova pollution

Authors: Tim Lichtenberg, Richard J. Parker, Michael R. Meyer

Abstract: Heating by short-lived radioisotopes (SLRs) such as aluminum-26 and iron-60 fundamentally shaped the thermal history and interior structure of Solar System planetesimals during the early stages of planetary formation. The subsequent thermo-mechanical evolution, such as internal differentiation or rapid volatile degassing, yields important implications for the final structure, composition and evolu… ▽ More Heating by short-lived radioisotopes (SLRs) such as aluminum-26 and iron-60 fundamentally shaped the thermal history and interior structure of Solar System planetesimals during the early stages of planetary formation. The subsequent thermo-mechanical evolution, such as internal differentiation or rapid volatile degassing, yields important implications for the final structure, composition and evolution of terrestrial planets. SLR-driven heating in the Solar System is sensitive to the absolute abundance and homogeneity of SLRs within the protoplanetary disk present during the condensation of the first solids. In order to explain the diverse compositions found for extrasolar planets, it is important to understand the distribution of SLRs in active planet formation regions (star clusters) during their first few Myr of evolution. By constraining the range of possible effects, we show how the imprint of SLRs can be extrapolated to exoplanetary systems and derive statistical predictions for the distribution of aluminum-26 and iron-60 based on N-body simulations of typical to large clusters (1000-10000 stars) with a range of initial conditions. We quantify the pollution of protoplanetary disks by supernova ejecta and show that the likelihood of enrichment levels similar to or higher than the Solar System can vary considerably, depending on the cluster morphology. Furthermore, many enriched systems show an excess in radiogenic heating compared to Solar System levels, which implies that the formation and evolution of planetesimals could vary significantly depending on the birth environment of their host stars. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 15 pages, 8 figures, 4 tables; accepted for publication in MNRAS; associated video files can be found at http://timlichtenberg.net/2016_enrichment.html

Investigating the structure and fragmentation of a highly filamentary IRDC

Authors: J. D. Henshaw, P. Caselli, F. Fontani, I. Jimenez-Serra, J. C. Tan, S. N. Longmore, J. E. Pineda, R. J. Parker, A. T. Barnes

Abstract: We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with th… ▽ More We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with that predicted by theoretical work describing the fragmentation of hydrodynamic fluid cylinders, finding a significant (a factor of ~8) discrepancy between the two. Our observations are consistent with the picture emerging from kinematic studies of molecular clouds suggesting that the cores are harboured within a complex network of independent sub-filaments. This result emphasizes the importance of considering the underlying physical structure, and potentially, dynamically important magnetic fields, in any fragmentation analysis. The identified cores exhibit a range in (peak) beam-averaged column density ($3.6{\rm x}10^{23}{\rm cm}^{-2}<N_{H,c}<8.0{\rm x}10^{23}{\rm cm}^{-2}$), mass ($8.1M_{\odot}<M_{c}<26.1M_{\odot}$), and number density ($6.1{\rm x}10^{5}{\rm cm}^{-3}<n_{H, c, eq}<14.7{\rm x}10^{5}{\rm cm}^{-3}$). Two of these cores, dark in the mid-infrared, centrally-concentrated, monolithic (with no traceable substructure at our PdBI resolution), and with estimated masses of the order ~20-25$M_{\odot}$, are good candidates for the progenitors of intermediate-to-high-mass stars. Virial parameters span a range $0.2<α_{\rm vir}<1.3$. Without additional support, possibly from dynamically important magnetic fields with strengths of the order 230$μ$G<B<670$μ$G, the cores are susceptible to gravitational collapse. These results may imply a multilayered fragmentation process, which incorporates the formation of sub-filaments, embedded cores, and the possibility of further fragmentation. △ Less

Submitted 15 August, 2016; v1 submitted 25 July, 2016; originally announced July 2016.

Comments: 14 pages (plus 10 page appendix), 6 figures (plus 11 in appendix), 3 tables. Accepted for publication in MNRAS. Abstract modified due to character limit. This updated version includes minor corrections (including typos and updated references)