Showing 1–50 of 150 results for author: Stone, C

Caustics: A Python Package for Accelerated Strong Gravitational Lensing Simulations

Authors: Connor Stone, Alexandre Adam, Adam Coogan, M. J. Yantovski-Barth, Andreas Filipp, Landung Setiawan, Cordero Core, Ronan Legin, Charles Wilson, Gabriel Missael Barco, Yashar Hezaveh, Laurence Perreault-Levasseur

Abstract: Gravitational lensing is the deflection of light rays due to the gravity of intervening masses. This phenomenon is observed in a variety of scales and configurations, involving any non-uniform mass such as planets, stars, galaxies, clusters of galaxies, and even the large scale structure of the universe. Strong lensing occurs when the distortions are significant and multiple images of the backgrou… ▽ More Gravitational lensing is the deflection of light rays due to the gravity of intervening masses. This phenomenon is observed in a variety of scales and configurations, involving any non-uniform mass such as planets, stars, galaxies, clusters of galaxies, and even the large scale structure of the universe. Strong lensing occurs when the distortions are significant and multiple images of the background source are observed. The lens objects must align on the sky of order ~1 arcsecond for galaxy-galaxy lensing, or 10's of arcseonds for cluster-galaxy lensing. As the discovery of lens systems has grown to the low thousands, these systems have become pivotal for precision measurements and addressing critical questions in astrophysics. Notably, they facilitate the measurement of the Universe's expansion rate, dark matter, supernovae, quasars, and the first stars among other topics. With future surveys expected to discover hundreds of thousands of lensing systems, the modelling and simulation of such systems must occur at orders of magnitude larger scale then ever before. Here we present `caustics`, a Python package designed to handle the extensive computational demands of modeling such a vast number of lensing systems. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 13 pages, 3 figures, submitted to JOSS

The Nearly Universal Disk Galaxy Rotation Curve

Authors: Raj Patel, Nikhil Arora, Stéphane Courteau, Connor Stone, Matthew Frosst, Lawrence Widrow

Abstract: The Universal Rotation Curve (URC) of disk galaxies was originally proposed to predict the shape and amplitude of any rotation curve (RC) based solely on photometric data. Here, the URC is investigated with an extensive set of spatially-resolved rotation curves drawn from the PROBES-I, PROBES-II, and MaNGA data bases with matching multi-band surface brightness profiles from the DESI-LIS and WISE s… ▽ More The Universal Rotation Curve (URC) of disk galaxies was originally proposed to predict the shape and amplitude of any rotation curve (RC) based solely on photometric data. Here, the URC is investigated with an extensive set of spatially-resolved rotation curves drawn from the PROBES-I, PROBES-II, and MaNGA data bases with matching multi-band surface brightness profiles from the DESI-LIS and WISE surveys for 3,846 disk galaxies. Common URC formulations fail to achieve an adequate level of accuracy to qualify as truly universal over fully sampled RCs. We develop neural network (NN) equivalents for the proposed URCs which predict RCs with higher accuracy, showing that URC inaccuracies are not due to insufficient data but rather non-optimal formulations or sampling effects. This conclusion remains even if the total RC sample is pruned for symmetry. The latest URC prescriptions and their NN equivalents trained on our sub-sample of 579 disk galaxies with symmetric RCs perform similarly to the URC/NN trained on the complete data sample. We conclude that a URC with an acceptable level of accuracy ($ΔV_{\rm circ} \lesssim15$ per cent) at all radii would require a detailed modelling of a galaxy's central regions and outskirts (e.g., for baryonic effects leading to contraction or expansion of any dark-matter-only halo). △ Less

Submitted 19 June, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: 36 full pages, 13 figures, 2 appendices. Accepted for publication in the Astrophysical Journal

Elevated Rates of Tidal Disruption Events in Active Galactic Nuclei

Authors: Karamveer Kaur, Nicholas C. Stone

Abstract: Advances in time domain astronomy have produced a growing population of flares from galactic nuclei, including both tidal disruption events (TDEs) and flares in active galactic nuclei (AGN). Because TDEs are uncommon and AGN variability is abundant, large-amplitude AGN flares are usually not categorized as TDEs. While TDEs are normally channelled by the collisional process of two-body scatterings… ▽ More Advances in time domain astronomy have produced a growing population of flares from galactic nuclei, including both tidal disruption events (TDEs) and flares in active galactic nuclei (AGN). Because TDEs are uncommon and AGN variability is abundant, large-amplitude AGN flares are usually not categorized as TDEs. While TDEs are normally channelled by the collisional process of two-body scatterings over relaxation timescale, the quadrupole moment of a gas disk alters the stellar orbits, allowing them to collisionlessly approach the central massive black hole (MBH). This leads to an effectively enlarged loss cone, the \emph{loss wedge}. Earlier studies found a moderate enhancement, up to a factor $\sim 2-3$, of TDE rates $\dot{N}_{\rm 2b} $ for a static axisymmetric perturbation. Here we study the loss wedge problem for an evolving AGN disk, which can capture large number of stars into the growing loss wedge over much shorter times. The rates $\dot{N}_{\rm cl}$ of collisionless TDEs produced by these time-evolving disks are much higher than the collisional rates $\dot{N}_{\rm 2b}$ in a static loss wedge. We calculate the response of a stellar population to the axisymmetric potential of an adiabatically growing AGN disk and find that the highest rates of collisionless TDEs are achieved for the largest (i) MBH masses $M_{\bullet}$ and (ii) disk masses $M_{\rm d}$. For $M_{\bullet}\sim 10^7 M_\odot$ and $M_{\rm d} \sim 0.1 M_{\bullet}$, the rate enhancement can be up to a factor $\dot{N}_{\rm cl}/\dot{N}_{\rm 2b} \sim 10$. The orbits of collisionless TDEs sometimes have a preferred orientation in apses, carrying implications for observational signatures of resulting flares. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: Submitted to ApJ; 24 pages, 8 figures; comments are welcome

Parker Solar Probe Observations of Energetic Particles in the Flank of a Coronal Mass Ejection Close to the Sun

Authors: N. A. Schwadron, Stuart D. Bale, J. Bonnell, A. Case, M. Shen, E. R. Christian, C. M. S. Cohen, A. J. Davis, M. I. Desai, K. Goetz, J. Giacalone, M. E. Hill, J. C. Kasper, K. Korreck, D. Larson, R. Livi, T. Lim, R. A. Leske, O. Malandraki, D. Malaspina, W. H. Matthaeus, D. J. McComas, R. L. McNutt Jr., R. A. Mewaldt, D. G. Mitchell , et al. (10 additional authors not shown)

Abstract: We present an event observed by Parker Solar Probe at $\sim$0.2 au on March 2, 2022 in which imaging and \emph{in situ} measurements coincide. During this event, PSP passed through structures on the flank of a streamer blowout CME including an isolated flux tube in front of the CME, a turbulent sheath, and the CME itself. Imaging observations and \emph{in situ} helicity and principal variance sign… ▽ More We present an event observed by Parker Solar Probe at $\sim$0.2 au on March 2, 2022 in which imaging and \emph{in situ} measurements coincide. During this event, PSP passed through structures on the flank of a streamer blowout CME including an isolated flux tube in front of the CME, a turbulent sheath, and the CME itself. Imaging observations and \emph{in situ} helicity and principal variance signatures consistently show the presence of flux ropes internal to the CME. In both the sheath, and the CME interval, the distributions are more isotropic, the spectra are softer, and the abundance ratios of Fe/O and He/H are lower than those in the isolated flux tube, and yet elevated relative to typical plasma and SEP abundances. These signatures in the sheath and the CME indicate that both flare populations and those from the plasma are accelerated to form the observed energetic particle enhancements. In contrast, the isolated flux tube shows large streaming, hard spectra and large Fe/O and He/H ratios, indicating flare sources. Energetic particle fluxes are most enhanced within the CME interval from suprathermal through energetic particle energies ($\sim$ keV to $>10$ MeV), indicating particle acceleration, and confinement local to the closed magnetic structure. The flux-rope morphology of the CME helps to enable local modulation and trapping of energetic particles, particularly along helicity channels and other plasma boundaries. Thus, the CME acts to build-up energetic particle populations, allowing them to be fed into subsequent higher energy particle acceleration throughout the inner heliosphere where a compression or shock forms on the CME front. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: 41 pages, 19 figures, In Press

Journal ref: Astrophysical Journal, 2024

Euclid: ERO -- NISP-only sources and the search for luminous $z=6-8$ galaxies

Authors: J. R. Weaver, S. Taamoli, C. J. R. McPartland, L. Zalesky, N. Allen, S. Toft, D. B. Sanders, H. Atek, R. A. A. Bowler, D. Stern, C. J. Conselice, B. Mobasher, I. Szapudi, P. R. M. Eisenhardt, G. Murphree, I. Valdes, K. Ito, S. Belladitta, P. A. Oesch, S. Serjeant, D. J. Mortlock, N. A. Hatch, M. Kluge, B. Milvang-Jensen, G. Rodighiero , et al. (163 additional authors not shown)

