-
Unsupervised anomaly detection in spatio-temporal stream network sensor data
Authors:
Edgar Santos-Fernandez,
Jay M. Ver Hoef,
Erin E. Peterson,
James McGree,
Cesar A. Villa,
Catherine Leigh,
Ryan Turner,
Cameron Roberts,
Kerrie Mengersen
Abstract:
The use of in-situ digital sensors for water quality monitoring is becoming increasingly common worldwide. While these sensors provide near real-time data for science, the data are prone to technical anomalies that can undermine the trustworthiness of the data and the accuracy of statistical inferences, particularly in spatial and temporal analyses. Here we propose a framework for detecting anomal…
▽ More
The use of in-situ digital sensors for water quality monitoring is becoming increasingly common worldwide. While these sensors provide near real-time data for science, the data are prone to technical anomalies that can undermine the trustworthiness of the data and the accuracy of statistical inferences, particularly in spatial and temporal analyses. Here we propose a framework for detecting anomalies in sensor data recorded in stream networks, which takes advantage of spatial and temporal autocorrelation to improve detection rates. The proposed framework involves the implementation of effective data imputation to handle missing data, alignment of time-series to address temporal disparities, and the identification of water quality events. We explore the effectiveness of a suite of state-of-the-art statistical methods including posterior predictive distributions, finite mixtures, and Hidden Markov Models (HMM). We showcase the practical implementation of automated anomaly detection in near-real time by employing a Bayesian recursive approach. This demonstration is conducted through a comprehensive simulation study and a practical application to a substantive case study situated in the Herbert River, located in Queensland, Australia, which flows into the Great Barrier Reef. We found that methods such as posterior predictive distributions and HMM produce the best performance in detecting multiple types of anomalies. Utilizing data from multiple sensors deployed relatively near one another enhances the ability to distinguish between water quality events and technical anomalies, thereby significantly improving the accuracy of anomaly detection. Thus, uncertainty and biases in water quality reporting, interpretation, and modelling are reduced, and the effectiveness of subsequent management actions improved.
△ Less
Submitted 11 September, 2024;
originally announced September 2024.
-
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.
-
New candidate hypervelocity red clump stars in the inner Galactic bulge
Authors:
A. Luna,
T. Marchetti,
M. Rejkuba,
N. W. C. Leigh,
J. Alonso-García,
A. Valenzuela Navarro,
D. Minniti,
L. C. Smith
Abstract:
We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the de…
▽ More
We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the detection of HVSs, ejection mechanism models can be constrained by exploring the stellar dynamics in the Galactic centre through a recent stellar interaction with the SMBH. Based on a previously developed methodology by our group, we searched with a sample of preliminary data from version 2 of the Vista Variables in the Via Lactea (VVV) Infrared Astrometric Catalogue (VIRAC2) and Gaia DR3 data, including accurate optical and NIR proper motions. This search resulted in a sample of 46 stars with transverse velocities larger than the local escape velocity within the Galactic bulge, of which 4 are prime candidate HVSs with high-proper motions consistent with being ejections from the Galactic centre. Adding to that, we studied a sample of reddened stars without a Gaia DR3 counterpart and found 481 stars with transverse velocities larger than the local escape velocity, from which 65 stars have proper motions pointing out of the Galactic centre and are candidate HVSs. In total, we found 69 candidate HVSs pointing away from the Galactic centre with transverse velocities larger than the local escape velocity.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
A catalogue of Galactic GEMS: Globular cluster Extra-tidal Mock Stars
Authors:
Steffani M. Grondin,
Jeremy J. Webb,
James M. M. Lane,
Joshua S. Speagle,
Nathan W. C. Leigh
Abstract:
This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are p…
▽ More
This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way's dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way.
△ Less
Submitted 16 January, 2024; v1 submitted 13 October, 2023;
originally announced October 2023.
-
In-plane Tidal Disruption of Stars in Disks of Active Galactic Nuclei
Authors:
Taeho Ryu,
Barry McKernan,
Saavik Ford,
Matteo Cantiello,
Matthew Graham,
Daniel Stern,
Nathan W. C Leigh
Abstract:
Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk g…
▽ More
Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk gas (pro- and retro-grade). Key findings are: 1) Debris experiences continuous perturbations from the disk gas, which can result in significant and continuous changes in debris energy and angular momentum compared to `naked' TDEs. 2) Above a critical density of a disk around a SMBH with mass $M_{\bullet}$ ($ρ_{\rm crit} \sim 10^{-8}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$) for retrograde stars, both bound and unbound debris is fully mixed into the disk. The density threshold for no bound debris return, inhibiting the accretion component of TDEs, is $ρ_{\rm crit,bound} \sim 10^{-9}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$. 3) Observationally, AGN-TDEs transition from resembling naked TDEs in the limit of $ρ_{\rm disk}\lesssim 10^{-2}ρ_{\rm crit,bound}$ to fully muffled TDEs with associated inner disk state changes at $ρ_{\rm disk}\gtrsimρ_{\rm crit,bound}$, with a superposition of AGN+TDE in between. Stellar or remnant passages themselves can significantly perturb the inner disk. This can lead to an immediate X-ray signature and optically detectable inner disk state changes, potentially contributing to the changing-look AGN phenomenon. 4) Debris mixing can enriches the average disk metallicity over time if the star's metallicity exceeds that of the disk gas.
△ Less
Submitted 1 October, 2023;
originally announced October 2023.
-
The dominant mechanism(s) for populating the outskirts of star clusters with neutron star binaries
Authors:
Nathan W. C. Leigh,
Claire S. Ye,
Steffani M. Grondin,
Giacomo Fragione,
Jeremy J. Webb,
Craig O. Heinke
Abstract:
It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of or…
▽ More
It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of orbital separations should a BH-BH binary be both dynamically hard and produce a sufficiently low recoil velocity to retain the NS binary in the cluster. Hence, BH binaries are in general likely to eject NSs from clusters. We explore several alternative mechanisms that would cause NS/MSP binaries to be observed in the outskirts of their host clusters after a Hubble time. The most likely mechanism is a three-body interaction involving the NS/MSP binary and a normal star. We compare to Monte Carlo simulations of cluster evolution for the globular clusters NGC 6752 and 47 Tuc, and show that the models not only confirm that normal three-body interactions involving all stellar-mass objects are the dominant mechanism for putting NS/MSP binaries into the cluster outskirts, they also reproduce the observed NS/MSP binary radial distributions without needing to invoke the presence of a massive BH binary. Higher central densities and an episode of core-collapse can broaden the radial distributions of NSs/MSPs and NS/MSP binaries due to three-body interactions, making these clusters more likely to host NSs in the cluster outskirts.
△ Less
Submitted 22 September, 2023;
originally announced September 2023.
-
The chaotic four-body problem in Newtonian gravity -- II. An ansatz-based approach to analytic solutions
Authors:
Carlos M. Barrera Retamal,
Nathan W. C. Leigh,
Nicholas C. Stone
Abstract:
In this paper, we continue our analysis of the chaotic four-body problem by presenting a general ansatz-based analytic treatment using statistical mechanics, where each outcome of the four-body problem is regarded as some variation of the three-body problem (e.g., when two single stars are produced, called the 2+1+1 outcome, each ejection event is modeled as its own three-body interaction by assum…
▽ More
In this paper, we continue our analysis of the chaotic four-body problem by presenting a general ansatz-based analytic treatment using statistical mechanics, where each outcome of the four-body problem is regarded as some variation of the three-body problem (e.g., when two single stars are produced, called the 2+1+1 outcome, each ejection event is modeled as its own three-body interaction by assuming that the ejections are well separated in time). This is a generalization of the statistical mechanics treatment of the three-body problem based on the density-of-states formalism. In our case, we focus on the interaction of two binary systems, after which we divide our results into three possible outcome scenarios (2+2, 2+1+1, and 3+1). For each outcome, we apply an ansatz-based approach to deriving analytic distribution functions that describe the properties of the products of chaotic four-body interactions involving point particles. To test our theoretical distributions, we perform a set of scattering simulations in the equal-mass point particle limit using FEWBODY. We compare our final theoretical distributions to the simulations for each particular scenario, finding consistently good agreement between the two. The highlights of our results include that binary-binary scatterings act to systematically destroy binaries producing instead a single binary and two ejected stars or a stable triple, the 2+2 outcome produces the widest binaries and the 2+1+1 outcome produces the most compact binaries.
△ Less
Submitted 23 August, 2023;
originally announced August 2023.
-
Conditional normalization in time series analysis
Authors:
Puwasala Gamakumara,
Edgar Santos-Fernandez,
Priyanga Dilini Talagala,
Rob J. Hyndman,
Kerrie Mengersen,
Catherine Leigh
Abstract:
Time series often reflect variation associated with other related variables. Controlling for the effect of these variables is useful when modeling or analysing the time series. We introduce a novel approach to normalize time series data conditional on a set of covariates. We do this by modeling the conditional mean and the conditional variance of the time series with generalized additive models us…
▽ More
Time series often reflect variation associated with other related variables. Controlling for the effect of these variables is useful when modeling or analysing the time series. We introduce a novel approach to normalize time series data conditional on a set of covariates. We do this by modeling the conditional mean and the conditional variance of the time series with generalized additive models using a set of covariates. The conditional mean and variance are then used to normalize the time series. We illustrate the use of conditionally normalized series using two applications involving river network data. First, we show how these normalized time series can be used to impute missing values in the data. Second, we show how the normalized series can be used to estimate the conditional autocorrelation function and conditional cross-correlation functions via additive models. Finally we use the conditional cross-correlations to estimate the time it takes water to flow between two locations in a river network.
△ Less
Submitted 21 May, 2023;
originally announced May 2023.
-
Short-term prediction of stream turbidity using surrogate data and a meta-model approach
Authors:
Bhargav Rele,
Caleb Hogan,
Sevvandi Kandanaarachchi,
Catherine Leigh
Abstract:
Many water-quality monitoring programs aim to measure turbidity to help guide effective management of waterways and catchments, yet distributing turbidity sensors throughout networks is typically cost prohibitive. To this end, we built and compared the ability of dynamic regression (ARIMA), long short-term memory neural nets (LSTM), and generalized additive models (GAM) to forecast stream turbidit…
▽ More
Many water-quality monitoring programs aim to measure turbidity to help guide effective management of waterways and catchments, yet distributing turbidity sensors throughout networks is typically cost prohibitive. To this end, we built and compared the ability of dynamic regression (ARIMA), long short-term memory neural nets (LSTM), and generalized additive models (GAM) to forecast stream turbidity one step ahead, using surrogate data from relatively low-cost in-situ sensors and publicly available databases. We iteratively trialled combinations of four surrogate covariates (rainfall, water level, air temperature and total global solar exposure) selecting a final model for each type that minimised the corrected Akaike Information Criterion. Cross-validation using a rolling time-window indicated that ARIMA, which included the rainfall and water-level covariates only, produced the most accurate predictions, followed closely by GAM, which included all four covariates. We constructed a meta-model, trained on time-series features of turbidity, to take advantage of the strengths of each model over different time points and predict the best model (that with the lowest forecast error one-step prior) for each time step. The meta-model outperformed all other models, indicating that this methodology can yield high accuracy and may be a viable alternative to using measurements sourced directly from turbidity-sensors where costs prohibit their deployment and maintenance, and when predicting turbidity across the short term. Our findings also indicated that temperature and light-associated variables, for example underwater illuminance, may hold promise as cost-effective, high-frequency surrogates of turbidity, especially when combined with other covariates, like rainfall, that are typically measured at coarse levels of spatial resolution.
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
Globular Cluster UVIT legacy Survey (GlobUleS) III. Omega Centauri in Far-Ultraviolet
Authors:
Deepthi S. Prabhu,
Annapurni Subramaniam,
Snehalata Sahu,
Chul Chung,
Nathan W. C. Leigh,
Emanuele Dalessandro,
Sourav Chatterjee,
N. Kameswara Rao,
Michael Shara,
Patrick Cote,
Samyaday Choudhury,
Gajendra Pandey,
Aldo A. R. Valcarce,
Gaurav Singh,
Joesph E. Postma,
Sharmila Rani,
Avrajit Bandyopadhyay,
Aaron M. Geller,
John Hutchings,
Thomas Puzia,
Mirko Simunovic,
Young-Jong Sohn,
Sivarani Thirupathi,
Ramakant Singh Yadav
Abstract:
We present the first comprehensive study of the most massive globular cluster Omega Centauri in the far-ultraviolet (FUV) extending from the center to ~ 28% of the tidal radius using the Ultraviolet Imaging Telescope aboard AstroSat. A comparison of the FUV-optical color-magnitude diagrams with available canonical models reveals that the horizontal branch (HB) stars bluer than the knee (hHBs) and…
▽ More
We present the first comprehensive study of the most massive globular cluster Omega Centauri in the far-ultraviolet (FUV) extending from the center to ~ 28% of the tidal radius using the Ultraviolet Imaging Telescope aboard AstroSat. A comparison of the FUV-optical color-magnitude diagrams with available canonical models reveals that the horizontal branch (HB) stars bluer than the knee (hHBs) and the white dwarfs (WDs) are fainter in the FUV by ~ 0.5 mag than model predictions. They are also fainter than their counterparts in M13, another massive cluster. We simulated HB with at least five subpopulations, including three He-rich populations with a substantial He enrichment of Y up to 0.43 dex, to reproduce the observed FUV distribution. We find the He-rich younger subpopulations to be radially more segregated than the He-normal older ones, suggesting an in-situ enrichment from older generations. The Omega Cen hHBs span the same effective temperature range as their M13 counterparts, but some have smaller radii and lower luminosities. This may suggest that a fraction of Omega Cen hHBs are less massive than those of M13, similar to the result derived from earlier spectroscopic studies of outer extreme HB stars. The WDs in Omega Cen and M13 have similar luminosity-radius-effective temperature parameters, and 0.44 - 0.46 M$_\odot$ He-core WD model tracks evolving from progenitors with Y = 0.4 dex are found to fit the majority of these. This study provides constraints on the formation models of Omega Cen based on the estimated range in age, [Fe/H] and Y (in particular), for the HB stars.
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
Searching for the extra-tidal stars of globular clusters using high-dimensional analysis and a core particle spray code
Authors:
Steffani M. Grondin,
Jeremy J. Webb,
Nathan W. C. Leigh,
Joshua S. Speagle,
Reem J. Khalifeh
Abstract:
Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data ana…
▽ More
Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data analysis and a particle spray code to identify extra-tidal stars of any Galactic globular cluster using M3 as a case study. Using the t-Stochastic Neighbour Embedding (t-SNE) and Uniform Manifold Approximation and Projection (UMAP) machine learning dimensionality reduction algorithms, we first identify a set of 103 extra-tidal candidates in the APOGEE DR17 data catalogue with chemical abundances similar to M3 stars. To confirm each candidate's extra-tidal nature, we introduce corespray - a new Python-based three-body particle spray code that simulates extra-tidal stars for any Galactic globular cluster. Using Gaia EDR3 proper motions and APOGEE DR17 radial velocities, we apply multivariate Gaussian modelling and an extreme deconvolution to identify the extra-tidal candidates that are more likely to be associated with a distribution of corespray-simulated M3 extra-tidal stars than the field. Through these methods, we identify 10 new high-probability extra-tidal stars produced via three-body interactions in M3. We also explore whether any of our extra-tidal candidates are consistent with being ejected from M3 through different dynamical processes. Future applications of corespray will yield better understandings of core dynamics, star formation histories and binary fractions in globular clusters.
△ Less
Submitted 16 November, 2022; v1 submitted 22 July, 2022;
originally announced July 2022.
-
Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc
Authors:
Syeda S. Nasim,
Gaia Fabj,
Freddy Caban,
Amy Secunda,
K. E. Saavik Ford,
Barry McKernan,
Jillian M. Bellovary,
Nathan W. C. Leigh,
Wladimir Lyra
Abstract:
Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population w…
▽ More
Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population which drives sBH-sBH mergers detectable in gravitational waves using LIGO-Virgo-KAGRA (LVK) or sBH-SMBH mergers detectable with LISA (Laser Interferometer Space Antenna). Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive 'immortal' stars. Here, we investigate the behaviour of polar and retrograde orbiters $(i \geq 90^{\circ})$ interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation $(i<90^{\circ})$ during capture, and decrease $a$ dramatically towards the SMBH. For sBH, we find a critical angle $i_{\rm ret} \sim 113^{\circ}$, below which retrograde sBH decay towards embedded prograde orbits $(i \to 0^{\circ})$, while for $i_{\rm o}>i_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution.
△ Less
Submitted 17 May, 2023; v1 submitted 19 July, 2022;
originally announced July 2022.
-
The thermodynamics of stellar multiplicity: dynamical evolution of binary star populations in dense stellar environments
Authors:
N. W. C. Leigh,
N. C. Stone,
J. J. Webb,
W. Lyra
Abstract:
We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant…
▽ More
We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single-binary interactions, that drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step toward the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions.
△ Less
Submitted 30 May, 2022;
originally announced May 2022.
-
Tidal Disruption Events by Compact Supermassive Black Hole Binaries
Authors:
Taeho Ryu,
Alessandro A. Trani,
Nathan W. C. Leigh
Abstract:
Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass rat…
▽ More
Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass ratio ($10^{-3}-1$), the ratio of the binary semimajor axis to the hardening radius ($10^{-4}-1$), the binary eccentricity ($0.0-0.9$) and the stellar mass ($0.3-3M_{\odot}$). This is a significant extension of the parameter space compared to previous work. We show that the encounter probabilities of all three events are well-described by the encounter cross section, which is proportional to the pericenter distance. The probability of full disruptions by supermassive black hole binaries can be enhanced by up to a factor of $40-50$ or suppressed by up to a factor of $10$, relative to that by single black holes, depending on the binary parameters. Relativistic effects are not important for the primary black hole mass $\leq 10^{7}M_{\odot}$, but can provide an additional enhancement of the full disruption probability by less than a factor of $2-3$ for higher primary black hole masses. We provide a fitting formula for the full disruption probability by the hard supermassive black hole binaries that works for a wide range of parameters. We also find that partial disruption events can occur multiple times before full disruptions or direct captures, and their probabilities can be greater than that of full disruption events by a factor of three. Because partial disruption events can induce stellar spins and mass loss and change the orbits, it can significantly affect the overall full disruption event rate and the shape of the light curves.
△ Less
Submitted 2 August, 2022; v1 submitted 15 February, 2022;
originally announced February 2022.
-
The mean free path approximation and stellar collisions in star clusters: Numerical exploration of the analytic rates and the role of perturbations on binary star mergers
Authors:
Bastián Reinoso,
Nathan W. C. Leigh,
Carlos M. Barrera-Retamal,
Dominik Schleicher,
Ralf S. Klessen,
Amelia M. Stutz
Abstract:
In this paper we compute predictions for the number of stellar collisions derived from analytic models based on the mean free path (MFP) approximation and compare them to the results of $N$-body simulations. Our goal is to identify the cluster conditions under which the MFP approximation remains valid. Adopting a range of particle numbers ($100\leq N\leq5000$) and different combinations of particl…
▽ More
In this paper we compute predictions for the number of stellar collisions derived from analytic models based on the mean free path (MFP) approximation and compare them to the results of $N$-body simulations. Our goal is to identify the cluster conditions under which the MFP approximation remains valid. Adopting a range of particle numbers ($100\leq N\leq5000$) and different combinations of particle masses and radii, we explore three different channels leading to stellar collisions, all of which are expected to occur in realistic stellar environments. At high densities, binaries form from isolated three-body interactions of single stars. Hence, we consider collisions between single stars and collisions involving binary stars, after they form in our simulations. For the latter, we consider two channels for mergers, namely direct stellar collisions during chaotic single-binary interactions and perturbation-driven mergers of binaries due to random walks in eccentricity approaching unity. In the densest systems considered here, a very massive object is formed at the cluster centre, causing local stellar orbits to become increasingly Keplerian and the assumptions going into our analytic model to break down. Before reaching this limit, we obtain excellent agreement between our theoretical predictions and the simulations: the analytic rates are typically accurate to within one standard deviation for the entire parameter space considered here, but the agreement is best for short integration times. Our results have direct implications for blue straggler formation in dense star clusters, and stellar mergers in galactic nuclei hosting massive black holes.
△ Less
Submitted 4 November, 2021;
originally announced November 2021.
-
Starfall: A heavy rain of stars in 'turning on' AGN
Authors:
B. McKernan,
K. E. S. Ford,
M. Cantiello,
M. J. Graham,
A. S. Jermyn,
N. W. C. Leigh,
T. Ryu,
D. Stern
Abstract:
As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and…
▽ More
As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and onto the SMBH. A large AGN loss cone ($θ_{\rm AGN,lc}$) can result from binary scatterings in the inner disk and yield tidal disruption events (TDEs). Signatures of starfall include optical/UV flares that rise in luminosity over time, particularly in the inner disk. If the SMBH mass is $M_{\rm SMBH} \ge 10^{8}M_{\odot}$, flares truncate abruptly and the star is swallowed. If $M_{\rm SMBH}<10^{8}M_{\odot}$, and if the infalling orbit lies within $θ_{\rm AGN,lc}$, the flare is followed by a TDE which can be prograde or retrograde relative to the AGN inner disk. Retrograde AGN TDEs are over-luminous and short-lived as in-plane ejecta collide with the inner disk and a lower AGN state follows. Prograde AGN TDEs add angular momentum to inner disk gas and so start off looking like regular TDEs but are followed by an AGN high state. Searches for such flare signatures test models of AGN `turn on', SMBH mass, as well as disk properties and the embedded population.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei
Authors:
Yihan Wang,
Barry McKernan,
Saavik Ford,
Rosalba Perna,
Nathan W. C. Leigh,
Mordecai-Mark Mac Low
Abstract:
Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-maj…
▽ More
Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries.
△ Less
Submitted 7 October, 2021;
originally announced October 2021.
-
Hot Jupiter formation in dense clusters: secular chaos in multi-planetary systems
Authors:
Yihan Wang,
Rosalba Perna,
Nathan W. C. Leigh,
Michael M. Shara
Abstract:
Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity ex…
▽ More
Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity excitation mechanisms such as the Zeipel-Lidov-Kozai (ZLK) effect and planet-planet scattering. Previous studies have shown that if there are two giant planets in the planetary system, close flybys can efficiently activate the ZLK mechanism, thus triggering high eccentricity tidal migration and ultimately form hot Jupiters in star clusters. Here we extend our previous study with a multi-planet (triple) system. We perform high precision, high-accuracy few-body simulations of stellar flybys and subsequent planetary migration within the perturbed planetary systems using the code {\tt SpaceHub}. Our simulations demonstrate that a single close flyby on a multi-planet system in a cluster can activate secular chaos and ultimately lead to hot Jupiter formation via high eccentricity migration. We find that the hot Jupiter formation rate per system increases with both the size of the planetary system as well as with the mass of the outer planet, and we quantify the relative formation fractions for a range of parameters. Hot Jupiters formed via secular chaos are expected to be accompanied by massive companions with very long periods. Our study further shows that this flyby-induced secular chaos is preferred in low-density clusters where multi-planet systems are more likely to survive, and that it contributes a significant fraction of the hot Jupiter formation in star clusters compared to the flyby-induced ZLK mechanism.
△ Less
Submitted 27 September, 2021;
originally announced September 2021.
-
Chaos in the vicinity of a singularity in the Three-Body Problem: The equilateral triangle experiment in the zero angular momentum limit
Authors:
Hugo D. Parischewsky,
Alessandro A. Trani,
Nathan W. C. Leigh
Abstract:
We present numerical simulations of the gravitational three-body problem, in which three particles lie at rest close to the vertices of an equilateral triangle. In the unperturbed problem, the three particles fall towards the center of mass of the system to form a three-body collision, or singularity, where the particles overlap in space and time. By perturbing the initial positions of the particl…
▽ More
We present numerical simulations of the gravitational three-body problem, in which three particles lie at rest close to the vertices of an equilateral triangle. In the unperturbed problem, the three particles fall towards the center of mass of the system to form a three-body collision, or singularity, where the particles overlap in space and time. By perturbing the initial positions of the particles, we are able to study chaos in the vicinity of the singularity. Here we cover both the singular region close to the unperturbed configuration and the binary-single scattering regime where one side of the triangle is very short compared to the other two. We make phase space plots to study the regular and ergodic subsets of our simulations and compare them with the outcomes expected from the statistical escape theory of the three-body problem. We further provide fits to the ergodic subset to characterize the properties of the left-over binaries. We identify the discrepancy between the statistical theory and the simulations in the regular subset of interactions, which only exhibits weak chaos. As we decrease the scale of the perturbations in the initial positions, the phase space becomes entirely dominated by regular interactions, according to our metric for chaos. Finally, we show the effect of general relativity corrections by simulating the same scenario with the inclusion of post-Newtonian corrections to the equations of motion.
△ Less
Submitted 7 November, 2023; v1 submitted 13 August, 2021;
originally announced August 2021.
-
Understanding links between water-quality variables and nitrate concentration in freshwater streams using high-frequency sensor data
Authors:
Claire Kermorvant,
Benoit Liquet,
Guy Litt,
Kerrie Mengersen,
Erin Peterson,
Rob Hyndman,
Jeremy B. Jones Jr.,
Catherine Leigh
Abstract:
Real time monitoring using in situ sensors is becoming a common approach for measuring water quality within watersheds. High frequency measurements produce big data sets that present opportunities to conduct new analyses for improved understanding of water quality dynamics and more effective management of rivers and streams. Of primary importance is enhancing knowledge of the relationships between…
▽ More
Real time monitoring using in situ sensors is becoming a common approach for measuring water quality within watersheds. High frequency measurements produce big data sets that present opportunities to conduct new analyses for improved understanding of water quality dynamics and more effective management of rivers and streams. Of primary importance is enhancing knowledge of the relationships between nitrate, one of the most reactive forms of inorganic nitrogen in the aquatic environment, and other water quality variables. We analysed high frequency water quality data from in situ sensors deployed in three sites from different watersheds and climate zones within the National Ecological Observatory Network, USA. We used generalised additive mixed models to explain the nonlinear relationships at each site between nitrate concentration and conductivity, turbidity, dissolved oxygen, water temperature, and elevation. Temporal auto correlation was modelled with an auto regressive moving average model and we examined the relative importance of the explanatory variables. Total deviance explained by the models was high for all sites. Although variable importance and the smooth regression parameters differed among sites, the models explaining the most variation in nitrate contained the same explanatory variables. This study demonstrates that building a model for nitrate using the same set of explanatory water quality variables is achievable, even for sites with vastly different environmental and climatic characteristics. Applying such models will assist managers to select cost effective water quality variables to monitor when the goals are to gain a spatially and temporally in depth understanding of nitrate dynamics and adapt management plans accordingly.
△ Less
Submitted 3 June, 2021;
originally announced June 2021.
-
Stellar collisions in flattened and rotating Pop. III star clusters
Authors:
M. Z. C. Vergara,
D. R. G. Schleicher,
T. C. N. Boekholt,
B. Reinoso,
M. Fellhauer,
R. S. Klessen,
N. W. C. Leigh
Abstract:
Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive c…
▽ More
Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive cluster, in which the most massive star moves to the center of the cluster due to dynamical friction. This increases the rate of collisions and mergers since massive stars have larger collisional cross sections. This can lead to runaway growth of a very massive star which may collapse to become an intermediate-mass black hole. Here we investigate the dynamical evolution of Miyamoto-Nagai models that allow us to describe dense stellar clusters, including flattening and different degrees of rotation. We find that the collisions in these clusters depend mostly on the number of stars and the initial stellar radii for a given radial size of the cluster. By comparison, rotation seems to affect the collision rate by at most $20\%$. For flatness, we compared spherical models with systems that have a scale height of about $10\%$ of their radial extent, in this case finding a change in the collision rate of less than $25\%$. Overall, we conclude that the parameters only have a minor effect on the number of collisions. Our results also suggest that rotation helps to retain more stars in the system, reducing the number of escapers by a factor of $2-3$ depending on the model and the specific realization. After two million years, a typical lifetime of a very massive star, we find that about $630$ collisions occur in typical models with $N=10^4$, $R=100$ $\rm~R_\odot$ and a half-mass radius of $0.1$ $\rm~pc$, leading to a mass of about $6.3\times10^3$ $\rm~M_\odot$ for the most massive object.
△ Less
Submitted 3 April, 2021;
originally announced April 2021.
-
Spin misalignment of black hole binaries from young star clusters: implications for the origin of gravitational waves events
Authors:
Alessandro Alberto Trani,
Ataru Tanikawa,
Michiko S. Fujii,
Nathan W. C. Leigh,
Jun Kumamoto
Abstract:
Recent studies indicate that the progenitors of merging black hole (BH) binaries from young star clusters can undergo a common envelope phase just like isolated binaries. If the stars emerge from the common envelope as naked cores, tidal interactions can efficiently synchronize their spins before they collapse into BHs. Contrary to the isolated case, these binary BHs can also undergo dynamical int…
▽ More
Recent studies indicate that the progenitors of merging black hole (BH) binaries from young star clusters can undergo a common envelope phase just like isolated binaries. If the stars emerge from the common envelope as naked cores, tidal interactions can efficiently synchronize their spins before they collapse into BHs. Contrary to the isolated case, these binary BHs can also undergo dynamical interactions with other BHs in the cluster before merging. The interactions can tilt the binary orbital plane, leading to spin-orbit misalignment. We estimate the spin properties of merging binary BHs undergoing this scenario by combining up-to-date binary population synthesis and accurate few-body simulations. We show that post-common envelope binary BHs are likely to undergo only a single encounter, due to the high binary recoil velocity and short coalescence times. Adopting conservative limits on the binary-single encounter rates, we obtain a local BH merger rate density of ~6.6 yr^-1 Gpc^-3. Assuming low (<0.2) natal BH spins, this scenario reproduces the trends in the distributions of effective spin Xeff and precession parameters Xp inferred from GWTC-2, including the peaks at (Xeff, Xp) ~ (0.1, 0.2) and the tail at negative Xeff values.
△ Less
Submitted 2 April, 2021; v1 submitted 2 February, 2021;
originally announced February 2021.
-
Testing the Flux-based statistical prediction of the Three-Body Problem
Authors:
Viraj Manwadkar,
Barak Kol,
Alessandro A. Trani,
Nathan W. C. Leigh
Abstract:
We present an extensive comparison between the statistical properties of non-hierarchical three-body systems and the corresponding three-body theoretical predictions. We perform and analyze 1 million realizations for each different initial condition considering equal and unequal mass three-body systems to provide high accuracy statistics. We measure 4 quantities characterizing the statistical dist…
▽ More
We present an extensive comparison between the statistical properties of non-hierarchical three-body systems and the corresponding three-body theoretical predictions. We perform and analyze 1 million realizations for each different initial condition considering equal and unequal mass three-body systems to provide high accuracy statistics. We measure 4 quantities characterizing the statistical distribution of ergodic disintegrations: escape probability of each body, the characteristic exponent for escapes by a narrow margin, predicted absorptivity as a function of binary energy and binary angular momentum, and, finally, the lifetime distribution. The escape probabilities are shown to be in agreement down to the 1% level with the emissivity-blind, flux-based theoretical prediction. This represents a leap in accuracy compared to previous three-body statistical theories. The characteristic exponent at the threshold for marginally unbound escapes is an emissivity-independent flux-based prediction, and the measured values are found to agree well with the prediction. We interpret both tests as strong evidence for the flux-based three-body statistical formalism. The predicted absorptivity and lifetime distributions are measured to enable future tests of statistical theories.
△ Less
Submitted 10 January, 2021;
originally announced January 2021.
-
Formation of supermassive black hole seeds in nuclear star clusters via gas accretion and runaway collisions
Authors:
Arpan Das,
Dominik R. G. Schleicher,
Nathan W. C. Leigh,
Tjarda C. N. Boekholt
Abstract:
More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. I…
▽ More
More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. In this paper, we explore how SMSs with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealised N-body simulations. We explore physically motivated accretion scenarios, e.g. Bondi-Hoyle-Lyttleton accretion and Eddington accretion, as well as simplified scenarios such as constant accretions. While gas is present, the accretion timescale remains considerably shorter than the timescale for collisions with the most massive object (MMO). However, overall the timescale for collisions between any two stars in the cluster can become comparable or shorter than the accretion timescale, hence collisions still play a crucial role in determining the final mass of the SMSs. We find that the problem is highly sensitive to the initial conditions and our assumed recipe for the accretion, due to the highly chaotic nature of the problem. The key variables that determine the mass growth mechanism are the mass of the MMO and the gas reservoir that is available for the accretion. Depending on different conditions, SMSs of masses $\sim10^{3-5} \,\mathrm{M_{\odot}}$ can form for all three accretion scenarios considered in this work.
△ Less
Submitted 14 May, 2021; v1 submitted 2 December, 2020;
originally announced December 2020.
-
A phylogenetic analysis of galaxies in the Coma Cluster and the field: a new approach to galaxy evolution
Authors:
M. Martínez-Marín,
R. Demarco,
G. Cabrera-Vives,
P. Cerulo,
N. W. C. Leigh,
R. Herrera-Camus
Abstract:
We propose a phylogenetic approach (PA) as a novel and robust tool to detect galaxy populations (GPs) based on their chemical composition. The branches of the tree are interpreted as different GPs and the length between nodes as the internal chemical variation along a branch. We apply the PA using 30 abundance indices from the Sloan Digital Sky Survey to 475 galaxies in the Coma Cluster and 438 ga…
▽ More
We propose a phylogenetic approach (PA) as a novel and robust tool to detect galaxy populations (GPs) based on their chemical composition. The branches of the tree are interpreted as different GPs and the length between nodes as the internal chemical variation along a branch. We apply the PA using 30 abundance indices from the Sloan Digital Sky Survey to 475 galaxies in the Coma Cluster and 438 galaxies in the field. We find that a dense environment, such as Coma, shows several GPs, which indicates that the environment is promoting galaxy evolution. Each population shares common properties that can be identified in colour magnitude space, in addition to minor structures inside the red sequence. The field is more homogeneous, presenting one main GP. We also apply a principal component analysis (PCA) to both samples, and find that the PCA does not have the same power in identifying GPs.
△ Less
Submitted 13 November, 2020;
originally announced November 2020.
-
Hot Jupiter and ultra-cold Saturn formation in dense star clusters
Authors:
Yi-Han Wang,
Nathan W. C. Leigh,
Rosalba Perna,
Michael M. Shara
Abstract:
The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequen…
▽ More
The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy.
△ Less
Submitted 2 November, 2020;
originally announced November 2020.
-
Small-N Collisional Dynamics V: Beyond the Realm of Not-So-Small-N
Authors:
Carlos Barrera,
Nathan W. C. Leigh,
Bastián Reinoso,
Amelia M. Stutz,
Dominik Schleicher
Abstract:
Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important application is stellar collisions, which occur commonly in dense star clusters, and their relevance for the formation of various types of stellar exotica. In this paper,…
▽ More
Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important application is stellar collisions, which occur commonly in dense star clusters, and their relevance for the formation of various types of stellar exotica. In this paper, we return to our study of the collision rates and probabilities during small-number chaotic gravitational interactions ($N$ $\lesssim$ 10), moving beyond the small-number particle limit and into the realm of larger particle numbers ($N$ $\gtrsim$ 10$^3$) to test the extent of validity of our analytic model as a function of the particle properties and the number of interacting particles. This is done using direct $N$-body simulations of stellar collisions in dense star clusters, by varying the relative numbers of particles with different particle masses and radii. We compute the predicted rate of collisions using the mean free path approximation, adopting the point-particle limit and using the sticky-star approximation as our collision criterion. We evaluate its efficacy in the regime where gravitational-focusing is important by comparing the theoretical rates to numerical simulations. Using the tools developed in previous papers in this series, in particular Collision Rate Diagrams, we illustrate that our predicted and simulated rates are in excellent agreement, typically consistent with each other to within one standard deviation.
△ Less
Submitted 30 October, 2020;
originally announced November 2020.
-
Massive Stellar Triples Leading to Sequential Binary Black-Hole Mergers in the Field
Authors:
Alejandro Vigna-Gómez,
Silvia Toonen,
Enrico Ramirez-Ruiz,
Nathan W. C. Leigh,
Jeff Riley,
Carl-Johan Haster
Abstract:
Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are…
▽ More
Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well.
△ Less
Submitted 22 December, 2020; v1 submitted 26 October, 2020;
originally announced October 2020.
-
Filament rotation in the California L1482 cloud
Authors:
Rodrigo H. Álvarez-Gutiérrez,
Amelia M. Stutz,
Chi Y. Law,
Stefan Reissl,
Ralf S. Klessen,
Nathan W. C. Leigh,
Hong-Li Liu,
Rodrigo A. Reeves
Abstract:
We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C…
▽ More
We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation.
△ Less
Submitted 17 December, 2020; v1 submitted 21 October, 2020;
originally announced October 2020.
-
AGN as Potential Factories for Eccentric Black Hole Mergers
Authors:
J. Samsing,
I. Bartos,
D. J. D'Orazio,
Z. Haiman,
B. Kocsis,
N. W. C. Leigh,
B. Liu,
M. E. Pessah,
H. Tagawa
Abstract:
There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to th…
▽ More
There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to the merger, as basic physics would argue for circularization. Here we show that AGN-disk environments can lead to an excess of eccentric mergers, if the interactions between single and binary black holes are frequent, and occur with mutual inclinations of less than a few degrees. We further illustrate that this eccentric population has a different distribution of the inclination between the spin vectors of the black holes and their orbital angular momentum at merger, referred to as the spin-orbit tilt, compared to the remaining circular mergers.
△ Less
Submitted 21 April, 2022; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Evolution of Retrograde Orbiters in an AGN Disk
Authors:
Amy Secunda,
Betsy Hernandez,
Jeremy Goodman,
Nathan W. C. Leigh,
Barry McKernan,
K. E. Saavik Ford,
Jose I. Adorno
Abstract:
AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifyin…
▽ More
AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifying the retrograde contribution is important, since roughly half of all orbiters should initially be on retrograde orbits when the disk forms. We perform an analytic calculation of the evolution of ROs in an AGN disk. Because this evolution could cause the orbits of ROs to cross those of prograde BBHs, we derive the collision rate between a given RO and a given BBH orbiting in the prograde direction. In the examples given here, ROs in the inner region of the disk experience a rapid decrease in the semimajor axis of their orbits while also becoming highly eccentric in less than a million years. This rapid orbital evolution could lead to extreme mass ratio inspirals detectable by the Laser Interferometer Space Antenna. The collision rates of our example ROs with prograde BBHs in the migration trap depend strongly on the volume of the inner radiation-pressure-dominated region which depends on the mass of the supermassive black hole (SMBH). Rates are lowest for larger mass SMBHs, which dominate the AGN merger channel, suggesting that merger rates for this channel may not be significantly altered by ROs.
△ Less
Submitted 4 February, 2021; v1 submitted 8 September, 2020;
originally announced September 2020.
-
The effects of a background potential in star cluster evolution: a delay in the relaxation time-scale and runaway collision processes
Authors:
B. Reinoso,
D. R. G. Schleicher,
M. Fellhauer,
N. W. C. Leigh,
R. S. Klessen
Abstract:
Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explic…
▽ More
Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times.
△ Less
Submitted 15 May, 2020;
originally announced May 2020.
-
Mergers of Equal-Mass Binaries with Compact Object Companions from Mass Transfer in Triple Star Systems
Authors:
N. W. C. Leigh,
S. Toonen,
S. F. Portegies Zwart,
R. Perna
Abstract:
In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analyti…
▽ More
In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion.
△ Less
Submitted 23 June, 2020; v1 submitted 13 May, 2020;
originally announced May 2020.
-
Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes II. Spins and Incoming Objects
Authors:
Amy Secunda,
Jillian Bellovary,
Mordecai-Mark Mac Low,
K. E. Saavik Ford,
Barry McKernan,
Nathan W. C. Leigh,
Wladimir Lyra,
Zsolt Sandor,
Jose I. Adorno
Abstract:
The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previou…
▽ More
The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. Orbiters entering the inner disk form BBHs with sBHs on resonant orbits near the migration trap. When these sBHs reach ~80 Msun, they form BBHs with sBHs in the migration trap, which over 10 Myr reach ~1000 Msun. We find 68% of the BBHs in our simulation orbit in the retrograde direction, which implies BBHs in our merger channel will have small dimensionless aligned spins, χ_eff. Overall, our models produce BBHs that resemble both the majority of BBH mergers detected thus far (0.66 to 120 Gpc^-3 yr^-1 ) and two recent unusual detections, GW190412 (~0.3 Gpc^-3 yr^-1 ) and GW190521 (~0.1 Gpc^-3 yr^-1 ).
△ Less
Submitted 25 September, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
Chaos and Lévy Flights in the Three-Body Problem
Authors:
Viraj Manwadkar,
Alessandro A. Trani,
Nathan W. C. Leigh
Abstract:
We study chaos and Lévy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the…
▽ More
We study chaos and Lévy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the ergodic interactions and Lévy flight interactions. This enables us to study the three-body lifetime distribution in more detail by decomposing it into the individual distributions from the different kinds of interactions. We observe that ergodic interactions follow an exponential decay distribution similar to that of radioactive decay. Meanwhile, Lévy flight interactions follow a power-law distribution. Lévy flights in fact dominate the tail of the general three-body lifetime distribution, providing conclusive evidence for the speculated connection between power-law tails and Lévy flight interactions. We propose a new physically-motivated model for the lifetime distribution of three-body systems and discuss how it can be used to extract information about the underlying ergodic and Lévy flight interactions. We discuss mass ejection probabilities in three-body systems in the ergodic limit and compare it to previous ergodic formalisms. We introduce a novel mechanism for a three-body relaxation process and discuss its relevance in general three-body systems.
△ Less
Submitted 19 June, 2020; v1 submitted 11 April, 2020;
originally announced April 2020.
-
Giant planet swaps during close stellar encounters
Authors:
Yi-Han Wang,
Rosalba Perna,
Nathan W. C. Leigh
Abstract:
The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not…
▽ More
The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $Γ\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (σ/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $σ$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
Planetary Architectures in Interacting Stellar Environments
Authors:
Yi-Han Wang,
Rosalba Perna,
Nathan W. C. Leigh
Abstract:
The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers…
▽ More
The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers (stars and binaries), we explore the reconfiguration to the planetary system due to the encounters. In particular, via an exploration focused on the strong scattering regime, we derive the fraction of encounters which result in planet ejections, planet transfers and collisions by the interloper star/binary, as a function of the characteristics of the environment (density, velocity dispersion), and for different masses of the flyby star/binary. We find that binary interlopers can significantly increase the cross section of planet ejections and collisions, while they only slightly change the cross section for planet transfers. Therefore, in environments with high binary fractions, floating planets are expected to be relatively common, while in environments with low binary fractions, where the cross sections of planet ejection and transfer are comparable, the rate of planet exchanges between two stars will be comparable to the rate of production of free-floating planets.
△ Less
Submitted 13 February, 2020;
originally announced February 2020.
-
The X-ray Emissivity of Low-Density Stellar Populations
Authors:
C. O. Heinke,
M. G. Ivanov,
E. W. Koch,
R. Andrews,
L. Chomiuk,
H. N. Cohn,
S. Crothers,
T. de Boer,
N. Ivanova,
A. K. H. Kong,
N. W. C. Leigh,
P. M. Lugger,
L. Nelson,
C. J. Parr,
E. W. Rosolowsky,
A. J. Ruiter,
C. L. Sarazin,
A. W. Shaw,
G. R. Sivakoff,
M. van den Berg
Abstract:
The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations,…
▽ More
The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.
△ Less
Submitted 26 January, 2020;
originally announced January 2020.
-
On the eccentricity evolution of massive black hole binaries in stellar backgrounds
Authors:
Matteo Bonetti,
Alexander Rasskazov,
Alberto Sesana,
Massimo Dotti,
Francesco Haardt,
Nathan W. C. Leigh,
Manuel Arca Sedda,
Giacomo Fragione,
Elena Rossi
Abstract:
We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrap…
▽ More
We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).
△ Less
Submitted 4 February, 2020; v1 submitted 7 January, 2020;
originally announced January 2020.
-
SSNdesign -- an R package for pseudo-Bayesian optimal and adaptive sampling designs on stream networks
Authors:
Alan R. Pearse,
James M. McGree,
Nicholas A. Som,
Catherine Leigh,
Jay M. Ver Hoef,
Paul Maxwell,
Erin E. Peterson
Abstract:
Streams and rivers are biodiverse and provide valuable ecosystem services. Maintaining these ecosystems is an important task, so organisations often monitor the status and trends in stream condition and biodiversity using field sampling and, more recently, autonomous in-situ sensors. However, data collection is often costly and so effective and efficient survey designs are crucial to maximise info…
▽ More
Streams and rivers are biodiverse and provide valuable ecosystem services. Maintaining these ecosystems is an important task, so organisations often monitor the status and trends in stream condition and biodiversity using field sampling and, more recently, autonomous in-situ sensors. However, data collection is often costly and so effective and efficient survey designs are crucial to maximise information while minimising costs. Geostatistics and optimal and adaptive design theory can be used to optimise the placement of sampling sites in freshwater studies and aquatic monitoring programs. Geostatistical modelling and experimental design on stream networks pose statistical challenges due to the branching structure of the network, flow connectivity and directionality, and differences in flow volume. Thus, unique challenges of geostatistics and experimental design on stream networks necessitates the development of new open-source software for implementing the theory. We present SSNdesign, an R package for solving optimal and adaptive design problems on stream networks that integrates with existing open-source software. We demonstrate the mathematical foundations of our approach, and illustrate the functionality of SSNdesign using two case studies involving real data from Queensland, Australia. In both case studies we demonstrate that the optimal or adaptive designs outperform random and spatially balanced survey designs. The SSNdesign package has the potential to boost the efficiency of freshwater monitoring efforts and provide much-needed information for freshwater conservation and management.
△ Less
Submitted 1 December, 2019;
originally announced December 2019.
-
A Statistical Solution to the Chaotic, Non-Hierarchical Three-Body Problem
Authors:
Nicholas C. Stone,
Nathan W. C. Leigh
Abstract:
The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost a…
▽ More
The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost always disintegrate into a single escaping star and a stable, bound binary, but the chaotic nature of the three-body problem prevents the derivation of tractable analytic formulae deterministically mapping initial conditions to final outcomes. However, chaos also motivates the assumption of ergodicity, suggesting that the distribution of outcomes is uniform across the accessible phase volume. Here, we use the ergodic hypothesis to derive a complete statistical solution to the non-hierarchical three-body problem, one which provides closed-form distributions of outcomes (e.g. binary orbital elements) given the conserved integrals of motion. We compare our outcome distributions to large ensembles of numerical three-body integrations, and find good agreement, so long as we restrict ourselves to "resonant" encounters (the ~50% of scatterings that undergo chaotic evolution). In analyzing our scattering experiments, we identify "scrambles" (periods in time where no pairwise binaries exist) as the key dynamical state that ergodicizes a non-hierarchical triple. The generally super-thermal distributions of survivor binary eccentricity that we predict have notable applications to many astrophysical scenarios. For example, non-hierarchical triples produced dynamically in globular clusters are a primary formation channel for black hole mergers, but the rates and properties of the resulting gravitational waves depend on the distribution of post-disintegration eccentricities.
△ Less
Submitted 11 September, 2019;
originally announced September 2019.
-
Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems
Authors:
Giacomo Fragione,
Brian D. Metzger,
Rosalba Perna,
Nathan W. C. Leigh,
Bence Kocsis
Abstract:
Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effect…
▽ More
Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.
△ Less
Submitted 27 April, 2020; v1 submitted 2 August, 2019;
originally announced August 2019.
-
Kinematic study of the association Cyg OB3 with Gaia DR2
Authors:
Anjali Rao,
Poshak Gandhi,
Christian Knigge,
John A. Paice,
Nathan W. C. Leigh,
Douglas Boubert
Abstract:
We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean pec…
▽ More
We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is ~22 km/s, dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds <20 km/s with respect to the association for a majority of the sample stars.
△ Less
Submitted 1 May, 2020; v1 submitted 2 August, 2019;
originally announced August 2019.
-
The evolution of kicked stellar-mass black holes in star cluster environments II. Rotating star clusters
Authors:
Jeremy J. Webb,
Nathan W. C. Leigh,
Roberto Serrano,
Jillian Bellovary,
K. E. Saavik Ford,
Barry McKernan,
Mario Spera,
Alessandro A. Trani
Abstract:
In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic fra…
▽ More
In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.
△ Less
Submitted 9 July, 2019;
originally announced July 2019.
-
The Keplerian three-body encounter II. Comparisons with isolated encounters and impact on gravitational wave merger timescales
Authors:
Alessandro A. Trani,
Mario Spera,
Nathan W. C. Leigh,
Michiko S. Fujii
Abstract:
We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects t…
▽ More
We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects the breakup of triples by tidally limiting the system, so that the triples break earlier with lower breakup velocity, leaving behind slightly larger binaries (relative to the isolated case). The breakup direction becomes anisotropic and tends to follow the shape of the Hill region of the triple, favouring breakups in the radial direction. Furthermore, the tidal field can torque the system, leading to angular momentum exchanges between the triple and its orbit about the SMBH. This process changes the properties of the final binary, depending on the initial angular momentum of the triple. Finally, the tidal field also affects binary-single encounters: binaries tend to become both harder and more eccentric with respect to encounters that occur in isolation. Consequently, single-binary scattering in a deep Keplerian potential produces binaries with shorter gravitational wave merger timescales.
△ Less
Submitted 24 September, 2019; v1 submitted 16 April, 2019;
originally announced April 2019.
-
Constraining Mass Transfer Histories of Blue Straggler Stars with COS Spectroscopy of White Dwarf Companions
Authors:
Natalie M. Gosnell,
Emily M. Leiner,
Robert D. Mathieu,
Aaron M. Geller,
Christian Knigge,
Alison Sills,
Nathan W. C. Leigh
Abstract:
Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with…
▽ More
Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with HST COS far-ultraviolet spectroscopy. Adopting a Gaia-based cluster distance of $1847\pm107$ pc, we determine that one system, WOCS 4540, formed through Case C mass transfer resulting in a CO-core white dwarf with $T_{\text{eff}}=17000^{+140}_{-200}$ K and a $\log g=7.80^{+0.06}_{-0.06}$, corresponding to a mass of $0.53^{+0.03}_{-0.03}$ $M_{\odot}$ and a cooling age of $105^{+6}_{-5}$ Myr. The other system, WOCS 5379, formed through Case B mass transfer resulting in a He-core white dwarf with $T_{\text{eff}}=15500^{+170}_{-150}$ K and a $\log g=7.50^{+0.06}_{-0.05}$, corresponding to a mass of $0.42^{+0.02}_{-0.02}$ $M_{\odot}$ and an age of $250^{+20}_{-20}$ Myr. The WD parameters are consistent across four different cluster distance assumptions. We determine possible progenitor binary systems with a grid of accretion models using MESA, and investigate whether these systems would lead to stable or unstable mass transfer. WOCS 4540 likely resulted from stable mass transfer during periastron passage in an eccentric binary system, while WOCS 5379 challenges our current understanding of the expected outcomes for mass transfer from red giant branch stars. Both systems are examples of the value in using detailed analyses to fine-tune our physical understanding of binary evolutionary processes.
△ Less
Submitted 19 August, 2019; v1 submitted 3 April, 2019;
originally announced April 2019.
-
A feature-based framework for detecting technical outliers in water-quality data from in situ sensors
Authors:
Priyanga Dilini Talagala,
Rob J. Hyndman,
Catherine Leigh,
Kerrie Mengersen,
Kate Smith-Miles
Abstract:
Outliers due to technical errors in water-quality data from in situ sensors can reduce data quality and have a direct impact on inference drawn from subsequent data analysis. However, outlier detection through manual monitoring is unfeasible given the volume and velocity of data the sensors produce. Here, we proposed an automated framework that provides early detection of outliers in water-quality…
▽ More
Outliers due to technical errors in water-quality data from in situ sensors can reduce data quality and have a direct impact on inference drawn from subsequent data analysis. However, outlier detection through manual monitoring is unfeasible given the volume and velocity of data the sensors produce. Here, we proposed an automated framework that provides early detection of outliers in water-quality data from in situ sensors caused by technical issues.The framework was used first to identify the data features that differentiate outlying instances from typical behaviours. Then statistical transformations were applied to make the outlying instances stand out in transformed data space. Unsupervised outlier scoring techniques were then applied to the transformed data space and an approach based on extreme value theory was used to calculate a threshold for each potential outlier. Using two data sets obtained from in situ sensors in rivers flowing into the Great Barrier Reef lagoon, Australia, we showed that the proposed framework successfully identified outliers involving abrupt changes in turbidity, conductivity and river level, including sudden spikes, sudden isolated drops and level shifts, while maintaining very low false detection rates. We implemented this framework in the open source R package oddwater.
△ Less
Submitted 17 February, 2019;
originally announced February 2019.
-
In Search of the Thermal Eccentricity Distribution
Authors:
Aaron M. Geller,
Nathan W. C. Leigh,
Mirek Giersz,
Kyle Kremer,
Frederic A. Rasio
Abstract:
About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary…
▽ More
About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.
△ Less
Submitted 31 January, 2019;
originally announced February 2019.
-
A Triple Origin for Twin Blue Stragglers in Close Binaries
Authors:
Simon Portegies Zwart,
Nathan W. C. Leigh
Abstract:
We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed…
▽ More
We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed model. We predict the following properties for twin BSs: (1) For the outer tertiary orbit, the initial orbital period should lie between 220 days $\lesssim$ P$_{\rm out}$ $\lesssim$ 1100 days, assuming initial masses for the inner binary components of $m_{\rm 1} = 1.1$ M$_{\odot}$ and $m_{\rm 2} =$ 0.9 M$_{\odot}$ and an outer tertiary mass of $m_{\rm 3} = 1.4$ M$_{\odot}$. After Roche-lobe overflow, the outer star turns into a white dwarf (WD) of mass 0.43 to 0.54\,\MSun. There is a correlation between the mass of this WD and the outer orbital period: more massive WDs will be on wider orbits. (3) The rotational axes of both BSs will be aligned with each other and the orbital plane of the outer tertiary WD. (4) The BSs will have roughly equal masses, independent of their initial masses (since the lower mass star accretes the most). The dominant accretor should, therefore, be more enriched by the accreted material. Hence, one of the BSs will appear to be more enriched by either He, C and O or by s-process elements, if the donor started Roche lobe overflow on, respectively, the red giant or asymptotic giant branch. (5) Relative to old clusters, twin BSs in close binaries formed from the proposed mechanism should be more frequent in the Galactic field and open clusters with ages $\lesssim$ 4-6 Gyr, since then the donor will have a radiative envelope. (6) The orbit of the binary BS will have a small semi-major axis (typically $\aplt 0.3$\,au) and be close to circular ($e \aplt 0.2$).
△ Less
Submitted 30 January, 2019;
originally announced January 2019.
-
Black hole and neutron star mergers in galactic nuclei
Authors:
Giacomo Fragione,
Evgeni Grishin,
Nathan W. C. Leigh,
Hagai. B. Perets,
Rosalba Perna
Abstract:
Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies,…
▽ More
Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.
△ Less
Submitted 12 June, 2019; v1 submitted 26 November, 2018;
originally announced November 2018.