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: 23 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 diving into the paper, Impact of the equation of state on core collapse supernovae I: the low-T/W instability.

Jocelyn: So we know this instability is robust across all nuclear equation of state models considered for this rapidly rotating progenitor.

Subrahmanyan: That means it happens regardless of what specific physics is involved in those different models, which is a big deal.

Vera: Exactly, but the details about when it starts and what gravitational waves it produces really change depending on the equation of state.

Jocelyn: The study looked at how variations in nuclear physics, like differences in stiffness at supranuclear densities characterized by kappa two affect this instability.

Subrahmanyan: It systematically explores how choosing between Skyrme and Relativistic Mean Field type approaches changes the outcome.

Vera: The onset time varies; for instance, the earliest onset of this spiral structure is observed around fifty ms after bounce in some models like BH at approximately fifty ms.

Jocelyn: But in other models, like BB, the pattern only becomes visible later around one hundred ms.

Subrahmanyan: Despite those timing differences and model specifics, the overall structure remains coherent; all models show a low-order spiral structure dominated by azimuthal modes m=one and m=two near the equatorial plane.

Vera: The most important part for multimessenger astronomy is that while the characteristic frequencies remain quasi-constant during its active phase, they are directly correlated with the equation of state properties.

Jocelyn: Specifically, the dominant gravitational wave frequency correlates approximately linearly with the effective stiffness and compactness of the proto-neutron star.

Subrahmanyan: Models characterized by a more compact or stiffer equation of state produce higher characteristic gravitational wave frequencies.

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

Jocelyn: The resulting quasi-periodic modulations are characterized by specific typical values: about three hundred thirty Hz for gravitational waves and about one hundred seventy Hz for neutrino emission.

Subrahmanyan: That places the low-T/W instability in a distinct spectral window, above the Standing Accretion Shock Instability range of one hundred Hz but below the dominant proto-neutron star oscillation frequencies of five hundred Hz.

Vera: So we’re looking at this as a potential unique signature because it sits right in that gap.

Jocelyn: The paper impacts how we probe the properties of nuclear matter in these extreme astrophysical environments through its multimessenger signatures.

Lucky paper: 2509.03802: Vera: So we’re moving onto Structure of Turbulent Radiative Mixing Layers.

Jocelyn: This paper looks at turbulent radiative mixing layers, or TRMLs, and how they form in the interstellar medium.

Subrahmanyan: They show that a positive vertical compressive stress Rzz drives the TRML pressure to dip, which establishes a universal steady-state structure across all nuclear equation of state models.

Vera: That is interesting because it suggests this structure doesn't depend on the specific physics in the equation of state.

Jocelyn: But they also found that details like when it starts and what kind of gravitational waves it produces change depending on the model.

Subrahmanyan: The dominant gravitational wave frequency is linked to how compact and stiff the proto-neutron star is, which implies that this instability's frequency might be a probe of the equation of state.

Vera: So we have universal structure but model-dependent dynamics.

Jocelyn: This paper, Structure of Turbulent Radiative Mixing Layers, explores these turbulent radiative mixing layers where shear drives turbulence between dense and diffuse gas.

Subrahmanyan: They found that in quasi-isobaric TRMLs, upstream gas cools and compresses before streamwise momentum is fully mixed, yielding a negative shear stress Rxz and a positive compressive stress Rzz.

Vera: That positive compressive stress Rzz causes a pronounced thermal-pressure dip within the TRML.

Jocelyn: And they found that radiative losses are offset by the divergences of enthalpy flux and turbulent heat flux Qz, which is important for modeling.

Subrahmanyan: The volume-averaged temperature follows a tanh profile, resulting in universal emissivity distributions that are consistent with simulations.

Vera: The cooling-rate surface density saturates and becomes independent of box size in the strong-cooling limit, which establishes the universal structure of TRMLs.

Jocelyn: How does this connect to how we model multiphase gas dynamics in galactic environments?

Subrahmanyan: These results establish TRMLs as universal building blocks for multiphase gas dynamics, and they enable numerical calibration of the entrainment rate as a function of key parameters of the setup.

Vera: So it’s about using this universal structure to calibrate transport coefficients.

Jocelyn: And what about the time evolution part?