Abstract: This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ (… ▽ More This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ ($M_{\rm UV}\lesssim-22$). Beyond their still uncertain role in reionisation, these UV-bright galaxies are ideal laboratories from which to study galaxy formation and constrain the bright-end of the UV luminosity function. Of the 501994 sources detected from a combined $Y_{\rm E}$, $J_{\rm E}$, and $H_{\rm E}$ NISP detection image, 168 do not have any appreciable VIS/$I_{\rm E}$ flux. These objects span a range in spectral colours, separated into two classes: 139 extremely red sources; and 29 Lyman-break galaxy candidates. Best-fit redshifts and spectral templates suggest the former is composed of both $z\gtrsim5$ dusty star-forming galaxies and $z\approx1-3$ quiescent systems. The latter is composed of more homogeneous Lyman break galaxies at $z\approx6-8$. In both cases, contamination by L- and T-type dwarfs cannot be ruled out with Euclid images alone. Additional contamination from instrumental persistence is investigated using a novel time series analysis. This work lays the foundation for future searches within the Euclid Deep Fields, where thousands more $z\gtrsim6$ Lyman break systems and extremely red sources will be identified. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures, paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- A preview of the Euclid era through a galaxy cluster magnifying lens

Authors: H. Atek, R. Gavazzi, J. R. Weaver, J. M. Diego, T. Schrabback, N. A. Hatch, N. Aghanim, H. Dole, W. G. Hartley, S. Taamoli, G. Congedo, Y. Jimenez-Teja, J. -C. Cuillandre, E. Bañados, S. Belladitta, R. A. A. Bowler, M. Franco, M. Jauzac, G. Mahler, J. Richard, P. -F. Rocci, S. Serjeant, S. Toft, D. Abriola, P. Bergamini , et al. (178 additional authors not shown)

Abstract: We present the first analysis of the Euclid Early Release Observations (ERO) program that targets fields around two lensing clusters, Abell 2390 and Abell 2764. We use VIS and NISP imaging to produce photometric catalogs for a total of $\sim 500\,000$ objects. The imaging data reach a $5\,σ$ typical depth in the range 25.1-25.4 AB in the NISP bands, and 27.1-27.3 AB in the VIS band. Using the Lyma… ▽ More We present the first analysis of the Euclid Early Release Observations (ERO) program that targets fields around two lensing clusters, Abell 2390 and Abell 2764. We use VIS and NISP imaging to produce photometric catalogs for a total of $\sim 500\,000$ objects. The imaging data reach a $5\,σ$ typical depth in the range 25.1-25.4 AB in the NISP bands, and 27.1-27.3 AB in the VIS band. Using the Lyman-break method in combination with photometric redshifts, we identify $30$ Lyman-break galaxy (LBG) candidates at $z>6$ and 139 extremely red sources (ERSs), most likely at lower redshift. The deeper VIS imaging compared to NISP means we can routinely identify high-redshift Lyman breaks of the order of $3$ magnitudes, which reduces contamination by brown dwarf stars and low-redshift galaxies. Spectroscopic follow-up campaigns of such bright sources will help constrain both the bright end of the ultraviolet galaxy luminosity function and the quasar luminosity function at $z>6$, and constrain the physical nature of these objects. Additionally, we have performed a combined strong lensing and weak lensing analysis of A2390, and demonstrate how Euclid will contribute to better constraining the virial mass of galaxy clusters. From these data, we also identify optical and near-infrared counterparts of known $z>0.6$ clusters, which exhibit strong lensing features, establishing the ability of Euclid to characterize high-redshift clusters. Finally, we provide a glimpse of Euclid's ability to map the intracluster light out to larger radii than current facilities, enabling a better understanding of the cluster assembly history and mapping of the dark matter distribution. This initial dataset illustrates the diverse spectrum of legacy science that will be enabled by the Euclid survey. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations. 17 pages, 12 figures

Euclid: Early Release Observations -- The intracluster light and intracluster globular clusters of the Perseus cluster

Authors: M. Kluge, N. A. Hatch, M. Montes, J. B. Golden-Marx, A. H. Gonzalez, J. -C. Cuillandre, M. Bolzonella, A. Lançon, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, A. Boselli, M. Cantiello, J. G. Sorce, F. R. Marleau, P. -A. Duc, E. Sola, M. Urbano, S. L. Ahad, Y. M. Bahé, S. P. Bamford, C. Bellhouse, F. Buitrago, P. Dimauro , et al. (163 additional authors not shown)

Abstract: We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus clu… ▽ More We study the intracluster light (ICL) and intracluster globular clusters (ICGCs) in the nearby Perseus galaxy cluster using Euclid's EROs. By modelling the isophotal and iso-density contours, we map the distributions and properties of the ICL and ICGCs out to a radius of 600 kpc (~1/3 of the virial radius) from the brightest cluster galaxy (BCG). We find that the central 500 kpc of the Perseus cluster hosts 70000$\pm$2800 GCs and $1.6\times10^{12}$ L$_\odot$ of diffuse light from the BCG+ICL in the near-infrared H$_E$. This accounts for 37$\pm$6% of the cluster's total stellar luminosity within this radius. The ICL and ICGCs share a coherent spatial distribution, suggesting a common origin or that a common potential governs their distribution. Their contours on the largest scales (>200 kpc) are offset from the BCG's core westwards by 60 kpc towards several luminous cluster galaxies. This offset is opposite to the displacement observed in the gaseous intracluster medium. The radial surface brightness profile of the BCG+ICL is best described by a double Sérsic model, with 68$\pm$4% of the H$_E$ light in the extended, outer component. The transition between these components occurs at ~50 kpc, beyond which the isophotes become increasingly elliptical and off-centred. The radial ICGC number density profile closely follows the BCG+ICL profile only beyond this 50 kpc radius, where we find an average of 60 GCs per $10^9$ M$_\odot$ of diffuse stellar mass. The BCG+ICL colour becomes increasingly blue with radius, consistent with the stellar populations in the ICL having subsolar metallicities [Fe/H]~-0.6. The colour of the ICL, and the specific frequency and luminosity function of the ICGCs suggest that the ICL+ICGCs were tidally stripped from the outskirts of massive satellites with masses of a few $\times10^{10}$ M$_\odot$, with an increasing contribution from dwarf galaxies at large radii. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations. 24 pages, 17 figures, 5 tables. Submitted to A&A

Euclid: Early Release Observations -- Dwarf galaxies in the Perseus galaxy cluster

Authors: F. R. Marleau, J. -C. Cuillandre, M. Cantiello, D. Carollo, P. -A. Duc, R. Habas, L. K. Hunt, P. Jablonka, M. Mirabile, M. Mondelin, M. Poulain, T. Saifollahi, R. Sánchez-Janssen, E. Sola, M. Urbano, R. Zöller, M. Bolzonella, A. Lançon, R. Laureijs, O. Marchal, M. Schirmer, C. Stone, A. Boselli, A. Ferré-Mateu, N. A. Hatch , et al. (171 additional authors not shown)

Abstract: We make use of the unprecedented depth, spatial resolution, and field of view of the Euclid Early Release Observations of the Perseus galaxy cluster to detect and characterise the dwarf galaxy population in this massive system. The Euclid high resolution VIS and combined VIS+NIR colour images were visually inspected and dwarf galaxy candidates were identified. Their morphologies, the presence of n… ▽ More We make use of the unprecedented depth, spatial resolution, and field of view of the Euclid Early Release Observations of the Perseus galaxy cluster to detect and characterise the dwarf galaxy population in this massive system. The Euclid high resolution VIS and combined VIS+NIR colour images were visually inspected and dwarf galaxy candidates were identified. Their morphologies, the presence of nuclei, and their globular cluster (GC) richness were visually assessed, complementing an automatic detection of the GC candidates. Structural and photometric parameters, including Euclid filter colours, were extracted from 2-dimensional fitting. Based on this analysis, a total of 1100 dwarf candidates were found across the image, with 638 appearing to be new identifications. The majority (96%) are classified as dwarf ellipticals, 53% are nucleated, 26% are GC-rich, and 6% show disturbed morphologies. A relatively high fraction of galaxies, 8%, are categorised as ultra-diffuse galaxies. The majority of the dwarfs follow the expected scaling relations. Globally, the GC specific frequency, S_N, of the Perseus dwarfs is intermediate between those measured in the Virgo and Coma clusters. While the dwarfs with the largest GC counts are found throughout the Euclid field of view, those located around the east-west strip, where most of the brightest cluster members are found, exhibit larger S_N values, on average. The spatial distribution of the dwarfs, GCs, and intracluster light show a main iso-density/isophotal centre displaced to the west of the bright galaxy light distribution. The ERO imaging of the Perseus cluster demonstrates the unique capability of Euclid to concurrently detect and characterise large samples of dwarfs, their nuclei, and their GC systems, allowing us to construct a detailed picture of the formation and evolution of galaxies over a wide range of mass scales and environments. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 44 pages, 24 figures, 5 tables, paper submitted to A&A as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions

Authors: J. -C. Cuillandre, M. Bolzonella, A. Boselli, F. R. Marleau, M. Mondelin, J. G. Sorce, C. Stone, F. Buitrago, Michele Cantiello, K. George, N. A. Hatch, L. Quilley, F. Mannucci, T. Saifollahi, R. Sánchez-Janssen, F. Tarsitano, C. Tortora, X. Xu, H. Bouy, S. Gwyn, M. Kluge, A. Lançon, R. Laureijs, M. Schirmer, Abdurro'uf , et al. (177 additional authors not shown)

Abstract: The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exception… ▽ More The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 35 figures, Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Globular clusters in the Fornax galaxy cluster, from dwarf galaxies to the intracluster field

Authors: T. Saifollahi, K. Voggel, A. Lançon, Michele Cantiello, M. A. Raj, J. -C. Cuillandre, S. S. Larsen, F. R. Marleau, A. Venhola, M. Schirmer, D. Carollo, P. -A. Duc, A. M. N. Ferguson, L. K. Hunt, M. Kümmel, R. Laureijs, O. Marchal, A. A. Nucita, R. F. Peletier, M. Poulain, M. Rejkuba, R. Sánchez-Janssen, M. Urbano, Abdurro'uf, B. Altieri , et al. (174 additional authors not shown)

Abstract: We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial… ▽ More We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial GCs injected into the data shows that Euclid's data in $I_{\rm E}$ band is 80% complete at about $I_{\rm E} \sim 26.0$ mag ($M_{V\rm } \sim -5.0$ mag), and resolves GCs as small as $r_{\rm h} = 2.5$ pc. In the $I_{\rm E}$ band, we detect more than 95% of the known GCs from previous spectroscopic surveys and GC candidates of the ACS Fornax Cluster Survey, of which more than 80% are resolved. We identify more than 5000 new GC candidates within the field of view down to $I_{\rm E}$ mag, about 1.5 mag fainter than the typical GC luminosity function turn-over magnitude, and investigate their spatial distribution within the intracluster field. We then focus on the GC candidates around dwarf galaxies and investigate their numbers, stacked luminosity distribution and stacked radial distribution. While the overall GC properties are consistent with those in the literature, an interesting over-representation of relatively bright candidates is found within a small number of relatively GC-rich dwarf galaxies. Our work confirms the capabilities of Euclid data in detecting GCs and separating them from foreground and background contaminants at a distance of 20 Mpc, particularly for low-GC count systems such as dwarf galaxies. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Deep anatomy of nearby galaxies

Authors: L. K. Hunt, F. Annibali, J. -C. Cuillandre, A. M. N. Ferguson, P. Jablonka, S. S. Larsen, F. R. Marleau, E. Schinnerer, M. Schirmer, C. Stone, C. Tortora, T. Saifollahi, A. Lançon, M. Bolzonella, S. Gwyn, M. Kluge, R. Laureijs, D. Carollo, M. L. M. Collins, P. Dimauro, P. -A. Duc, D. Erkal, J. M. Howell, C. Nally, E. Saremi , et al. (174 additional authors not shown)

Abstract: Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from… ▽ More Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from about 0.5 Mpc to 8.8 Mpc. Our assessment of the surface brightness depths in the stacked Euclid images confirms previous estimates in 100 arcsec^2 regions of 1sigma=30.5 mag/arcsec^2 for VIS, but slightly deeper than previous estimates for NISP with 1sigma=29.2-29.4 mag/arcsec^2. By combining Euclid HE, YE, and IE into RGB images, we illustrate the large field-of-view covered by a single Reference Observing Sequence, together with exquisite detail on parsec scales in these nearby galaxies. Radial surface brightness and color profiles demonstrate galaxy colors in agreement with stellar population synthesis models. Standard stellar photometry selection techniques find approximately 1.3 million stars across the 6 galaxy fields. Euclid's resolved stellar photometry allows us to constrain the star-formation histories of these galaxies, by disentangling the distributions of young stars, as well as asymptotic giant branch and red giant branch stellar populations. We finally examine 2 galaxies individually for surrounding satellite systems. Our analysis of the ensemble of dwarf satellites around NGC6744 reveals a new galaxy, EDwC1, a nucleated dwarf spheroidal at the end of a spiral arm. Our new census of the globular clusters around NGC2403 yields 9 new star-cluster candidates, 8 of which with colors indicative of evolved stellar populations. In summary, our investigation of the 6 Showcase galaxies demonstrates that Euclid is a powerful probe of the anatomy of nearby galaxies [abridged]. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 36 pages; 20 figures in main text; 4 Appendices. Submitted to A&A, as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry

Authors: J. -C. Cuillandre, E. Bertin, M. Bolzonella, H. Bouy, S. Gwyn, S. Isani, M. Kluge, O. Lai, A. Lançon, D. A. Lang, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, Abdurro'uf, N. Aghanim, B. Altieri, F. Annibali, H. Atek, P. Awad, M. Baes, E. Bañados, D. Barrado, S. Belladitta, V. Belokurov , et al. (240 additional authors not shown)

Abstract: The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline t… ▽ More The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Tidal disruption event AT2020ocn: early-time X-ray flares caused by a possible disc alignment process

Authors: Z. Cao, P. G. Jonker, D. R. Pasham, S. Wen, N. C. Stone, A. I. Zabludoff

Abstract: A tidal disruption event (TDE) may occur when a star is torn apart by the tidal force of a black hole (BH). Eventually, an accretion disc is thought to form out of stellar debris falling back towards the BH. If the star's orbital angular momentum vector prior to disruption is not aligned with the BH spin angular momentum vector, the disc will be tilted with respect to the BH equatorial plane. The… ▽ More A tidal disruption event (TDE) may occur when a star is torn apart by the tidal force of a black hole (BH). Eventually, an accretion disc is thought to form out of stellar debris falling back towards the BH. If the star's orbital angular momentum vector prior to disruption is not aligned with the BH spin angular momentum vector, the disc will be tilted with respect to the BH equatorial plane. The disc will eventually be drawn into the BH equatorial plane due to a combination of the Bardeen-Petterson effect and internal torques. Here, we analyse the X-ray and UV observations of the TDE AT2020ocn obtained by Swift, XMM-Newton, and NICER. The X-ray light curve shows strong flares during the first $\approx100$ days, while, over the same period, the UV emission decays gradually. We find that the X-ray flares can be explained by a model that also explains the spectral evolution. This model includes a slim disc viewed under a variable inclination plus an inverse-Comptonisation component processing the slim disc emission. A scenario where the ongoing Lense-Thirring precession during the disc alignment process is responsible for the observed inclination variations is consistent with the data. In later observations, we find that the X-ray spectrum of AT2020ocn becomes harder, while the mass accretion rate remains at super-Eddington levels, suggesting the formation of a corona in line with accretion onto other compact objects. We constrain the BH mass to be $(7^{+13}_{-3})\times10^{5}$ M$_\odot$ at the 1$σ$ (68%) confidence level. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures, 5 tables, with supplementary materials. Accepted for publication in ApJ

AT2018fyk: Candidate Tidal Disruption Event by a (Super)massive Black Hole Binary

Authors: S. Wen, P. G. Jonker, A. J. Levan, D. Li, N. C. Stone, A. I. Zabludoff, Z. Cao, T. Wevers, D. R. Pasham, C. Lewin, E. Kara

Abstract: The tidal disruption event (TDE) AT2018fyk has unusual X-ray, UV, and optical light curves that decay over the first $\sim$600d, rebrighten, and decay again around 1200d. We explain this behavior as a one-off TDE associated with a massive black hole (BH) \emph{binary}. The sharp drop-offs from $t^{-5/3}$ power laws at around 600d naturally arise when one BH interrupts the debris fallback onto the… ▽ More The tidal disruption event (TDE) AT2018fyk has unusual X-ray, UV, and optical light curves that decay over the first $\sim$600d, rebrighten, and decay again around 1200d. We explain this behavior as a one-off TDE associated with a massive black hole (BH) \emph{binary}. The sharp drop-offs from $t^{-5/3}$ power laws at around 600d naturally arise when one BH interrupts the debris fallback onto the other BH. The BH mass $M_\bullet$ derived from fitting X-ray spectra with a slim disk accretion model and, independently, from fitting the early UV/optical light curves, is smaller by two orders of magnitude than predicted from the $M_\bullet$--$σ_*$ host galaxy relation, suggesting that the debris is accreted onto the secondary, with fallback cut off by the primary. Furthermore, if the rebrightening were associated with the primary, it should occur around 5000d, not the observed 1200d. The secondary's mass and dimensionless spin is $M_{\bullet,{\rm s}}=2.7^{+0.5}_{-1.5} \times 10^5 M_\odot$ and $a_{\bullet,{\rm s}}>0.3$ (X-ray spectral fitting), while the primary's mass is $M_{\bullet,{\rm p}}=10^{7.7\pm0.4}M_\odot$ ($M_\bullet$-$σ_*$ relation). An intermediate mass BH secondary is consistent with the observed UV/optical light curve decay, i.e., the secondary's outer accretion disk is too faint to produce a detectable emission floor. The time of the first accretion cutoff constrains the binary separation to be $(6.7\pm 1.2) \times 10^{-3}~{\rm pc}$. X-ray spectral fitting and timing analysis indicate that the hard X-rays arise from a corona above the secondary's disk. The early UV/optical emission, suggesting a super-Eddington phase for the secondary, possibly originates from shocks arising from debris circularization. △ Less

Submitted 18 May, 2024; v1 submitted 1 May, 2024; originally announced May 2024.

Comments: 5 tables, 7 figures. Accepted for publication in APJ

Repeating partial disruptions and two-body relaxation

Authors: Luca Broggi, Nicholas C. Stone, Taeho Ryu, Elisa Bortolas, Massimo Dotti, Matteo Bonetti, Alberto Sesana

Abstract: Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only… ▽ More Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitation. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruptions that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates. △ Less

Submitted 19 June, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: Open Journal of Astrophysics, accepted for publication

The dark balance: quantifying the inner halo response to active galactic nuclei feedback in galaxies

Authors: Nikhil Arora, Stephane Courteau, Andrea V. Maccio, Changhyun Cho, Raj Patel, Connor Stone

Abstract: This paper presents a study of the impact of supermassive black hole (SMBH) feedback on dark matter (DM) halos in numerical NIHAO simulations of galaxies. In particular, the amount of DM displaced via active galactic nuclei (AGN) feedback and the physical scale over which AGN feedback affects the DM halo are quantified by comparing NIHAO simulations with and without AGN feedback. NIHAO galaxies wi… ▽ More This paper presents a study of the impact of supermassive black hole (SMBH) feedback on dark matter (DM) halos in numerical NIHAO simulations of galaxies. In particular, the amount of DM displaced via active galactic nuclei (AGN) feedback and the physical scale over which AGN feedback affects the DM halo are quantified by comparing NIHAO simulations with and without AGN feedback. NIHAO galaxies with $\log(M_*/M_{\rm \odot})\geq 10.0$ show a growing central DM suppression of 0.2 dex (~40%) from z = 1.5 to the present relative to noAGN feedback simulations. The growth of the DM suppression is related to the mass evolution of the SMBH and the gas mass in the central regions. For the most massive NIHAO galaxies with $\log(M_*/M_{\rm \odot}) > 10.5$, partially affected by numerical resolution, the central DM suppression peaks at z = 0.5, after which halo contraction overpowers AGN feedback due a shortage of gas and, thus, SMBH growth. The spatial scale, or ``sphere of influence,'' over which AGN feedback affects the DM distribution decreases as a function of time for MW-mass galaxies (from ~16 kpc at z = 1.5 to ~7.8 kpc at z = 0) as a result of halo contraction due to stellar growth. For the most massive NIHAO galaxies, the size of the sphere of influence remains constant (~16 kpc) for z > 0.5 owing to the balance between AGN feedback and halo contraction. △ Less

Submitted 6 March, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

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

XMM-Newton-discovered Fast X-ray Transients: Host galaxies and limits on contemporaneous detections of optical counterparts

Authors: D. Eappachen, P. G. Jonker, J. Quirola-Vásquez, D. Mata Sánchez, A. Inkenhaag, A. J. Levan, M. Fraser, M. A. P. Torres, F. E. Bauer, A. A. Chrimes, D. Stern, M. J. Graham, S. J. Smartt, K. W. Smith, M. E. Ravasio, A. I. Zabludoff, M. Yue, F. Stoppa, D. B. Malesani, N. C. Stone, S. Wen

Abstract: Extragalactic fast X-ray transients (FXTs) are a class of soft (0.3-10 keV) X-ray transients lasting a few hundred seconds to several hours. Several progenitor mechanisms have been suggested to produce FXTs, including supernova shock breakouts, binary neutron star mergers, or tidal disruptions involving an intermediate-mass black hole and a white dwarf. We present detailed host studies, including… ▽ More Extragalactic fast X-ray transients (FXTs) are a class of soft (0.3-10 keV) X-ray transients lasting a few hundred seconds to several hours. Several progenitor mechanisms have been suggested to produce FXTs, including supernova shock breakouts, binary neutron star mergers, or tidal disruptions involving an intermediate-mass black hole and a white dwarf. We present detailed host studies, including spectroscopic observations of the host galaxies of 7 XMM-Newton-discovered FXTs. The candidate hosts lie at redshifts 0.0928 $< z <$ 0.645 implying peak X-ray luminosities of 10$^{43}$ erg s$^{-1}$ $< L_X <$ 10$^{45}$ erg s$^{-1}$,and physical offsets of 1 kpc < $r_\mathrm{proj}$ < 22 kpc. These observations increase the number of FXTs with a spectroscopic redshift measurement by a factor of 2, although we note that one event is re-identified as a Galactic flare star. We infer host star formation rates and stellar masses by fitting the combined spectroscopic and archival photometric data. We also report on a contemporaneous optical counterpart search to the FXTs in Pan-STARRS and ATLAS by performing forced photometry at the position of the FXTs. We do not find any counterpart in our search. Given our constraints, including peak X-ray luminosities, optical limits, and host properties, we find that XRT 110621 is consistent with a SN SBO event. Spectroscopic redshifts of likely host galaxies for four events imply peak X-ray luminosities that are too high to be consistent with SN SBOs, but we are unable to discard either the BNS or WD-IMBH TDE scenarios for these FXTs. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: 17 pages, 10 figures, 6 tables. Accepted for publication in MNRAS

Echoes in the Noise: Posterior Samples of Faint Galaxy Surface Brightness Profiles with Score-Based Likelihoods and Priors

Authors: Alexandre Adam, Connor Stone, Connor Bottrell, Ronan Legin, Yashar Hezaveh, Laurence Perreault-Levasseur

Abstract: Examining the detailed structure of galaxy populations provides valuable insights into their formation and evolution mechanisms. Significant barriers to such analysis are the non-trivial noise properties of real astronomical images and the point spread function (PSF) which blurs structure. Here we present a framework which combines recent advances in score-based likelihood characterization and dif… ▽ More Examining the detailed structure of galaxy populations provides valuable insights into their formation and evolution mechanisms. Significant barriers to such analysis are the non-trivial noise properties of real astronomical images and the point spread function (PSF) which blurs structure. Here we present a framework which combines recent advances in score-based likelihood characterization and diffusion model priors to perform a Bayesian analysis of image deconvolution. The method, when applied to minimally processed \emph{Hubble Space Telescope} (\emph{HST}) data, recovers structures which have otherwise only become visible in next-generation \emph{James Webb Space Telescope} (\emph{JWST}) imaging. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 5+5 pages, 10 figures, Machine Learning and the Physical Sciences Workshop, NeurIPS 2023

Towards Exploratory Reformulation of Constraint Models

Authors: Ian Miguel, András Z. Salamon, Christopher Stone

Abstract: It is well established that formulating an effective constraint model of a problem of interest is crucial to the efficiency with which it can subsequently be solved. Following from the observation that it is difficult, if not impossible, to know a priori which of a set of candidate models will perform best in practice, we envisage a system that explores the space of models through a process of ref… ▽ More It is well established that formulating an effective constraint model of a problem of interest is crucial to the efficiency with which it can subsequently be solved. Following from the observation that it is difficult, if not impossible, to know a priori which of a set of candidate models will perform best in practice, we envisage a system that explores the space of models through a process of reformulation from an initial model, guided by performance on a set of training instances from the problem class under consideration. We plan to situate this system in a refinement-based approach, where a user writes a constraint specification describing a problem above the level of abstraction at which many modelling decisions are made. In this position paper we set out our plan for an exploratory reformulation system, and discuss progress made so far. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 13 pages, 6 figures

ACM Class: F.4.2

Exploring Benchmarks for Self-Driving Labs using Color Matching

Authors: Tobias Ginsburg, Kyle Hippe, Ryan Lewis, Doga Ozgulbas, Aileen Cleary, Rory Butler, Casey Stone, Abraham Stroka, Ian Foster

Abstract: Self Driving Labs (SDLs) that combine automation of experimental procedures with autonomous decision making are gaining popularity as a means of increasing the throughput of scientific workflows. The task of identifying quantities of supplied colored pigments that match a target color, the color matching problem, provides a simple and flexible SDL test case, as it requires experiment proposal, sam… ▽ More Self Driving Labs (SDLs) that combine automation of experimental procedures with autonomous decision making are gaining popularity as a means of increasing the throughput of scientific workflows. The task of identifying quantities of supplied colored pigments that match a target color, the color matching problem, provides a simple and flexible SDL test case, as it requires experiment proposal, sample creation, and sample analysis, three common components in autonomous discovery applications. We present a robotic solution to the color matching problem that allows for fully autonomous execution of a color matching protocol. Our solution leverages the WEI science factory platform to enable portability across different robotic hardware, the use of alternative optimization methods for continuous refinement, and automated publication of results for experiment tracking and post-hoc analysis. △ Less

Submitted 30 September, 2023; originally announced October 2023.

A Large Language Model Approach to Educational Survey Feedback Analysis

Authors: Michael J. Parker, Caitlin Anderson, Claire Stone, YeaRim Oh

Abstract: This paper assesses the potential for the large language models (LLMs) GPT-4 and GPT-3.5 to aid in deriving insight from education feedback surveys. Exploration of LLM use cases in education has focused on teaching and learning, with less exploration of capabilities in education feedback analysis. Survey analysis in education involves goals such as finding gaps in curricula or evaluating teachers,… ▽ More This paper assesses the potential for the large language models (LLMs) GPT-4 and GPT-3.5 to aid in deriving insight from education feedback surveys. Exploration of LLM use cases in education has focused on teaching and learning, with less exploration of capabilities in education feedback analysis. Survey analysis in education involves goals such as finding gaps in curricula or evaluating teachers, often requiring time-consuming manual processing of textual responses. LLMs have the potential to provide a flexible means of achieving these goals without specialized machine learning models or fine-tuning. We demonstrate a versatile approach to such goals by treating them as sequences of natural language processing (NLP) tasks including classification (multi-label, multi-class, and binary), extraction, thematic analysis, and sentiment analysis, each performed by LLM. We apply these workflows to a real-world dataset of 2500 end-of-course survey comments from biomedical science courses, and evaluate a zero-shot approach (i.e., requiring no examples or labeled training data) across all tasks, reflecting education settings, where labeled data is often scarce. By applying effective prompting practices, we achieve human-level performance on multiple tasks with GPT-4, enabling workflows necessary to achieve typical goals. We also show the potential of inspecting LLMs' chain-of-thought (CoT) reasoning for providing insight that may foster confidence in practice. Moreover, this study features development of a versatile set of classification categories, suitable for various course types (online, hybrid, or in-person) and amenable to customization. Our results suggest that LLMs can be used to derive a range of insights from survey text. △ Less

Submitted 26 June, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Journal ref: Int J Artif Intell Educ (2024)

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.

Comments: 10 pages, 8 figures. Comments are welcome

Towards a Modular Architecture for Science Factories

Authors: Rafael Vescovi, Tobias Ginsburg, Kyle Hippe, Doga Ozgulbas, Casey Stone, Abraham Stroka, Rory Butler, Ben Blaiszik, Tom Brettin, Kyle Chard, Mark Hereld, Arvind Ramanathan, Rick Stevens, Aikaterini Vriza, Jie Xu, Qingteng Zhang, Ian Foster

Abstract: Advances in robotic automation, high-performance computing (HPC), and artificial intelligence (AI) encourage us to conceive of science factories: large, general-purpose computation- and AI-enabled self-driving laboratories (SDLs) with the generality and scale needed both to tackle large discovery problems and to support thousands of scientists. Science factories require modular hardware and softwa… ▽ More Advances in robotic automation, high-performance computing (HPC), and artificial intelligence (AI) encourage us to conceive of science factories: large, general-purpose computation- and AI-enabled self-driving laboratories (SDLs) with the generality and scale needed both to tackle large discovery problems and to support thousands of scientists. Science factories require modular hardware and software that can be replicated for scale and (re)configured to support many applications. To this end, we propose a prototype modular science factory architecture in which reconfigurable modules encapsulating scientific instruments are linked with manipulators to form workcells, that can themselves be combined to form larger assemblages, and linked with distributed computing for simulation, AI model training and inference, and related tasks. Workflows that perform sets of actions on modules can be specified, and various applications, comprising workflows plus associated computational and data manipulation steps, can be run concurrently. We report on our experiences prototyping this architecture and applying it in experiments involving 15 different robotic apparatus, five applications (one in education, two in biology, two in materials), and a variety of workflows, across four laboratories. We describe the reuse of modules, workcells, and workflows in different applications, the migration of applications between workcells, and the use of digital twins, and suggest directions for future work aimed at yet more generality and scalability. Code and data are available at https://ad-sdl.github.io/wei2023 and in the Supplementary Information △ Less

Submitted 17 October, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

AstroPhot: Fitting Everything Everywhere All at Once in Astronomical Images

Authors: Connor Stone, Stephane Courteau, Jean-Charles Cuillandre, Yashar Hezaveh, Laurence Perreault-Levasseur, Nikhil Arora

Abstract: We present AstroPhot, a fast, powerful, and user-friendly Python based astronomical image photometry solver. AstroPhot incorporates automatic differentiation and GPU (or parallel CPU) acceleration, powered by the machine learning library PyTorch. Everything: AstroPhot can fit models for sky, stars, galaxies, PSFs, and more in a principled Chi^2 forward optimization, recovering Bayesian posterior i… ▽ More We present AstroPhot, a fast, powerful, and user-friendly Python based astronomical image photometry solver. AstroPhot incorporates automatic differentiation and GPU (or parallel CPU) acceleration, powered by the machine learning library PyTorch. Everything: AstroPhot can fit models for sky, stars, galaxies, PSFs, and more in a principled Chi^2 forward optimization, recovering Bayesian posterior information and covariance of all parameters. Everywhere: AstroPhot can optimize forward models on CPU or GPU; across images that are large, multi-band, multi-epoch, rotated, dithered, and more. All at once: The models are optimized together, thus handling overlapping objects and including the covariance between parameters (including PSF and galaxy parameters). A number of optimization algorithms are available including Levenberg-Marquardt, Gradient descent, and No-U-Turn MCMC sampling. With an object-oriented user interface, AstroPhot makes it easy to quickly extract detailed information from complex astronomical data for individual images or large survey programs. This paper outlines novel features of the AstroPhot code and compares it to other popular astronomical image modeling software. AstroPhot is open-source, fully Python based, and freely accessible here: https://github.com/Autostronomy/AstroPhot △ Less

Submitted 6 September, 2023; v1 submitted 3 August, 2023; originally announced August 2023.

Comments: 18 Pages, 9 figures, published in MNRAS

The Effect of Thermal Torques on AGN Disc Migration Traps and Gravitational Wave Populations

Authors: Evgeni Grishin, Shmuel Gilbaum, Nicholas C. Stone

Abstract: Accretion discs in active galactic nuclei (AGN) foster black hole (BH) formation, growth, and mergers. Stellar mass BHs migrate inwards under the influence of hydrodynamical torques unless they encounter a region where the torque flips sign. At these migration traps, BHs accumulate and merge via dynamical or gas-assisted interactions, producing high-frequency LIGO/Virgo/KAGRA (LVK) gravitational w… ▽ More Accretion discs in active galactic nuclei (AGN) foster black hole (BH) formation, growth, and mergers. Stellar mass BHs migrate inwards under the influence of hydrodynamical torques unless they encounter a region where the torque flips sign. At these migration traps, BHs accumulate and merge via dynamical or gas-assisted interactions, producing high-frequency LIGO/Virgo/KAGRA (LVK) gravitational wave (GW) sources and potentially cutting off the supply of extreme mass ratio inspirals that would otherwise make low-frequency, {\it LISA}-band GWs. In this paper, we study the interplay between different types of migration torques, focusing especially on the ``thermal torques'' generated by the thermal response of the AGN to embedded stellar-mass BHs that accrete through their own mini-discs.In contrast to previous work, we find that Type I torques cannot produce migration traps on their own, but thermal torques often do, particularly in low-mass AGN. The migration traps produced by thermal torques exist at much larger radii ($\sim 10^{3-5}$ gravitational radii) than do previously identified Type I traps, carrying implications for GW populations at multiple frequencies. Finally, we identify a bifurcation of AGN discs into two regimes: migration traps exist below a critical AGN luminosity, and do not at higher luminosities. This critical luminosity is fit as $\log_{10} L_{\rm AGN}^c = 45 - 0.32 \log_{10}{(α/0.01)}$ where $α$ is the AGN alpha viscosity parameter, a range compatible with recent claims that LVK GWs are not preferentially associated with high-luminosity AGN. △ Less

Submitted 19 March, 2024; v1 submitted 14 July, 2023; originally announced July 2023.

Comments: Accepted to MNRAS

Heating Galaxy Clusters with Interacting Dark Matter

Authors: Yutaro Shoji, Eric Kuflik, Yuval Birnboim, Nicholas C. Stone

Abstract: The overcooling of cool core clusters is a persistent puzzle in the astrophysics of galaxy clusters. We propose that it may naturally be resolved via interactions between the baryons of the intracluster medium (ICM) and its dark matter (DM). DM-baryon interactions can inject heat into the ICM to offset bremmstrahlung cooling, but these interactions are also strongly constrained by existing experim… ▽ More The overcooling of cool core clusters is a persistent puzzle in the astrophysics of galaxy clusters. We propose that it may naturally be resolved via interactions between the baryons of the intracluster medium (ICM) and its dark matter (DM). DM-baryon interactions can inject heat into the ICM to offset bremmstrahlung cooling, but these interactions are also strongly constrained by existing experiments and astrophysical observations. We survey existing constraints and combine these with the energetic needs of an observed sample of cool core clusters. We find that a robust parameter space exists for baryon-DM scattering solutions to the cooling flow problem, provided that only a sub-component of DM interacts strongly with the baryons. Interestingly, baryon-DM scattering is a thermally stable heating source so long as the baryon temperature is greater than $1/3-1/2$ the DM temperature, a condition that seems to be satisfied observationally. △ Less

Submitted 4 March, 2024; v1 submitted 14 June, 2023; originally announced June 2023.

Comments: 11 pages, 2 figures; corrected figure 1

Enhanced Extreme Mass Ratio Inspiral Rates into Intermediate Mass Black Holes

Authors: Ismail Qunbar, Nicholas C. Stone

Abstract: Extreme mass ratio inspirals (EMRIs) occur when stellar-mass compact objects begin a gravitational wave (GW) driven inspiral into massive black holes. EMRI waveforms can precisely map the surrounding spacetime, making them a key target for future space-based GW interferometers such as {\it LISA}, but their event rates and parameters are massively uncertain. One of the largest uncertainties is the… ▽ More Extreme mass ratio inspirals (EMRIs) occur when stellar-mass compact objects begin a gravitational wave (GW) driven inspiral into massive black holes. EMRI waveforms can precisely map the surrounding spacetime, making them a key target for future space-based GW interferometers such as {\it LISA}, but their event rates and parameters are massively uncertain. One of the largest uncertainties is the ratio of true EMRIs (which spend at least thousands of orbits in the {\it LISA} band) and direct plunges, which are in-band for at most a handful of orbits and are not detectable in practice. In this paper, we show that the traditional dichotomy between EMRIs and plunges -- EMRIs originate from small semimajor axes, plunges from large -- does not hold for intermediate-mass black holes with masses $M_\bullet \lesssim 10^5 M_\odot$. In this low-mass regime, a plunge always has an $\mathcal{O}(1)$ probability of failing and transitioning into a novel ``cliffhanger'' EMRI. Cliffhanger EMRIs are more easily produced for larger stellar-mass compact objects, and are less likely for smaller ones. This new EMRI production channel can dominate volumetric EMRI rates $\dot{n}_{\rm EMRI}$ if intermediate-mass black holes are common in dwarf galactic nuclei, potentially increasing $\dot{n}_{\rm EMRI}$ by an order of magnitude. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 8 pages, 6 figures, comments welcome

MaNGA galaxy properties -- II. A detailed comparison of observed and simulated spiral galaxy scaling relations

Authors: Nikhil Arora, Stéphane Courteau, Connor Stone, Andrea V. Macció

Abstract: We present a catalogue of dynamical properties for 2368 late-type galaxies from the MaNGA survey. The latter complements the catalogue of photometric properties for the same sample based on deep optical DESI photometry processed with AutoProf. Rotation curves (RCs), extracted by model fitting H$α$ velocity maps from the MaNGA Data Analysis Pipeline, extend out to 1.4 (1.9) R$_{e}$ for the primary… ▽ More We present a catalogue of dynamical properties for 2368 late-type galaxies from the MaNGA survey. The latter complements the catalogue of photometric properties for the same sample based on deep optical DESI photometry processed with AutoProf. Rotation curves (RCs), extracted by model fitting H$α$ velocity maps from the MaNGA Data Analysis Pipeline, extend out to 1.4 (1.9) R$_{e}$ for the primary (secondary) MaNGA samples. The RCs and ancillary MaNGA Pipe3D data products were used to construct various fundamental galaxy scaling relations that are also compared uniformly with similar relations from NIHAO zoom-in simulations. Simulated NIHAO galaxies were found to broadly reproduce the observed MaNGA galaxy population for $\log (M_*/{\rm M_{\odot}) > 8.5}$. Some discrepancies remain, such as those pertaining to central stellar densities and the diversity of RCs due to strong feedback schemes. Also presented are spatially-resolved scatters for the velocity-size-stellar mass (VRM$_*$) structural relations using MaNGA and NIHAO samples. The scatter for these relations in the galaxian interiors is a consequence of the diversity of inner RC shapes, while scatter in the outskirts is dictated by the large range of stellar surface densities which itself is driven by sporadic star formation. The detailed spatially-resolved scatter analysis highlights the complex interplay between local and global astrophysical processes and provides a strong constraint to numerical simulations. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 21 pages, 13 Figures, Accepted for publication in MNRAS

Optical/UV Emission in the Tidal Disruption Event ASASSN-14li: Implications of Disc Modeling

Authors: Sixiang Wen, Peter G. Jonker, Nicholas C. Stone, Sjoert Van Velzen, Ann I. Zabludoff

Abstract: We predict late-time optical/UV emission from tidal disruption events (TDEs) from our slim accretion disc model \citep{Wen20} and explore the impact of the black hole mass $M_\bullet$, black hole spin $a_\bullet$, and accretion disc size. We use these synthetic spectra to successfully fit the multi-band \emph{Swift} observations of ASASSN-14li at >350 days, setting only the host galaxy extinction… ▽ More We predict late-time optical/UV emission from tidal disruption events (TDEs) from our slim accretion disc model \citep{Wen20} and explore the impact of the black hole mass $M_\bullet$, black hole spin $a_\bullet$, and accretion disc size. We use these synthetic spectra to successfully fit the multi-band \emph{Swift} observations of ASASSN-14li at >350 days, setting only the host galaxy extinction and outer disc radius as free parameters and employing the $M_\bullet$, $a_\bullet$, disc inclination, and disc accretion rates derived from fitting 10 epochs of ASASSN-14li's X-ray spectra with the slim disc. To address the nature of the \emph{early}-time optical/UV emission, we consider two models: shock dissipation and reprocessing. We find that (1) the predicted late-time optical/UV colour (e.g., $u-w2$) is insensitive to black hole and disc parameters unless the disc spreads quickly; (2) a starburst galaxy extinction model is required to fit the data, consistent with ASASSN-14li's post-starburst host; (3) surprisingly, the outer disc radius is $\approx$2$\times$ the tidal radius and $\sim$constant at late times, showing that viscous spreading is slow or non-existent; (4) the shock model can be self-consistent if $M_\bullet \lesssim 10^{6.75}$M$_\odot$, i.e., on the low end of ASASSN-14li's $M_\bullet$ range ($10^{6.5-7.1}$M$_\odot$; 1$σ$ CL); larger black hole masses require disruption of an unrealistically massive progenitor star; (5) the gas mass needed for reprocessing, whether by a quasi-static or an outflowing layer, can be $<0.5$M$_\odot$, consistent with a (plausible) disruption of a solar-mass star. △ Less

Submitted 4 April, 2023; v1 submitted 1 April, 2023; originally announced April 2023.

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

Precision Measurement of the Specific Activity of $^{39}$Ar in Atmospheric Argon with the DEAP-3600 Detector

Authors: P. Adhikari, R. Ajaj, M. Alpízar-Venegas, P. -A. Amaudruz, J. Anstey, G. R. Araujo, D. J. Auty, M. Baldwin, M. Batygov, B. Beltran, H. Benmansour, C. E. Bina, J. Bonatt, W. Bonivento, M. G. Boulay, B. Broerman, J. F. Bueno, P. M. Burghardt, A. Butcher, M. Cadeddu, B. Cai, M. Cárdenas-Montes, S. Cavuoti, M. Chen, Y. Chen , et al. (125 additional authors not shown)

Abstract: The specific activity of the beta decay of $^{39}$Ar in atmospheric argon is measured using the DEAP-3600 detector. DEAP-3600, located 2 km underground at SNOLAB, uses a total of (3269 $\pm$ 24) kg of liquid argon distilled from the atmosphere to search for dark matter. This detector with very low background uses pulseshape discrimination to differentiate between nuclear recoils and electron recoi… ▽ More The specific activity of the beta decay of $^{39}$Ar in atmospheric argon is measured using the DEAP-3600 detector. DEAP-3600, located 2 km underground at SNOLAB, uses a total of (3269 $\pm$ 24) kg of liquid argon distilled from the atmosphere to search for dark matter. This detector with very low background uses pulseshape discrimination to differentiate between nuclear recoils and electron recoils and is well-suited to measure the decay of $^{39}$Ar. With 167 live-days of data, the measured specific activity at the time of atmospheric extraction is [0.964 $\pm$ 0.001 (stat) $\pm$ 0.024 (sys)] Bq/kg$_{\rm atmAr}$ which is consistent with results from other experiments. A cross-check analysis using different event selection criteria provides a consistent result. △ Less

Submitted 10 October, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Journal ref: Eur. Phys. J. C 83, 642 (2023)

The Rapidly Spinning Intermediate-Mass Black Hole 3XMM J150052.0+015452

Authors: Z. Cao, P. G. Jonker, S. Wen, N. C. Stone, A. I. Zabludoff

Abstract: A star tidally disrupted by a black hole can form an accretion disc with a super-Eddington mass accretion rate; the X-ray emission produced by the inner disc provides constraints on the black hole mass $M_\bullet$ and dimensionless spin parameter $a_\bullet$. Previous studies have suggested that the $M_\bullet$ responsible for the tidal disruption event 3XMM J150052.0+015452 (hereafter J150052) is… ▽ More A star tidally disrupted by a black hole can form an accretion disc with a super-Eddington mass accretion rate; the X-ray emission produced by the inner disc provides constraints on the black hole mass $M_\bullet$ and dimensionless spin parameter $a_\bullet$. Previous studies have suggested that the $M_\bullet$ responsible for the tidal disruption event 3XMM J150052.0+015452 (hereafter J150052) is $\sim$10$^{5} M_{\odot}$, in the intermediate black hole (IMBH) regime. Fitting multi-epoch XMM-Newton and Chandra X-ray spectra obtained after 2008 during the source's decade-long decay, with our latest slim accretion disc model gives $M_\bullet = 2.0^{+1.0}_{-0.3}\times10^{5} M_{\odot}$ (at 68% confidence) and $a_\bullet > 0.97$ (a 84.1% confidence lower limit). The spectra obtained between 2008-2014 are significantly harder than those after 2014, an evolution that can be well explained by including the effects of inverse-Comptonisation by a corona on the early-time spectra. The corona is present when the source accretion rate is super-Eddington, while there is no evidence for its effect in data obtained after 2014, when the mass accretion rate is around the Eddington-limit. Based on our spectral study, we infer that the corona is optically thick and warm ($kT_e=2.3^{+2.7}_{-0.8}$ keV). Our mass and spin measurements of J150052 confirm it as an IMBH and point to a rapid, near extremal, spin. These $M_\bullet$ and $a_\bullet$ values rule out both vector bosons and axions of masses $\sim10^{-16}$ eV. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 11 pages, 6 figures. Accepted for publication in MNRAS

The growth of intermediate mass black holes through tidal captures and tidal disruption events

Authors: Francesco Paolo Rizzuto, Thorsten Naab, Antti Rantala, Peter H. Johansson, Jeremiah P. Ostriker, Nicholas C. Stone, Shihong Liao, Dimitrios Irodotou

Abstract: We present $N\mathrm{-body} $ simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 $\mathrm{M_{\odot}}$. The models are evolved with the $N\mathrm{-body} $ code \textsc{bifrost} to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruptio… ▽ More We present $N\mathrm{-body} $ simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 $\mathrm{M_{\odot}}$. The models are evolved with the $N\mathrm{-body} $ code \textsc{bifrost} to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruption events (TDEs). We model star-BH tidal interactions using a velocity-dependent drag force, which causes orbital energy and angular momentum loss near the BH. About $\sim 20-30$ per cent of the stars within the spheres of influence of the black holes form Bahcall-Wolf cusps and prevent the systems from core collapse. Within the first 40 Myr of evolution, the systems experience 500 up to 1300 TDEs, depending on the initial cluster structure. Most ($> 95$ per cent) of the TDEs originate from stars in the Bahcall-Wolf cusp. We derive an analytical formula for the TDE rate as a function of the central BH mass, density and velocity dispersion of the clusters ($\dot{N}_{\mathrm{TDE}} \propto M\mathrm{_{BH}} ρσ^{-3}$). We find that TDEs can lead a 300 $\mathrm{M_{\odot}}$ BH to reach $\sim 7000 \mathrm{M_{\odot}}$ within a Gyr. This indicates that TDEs can drive the formation and growth of massive BHs in sufficiently dense environments, which might be present in the central regions of nuclear star clusters. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 17 pages, 18 figures

Loss Cone Shielding

Authors: Odelia Teboul, Nicholas C. Stone, Jeremiah P. Ostriker

Abstract: A star wandering close enough to a massive black hole (MBH) can be ripped apart by the tidal forces of the black hole. The advent of wide-field surveys at many wavelengths has quickly increased the number of tidal disruption events (TDEs) observed, and has revealed that i) observed TDE rates are lower than theoretical predictions and ii) E+A galaxies are significantly overrepresented. This overrep… ▽ More A star wandering close enough to a massive black hole (MBH) can be ripped apart by the tidal forces of the black hole. The advent of wide-field surveys at many wavelengths has quickly increased the number of tidal disruption events (TDEs) observed, and has revealed that i) observed TDE rates are lower than theoretical predictions and ii) E+A galaxies are significantly overrepresented. This overrepresentation further worsens the tension between observed and theoretically predicted TDEs for non-E+A galaxies. Classical loss cone theory focuses on the cumulative effect of many weak scatterings. However, a strong scattering can remove a star from the distribution before it can get tidally disrupted. Most stars undergoing TDEs come from within the radius of influence, the densest environments of the universe. In such environments, close encounters rare elsewhere become non-negligible. We revise the standard loss cone theory to take into account classical two-body interactions as well as strong scattering, collisions, tidal captures, and study under which conditions close encounters can shield the loss cone. We i) analytically derive the impact of strong scattering and other close encounters, ii) compute time-dependent loss cone dynamics including both weak and strong encounters, and iii) derive analytical solutions to the Fokker-Planck equation with strong scattering. We find that i) TDE rates can be reduced to up to an order of magnitude and ii) strong shielding preferentially reduces deeply plunging stars. We also show that stellar overdensities, one possible explanation for the E+A preference, can fail to increase TDE rates when taking into account strong scattering. △ Less

