Showing 1–50 of 247 results for author: Zwart, S P

Binary-single interactions with different mass ratios

Authors: Bruno Rando Forastier, Daniel Marín Pina, Mark Gieles, Simon Portegies Zwart, Fabio Antonini

Abstract: Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass… ▽ More Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass ratios of the binary components ($q_2\equiv m_2/m_1\le1$) and the incoming single ($q_3\equiv m_3/m_1$). We extract cross sections and rates for (i) GW capture during resonant interactions; (ii) GW inspiral in between resonant interactions and apply the results to different globular cluster conditions. We find that GW capture during resonant interactions is most efficient if $q_2\simeq q_3$ and that the mass-ratio distribution of BBH coalescence due to inspirals is $\propto m_1^{-1}q^{2.9+α}$, where $α$ is the exponent of the BH mass function. The total rate of GW captures and inspirals depends mostly on $m_1$ and is relatively insensitive to $q_2$ and $q_3$. We show that eccentricity increase by non-resonant encounters approximately doubles the rate of BBH inspiral in between resonant encounters. For a given GC mass and radius, the BBH merger rate in metal-rich GCs is approximately double that of metal-poor GCs, because of their (on average) lower BH masses ($m_1$) and steeper BH mass function, yielding binaries with lower $q$. Our results enable the translating from the mass-ratio distribution of dynamically formed BBH mergers to the underlying BH mass function. The additional mechanism that leads to a doubling of the inspirals provides an explanation for the reported high fraction of in-cluster inspirals in $N$-body models of clusters. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 14 pages, 7 figures

Massive Star Cluster Formation II. Runaway Stars as Fossils of Sub-Cluster Mergers

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Simon Portegies Zwart, Eric P. Andersson, Sabrina M. Appel, Claude Cournoyer-Cloutier, Simon C. O. Glover, Stephen L. W. McMillan

Abstract: Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway star… ▽ More Two main mechanisms have classically been proposed for the formation of runaway stars. The binary supernova scenario (BSS) suggests that a massive star in a binary explodes as a supernova, ejecting its companion. The dynamical ejection scenario suggests that a star is ejected through a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the sub-cluster ejection scenario (SCES), where an infalling sub-cluster of stars is ejected out of the cluster by a slingshot interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of a young massive cluster forming from a $10^6\rm~M_\odot$ gas cloud using the Torch framework. This star cluster forms hierarchically through a sequence of sub-cluster mergers, determined by the initial turbulent, spherical initial conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming sub-clusters fall into the central potential, where they are ejected on a slingshot trajectory, forming groups of runaway stars that are distributed highly anisotropically. Runaways formed by the same SCES share similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling sub-clusters were tightened through ejection. This allows for another velocity kick via subsequent BSS, which is a promising avenue for producing hypervelocity stars unbound to the Galaxy. The SCES occurs when sub-cluster formation is resolved. We expect non-spherical initial gas distributions to increase runaway star numbers. The observation of groups of runaway stars formed via SCES thus reveals the assembly history of their natal clusters. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 11 pages, 10 figures. Submitted to A&A. Comments welcome!

Isles of regularity in a sea of chaos amid the gravitational three-body problem

Authors: Alessandro Alberto Trani, Nathan W. C. Leigh, Tjarda C. N. Boekholt, Simon Portegies Zwart

