Daily Summary for 2026-10-21

Introduction to the show: ident: Astrophysics Radio. The week's best astrophysics papers, unpacked for curious ears.

Vera: It's the twenty-first of October, twenty twenty-six, and this is the day's research.

Jocelyn: 25 new papers came out today.

Vera: Hi, I'm Vera.

Jocelyn: And I'm Jocelyn, here with Subrahmanyan.

Vera: We'll take the day in one pass, then pull out the papers we're staying with.

Jocelyn: Start us off with Impact of the equation of state on core collapse supernovae I.

The summary: Vera: Today is the twenty-first of October twenty twenty six.

Jocelyn: We mapped tidal tails in open clusters using Gaia DR3 data.

Subrahmanyan: These features show how clusters evolve and interact with their galactic environment.

Vera: The primary focus was detecting these tails in five specific open clusters.

Jocelyn: Which ones were those five clusters?

Subrahmanyan: BH 164, Alessi 2, NGC 2281, NGC 2354 and M67.

Vera: So we focused on those specific objects today.

Vera: We used unsupervised machine learning for detection.

Jocelyn: And statistical techniques like color-magnitude diagram matching?

Subrahmanyan: Plus principal component analysis, PCA.

Vera: We also used density-based clustering to pinpoint tails.

Jocelyn: What did those analyses reveal about the clusters?

Subrahmanyan: Tidal tails exist in every cluster studied.

Vera: These tails stretch forty to one hundred parsecs.

Jocelyn: And how many stars are in those tails?

Subrahmanyan: They contain between one hundred and two hundred stars.

Vera: Further investigation examined kinematics, photometry, and spatial distribution of the detected features.

Jocelyn: The luminosity functions of these tails resembled those of their parent clusters.

Subrahmanyan: They generally lacked massive stars compared to the main bodies and showed higher binary fractions than expected.

Vera: Significant rotation was first detected in M67 and NGC 2281 through this analysis.

Jocelyn: That finding is important for understanding cluster evolution.

Subrahmanyan: Indeed, it suggests a dynamic process at play within those systems.

Vera: We need to look closely at the stellar kinematics next.

Jocelyn: The photometric data provides constraints on the physical extent of these features.

Subrahmanyan: It seems the distribution confirms tidal stripping is active there.

Vera: This study also found that the luminosity functions resembled those of their parent clusters.

Jocelyn: That comparison helps us gauge how much material was lost.

Subrahmanyan: They generally lacked massive stars compared to the main bodies in these tails.

Vera: Higher binary fractions were also exhibited than we anticipated for these structures.

Jocelyn: We should focus on the rotational signatures detected in M67 and NGC 2281.

Subrahmanyan: That significant rotation was a key discovery from this analysis.

Vera: It opens new avenues for modeling cluster dynamics based on this data.

Jocelyn: The results provide concrete evidence for these tidal interactions in action.

Subrahmanyan: Agreed, the kinematic and photometric data are very illuminating here.

Vera: Today's lucky papers are Tidal Tails and Their Dynamics in Open Clusters Using Gaia DR3 because it uses machine learning to find and study tidal tails in five open clusters using Gaia data.

Jocelyn: Significant rotation was also detected in M67 and NGC 2281 for the first time.

Subrahmanyan: This paper is worth hearing because it detects significant rotation in M67 and NGC 2281 for the first time.

Vera: Agreed, that detection of rotation is a major finding.

Jocelyn: That research also analyzes stellar properties like kinematics and binaries to find these tails.

Subrahmanyan: Indeed, analyzing those specific properties reveals the underlying dynamics clearly.

Lucky paper: 2605.28944: Vera: We’re looking at the Impact of the equation of state on core collapse supernovae I right now.

Jocelyn: This paper looks at how different nuclear equation of state models affect a low-temperature W instability during core collapse in rapidly rotating stars.

Subrahmanyan: So, it uses three-dimensional neutrino magnetohydrodynamics simulations to explore the relationship between dense matter physics and these multimessenger signatures.

Vera: The main thing they found is that this low-T/W instability is robust across a broad spectrum of nuclear equation of state models tested for this rapidly rotating progenitor.

Jocelyn: It means the instability happens in every model they ran, which confirms it’s a fundamental feature for these kinds of progenitors.

Subrahmanyan: But the specific dynamics and signatures change depending on the underlying dense-matter physics encoded in those equation of state models.

Vera: They found that the onset time varies; for instance, the earliest onset of this spiral structure is around fifty ms after bounce in one model, like BH.

Jocelyn: That timing difference suggests that different physics leads to different initial conditions for the instability.

Subrahmanyan: Yes, while it's consistent that all models show a low-order spiral structure dominated by azimuthal modes m=one and m=two near the equatorial plane.

