Stars In M31 Research Articles

Uncertainty of Line-of-sight Velocity Measurements of Faint Stars from Low- and Medium-resolution Optical Spectra

Massively multiplexed spectrographs will soon gather large statistical samples of stellar spectra. The accurate estimation of uncertainties on derived parameters, such as the line-of-sight velocity v los, especially for spectra with low signal-to-noise ratios (S/Ns), is paramount. We generated an ensemble of simulated optical spectra of stars as if they were observed with low- and medium-resolution fiber-fed instruments on an 8 m class telescope, similar to the Subaru Prime Focus Spectrograph, and determined v los by fitting stellar templates to the simulated spectra. We compared the empirical errors of the derived parameters—calculated from an ensemble of simulations—to the asymptotic errors determined from the Fisher matrix, as well as from Monte Carlo sampling of the posterior probability. We confirm that the uncertainty of v los scales with the inverse square root of the S/N, but also show how this scaling breaks down at low S/N and analyze the error and bias caused by template mismatch. We outline a computationally optimized algorithm to fit multiexposure data and provide a mathematical model of stellar spectrum fitting that maximizes the so called significance, which allows for calculating the error from the Fisher matrix analytically. We also introduce the effective line count, and provide a scaling relation to estimate the errors of v los measurements based on stellar type. Our analysis covers a range of stellar types with parameters that are typical of the Galactic outer disk and halo, together with analogs of stars in M31 and in satellite dwarf spheroidal galaxies around the Milky Way.

TOI-2015 b: A Warm Neptune with Transit Timing Variations Orbiting an Active Mid-type M Dwarf

We report the discovery of a close-in (P orb = 3.349 days) warm Neptune with clear transit timing variations (TTVs) orbiting the nearby (d = 47.3 pc) active M4 star, TOI-2015. We characterize the planet's properties using Transiting Exoplanet Survey Satellite (TESS) photometry, precise near-infrared radial velocities (RVs) with the Habitable-zone Planet Finder Spectrograph, ground-based photometry, and high-contrast imaging. A joint photometry and RV fit yields a radius Rp=3.37−0.20+0.15R⊕ , mass mp=16.4−4.1+4.1M⊕ , and density ρp=2.32−0.37+0.38gcm−3 for TOI-2015 b, suggesting a likely volatile-rich planet. The young, active host star has a rotation period of P rot = 8.7 ± 0.9 days and associated rotation-based age estimate of 1.1 ± 0.1 Gyr. Though no other transiting planets are seen in the TESS data, the system shows clear TTVs of super-period Psup≈430days and amplitude ∼100 minutes. After considering multiple likely period-ratio models, we show an outer planet candidate near a 2:1 resonance can explain the observed TTVs while offering a dynamically stable solution. However, other possible two-planet solutions—including 3:2 and 4:3 resonances—cannot be conclusively excluded without further observations. Assuming a 2:1 resonance in the joint TTV-RV modeling suggests a mass of mb=13.3−4.5+4.7M⊕ for TOI-2015 b and mc=6.8−2.3+3.5M⊕ for the outer candidate. Additional transit and RV observations will be beneficial to explicitly identify the resonance and further characterize the properties of the system.

Abundances of Neutron-capture Elements in 62 Stars in the Globular Cluster Messier 15

M15 is a globular cluster with a known spread in neutron-capture elements. This paper presents abundances of neutron-capture elements for 62 stars in M15. Spectra were obtained with the Michigan/Magellan Fiber System spectrograph, covering a wavelength range from ∼4430 to 4630 Å. Spectral lines from Fe i, Fe ii, Sr i, Zr ii, Ba ii, La ii, Ce ii, Nd ii, Sm ii, Eu ii, and Dy ii were measured, enabling classifications and neutron-capture abundance patterns for the stars. Of the 62 targets, 44 are found to be highly Eu-enhanced r-II stars, another 17 are moderately Eu-enhanced r-I stars, and one star is found to have an s-process signature. The neutron-capture patterns indicate that the majority of the stars are consistent with enrichment by the r-process. The 62 target stars are found to show significant star-to-star spreads in Sr, Zr, Ba, La, Ce, Nd, Sm, Eu, and Dy, but no significant spread in Fe. The neutron-capture abundances are further found to have slight correlations with sodium abundances from the literature, unlike what has been previously found; follow-up studies are needed to verify this result. The findings in this paper suggest that the Eu-enhanced stars in M15 were enhanced by the same process, that the nucleosynthetic source of this Eu pollution was the r-process, and that the r-process source occurred as the first generation of cluster stars was forming.

