Source author record

Jarrod Hurley

Jarrod Hurley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.HE astro-ph astro-ph.IM astro-ph.CO gr-qc physics.comp-ph

Catalog footprint

What is connected

11works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

The Secret Lives of Open Clusters: a Multiwavelength Examination of Three Open Clusters

Star clusters are well known for their dynamical interactions, an outcome of their high stellar densities; in this paper we use multiwavelength observations to search for the unique outcomes of these interactions in three nearby Galactic open clusters: IC 2602 (30 Myr), NGC 2632 (750 Myr) and M67 (4 Gyr). We compared X-ray observations from all-sky surveys like eROSITA, plus archival observations from Chandra X-ray Observatory, survey radio observations from ASKAP's Evolutionary Map of the Universe survey plus archival VLA observations, in conjunction with new cluster catalogs with Gaia. From X-ray, we found 77 X-ray sources likely associated with IC 2602, 31 X-ray sources in NGC 2632, and 31 near M67's central regions. We were further able to classify these X-ray sources based on their optical variability and any radio emission. Three IC 2602 X-ray sources had radio counterparts, which are likely all chromospherically active binary stars. We also identified luminous radio and X-ray variability from a spectroscopic triple system in M67, WOCS 3012/S1077, which is either consistent with a quiescent black hole binary, or due to an active binary stellar system. A recent population study of optical variables by Anderson & Hunt 2025 shows that the population of optical variables in open clusters clearly changes over cluster age; this pilot study gives evidence that the X-ray population also changes with time, and demonstrates the need for a broader multiwavelength study of Galactic open clusters.

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

preprint2022arXiv

Dynamical double black holes and their host cluster properties

We investigate the relationship between the global properties of star clusters and their double black hole (DBH) populations. We use the code {\tt NBODY6} to evolve a suite of star cluster models with an initial mass of $\mathcal{O}(10^4)$M$_\odot$ and varying initial parameters. We conclude that cluster metallicity plays the most significant role in determining the lifespan of a cluster, while the initial half-mass radius is dominant in setting the rate of BH exchange interactions in the central cluster regions. We find that the mass of interacting BHs, rather than how frequently their interactions with other BHs occur, is more crucial in the thermal expansion and eventual evaporation of the cluster. We formulate a novel approach to easily quantify the degree of BH-BH dynamical activity in each model. We report 12 in-cluster and three out-of-cluster (after ejection from the cluster) DBH mergers, of different types (inspiral, eccentric, hierarchical) across the ten $N$-body models presented. Our DBH merger efficiency is 3--4$\times10^{-5}$ mergers per M$_\odot$. We note the cluster initial density plays the most crucial role in determining the number of DBH mergers, with the potential presence of a transitional density point (between 1.2-3.8$\times10^3$M$_\odot$/pc$^3$) below which the number of in-cluster mergers increases with cluster density and above which the increased stellar density acts to prevent in-cluster BH mergers. The importance of the history of dynamical interactions within the cluster in setting up the pathways to ejected DBH mergers is also discussed. Our findings show a broad match with observed LIGO-Virgo DBH mergers.

preprint2022arXiv

Explaining the differences in massive star models from various simulations

The evolution of massive stars is the basis of several astrophysical investigations, from predicting gravitational-wave event rates to studying star-formation and stellar populations in clusters. However, uncertainties in massive star evolution present a significant challenge when accounting for these models' behaviour in stellar population studies. In this work, we present a comparison between five published sets of stellar models from the BPASS, BoOST, Geneva, MIST, and PARSEC simulations at near-solar metallicity. The different sets of stellar models have been computed using slightly different physical inputs in terms of mass-loss rates and internal mixing properties. Moreover, these models also employ various pragmatic methods to overcome the numerical difficulties that arise due to the presence of density inversions in the outer layers of stars more massive than 40 M$_\odot$. These density inversions result from the combination of inefficient convection in the low-density envelopes of massive stars and the excess of radiative luminosity to the Eddington luminosity. We find that the ionizing radiation released by the stellar populations can change by up to 18 percent, the maximum radial expansion of a star can differ between 100-1600 R$_\odot$, and the mass of the stellar remnant can vary up to 20 M$_\odot$ between the five sets of simulations. We conclude that any attempts to explain observations that rely on the use of models of stars more massive than 40 M$_\odot$ should be made with caution.