Vera: The key diagnostic here is how those characteristic frequencies correlate with the EOS properties of the proto-neutron star.

Jocelyn: The paper says that while the characteristic frequencies stay quasi-constant during its active phase, they are directly correlated with the effective stiffness and compactness of the proto-neutron star.

Subrahmanyan: They specifically say models characterized by a more compact or stiffer equation of state produce higher characteristic gravitational wave frequencies.

Vera: That suggests measuring this frequency could be a crucial diagnostic tool for constraining the dense-matter equation of state, complementing neutrino data.

Jocelyn: The resulting quasi-periodic modulations in both gravitational waves and neutrinos are characterized by specific typical values.

Subrahmanyan: They cite about three hundred thirty Hz for gravitational waves and about one hundred seventy Hz for neutrino emission, placing the low-T/W instability in a distinct spectral window.

Vera: That separation is important because it lies above the typical standing accretion shock instability range of one hundred Hz but below the dominant proto-neutron star oscillation frequencies at five hundred Hz.

Jocelyn: That distinction makes the low-T/W instability signature potentially unique and measurable in future multimessenger observations.

Subrahmanyan: Overall, this research concludes that the low-T/W instability is a robust physical feature of rapidly rotating core-collapse supernovae across a wide range of nuclear equation of state models.

Vera: Crucially, while its existence is universal, its characteristic multimessenger signatures retain measurable sensitivity to the stiffness and compactness parameters kappa two.

Jocelyn: So the frequency of this instability-driven gravitational wave signal offers a promising avenue for probing the properties of nuclear matter in these extreme astrophysical environments.

Lucky paper: 2403.05128: Vera: We’re moving on to the paper on High-energy extension of the gamma-ray band observable with an electron-tracking Compton camera because it tackles a real sensitivity gap in MeV gamma-ray astronomy.

Jocelyn: So, the core issue is that conventional Compton cameras struggle to determine the arrival directions of gamma rays uniquely.

Subrahmanyan: They can only constrain the incoming direction within an annulus region in the sky, which limits their practical angular resolution to worse than twenty degrees at six hundred sixty-two keV.

Vera: This new method uses an electron-tracking Compton camera, or ETCC, and looks for double-hit events where one hit comes from scattered gamma rays and the other from a recoil electron that escaped the tracker.

Jocelyn: How does geometrically identifying those scattering points help improve things?

Subrahmanyan: The reconstruction is achieved by geometrically identifying the Compton scattering point for these double-hit events.

Vera: The paper shows that this geometric identification leads to an angular resolution of eight point four degrees at the half-power radius at one point five MeV, which is twice as good as in the previous study.

Jocelyn: So, we’re talking about a factor of two improvement in resolution for these higher energy events?

Subrahmanyan: Exactly, and they managed to extend the observable energy range with the ETCC from zero point two–two point one MeV up to three point five MeV.

Vera: That extension is huge because it allows the detector to target things like the 26Al line emission at one point eight MeV or even molecular cloud de-excitation lines at four–six MeV and the COMPTEL excess.

Jocelyn: So, this method directly fills that sensitivity gap in the MeV band for future observations.

Subrahmanyan: The paper confirms this by showing that the total effective area for these double-hit events is almost constant at (five point zero–five point five) times ten to the minus two cm squared for the energy range of one point five to three point zero MeV, indicating a real extension of the observation band to higher energies.

Vera: It's interesting that this new analysis method was tested on laboratory data using a source like Carbon-sixty and the simulation result agreed with actual data within a factor of one point two.

Jocelyn: That agreement between simulation and real data is important for validating the technique before it goes into space.

Subrahmanyan: It shows that the reconstruction method works well in practice, even when compared to previous single-hit event analyses which had a worse resolution.

Lucky paper: 2601.04294: Vera: We're looking at a paper called A New Bayesian Framework with Mass-Radius-Space Priors to Constrain the Neutron Star Equation of State.

Jocelyn: This work proposes a new way to infer the neutron star equation of state by using mass and radius observations directly, instead of starting in pressure energy density space.

Subrahmanyan: They call this approach the Parameterized Mass-Radius, or PMR method, which uses an analytical scheme to invert the Tolman-Oppenheimer-Volkoff equations.

Vera: So they are taking a mass and radius curve and trying to map it back to the pressure energy density relation.

Jocelyn: They establish boundaries for this mass-radius space using physics informed priors, which is what they call establishing the M-R boundary.

Subrahmanyan: For low densities, they use a piecewise polytropic representation below a density of zero point zero four fm minus three cubic femtometers.

Vera: Above that density, they use the central value and one sigma uncertainty band from N3LO χEFT calculations of neutron star matter.