JWST/NIRCam Detection of the Fomalhaut C Debris Disk in Scattered Light

Observations of debris disks offer important insights into the formation and evolution of planetary systems. Though M dwarfs make up approximately 80% of nearby stars, very few M dwarf debris disks have been studied in detail—making it unclear how or if the information gleaned from studying debris disks around more massive stars extends to the more abundant M dwarf systems. We report the first scattered-light detection of the debris disk around the M4 star Fomalhaut C using JWST's Near Infrared Camera (NIRCam; 3.6 and 4.4 μm). This result adds to the prior sample of only four M dwarf debris disks with detections in scattered light and marks the latest spectral type and oldest star among them. The size and orientation of the disk in these data are generally consistent with the prior Atacama Large Millimeter/submillimeter Array submillimeter detection. Though no companions are identified, these data provide strong constraints on their presence—with sensitivity sufficient to recover sub-Saturn mass objects in the vicinity of the disk. This result illustrates the unique capability of JWST to uncover elusive M dwarf debris disks in scattered light and lays the groundwork for deeper studies of such objects in the 2–5 μm regime.

Asteroseismology of three evolved stars in M67: testing systematic biases in seismic masses and ages

ABSTRACT Asteroseismology can precisely characterize stars, but the accuracy of seismic inference is still an open question. In this work, we use three evolved stars, including one late subgiant (EPIC 211411922) and two young red giants (EPIC 211409560 and EPIC 211416749) in the open cluster M67 to validate seismic inferences of stellar mass and age. We start by measuring oscillation mode frequencies and identifying acoustic and mixed modes in all three stars. We then fit the stars to determine their masses and ages with a detailed modelling approach. Our age estimates for all stars agree well with the literature consensus age range of M67 (3.0–5.0 Gyr). Mass estimates for EPIC 211411922 and EPIC 211409560 are sensible because they are slightly higher than the typical turn-off mass of the cluster. The inferred mass for EPIC 211416749 is underestimated by 5–10 per cent compared with the turn-off mass, but this could be caused by poor data quality. Our results indicate that the detailed modelling can determine sensible masses and ages for evolved stars near the base of the red giant branch.

Ring nebulae around Wolf–Rayet stars in M33 as seen by SITELLE

ABSTRACT We have conducted an analysis of nebulae around Wolf–Rayet (WR) stars in M33 using data collected by the imaging Fourier transform spectrometer SITELLE at the Canada–France–Hawaii telescope as part of the SIGNALS Large Program. Of the 211 known WR stars in M33, 178 are located in the fields observed in this study. We present the results of this analysis in the form of a comprehensive summary of all nebulae found around the observed WR stars. Based on three criteria we find to be the most effective for their detection, we detect a clear association with a circumstellar bubble around 33 of them (19 per cent). Our results show that the presence of bubbles does not correlate with the spectral type of the central star. The mean diameter of the WR nebulae we have found is 21 parsec.

Prospects for Detecting Gaps in Globular Cluster Stellar Streams in External Galaxies with the Nancy Grace Roman Space Telescope

Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming Nancy Grace Roman Space Telescope (Roman) presents a unique opportunity to increase the number of observed streams and gaps significantly. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar 5–like stream and generate mock observations of these gaps with background stars in M31 and foreground Milky Way stellar fields. We assess Roman's ability to detect gaps out to 10 Mpc through visual inspection and with the gap-finding tool FindTheGap. We conclude that gaps of ≈1.5 kpc in streams that are created from subhalos of masses ≥5 × 106 M ⊙ are detectable within a 2–3 Mpc volume in exposure times of 1000 s to 1 hr. This volume contains ≈150 galaxies, including ≈eight galaxies with luminosities >109 L ⊙. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models.

A Comparative Study between M30 and M92: M92 is a Merger Remnant with a Large Helium Enhancement