Subrahmanyan: At early times, a radiative shear layer evolves similarly to its adiabatic counterpart, but after roughly one cooling time t0, the system settles into a statistically steady turbulent radiative mixing layer.

Vera: They define t0 as tcool(T0 equals pTcTh) equals p0(gamma minus one)n2 over twenty lambda T0.

Jocelyn: Once steady state is established, the mean density profile is stationary in the TRML frame, and the net vertical mass flux ρuz is uniform across the layer.

Subrahmanyan: Within the TRML, rising dense plumes and sinking hot clumps generate a strong positive turbulent mass flux δρ deltauz.

Vera: The vertical momentum equation attains steady state: ddz ρuz squared plus Rzz plus p equals zero.

Jocelyn: Locally, cooling is balanced by the divergence of enthalpy flux, B times ρuz, which is pdV work, and turbulent heat flux Qt equals δB delta(ρuz).

Subrahmanyan: The mean TRML temperature profile is well described by T z equals Th minus Tc2 tanh z over z0 plus Th plus Tc2, with zero defined at T equal to (Th plus Tc) over two.

Vera: This paper, Structure of Turbulent Radiative Mixing Layers, shows that shear-driven turbulence fundamentally changes radiative cooling layers between dense and diffuse CGM phases.

Jocelyn: So the pressure dip caused by the compressive Rzz greater than zero stress allows a steady, subsonic cooling flow instead of zero stress allowing a steady, subsonic cooling flow.

Subrahmanyan: For sufficiently large shear to cooling time ratio ξ equals L perpendicular over delta ut0, the transport coefficients and the entrainment or cooling rate density Σdot cool become independent of the box size.

Vera: That independence from box size is a defining hallmark of TRMLs.

Jocelyn: This really means we can use these simulations to calibrate rates in real astrophysical systems like galactic environments.

Subrahmanyan: This research establishes TRMLs as universal building blocks for multiphase gas dynamics, consistent with the simulation results they presented.

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

Jocelyn: This work uses a machine learning convolutional neural network to find Lyman alpha emitters hidden in the Dark Energy Spectroscopic Instrument data release one dataset.

Subrahmanyan: It addresses a problem because current DESI pipelines aren't designed to measure redshifts for galaxies at redshift greater than two.

Vera: The architecture has two main parts, starting with fully-convolutional layers without any fully-connected layers, where each convolution layer is followed by a ReLU activation function and a max pooling layer to downscale the size of the input spectrum.

Jocelyn: 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 it outputs a lower resolution spectrum with seventy-six pixels and one pixel being about twenty-five arcseconds.

Subrahmanyan: Each output pixel value is between zero and one, indicating the probability of that region in the original spectrum containing a Ly alpha emission line.

Vera: Then there's the second part, which is a combination of a slicer and a fitter where the slicer identifies pixels with a probability greater than zero point five, and then extracts that region from the original spectrum.

Jocelyn: That extracted spectrum is passed to a fitter consisting of two convolution layers followed by three fully-connected layers to fit the red peak of the Ly alpha emission line with a peak-normalized Gaussian.

Subrahmanyan: The training dataset has six thousand five hundred sixty spectra with a one to one ratio between LAEs and non-LAEs, and they added one thousand two hundred thirteen quasars and one thousand two hundred thirteen low-redshift ELGs to balance the numbers.

Vera: The model is trained using the dimensionless signal-to-noise ratio spectra as input, normalizing the flux by its noise.

Jocelyn: The final model achieved a purity of ninety-five point two percent and a completeness of ninety-five point nine percent for detecting LAEs in the test dataset.

Subrahmanyan: Applying this model to about two × ten six spectra targeted as emission-line galaxies, they detected nineteen thousand six hundred eighty-five LAEs from z approximately two to three point five within twelve minutes on a single GPU.

Vera: The detected LAEs are mostly at the bright end of the luminosity function with a Ly alpha luminosity greater than or equal to ten forty-three erg per second.

Jocelyn: The paper also identifies several misidentification patterns by the DESI pipeline, Redrock, among the VI-confirmed LAEs.

Subrahmanyan: They found that twenty-three out of sixty-two LAEs, which is about thirty-seven percent, do not detect Ly alpha lines and therefore no robust redshifts are assigned.