Jocelyn: Then they build a discrete mesh with horizontal lines separated by zero point one M solar masses to specify mass and angled lines separated by zero point two kilometers for radius.

Subrahmanyan: They anchor this grid to the fixed point (Mt, Rt) and determine the slope using the largest value of dM over dR between intersections of the right-hand boundary with grid points.

Vera: That resulting slope for that N3LO equation of state is found to be two point two M solar masses per kilometer minus one.

Jocelyn: They then construct all possible continuous mass-radius curves by connecting these grid points starting from the fixed point and ending at a mass greater than or equal to two point zero M solar masses.

Subrahmanyan: Each curve is checked for causality and thermodynamic stability at every density, and the probability of a given M-R curve is computed using Bayes’ theorem.

Vera: The analysis showed that the relative importance of prior modeling uncertainties versus observational uncertainties was a key issue here.

Jocelyn: They decomposed the total uncertainty into model uncertainty plus prior uncertainty, sigma two tot(E) equals sigma two model(E) plus sigma two prior(E).

Subrahmanyan: The PMR method is generally much faster because it avoids EOS parameterization and repeated TOV solutions.

Vera: And they found that the sensitivity to prior assumptions in the PMR method stays below fifty percent, which is lower than traditional frameworks.

Jocelyn: Traditional methods showed mean values around seventy percent sensitivity with typical variations of fifty to eighty percent.

Subrahmanyan: The conclusion is that the PMR method offers advantages because both natural priors and observational data are defined in the M-R space.

Vera: This suggests a need for alternative inference frameworks that offer broader coverage of prior space and reduced dependence on subjective prior assumptions.

Jocelyn: They also mentioned using spectral representations, like the Spectral4 parameterization, to cover the vast majority of faithful equations of state efficiently.

Subrahmanyan: Future work will explore more complex M-R relations, including those with multiple phase transitions, using improved inversion techniques to map those curves to corresponding equations of state.

Lucky paper: 2605.12503: Vera: We're looking at Unveiling Hidden Lyman Alpha Emitters in the DESI DR1 Data today.

Jocelyn: This paper is about using a convolutional neural network architecture to find Lyman alpha emitters hidden in the DESI data release one spectroscopic dataset.

Subrahmanyan: It focuses on solving a problem because current pipelines are not designed to detect and measure the redshifts of galaxies at z greater than two.

Vera: So how does this machine learning model actually work to do that?

Jocelyn: The architecture has two parts, first, it uses fully convolutional layers without any fully-connected layers, each one followed by a ReLU activation function and a max pooling layer to downscale the size of the input spectrum.

Subrahmanyan: That first part takes the blue arm signal-to-noise ratio spectrum, which is two thousand four hundred thirty-eight pixels with one pixel being zero point eight arcseconds per pixel, and outputs a lower resolution spectrum of seventy-six pixels with one pixel being approximately twenty-five arcseconds per pixel.

Vera: What does that lower resolution output actually represent for the model?

Jocelyn: Each output pixel value is between zero and one indicating the probability that the corresponding region in the original spectrum contains a Ly alpha emission line.

Subrahmanyan: Then there's a second part which is a combination of a slicer and a fitter. The slicer first identifies the pixel from that first part output with a probability greater than zero point five, and then extracts that corresponding region from the original spectrum.

Vera: After extracting that specific region, what happens to it next?

Jocelyn: That extracted spectrum is then passed to a fitter which consists of two convolution layers followed by three fully-connected layers. This fitter fits the red peak of the LAE Ly alpha emission line with a peak-normalized Gaussian function.

Subrahmanyan: The training dataset used six thousand five hundred sixty spectra with a one point one ratio between Lyman alpha emitters and non-Lyman alpha emitters.

Vera: And how did they balance out that imbalance for the model training?

Jocelyn: They added supplementary spectra consisting of one thousand two hundred thirteen quasars and another one thousand two hundred thirteen low-redshift emission line galaxies to balance the number of LAEs and non-LAEs.

Subrahmanyan: The performance results show a final model yields ninety-five point two percent purity and ninety-five point nine percent completeness for detecting LAEs in the test dataset.

Vera: How many actual Lyman alpha emitters did this model find when applied to the larger test set?

Jocelyn: Applying this model to about two point zero times ten six spectra targeted as emission-line galaxies, the study detected nineteen thousand six hundred eighty-five LAEs from z approximately two point zero to three point five within twelve minutes with a single GPU.

Subrahmanyan: Those detected LAEs are mostly at the bright end of the luminosity function with a Lyman alpha luminosity of greater than or equal to ten to the power of four plus thirty erg per second.

Vera: The paper also points out some specific issues that happened because of how the DESI pipeline works.

Jocelyn: They identified several misidentification patterns by the Redrock pipeline among the VI confirmed LAEs, specifically a lack of robust redshift detection for twenty-three out of sixty-two LAEs, which is about thirty-seven percent.