We perform a comparative study of the ex situ second-parameter pair globular clusters (GCs) M30 and M92, having similar metallicities but different horizontal branch morphologies. We obtain similar mean primordial carbon abundances for both clusters. However, M92 shows a large dispersion in carbon due to a more extended C–N anticorrelation, while M30 exhibits a higher primordial nitrogen abundance, suggesting that they have different chemical enrichment histories. Our new results confirm our previous result that M92 is a metal-complex GC showing a bimodal metallicity distribution. We also find that the metal-rich group of stars in M92 shows a helium enhancement as large as ΔY ∼ 0.05 from the red giant branch bump V magnitudes, which can also be supported by (i) a lack of bright red giant branch stars, (ii) synthetic evolutionary horizontal branch population models and (iii) the more extended spatial distribution due to different degree of the diffusion process from their lower masses. We reinterpret the [Eu/Fe] measurements by others, finding that the two metallicity groups of stars in M92 have significantly different [Eu/Fe] abundances with small scatters. This strongly suggests that they formed independently out of well-mixed interstellar media in different environments. We suggest that M92 is a more complex system than a normal GC, most likely a merger remnant of two GCs or an even more complex system. In the Appendix, we address the problems with the recently developed color–temperature relations and the usage of broadband photometry in the populational taggings.

Water-rich Disks around Late M Stars Unveiled: Exploring the Remarkable Case of Sz 114

We present an analysis of the JDISCS JWST/MIRI-MRS spectrum of Sz 114, an accreting M5 star surrounded by a large dust disk with a shallow gap at ∼39 au. The spectrum is molecule-rich; we report the detection of water, CO, CO2, HCN, C2H2, and H2. The only identified atomic/ionic transition is from [Ne ii] at 12.81 μm. A distinct feature of this spectrum is the forest of water lines with the 17.22 μm emission surpassing that of most mid-to-late M star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, the flux ratios of C2H2/H2O and HCN/H2O in Sz 114 also resemble those of earlier-type disks, with a slightly elevated CO2/H2O ratio. While accretional heating can boost all infrared lines, the unusual properties of Sz 114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and the presence of dust substructures. The latter delays the inward drift of icy pebbles and helps preserve a lower C/O ratio over an extended period. In contrast, mid-to-late M-star disks—which are typically faint, small in size, and likely lack significant substructures—may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of mid-to-late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models.

R-process Abundance Patterns in the Globular Cluster M92

Whereas light-element abundance variations are a hallmark of globular clusters, there is little evidence for variations in neutron-capture elements. A significant exception is M15, which shows a star-to-star dispersion in neutron-capture abundances of at least one order of magnitude. The literature contains evidence both for and against a neutron-capture dispersion in M92. We conducted an analysis of archival Keck/HIRES spectra of 35 stars in M92, 29 of which are giants, which we use exclusively for our conclusions. M92 conforms to the abundance variations typical of massive clusters. Like other globular clusters, its neutron-capture abundances were generated by the r-process. We confirm a star-to-star dispersion in r-process abundances. Unlike M15, the dispersion is limited to “first-generation” (low-Na, high-Mg) stars, and the dispersion is smaller for Sr, Y, and Zr than for Ba and the lanthanides. This is the first detection of a relation between light-element and neutron-capture abundances in a globular cluster. We propose that a source of the main r-process polluted the cluster shortly before or concurrently with the first generation of star formation. The heavier r-process abundances were inhomogeneously distributed while the first-generation stars were forming. The second-generation stars formed after several crossing times (∼0.8 Myr); hence, the second generation shows no r-process dispersion. This scenario imposes a minimum temporal separation of 0.8 Myr between the first and second generations.

M67 Blue Stragglers with High-resolution Infrared Spectroscopy

We report on the first detailed infrared chemical analysis of five binary members (S277, S997, S975, S1031, and S1195) in the open cluster M67 (NGC 2682). These stars are located outside (bluer and/or brighter than) the main-sequence turnoff region in M67. High-resolution (R ∼ 45,000) near-infrared spectra were obtained with the Immersion GRating INfrared Spectrograph (IGRINS) at the McDonald Observatory 2.7 m Harlan J. Smith Telescope, providing full spectral coverage of the H and K bands. The abundances of C, Na, Mg, Al, Si, S, Ca, Fe, and Ni are measured using neutral atomic absorption lines. We detect v sin i ≥ 25 km s−1 in three of our program stars: S1031, S975, and S1195. We find our derived abundances to be in good agreement with turnoff star abundances, similar to published analyses of blue straggler stars in M67 from optical spectra. Detection of a carbon enhancement or depletion resulting from mass transfer is difficult due to the uncertainties in the carbon abundance and the relatively modest changes that may occur through red giant and asymptotic giant branch evolution.