Vera: Another issue was when twenty-one out of those sixty-two LAEs, about thirty-four percent, were mis-identified as QSOs with correct redshifts because the best-fit models showed broader Ly alpha line profiles with lower amplitude than the observed ones.

Jocelyn: Then there was another group, about twenty-one percent of them, or thirteen out of those sixty-two where the pipeline misidentified Ly alpha lines as O II doublet lines and yielded incorrect redshifts.

Subrahmanyan: Most of those misidentified lines had assigned redshifts between zero point one and zero point five.

Vera: The final catalog consists of two samples, one with nineteen thousand six hundred eighty-five LAEs selected by a delta chi squared less than twenty and another with four thousand seven hundred four from the greater than or equal to delta chi squared twenty selection.

Jocelyn: Looking at the results, the LAEs are blue galaxies at the bright end of their Ly alpha luminosity function with a luminosity greater than or equal to ten forty-three erg per second and redshift between two and three point five.

Subrahmanyan: The composite spectrum also shows rich spectral information including P-Cygni features like N v and C iv along with broad He ii emission and Fe II absorption/emission, as well as MgII emission lines at longer wavelengths not covered in previous studies.

Vera: They estimated the average velocity offset being about two hundred fifty kilometers per second between the redshifts based on Ly alpha emission lines and the systematic redshifts of galaxies.

Jocelyn: The paper also plans to improve their model by using the entire spectrum as model input to reduce false-positive rates using other spectral lines, and they will explore attention-based Transformer architectures to capture correlations between multiple emission lines like Ly alpha and CIV.

Lucky paper: 2509.09279: Vera: We're looking at Tidal Tails and Their Dynamics in Open Clusters Using Gaia DR3 today.

Jocelyn: It uses unsupervised machine learning to find and analyze tidal tails in five open clusters using Gaia DR3 data, Vera.

Subrahmanyan: The core finding is that tidal tails were present in all five clusters, extending from forty pc up to one hundred pc and containing between one hundred and two hundred stars.

Vera: That's a pretty consistent result across the sample.

Jocelyn: And how did they go about filtering the data initially?

Subrahmanyan: They started with rigorous data filtering using proper motion cuts based on cluster-specific bounds, resulting in a Filtered Sample.

Vera: Proper motion errors were kept below zero point five mas/yr, plus positive parallaxes and relative parallax errors less than zero point one.

Jocelyn: So they were being very picky about the quality of the astrometric measurements before starting the detection steps?

Subrahmanyan: Exactly, and they identified the core cluster using five normalized astrometric features: ra, dec, pmra, pmdec and distance.

Vera: Then they used a multi-step approach to actually detect those tails.

Jocelyn: What was the first step in that detection methodology?

Subrahmanyan: The first step was Radial Velocity filtering with a plus minus ten km/s window around the core cluster’s mean radial velocity as a loose filter.

Vera: That sounds like a good way to isolate the cluster members before looking for tails.

Jocelyn: After that, what did they do to identify potential tail stars?

Subrahmanyan: Then they used CMD Matching with a KD-Tree structure for efficient nearest-neighbor searches, retaining only those with dCMD less than or equal to zero point zero five.

Vera: So they were filtering based on color magnitude diagram similarity first.

Jocelyn: After that, PCA and clustering determined the primary elongation axis and aligned the coordinate frame before they used DBSCAN to separate tail stars from field noise.

Subrahmanyan: That whole process is what allowed them to analyze the tidal tail kinematics using proper motion corrections, orbital modeling, luminosity functions, binary statistics and rotational dynamics.

Vera: The kinematic findings were quite specific regarding rotation.

Jocelyn: Can you tell us about the key rotational results?

Subrahmanyan: They found significant rotation in M67 and NGC two thousand two hundred eighty-one for the first time, with M67 showing a rotational velocity of vphi equal to minus zero point eight four plus or minus zero point zero seven km/s.

Vera: And how did they estimate the mass for those clusters?

Jocelyn: They used Plummer model fits which provided mass estimates of eight hundred thirty-two to eight hundred seventy-six solar masses for M67 and seven hundred forty-one to eight hundred three solar masses for NGC two thousand two hundred eighty-one.

