Daily Summary for 2026-10-22

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

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

Jocelyn: 11 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 VAPOLA - A multiyear, multiband polarization survey of AGNs and Sgr A* at millimeter wavelengths with ALMA II. Spectropolarimetric properties and their evolution from 2017 to 2023.

The summary: Vera: It is the twenty-second of October, twenty twenty six.

Jocelyn: The VAPOLA survey analyzed spectropolarimetric properties over six years.

Subrahmanyan: This helps understand polarization evolution in active galactic nuclei and Sagittarius A star observations.

Vera: Total intensity stays stable, but linear polarization fraction and electric vector position angle show strong day-to-day changes.

Jocelyn: This variability is often linked to flaring activity in sources like M87 and Sgr A star.

Subrahmanyan: This variability is often linked to flaring activity as seen in VAPOLA.

Vera: Further investigation into jet physics showed Faraday rotation measure increases with observing frequency for all sources with reliable multi-band measurements.

Jocelyn: This suggests a magnetized sheath around the jet, possibly originating in the accretion flow for Sgr A star.

Subrahmanyan: This connects to VAPOLA because it probes the magnetic environment surrounding these central objects.

Vera: Researchers are developing an effective field theory for Lyman alpha forest fluctuations at the one-loop order.

Jocelyn: This models how small density variations affect the power spectrum of light in the early universe.

Subrahmanyan: This expansion matches previous bottom-up models and fits hydrodynamic simulation data with sub percent accuracy up to k= 3 h Mpc-1 at z= 2.8, as described in Don't miss the forest for the trees: the Lyman alpha forest power spectrum in effective field theory.

Vera: Primordial black holes of asteroid mass are viable dark matter candidates.

Jocelyn: Their capture depends on host star properties and three-body interactions.

Subrahmanyan: This framework combines analytic calculations with simulations to predict outcomes.

Vera: That detail is in The Life and Death of Stars That Capture Primordial Black Holes.

Jocelyn: Researchers derived an empirical constraint on dust attenuation law using JWST spectroscopy.

Subrahmanyan: The ensemble-average behavior aligns with local relations despite UV slope differences.

Vera: That result comes from An empirical derivation of the average dust attenuation law at 2z<7.

Jocelyn: J1244-LyC1 is the first confirmed high-redshift LyC-leaking merger.

Subrahmanyan: This shows galaxy mergers create conditions where ionizing radiation escapes more easily.

Vera: That means when galaxies collide, ionizing radiation escapes more easily than thought.

Vera: This discovery links galaxy mergers to reionization through early universe ionization history mechanisms.

Jocelyn: The MAMMOTH-LyC survey found J1244-LyC1, a system with 10 point 02 solar masses and a significant six sigma LyC detection.

Vera: Its spatial resolution resolved the emission into multiple peaks matching the disturbed shape, confirming multi-site leakage.

Jocelyn: This confirms genuine leakage, building on Lagrangian perturbation theory showing nonlinearities lead to force screening at large scales.

Vera: MAMMOTH-LyC suggests tidal interactions disrupt the interstellar medium in merging systems creating low-density escape routes.

Jocelyn: Observations of our galaxy's ISM, from pulsar scintillation arc curvature measurements, reveal rich plasma density fluctuation structures.

Vera: VAPOLA is a multiyear survey analyzing AGN polarization properties over several years to see how they change.

Jocelyn: Don't miss the forest for the trees develops an effective field theory to describe fluctuations in the Lyman alpha forest.

Vera: The Life and Death of Stars That Capture Primordial Black Holes models how stars can capture primordial black holes and evolve these systems.

Jocelyn: An empirical derivation of the average dust attenuation law at 2<z<7 used JWST data to determine how dust affects light from massive star-forming galaxies.

Vera: Probing Primordial Black Hole Mergers in Clusters with Pulsar Timing Data searches for evidence that primordial black holes cause the observed gravitational wave background.

Jocelyn: Formation and Evolution of Wolf-Rayet Planetary Nebulae through a Late Thermal Pulse simulates how a late thermal pulse causes a star to form specific planetary nebulae.

Vera: Silicate cosmic dust grain collisions in the interstellar medium is a molecular dynamics study predicting how dust grains shatter when they collide.