preprint2020arXiv

The fates of massive stars: exploring uncertainties in stellar evolution with METISSE

In the era of advanced electromagnetic and gravitational wave detectors, it has become increasingly important to effectively combine and study the impact of stellar evolution on binaries and dynamical systems of stars. Systematic studies dedicated to exploring uncertain parameters in stellar evolution are required to account for the recent observations of the stellar populations. We present a new approach to the commonly used Single-Star Evolution (SSE) fitting formulae, one that is more adaptable: Method of Interpolation for Single Star Evolution (METISSE). It makes use of interpolation between sets of pre-computed stellar tracks to approximate evolution parameters for a population of stars. We have used METISSE with detailed stellar tracks computed by the Modules for Experiments in Stellar Astrophysics (MESA), Bonn Evolutionary Code (BEC) and Cambridge STARS code. METISSE better reproduces stellar tracks computed using the STARS code compared to SSE, and is on average three times faster. Using stellar tracks computed with MESA and BEC, we apply METISSE to explore the differences in the remnant masses, the maximum radial expansion, and the main-sequence lifetime of massive stars. We find that different physical ingredients used in the evolution of stars, such as the treatment of radiation dominated envelopes, can impact their evolutionary outcome. For stars in the mass range 9 to 100 M$_\odot$, the predictions of remnant masses can vary by up to 20 M$_\odot$, while the maximum radial expansion achieved by a star can differ by an order of magnitude between different stellar models.

preprint2013arXiv

Comparing the White Dwarf Cooling Sequences in 47 Tuc and NGC 6397

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

preprint2013arXiv

MOCCA Code for Star Cluster Simulations - II. Comparison with N-body Simulations

We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries.

preprint2011arXiv

A Deep, Wide-Field, and Panchromatic View of 47 Tuc and the SMC with HST: Observations and Data Analysis Methods

In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.

preprint2010arXiv

The long and the short of it: modelling double neutron star and collapsar Galactic dynamics

The work presented here examines populations of double compact binary systems and tidally enhanced collapsars. We make use of BINPOP and BINKIN, two components of a recently developed population synthesis package. Results focus on correlations of both binary and spatial evolutionary population characteristics. Pulsar and long duration gamma-ray burst observations are used in concert with our models to draw the conclusions that: double neutron star binaries can merge rapidly on timescales of a few million years (much less than that found for the observed double neutron star population), common envelope evolution within these models is a very important phase in double neutron star formation, and observations of long gamma-ray burst projected distances are more centrally concentrated than our simulated coalescing double neutron star and collapsar Galactic populations. Better agreement is found with dwarf galaxy models although the outcome is strongly linked to the assumed birth radial distribution. The birth rate of the double neutron star population in our models range from 4-160 Myr^-1 and the merger rate ranges from 3-150 Myr^-1. The upper and lower limits of the rates results from including electron capture supernova kicks to neutron stars and decreasing the common envelope efficiency respectively. Our double black hole merger rates suggest that black holes should receive an asymmetric kick at birth.

preprint2008arXiv

A Deep HST Study of the Globular Cluster NGC 6397: Reduction Methods

We describe here the reduction methods that we developed to study the faintest red dwarfs and white dwarfs in an outer field of NGC6397, which was observed by \hst for 126 orbits in 2005. The particular challenge of this data set is that the faintest stars are not readily visible in individual exposures, so special care must be taken to combine the information in all the exposures in order to identify and measure them. Unfortunately, it is hard to find the faintest stars without also finding a large number of faint galaxies, so we developed specialized tools to distinguish between the point-like stars and the barely resolved galaxies. We found that artificial-star tests, while obviously necessary for completeness determination, can also play an important role in helping us optimize our finding and measuring algorithms. Although this paper focuses on this data set specifically, many of the techniques are new and might find application in other work, particularly when a large number of images is available for a single field.

preprint1998arXiv

The Binary Second Sequence in Cluster Colour--Magnitude Diagrams

We show how the second sequence seen lying above the main sequence in cluster colour magnitude diagrams results from binaries with a large range of mass ratios and not just from those with equal masses. We conclude that the presence of a densely populated second sequence, with only sparse filling in between it and the single star main sequence, does not necessarily imply that binary mass ratios are close to unity.