Subrahmanyan: Then there's another issue where twenty-one out of those same sixty-two LAEs, about thirty-four percent, are misidentified as QSOs with correct redshifts because the best-fit models show broader Ly alpha line profiles with lower amplitude than what was observed.

Vera: And what about the third category of misidentification mentioned in that paper?

Jocelyn: Another twenty-one percent of LAEs, which is thirteen out of sixty-two are misidentified as O II lines and other lines with incorrect redshifts, most having assigned redshifts between z equals zero point one and z equals zero point five.

Subrahmanyan: The final catalog consists of two samples, including those LAEs with a selection based on delta chi squared less than or equal to twenty.

Vera: Looking at the properties of this final sample, what kind of galaxies are these LAEs?

Jocelyn: The exploration shows that the LAEs are blue galaxies at the bright end of the Lyman alpha luminosity function with a Lyman alpha luminosity greater than or equal to ten to the power of four plus thirty erg per second from redshift two point zero to three point five.

Subrahmanyan: The composite spectrum also shows rich spectral information, including P-Cygni features for N v and C iv along with broad He ii emission and Fe II absorption/emission lines, and MgII emission lines in the longer wavelength that wasn't covered in previous studies.

Vera: What is the average velocity offset they estimated between the redshifts?

Jocelyn: The study estimates an average velocity offset being about two hundred fifty kilometers per second between those redshifts based on the Ly alpha emission lines and systematic redshifts of galaxies.

Subrahmanyan: The authors plan to keep improving their model by using the entire spectrum as the model input, which will allow their model to reduce the false positive rate using other spectral lines and also recover LAEs at z greater than three point five.

Vera: They are also planning to explore attention-based Transformer architectures to capture correlations between multiple emission lines like Ly alpha and CIV.

Jocelyn: That should further enhance the model's capability to distinguish LAEs associated with Active Galactic Nuclei, or AGNs.

Lucky paper: 2601.12315: Vera: We're looking at the MAGPI Survey today, which looks at the co-evolution of baryons and dark matter in star-forming disk-like galaxies at intermediate redshifts between zero point one and zero point eight five z.

Jocelyn: That paper uses spatially resolved kinematic data to figure out how dark matter content changes with galaxy structure.

Subrahmanyan: It uses threeDBarolo to model the gas kinematics, getting moment maps, surface brightness profiles, rotation velocity profiles, and velocity dispersion profiles.

Vera: They corrected the observed rotation velocity for inclination and gas pressure support to get intrinsic rotation curves.

Jocelyn: Then they decompose that total dynamical mass profile into stellar disk mass MDisk(R), bulge mass Mbulge(R), atomic gas HI(R), molecular gas H2(R), and dark matter MDM(R).

Subrahmanyan: The DM fraction fDM is calculated by comparing the dark matter component to the total dynamical mass within a given radius.

Vera: They defined two regions for this calculation: the inner region, which is within Re, and the outskirts, from Re to Ropt.

Jocelyn: The key finding here is that low-mass galaxies with a stellar mass Mstar less than one hundred nine point five Msun are dark matter dominated in both these regions.

Subrahmanyan: Specifically, they found an average fDM of about zero point eight five within both the inner region Re and the outskirts Ropt for these low-mass galaxies.

Vera: But then they looked at high-mass galaxies, Mstar greater than or equal to one thousand ten Msun.

Jocelyn: In those high-mass systems, the DM fraction is lower in the inner regions compared to the outskirts.

Subrahmanyan: That suggests an inside-out disk growth pattern for those larger galaxies.

Vera: They also showed that at a fixed circular velocity, low-mass galaxies have about forty percent higher DM fraction within Re than their high-mass counterparts.

Jocelyn: So the DM fraction decreases as the baryon surface density increases, as shown in Figure eight.

Subrahmanyan: That inverse correlation between dark matter fraction and baryon surface density is what they highlight.

Vera: They also found that galaxies with a central baryon surface density Sigma bar greater than or equal to one hundred eight point zero Msun per square kiloparsec cannot have a DM fraction less than the total observed baryonic mass.

Jocelyn: That’s a strong constraint on the distribution of matter in these systems.

Subrahmanyan: The study also showed that both high- and low-mass systems show no significant change, less than ten percent, in their DM fractions when looking out to Ropt.

Vera: It really points toward that inside-out growth model being consistent across different galaxy sizes.

Jocelyn: The paper title is The MAGPI Survey: co-evolution of baryons and dark matter in star-forming disk-like galaxies at zero point one z zero point eight five.

Subrahmanyan: That analysis provides a detailed look at how baryons and dark matter evolve together in these disk-like galaxies across that redshift range.