Abstract: The three-body problem (3BP) poses a longstanding challenge in physics and celestial mechanics. Despite the impossibility of obtaining general analytical solutions, statistical theories have been developed based on the ergodic principle. This assumption is justified by chaos, which is expected to fully mix the accessible phase space of the 3BP. This study probes the presence of regular (i.e. non c… ▽ More The three-body problem (3BP) poses a longstanding challenge in physics and celestial mechanics. Despite the impossibility of obtaining general analytical solutions, statistical theories have been developed based on the ergodic principle. This assumption is justified by chaos, which is expected to fully mix the accessible phase space of the 3BP. This study probes the presence of regular (i.e. non chaotic) trajectories within the 3BP and assesses their impact on statistical escape theories. Using numerical simulations, we establish criteria for identifying regular trajectories and analyse their impact on statistical outcomes. Our analysis reveals that regular trajectories occupy up to 32% of the phase space, and their outcomes defy the predictions of statistical escape theories. The coexistence of regular and chaotic regions at all scales is characterized by a multi-fractal behaviour. Integration errors manifest as numerical chaos, artificially enhancing the mixing of the phase space and affecting the reliability of individual simulations, yet preserving the statistical correctness of an ensemble of realizations. Our findings underscore the challenges in applying statistical escape theories to astrophysical problems, as they may bias results by excluding the outcome of regular trajectories. This is particularly important in the context of formation scenarios of gravitational wave mergers, where biased estimates of binary eccentricity can significantly impact estimates of coalescence efficiency and detectable eccentricity. △ Less

Submitted 3 July, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 15 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Magnetospheric Flows in X-ray Pulsars I: Instability at super-Eddington regime of accretion

Authors: A. A. Mushtukov, A. Ingram, V. F. Suleimanov, N. DiLullo, M. Middleton, S. S. Tsygankov, M. van der Klis, S. Portegies Zwart

Abstract: Within the magnetospheric radius, the geometry of accretion flow in X-ray pulsars is shaped by a strong magnetic field of a neutron star. Starting at the magnetospheric radius, accretion flow follows field lines and reaches the stellar surface in small regions located close to the magnetic poles of a star. At low mass accretion rates, the dynamic of the flow is determined by gravitational attracti… ▽ More Within the magnetospheric radius, the geometry of accretion flow in X-ray pulsars is shaped by a strong magnetic field of a neutron star. Starting at the magnetospheric radius, accretion flow follows field lines and reaches the stellar surface in small regions located close to the magnetic poles of a star. At low mass accretion rates, the dynamic of the flow is determined by gravitational attraction and rotation of the magnetosphere due to the centrifugal force. At the luminosity range close to the Eddington limit and above it, the flow is additionally affected by the radiative force. We construct a model simulating accretion flow dynamics over the magnetosphere, assuming that the flow strictly follows field lines and is affected by gravity, radiative and centrifugal forces only. The magnetic field of a NS is taken to be dominated by the dipole component of arbitrary inclination with respect to the accretion disc plane. We show that accretion flow becomes unstable at high mass accretion rates and tends to fluctuate quasi-periodically with a typical period comparable to the free-fall time from the inner disc radius. The inclination of a magnetic dipole with respect to the disc plane and strong anisotropy of X-ray radiation stabilise the mass accretion rate at the poles of a star, but the surface density of material covering the magnetosphere fluctuates even in this case. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 14 pages, 13 figures, submitted to MNRAS

Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars

Authors: Brooke Polak, Mordecai-Mark Mac Low, Ralf S. Klessen, Jia Wei Teh, Claude Cournoyer-Cloutier, Eric P. Andersson, Sabrina M. Appel, Aaron Tran, Sean C. Lewis, Maite J. C. Wilhelm, Simon Portegies Zwart, Simon C. O. Glover, Long Wang, Stephen L. W. McMillan

Abstract: The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg… ▽ More The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Submitted to A&A

The Steady State of Intermediate-Mass Black Holes Near a Supermassive Black Hole

Authors: Erwan Hochart, Simon Portegies Zwart

Abstract: Aims: Investigate properties of a cluster of intermediate-mass black holes surrounding a supermassive black hole. Methods: We simulate clusters of equal-mass intermediate-mass black holes ($m_{\rm{IMBH}} = 10^{3}$ ${\rm{M_\odot}}$) initialised in a shell between $0.15\leq r$ [pc] $\leq 0.25$ centered about a supermassive black hole. We explore the influence of the cluster population and supermassi… ▽ More Aims: Investigate properties of a cluster of intermediate-mass black holes surrounding a supermassive black hole. Methods: We simulate clusters of equal-mass intermediate-mass black holes ($m_{\rm{IMBH}} = 10^{3}$ ${\rm{M_\odot}}$) initialised in a shell between $0.15\leq r$ [pc] $\leq 0.25$ centered about a supermassive black hole. We explore the influence of the cluster population and supermassive black hole mass on the merger rate, the ejection rate and the escape velocity. For $M_{\text{SMBH}} = 4\times10^{6}$ ${\rm {M}_\odot}$, we use both a Newtonian and post-Newtonian formalism, going up to the 2.5th order and including cross-terms. For the other two SMBH masses ($M_{\rm{SMBH}} = 4\times10^{5}$ ${\rm{M_\odot}}$ and $M_{\rm{SMBH}} = 4\times10^{7}$ $\rm{M_\odot}$), we model the system only taking into account relativistic effects. The simulations end once a black hole escapes the cluster, a merger occurs, or the system has evolved till $100$ Myr. Results: The post-Newtonian formalism accelerates the loss rate of intermediate-mass black holes. Ejections occur more often for lower supermassive black hole masses while more massive ones increase the rate of mergers. Although relativistic effects allow for circularisation, all merging binaries have $e \gtrsim 0.97$. Strong gravitational wave signals are suppressed during our Newtonian calculations. Weaker and more frequent signals are expected from gravitational wave radiation emitted in a fly-by. In our post-Newtonian calculations, $30/406$ of the gravitational wave events capable of being observed with LISA and $μ$Ares are detected as gravitational wave capture binaries with the remaining being in-cluster mergers. Throughout our investigation, no IMBH-IMBH binaries were detected. △ Less

Submitted 26 February, 2024; v1 submitted 8 December, 2023; originally announced December 2023.

Comments: 18 pages, 18 figures

Journal ref: A&A 685, A123 (2024)

The origin and evolution of wide Jupiter Mass Binary Objects in young stellar clusters

Authors: Simon Portegies Zwart, Erwan Hochart

Abstract: The recently observed population of 540 free-floating Jupiter-mass objects, including 40 dynamically soft pairs in the Trapezium cluster have raised interesting questions on their formation and evolution. We test various scenarios for the origin and survivability of these free floating Jupiter-mass objects and Jupiter-mass Binary Objects (JuMBOs) in the Trapezium cluster. The numerical calculation… ▽ More The recently observed population of 540 free-floating Jupiter-mass objects, including 40 dynamically soft pairs in the Trapezium cluster have raised interesting questions on their formation and evolution. We test various scenarios for the origin and survivability of these free floating Jupiter-mass objects and Jupiter-mass Binary Objects (JuMBOs) in the Trapezium cluster. The numerical calculations are performed by direct N-body integration of the stars and planets in the Trapezium cluster starting with a wide variety of planets in various configurations. We discuss four models: SPP, in which selected stars have two outer orbiting Jupiter-mass planets; SPM, where selected stars are orbited by Jupiter-mass planet-moon pairs; ISF in which JuMBOs form in situ with the stars, and FFC, where we introduce a population of free-floating single Jupiter-mass objects, but no initialized binaries. Models FFC and SPP fail to produce enough JuMBOs. Models SPM can produce sufficient JuMBOs, but requires unusually wide orbits for the planet-moon system around the star. The observed JuMBOs and free-floating Jupiter-mass objects in the Trapezium cluster are best reproduced if they formed in pairs and as free-floaters together with the other stars in a smooth (Plummer) density profile with a virial radius of 0.5pc. A fractal stellar distribution also works, but requires relatively recent formations (>0.2Myr after the other stars formed) or a high (50%) initial binary fraction. This would make the primordial binary fraction of JuMBOs even higher than the already large observation fraction of 8%. The fraction of JuMBOs will continue to drop with time, and the lack of JuMBOs in Upper Scorpius could then result in its higher age, causing more JuMBOs to be ionized. We then also predict that the interstellar density of Jupiter-mass objects (mostly singles with 2% lucky surviving binaries) is 0.05/pc$^{3}$. △ Less

Submitted 12 March, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted for publiation in golden open-access journal SciPost Astronomy, see https://scipost.org/submission/scipost_202312_00003v1/ Software available at: https://github.com/spzwart/JuMBOs.git Simulation data available at: https://zenodo.org/records/10149241

A hybrid approach for solving the gravitational N-body problem with Artificial Neural Networks

Authors: Veronica Saz Ulibarrena, Philipp Horn, Simon Portegies Zwart, Elena Sellentin, Barry Koren, Maxwell X. Cai

Abstract: Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural Networks (ANNs) to replace expensive parts of the integration of planetary systems. Neural networks that include physical knowledge have grown in popularity in the l… ▽ More Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural Networks (ANNs) to replace expensive parts of the integration of planetary systems. Neural networks that include physical knowledge have grown in popularity in the last few years, although few attempts have been made to use them to speed up the simulation of the motion of celestial bodies. We study the advantages and limitations of using Hamiltonian Neural Networks to replace computationally expensive parts of the numerical simulation. We compare the results of the numerical integration of a planetary system with asteroids with those obtained by a Hamiltonian Neural Network and a conventional Deep Neural Network, with special attention to understanding the challenges of this problem. Due to the non-linear nature of the gravitational equations of motion, errors in the integration propagate. To increase the robustness of a method that uses neural networks, we propose a hybrid integrator that evaluates the prediction of the network and replaces it with the numerical solution if considered inaccurate. Hamiltonian Neural Networks can make predictions that resemble the behavior of symplectic integrators but are challenging to train and in our case fail when the inputs differ ~7 orders of magnitude. In contrast, Deep Neural Networks are easy to train but fail to conserve energy, leading to fast divergence from the reference solution. The hybrid integrator designed to include the neural networks increases the reliability of the method and prevents large energy errors without increasing the computing cost significantly. For this problem, the use of neural networks results in faster simulations when the number of asteroids is >70. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: Accepted for publication in the Journal of Computational Physics

Punctuated chaos and the unpredictability of the Galactic center S-star orbital evolution

Authors: Simon Portegies Zwart, Tjarda Boekholt, Douglas Heggie

Abstract: We investigate the chaotic behavior of the S-star cluster in the Galactic center using precise $N$-body calculations, free from round-off or discretization errors. Our findings reveal that chaos among the Galactic center S-stars arises from close encounters, particularly among pairs and near the massive central body. These encounters induce perturbations, causing sudden changes in the orbital ener… ▽ More We investigate the chaotic behavior of the S-star cluster in the Galactic center using precise $N$-body calculations, free from round-off or discretization errors. Our findings reveal that chaos among the Galactic center S-stars arises from close encounters, particularly among pairs and near the massive central body. These encounters induce perturbations, causing sudden changes in the orbital energies of the interacting stars. Consequently, neighboring solutions experience roughly exponential growth in separation. We propose a theory of "punctuated chaos" that describes the S-star cluster's chaotic behavior. This phenomenon results from nearly linear growth in the separation between neighboring orbits after repeated finite perturbations. Each participating star's orbit experiences discrete, abrupt changes in energy due to the perturbations. The cumulative effect of these events is further amplified by the steady drift in orbital phase. In the Galactic center, perturbations originate from coincidental encounters occurring within a distance of $\aplt 100$\,au between at least two stars (in some cases, three stars). Our model satisfactorily explains the observed exponential growth in the 27 S-star cluster. We determine that the S-star system has a Lyapunov time scale of approximately 462 +/-74 years. For the coming millennium, chaos in the S-star cluster will be driven mainly by a few of the closest orbiting stars: S2, S5, S6, S8, S9, S14, S18, S31, S21, S24, S27, S29, and S38. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS

Modeling the chemical enrichment history of the Bulge Fossil Fragment Terzan 5

Authors: Donatella Romano, Francesco R. Ferraro, Livia Origlia, Simon Portegies Zwart, Barbara Lanzoni, Chiara Crociati, Davide Massari, Emanuele Dalessandro, Alessio Mucciarelli, R. Michael Rich, Francesco Calura, Francesca Matteucci

Abstract: Terzan 5 is a heavily obscured stellar system located in the inner Galaxy. It has been postulated to be a stellar relic, a Bulge Fossil Fragment witnessing the complex history of the assembly of the Milky Way bulge. In this paper, we follow the chemical enrichment of a set of putative progenitors of Terzan 5 to assess whether the chemical properties of this cluster fit within a formation scenario… ▽ More Terzan 5 is a heavily obscured stellar system located in the inner Galaxy. It has been postulated to be a stellar relic, a Bulge Fossil Fragment witnessing the complex history of the assembly of the Milky Way bulge. In this paper, we follow the chemical enrichment of a set of putative progenitors of Terzan 5 to assess whether the chemical properties of this cluster fit within a formation scenario in which it is the remnant of a primordial building block of the bulge. We can explain the metallicity distribution function and the runs of different element-to-iron abundance ratios as functions of [Fe/H] derived from optical-infrared spectroscopy of giant stars in Terzan 5, by assuming that the cluster experienced two major star formation bursts separated by a long quiescent phase. We further predict that the most metal-rich stars in Terzan 5 are moderately He-enhanced and a large spread of He abundances in the cluster, Y = 0.26-0.335. We conclude that current observations fit within a formation scenario in which Terzan 5 originated from a pristine, or slightly metal-enriched, gas clump about one order of magnitude more massive than its present-day mass. Losses of gas and stars played a major role in shaping Terzan 5 the way we see it now. The iron content of the youngest stellar population is better explained if the white dwarfs that give rise to type Ia supernovae (the main Fe factories) sink towards the cluster center, rather than being stripped by the strong tidal forces exerted by the Milky Way in the outer regions. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 14 pages, 2 tables, 8 figures, accepted for publication in ApJ

Early Evolution and 3D Structure of Embedded Star Clusters

Authors: Claude Cournoyer-Cloutier, Alison Sills, William E. Harris, Sabrina M. Appel, Sean C. Lewis, Brooke Polak, Aaron Tran, Martijn J. C. Wilhelm, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart

Abstract: We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at… ▽ More We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: 15 pages, 10 figures, to be published in MNRAS

Radiation shielding of protoplanetary discs in young star-forming regions

Authors: Martijn J. C. Wilhelm, Simon Portegies Zwart, Claude Cournoyer-Cloutier, Sean C. Lewis, Brooke Polak, Aaron Tran, Mordecai-Mark Mac Low

Abstract: Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two en… ▽ More Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 23 pages, 22 figures, 1 table, accepted for publication in MNRAS

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.

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

Early-Forming Massive Stars Suppress Star Formation and Hierarchical Cluster Assembly

Authors: Sean C. Lewis, Stephen L. W. McMillan, Mordecai-Mark Mac Low, Claude Cournoyer-Cloutier, Brooke Polak, Martijn J. C. Wilhelm, Aaron Tran, Alison Sills, Simon Portegies Zwart, Ralf S. Klessen, Joshua E. Wall

Abstract: Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation… ▽ More Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite \texttt{Torch}, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer, $N$-body dynamics, and stellar feedback to model four initially identical $10^4$ M$_\odot$ giant molecular clouds with a Gaussian density profile peaking at $521.5 \mbox{ cm}^{-3}$. Using the \texttt{Torch} software suite through the \texttt{AMUSE} framework we modify three of the models to ensure that the first star that forms is very massive (50, 70, 100 M$_\odot$). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star forming region. The star formation rate is suppressed, reducing the total mass of stars formed. Our fiducial control model without an early massive star has a larger star formation rate and total efficiency by up to a factor of three and a higher average star formation efficiency per free-fall time by up to a factor of seven. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster. △ Less

Submitted 28 February, 2023; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 17 pages, 7 figures. Published in ApJ

Journal ref: ApJ 944 211 (2023)

Bright X-ray Pulsars: how outflows influence beaming, pulsations and pulse phase lags

Authors: Alexander A. Mushtukov, Simon Portegies Zwart

Abstract: Extreme accretion in X-ray pulsars (XRPs) results in radiation-driven outflows launched from the inner parts of the accretion disc. The outflows affect the apparent luminosity of the XRPs and their pulsations through the geometrical beaming. We model processes of geometrical beaming and pulse formation using Monte Carlo simulations. We confirm our earlier statement that strong amplification of lum… ▽ More Extreme accretion in X-ray pulsars (XRPs) results in radiation-driven outflows launched from the inner parts of the accretion disc. The outflows affect the apparent luminosity of the XRPs and their pulsations through the geometrical beaming. We model processes of geometrical beaming and pulse formation using Monte Carlo simulations. We confirm our earlier statement that strong amplification of luminosity due to the collimation of X-ray photons is inconsistent with a large pulsed fraction. Accounting for relativistic aberration due to possibly high outflow velocity ($\sim 0.2c$) does not affect this conclusion. We demonstrate that the beaming causes phase lags of pulsations. Within the opening angle of the accretion cavity formed by the outflows, phase lags tend to be sensitive to observers viewing angles. Variations in outflow geometry and corresponding changes of the phase lags might influence the detectability of pulsation in bright X-ray pulsars and ULXs. We speculate that the strong geometrical beaming is associated with large radiation pressure on the walls of accretion cavity due to multiple photons reflections. We expect that the mass loss rate limits geometrical beaming: strong beaming becomes possible only under sufficiently large fractional mass loss rate from the disc. △ Less

Submitted 22 November, 2022; v1 submitted 16 November, 2022; originally announced November 2022.

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

The paradox of infinitesimal granularity: Chaos and the reversibility of time in Newton's theory of gravity

Authors: Simon Portegies Zwart, Tjarda Boekholt

Abstract: The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction… ▽ More The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction postulate then introduces fundamental uncertainties in the Cartesian phase space coordinates that propagate through classical three-body dynamics to macroscopic scales within the triple's lifetime. The spontaneous collapse of the wave function then drives the global chaotic behavior of the Universe through the superposition of triple systems (and probably multi-body systems). The paradox of infinitesimal granularity then arises from the superposition principle, which states that any multi-body system is composed of an ensemble of three-body problems. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: Accepted for publication in IC M2 Physical Sciences AIP Conference Proceedings

Report number: C01-Y22-P276

A Monte Carlo study of early gas expulsion and evolution of star clusters: new simulations with the MOCCA code in the AMUSE framework

Authors: A. Leveque, M. Giersz, S. Banerjee, E. Vesperini, J. Hong, S. Portegies Zwart

Abstract: We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in… ▽ More We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 13 pages, 6 figures

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.

Comments: 14 Pages, 8 Figures, accepted for publication in PRD

Statistical features of multiple Compton scattering in a strong magnetic field

Authors: Alexander A. Mushtukov, Ivan D. Markozov, Valery F. Suleimanov, Dmitrij I. Nagirner, Alexander D. Kaminker, Alexander. Y. Potekhin, Simon Portegies Zwart

Abstract: Compton scattering is a key process shaping spectra formation and accretion flow dynamics in accreting strongly magnetized neutron stars. A strong magnetic field affects the scattering cross section and makes it dependent on photon energy, momentum, and polarization state. Using Monte Carlo simulations, we investigate statistical features of Compton scattering of polarized X-ray radiation in a str… ▽ More Compton scattering is a key process shaping spectra formation and accretion flow dynamics in accreting strongly magnetized neutron stars. A strong magnetic field affects the scattering cross section and makes it dependent on photon energy, momentum, and polarization state. Using Monte Carlo simulations, we investigate statistical features of Compton scattering of polarized X-ray radiation in a strong magnetic field. Our analysis is focused on photon gas behaviour well inside the scattering region. We take into account the resonant scattering at the fundamental cyclotron frequency, thermal distribution of electrons at the ground Landau level, and bulk velocity of the electron gas. We show that (i) the photons scattered around the cyclotron energy by the electron gas at rest tend to acquire the final energy close to the cyclotron one with a very small dispersion measure; (ii) the redistribution of photons within the Doppler core of cyclotron resonance differs significantly from the complete redistribution; (iii) the efficiency of momentum transfer from photons to the electron gas is affected by the temperature of electron gas both for photons at cyclotron energy and below it; (iv) the momentum transfer from photons to the electron gas of non-zero bulk velocity is more efficient in the case of magnetic scattering. △ Less

Submitted 26 May, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 18 pages, 13 figures, accepted for publication in Phys.Rev.D

Journal ref: Phys. Rev. D 105, 103027 (2022)

Cosmology with the Laser Interferometer Space Antenna

Authors: Pierre Auclair, David Bacon, Tessa Baker, Tiago Barreiro, Nicola Bartolo, Enis Belgacem, Nicola Bellomo, Ido Ben-Dayan, Daniele Bertacca, Marc Besancon, Jose J. Blanco-Pillado, Diego Blas, Guillaume Boileau, Gianluca Calcagni, Robert Caldwell, Chiara Caprini, Carmelita Carbone, Chia-Feng Chang, Hsin-Yu Chen, Nelson Christensen, Sebastien Clesse, Denis Comelli, Giuseppe Congedo, Carlo Contaldi, Marco Crisostomi , et al. (155 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations exten… ▽ More The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Report number: LISA CosWG-22-03

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Authors: Katja Stock, Dimitri Veras, Maxwell X. Cai, Rainer Spurzem, Simon Portegies Zwart

Abstract: Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of th… ▽ More Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations. △ Less

Submitted 2 March, 2022; originally announced March 2022.

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

Survey of Orion Disks with ALMA (SODA) I: Cloud-level demographics of 873 protoplanetary disks

Authors: S. E. van Terwisga, A. Hacar, E. F. van Dishoeck, R. Oonk, S. Portegies Zwart

Abstract: (Abridged) Surveys of protoplanetary disks in nearby star-forming regions (SFRs) have provided important information on their demographics, but due to their sample sizes, they cannot be used to study how disk properties vary with the environment. We conduct a survey of the unresolved millimeter continuum emission of 873 protoplanetary disks identified by Spitzer in the L1641 and L1647 regions of t… ▽ More (Abridged) Surveys of protoplanetary disks in nearby star-forming regions (SFRs) have provided important information on their demographics, but due to their sample sizes, they cannot be used to study how disk properties vary with the environment. We conduct a survey of the unresolved millimeter continuum emission of 873 protoplanetary disks identified by Spitzer in the L1641 and L1647 regions of the Orion A cloud. This is the largest such survey yet, allowing us to identify even weak trends in the median disk mass as a function of position in the cloud and cluster membership. The sample detection rates and median masses are also compared to those of nearby (<300 pc) SFRs. The sample was observed with ALMA at 225 GHz, with a median rms of 0.08 mJy/beam (1.5 M$_{\oplus}$). The data were reduced and imaged using an innovative parallel data processing approach. We detect 58% (502/873) of the observed disks. This includes 20 disks with dust masses >100 M$_{\oplus}$, and two objects associated with extended dust emission. We infer a median disk dust mass in the full sample of $2.2^{+0.2}_{-0.2}$ M$_{\oplus}$. In L1641 and L1647 median dust masses are $2.1^{+0.2}_{-0.2}$ M$_{\oplus}$ and $2.6^{+0.4}_{-0.5}$ M$_{\oplus}$, respectively. The disk mass distribution of the full sample is similar to that of nearby low-mass SFRs at similar ages of 1-3 Myr. We find only weak trends in disk (dust) masses with galactic longitude and between the YSO clusters identified in the sample, with median masses varying by $\lesssim$ 50%. Age differences may explain these median disk mass variations. Apart from this, disk masses are essentially constant at scales of ~100 pc. This also suggests that the majority of disks, even in different SFRs, are formed with similar initial masses and evolve at similar rates, assuming no external irradiation, with disk mass loss rates of $\sim 10^{-8}$ M$_{\odot}$/yr. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted in A&A. 19 pages, 13 figures, 5 tables in the main text; 5 pages, 1 table and 5 additional figures in appendices

Journal ref: A&A 661, A53 (2022)

Impact of bar resonances in the velocity-space distribution of the solar neighbourhood stars in a self-consistent $N$-body Galactic disc simulation

Authors: Tetsuro Asano, Michiko S. Fujii, Junichi Baba, Jeroen Bédorf, Elena Sellentin, Simon Portegies Zwart

Abstract: The velocity-space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent $N$-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity-space distributions in the $N$-… ▽ More The velocity-space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent $N$-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity-space distributions in the $N$-body model and that of the solar neighbourhood, using Kullback-Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity-space distribution. The KLD fluctuates with time, which indicates the velocity-space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity-space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,φ)=(8.2\;\mathrm{kpc}, 30^{\circ})$, where $R$ and $φ$ are the distance from the galactic centre and the angle with respect to the bar's major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity-space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity-space distribution even though the galactic potential is not static. △ Less

Submitted 23 May, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 9 pages, 11 figures. Accepted by MNRAS

Neural Symplectic Integrator with Hamiltonian Inductive Bias for the Gravitational $N$-body Problem

Authors: Maxwell X. Cai, Simon Portegies Zwart, Damian Podareanu

Abstract: The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general analytical solution for $N>2$. Can an $N$-body problem be solved accurately by a neural network (NN)? Can a NN observe long-term conservation of energy and orbital a… ▽ More The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general analytical solution for $N>2$. Can an $N$-body problem be solved accurately by a neural network (NN)? Can a NN observe long-term conservation of energy and orbital angular momentum? Inspired by Wistom & Holman (1991)'s symplectic map, we present a neural $N$-body integrator for splitting the Hamiltonian into a two-body part, solvable analytically, and an interaction part that we approximate with a NN. Our neural symplectic $N$-body code integrates a general three-body system for $10^{5}$ steps without diverting from the ground truth dynamics obtained from a traditional $N$-body integrator. Moreover, it exhibits good inductive bias by successfully predicting the evolution of $N$-body systems that are no part of the training set. △ Less

Submitted 28 November, 2021; originally announced November 2021.

Comments: 7 pages, 2 figures, accepted for publication at the NeurIPS 2021 workshop "Machine Learning and the Physical Sciences"

Exploring the possibility of Peter Pan discs across stellar mass

Authors: Martijn J. C. Wilhelm, Simon Portegies Zwart

Abstract: Recently, several accreting M dwarf stars have been discovered with ages far exceeding the typical protoplanetary disc lifetime. These `Peter Pan discs' can be explained as primordial discs that evolve in a low-radiation environment. The persistently low masses of the host stars raise the question whether primordial discs can survive up to these ages around stars of higher mass. In this work we ex… ▽ More Recently, several accreting M dwarf stars have been discovered with ages far exceeding the typical protoplanetary disc lifetime. These `Peter Pan discs' can be explained as primordial discs that evolve in a low-radiation environment. The persistently low masses of the host stars raise the question whether primordial discs can survive up to these ages around stars of higher mass. In this work we explore the way in which different mass loss processes in protoplanetary discs limit their maximum lifetimes, and how this depends on host star mass. We find that stars with masses $\lesssim$ 0.6 M$_\odot$ can retain primordial discs for $\sim$50 Myr. At stellar masses $\gtrsim$ 0.8 M$_\odot$, the maximum disc lifetime decreases strongly to below 50 Myr due to relatively more efficient accretion and photoevaporation by the host star. Lifetimes up to 15 Myr are still possible for all host star masses up to $\sim$2 M$_\odot$. For host star masses between 0.6 and 0.8 M$_\odot$, accretion ceases and an inner gap forms before 50 Myr in our models. Observations suggest that such a configuration is rapidly dispersed. We conclude that Peter Pan discs can only occur around M dwarf stars. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Comments: Accepted by MNRAS, 8 pages, 5 figures

Stellar triples on the edge; Comprehensive overview of the evolution of destabilised triples leading to stellar and binary exotica

Authors: S. Toonen, T. C. N. Boekholt, S. Portegies Zwart

Abstract: Hierarchical triple stars are ideal laboratories for studying the interplay between orbital dynamics and stellar evolution. Both stellar wind mass loss and three-body dynamics cooperate to destabilise triples, which can lead to a variety of astrophysical exotica. So far our understanding of their evolution was mainly built upon results from extensive binary-single scattering experiments. Starting… ▽ More Hierarchical triple stars are ideal laboratories for studying the interplay between orbital dynamics and stellar evolution. Both stellar wind mass loss and three-body dynamics cooperate to destabilise triples, which can lead to a variety of astrophysical exotica. So far our understanding of their evolution was mainly built upon results from extensive binary-single scattering experiments. Starting from generic initial conditions, we evolve an extensive set of hierarchical triples using a combination of the triple evolution code TRES and an N-body code. We find that the majority of triples preserve their hierarchy throughout their evolution, which is in contradiction with the commonly adopted picture that unstable triples always experience a chaotic, democratic resonant interaction. The duration of the unstable phase is much longer than expected, so that stellar evolution cannot be neglected. Typically an unstable triple dissolve into a single star and a binary; sometimes democratically (the initial hierarchy is lost and the lightest body usually escapes), but also in a hierarchical way (the tertiary is ejected in a slingshot, independent of its mass). Collisions are common, and mostly involve the two original inner binary components still on the main-sequence. This contradicts the idea that collisions with a giant during democratic encounters dominate. Together with collisions in stable triples, we find that triple evolution is the dominant mechanism for stellar collisions in the Milky Way. Furthermore, our simulations produce runaway and walk-away stars with speeds up to several tens km/s, with a maximum of a few 100km/s. We suggest that destabilised triples can alleviate the tension behind the origin of the observed run-away stars. Lastly, we present a promising indicator to make general predictions for the fate of a specific triple, based on the initial inclination of the system. △ Less

Submitted 1 March, 2022; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 25 pages, 15 figures, 3 tables, accepted by A&A. Reproduction package at DOI: 10.5281/zenodo.6319686

Interstellar objects follow the collapse of molecular clouds

Authors: Susanne Pfalzner, Dylan Paterson, Michele T. Bannister, Simon Portegies Zwart

Abstract: Interstellar objects (ISOs), the parent population of 1I/Oumuamua and 2I/Borisov, are abundant in the interstellar medium of the Milky Way. This means that the interstellar medium, including molecular cloud regions, has three components: gas, dust, and ISOs. From the observational constraints for the field density of ISOs drifting in the solar neighbourhood, we infer a typical molecular cloud of 1… ▽ More Interstellar objects (ISOs), the parent population of 1I/Oumuamua and 2I/Borisov, are abundant in the interstellar medium of the Milky Way. This means that the interstellar medium, including molecular cloud regions, has three components: gas, dust, and ISOs. From the observational constraints for the field density of ISOs drifting in the solar neighbourhood, we infer a typical molecular cloud of 10 pc diameter contains some 10$^{18}$ ISOs. At typical sizes ranging from hundreds of metres to tens of km, ISOs are entirely decoupled from the gas dynamics in these molecular clouds. Here we address the question of whether ISOs can follow the collapse of molecular clouds. We perform low-resolution simulations of the collapse of molecular clouds containing initially static ISO populations toward the point where stars form. In this proof-of-principle study, we find that the interstellar objects definitely follow the collapse of the gas -- and many become bound to the new-forming numerical approximations to future stars (sinks). At minimum, 40\% of all sinks have one or more ISO test particles gravitationally bound to them for the initial ISO distributions tested here. This value corresponds to at least $10^{10}$ actual interstellar objects being bound after three initial free-fall times. Thus, ISOs are a relevant component of star formation. We find that more massive sinks bind disproportionately large fractions of the initial ISO population, implying competitive capture of ISOs. Sinks can also be solitary, as their ISOs can become unbound again -- particularly if sinks are ejected from the system. Emerging planetary systems will thus develop in remarkably varied environments, ranging from solitary to richly populated with bound ISOs. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Comments: 15 pages, 7 figures, accepted for ApJ

Oort cloud Ecology II: The chronology of the formation of the Oort cloud

Authors: Simon Portegies Zwart, Santiago Torres, Maxwell X. Cai, Anthony Brown

Abstract: We present a chronology of the formation and early evolution of the Oort cloud by simulations. These simulations start with the Solar System being born with planets and asteroids in a stellar cluster orbiting the Galactic center. Upon ejection from its birth environment, we continue to follow the evolution of the Solar System while it navigates the Galaxy as an isolated planetary system. We conclu… ▽ More We present a chronology of the formation and early evolution of the Oort cloud by simulations. These simulations start with the Solar System being born with planets and asteroids in a stellar cluster orbiting the Galactic center. Upon ejection from its birth environment, we continue to follow the evolution of the Solar System while it navigates the Galaxy as an isolated planetary system. We conclude that the range in semi-major axis between 100au and several 10$^3$\,au still bears the signatures of the Sun being born in a 1000MSun/pc$^3$ star cluster, and that most of the outer Oort cloud formed after the Solar System was ejected. The ejection of the Solar System, we argue, happened between 20Myr and 50Myr after its birth. Trailing and leading trails of asteroids and comets along the Sun's orbit in the Galactic potential are the by-product of the formation of the Oort cloud. These arms are composed of material that became unbound from the Solar System when the Oort cloud formed. Today, the bulk of the material in the Oort cloud ($\sim 70$\%) originates from the region in the circumstellar disk that was located between $\sim 15$\,au and $\sim 35$\,au, near the current location of the ice giants and the Centaur family of asteroids. According to our simulations, this population is eradicated if the ice-giant planets are born in orbital resonance. Planet migration or chaotic orbital reorganization occurring while the Solar System is still a cluster member is, according to our model, inconsistent with the presence of the Oort cloud. About half the inner Oort cloud, between 100 and $10^4$\,au, and a quarter of the material in the outer Oort cloud, $\apgt 10^4$\,au, could be non-native to the Solar System but was captured from free-floating debris in the cluster or from the circumstellar disk of other stars in the birth cluster. △ Less

Submitted 11 June, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: Corrected proofs for publication in A&A

Journal ref: A&A 652, A144 (2021)

Inside-Out Planet Formation: VI. Oligarchic Coagulation of Planetesimals from a Pebble Ring?

Authors: Maxwell X. Cai, Jonathan C. Tan, Simon Portegies Zwart

Abstract: Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations… ▽ More Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations, we study the collisional evolution of planetesimals formed from such a pebble ring, in particular examining whether a single dominant planet emerges. We consider a variety of models, including some in which the planetesimals are continuing to grow via pebble accretion. We find that the planetesimal ring undergoes oligarchic evolution, and typically turns into 2 or 3 surviving oligarchs on nearly coplanar and circular orbits, independent of the explored initial conditions or form of pebble accretion. The most massive oligarchs typically consist of about $70\%$ of the total mass, with the building-up process typically finishing within $\sim 10^5$ years. However, a relatively massive secondary planet always remains with $\sim30-65\%$ of the mass of the primary. Such secondary planets have properties that are inconsistent with the observed properties of the innermost pairs of planets in STIPs. Thus, for IOPF to be a viable theory for STIP formation, it needs to be shown how oligarchic growth of a relatively massive secondary from the initial pebble ring can be avoided. We discuss some potential additional physical processes that should be included in the modeling and explored as next steps. △ Less

Submitted 13 December, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

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

Lucky planets: how circum-binary planets survive the supernova in one of the inner-binary components

Authors: Fedde Fagginger Auer, Simon Portegies Zwart

Abstract: A planet hardly ever survives the supernova of the host star in a bound orbit, because mass loss in the supernova and the natal kick imparted to the newly formed compact object cause the planet to be ejected. A planet in orbit around a binary has a considerably higher probability to survive the supernova explosion of one of the inner binary stars. In those cases, the planet most likely remains bou… ▽ More A planet hardly ever survives the supernova of the host star in a bound orbit, because mass loss in the supernova and the natal kick imparted to the newly formed compact object cause the planet to be ejected. A planet in orbit around a binary has a considerably higher probability to survive the supernova explosion of one of the inner binary stars. In those cases, the planet most likely remains bound to the companion of the exploding star, whereas the compact object is ejected. We estimate this to happen to $\sim 1/33$ the circum-binary planetary systems. These planetary orbits tend to be highly eccentric ($e \apgt 0.9$), and $\sim 20$\,\% of these planets have retrograde orbits compared to their former binary. The probability that the planet as well as the binary (now with a compact object) remains bound is about ten times smaller ($\sim 3\cdot 10^{-3}$). We then expect the Milky way Galaxy to host $\aplt 10$ x-ray binaries that are still orbited by a planet, and $\aplt 150$ planets that survived in orbit around the compact object's companion. These numbers should be convolved with the fraction of massive binaries that is orbited by a planet. △ Less

Submitted 4 April, 2022; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: published in SciPost Astronomy, see https://scipost.org/SciPostAstro.2.1.002

Journal ref: SciPost Astro. 2, 002 (2022)

Evolution of circumstellar discs in young star-forming regions

Authors: Francisca Concha-Ramírez, Martijn J. C Wilhelm, Simon Portegies Zwart

Abstract: The evolution of circumstellar discs is influenced by their surroundings. The relevant processes include external photoevaporation due to nearby stars, and dynamical truncations. The impact of these processes on disc populations depends on the star-formation history and on the dynamical evolution of the region. Since star formation history and the phase-space characteristics of the stars are impor… ▽ More The evolution of circumstellar discs is influenced by their surroundings. The relevant processes include external photoevaporation due to nearby stars, and dynamical truncations. The impact of these processes on disc populations depends on the star-formation history and on the dynamical evolution of the region. Since star formation history and the phase-space characteristics of the stars are important for the evolution of the discs, we start simulating the evolution of the star cluster with the results of molecular cloud collapse simulations. In the simulation we form stars with circumstellar discs, which can be affected by different processes. Our models account for the viscous evolution of the discs, internal and external photoevaporation of gas, external photoevaporation of dust, and dynamical truncations. All these processes are resolved together with the dynamical evolution of the cluster, and the evolution of the stars. An extended period of star formation, lasting for at least 2 Myr, results in some discs being formed late. These late formed discs have a better chance of survival because the cluster gradually expands with time, and a lower local stellar density reduces the effects of photoevaporation and dynamical truncation. Late formed discs can then be present in regions of high UV radiation, solving the proplyd lifetime problem. We also find a considerable fraction of discs that lose their gas content, but remain sufficiently rich in solids to be able to form a rocky planetary system. △ Less

Submitted 26 May, 2022; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS. 15 pages, 12 figures

Pulsating ULXs: large pulsed fraction excludes strong beaming

Authors: Alexander A. Mushtukov, Simon Portegies Zwart, Sergey S. Tsygankov, Dmitrij I. Nagirner, Juri Poutanen

Abstract: The recent discovery of pulsating ultra-luminous X-ray sources (ULXs) shows that the apparent luminosity of accreting neutron stars (NSs) can exceed the Eddington luminosity by a factor of hundreds. The relation between the actual and apparent luminosity is a key ingredient in theoretical models of ULXs but it is still under debate. A typical feature of the discovered pulsating ULXs is a large pul… ▽ More The recent discovery of pulsating ultra-luminous X-ray sources (ULXs) shows that the apparent luminosity of accreting neutron stars (NSs) can exceed the Eddington luminosity by a factor of hundreds. The relation between the actual and apparent luminosity is a key ingredient in theoretical models of ULXs but it is still under debate. A typical feature of the discovered pulsating ULXs is a large pulsed fraction (PF). Using Monte Carlo simulations, we consider a simple geometry of accretion flow and test the possibility of simultaneous presence of a large luminosity amplification due the geometrical beaming and a high PF. We argue that these factors largely exclude each other and only a negligible fraction of strongly beamed ULX pulsars can show PF above 10 per cent. Discrepancy between this conclusion and current observations indicate that pulsating ULXs are not strongly beamed and their apparent luminosity is close to the actual one. △ Less

Submitted 9 December, 2020; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 6 pages, 5 figures; accepted for publication in MNRAS

Oort cloud Ecology II: Extra-solar Oort clouds and the origin of asteroidal interlopers

Authors: Simon Portegies Zwart

Abstract: We simulate the formation and evolution of Oort clouds around the 200 nearest stars (within 16pc according to the Gaia DR2) database. This study is performed by numerically integrating the planets and minor bodies in orbit around the parent star and in the Galactic potential. The calculations start 1\,Gyr ago and continue for 100Myr into the future. In this time frame, we simulate how asteroids (a… ▽ More We simulate the formation and evolution of Oort clouds around the 200 nearest stars (within 16pc according to the Gaia DR2) database. This study is performed by numerically integrating the planets and minor bodies in orbit around the parent star and in the Galactic potential. The calculations start 1\,Gyr ago and continue for 100Myr into the future. In this time frame, we simulate how asteroids (and planets) are ejected from the star's vicinity and settle in an Oort cloud and how they escape the local stellar gravity to form tidal steams. A fraction of 0.0098 to 0.026 of the asteroids remain bound to their parent star. The orbits of these asteroids isotropizes and circularizes due to the influence of the Galactic tidal field to eventually form an Oort cloud between 10^4 and 2 10^5au. We estimate that 6% of the nearby stars may have a planet in its Oort cloud. The majority of asteroids (and some of the planets) become unbound from the parent star to become free floating in the Galactic potential. These soli lapides remain in a similar orbit around the Galactic center as their host star, forming dense streams of rogue interstellar asteroids and planets. The Solar system occasionally passes through such tidal streams, potentially giving rise to occasional close encounters with object in this stream. The two recently discovered objects, 1I/(2017 Q3) 'Oumuamua and 2I/(2019 Q4) Borisov, may be such objects. Although the direction from which an individual solus lapis originated cannot easily be traced back to the original host, multiple such objects coming from the same source might help to identify their origin. At the moment the Solar system is in the bow or wake of the tidal stream of 10 of the nearby stars which might contribute considerably to the interaction rate. (abridged) △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: accepted for publication in A&A, see animation https://youtu.be/0fYeAW3e9bQ

Journal ref: A&A 647, A136 (2021)

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

Authors: Claude Cournoyer-Cloutier, Aaron Tran, Sean Lewis, Joshua E. Wall, William E. Harris, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Simon Portegies Zwart, Alison Sills

Abstract: The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled M… ▽ More The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 15 pages, 14 figures, submitted to MNRAS and edited to address positive referee's report

APEX-SEPIA660 Early Science: Gas at densities above $10^7$ cm$^{-3}$ towards OMC-1

Authors: A. Hacar, M. R. Hogerheijde, D. Harsono, S. Portegies Zwart, C. De Breuck, K. Torstensson, W. Boland, A. M. Baryshev, R. Hesper, J. Barkhof, J. Adema, M. E. Bekema, A. Koops, A. Khudchenko, R. Stark

Abstract: Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1… ▽ More Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1 region located in the vicinity of the Orion Nebula Cluster (ONC). Methods. We mapped the molecular emission at 652 GHz in OMC-1 as part of the APEX-SEPIA660 Early Science. Results. We detect bright and extended N$_2$H$^+$ (J=7-6) line emission along the entire OMC-1 region. Comparisons with previous ALMA data of the (J=1-0) transition and radiative transfer models indicate that the line intensities observed in this N$_2$H$^+$ (7-6) line are produced by large mass reservoirs of gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Conclusions. The first detection of this N$_2$H$^+$ (7-6) line at parsec-scales demonstrates the extreme density conditions of the star-forming gas in young massive clusters such as the ONC. Our results highlight the unique combination of sensitivity and mapping capabilities of the new SEPIA660 receiver for the study of the ISM properties at high frequencies. △ Less

Submitted 26 September, 2020; originally announced September 2020.

Comments: 8 pages, 5 Figures. Accepted for publication by A&A

Journal ref: A&A 644, A133 (2020)

The Ecological Impact of High-performance Computing in Astrophysics

Authors: Simon Portegies Zwart

Abstract: The importance of computing in astronomy continues to increase, and so is its impact on the environment. When analyzing data or performing simulations, most researchers raise concerns about the time to reach a solution rather than its impact on the environment. Luckily, a reduced time-to-solution due to faster hardware or optimizations in the software generally also leads to a smaller carbon footp… ▽ More The importance of computing in astronomy continues to increase, and so is its impact on the environment. When analyzing data or performing simulations, most researchers raise concerns about the time to reach a solution rather than its impact on the environment. Luckily, a reduced time-to-solution due to faster hardware or optimizations in the software generally also leads to a smaller carbon footprint. This is not the case when the reduced wall-clock time is achieved by overclocking the processor, or when using supercomputers. The increase in the popularity of interpreted scripting languages, and the general availability of high-performance workstations form a considerable threat to the environment. A similar concern can be raised about the trend of running single-core instead of adopting efficient many-core programming paradigms. In astronomy, computing is among the top producers of green-house gasses, surpassing telescope operations. Here I hope to raise the awareness of the environmental impact of running non-optimized code on overpowered computer hardware. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Originated at EAS 2020 conference, sustainability session by https://astronomersforplanet.earth - published in Nature Astronomy, September 2020

On the survival of resonant and non-resonant planetary systems in star clusters

Authors: Katja Stock, Maxwell X. Cai, Rainer Spurzem, M. B. N. Kouwenhoven, Simon Portegies Zwart

Abstract: Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing… ▽ More Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 21 pages, 21 figures, accepted by MNRAS

Spectrum formation in X-ray pulsars at very low mass accretion rate: Monte-Carlo approach

Authors: Alexander A. Mushtukov, Valery F. Suleimanov, Sergey S. Tsygankov, Simon Portegies Zwart

Abstract: It has been recently discovered that the transition of X-ray pulsars to the low luminosity state ($L\lesssim10^{35}\,{\rm erg\,s^{-1}}$) is accompanied by a dramatic spectral change. Namely, the typical power-law-like spectrum with high energy cutoff transforms into a two-component structure with a possible cyclotron absorption feature on top of it. It was proposed that these spectral characterist… ▽ More It has been recently discovered that the transition of X-ray pulsars to the low luminosity state ($L\lesssim10^{35}\,{\rm erg\,s^{-1}}$) is accompanied by a dramatic spectral change. Namely, the typical power-law-like spectrum with high energy cutoff transforms into a two-component structure with a possible cyclotron absorption feature on top of it. It was proposed that these spectral characteristics can be explained qualitatively by the emission of cyclotron photons in the atmosphere of the neutron star caused by collisional excitation of electrons to upper Landau levels and further comptonization of the photons by electron gas. The latter is expected to be overheated in a thin top layer of the atmosphere. In this paper, we perform Monte Carlo simulations of the radiative transfer in the atmosphere of an accreting neutron star while accounting for a resonant scattering of polarized X-ray photons by thermally distributed electrons. The spectral shape is shown to be strongly polarization-dependent in soft X-rays ($\lesssim 10\,{\rm keV}$) and near the cyclotron scattering feature. The results of our numerical simulations are tested against the observational data of X-ray pulsar A 0535+262 in the low luminosity state. We show that the spectral shape of the pulsar can be reproduced by the proposed theoretical model. The applications of the discovery to the observational studies of accreting neutron stars are discussed. △ Less

Submitted 17 March, 2021; v1 submitted 24 June, 2020; originally announced June 2020.

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

Effects of stellar density on the photoevaporation of circumstellar discs

Authors: Francisca Concha-Ramírez, Martijn J. C. Wilhelm, Simon Portegies Zwart, Sierk E. van Terwisga, Alvaro Hacar

Abstract: Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the… ▽ More Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the most important ones. Theoretical and observational evidence point to circumstellar discs losing mass quickly when in the vicinity of massive, bright stars. In this work we simulate circumstellar discs in clustered environments in a range of stellar densities, where the photoevaporation mass-loss process is resolved simultaneously with the stellar dynamics, stellar evolution, and the viscous evolution of the discs. Our results indicate that external photoevaporation is efficient in depleting disc masses and that the degree of its effect is related to stellar density. We find that a local stellar density lower than 100 stars pc$^{-2}$ is necessary for discs massive enough to form planets to survive for \SI{2.0}{Myr}. There is an order of magnitude difference in the disc masses in regions of projected density 100 stars pc$^{-2}$ versus $10^4$ stars pc$^{-2}$. We compare our results to observations of the Lupus clouds, the Orion Nebula Cluster, the Orion Molecular Cloud-2, Taurus, and NGC 2024, and find that the trends observed between region density and disc masses are similar to those in our simulations. △ Less

Submitted 20 November, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS. 9 pages, 4 figures

Trimodal structure of Hercules stream explained by originating from bar resonances

Authors: Tetsuro Asano, Michiko S. Fujii, Junichi Baba, Jeroen Bédorf, Elena Sellentin, Simon Portegies Zwart

Abstract: Gaia Data Release 2 revealed detailed structures of nearby stars in phase space. These include the Hercules stream, whose origin is still debated. Most of the previous numerical studies conjectured that the observed structures originate from orbits in resonance with the bar, based on static potential models for the Milky Way. We, in contrast, approach the problem via a self-consistent, dynamic, an… ▽ More Gaia Data Release 2 revealed detailed structures of nearby stars in phase space. These include the Hercules stream, whose origin is still debated. Most of the previous numerical studies conjectured that the observed structures originate from orbits in resonance with the bar, based on static potential models for the Milky Way. We, in contrast, approach the problem via a self-consistent, dynamic, and morphologically well-resolved model, namely a full $N$-body simulation of the Milky Way. Our simulation comprises about 5.1 billion particles in the galactic stellar bulge, bar, disk, and dark-matter halo and is evolved to 10 Gyr. Our model's disk component is composed of 200 million particles, and its simulation snapshots are stored every 10 Myr, enabling us to resolve and classify resonant orbits of representative samples of stars. After choosing the Sun's position in the simulation, we compare the distribution of stars in its neighborhood with Gaia's astrometric data, thereby establishing the role of identified resonantly trapped stars in the formation of Hercules-like structures. From our orbital spectral-analysis we identify multiple, especially higher order resonances. Our results suggest that the Hercules stream is dominated by the 4:1 and 5:1 outer Lindblad and corotation resonances. In total, this yields a trimodal structure of the Hercules stream. From the relation between resonances and ridges in phase space, our model favored a slow pattern speed of the Milky-Way bar (40--45 $\mathrm{km \; s^{-1} \; kpc^{-1}}$). △ Less

Submitted 15 September, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 11 pages, 9 figures, MNRAS accepted

The evolution of stellar triples: The most common evolutionary pathways

Authors: S. Toonen, S. Portegies Zwart, A. Hamers, D. Bandopadhyay

Abstract: Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples and higher-order multiples are common. Whereas the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. To fill this gap in our general understanding of stellar lives, we aim to systematically explo… ▽ More Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples and higher-order multiples are common. Whereas the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. To fill this gap in our general understanding of stellar lives, we aim to systematically explore the long-term evolution of triples and to map out the most common evolutionary pathways that triples go through. We quantitatively study how triples evolve, which processes are most relevant, and how this differs from binary evolution. We simulate the evolution of several large populations of triples with a population synthesis approach. We make use of the triple evolution code TRES to simulate the evolution of each triple in a consistent way; including three-body dynamics (based on the secular approach), stellar evolution and their mutual influences. We simulate the evolution of the system up until mass transfer starts, the system becomes dynamically unstable, or a Hubble time has passed. We find that stellar interactions are common in triples. Compared to a binary population, we find that the fraction of systems that can undergo mass transfer is about 2 to 3 times larger in triples. Moreover, whereas in binaries the orbits typically reach circularisation before Roche-lobe overflow, this is not true anymore in triples. In our simulations, about 40% of systems retain an eccentric orbit. Additionally, we discuss various channels of triple evolution in detail such as those where the secondary or the tertiary is the first star to initiate a mass transfer event. △ Less

Submitted 6 July, 2020; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: updated version: accepted for publication in A&A 18 pages, 16 figures, 2 tables

Modelling of the Effects of Stellar Feedback during Star Cluster Formation Using a Hybrid Gas and N-Body Method

Authors: Joshua E. Wall, Mordecai-Mark Mac Low, Stephen L. W. McMillan, Ralf S. Klessen, Simon Portegies Zwart, Andrew Pellegrino

Abstract: Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of… ▽ More Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ~1000 stars formed from clouds with a mass range of 10^3-10^5 M_sun. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations. △ Less

Submitted 10 October, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted to ApJ. 29 pages, 18 figures. Source code at https://bitbucket.org/torch-sf/torch/ and documentation at https://torch-sf.bitbucket.io/

Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics

Authors: Simon Portegies Zwart, Inti Pelupessy, Carmen Martinez-Barbosa, Arjen van Elteren, Steve McMillan

Abstract: Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method… ▽ More Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset. The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. ... abridged ... △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: Accepted for publication in Communications in Nonlinear Science and Numerical Simulation (CNSNS) The associated software is part of the AMUSE framework and can be downloaded from http:www.amusecode.org

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.

Comments: Accepted for publication in MNRAS

Gargantuan chaotic gravitational three-body systems and their irreversibility to the Planck length

Authors: Tjarda Boekholt, Simon Portegies Zwart, Mauri Valtonen

Abstract: Chaos is present in most stellar dynamical systems and manifests itself through the exponential growth of small perturbations. Exponential divergence drives time irreversibility and increases the entropy in the system. A numerical consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation betwee… ▽ More Chaos is present in most stellar dynamical systems and manifests itself through the exponential growth of small perturbations. Exponential divergence drives time irreversibility and increases the entropy in the system. A numerical consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation between chaos in stellar dynamical systems and the level of irreversibility remained undetermined. In this work we study chaotic three-body systems in free fall initially using the accurate and precise N-body code Brutus, which goes beyond standard double-precision arithmetic. We demonstrate that the fraction of irreversible solutions decreases as a power law with numerical accuracy. This can be derived from the distribution of amplification factors of small initial perturbations. Applying this result to systems consisting of three massive black holes with zero total angular momentum, we conclude that up to five percent of such triples would require an accuracy of smaller than the Planck length in order to produce a time-reversible solution, thus rendering them fundamentally unpredictable. △ Less

Submitted 10 February, 2020; originally announced February 2020.

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

Identifying Earth-impacting asteroids using an artificial neural network

Authors: John D. Hefele, Francesco Bortolussi, Simon Portegies Zwart

Abstract: By means of a fully connected artificial neural network, we identified asteroids with the potential to impact Earth. The resulting instrument, named the Hazardous Object Identifier (HOI), was trained on the basis of an artificial set of known impactors which were generated by launching objects from Earth's surface and integrating them backward in time. HOI was able to identify 95.25% of the known… ▽ More By means of a fully connected artificial neural network, we identified asteroids with the potential to impact Earth. The resulting instrument, named the Hazardous Object Identifier (HOI), was trained on the basis of an artificial set of known impactors which were generated by launching objects from Earth's surface and integrating them backward in time. HOI was able to identify 95.25% of the known impactors simulated that were present in the test set as potential impactors. In addition, HOI was able to identify 90.99% of the potentially hazardous objects identified by NASA, without being trained on them directly. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Journal ref: A&A 634, A45 (2020)

Newton vs the machine: solving the chaotic three-body problem using deep neural networks

Authors: Philip G. Breen, Christopher N. Foley, Tjarda Boekholt, Simon Portegies Zwart

Abstract: Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. W… ▽ More Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. We show that an ensemble of solutions obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at fixed computational cost and up to 100 million times faster than a state-of-the-art solver. Our results provide evidence that, for computationally challenging regions of phase-space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black-hole binary systems or the origin of the core collapse in dense star clusters. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 6 pages, 6 figures

Searching for Solar Siblings in APOGEE and $Gaia$ DR2 with N-body Simulations

Authors: Jeremy J. Webb, Natalie Price-Jones, Jo Bovy, Simon Portegies Zwart, Jason A. S. Hunt, J. Ted Mackereth, Henry W. Leung

Abstract: We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of si… ▽ More We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67. △ Less

Submitted 23 March, 2020; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 12 pages, 9 figures, 3 tables, Accepted for publication in MNRAS

Bonsai-SPH: A GPU accelerated astrophysical Smoothed Particle Hydrodynamics code

Authors: Jeroen Bédorf, Simon Portegies Zwart

Abstract: We present the smoothed-particle hydrodynamics simulation code, Bonsai-SPH, which is a continuation of our previously developed gravity-only hierarchical $N$-body code (called Bonsai). The code is optimized for Graphics Processing Unit (GPU) accelerators which enables researchers to take advantage of these powerful computational resources. Bonsa-SPH produces simulation results comparable with stat… ▽ More We present the smoothed-particle hydrodynamics simulation code, Bonsai-SPH, which is a continuation of our previously developed gravity-only hierarchical $N$-body code (called Bonsai). The code is optimized for Graphics Processing Unit (GPU) accelerators which enables researchers to take advantage of these powerful computational resources. Bonsa-SPH produces simulation results comparable with state-of-the-art, CPU based, codes, but using an order of magnitude less computation time. The code is freely available online and the details are described in this work. △ Less

Submitted 14 February, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: Updated intro and multi-GPU sections. 23 pages, 9 figures. Submission to SciPost