Subrahmanyan: Both clusters exhibited tidal tails extending one to four times their Jacobi radii with tail lengths between forty pc and one hundred pc.

Vera: That length range seems substantial for these structures.

Jocelyn: The study concluded that optimized clustering methods can detect extended structures up to a kiloparsec, which is much larger than the convergent point method which is typically limited to within approximately four hundred pc.

Subrahmanyan: They also noted that detecting rotation in outer regions can be used as a test for validating the presence of tidal tails because six-dimensional information is poor at larger distances.

Vera: So it's using the rotation data to confirm the physical presence of these tails at greater distances.

Jocelyn: The paper on Tidal Tails and Their Dynamics in Open Clusters Using Gaia DR3 gives us a good idea of how structure evolves within dense stellar systems.

Subrahmanyan: It really shows that optimized clustering methods can detect extended structures up to a kiloparsec, increasing the range of clusters that can be studied compared to other techniques, such as the convergent point method which is typically limited to within approximately four hundred pc.

Lucky paper: 2504.17506: Vera: We're looking at Supernova dust synthesis Carbonaceous dust in the very early universe today.

Jocelyn: This paper is about how molecular formation reactions affect grain formation in supernova ejecta and how mixing and fallback change the dust composition.

Subrahmanyan: The core finding is that including molecule formation reactions significantly affects grain production, specifically showing that carbon grain mass is underestimated by about forty percent compared to models where both are solved self-consistently.

Vera: That's a big correction for how we calculate the dust content in these explosions.

Jocelyn: They developed a state-of-the-art dust synthesis code that solves molecule and dust formation self-consistently within the supernova ejecta, incorporating hydrodynamics and radiative transfer as well.

Subrahmanyan: The model starts by calculating preSN stars using Umeda et al two thousand with their stellar evolution code, then proceeds to explosive nucleosynthesis from Tominaga et al two thousand seven to get the radial distribution of elemental abundances.

Vera: So they are tying the initial stellar state directly into the subsequent dust formation calculations.

Jocelyn: They also looked at how mixing and fallback influence this composition, considering two models for C-normal and C-enhanced metal-poor stars.

Subrahmanyan: In the faint SN model M25E10F, while the total dust amount is about eighty-five percent smaller, the carbon grain mass is only about twenty percent smaller than in the hypernova model.

Vera: And that leads to a carbonaceous to silicate grains ratio of about zero point three in that faint SN case.

Jocelyn: The conclusion drawn from this is that faint supernovae can be a significant carbonaceous dust factory in the early universe, and both faint SNe and faint hypernovae can cause the UV bump we see.

Subrahmanyan: They found that the high mass ratio of C grains to silicates in the faint SN model is about seven which they suggest might explain the observed spectral feature of a galaxy called JADES-GS-z6-zero.

Vera: So this work is presented as a very first step to uncover the cosmic dust formation history.

Jocelyn: They also mentioned that higher MPopIII means the dust mass is larger because the synthesized metal mass itself is larger, and they suggest we might need to consider dust destruction from reverse shocks propagating in ejecta.

Subrahmanyan: The total dust mass produced in model M25E10 is about one solar mass of dust, and the dominant species identified are metallic Si, Mg2SiO4, and SiO2.

Vera: The results show that the C grain production isn't significantly changed because the C layer at Mr around six solar masses isn't affected by mixing or fallback in those specific scenarios.

Jocelyn: The time evolution of gas temperature and chemical species abundance at an enclosed mass of Mr equals six point two solar masses was also analyzed, showing that the M25E10mol.fix curve overlapped with the M25E10 model because they used the same SN model.

Subrahmanyan: The paper also tested a lower explosion energy case with ESN equals one B, and found that the carbon grain mass for M25E1F is almost unchanged from M25E.

Vera: This research is important because it provides a concrete mechanism for how carbonaceous dust forms in the early universe via faint supernovae.

Jocelyn: The study was funded by the UK Science and Technology Facility Council through grants ST/R000905/one ST/V000632/one and ST/Y001443/one.

Subrahmanyan: And they used matplotlib from Hunter twenty-five for the software in this work.