Jocelyn: Unveiling dark forces with measurements of the Large Scale Structure of the Universe studies how long-range dark forces affect galaxy distribution to constrain dark matter.

Vera: Lagrangian perturbation theory for modified gravity creates a mathematical method to calculate gravity behavior in different modified gravity theories on large scales.

Jocelyn: Precision Determination of Scintillation Screen Parameters from Annual Modulation Measurement of Pulsar Scintillation Arc Curvature with the FAST Telescope uses pulsar timing data to measure plasma screens.

Vera: These papers are VAPOLA for AGN polarization evolution, Don't miss the forest for the trees for Lyman alpha forest power spectrum, The Life and Death of Stars That Capture Primordial Black Holes, An empirical derivation of the average dust attenuation law at 2<z<7, Probing Primordial Black Hole Mergers in Clusters with Pulsar Timing Data, Formation and Evolution of Wolf-Rayet Planetary Nebulae through a Late Thermal Pulse, Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study.

Vera: Agreed. The show is over.

Jocelyn: Indeed. This is the end of today's review.

Lucky paper: 2608.14900: Vera: We’re looking at VAPOLA today, a multiyear survey of AGNs and Sagittarius A* using ALMA II to track spectropolarimetric properties from two thousand seventeen to two thousand twenty-three.

Jocelyn: So we know the total intensity stays stable on weekly timescales, but the polarization fraction and electric vector position angle show strong day-to-day changes.

Subrahmanyan: That variability is often linked to flaring activity, as seen in several objects examined by VAPOLA. The paper looks at how these properties evolve across different source classes too.

Vera: We see that BLLac objects are systematically more polarized than FSRQs at millimeter and sub-millimeter wavelengths in this study.

Jocelyn: That’s an important distinction, because it tells us something about the underlying physics of the jets themselves.

Subrahmanyan: The study provides strong evidence regarding the origin of Faraday rotation measure, RM. A key result is that the magnitude of RM increases with observing frequency for all sources with reliable multi-band measurements available.

Vera: That frequency dependence points toward Faraday rotation arising from a magnetized sheath surrounding the relativistic jet.

Jocelyn: And polarization fraction also increases with observing frequency, which reinforces that idea about magnetic fields in those regions.

Subrahmanyan: The relatively large RM values observed, approximately ten four to ten five rad m-two strongly suggest ordered magnetic field structures in close proximity to the central engine.

Vera: But the authors are careful and they mention that the possibility of RM originating in an accretion flow for Sgr A*, where a collimated jet isn't fully established yet, remains a viable alternative hypothesis.

Jocelyn: That’s fair, because Sgr A* is such a complex system. The VAPOLA analysis successfully maps the evolution of these parameters for key sources over several years.

Lucky paper: 2309.10133: Vera: We are looking at Don't miss the forest for the trees today, Jocelyn. This paper deals with modeling fluctuations in the Lyman alpha forest using an effective field theory framework.

Jocelyn: It’s about taking those small density variations and describing how they affect the power spectrum of light in the early universe. The core idea is establishing a self-consistent, one-loop order EFT expansion for these fluctuations.

Subrahmanyan: The authors prove this by showing mathematical identity between two methods: a "bottom-up" expansion based on transmitted flux and a "top-down" derivation using the exponential map of the optical depth field.

Vera: They also point out inconsistencies in tree-level perturbative expansions when counterterms are omitted, which is something they are fixing by introducing a generalization of the FFTLog method.

Jocelyn: That generalization lets them efficiently calculate one-loop integrals involving new line-of-sight dependent operators, and then they rigorously test this resulting one-loop EFT model against hydrodynamic simulations.

Subrahmanyan: The paper demonstrates exceptional predictive power, fitting the simulation data with sub percent accuracy across both three dee and 1D correlation analyses up to a wavenumber of k=three h Mpc-one at a redshift of z=two point eight.

Vera: So, this means we have a robust one-loop EFT model that is validated by internal consistency checks and tested against actual simulation data.

Jocelyn: It gives cosmologists a practical tool for full-shape analyses of Lyman alpha forest data, which was the main goal when they developed this model.

Subrahmanyan: The implication is that we can now use this EFT expansion effectively for modeling those early universe density fluctuations at high redshift.

Lucky paper: 2412.15989: Vera: We're moving into segment five now, and we're looking at a paper titled Probing Primordial Black Hole Mergers in Clusters with Pulsar Timing Data.