Submitted 26 October, 2023; v1 submitted 10 November, 2022; originally announced November 2022.

Comments: 14 pages, 10 figures, Accepted for publication MNRAS

Magnetically Dominated Disks in Tidal Disruption Events and Quasi-Periodic Eruptions

Authors: Karamveer Kaur, Nicholas C. Stone, Shmuel Gilbaum

Abstract: The classical radiation pressure instability has been a persistent theoretical feature of thin, radiatively efficient accretion disks with accretion rates 1 to 100 per cent of the Eddington rate. But there is only limited evidence of its occurrence in nature: rapid heartbeat oscillations of a few X-ray binaries and now, perhaps, the new class of hourly X-ray transients called quasi-periodic erupti… ▽ More The classical radiation pressure instability has been a persistent theoretical feature of thin, radiatively efficient accretion disks with accretion rates 1 to 100 per cent of the Eddington rate. But there is only limited evidence of its occurrence in nature: rapid heartbeat oscillations of a few X-ray binaries and now, perhaps, the new class of hourly X-ray transients called quasi-periodic eruptions (QPEs). The accretion disks formed in tidal disruption events (TDEs) have been observed to peacefully trespass through the range of unstable accretion rates without exhibiting any clear sign of the instability. We try to explain the occurrence or otherwise of this instability in these systems, by constructing steady state 1D models of thin magnetic accretion disks. The local magnetic pressure in the disk is assumed to be dominated by toroidal fields arising from a dynamo sourced by magneto-rotational instability (MRI). We choose a physically motivated criterion of MRI saturation, validated by recent magnetohydrodynamic simulations, to determine the strength of magnetic pressure in the disk. The resulting magnetic pressure support efficiently shrinks: (1) the parameter space of unstable mass accretion rates, explaining the absence of instability in systems such as TDEs and (2) the range of unstable radii in the inner accretion disk, which can shorten the quasi-periods of instability limit-cycles by more than three orders of magnitude, explaining the observed periods ( a few hrs) of QPEs. In addition to examining stability properties of strongly magnetized disks, we predict other observational signatures such as spectral hardening factors and jet luminosities to test the compatibility of our disk models with observations of apparently stable TDE disks. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: Submitted to MNRAS

