-
Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations
Authors:
M. C. Vergara,
D. R. G. Schleicher,
A. Escala,
B. Reinoso,
F. Flammini Dotti,
A. W. H. Kamlah,
M. Liempi,
N. Hoyer,
N. Neumayer,
R. Spurzem
Abstract:
This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar syst…
▽ More
This paper explores the theoretical relation between star clusters and black holes within, focusing on the potential role of nuclear star clusters (NSCs), globular clusters (GCs), and ultra compact dwarf galaxies (UCDs) as environments that allow for black hole formation via stellar collisions. This study aims to identify the optimal conditions for stellar collisions across a range of stellar systems leading to the formation of very massive stars that subsequently collapse into black holes. We analyze data from numerical simulations and observations of diverse stellar systems, encompassing various initial conditions, initial mass functions, and evolution scenarios. We computed a critical mass, determined by the interplay of collision time, system age, and initial properties of the star cluster. The efficiency of black hole formation ($ε_{\mathrm{BH}}$) is defined as the ratio of initial stellar mass divided by critical mass. We find that stellar systems with a ratio of initial stellar mass over critical mass above 1 exhibit a high efficiencies in terms of black hole formation, ranging from $30-100\%$. While there is some scatter, potentially attributed to complex system histories and the presence of gas, the results highlight the potential for achieving high efficiencies via a purely collisional channel in black hole formation. In conclusion, this theoretical exploration elucidates the connection between star clusters and black hole formation. The study underscores the significance of UCDs, GCs, and NSCs as environments conducive to the black hole formation scenario via stellar collisions. The defined black hole formation efficiency ($ε_{\mathrm{BH}}$) is shown to be influenced by the ratio of the initial stellar mass to the critical mass.
△ Less
Submitted 11 June, 2024; v1 submitted 20 May, 2024;
originally announced May 2024.
-
Hot Jupiter Formation in Dense Star Clusters
Authors:
Leonard Benkendorff,
Francesco Flammini Dotti,
Katja Stock,
Maxwell Xu Cai,
Rainer Spurzem
Abstract:
Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-…
▽ More
Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-body codes: NBODY6++GPU for the dynamics of dense star clusters with 32 000 and 64 000 stellar members and LonelyPlanets used to follow 200 identical planetary systems around solar mass stars in those star clusters. We use different sets with 3, 4, or 5 planets and with the innermost planet at a semi-major axis of 5 au or 1 au and follow them for 100 Myr in our simulations. The results indicate that HJs are generated with high efficiency in dense star clusters if the innermost planet is already close to the host star at a semi-major axis of 1 au. If the innermost planet is initially beyond a semi-major axis of 5 au, the probability of a potential HJ ranges between $1.5-4.5$ percent. Very dense stellar neighborhoods tend to eject planets rather than forming HJs. A correlation between HJ formation and angular momentum deficit (AMD) is not witnessed. Young Hot Jupiters ($t_{\rm age} < 100$ Myrs) have only been found, in our simulations, in planetary systems with the innermost planet at a semi-major axis of 1 au.
△ Less
Submitted 21 January, 2024;
originally announced January 2024.
-
The dynamical evolution of protoplanetary disks and planets in dense star clusters
Authors:
Francesco Flammini Dotti,
Roberto Capuzzo-Dolcetta,
M. B. N. Kouwenhoven
Abstract:
Most stars are born in dense stellar environments where the formation and early evolution of planetary systems may be significantly perturbed by encounters with neighbouring stars. To investigate on the fate of circumstellar gas disks and planets around young stars dense stellar environments, we numerically evolve star-disk-planet systems. We use the $N$-body codes NBODY6++GPU and SnIPES for the d…
▽ More
Most stars are born in dense stellar environments where the formation and early evolution of planetary systems may be significantly perturbed by encounters with neighbouring stars. To investigate on the fate of circumstellar gas disks and planets around young stars dense stellar environments, we numerically evolve star-disk-planet systems. We use the $N$-body codes NBODY6++GPU and SnIPES for the dynamical evolution of the stellar population, and the SPH-based code GaSPH for the dynamical evolution of protoplanetary disks. The secular evolution of a planetary system in a cluster differs from that of a field star. Most stellar encounters are tidal, adiabatic and nearly-parabolic. The parameters that characterize the impact of an encounter include the orientation of the protoplanetary disk and planet relative to the orbit of the encountering star, and the orbital phase and the semi-major axis of the planet. We investigate this dependence for close encounters ($r_p/a\leq 100$, where $r_p$ is the periastron distance of the encountering star and $a$ is the semi-major axis of the planet). We also investigate distant perturbers ($r_p/a\gg 100$), which have a moderate effect on the dynamical evolution of the planet and the protoplanetary disk. We find that the evolution of protoplanetary disks in star clusters differs significantly from that of isolated systems. When interpreting the outcome of the planet formation process, it is thus important to consider their birth environments.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
3-D selection of 167 sub-stellar companions to nearby stars
Authors:
Fabo Feng,
R. Paul Butler,
Steven S. Vogt,
Matthew S. Clement,
C. G. Tinney,
Kaiming Cui,
Masataka Aizawa,
Hugh R. A. Jones,
J. Bailey,
Jennifer Burt,
B. D. Carter,
Jeffrey D. Crane,
Francesco Flammini Dotti,
Bradford Holden,
Bo Ma,
Masahiro Ogihara,
Rebecca Oppenheimer,
S. J. O'Toole,
Stephen A. Shectman,
Robert A. Wittenmyer,
Sharon X. Wang,
D. J. Wright,
Yifan Xuan
Abstract:
We analyze 5108 AFGKM stars with at least five high precision radial velocity points as well as Gaia and Hipparcos astrometric data utilizing a novel pipeline developed in previous work. We find 914 radial velocity signals with periods longer than 1000\,d. Around these signals, 167 cold giants and 68 other types of companions are identified by combined analyses of radial velocity, astrometry, and…
▽ More
We analyze 5108 AFGKM stars with at least five high precision radial velocity points as well as Gaia and Hipparcos astrometric data utilizing a novel pipeline developed in previous work. We find 914 radial velocity signals with periods longer than 1000\,d. Around these signals, 167 cold giants and 68 other types of companions are identified by combined analyses of radial velocity, astrometry, and imaging data. Without correcting for detection bias, we estimate the minimum occurrence rate of the wide-orbit brown dwarfs to be 1.3\%, and find a significant brown dwarf valley around 40 $M_{\rm Jup}$. We also find a power-law distribution in the host binary fraction beyond 3 au similar to that found for single stars, indicating no preference of multiplicity for brown dwarfs. Our work also reveals nine sub-stellar systems (GJ 234 B, GJ 494 B, HD 13724 b, HD 182488 b, HD 39060 b and c, HD 4113 C, HD 42581 d, HD 7449 B, and HD 984 b) that have previously been directly imaged, and many others that are observable at existing facilities. Depending on their ages we estimate that an additional 10-57 sub-stellar objects within our sample can be detected with current imaging facilities, extending the imaged cold (or old) giants by an order of magnitude.
△ Less
Submitted 13 September, 2022; v1 submitted 26 August, 2022;
originally announced August 2022.
-
The impact of stellar evolution on rotating star clusters: the gravothermal-gravogyro catastrophe and the formation of a bar of black holes
Authors:
A. W. H. Kamlah,
R. Spurzem,
P. Berczik,
M. Arca Sedda,
F. Flammini Dotti,
N. Neumayer,
X. Pang,
Q. Shu,
A. Tanikawa,
M. Giersz
Abstract:
We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact o…
▽ More
We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact of rotation and stellar evolution on the star cluster dynamics. In all runs for rotating star clusters we detect a previously predicted mechanism: an initial phase of violent relaxation followed by the so-called gravogyro catastrophe. We find that the gravogyro catastrophe reaches a finite amplitude, which depends in strength on the level of the bulk rotation, and then levels off. After this phase the angular momentum is transferred from high-mass to low-mass particles in the cluster (both stars and compact objects). Simultaneously, the system becomes gravothermally unstable and collapses, thus undergoing the so-called gravothermal-gravogyro catastrophe. Comparing models with and without stellar evolution, we find an interesting difference. When stellar evolution is not taken into account, the whole process proceeds at a faster pace. The population of heavy objects tend to form a triaxial structure that rotates in the cluster centre. When stellar evolution is taken into account, we find that such a {\it rotating bar} is populated by stellar black holes and their progenitors. The triaxial structure becomes axisymmetric over time, but we also find that the models without stellar evolution suffer repeated gravogyro catastrophes as sufficient angular momentum and mass are removed by the tidal field.
△ Less
Submitted 31 August, 2022; v1 submitted 9 May, 2022;
originally announced May 2022.
-
A highly mutually-inclined, compact warm-Jupiter system KOI-984 ?
Authors:
Leilei Sun,
Panagiotis Ioannidis,
Shenghong Gu,
Jürgen H. M. M. Schmitt,
Xiaobin Wang,
M. B. N. Kouwenhoven,
Volker Perdelwitz,
Francesco Flammini Dotti,
Stefan Czesla
Abstract:
The discovery of a population of close-orbiting giant planets ($\le$ 1 au) has raised a number of questions about their origins and dynamical histories. These issues have still not yet been fully resolved, despite over 20 years of exoplanet detections and a large number of discovered exoplanets. In particular, it is unclear whether warm Jupiters (WJs) form in situ, or whether they migrate from fur…
▽ More
The discovery of a population of close-orbiting giant planets ($\le$ 1 au) has raised a number of questions about their origins and dynamical histories. These issues have still not yet been fully resolved, despite over 20 years of exoplanet detections and a large number of discovered exoplanets. In particular, it is unclear whether warm Jupiters (WJs) form in situ, or whether they migrate from further outside and are even currently migrating to form hot Jupiters (HJs). Here, we report the possible discovery and characterization of the planets in a highly mutually-inclined ($I_{\rm mut}\simeq 45^\circ$), compact two-planet system (KOI-984), in which the newly discovered warm Jupiter KOI-984$c$ is on a 21.5-day, moderately eccentric ($e\simeq 0.4$) orbit, in addition to a previously known 4.3-day planet candidate KOI-984$b$. Meanwhile, the orbital configuration of a moderately inclined ($I_{\rm mut}\simeq 15^\circ$), low-mass ($m_{c}\simeq 24 M_{\oplus}$;$P_b\simeq 8.6$ days) perturbing planet near 1:2 mean motion resonace with KOI-984$b$ could also well reproduce observed transit timing variations and transit duration variations of KOI-984$b$. Such an eccentric WJ with a close-in sibling would pose a challenge to the proposed formation and migration mechanisms of WJs, if the first scenario is supported with more evidences in near future; this system with several other well-measured inclined WJ systems (e.g., Kepler-419 and Kepler-108) may provide additional clues for the origin and dynamical histories of WJs.
△ Less
Submitted 22 November, 2021; v1 submitted 18 November, 2021;
originally announced November 2021.
-
The long-term evolution of main-sequence binaries in DRAGON simulations
Authors:
Qi Shu,
Xiaoying Pang,
Francesco Flammini Dotti,
M. B. N. Kouwenhoven,
Manuel Arca Sedda,
Rainer Spurzem
Abstract:
We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have…
▽ More
We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have different initial stellar mass functions: Kroupa 1993 (D1-R7-IMF93) and Kroupa 2001 (D2-R7-IMF01); and different initial mass ratio distributions: random paring (D1-R7-IMF93) and a power-law (D1-R7-IMF93). In general, the mass ratio of a population of binaries increases over time due to stellar evolution, which is less significant in D2-R7-IMF01. In D1-R7-IMF93, primordial binaries with mass ratio $q \approx$ 0.2 are most common, and the frequency linearly declines with increasing $q$ at all times. Dynamical binaries of both models have higher eccentricities and larger semi-major axes than primordial binaries. They are preferentially located in the inner part of the star cluster. Secular evolution of binary orbital parameters does not depend on the initial mass-ratio distribution, but is sensitive to the initial binary distribution of the system. At t = 12 Gyr, the binary fraction decreases radially outwards, and mass segregation is present. A color difference of 0.1 mag in $F330W-F814W$ and 0.2 mag in $NUV-y$ between the core and the outskirts of both clusters is seen, which is a reflection of the binary radial distribution and the mass segregation in the cluster. The complete set of data for primordial and dynamical binary systems at all snapshot intervals is made publicly available.
△ Less
Submitted 30 November, 2020;
originally announced November 2020.
-
Planetary Systems in a Star Cluster II: intermediate-mass black holes and planetary systems
Authors:
Francesco Flammini Dotti,
M. B. N. Kouwenhoven,
Qi Shu,
Wei Hao,
Rainer Spurzem
Abstract:
Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of si…
▽ More
Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of single-planet systems and free-floating planets, over a period of 100~Myr, in star clusters without an IMBH, and in clusters with a central IMBH of mass $100~M_\odot$ or $200~M_\odot$. In the central region ($r\lesssim 0.2$~pc), the IMBH's tidal influence on planetary systems is typically 10~times stronger than the average neighbour star. For a star cluster with a $200M_\odot$ IMBH, the region in which the IMBH's influence is stronger within the virial radius ($\sim 1$~pc). The IMBH quenches mass segregation, and the stars in the core tend to move towards intermediate regions. The ejection rate of both stars and planets is higher when an IMBH is present. The rate at which planets are expelled from their host star rate is higher for clusters with higher IMBH masses, for $t<0.5 t_{rh}$, while remains mostly constant while the star cluster fills its Roche lobe, similar to a star cluster without an IMBH. The disruption rate of planetary systems is higher in initially denser clusters, and for wider planetary orbits, but this rate is substantially enhanced by the presence of a central IMBH.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
Linking the formation and fate of exo-Kuiper belts within solar system analogues
Authors:
Dimitri Veras,
Katja Reichert,
Francesco Flammini Dotti,
Maxwell X. Cai,
Alexander J. Mustill,
Andrew Shannon,
Catriona H. McDonald,
Simon Portegies Zwart,
M. B. N. Kouwenhoven,
Rainer Spurzem
Abstract:
Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially em…
▽ More
Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
Planetary systems in a star cluster I: the Solar system scenario
Authors:
Francesco Flammini Dotti,
M. B. N. Kouwenhoven,
Maxwell Xu Cai,
Rainer Spurzem
Abstract:
Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planeta…
▽ More
Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planetary systems escape from the star cluster within 50 Myr; those with low escape speeds often remain intact during and after the escape process. While most planetary systems inside the star cluster remain intact, a subset is strongly perturbed during the first 50 Myr. Over the course of time, 0.3 % - 5.3 % of the planets escape, sometimes up to tens of millions of years after a stellar encounter occurred. Survival rates are highest for Jupiter, while Uranus and Neptune have the highest escape rates. Unless directly affected by a stellar encounter itself, Jupiter frequently serves as a barrier that protects the terrestrial planets from perturbations in the outer planetary system. In low-density environments, Jupiter provides protection from perturbations in the outer planetary system, while in high-density environments, direct perturbations of Jupiter by neighbouring stars is disruptive to habitable-zone planets. The diversity amongst planetary systems that is present in the star clusters at 50 Myr, and amongst the escaping planetary systems, is high, which contributes to explaining the high diversity of observed exoplanet systems in star clusters and in the Galactic field
△ Less
Submitted 21 August, 2019;
originally announced August 2019.
-
SPHERE dynamical and spectroscopic characterization of HD142527B
Authors:
R. Claudi,
A. -L. Maire,
D. Mesa,
A. Cheetham,
C. Fontanive,
R. Gratton,
A. Zurlo,
H. Avenhaus,
T. Bhowmik,
B. Biller,
A. Boccaletti,
M. Bonavita,
M. Bonnefoy,
E. Cascone,
G. Chauvin,
A. Delboulbè,
S. Desidera,
V. D'Orazi,
P. Feautrier,
M. Feldt,
F. Flammini Dotti,
J. H. Girard,
E. Giro,
M. Janson,
J. Hagelberg
, et al. (28 additional authors not shown)
Abstract:
We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty o…
▽ More
We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty of one spectral subtype, compatible with an object of M=0.11 \pm 0.06 MSun and R=0.15 \pm 0.07 RSun. The determination of the mass remains a challenge using contemporary evolutionary models, as they do not account for the energy input due to accretion from infalling material. We consider that the spectral type of the secondary may also be earlier than the type we derived from IFS spectra. From dynamical considerations, we further constrain the mass to 0.26^{+0.16}_{-0.14} MSun , which is consistent with both our spectroscopic analysis and the values reported in the literature. Following previous methods, the lower and upper dynamical mass values correspond to a spectral type between M2.5 and M5.5 for the companion. By fitting the astrometric points, we find the following orbital parameters: a period of P=35-137 yr; an inclination of i=121-130 deg.; , a value of Omega=124-135 deg for the longitude of node, and an 68% confidence interval of \sim 18 - 57 au for the separation at periapsis. Eccentricity and time at periapsis passage exhibit two groups of values: \sim0.2-0.45 and \sim0.45-0.7 for e, and \sim 2015-2020 and \sim2020-2022 for T_0. While these orbital parameters might at first suggest that HD142527B is not the companion responsible for the outer disk truncation, a previous hydrodynamical analysis of this system showed that they are compatible with a companion that is able to produce the large cavity and other observed features.
△ Less
Submitted 19 December, 2018;
originally announced December 2018.
-
Planetary Systems in Star Clusters: the dynamical evolution and survival
Authors:
Francesco Flammini Dotti,
Maxwell Xu Cai,
Rainer Spurzem,
M. B. N. Kouwenhoven
Abstract:
Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary co…
▽ More
Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary companions. To understand the origin and dynamical evolution of such exoplanet systems, it is necessary to carefully study the effect of their environments. Here, we combine theoretical estimates with state-of-the-art numerical simulations of evolving planetary systems similar to our own solar system in different star cluster environments. We combine the REBOUND planetary system evolution code, and the NBODY6++GPU star cluster evolution code, integrated in the AMUSE multi-physics environment. With our study we can constrain the effect of external perturbations of different environments on the planets and debris structures of a wide variety of planetary systems, which may play a key role for the habitability of exoplanets in the Universe.
△ Less
Submitted 30 November, 2018;
originally announced November 2018.