Daily Summary for 2026-10-23

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

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

Jocelyn: 18 new papers came out today.

Vera: 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 Stardust Galaxies at z>9.

The day's research: Vera: Burgarella et al.'s AKMPC investigates dust production in early galaxies.

Jocelyn: It suggests supernova remnants make dust before grain growth starts.

Vera: This shift explains UV-bright galaxy abundance at high redshifts.

Subrahmanyan: The PAHs work models polycyclic aromatic hydrocarbon lifecycles.

Jocelyn: Grain-grain shattering in the medium drives PAH abundance evolution.

Vera: That links metallicity to dust growth systematically across galaxies.

Subrahmanyan: Searching for signatures of self-interacting dark matter alters cluster density profiles.

Jocelyn: Baryonic and projection effects usually reduce those deviations significantly.

Vera: Calibrating optical galaxy cluster projection effects corrects line-of-sight contamination.

Subrahmanyan: This technique quantifies projection effects needed for accurate cluster measurements.

Vera: Io's tidal response links interior, partial melting, and Juno observations quantitatively.

Jocelyn: It models Io as three layers and finds heating when melt fractions are below the critical threshold.

Subrahmanyan: This suggests a heterogeneous mantle structure, pointing to a magmatic sponge.

Vera: Mass flux analysis confirms melt percolation exceeds thermodynamic production, implying efficient drainage.

Jocelyn: That contrasts with jet power estimates for Flat-Spectrum Radio Quasars being insensitive to particle energy distribution.

Subrahmanyan: The analysis shows melt presence decreases viscosity and increases anelasticity in the upper mantle.

Vera: This amplifies upper-mantle dissipation, connecting to modeling X-ray reflection spectra from disk returning radiation.

Jocelyn: That suggests returning radiation from a disk could explain reflection features in black-hole X-ray binaries.

Subrahmanyan: We see connections between turbulence and structure in predicting prestellar cores using machine learning.

Vera: That work highlights the need for spatial architectures beyond point-wise models to capture differential operators.

Jocelyn: It uses machine learning to forecast core collapse based on local density and magnetic field magnitude.

Lucky paper: 2605.09829: Vera: Let's talk about Stardust Galaxies at z>nine and how it explains the excess of UV-bright galaxies we see at those high redshifts.

Jocelyn: So, the core idea here is that before efficient interstellar medium grain growth kicks in, dust is primarily made by supernova remnants.

Subrahmanyan: That initial regime means the surviving dust population is biased toward larger grains, which gives it a much lower intrinsic far-UV dust opacity.

Vera: It says this intrinsic far-UV dust opacity is about ten to the power of-four cm squared g-one which is way lower than standard interstellar medium values.

Jocelyn: And when metallicity reaches a critical value, say zero point one Z, that's when metal accretion onto grain surfaces becomes efficient, leading to ISM grain growth.

Subrahmanyan: This defines a hybrid model where the dust mass-to-stellar mass relation is blended between the SN-dominated branch and the ISM grain-growth branch.

Vera: The transition boundaries are set at x lo about seven point eight and x hi about nine point zero, with a gas-fraction transition midpoint at x b around eight point three nine.

Jocelyn: They use mixed radiative-transfer geometries to model photon escape, finding the best fit is achieved at the lower search boundary for the active-dust fraction tau lambda.

Subrahmanyan: The effective transport opacity kappa eff is calculated as zero point three five and g FUV is zero point six zero in their fiducial calculation, which leads to a far-UV attenuation of AFUV of zero point three five.

Vera: That low intrinsic opacity is what allows them to match the JWST measurements that require modest dust corrections for their theoretical predictions to align with observed UV luminosity functions.

Jocelyn: The paper also incorporates a sophisticated mass-metallicity relation that distinguishes between those low-metallicity, SN-dominated regimes and the high-metallicity ISM growth regimes.

Subrahmanyan: They define the Dust-to-Metal ratio using an empirical broken DTM–metallicity relation, suggesting lower DTM ratios in the early universe.

Lucky paper: 2508.15337: Vera: We're looking at the paper on Hierarchical Black Hole Mergers in AGN Disks, which is a pretty dense model by Lee et al.

Jocelyn: This work looks at how dense environments around active galactic nuclei can trigger repeated mergers of stellar-mass black holes.

Vera: It uses a framework based on the fastcluster code extended with AGN specific prescriptions to describe this process.