Jocelyn: This paper dives into whether the stochastic gravitational wave background suggested by Pulsar Timing Array datasets can compete with scalar induced GWs.

Subrahmanyan: The analysis finds that scalar-induced GWs dominate in a broad primordial black hole mass distribution, although specific models with enhanced clustering could yield a background dominated by PBH mergers.

Vera: They derive the PBH mass function from a power law primordial spectrum of curvature perturbations using the Press-Schechter formalism.

Jocelyn: The mass is approximated by the horizon mass at re-entry, given by MH(k) equals four pi M2 p over H(one) equal to twenty solar masses times k106 Mpc minus one g four star s of Tcr g minus three star of Tcr.

Subrahmanyan: The PBH mass function fPBH(m) is then expressed in terms of the primordial power spectrum P zeta (k) and the critical threshold delta c.

Vera: For a broad mass distribution obtained for values of ns equals one the behavior of the equation of state w(T) with temperature, especially during the QCD epoch, can play an important role.

Jocelyn: They also look at scalar induced gravitational waves which are induced by scalar perturbations.

Subrahmanyan: For Gaussian perturbations, this induced GW spectrum can be sizeable for enhanced primordial curvature perturbations at small scales around P zeta equals ten minus two.

Vera: The observable induced GW spectrum today is defined as omega g w zero equal to one over rho c times d rho zero g w over d log f.

Jocelyn: This can be expressed in terms of the induced GW spectrum in the radiation era through equations twenty-four eighty-six.

Subrahmanyan: For the spectrum defined in equation (three), one can derive an approximate analytical expression for the induced GW spectrum today given by omega SIGW g w zero equals seven point five times ten minus five g star of T g star s of T minus four over three Q ns f f star squared ns minus one theta of f minus f min theta of f max minus f.

Vera: Then there's the GW background from PBH binaries considered through two main channels, early universe binaries and late-universe binaries that can dynamically form in clusters via capture.

Jocelyn: The merger rate from both these channels has many uncertainties as discussed in section IV.

Subrahmanyan: For the SGWB from unresolved PBH binaries, the spectrum follows the standard relationship to the comoving number density of compact binary sources.

Vera: The differential GW spectrum emitted in equation (twenty-five) can be expressed by plugging equations twenty-four and fourteen into equation twenty-five giving omega PBHB g w zero equals f two over three pi G squared over three rho c z d log m1 d log m2 M five over three c with V equal to V early plus V late.

Jocelyn: The total GWB is given by summing equations thirteen and twenty-nine meaning h two omega g w equals h two omega SIGW g w plus h two omega PBHB g w.

Subrahmanyan: In the fixed clustering analysis, the PTA signal is always dominated by the scalar induced GWB with the PBH binaries giving a much sub-leading contribution to the total GW spectrum.

Vera: To explain that common-spectrum process, a spectral index ns equals one point three would yield the required increasing frequency dependence of the GW spectrum from SIGWs, omega SIGW g w proportional to f two times ns minus one.

Jocelyn: The analysis imposes the CMB µ-distortion constraints on the curvature power spectrum at scales k less than or equal to one hundred five Mpc minus one.

Subrahmanyan: Since gravitational waves act as radiation, they are thus constrained by current bounds on the effective number of neutrino species limit, delta N eff less than zero point twenty-eight.

Vera: The LVK constraint imposes that at f LVK equals twenty-five Hz, omega g w is less than or equal to one point seven times ten minus eight.

Jocelyn: The posterior values for the fixed clustering analysis in figure eight do not violate this constraint.

Subrahmanyan: The main results from the Bayesian inference show that the contribution to the GWB from PBH binaries is largely subdominant with respect to the GWB from scalar induced perturbations for power law PBH mass distributions.

Vera: It is not possible to explain the PTA signal with only PBH binaries in late time PBH clusters under standard PBH formation scenario with Poissonian fluctuations.

Jocelyn: The analysis reveals that the required clustering factor Rcl to explain the PTA signal with late PBH binaries is about one hundred four Gpc minus three year minus one larger than the Rcl values arising from Poissonian fluctuations.

Subrahmanyan: They also show that PBH mass distributions with ns less than or equal to one are statistically preferred because they naturally suppress low-mass PBHs from overclosing the Universe and do not require a small scale cut-off k max.