The Panchromatic Hubble Andromeda Treasury. XX. The Disk of M31 is Thick

We present a new approach to measuring the thickness of a partially face-on stellar disk, using dust geometry. In a moderately-inclined disk galaxy, the fraction of reddened stars is expected to be 50% everywhere, assuming that dust lies in a thin midplane. In a thickened disk, however, a wide range of radii project onto the line of sight. Assuming stellar density declines with radius, this geometrical projection leads to differences in the numbers of stars on the near and far sides of the thin dust layer. The fraction of reddened stars will thus differ from the 50% prediction, with a deviation that becomes larger for puffier disks. We map the fraction of reddened red giant branch (RGB) stars across M31, which shows prominent dust lanes on only one side of the major axis. The fraction of reddened stars varies systematically from 20% to 80%, which requires that these stars have an exponential scale height h z that is 0.14 ± 0.015 times the exponential scale length (h r ≈ 5.5 kpc). M31's RGB stars must therefore have h z = 770 ± 80 pc, which is far thicker than the Milky Way’s thin disk, but comparable to its thick disk. The lack of a significant thin disk in M31 is unexpected, but consistent with its interaction history and high disk velocity dispersion. We suggest that asymmetric reddening be used as a generic criterion for identifying “thick disk”-dominated systems, and discuss prospects for future 3D tomographic mapping of the gas and stars in M31.

Newly Discovered Wolf–Rayet Stars in M31

The evolved massive star populations of the Local Group galaxies are generally thought to be well understood. However, recent work has suggested that the Wolf–Rayet (WR) content of M31 may have been underestimated. We therefore began a pilot project to search for new WRs in M31 and to reexamine the completeness of our previous WR survey, finished almost a decade prior. Our improved imaging data and spectroscopic follow-up confirmed 19 new WRs across three small fields in M31. These newly discovered WRs are generally fainter than the previously known sample due to slightly increased reddening as opposed to intrinsic faintness. From these findings, we estimate that there are another ∼60 WRs left to be discovered in M31; however, the overall ratio of WN-type (nitrogen-rich) to WC-type (carbon-rich) WRs remains unchanged with our latest additions to the M31 WR census. We are in the process of extending this pilot WR survey to include the rest of M31, and a more complete population will be detailed in our future work.

Effects of Full-sphere Convection on M-dwarf Dynamo Action, Flux Emergence, and Spin-down

We present an analysis of 21 global MHD simulations of the dynamo action taking place in the interiors of fully convective M4 stars. We leverage our past work involving 45 simulations of shell-convecting-type M2 stars in a very similar parameter space to identify properties that their dynamos share and ways in which they differ. We find that the parameter scalings of many fundamental dynamo measures including the ratio of magnetic to kinetic energy and the cycle period are shared not only between the two sets of M-dwarf models but also with the broader literature on K and G dwarfs, suggesting that their convective dynamos are similar in nature. We find that the surface fields of the fully convective stars are stronger and more organized than those in the shell-convecting models and indicate faster spin-down than is observed, for which we suggest possible remediations. We note an abundance of dynamo states that are localized within one hemisphere at a time and find that such states are well described by a quasi-linear superposition of even- and odd-parity eigenmodes with similar amplitudes. We observe that the turbulent stirring of the deep fields of the fully convective models allows them to attain far higher peak amplitudes than those in the tachoclines of the partially convective models. We employ our AI-enabled pipeline, LoopNet, to search for buoyant, fibril magnetic fields in the simulated stars and find nearly twice as many such structures forming in the M4 models, in good agreement with observed differences in their flaring rates.

Elemental Abundances in M31: Individual and Coadded Spectroscopic [Fe/H] and [α/Fe] throughout the M31 Halo with SPLASH