Subrahmanyan: The model describes the AGN disk using the steady-state model of Sirko and Goodman, which gives us radial profiles for gas surface density sigma g and aspect ratio h.

Jocelyn: So they are numerically computing those profiles for given values of the supermassive black hole mass M• and the disk viscosity parameter alpha.

Vera: They start with stellar-origin black holes formed in a nuclear star cluster at solar metallicity using the SEvN code, distributed uniformly in radius.

Subrahmanyan: These black holes are captured by the disk via gas drag, and once embedded, they undergo migration toward regions of higher numerical density where binary formation is favored.

Jocelyn: After formation, these binaries evolve through three main processes: gas hardening, three-body encounters with single black holes, and gravitational-wave emission driving the final inspiral.

Vera: They compute remnant mass mrem and spin chi rem from numerical relativity fits to see if they can undergo further capture and hierarchical mergers.

Subrahmanyan: For those hierarchical mergers to happen, three conditions must be met: the recoil velocity vkick has to be below the local escape velocity vesc in ninety-nine point four percent of cases, the cumulative accreted mass Macc does not exceed the total captured black hole mass in seventy-two point nine percent, and the remnant formation timescale is shorter than tau at forty-one point three percent.

Jocelyn: If those conditions are met, this process creates a chain of progressively more massive binary black hole mergers with low mass ratios q less than or equal to zero point one and effective spins spanning extreme values up to chi eff of plus or minus zero point eight.

Vera: The inferred fraction of AGN contribution to observed binary black hole mergers has a median value around eight percent, with a credible interval from less than one percent up to about twenty percent.

Subrahmanyan: They also look at how the redshift evolution of AGN-driven merger rates can inform SMBH growth across cosmic time.

Lucky paper: 2608.21963: Vera: Now we're looking at a paper on searching for signatures of self-interacting dark matter in halos from full physics simulations.

Jocelyn: It dives into how this self-interacting dark matter, or SIDM, affects the structure of galaxy groups and clusters when you look ahead to future optical surveys.

Subrahmanyan: The study focuses on two types of SIDM models: rare SIDM and frequent SIDM, defined by how the scattering cross-section per unit mass changes with relative velocity.

Vera: They set the normalization for this transfer cross-section as sigma T zero over mDM equals forty cm2 g-one and the transition velocity vc is at two hundred km s-.

Jocelyn: The simulation setup used Dianoga one times SIDM simulations, zooming in on regions of comoving length one Gpc with cosmological parameters omegam = zero point two four and ns = zero point nine six.

Subrahmanyan: They re-simulated six zoom-in regions, using both a CDM model and both the rare and frequent SIDM models, including baryons in the full-physics runs.

Vera: In lower-mass halos, SIDM effectively produces more cored internal three dee profiles than CDM with relative differences of about twenty percent in the DMO case.

Jocelyn: The full-physics case results show smaller deviations because baryons counteract those effects from the SIDM.

Subrahmanyan: Both DMO and FP SIDM profiles also have a mass excess of about five percent at intermediate radii, which they attribute to heat transfer from collisions.

Vera: The most massive SIDM halos are more affected by baryons, leading to cuspier inner profiles than CDM and a much less pronounced mass excess at larger radii.

Jocelyn: They found that projecting the mass and density profiles along the line of sight alleviates the effects of SIDM at all halo masses, which leads to cuspier profiles than the three dee case.

Subrahmanyan: But weak lensing shear profiles still reach deviations up to about twenty percent above resolution limit for M200 halos around one thousand thirteen solar masses, across all lines of sight.

Vera: That means weak lensing inference of surface density excess profiles could prove particularly useful for detecting SIDM effects in those larger clusters.

Jocelyn: The median NFW concentration remains only slightly affected by SIDM at all mass scales, with deviations at the ≲ five percent level in the projected full-physics case.

Subrahmanyan: The study concludes that while SIDM imprints on halo mass and density profiles are washed out from both baryons and two-dimensional projection, there are certain observationally promising radial scales that may act as robust probes through weak lensing.

Vera: They point to the internal regions of group-sized halos and the cores of very massive clusters as places where this detection might be possible.

Jocelyn: The future work will focus on studying the halo velocity dispersion profiles to see if that can represent a viable probe for detecting SIDM.

Lucky paper: 2603.25331: Vera: We’re moving into a paper about jet power estimates from broadband SED modeling of Flat Spectrum Radio Quasars PKS one thousand four hundred forty-one plustwenty-five and Ton five hundred ninety-nine.