Vera: The resulting maximum likelihood GW spectrum for this analysis is shown in figure five right where the total GW background is a combination from scalar induced shown in blue solid till the cut-off k max of f max equals one hundred seven Mpc minus one, which is ten minus eight hertz.

Jocelyn: The best fit values used for the GW spectrum shown in figure five are in violation of this LVK bound.

Subrahmanyan: The analysis uses a broad PBH mass function including features from the QCD epoch as detailed in section II.

Lucky paper: 2603.26536: Vera: Let's move on to our fourth paper of five, which is Formation and Evolution of Wolf-Rayet Planetary Nebulae through a Late Thermal Pulse.

Jocelyn: This paper looks at the first radiation-hydrodynamical simulations for a born-again planetary nebula triggered by a late thermal pulse.

Subrahmanyan: The research uses stellar evolution models, specifically with the mesa code, to simulate a low-mass star experiencing this born-again event.

Vera: They activated updated opacity tables for H-deficient, C-rich material used for R Coronae Borealis stars and adopted a massloss prescription tailored for H-deficient WR winds.

Jocelyn: The simulations use the pluto code with a finite-volume formulation covering a radial extent of zero point six pc, starting from an initial medium density of nISM = one cm−three and temperature of TISM = one hundred K.

Subrahmanyan: The simulations show that after roughly one thousand yr of post-LTP evolution, the H-deficient ejecta will mix with the clumps and filaments formed in the innermost nebular rim.

Vera: The total mass ejected by the born-again event is M LTP = three point four × ten−four solar masses, which is relatively small compared to the ionized nebula mass.

Jocelyn: That mixing means once it happens, you can't really tell if the H-deficient ejecta are even there anymore.

Subrahmanyan: This explains why PNe hosting WR CSPN don't show significant abundance differences when compared with PNe with H-rich CSPN, both cases looking consistent with models of H-rich post-AGB stars.

Vera: The high degree of turbulent mixing can also explain why WR PNe are more easily detected by X-ray instruments.

Jocelyn: A synthetic optical spectrum computed four thousand five hundred yr after the LTP event predicts that the H-deficient CSPN will evolve extremely fast into a WO-type star.

Subrahmanyan: The model predicts a carbon-to-oxygen abundance ratio of X C/X O greater than or equal to one point zero at late stages of the LTP evolution, which should lead to large amounts of C-rich dust during this short phase.

Vera: The nitrogen abundance in NGC five thousand one hundred eighty-nine stands out because its XN is zero point zero one compared to other WR CSPN values where XN is less than or equal to three × ten−three which aligns with LTP models.

Jocelyn: This difference in nitrogen abundance seems to suggest that XN can be used to tell the difference between an LTP and a VLTP event.

Subrahmanyan: Our stellar evolution model predicts abundances slightly different from previous LTP models, but these new predictions are very similar to those of WR CSPN reported in the literature.

Vera: Overall, our model reproduces well the observed abundances of H-deficient WR CSPN.

Jocelyn: The radiation-hydrodynamical simulations show that a star forming a PN and undergoing an LTP event develops a short-lived double-shell structure.

Subrahmanyan: Our model indicates that after roughly one thousand yr of post-LTP evolution, the H-deficient ejecta will mix with the clumps and filaments formed in the innermost nebular rim.

Vera: The LTP event produces a layer of enhanced turbulence inside the H-rich nebula, which is consistent with previous evidence that WR PNe have larger expansion velocities and turbulence.

Jocelyn: The combination of born-again ejecta and high mechanical luminosity from the post-LTP wind produces a high degree of turbulence.

Subrahmanyan: We compare our LTP model predictions with a grid of born-again stellar models, including both LTP and VLTP events, to probe the origins of PNe surrounding WR CSPN.

Vera: The paper confirms that WR CSPNe can be produced through a born-again event based on these results.

Jocelyn: This work is Formation and Evolution of Wolf-Rayet Planetary Nebulae through a Late Thermal Pulse.

Lucky paper: 2605.22225: Vera: We're coming to paper number five of five which is Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study.

Jocelyn: This paper looks at how silicate grains shatter and vaporize when they collide in space using molecular dynamics simulations.