We present spectroscopic chemical abundances of red giant branch stars in Andromeda (M31), using medium-resolution (R ∼ 6000) spectra obtained via the Spectroscopic and Photometric Landscape of Andromeda’s Stellar Halo survey. In addition to individual chemical abundances, we coadd low signal-to-noise ratio spectra of stars to obtain a high enough signal to measure average [Fe/H] and [α/Fe] abundances. We obtain individual and coadded measurements for [Fe/H] and [α/Fe] for M31 halo stars, covering a range of 9–180 kpc in projected radius from the center of M31. With these measurements, we greatly increase the number of outer halo (R proj > 50 kpc) M31 stars with spectroscopic [Fe/H] and [α/Fe], adding abundance measurements for 45 individual stars and 33 coadds from a pool of an additional 174 stars. We measure the spectroscopic metallicity ([Fe/H]) gradient, finding a negative radial gradient of −0.0084 ± 0.0008 for all stars in the halo, consistent with gradient measurements obtained using photometric metallicities. Using the first measurements of [α/Fe] for M31 halo stars covering a large range of projected radii, we find a positive gradient (+0.0027 ± 0.0005) in [α/Fe] as a function of projected radius. We also explore the distribution in [Fe/H]–[α/Fe] space as a function of projected radius for both individual and coadded measurements in the smooth halo, and compare these measurements to those stars potentially associated with substructure. These spectroscopic abundance distributions add to existing evidence that M31 has had an appreciably different formation and merger history compared to our own Galaxy.

Photometric variability of blue straggler stars in M67 with TESS and K2

ABSTRACT Blue straggler stars (BSSs) are formed through mass transfers or mergers in binaries. The recent detections of white dwarf (WD) companions to BSSs in M67 suggested a mass transfer pathway of formation. In search of a close companion to five BSSs in M67 that are known to be spectroscopic binaries, we study the light curves from K2 and TESS data. We use PHOEBE to analyse the light curves and estimate the properties of the companions. We detect variability in WOCS 1007, and the light curve is dominated by ellipsoidal variation. Using the light curve and radial velocity measurements, we estimate its orbital period to be 4.212 ± 0.041 d and e = 0.206 ± 002. The mass of the companion is estimated to be 0.22 ± 0.05 M⊙ with a radius of 0.078 ± 0.027 R⊙, confirming it to be an LM WD with Teff = 14300 ± 1100 K. The estimated mass of the BSS, 1.95 ± 0.26 M⊙, is similar to that estimated from isochrones. The BSS in WOCS 1007 shows δ Scuti pulsations, although it is slightly deformed and likely to be formed through an efficient mass transfer. Though we detect a light curve for WOCS 4003 showing grazing eclipse with ellipsoidal variation, the estimated parameters are inconclusive. Apart from the 0.44 d period, we found smaller eclipses with a period of 1.1 d, suggesting a compact triple system. In the cases of WOCS 4003, WOCS 5005, and WOCS 1025, no eclipses or pulsations are detected, confirming the absence of any short-period inner binary with high inclination in these BSSs.

A 1.55 R⊕ habitable-zone planet hosted by TOI-715, an M4 star near the ecliptic South Pole

ABSTRACT A new generation of observatories is enabling detailed study of exoplanetary atmospheres and the diversity of alien climates, allowing us to seek evidence for extraterrestrial biological and geological processes. Now is therefore the time to identify the most unique planets to be characterized with these instruments. In this context, we report on the discovery and validation of TOI-715 b, a $R_{\rm b}=1.55\pm 0.06\rm R_{\oplus }$ planet orbiting its nearby (42 pc) M4 host (TOI-715/TIC 271971130) with a period $P_{\rm b} = 19.288004_{-0.000024}^{+0.000027}$ d. TOI-715 b was first identified by TESS and validated using ground-based photometry, high-resolution imaging and statistical validation. The planet’s orbital period combined with the stellar effective temperature $T_{\rm eff}=3075\pm 75~\rm K$ give this planet an installation $S_{\rm b} = 0.67_{-0.20}^{+0.15}~\rm S_\oplus$, placing it within the most conservative definitions of the habitable zone for rocky planets. TOI-715 b’s radius falls exactly between two measured locations of the M-dwarf radius valley; characterizing its mass and composition will help understand the true nature of the radius valley for low-mass stars. We demonstrate TOI-715 b is amenable for characterization using precise radial velocities and transmission spectroscopy. Additionally, we reveal a second candidate planet in the system, TIC 271971130.02, with a potential orbital period of $P_{02} = 25.60712_{-0.00036}^{+0.00031}$ d and a radius of $R_{02} = 1.066\pm 0.092\, \rm R_{\oplus }$, just inside the outer boundary of the habitable zone, and near a 4:3 orbital period commensurability. Should this second planet be confirmed, it would represent the smallest habitable zone planet discovered by TESS to date.