Jocelyn: That paper looked at how jet power estimates change depending on the assumed particle energy distribution.

Vera: They tested four models: a broken power law, a log-parabola, and two energy-dependent models for diffusion or acceleration timescales.

Jocelyn: The key finding is that the jet power estimates from models with intrinsic curvature are of the same order as those from the broken power law distribution.

Vera: This contrasts sharply with High Synchrotron Peaked Blazars where estimates can differ by nearly two orders of magnitude between models.

Subrahmanyan: They attribute that difference to lower electron break energies typically seen in Flat Spectrum Radio Quasars compared to HBLs.

Jocelyn: So, unlike in HBLs, the estimated jet powers for FSRQs seem relatively stable regardless of the particle energy distribution assumed.

Vera: For PKS one thousand four hundred forty-one plustwenty-five with a Lorentz factor Gamma equal to twenty and a gamma minimum of ten the power was around ten forty-six erg s-one.

Jocelyn: And for Ton five hundred ninety-nine all models showed jet powers in the range of about ten forty-five point five to ten forty-six point five erg s-one.

Subrahmanyan: The calculated proton power component reported in Table three and Table four was found to be less than that of HBL like Mkn five hundred one.

Vera: That suggests the proton power didn't seem to dominate in FSRQs, which might explain why the jet power isn't sensitive to gamma minimum variation.

Jocelyn: The researchers analyzed the variation of chi squared with gamma minimum to ensure physically acceptable solutions.

Subrahmanyan: Although X-ray spectra could fit with either SSC or EC components, EC fits yielded unphysical gamma minimum and equipartition values, so they kept only SSC fits for the X-ray band.

Vera: They concluded that external Compton processes and weaker dependence on spectral curvature dominate the estimation of jet power in these sources.

Lucky paper: 2602.16043: Vera: Let's talk about the paper on Diffusive Instabilities in Dusty Disks, which is paper five of five this show.

Jocelyn: The paper revisits the diffusive instability in dusty disks when dust mass diffusivity or viscosity decreases steeply with increasing dust density.

Vera: The basic criterion for diffusion-slope-driven growth is βdiff less than or equal to minus two for small dust stopping times.

Subrahmanyan: This is driven by the density dependence of the transport coefficients themselves, not just background gradients.

Jocelyn: They linearize around a constant background state and get a fourth-order dispersion relation n4 plus a3n3 plus a2n2 plus a1n plus a0 equals zero.

Vera: A sufficient criterion for an unstable static mode is that the coefficient of the leading term, which is n0, must be less than zero.

Subrahmanyan: This immediately implies that βdiff less than minus two and βvisc less than one are approximate criteria for this instability.

Jocelyn: The growth rates for the fastest growing mode are analyzed using quantities like tau s equal to ts omega and sigma equal to n over omega.

Vera: The analysis confirms that beta diff less than or equal to minus two leads to an instability when sigma is less than one.

Subrahmanyan: They then look at nonlinear evolution numerically and find that for power-law closures D proportional to sigma beta diff with beta diff less than zero, the nonlinear evolution does not saturate.

Jocelyn: Instead, steepening gradients amplify the nonlinear dust-pressure term and drive finite-time collapse into increasingly sharp spikes.

Vera: To handle this, they test a simple piecewise closure where the negative diffusion slope operates only over a finite density interval to eliminate blowup.

Subrahmanyan: This produces peak densities controlled by an imposed saturation scale sigma sat, but they admit their model lacks a radial gas pressure gradient which is usually in simulations.

Jocelyn: The model doesn't reproduce multidimensional streaming-instability simulations because it misses that background equilibrium drift.

Vera: They propose a piecewise diffusion law D equal to D0 times fpw sigma with fpw defined piecewise to enforce saturation at high dust densities.

Subrahmanyan: This modification eliminates blowup and produces narrow filaments whose peak densities are controlled by the saturation scale sigma sat.

Jocelyn: The authors support the robustness of diffusive instabilities as a linear organizing mechanism whenever dust diffusion decreases sufficiently steeply with density.

Vera: But they stress that realistic nonlinear outcomes depend sensitively on how turbulent diffusion behaves at high dust loading and on multidimensional saturation effects not captured by their one-dimensional framework.

Subrahmanyan: Future work should focus on calibrating density-dependent transport closures from streaming instability simulations and extending the model to include a background pressure gradient.