PROBES-I: A Compendium of Deep Rotation Curves and Matched multiband Photometry

Authors: Connor Stone, Stephane Courteau, Nikhil Arora, Matthew Frosst, Thomas Jarrett

Abstract: We present the Photometry and Rotation Curve Observations from Extragalactic Surveys (PROBES) compendium of extended rotation curves for 3163 late-type spirals, with matching homogeneous multiband photometry for 1677 of them. PROBES rotation curves originally extracted from Halpha long-slit spectra and aperture synthesis HI (21cm) velocity maps typically extend out to a median 2R_e (or 1R_{23.5, r… ▽ More We present the Photometry and Rotation Curve Observations from Extragalactic Surveys (PROBES) compendium of extended rotation curves for 3163 late-type spirals, with matching homogeneous multiband photometry for 1677 of them. PROBES rotation curves originally extracted from Halpha long-slit spectra and aperture synthesis HI (21cm) velocity maps typically extend out to a median 2R_e (or 1R_{23.5, r}). Our uniform photometry takes advantage of GALEX, DESI-LIS, and WISE images and the software AutoProf to yield multiband azimuthally averaged surface brightness profiles that achieve depths greater than 25 mag/arcsec^2 (FUV, NUV), 27 mag/arcsec^2 (g, r), and 26 mag/arcsec^2 (z, W1, W2). With its library of spatially resolved profiles and an extensive table of structural parameters, the versatile PROBES data set will benefit studies of galaxy structure and formation. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 18 pages, 12 figures, Published in ApJS

Journal ref: ApJS 262 33 (2022)

A Tale of a Tail: A Tidally-Disrupting Ultra-Diffuse Galaxy in the M81 Group

Authors: Rokas Žemaitis, Annette M. N. Ferguson, Sakurako Okamoto, Jean-Charles Cuillandre, Connor J. Stone, Nobuo Arimoto, Mike J. Irwin

Abstract: We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper… ▽ More We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper Suprime-Cam survey of the M81 Group. It has an average surface brightness of $μ_g \sim 32$ mag arcsec$^{-2}$ and can be traced for over a degree on the sky (60 kpc at the distance of F8D1) with our current imagery. We revisit the main body properties of F8D1 using deep multiband imagery acquired with MegaCam on CFHT and measure effective radii of $1.7-1.9$ kpc, central surface brightnesses of $24.7-25.7$ mag and a stellar mass of $\sim7 \times 10^7 M_{\odot}$. Assuming a symmetric feature on the other side of the galaxy, we calculate that $30-36$% of F8D1's present-day luminosity is contained in the tail. We argue that the most likely origin of F8D1's disruption is a recent close passage to M81, which would have stripped its gas and quenched its star formation. As the only UDG that has so far been studied to such faint surface brightness depths, the unveiling of F8D1's tidal disruption is important. It leaves open the possibility that many other UDGs could be the result of similar processes, with the most telling signatures of this lurking below current detection limits. △ Less

Submitted 27 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

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

The scientific payload of the Ultraviolet Transient Astronomy Satellite (ULTRASAT)

Authors: Sagi Ben-Ami, Yossi Shvartzvald, Eli Waxman, Udi Netzer, Yoram Yaniv, Viktor M. Algranatti, Avishay Gal-Yam, Ofer Lapid, Eran Ofek, Jeremy Topaz, Iair Arcavi, Arooj Asif, Shlomi Azaria, Eran Bahalul, Merlin F. Barschke, Benjamin Bastian-Querner, David Berge, Vlad D. Berlea, Rolf Buhler, Louise Dittmar, Anatoly Gelman, Gianluca Giavitto, Or Guttman, Juan M. Haces Crespo, Daniel Heilbrunn , et al. (23 additional authors not shown)

Abstract: The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is a space-borne near UV telescope with an unprecedented large field of view (200 sq. deg.). The mission, led by the Weizmann Institute of Science and the Israel Space Agency in collaboration with DESY (Helmholtz association, Germany) and NASA (USA), is fully funded and expected to be launched to a geostationary transfer orbit in Q2/3 of 202… ▽ More The Ultraviolet Transient Astronomy Satellite (ULTRASAT) is a space-borne near UV telescope with an unprecedented large field of view (200 sq. deg.). The mission, led by the Weizmann Institute of Science and the Israel Space Agency in collaboration with DESY (Helmholtz association, Germany) and NASA (USA), is fully funded and expected to be launched to a geostationary transfer orbit in Q2/3 of 2025. With a grasp 300 times larger than GALEX, the most sensitive UV satellite to date, ULTRASAT will revolutionize our understanding of the hot transient universe, as well as of flaring galactic sources. We describe the mission payload, the optical design and the choice of materials allowing us to achieve a point spread function of ~10arcsec across the FoV, and the detector assembly. We detail the mitigation techniques implemented to suppress out-of-band flux and reduce stray light, detector properties including measured quantum efficiency of scout (prototype) detectors, and expected performance (limiting magnitude) for various objects. △ Less

Submitted 11 March, 2023; v1 submitted 30 July, 2022; originally announced August 2022.

Comments: Presented in the SPIE Astronomical Telescopes + Instrumentation 2022

Population-Level Inference of Strong Gravitational Lenses with Neural Network-Based Selection Correction

Authors: Ronan Legin, Connor Stone, Yashar Hezaveh, Laurence Perreault-Levasseur

Abstract: A new generation of sky surveys is poised to provide unprecedented volumes of data containing hundreds of thousands of new strong lensing systems in the coming years. Convolutional neural networks are currently the only state-of-the-art method that can handle the onslaught of data to discover and infer the parameters of individual systems. However, many important measurements that involve strong l… ▽ More A new generation of sky surveys is poised to provide unprecedented volumes of data containing hundreds of thousands of new strong lensing systems in the coming years. Convolutional neural networks are currently the only state-of-the-art method that can handle the onslaught of data to discover and infer the parameters of individual systems. However, many important measurements that involve strong lensing require population-level inference of these systems. In this work, we propose a hierarchical inference framework that uses the inference of individual lensing systems in combination with the selection function to estimate population-level parameters. In particular, we show that it is possible to model the selection function of a CNN-based lens finder with a neural network classifier, enabling fast inference of population-level parameters without the need for expensive Monte Carlo simulations. △ Less

Submitted 8 July, 2022; originally announced July 2022.

Comments: 8 pages, 5 figures, accepted at the ICML 2022 Workshop on Machine Learning for Astrophysics

Revisiting Stellar Orbits and the Sgr A$^*$ Quadrupole Moment

Authors: Yael Alush, Nicholas Chamberlain Stone

Abstract: The "no-hair" theorem can, in principle, be tested at the center of the Milky Way by measuring the spin and the quadrupole moment of Sgr A$^*$ with the orbital precession of S-stars, measured over their full periods. Contrary to the original method, we show why it is possible to test the no-hair theorem using observations from only a single star, by measuring precession angles over a half-orbit. T… ▽ More The "no-hair" theorem can, in principle, be tested at the center of the Milky Way by measuring the spin and the quadrupole moment of Sgr A$^*$ with the orbital precession of S-stars, measured over their full periods. Contrary to the original method, we show why it is possible to test the no-hair theorem using observations from only a single star, by measuring precession angles over a half-orbit. There are observational and theoretical reasons to expect S-stars to spin rapidly, and we have quantified the effect of stellar spin, via spin-curvature coupling (the leading-order manifestation of the Mathisson-Papapetrou-Dixon equations), on future quadrupole measurements. We find that they will typically suffer from errors of order a few percentage points, but for some orbital parameters, the error can be much higher. We re-examine the more general problem of astrophysical noise sources that may impede future quadrupole measurements, and find that a judicious choice of measurable precession angles can often eliminate individual noise sources. We have derived optimal combinations of observables to eliminate the large noise source of mass precession, the novel noise of spin-curvature coupling due to stellar spin, and the more complicated noise source arising from transient quadrupole moments in the stellar potential. △ Less

Submitted 6 December, 2022; v1 submitted 5 July, 2022; originally announced July 2022.

Comments: Minor changes made for consistency with PRD version

The Spectroscopy and H-band Imaging of Virgo cluster galaxies (SHIVir) Survey: Data Catalogue and Kinematic Profiles

Authors: Nathalie N. -Q. Ouellette, Stéphane Courteau, Jon A. Holtzman, Michael McDonald, Michele Cappellari, Joel C. Roediger, Patrick Côté, Julianne J. Dalcanton, Elena Dalla Bontà, Laura Ferrarese, R. Brent Tully, Connor Stone, Eric W. Peng

Abstract: The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reductio… ▽ More The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reduction, and analysis for this SHIVir sample. We have used and optimised the \texttt{pPXF} routine to extract stellar kinematics from our data. Ultimately, resolved kinematic profiles (rotation curves and velocity dispersion profiles) are available for 133 SHIVir galaxies. A comprehensive database of photometric and kinematic parameters for the SHIVir sample is presented with: grizH magnitudes, effective surface brightnesses, effective and isophotal radii, rotational velocities, velocity dispersions, and stellar and dynamical masses. Parameter distributions highlight some bimodal distributions and possible sample biases. A qualitative study of resolved extended velocity dispersion profiles suggests a link between the so-called ``sigma-drop'' kinematic profile and the presence of rings in lenticular S0 galaxies. Rising dispersion profiles are linked to early-type spirals or dwarf ellipticals for which a rotational component is significant, whereas peaked profiles are tied to featureless giant ellipticals. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 21 pages, 16 figures, accepted for publication in MNRAS (doi: 10.1093/mnras/stac1347). Supplemental material available at https://www.queensu.ca/academia/courteau/links

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 2356-2375

Stream-Disk Shocks as the Origins of Peak Light in Tidal Disruption Events

Authors: Elad Steinberg, Nicholas C. Stone

Abstract: Tidal disruption events occur when stars are ripped apart by massive black holes, and result in highly luminous, multi-wavelength flares. Optical/UV observations of tidal disruption events (TDEs) contradict simple models of TDE emission, but the debate between alternative models (e.g. shock power or reprocessed accretion power remains unsettled, as the dynamic range of the problem has so far preve… ▽ More Tidal disruption events occur when stars are ripped apart by massive black holes, and result in highly luminous, multi-wavelength flares. Optical/UV observations of tidal disruption events (TDEs) contradict simple models of TDE emission, but the debate between alternative models (e.g. shock power or reprocessed accretion power remains unsettled, as the dynamic range of the problem has so far prevented ab initio hydrodynamical simulations. Consequently, past simulations have resorted to unrealistic parameter choices, artificial mass injection schemes or very short run-times. Here we present a 3D radiation-hydrodynamic simulation of a TDE flare from disruption to peak emission, with typical astrophysical parameters. At early times, shocks near pericenter power the light curve and a novel source of X-ray emission, but circularization and outflows are inefficient. Near peak light, stream-disk shocks efficiently circularize returning debris, power stronger outflows, and reproduce observed peak optical/UV luminosities. Peak emission in this simulation is shock-powered, but upper limits on accretion power become competitive near peak light as circularization runs away. This simulation shows how deterministic predictions of TDE light curves and spectra can be calculated using moving-mesh hydrodynamics algorithms. △ Less

Submitted 25 December, 2023; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: Updated to match final version accepted for publication. Movies at https://www.youtube.com/watch?v=O3IWCPO_Thk https://www.youtube.com/watch?v=_kpBQjwf3Z4

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.

Comments: 12 pages, 8 figures, 2 appendices

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

Authors: Diego J. Muñoz, Nicholas C. Stone, Cristobal Petrovich, Frederic A. Rasio

Abstract: We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) si… ▽ More We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 15 pages (before references and appendices), 11 figures, comments welcome

A Library of Synthetic X-ray Spectra for Fitting Tidal Disruption Events

Authors: Sixiang Wen, Peter G. Jonker, Nicholas C. Stone, Ann I. Zabludoff, Zheng Cao

Abstract: We present a tabulated version of our slim disk model for fitting tidal disruption events (TDEs). We create a synthetic X-ray spectral library by ray-tracing stationary general relativistic slim disks and including gravitational redshift, Doppler, and lensing effects self-consistently. We introduce the library to reduce computational expense and increase access for fitting future events. Fitting r… ▽ More We present a tabulated version of our slim disk model for fitting tidal disruption events (TDEs). We create a synthetic X-ray spectral library by ray-tracing stationary general relativistic slim disks and including gravitational redshift, Doppler, and lensing effects self-consistently. We introduce the library to reduce computational expense and increase access for fitting future events. Fitting requires interpolation between the library spectra; the interpolation error in the synthetic flux is generally $<10\%$ (it can rise to $40\%$ when the disk is nearly edge-on). We fit the X-ray spectra of the TDEs ASASSN-14li and ASASSN-15oi, successfully reproducing our earlier constraints on black hole mass $M_\bullet$ and spin $a_\bullet$ from full on-the-fly ray-tracing. We use the library to fit mock observational data to explore the degeneracies among parameters, finding that 1) hotter thermal disk and edge-on inclination angle spectra offer tighter constraints on $M_\bullet$ and $a_\bullet$; 2) the constraining power of spectra on $M_\bullet$ and $a_\bullet$ increases as a power-law with the number of X-ray counts, and the index of the power law is higher for hotter thermal disk spectra; 3) multi-epoch X-ray spectra partially break the degeneracy between $M_\bullet$ and $a_\bullet$; 4) the time-dependent level of X-ray absorption can be constrained from spectral fitting. The tabulated model and slim disk model are {\href{https://doi.org/10.25739/hfhz-xn60}{here.} △ Less

Submitted 16 May, 2022; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted for Publication in ApJ

The diversity of spiral galaxies explained

Authors: Matthew Frosst, Stéphane Courteau, Nikhil Arora, Connor Stone, Andrea V. Macciò, Marvin Blank

Abstract: An extensive catalog of spatially-resolved galaxy rotation curves and multi-band optical light profiles for 1752 observed spiral galaxies is assembled to explore the drivers of diversity in galaxy structural parameters, rotation curve shapes, and stellar mass profiles. Similar data were extracted from the NIHAO galaxy simulations to identify any differences between observations and simulations. Se… ▽ More An extensive catalog of spatially-resolved galaxy rotation curves and multi-band optical light profiles for 1752 observed spiral galaxies is assembled to explore the drivers of diversity in galaxy structural parameters, rotation curve shapes, and stellar mass profiles. Similar data were extracted from the NIHAO galaxy simulations to identify any differences between observations and simulations. Several parameters, including the inner slope "S" of a rotation curve (RC), were tested for diversity. Two distinct populations are found in observed and simulated galaxies; (i) blue, low mass spirals with stellar mass M* < 10^9.3 Msol and roughly constant "S", and (ii) redder, more massive and more diverse spirals with rapidly increasing "S". In all cases, the value of "S" seems equally contributed by the baryonic and non-baryonic (dark) matter. Diversity is shown to increase mildly with mass. Numerical simulations reproduce well most baryon-dominated galaxy parameter distributions, such as the inner stellar mass profile slope and baryonic scaling relations, but they struggle to match the full diversity of observed galaxy rotation curves (through "S") and most dark-matter-dominated parameters. To reproduce observations, the error broadening of the simulation's intrinsic spread of RC metrics would have to be tripled. The differences in various projections of observed and simulated scaling relations may reflect limitations of current sub-grid physics models to fully capture the complex nature of galaxies. For instance, AGNs are shown to have a significant effect on the shapes of simulated RCs. The inclusion of AGN feedback brings simulated and observed inner RC shapes into closer agreement. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 22 pages, 14 figures, submitted to MNRAS, comments welcome

Astrophysics with the Laser Interferometer Space Antenna

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

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

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

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

Massive black hole formation in dense stellar environments: Enhanced X-ray detection rates in high velocity dispersion nuclear star clusters

Authors: Vivienne F. Baldassare, Nicholas C. Stone, Adi Foord, Elena Gallo, Jeremiah P. Ostriker

Abstract: We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass… ▽ More We analyze Chandra X-ray Observatory imaging of 108 galaxies hosting nuclear star clusters (NSCs) to search for signatures of massive black holes (BHs). NSCs are extremely dense stellar environments with conditions that can theoretically facilitate massive BH formation. Recent work by Stone et al. (2017) finds that sufficiently dense NSCs should be unstable to the runaway growth of a stellar mass BH into a massive BH via tidal captures. Furthermore, there is a velocity dispersion threshold ($40\;\rm{km\;s^{-1}}$) above which NSCs should inevitably form a massive BH. To provide an observational test of these theories, we measure X-ray emission from NSCs and compare to the measured velocity dispersion and tidal capture runaway timescale. We find that NSCs above the $40\;\rm{km\;s^{-1}}$ threshold are X-ray detected at roughly twice the rate of those below (after accounting for contamination from X-ray binaries). These results are consistent with a scenario in which dense, high-velocity NSCs can form massive BHs, providing a formation pathway that does not rely on conditions found only at high redshift. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: Resubmitted to AAS Journals on March 1, 2022

Adelie: Continuous Address Space Layout Re-randomization for Linux Drivers

Authors: Ruslan Nikolaev, Hassan Nadeem, Cathlyn Stone, Binoy Ravindran

Abstract: While address space layout randomization (ASLR) has been extensively studied for user-space programs, the corresponding OS kernel's KASLR support remains very limited, making the kernel vulnerable to just-in-time (JIT) return-oriented programming (ROP) attacks. Furthermore, commodity OSs such as Linux restrict their KASLR range to 32 bits due to architectural constraints (e.g., x86-64 only support… ▽ More While address space layout randomization (ASLR) has been extensively studied for user-space programs, the corresponding OS kernel's KASLR support remains very limited, making the kernel vulnerable to just-in-time (JIT) return-oriented programming (ROP) attacks. Furthermore, commodity OSs such as Linux restrict their KASLR range to 32 bits due to architectural constraints (e.g., x86-64 only supports 32-bit immediate operands for most instructions), which makes them vulnerable to even unsophisticated brute-force ROP attacks due to low entropy. Most in-kernel pointers remain static, exacerbating the problem when pointers are leaked. Adelie, our kernel defense mechanism, overcomes KASLR limitations, increases KASLR entropy, and makes successful ROP attacks on the Linux kernel much harder to achieve. First, Adelie enables the position-independent code (PIC) model so that the kernel and its modules can be placed anywhere in the 64-bit virtual address space, at any distance apart from each other. Second, Adelie implements stack re-randomization and address encryption on modules. Finally, Adelie enables efficient continuous KASLR for modules by using the PIC model to make it (almost) impossible to inject ROP gadgets through these modules regardless of gadget's origin. Since device drivers (typically compiled as modules) are often developed by third parties and are typically less tested than core OS parts, they are also often more vulnerable. By fully re-randomizing device drivers, the last two contributions together prevent most JIT ROP attacks since vulnerable modules are very likely to be a starting point of an attack. Furthermore, some OS instances in virtualized environments are specifically designated to run device drivers, where drivers are the primary target of JIT ROP attacks. Our evaluation shows high efficiency of Adelie's approach. [full abstract is in the paper] △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 27th ACM International Conference on Architectural Support for Programming Languages and Operating Systems (ASPLOS '22), February 28 - March 4, 2022, Lausanne, Switzerland

Density Wakes due to Dynamical Friction in Cored Potentials

Authors: Karamveer Kaur, Nicholas C. Stone

Abstract: Dynamical friction is often modeled with reasonable accuracy by the widely used Chandrasekhar formula. However, in some circumstances, Chandrasekhar's local and uniform approximations can break down severely. An astrophysically important example is the "core stalling" phenomenon seen in N-body simulations of massive perturber inspiralling into the near-harmonic potential of a stellar system's cons… ▽ More Dynamical friction is often modeled with reasonable accuracy by the widely used Chandrasekhar formula. However, in some circumstances, Chandrasekhar's local and uniform approximations can break down severely. An astrophysically important example is the "core stalling" phenomenon seen in N-body simulations of massive perturber inspiralling into the near-harmonic potential of a stellar system's constant-density core (and possibly also in direct observations of dwarf galaxies with globular clusters). In this paper we use the linearized collisionless Boltzmann equation to calculate the global response of a cored galaxy to the presence of a massive perturber. We evaluate the density deformation, or wake, due to the perturber and study its geometrical structure to better understand the phenomenon of core stalling. We also evaluate the dynamical friction torque acting on perturber from the Lynden-Bell--Kalnajs (LBK) formula. In agreement with past work, we find that the dynamical friction force arising from corotating resonances is greatly weakened, relative to the Chandrasekhar formula, inside a constant density core. In contrast to past work, however, we find that a population of previously neglected high-order and non-corotating resonances sustain a minimum level of frictional torque at ~10 % of the torque from Chandrasekhar formula. This suggests that complete core stalling likely requires phenomena beyond the LBK approach; we discuss several possible explanations. Additionally, to study core stalling for multiple perturbers, we investigate approximate secular dynamical interactions (akin to Lidov-Kozai dynamics) between two perturbers orbiting a cored stellar system and derive a criterion for instability arising due to their close encounters. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: Submitted for publication to MNRAS; comments are welcome