-
A Systematic Review of Edge Case Detection in Automated Driving: Methods, Challenges and Future Directions
Authors:
Saeed Rahmani,
Sabine Rieder,
Erwin de Gelder,
Marcel Sonntag,
Jorge Lorente Mallada,
Sytze Kalisvaart,
Vahid Hashemi,
Simeon C. Calvert
Abstract:
The rapid development of automated vehicles (AVs) promises to revolutionize transportation by enhancing safety and efficiency. However, ensuring their reliability in diverse real-world conditions remains a significant challenge, particularly due to rare and unexpected situations known as edge cases. Although numerous approaches exist for detecting edge cases, there is a notable lack of a comprehen…
▽ More
The rapid development of automated vehicles (AVs) promises to revolutionize transportation by enhancing safety and efficiency. However, ensuring their reliability in diverse real-world conditions remains a significant challenge, particularly due to rare and unexpected situations known as edge cases. Although numerous approaches exist for detecting edge cases, there is a notable lack of a comprehensive survey that systematically reviews these techniques. This paper fills this gap by presenting a practical, hierarchical review and systematic classification of edge case detection and assessment methodologies. Our classification is structured on two levels: first, categorizing detection approaches according to AV modules, including perception-related and trajectory-related edge cases; and second, based on underlying methodologies and theories guiding these techniques. We extend this taxonomy by introducing a new class called "knowledge-driven" approaches, which is largely overlooked in the literature. Additionally, we review the techniques and metrics for the evaluation of edge case detection methods and identified edge cases. To our knowledge, this is the first survey to comprehensively cover edge case detection methods across all AV subsystems, discuss knowledge-driven edge cases, and explore evaluation techniques for detection methods. This structured and multi-faceted analysis aims to facilitate targeted research and modular testing of AVs. Moreover, by identifying the strengths and weaknesses of various approaches and discussing the challenges and future directions, this survey intends to assist AV developers, researchers, and policymakers in enhancing the safety and reliability of automated driving (AD) systems through effective edge case detection.
△ Less
Submitted 10 October, 2024;
originally announced October 2024.
-
Massive Star Cluster Formation with Binaries. I. Evolution of Binary Populations
Authors:
Claude Cournoyer-Cloutier,
Alison Sills,
William E. Harris,
Brooke Polak,
Steven Rieder,
Eric P. Andersson,
Sabrina M. Appel,
Mordecai-Mark Mac Low,
Stephen McMillan,
Simon Portegies Zwart
Abstract:
We study the evolution of populations of binary stars within massive cluster-forming regions. We simulate the formation of young massive star clusters within giant molecular clouds with masses ranging from 2 x 10$^{4}$ to 3.2 x 10$^{5}$ M$_{\odot}$. We use Torch, which couples stellar dynamics, magnetohydrodynamics, star and binary formation, stellar evolution, and stellar feedback through the AMU…
▽ More
We study the evolution of populations of binary stars within massive cluster-forming regions. We simulate the formation of young massive star clusters within giant molecular clouds with masses ranging from 2 x 10$^{4}$ to 3.2 x 10$^{5}$ M$_{\odot}$. We use Torch, which couples stellar dynamics, magnetohydrodynamics, star and binary formation, stellar evolution, and stellar feedback through the AMUSE framework. We find that the binary fraction decreases during cluster formation at all molecular cloud masses. The binaries' orbital properties also change, with stronger and quicker changes in denser, more massive clouds. Most of the changes we see can be attributed to the disruption of binaries wider than 100 au, although the close binary fraction also decreases in the densest cluster-forming region. The binary fraction for O stars remains above 90%, but exchanges and dynamical hardening are ubiquitous, indicating that O stars undergo frequent few-body interactions early during the cluster formation process. Changes to the populations of binaries are a by-product of hierarchical cluster assembly: most changes to the binary population take place when the star formation rate is high and there are frequent mergers between sub-clusters in the cluster-forming region. A universal primordial binary distribution based on observed inner companions in the Galactic field is consistent with the binary populations of young clusters with resolved stellar populations, and the scatter between clusters of similar masses could be explained by differences in their formation history.
△ Less
Submitted 9 October, 2024;
originally announced October 2024.
-
Two waves of massive stars running away from the young cluster R136
Authors:
Mitchel Stoop,
Alex de Koter,
Lex Kaper,
Sarah Brands,
Simon Portegies Zwart,
Hugues Sana,
Fiorenzo Stoppa,
Mark Gieles,
Laurent Mahy,
Tomer Shenar,
Difeng Guo,
Gijs Nelemans,
Steven Rieder
Abstract:
Massive stars are predominantly born in stellar associations or clusters. Their radiation fields, stellar winds, and supernovae strongly impact their local environment. In the first few million years of a cluster's life, massive stars are dynamically ejected running away from the cluster at high speed. However, the production rate of dynamically ejected runaways is poorly constrained. Here we repo…
▽ More
Massive stars are predominantly born in stellar associations or clusters. Their radiation fields, stellar winds, and supernovae strongly impact their local environment. In the first few million years of a cluster's life, massive stars are dynamically ejected running away from the cluster at high speed. However, the production rate of dynamically ejected runaways is poorly constrained. Here we report on a sample of 55 massive runaway stars ejected from the young cluster R136 in the Large Magellanic Cloud. Astrometric analysis with Gaia reveals two channels of dynamically ejected runaways. The first channel ejects massive stars in all directions and is consistent with dynamical interactions during and after the birth of R136. The second channel launches stars in a preferred direction and may be related to a cluster interaction. We find that 23-33% of the most luminous stars initially born in R136 are runaways. Model predictions have significantly underestimated the dynamical escape fraction of massive stars. Consequently, their role in shaping and heating the interstellar and galactic medium, along with their role in driving galactic outflows, is far more important than previously thought.
△ Less
Submitted 8 October, 2024;
originally announced October 2024.
-
Gaussian-Based and Outside-the-Box Runtime Monitoring Join Forces
Authors:
Vahid Hashemi,
Jan Křetínský,
Sabine Rieder,
Torsten Schön,
Jan Vorhoff
Abstract:
Since neural networks can make wrong predictions even with high confidence, monitoring their behavior at runtime is important, especially in safety-critical domains like autonomous driving. In this paper, we combine ideas from previous monitoring approaches based on observing the activation values of hidden neurons. In particular, we combine the Gaussian-based approach, which observes whether the…
▽ More
Since neural networks can make wrong predictions even with high confidence, monitoring their behavior at runtime is important, especially in safety-critical domains like autonomous driving. In this paper, we combine ideas from previous monitoring approaches based on observing the activation values of hidden neurons. In particular, we combine the Gaussian-based approach, which observes whether the current value of each monitored neuron is similar to typical values observed during training, and the Outside-the-Box monitor, which creates clusters of the acceptable activation values, and, thus, considers the correlations of the neurons' values. Our experiments evaluate the achieved improvement.
△ Less
Submitted 8 October, 2024;
originally announced October 2024.
-
Monitizer: Automating Design and Evaluation of Neural Network Monitors
Authors:
Muqsit Azeem,
Marta Grobelna,
Sudeep Kanav,
Jan Kretinsky,
Stefanie Mohr,
Sabine Rieder
Abstract:
The behavior of neural networks (NNs) on previously unseen types of data (out-of-distribution or OOD) is typically unpredictable. This can be dangerous if the network's output is used for decision-making in a safety-critical system. Hence, detecting that an input is OOD is crucial for the safe application of the NN. Verification approaches do not scale to practical NNs, making runtime monitoring m…
▽ More
The behavior of neural networks (NNs) on previously unseen types of data (out-of-distribution or OOD) is typically unpredictable. This can be dangerous if the network's output is used for decision-making in a safety-critical system. Hence, detecting that an input is OOD is crucial for the safe application of the NN. Verification approaches do not scale to practical NNs, making runtime monitoring more appealing for practical use. While various monitors have been suggested recently, their optimization for a given problem, as well as comparison with each other and reproduction of results, remain challenging. We present a tool for users and developers of NN monitors. It allows for (i) application of various types of monitors from the literature to a given input NN, (ii) optimization of the monitor's hyperparameters, and (iii) experimental evaluation and comparison to other approaches. Besides, it facilitates the development of new monitoring approaches. We demonstrate the tool's usability on several use cases of different types of users as well as on a case study comparing different approaches from recent literature.
△ Less
Submitted 16 May, 2024;
originally announced May 2024.
-
Observational Bias and Young Massive Cluster Characterisation II. Can Gaia accurately observe young clusters and associations?
Authors:
Anne S. M. Buckner,
Tim Naylor,
Clare L. Dobbs,
Steven Rieder,
Thomas J. R. Bending
Abstract:
Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars and extinction. As such, how accurate are measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion- type), synthetically observe them to obtain luminosities, accounting for extinction and the inherent limits of Gaia, then place them within…
▽ More
Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars and extinction. As such, how accurate are measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion- type), synthetically observe them to obtain luminosities, accounting for extinction and the inherent limits of Gaia, then place them within the real Gaia DR3 catalogue. We then attempt to rediscover the clusters at distances of between 500pc and 4300pc. We show the spatial and kinematic criteria which are best able to pick out the simulated clusters, maximising completeness and minimising contamination. We then compare the properties of the 'observed' clusters with the original simulations. We looked at the degree of clustering, the identification of clusters and subclusters within the datasets, and whether the clusters are expanding or contracting. Even with a high level of incompleteness (e.g. $<2\%$ stellar members identified), similar qualitative conclusions tend to be reached compared to the original dataset, but most quantitative conclusions are likely to be inaccurate. Accurate determination of the number, stellar membership and kinematic properties of subclusters, are the most problematic to correctly determine, particularly at larger distances due to the disappearance of cluster substructure as the data become more incomplete, but also at smaller distances where the misidentification of asterisms as true structure can be problematic. Unsurprisingly, we tend to obtain better quantitative agreement of properties for our more massive Westerlund2-type cluster. We also make optical style images of the clusters over our range of distances.
△ Less
Submitted 13 November, 2023; v1 submitted 31 October, 2023;
originally announced October 2023.
-
Guessing Winning Policies in LTL Synthesis by Semantic Learning
Authors:
Jan Kretinsky,
Tobias Meggendorfer,
Maximilian Prokop,
Sabine Rieder
Abstract:
We provide a learning-based technique for guessing a winning strategy in a parity game originating from an LTL synthesis problem. A cheaply obtained guess can be useful in several applications. Not only can the guessed strategy be applied as best-effort in cases where the game's huge size prohibits rigorous approaches, but it can also increase the scalability of rigorous LTL synthesis in several w…
▽ More
We provide a learning-based technique for guessing a winning strategy in a parity game originating from an LTL synthesis problem. A cheaply obtained guess can be useful in several applications. Not only can the guessed strategy be applied as best-effort in cases where the game's huge size prohibits rigorous approaches, but it can also increase the scalability of rigorous LTL synthesis in several ways. Firstly, checking whether a guessed strategy is winning is easier than constructing one. Secondly, even if the guess is wrong in some places, it can be fixed by strategy iteration faster than constructing one from scratch. Thirdly, the guess can be used in on-the-fly approaches to prioritize exploration in the most fruitful directions.
In contrast to previous works, we (i)~reflect the highly structured logical information in game's states, the so-called semantic labelling, coming from the recent LTL-to-automata translations, and (ii)~learn to reflect it properly by learning from previously solved games, bringing the solving process closer to human-like reasoning.
△ Less
Submitted 24 May, 2023;
originally announced May 2023.
-
Bringing Stellar Evolution & Feedback Together: Summary of proposals from the Lorentz Center Workshop, 2022
Authors:
Sam Geen,
Poojan Agrawal,
Paul A. Crowther,
B. W. Keller,
Alex de Koter,
Zsolt Keszthelyi,
Freeke van de Voort,
Ahmad A. Ali,
Frank Backs,
Lars Bonne,
Vittoria Brugaletta,
Annelotte Derkink,
Sylvia Ekström,
Yvonne A. Fichtner,
Luca Grassitelli,
Ylva Götberg,
Erin R. Higgins,
Eva Laplace,
Kong You Liow,
Marta Lorenzo,
Anna F. McLeod,
Georges Meynet,
Megan Newsome,
G. André Oliva,
Varsha Ramachandran
, et al. (12 additional authors not shown)
Abstract:
Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this whit…
▽ More
Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as ``feedback''. Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates, and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting `Bringing Stellar Evolution and Feedback Together' in April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in.
△ Less
Submitted 31 January, 2023;
originally announced January 2023.
-
Using Molecular Gas Observations to Guide Initial Conditions for Star Cluster Simulations
Authors:
Alison Sills,
Steven Rieder,
Anne S. M. Buckner,
Alvaro Hacar,
Simon Portegies Zwart,
Paula S. Teixeira
Abstract:
The earliest evolution of star clusters involves a phase of co-existence of both newly-formed stars, and the gas from which they are forming. Observations of the gas in such regions provide a wealth of data that can inform the simulations which are needed to follow the evolution of such objects forward in time. We present a method for transforming the observed gas properties into initial condition…
▽ More
The earliest evolution of star clusters involves a phase of co-existence of both newly-formed stars, and the gas from which they are forming. Observations of the gas in such regions provide a wealth of data that can inform the simulations which are needed to follow the evolution of such objects forward in time. We present a method for transforming the observed gas properties into initial conditions for simulations that include gas, stars, and ongoing star formation. We demonstrate our technique using the Orion Nebula Cluster. Since the observations cannot provide all the necessary information for our simulations, we make choices for the missing data and assess the impact of those choices. We find that the results are insensitive to the adopted choices of the gas velocity in the plane of the sky. The properties of the surrounding gas cloud (e.g. overall density and size), however, have an effect on the star formation rate and pace of assembly of the resultant star cluster. We also analyze the stellar properties of the cluster and find that the stars become more tightly clustered and in a stronger radial distribution even as new stars form in the filament.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Runtime Monitoring for Out-of-Distribution Detection in Object Detection Neural Networks
Authors:
Vahid Hashemi,
Jan Křetínsky,
Sabine Rieder,
Jessica Schmidt
Abstract:
Runtime monitoring provides a more realistic and applicable alternative to verification in the setting of real neural networks used in industry. It is particularly useful for detecting out-of-distribution (OOD) inputs, for which the network was not trained and can yield erroneous results. We extend a runtime-monitoring approach previously proposed for classification networks to perception systems…
▽ More
Runtime monitoring provides a more realistic and applicable alternative to verification in the setting of real neural networks used in industry. It is particularly useful for detecting out-of-distribution (OOD) inputs, for which the network was not trained and can yield erroneous results. We extend a runtime-monitoring approach previously proposed for classification networks to perception systems capable of identification and localization of multiple objects. Furthermore, we analyze its adequacy experimentally on different kinds of OOD settings, documenting the overall efficacy of our approach.
△ Less
Submitted 15 December, 2022;
originally announced December 2022.
-
The early evolution of young massive clusters: The kinematic history of NGC6611 / M16
Authors:
Mitchel Stoop,
Lex Kaper,
Alex de Koter,
Difeng Guo,
Henny J. G. L. M. Lamers,
Steven Rieder
Abstract:
In the first few Myr the massive stars dynamically interact, produce runaways and affect the initial binary population. Observing and interpreting the dynamics of young massive clusters is key to our understanding of the star formation process and predicting the outcome of stellar evolution. We have studied NGC6611 in the Eagle Nebula (M16), a young massive cluster hosting 19 O stars. We used Gaia…
▽ More
In the first few Myr the massive stars dynamically interact, produce runaways and affect the initial binary population. Observing and interpreting the dynamics of young massive clusters is key to our understanding of the star formation process and predicting the outcome of stellar evolution. We have studied NGC6611 in the Eagle Nebula (M16), a young massive cluster hosting 19 O stars. We used Gaia EDR3 data to determine the membership, age, cluster dynamics and the kinematics of the massive stars including runaways. The membership analysis yields 137 members located at a mean distance of 1706 $\pm$ 7 pc. The colour - absolute magnitude diagram reveals a blue and a red population of pre-main-sequence stars, consistent with two distinct populations of stars. In line with earlier studies, the youngest population has a mean extinction $A_V$ = 3.6 $\pm$ 0.1 mag and an age = 1.3 $\pm$ 0.2 Myr, while the older population of stars has a mean extinction $A_V$ = 2.0 $\pm$ 0.1 mag and an age = 7.5 $\pm$ 0.4 Myr. The latter population is more spatially extended than the younger generation of stars. We argue that most of the OB stars belong to the younger population. We identify 8 runaways originating from the center of NGC6611, consistent with the dynamical ejection scenario. We show that ~ 50% of the O stars have velocities comparable to or greater than the escape velocity. These O stars can be traced back to the center of NGC6611 with kinematic ages ranging from 0 to 2 Myr. This suggests that dynamical interactions played an important role in the early evolution of NGC6611, which is surprising considering the low current stellar density. Comparing this to simulations of young massive clusters, the required initial radius of 0.1-0.5 pc is not consistent with that of NGC6611. The O stars could have initially formed in wide binaries and possibly harden through dynamical interactions.
△ Less
Submitted 6 December, 2022; v1 submitted 18 July, 2022;
originally announced July 2022.
-
Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes
Authors:
Alessandro Alberto Trani,
Steven Rieder,
Ataru Tanikawa,
Giuliano Iorio,
Riccardo Martini,
Georgii Karelin,
Hila Glanz,
Simon Portegies Zwart
Abstract:
We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory.…
▽ More
We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory. Our formalism can model the eccentricity during the common envelope inspiral, and it gives results consistent with smoothed particles hydrodynamical simulations. We apply our formalism to common envelope events from binary population synthesis models and find that the final eccentricity distribution resembles the observed distribution of post-common-envelope binaries. Our model can be used for time-resolved common-envelope evolution in population synthesis calculations or as part of binary interactions in direct N-body simulations of star clusters.
△ Less
Submitted 23 August, 2022; v1 submitted 26 May, 2022;
originally announced May 2022.
-
Observational Bias and Young Massive Cluster Characterisation I. 2D Perspective Effects
Authors:
Anne S. M. Buckner,
Kong You Liow,
Clare L. Dobbs,
Tim Naylor,
Steven Rieder
Abstract:
Understanding the formation and evolution of high mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are re…
▽ More
Understanding the formation and evolution of high mass star clusters requires comparisons between theoretical and observational data to be made. Unfortunately, while the full phase space of simulated regions is available, often only partial 2D spatial and kinematic data is available for observed regions. This raises the question as to whether cluster parameters determined from 2D data alone are reliable and representative of clusters real parameters and the impact of line-of-sight orientation. In this paper we derive parameters for a simulated cluster formed from a cloud-cloud collision with the full 6D phase space, and compare them with those derived from three different 2D line-of-sight orientations for the cluster. We show the same qualitative conclusions can be reached when viewing clusters in 2D versus 3D, but that drawing quantitative conclusions when viewing in 2D is likely to be inaccurate. The greatest divergence occurs in the perceived kinematics of the cluster, which in some orientations appears to be expanding when the cluster is actually contracting. Increases in the cluster density compounds pre-existing perspective issues, reducing the relative accuracy and consistency of properties derived from different orientations. This is particularly problematic for determination of the number, and membership, of subclusters present in the cluster. We find the fraction of subclusters correctly identified in 2D decreases as the cluster evolves, reaching less than 3.4% at the evolutionary end point for our cluster.
△ Less
Submitted 6 May, 2022;
originally announced May 2022.
-
RE-EDS Using GAFF Topologies: Application to Relative Hydration Free-Energy Calculations for Large Sets of Molecules
Authors:
Salome R. Rieder,
Benjamin Ries,
Kay Schaller,
Candide Champion,
Emilia P. Barros,
Philippe H. Huenenberger,
Sereina Riniker
Abstract:
Free-energy differences between pairs of end-states can be estimated based on molecular dynamics (MD) simulations using standard pathway-dependent methods such as thermodynamic integration (TI), free-energy perturbation, or Bennett's acceptance ratio. Replica-exchange enveloping distribution sampling (RE-EDS), on the other hand, allows for the sampling of multiple end-states in a single simulation…
▽ More
Free-energy differences between pairs of end-states can be estimated based on molecular dynamics (MD) simulations using standard pathway-dependent methods such as thermodynamic integration (TI), free-energy perturbation, or Bennett's acceptance ratio. Replica-exchange enveloping distribution sampling (RE-EDS), on the other hand, allows for the sampling of multiple end-states in a single simulation without the specification of any pathways. In this work, we use the RE-EDS method as implemented in GROMOS together with generalized AMBER force field (GAFF) topologies, converted to a GROMOS-compatible format with a newly developed GROMOS++ program amber2gromos, to compute relative hydration free energies for a series of benzene derivatives. The results obtained with RE-EDS are compared to the experimental data as well as calculated values from the literature. In addition, the estimated free-energy differences in water and in vacuum are compared to values from TI calculations carried out with GROMACS. The hydration free energies obtained using RE-EDS for multiple molecules are found to be in good agreement with both the experimental data and the results calculated using other free-energy methods. While all considered free-energy methods delivered accurate results, the RE-EDS calculations required the least amount of total simulation time. This work serves as a validation for the use of GAFF topologies with the GROMOS simulation package and the RE-EDS approach. Furthermore, the performance of RE-EDS for a large set of 28 end-states is assessed with promising results.
△ Less
Submitted 10 June, 2022; v1 submitted 4 April, 2022;
originally announced April 2022.
-
Grouped star formation: converting sink particles to stars in hydrodynamical simulations
Authors:
K. Y. Liow,
S. Rieder,
C. L. Dobbs,
S. E. Jaffa
Abstract:
Modelling star formation and resolving individual stars in numerical simulations of molecular clouds and galaxies is highly challenging. Simulations on very small scales can be sufficiently well resolved to consistently follow the formation of individual stars, whilst on larger scales sinks that have masses sufficient to fully sample the IMF can be converted into realistic stellar populations. How…
▽ More
Modelling star formation and resolving individual stars in numerical simulations of molecular clouds and galaxies is highly challenging. Simulations on very small scales can be sufficiently well resolved to consistently follow the formation of individual stars, whilst on larger scales sinks that have masses sufficient to fully sample the IMF can be converted into realistic stellar populations. However, as yet, these methods do not work for intermediate scale resolutions whereby sinks are more massive compared to individual stars but do not fully sample the IMF. In this paper, we introduce the grouped star formation prescription, whereby sinks are first grouped according to their positions, velocities, and ages, then stars are formed by sampling the IMF using the mass of the groups. We test our grouped star formation method in simulations of various physical scales, from sub-parsec to kilo-parsec, and from static isolated clouds to colliding clouds. With suitable grouping parameters, this star formation prescription can form stars that follow the IMF and approximately mimic the original stellar distribution and velocity dispersion. Each group has properties that are consistent with a star-forming region. We show that our grouped star formation prescription is robust and can be adapted in simulations with varying physical scales and resolution. Such methods are likely to become more important as galactic or even cosmological scale simulations begin to probe sub-parsec scales.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
The formation and early evolution of embedded star clusters in spiral galaxies
Authors:
Steven Rieder,
Clare Dobbs,
Thomas Bending,
Kong You Liow,
James Wurster
Abstract:
We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restri…
▽ More
We present Ekster, a new method for simulating star clusters from birth in a live galaxy simulation that combines the smoothed-particle hydrodynamics (SPH) method Phantom with the $N$-body method PeTar. With Ekster, it becomes possible to simulate individual stars in a simulation with only moderately high resolution for the gas, allowing us to study whole sections of a galaxy rather than be restricted to individual clouds. We use this method to simulate star and star cluster formation in spiral arms, investigating massive GMCs and spiral arm regions with lower mass clouds, from two galaxy models with different spiral potentials. After selecting these regions from pre-run galaxy simulations, we re-sample the particles to obtain a higher resolution. We then re-simulate these regions for 3 Myr to study where and how star clusters form. We analyse the early evolution of the embedded star clusters in these regions. We find that the massive GMC regions, which are more common with stronger spiral arms, form more massive clusters than the sections of spiral arms containing lower mass clouds. Clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.
△ Less
Submitted 18 November, 2021;
originally announced November 2021.
-
Two Directly Imaged, Wide-orbit Giant Planets around the Young, Solar Analog TYC 8998-760-1
Authors:
Alexander J. Bohn,
Matthew A. Kenworthy,
Christian Ginski,
Steven Rieder,
Eric E. Mamajek,
Tiffany Meshkat,
Mark J. Pecaut,
Maddalena Reggiani,
Jozua de Boer,
Christoph U. Keller,
Frans Snik,
John Southworth
Abstract:
Even though tens of directly imaged companions have been discovered in the past decades, the number of directly confirmed multiplanet systems is still small. Dynamical analysis of these systems imposes important constraints on formation mechanisms of these wide-orbit companions. As part of the Young Suns Exoplanet Survey (YSES) we report the detection of a second planetary-mass companion around th…
▽ More
Even though tens of directly imaged companions have been discovered in the past decades, the number of directly confirmed multiplanet systems is still small. Dynamical analysis of these systems imposes important constraints on formation mechanisms of these wide-orbit companions. As part of the Young Suns Exoplanet Survey (YSES) we report the detection of a second planetary-mass companion around the 17 Myr-old, solar-type star TYC 8998-760-1 that is located in the Lower Centaurus Crux subgroup of the Scorpius-Centaurus association. The companion has a projected physical separation of 320 au and several individual photometric measurements from 1.1 to 3.8 microns constrain a companion mass of $6\pm1\,M_\mathrm{Jup}$, which is equivalent to a mass ratio of $q=0.57\pm0.10\%$ with respect to the primary. With the previously detected $14\pm3\,M_\mathrm{Jup}$ companion that is orbiting the primary at 160 au, TYC 8998-760-1 is the first directly imaged multiplanet system that is detected around a young, solar analog. We show that circular orbits are stable, but that mildly eccentric orbits for either/both components ($e > 0.1$) are chaotic on Gyr timescales, implying in-situ formation or a very specific ejection by an unseen third companion. Due to the wide separations of the companions TYC 8998-760-1 is an excellent system for spectroscopic and photometric follow-up with space-based observatories such as the James Webb Space Telescope.
△ Less
Submitted 21 July, 2020;
originally announced July 2020.
-
The formation of young massive clusters by colliding flows
Authors:
C. L. Dobbs,
K. Y. Liow,
S. Rieder
Abstract:
Young massive clusters (YMCs) are the most intense regions of star formation in galaxies. Formulating a model for YMC formation whilst at the same time meeting the constraints from observations is highly challenging however. We show that forming YMCs requires clouds with densities $\gtrsim$ 100 cm$^{-3}$ to collide with high velocities ($\gtrsim$ 20 km s$^{-1}$). We present the first simulations w…
▽ More
Young massive clusters (YMCs) are the most intense regions of star formation in galaxies. Formulating a model for YMC formation whilst at the same time meeting the constraints from observations is highly challenging however. We show that forming YMCs requires clouds with densities $\gtrsim$ 100 cm$^{-3}$ to collide with high velocities ($\gtrsim$ 20 km s$^{-1}$). We present the first simulations which, starting from moderate cloud densities of $\sim100$ cm$^{-3}$, are able to convert a large amount of mass into stars over a time period of around 1 Myr, to produce dense massive clusters similar to those observed. Such conditions are commonplace in more extreme environments, where YMCs are common, but atypical for our Galaxy, where YMCs are rare.
△ Less
Submitted 20 April, 2020;
originally announced April 2020.
-
The Cosmic Ballet: spin and shape alignments of haloes in the cosmic web
Authors:
Punyakoti Ganeshaiah Veena,
Marius Cautun,
Rien van de Weygaert,
Elmo Tempel,
Bernard J. T. Jones,
Steven Rieder,
Carlos S. Frenk
Abstract:
We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $Λ$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, wi…
▽ More
We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $Λ$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.
△ Less
Submitted 19 September, 2018; v1 submitted 30 April, 2018;
originally announced May 2018.
-
Dynamical evolution of stars and gas of young embedded stellar sub-clusters
Authors:
Alison Sills,
Steven Rieder,
Jennifer Scora,
Jessica McCloskey,
Sarah Jaffa
Abstract:
We present simulations of the dynamical evolution of young embedded star clusters. Our initial conditions are directly derived from X-ray, infrared, and radio observations of local systems, and our models evolve both gas and stars simultaneously. Our regions begin with both clustered and extended distributions of stars, and a gas distribution which can include a filamentary structure in addition t…
▽ More
We present simulations of the dynamical evolution of young embedded star clusters. Our initial conditions are directly derived from X-ray, infrared, and radio observations of local systems, and our models evolve both gas and stars simultaneously. Our regions begin with both clustered and extended distributions of stars, and a gas distribution which can include a filamentary structure in addition to gas surrounding the stellar subclusters. We find that the regions become spherical, monolithic, and smooth quite quickly, and that the dynamical evolution is dominated by the gravitational interactions between the stars. In the absence of stellar feedback, the gas moves gently out of the centre of our regions but does not have a significant impact on the motions of the stars at the earliest stages of cluster formation. Our models at later times are consistent with observations of similar regions in the local neighbourhood. We conclude that the evolution of young proto-star clusters is relatively insensitive to reasonable choices of initial conditions. Models with more realism, such as an initial population of binary and multiple stars and ongoing star formation, are the next step needed to confirm these findings.
△ Less
Submitted 12 March, 2018;
originally announced March 2018.
-
Constraints on the size and dynamics of the J1407b ring system
Authors:
Steven Rieder,
Matthew A. Kenworthy
Abstract:
Context. J1407 (1SWASP J140747.93-394542.6 in full) is a young star in the Scorpius-Centaurus OB association that underwent a series of complex eclipses over 56 days in 2007. To explain these, it was hypothesised that a secondary substellar companion, J1407b, has a giant ring system filling a large fraction of the Hill sphere, causing the eclipses. Observations have not successfully detected J1407…
▽ More
Context. J1407 (1SWASP J140747.93-394542.6 in full) is a young star in the Scorpius-Centaurus OB association that underwent a series of complex eclipses over 56 days in 2007. To explain these, it was hypothesised that a secondary substellar companion, J1407b, has a giant ring system filling a large fraction of the Hill sphere, causing the eclipses. Observations have not successfully detected J1407b, but do rule out circular orbits for the companion around the primary star.
Aims. We test to what degree the ring model of J1407b could survive in an eccentric orbit required to fit the observations.
Methods. We run N-body simulations under the AMUSE framework to test the stability of Hill radius-filling systems where the companion is on an eccentric orbit.
Results. We strongly rule out prograde ring systems and find that a secondary of 60 to 100MJup with an 11 year orbital period and retrograde orbiting material can survive for at least $10^4$ orbits and produce eclipses with similar durations as the observed one.
△ Less
Submitted 21 November, 2016; v1 submitted 27 September, 2016;
originally announced September 2016.
-
Felix: A Topology based Framework for Visual Exploration of Cosmic Filaments
Authors:
Nithin Shivshankar,
Pratyush Pranav,
Vijay Natarajan,
Rien van de Weygaert,
E G Patrick Bos,
Steven Rieder
Abstract:
The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which s…
▽ More
The large-scale structure of the universe is comprised of virialized blob-like clusters, linear filaments, sheet-like walls and huge near empty three-dimensional voids. Characterizing the large scale universe is essential to our understanding of the formation and evolution of galaxies. The density range of clusters, walls and voids are relatively well separated, when compared to filaments, which span a relatively larger range. The large scale filamentary network thus forms an intricate part of the cosmic web.
In this paper, we describe Felix, a topology based framework for visual exploration of filaments in the cosmic web. The filamentary structure is represented by the ascending manifold geometry of the 2-saddles in the Morse-Smale complex of the density field. We generate a hierarchy of Morse-Smale complexes and query for filaments based on the density ranges at the end points of the filaments. The query is processed efficiently over the entire hierarchical Morse-Smale complex, allowing for interactive visualization.
We apply Felix to computer simulations based on the heuristic Voronoi kinematic model and the standard $Λ$CDM cosmology, and demonstrate its usefulness through two case studies. First, we extract cosmic filaments within and across cluster like regions in Voronoi kinematic simulation datasets. We demonstrate that we produce similar results to existing structure finders. Filaments that form the spine of the cosmic web, which exist in high density regions in the current epoch, are isolated using Felix. Also, filaments present in void-like regions are isolated and visualized. These filamentary structures are often over shadowed by higher density range filaments and are not easily characterizable and extractable using other filament extraction methodologies.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
-
The cosmic web in CosmoGrid void regions
Authors:
Steven Rieder,
Rien van de Weygaert,
Marius Cautun,
Burcu Beygu,
Simon Portegies Zwart
Abstract:
We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$Λ\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous…
▽ More
We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$Λ\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.
△ Less
Submitted 5 November, 2014;
originally announced November 2014.
-
MPWide: a light-weight library for efficient message passing over wide area networks
Authors:
Derek Groen,
Steven Rieder,
Simon Portegies Zwart
Abstract:
We present MPWide, a light weight communication library which allows efficient message passing over a distributed network. MPWide has been designed to connect application running on distributed (super)computing resources, and to maximize the communication performance on wide area networks for those without administrative privileges. It can be used to provide message-passing between application, mo…
▽ More
We present MPWide, a light weight communication library which allows efficient message passing over a distributed network. MPWide has been designed to connect application running on distributed (super)computing resources, and to maximize the communication performance on wide area networks for those without administrative privileges. It can be used to provide message-passing between application, move files, and make very fast connections in client-server environments. MPWide has already been applied to enable distributed cosmological simulations across up to four supercomputers on two continents, and to couple two different bloodflow simulations to form a multiscale simulation.
△ Less
Submitted 3 December, 2013;
originally announced December 2013.
-
Evolution of star clusters in a cosmological tidal field
Authors:
Steven Rieder,
Tomoaki Ishiyama,
Paul Langelaan,
Junichiro Makino,
Stephen L. W. McMillan,
Simon Portegies Zwart
Abstract:
We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializ…
▽ More
We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment.
We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.
△ Less
Submitted 27 September, 2013;
originally announced September 2013.
-
Assembly of filamentary void galaxy configurations
Authors:
Steven Rieder,
Rien van de Weygaert,
Marius Cautun,
Burcu Beygu,
Simon Portegies Zwart
Abstract:
We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed…
▽ More
We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.
△ Less
Submitted 30 July, 2013; v1 submitted 26 July, 2013;
originally announced July 2013.
-
High performance cosmological simulations on a grid of supercomputers
Authors:
Derek Groen,
Steven Rieder,
Simon Portegies Zwart
Abstract:
We present results from our cosmological N-body simulation which consisted of 2048x2048x2048 particles and ran distributed across three supercomputers throughout Europe. The run, which was performed as the concluding phase of the Gravitational Billion Body Problem DEISA project, integrated a 30 Mpc box of dark matter using an optimized Tree/Particle Mesh N-body integrator. We ran the simulation up…
▽ More
We present results from our cosmological N-body simulation which consisted of 2048x2048x2048 particles and ran distributed across three supercomputers throughout Europe. The run, which was performed as the concluding phase of the Gravitational Billion Body Problem DEISA project, integrated a 30 Mpc box of dark matter using an optimized Tree/Particle Mesh N-body integrator. We ran the simulation up to the present day (z=0), and obtained an efficiency of about 0.93 over 2048 cores compared to a single supercomputer run. In addition, we share our experiences on using multiple supercomputers for high performance computing and provide several recommendations for future projects.
△ Less
Submitted 26 September, 2011;
originally announced September 2011.
-
The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos
Authors:
Tomoaki Ishiyama,
Steven Rieder,
Junichiro Makino,
Simon Portegies Zwart,
Derek Groen,
Keigo Nitadori,
Cees de Laat,
Stephen McMillan,
Kei Hiraki,
Stefan Harfst
Abstract:
We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10…
▽ More
We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10^7 Msun. We have analyzed the spherically averaged density profiles of the three most massive halos which are of galaxy group size and contain at least 170 million particles. The slopes of these density profiles become shallower than -1 at the inner most radius. We also find a clear correlation of halo concentration with mass. The mass dependence of the concentration parameter cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory proposed by Navarro et al. gives reasonable agreement with this dependence. The spin parameter does not show a correlation with the halo mass. The probability distribution functions for both concentration and spin are well fitted by the log-normal distribution for halos with the masses larger than ~10^8 Msun. The subhalo abundance depends on the halo mass. Galaxy-sized halos have 50% more subhalos than ~10^{11} Msun halos have.
△ Less
Submitted 8 April, 2013; v1 submitted 10 January, 2011;
originally announced January 2011.
-
A Light-Weight Communication Library for Distributed Computing
Authors:
Derek Groen,
Steven Rieder,
Paola Grosso,
Cees de Laat,
Simon Portegies Zwart
Abstract:
We present MPWide, a platform independent communication library for performing message passing between computers. Our library allows coupling of several local MPI applications through a long distance network and is specifically optimized for such communications. The implementation is deliberately kept light-weight, platform independent and the library can be installed and used without administrati…
▽ More
We present MPWide, a platform independent communication library for performing message passing between computers. Our library allows coupling of several local MPI applications through a long distance network and is specifically optimized for such communications. The implementation is deliberately kept light-weight, platform independent and the library can be installed and used without administrative privileges. The only requirements are a C++ compiler and at least one open port to a wide area network on each site. In this paper we present the library, describe the user interface, present performance tests and apply MPWide in a large scale cosmological N-body simulation on a network of two computers, one in Amsterdam and the other in Tokyo.
△ Less
Submitted 16 August, 2010;
originally announced August 2010.
-
Probability distributions and Gleason's Theorem
Authors:
Sebastian Rieder,
Karl Svozil
Abstract:
We discuss concrete examples for frame functions and their associated density operators, as well as for non-Gleason type probability measures.
We discuss concrete examples for frame functions and their associated density operators, as well as for non-Gleason type probability measures.
△ Less
Submitted 8 August, 2006;
originally announced August 2006.