A Census of the Low Accretors. II. Accretion Properties

Much is known about the processes driving accretion from protoplanetary disks onto low-mass pre-main-sequence stars (T Tauri stars). Nevertheless, it is unclear how accretion stops. To determine the accretion properties and their relation to stellar properties and gain insight into the last stages of accretion, we present a detailed analysis of 24 low and possible accretors, previously identified using the He i λ10830 line. We model moderate-resolution Hα profiles of these stars using magnetospheric accretion flow models that account for the chromospheric contribution at the line center. Based on parameters derived from the fits of 20 stars that can be reproduced with the models, we find a power-law relation between the disk truncation radius and the mass accretion rate consistent with predictions from theory and simulations. Comparing the corotation and truncation radii, we find that most of our targets are accreting in the unstable regime and rule out the propeller as the main process stopping accretion. For the truncation radius to be the same as the magnetic radius, the dipole magnetic field and/or the efficiency parameter ξ need to be smaller than previously determined, suggesting that higher-order fields dominate in low accretion rates. Lastly, we determine that the lowest accretion rates that can be detected by Hα line modeling are (1–3) × 10−11 M ☉ yr−1 for M3 stars and (3–5) × 10−11 M ☉ yr−1 for K5 stars. These limits are lower than the observed accretion rates in our sample, suggesting that we have reached a physical lower limit. This limit, , is consistent with EUV-dominated photoevaporation.

The Prevalence of the α-bimodality: First JWST α-abundance Results in M31

AbstractWe present initial results from our JWST NIRSpec program to study the α-abundances in the M31 disk. The Milky Way has two chemically-defined disks, the low-α and high-α disks, which are closely related to the thin and thick disks, respectively. The origin of the two populations and the α-bimodality between them is not entirely clear, although there are now several models that can reproduce the observed features. To help constrain the models and discern the origin, we have undertaken a study of the chemical abundances of the M31 disk using JWST NIRSpec, in order to determine whether stars in M31’s disk also show an α-abundance bimodality. Approximately 100 stars were observed in our single NIRSpec field at a projected distance of 18 kpc from the M31 center. The 1-D extracted spectra have an average signal-to-noise ratio of 85 leading to statistical metallicity precision of 0.016 dex, α-abundance precision of 0.012 dex, and a radial velocity precision 8 km s-1 (mostly from systematics). The initial results indicate that, in contrast to the Milky Way, there is no α-bimodality in the M31 disk, and no low-α sequence. The entire stellar population falls along a single chemical sequence very similar to the MW’s high-α component which had a high star formation rate. While this is somewhat unexpected, the result is not that surprising based on other studies that found the M31 disk has a larger velocity dispersion than the MW and is dominated by a thick component. M31 has had a more active accretion and merger history than the MW which might explain the chemical differences.

Searching for the extra-tidal stars of globular clusters using high-dimensional analysis and a core particle spray code

ABSTRACT Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data analysis and a particle spray code to identify extra-tidal stars of any Galactic globular cluster using M3 as a case study. Using the t-Stochastic Neighbour Embedding and Uniform Manifold Approximation and Projection machine learning dimensionality reduction algorithms, we first identify a set of 103 extra-tidal candidates in the APOGEE DR17 data catalogue with chemical abundances similar to M3 stars. To confirm each candidate’s extra-tidal nature, we introduce corespray – a new python-based three-body particle spray code that simulates extra-tidal stars for any Galactic globular cluster. Using Gaia EDR3 proper motions and APOGEE DR17 radial velocities, we apply multivariate Gaussian modelling and an extreme deconvolution to identify the extra-tidal candidates that are more likely to be associated with a distribution of corespray-simulated M3 extra-tidal stars than the field. Through these methods, we identify 10 new high-probability extra-tidal stars produced via three-body interactions in M3. We also explore whether any of our extra-tidal candidates are consistent with being ejected from M3 through different dynamical processes. Future applications of corespray will yield better understandings of core dynamics, star formation histories, and binary fractions in globular clusters.