Subrahmanyan: The simulation methodology uses classical molecular dynamics to explore collisions between silicate grains with radii a from five Å up to fifty Å at velocities between zero point one km/s and twenty km/s.

Vera: They used two materials: pure amorphous silica SiO2 and an amorphous silicate of composition Mg1 point 3Fe0 point 3Si1O3 point 6, which they inferred from the astrodust model by Draine and Hensley in two thousand twenty-one.

Jocelyn: The chemical bond formation and breaking are modeled with the ReaxFF force field, and for the geometry, frequency, noncovalent methods they use GFN1-xTB.

Subrahmanyan: For the SiO2 grains, they start with spherical grains extracted from a cubic silicate bulk in the AMS molecular builder and run MD simulations to minimize energy up to two thousand K before relaxing down to ten K for an amorphous structure.

Vera: They estimated radii of about six point five, sixteen point five, thirty-five and fifty Å for the SiO2 grains, and structures for the astrodust grains were constructed using the multi-component artificial force-induced reaction method by Sameera et al. in two thousand sixteen with radii of five and seven point five Å respectively.

Jocelyn: Looking at collision outcomes as a function of velocity, they found that SiO2 grains coagulate for all collision velocities lower than about five km/s, without any bouncing behavior observed in their simulations.

Subrahmanyan: There is a rapid transition between five and seven km/s where the combined shattered and vaporized fractions increase to approximately ninety percent.

Vera: For increasingly higher velocities, most of the colliding grains are disrupted across all cases.

Jocelyn: In all cases, any reasonable choice for a threshold shattering velocity for either composition and any grain size is greater than the canonical value of two point seven km/s from Jones et al. in one thousand nine hundred ninety-six by a factor of approximately two.

Subrahmanyan: The size distributions of shattering products are generally not well-described by a single power law, displaying a complicated dependence on collision velocity, and they appear to become steeper with increased collision velocity.

Vera: The shapes of the grains themselves affect their optical and dynamical properties, and at low velocities, the resulting coagulated grain has a geometrically defined size of about ninety Å.

Jocelyn: That means it suggests the grain is composed of two spheres that largely retain their shape but stick together with some overlap, which they call a prolate grain.

Subrahmanyan: Higher collision velocities produce smaller grains, suggesting the constituent grains stick together with larger overlap in those cases.

Vera: Theoretically, the corrected Tielens et al. one thousand nine hundred ninety-four expression predicts a velocity threshold of seven point nine km/s for silicate materials, which they say is in moderately better agreement with their findings for SiO2 when compared to velocities at which half the grain is disrupted.

Jocelyn: However, the model by Hirashita and Kobayashi in two thousand thirteen fails to match these simulations even when modified to best match the correct threshold velocity.

Subrahmanyan: The new, higher velocity thresholds they obtain suggest that interstellar dust composed of silicate material is more resilient to shattering and vaporization than previously assumed.

Vera: They recommend adopting Equation five for the shattering threshold, which includes a step function in disrupted mass fraction as a function of velocity because the rapid transition shown in Figures eleven and twelve is justified.

Jocelyn: Their results also suggest that shattering and vaporization only need to be considered for collisions between grains with mass ratios less than ten.

Subrahmanyan: Further simulations with a larger range of grain sizes and collision parameters will be needed to understand how grain collision outcomes transition to the Tielens et al. one thousand nine hundred ninety-four results for grains large enough to be modeled with continuum dynamics.

Vera: We find that grain shattering thresholds for the silicate materials we simulate are about six km/s, which is a factor of approximately two higher than the canonical value of two point seven km/s from Jones et al. in one thousand nine hundred ninety-six.

Jocelyn: This discrepancy can be explained in part by correcting an error in the expression for the velocity thresholds defined by grain material properties derived in Tielens et al. one thousand nine hundred ninety-four and used for the numerical values quoted in Jones et al. one thousand nine hundred ninety-six.

Subrahmanyan: The prescription from Tielens et al. one thousand nine hundred ninety-four for the fraction of grain mass shocked to a given pressure, Equation seven as a function of collision velocity fails to predict the fraction of shattered or vaporized material observed in their numerical simulations.

Vera: We provide updated shattering velocity thresholds for the most likely candidate grain materials in Table two for use in future models of the dust population evolution on astrophysical scales.

Jocelyn: Indeed, this is the conclusion of Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study.