Papers for Thursday, Feb 26 2026

We present a deep-learning-based approach for identifying dark matter haloes in cosmological N-body simulations. Our framework consists of a volumetric Convolutional Neural Network to classify individual simulation particles as either halo or non-halo members, followed by a highly optimised and parallelised Friends-of-Friends clustering algorithm that groups the classified halo members into distinct haloes. The training data comprise simulations generated using GADGET-4, with labels obtained with the ROCKSTAR halo finder. Our models incorporate two main halo mass definitions, $M_{200\mathrm{b}}$ and $M_{\text{vir}}$, with similar performance. For haloes defined by the ROCKSTAR $M_{200\mathrm{b}}$ criterion, the classification network demonstrated stable performance across multiple simulation resolutions. For the highest resolution, it achieved over $98\%$ across all primary performance metrics when identifying halo particles. Furthermore, the FoF algorithm yielded halo catalogues with a purity generally exceeding $95\%$ and a stable completeness of $93\%$ for masses above $5\times10^{11} \, M_\odot$. Our pipeline recovered the centre-of-mass positions, velocities and halo masses with high fidelity, yielding a halo mass function consistent to within $5\%$ of the reference while faithfully reconstructing the internal density profiles. The primary objective of this study is to offer a faster and scalable alternative to conventional halo finders, achieving a speed-up of approximately one order of magnitude relative to ROCKSTAR, offering a promising pathway for modern simulation-based inference methods that rely on rapid and accurate structure identification.

Crater II (CraII), a large and low-density dwarf spheroidal galaxy, has unusual observed properties that are difficult to reproduce in cold dark matter simulations. Ongoing tidal disruption may help explain the discrepancies, as evidenced by the recent discovery of tidal tails. Here we present metallicity-sensitive narrowband photometry of the Ca II H and K lines from the Dark Energy Camera, covering $128$ deg$^2$ across the center and identified tidal tails of CraII as part of the Mapping the Ancient Galaxy in CaHK (MAGIC) survey. Our combined photometric metallicity, color-magnitude, proper motion, and parallax selections identify 162 CraII candidates. Of these, 37 candidates are located in the tidal tails which extend at least $7^\circ$ ($\sim 95$ kpc) from the center of CraII, suggesting it has lost $\gtrsim 25$% of its initial stellar mass. We confirm low contamination rates with dedicated control fields and highlight the extremely low surface brightness stellar features that can be uncovered with CaHK data, as faint as $\sim 36$ mag arcsec$^{-2}$. We also make the first detection of a metallicity gradient ($-0.34\pm0.17~{\rm dex}~{\rm deg}^{-1}$) in the center of the galaxy and infer a stream width of $w\sim 0.8^\circ$, roughly 50% larger than the CraII half-light radius. The detection of candidates in the most distant CraII pointings from its center implies that the tidal tails extend beyond our footprint. We compare the CraII stream to $N$-body models with "cored" and "cuspy" dark matter halo progenitors, determining that CraII's density profile is still ambiguous and warrants further modeling.

Historically, the study of the central regions of Type 1 AGN has been limited by the combination of the host galaxy spectrum with strong emission from the accretion disk and NLR/BLR, which prevented us from accurately probing the galactic and AGN properties in the central regions. Integral field spectroscopy allows us to correct for this effect and study both the unobscured cores of AGN host galaxies as well as the uncontaminated spectra of their central engines with unprecedented precision. Using MUSE WFM observations, in this work, we present a combined method for modelling and subtracting QSO light in type-1 AGN alongside results for one such source, PKS 2153-68 (z=0.028), both jetted and gamma-ray emitting. After separating the host galaxy and AGN spectra, we discuss the discovery of an unresolved and yet-to-be-identified high-velocity ($\sim$25000 km/s) short timescale ($\leq$1 yr) variable emission, unlike anything observed in other variable AGN.

When low-mass ($\lesssim 2$ $M_{\odot}$) red giants lose their envelopes to a companion just before the helium flash, the resulting mass transfer can produce binaries hosting hot subdwarfs, horizontal branch stars, and undermassive red clump stars. Recent work predicts a continuum of such products, from fully stripped hot subdwarfs to partially stripped horizontal branch and red clump stars, and suggests that young, metal-rich RR Lyrae can form when partial stripping leaves a helium-burning star in the instability strip. To enable direct comparison with observations, we model these binaries in a simulated Milky Way-like galaxy with a realistic metallicity-dependent star formation history and 3D dust map, generate epoch astrometry using Gaia's scanning law, and fit it with the cascade of astrometric models applied in Gaia DR3. We compare the simulated population to DR3 observations of hot subdwarfs, RR Lyrae, and red giants with high astrometric mass functions. The model significantly overpredicts the number of hot subdwarfs with astrometric binary solutions, partly because the predicted flux ratios are more unequal than observed. It also predicts $\gtrsim 100$ RR Lyrae with DR3 astrometric orbital solutions, while none are observed. We conclude that RR Lyrae in au-scale binaries may be substantially rarer than predicted. In contrast, the model plausibly explains the population of red clump stars with high astrometric mass functions, which we interpret as potential black hole impostors. We predict that $\sim 10 \times$ more stripped-star binaries will be detectable in DR4, whose sensitivity to longer periods will more strongly test wide-orbit systems.

Waves launched from the neutron star surface or inner magnetosphere propagate through the magnetosphere as small perturbations, but can grow relative to the background magnetic field and steepen into ``monster shocks'' -- ultra-relativistic magnetized shocks which can power high-energy emission from magnetars, neutron star mergers and collapse. They occur in magnetically dominated plasma and are described by relativistic magnetohydrodynamics (MHD). We present global relativistic MHD simulations of monster shocks in unperturbed and perturbed (``wrinkled'') backgrounds with a global dipolar geometry. Our simulations confirm analytical predictions for equatorial shocks and provide new insight into the behavior of oblique shocks off the equator. Simulations where the shock is formed through Alfvén mode to fast mode conversion are also presented, demonstrating the generic nature of the monster shock mechanism. We explore how the presence of additional modes in the magnetosphere modifies the shock behavior. Modes of comparable amplitude can fragment the shock front, substantially reduce the magnetization, produce localized enhancements in the Lorentz factor relative to an unperturbed dipole background, and intermittently generate additional shocks along a line of sight.

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining 2D general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models - which predicted proton-rich, $^{56}$Ni-dominated ejecta - strong magnetic fields eject ${\sim 0.2 M_\odot}$ of neutron-rich material $(\langle Y_e \rangle \sim 0.24)$ on dynamical timescales, before neutrino irradiation can raise the electron fraction, enabling strong $r$-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich $(X_{\rm lan} \approx 6\%)$ and dominated by near-infrared emission. We compute synthetic light curves in the LSST and JWST bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy $r$-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

We investigate the structural, photometric, and dynamical evolution of cosmic voids and their galaxy populations from $z=2.09$ to the present, focusing on void size as a key evolutionary parameter. Using void catalogs from four Millennium Simulation snapshots and SDSS data at $z<0.04$, we perform a unified analysis of void demographics, galaxy properties, and internal kinematics. Our analysis reveals clear evidence that cosmic voids exhibit a significant evolutionary trend of becoming progressively emptier toward low redshift, accompanied by a marked decline in the brightness and clustering of their galaxy populations. The void galaxy luminosity function evolves significantly: $M^{*}$ fades and $\alpha$ flattens with time, with large voids hosting brighter, more rapidly evolving galaxies than small voids. Stacked density profiles exhibit a universal shape when scaled by void radius, deepening and building more pronounced walls toward $z=0$. Galaxy spatial distributions reveal persistent size-dependent segregation, with galaxies in large voids lying farther from the center and more strongly clustered. Dynamical analysis of simulations shows coherent outward flows in all voids, with amplitudes decreasing toward $z=0$, providing a physical basis for observed redshift-space distortions. Comparison with SDSS broadly confirms these evolutionary trends but uncovers a non-zero central galaxy population in observed voids -- absent in $\Lambda$CDM predictions -- that may challenge current galaxy formation models in extreme underdensities. Future comparisons with additional simulations and deeper high-redshift surveys will provide stronger tests of $\Lambda$CDM in the most underdense regions.

Weak-interaction-mediated chiral imbalance generation in idealized massless electrons during core-collapse supernovae was once proposed to be the source of strong magnetic fields found in neutron stars. The effect goes by the name of chiral plasma instability (CPI). However, it was found that a finite electron mass damps out this process, inactivating the instability and preventing magnetic field growth. In this work we show that the instability can survive in the presence of abrupt density and temperature perturbation that drives the system sufficiently far out of weak equilibrium. As an example, we work with such perturbations generated by shockwaves which are common during both core collapse as well as neutron star mergers. We find that the chiral imbalance resulting from shock waves, under the right conditions of density and temperature, can sustain the chiral plasma instability despite the damping from the electron mass. Additionally, in an already magnetized medium, the chiral magnetic effect resulting from shock wave density and temperature perturbation can generate substantial ohmic heating. Our results imply that shockwaves during core-collapse supernovae and merging neutron stars can act as a source of strong heating in a magnetized medium as well as CPI.

The detection of sub-solar mass black holes is a milestone of modern astrophysics as it would open a window either onto new stellar physics or could potentially unveil the nature of Dark Matter as Primordial Black Holes. On November 12, 2025, the LIGO-Virgo-KAGRA (LVK) collaboration reported the compact binary merger candidate S251112cm, a system with no obvious EM counterpart, consistent with binary black hole merger with a chirp mass in the range $0.1-0.87 \, M_\odot$. The probability that at least one component has mass $<$1 $M_{\odot}$ is $>99\%$. Inspired by this trigger, we tested if a population of PBHs formed at Quantum Chromodynamics epoch with a broad mass function could account for a signal of this type. Our results, corresponding to a predicted event rate of $0.8 \,\text{yr}^{-1}$ as seen by LVK O3b, suggest that the observed merger rate of $0.23^{+0.86}_{-0.218}\,\text{yr}^{-1}\;(95\%\;\text{C.L.})$ if the trigger is confirmed as an astrophysical event would be compatible with such a model. Our predicted detection rate is also in agreement with current LVK expectations for stellar-mass binaries, remaining consistent with a scenario in which a non-negligible fraction of the $3-200 \;M_\odot$ mergers observed by LVK originate from Primordial Black Holes. If confirmed, this detection would place a lower limit to the PBH abundance $f_{PBH}>0.04$ for our adopted model.

Star formation in galaxies is governed by internal and environmental processes, yet their relative roles are not well understood. In particular, uncertainties in measurements of active galactic nuclei (AGN) host galaxies, combined with modeling limitations, obfuscate the impact of supermassive black hole feedback across environments and over time. Here we address this with a comprehensive analysis of ~60,000 nearby AGNs (z < 0.15 and new environment and halo-mass measurements for ~500,000 AGN and non-AGN host galaxies. This benchmark enables unified comparisons with three prominent cosmological simulations--SIMBA, TNG, and EAGLE--and reveals major, contrasting shortcomings. Simulations fail to reproduce observed trends linking star formation, quiescence, AGN luminosity, stellar mass, and halo mass. While simulations qualitatively capture that AGNs are more common in low-mass halos than in rich groups or clusters, detailed host demographics diverge strongly from observations. Partial agreement exists in the stellar mass distribution within large-scale structures, yet all simulations overproduce quenched low-mass satellites in massive halos, while misrepresenting quenched fractions of massive central galaxies and those in low-density environments, which are sensitive to feedback implementation. Improved AGN physics and modeling of multi-phase gas cooling and flows are required to capture the observed interplay between black holes, galaxies, and halos.

We present new constraints on dust attenuation and star-formation rates (SFRs) for 77 galaxies at redshifts $z=1.43-2.65$, using Paschen emission line detections from the JWST FRESCO survey and ground-based Balmer line measurements from the MOSDEF survey. Using nebular and continuum emission maps, we find that Paschen emission covers a smaller area than continuum emission observed in the F210M (2.1 $\mu$m; rest-frame optical) and F444W (4.4 $\mu$m; rest-frame near-IR) bands, and is preferentially located toward galaxy outskirts. These results suggest that current star formation is concentrated in regions farther from galaxy centers than older stellar populations traced by the continuum, indicative of inside-out star formation. With a careful accounting of slit-loss corrections for ground-based measurements, we calculate nebular reddening and dust-corrected SFRs using the Balmer decrement (H$\alpha$/H$\beta$) and Paschen-to-Balmer line ratios (Pa$\alpha$/H$\alpha$ and Pa$\beta$/H$\alpha$), assuming the Milky Way extinction curve. On average, Paschen-derived reddening and SFRs agree with Balmer-derived values; however, two galaxies exhibit significantly higher Paschen reddening and four show significantly higher Paschen SFRs. We find that non-unity dust covering fractions bias the Balmer decrement toward less reddened OB associations, while decrements involving the Paschen lines are less affected by this bias. These results highlight the enhanced sensitivity of the Paschen lines to the most heavily obscured OB associations in $z\sim2$ galaxies, particularly in galaxies with patchy dust geometries. Future studies using Paschen lines exclusively to measure nebular reddening will yield more robust constraints on the dustiest star-forming regions.

We present Near-Infrared (NIR) and Optical linear polarimetry towards the partially embedded cluster AFGL 6366S. The polarization ranges from 0.44-10.3 per cent in NIR and 0.16-11.22 per cent in Optical bands. The position angle spans $1^\circ - 179^\circ$ in both NIR and Optical bands. About 22 stars exhibit intrinsic polarization signatures. A polarization hole is evident towards the densest ($\sim 3.4 \times 10^{23} \mathrm{cm}^{-2}$) and warmest ($\sim 28.8 \mathrm{K}$) central cluster region. It is attributable to depolarization induced by Radiative Torque Disruption (RAT-D) of large grains and a modest contribution from magnetic-field tangling. The local magnetic field towards the cluster's central region is significantly misaligned with both the large-scale Galactic field and the long axis of the filament present in the region. The field morphology wraps around two dense molecular clumps of radii 0.34 pc and 0.22 pc and $\mathrm N(\mathrm H_2)$ = $(7.9 \pm 1.1) \times 10^{22}$ cm$^{-2}$ and $(4.3 \pm 0.5) \times 10^{22}$ cm$^{-2}$, respectively. The clumps are embedded in the filamentary structure and represent locally accelerated stages of mass accumulation. Gravitationally driven mass flows, largely perpendicular to the local magnetic field, produce a U-shaped field curvature across the filament axis. The plane-of-sky magnetic field strengths towards the two clumps are $ 447.91 \pm 83.81 \mu\mathrm{G}$ and $396.66 \pm 73.64 \mu\mathrm{G}$. The corresponding mass-to-flux ratios ($\lambda \sim 1.34$ and $0.82$) indicate that one clump is magnetically supercritical and the other is subcritical. The Alfven Mach numbers ($\mathcal{M}_A$) $\sim$ 0.395 and 0.393 indicate that both the clumps are in sub-Alfvénic state.

The spacings of super-Earths in multi-transiting systems exhibit a distribution that is broad and mostly featureless, with the exception of notable excesses of planet pairs situated a few percent wide of first-order mean motion resonances (MMRs). In this work, we extend the so-called "breaking-the-chains" model to account for both of these characteristics. Assuming that super-Earths are settled into stable chains of resonances after disk-driven migration, we show that scattering a planetesimal population that contains only a few percent of a system's mass can reorganize primordial chains in remarkable ways. The planetesimal scattering "rattles" the chains by repelling adjacent planet pairs wide of their initial MMRs. Some chains remain rattled but otherwise intact and make up the observed excesses wide of MMRs. In other systems, however, this initial rattling sows the seeds of later orbital instabilities that break the chains entirely. If individual planetesimals' masses are of order a Pluto mass or so, the onset of these instabilities can occur tens or hundreds of Myr after birth, naturally explaining the apparent disappearance of near-resonant pairs on this timescale. The origin of such Pluto-mass debris is currently unknown.

We present a radio-continuum survey of the LMC using the MeerKAT telescope, describe the full-Stokes products included in the first data release, and highlight some initial results. The observations are centred at 1.3 GHz with a bandwidth of 0.8 GHz. The imaging products comprise six fields of view, each encompassing $\sim$5$^\circ$ $\times$ 5$^\circ$ with the resulting images achieving a resolution of 8". The median broad-band Stokes~I image root-mean-square noise value is $\sim$11 $\mu$Jy beam$^{-1}$. The survey enables a variety of astrophysical studies, which we showcase with the presentation of a few findings. Within the LMC we identify a new supernova remnant candidate; present planetary nebulae and Wolf-Rayet stars without previous radio detections; and show the MeerKAT view of the well-known star-forming region 30 Doradus. We also present some examples of interesting foreground and background sources in the field, including the AB~Dor multiple-star system, a radio ring galaxy, a possible Odd Radio Circle, and a remarkable bent-tail radio galaxy.

The volatile budgets of giant planet satellites are critical to unraveling the origin of their building blocks within the circumplanetary disks that hosted them. The Galilean moons Ganymede and Callisto, as well as the Saturnian moon Titan, are known to be anomalously water rich on the basis of their mean densities and interior models informed by gravity data from Galileo and Cassini, characterized by ice-to-rock ratios around unity. Here, we show that the water-ice sublimation line in a decreting circumplanetary disk lends itself to the formation of a water-rich solid reservoir, serving as a natural site for the birthplace of icy satellites. Fundamentally, this reflects how interior to the ice line, water vapor is advected outward, while beyond it, water ice drifts inward as pebbles. Using a semi-analytic model for dust and vapor evolution, we simulate vapor and ice accumulation at the ice line, showing that solids just beyond it achieve steady-state ice-to-rock ratios a factor of a few higher than elsewhere in the disk. For typical disk parameters, this ice buildup occurs within a timescale of a few thousand years. We propose this as a first-order process that explains, at least to some extent, the compositions of three aforementioned satellites. We explore the impact of uncertain turbulence parameters on our results, namely the turbulent Schmidt number and Shakura-Sunyaev alpha, before discussing them in the context of icy satellite D/H ratios. We conclude by evaluating alternative scenarios for explaining water-rich satellites, based on the conversion of CO to CH4, with water as a by-product.

We present spatially resolved measurements of turbulence driving modes across entire extragalactic discs of NGC7793 and NGC1313, using Atacama Large Millimetre/submillimetre Array (ALMA) CO(J=2-1) observations at 13pc resolution. By applying a kernel-based analysis of density and velocity fluctuations, we map the turbulence driving parameter, b, which characterises the balance between solenoidal ($b\sim0.3$) and compressive ($b\sim1$) turbulent driving regimes. b is quantified as the ratio of the turbulent density fluctuations relative to the turbulent sonic Mach number, M. Both galaxies show predominantly solenoidal driving on average for the regions where we find valid results ($b\geq 0.33(\pm 0.05)^{+0.14}_{-0.10}$ in NGC7793; $b\geq 0.24(\pm 0.03)^{+0.10}_{-0.07}$ in NGC1313), noting that this is without including the influences of magnetic fields, making these measurements lower limits. We find substantial spatial variation of b, including localised regions of strongly compressive driving. NGC1313 exhibits higher turbulent Mach numbers and density dispersions than NGC7793, consistent with the disturbed morphology and recent satellite interaction in NGC1313. The turbulence in both NGC7793 and NGC1313 is supersonic ($3\lesssim M\lesssim 20$), and NGC1313 shows a radially decreasing trend of M with galactocentric radius. Radial trends indicate more solenoidal driving in the galaxy centres, potentially reflecting enhanced shear, and increasingly compressive modes in the outskirts. These results demonstrate that turbulence driving varies systematically with galactic environment and cannot be assumed uniform across discs. Our study applies a previously established method to larger scales and new data, linking local turbulence physics to global star formation regulation in galaxies, providing a new avenue for testing theoretical models with future integral field units and ALMA surveys.

The detection of ultra-high energy neutrinos (UHE-$\nu$), with enegies above 10 PeV, has been a long-time goal in astroparticle physics. Autonomous, radio-frequency (RF) UHE-$\nu$ detetectors have been deployed in polar regions that rely on the Askaryan effect in ice for the neutrino signal. The Askaryan effect occurs when the excess negative charge within a UHE-$\nu$ cascade radiates in a dense medium. UHE-$\nu$ can induce cascades that radiate in the RF bandwidth above thermal backgrounds. To identify UHE-$\nu$ signals in data from Askaryan-class detectors, analytic models of the Askaryan electromagnetic field have been created and matched to simulations and laboratory measurements. These models describe the Askaryan electromagnetic field, but leave the effects of signal propagation through polar ice and RF channel response to simulations. In this work, a fully analytic Askaryan model that accounts for these effects is presented. First, formulas for the observed voltage trace and its Hilbert envelope are calculated. Second, the analytic model is compared to UHE-$\nu$ signals at 100 PeV from NuRadioMC, a key Monte Carlo toolset in the field. Correlation coefficients between the analytic signal envelope and MC data in excess of $0.94$ are found, and 99.99% of UHE-$\nu$ signals pass a correlation threshold of $\rho\geq 0.4$. Analysis of RF thermal noise reveals that just 0.2 background events have $\rho\geq 0.4$ in 5 years at a 1 Hz thermal trigger rate. Finally, we describe future work related to the measurement of the logarithm of the UHE-$\nu$ cascade energy.

Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.

Galaxy clusters form through hierarchical assembly, where smaller substructures merge to build the largest gravitationally bound objects in the universe. These mergers, combined with feedback from AGN, filamentary accretion, and other energy injection processes, generate turbulence and perturbations within the intra-cluster medium (ICM). X-ray and Sunyaev-Zel'dovich (SZ) observations can be utilized to measure these ICM density and pressure inhomogeneities, in turn providing constraints on the effective Equation of State (EOS) of the perturbations and ICM velocities. In this work, we analyze deep SRG-eROSITA and Planck/SPT observations of Abell 3266 (A3266), a dynamically complex merging cluster with elongated morphology and significant substructure. We measure pressure and density fluctuations, and compute the power spectra and deprojected 3D amplitudes of these perturbations. We estimate the ratio of pressure-to-density fluctuation amplitudes as $1.00 \pm 0.55$ and non-thermal pressure support $0.068 \pm 0.050$. Density fluctuations are found to be stronger in the northern sector of the cluster compared to the south, consistent with ongoing accretion along a filamentary structure revealed by eROSITA. Further, we find the amplitude of density fluctuations increases with radius, qualitatively consistent with the trend found in cosmological simulations. Uncertainties in our results are dominated by the relatively low sensitivity of current Planck/SPT data, suggesting that improvements in SZ data quality could substantially improve our understanding of ICM energy injection, transport, and dissipation from this technique.

With the discovery of gravitational waves and fast radio bursts, wave optics has become increasingly relevant in astrophysics. This paper studies the behaviour of random gravitational and plasma lenses, presenting the refractive and diffractive expansions, with higher-order terms that allow error estimates and embody the counterintuitive resurgence phenomenon. Specifically, we show that the diffractive expansion converges for a broad class of bounded lens models and provides an efficient description of interference patterns across frequency regimes. Next, building on Picard-Lefschetz techniques, we derive the full refractive expansion to arbitrary order, organising it into a transseries. Near caustics, the standard transseries is supplemented with uniform asymptotics. We study this transseries, with both Borel and hyperasymptotic resummation yielding systematic approximations to lensing integrals at all frequencies. Our results give a framework for modelling wave optics lensing near caustics and beyond the geometric optics approximation and thereby illustrate how tools from resurgence and asymptotic analysis can be applied to practical problems in astrophysics. Near caustic singularities, the post-refractive corrections diverge, while the uniform asymptotic expansion becomes accurate. We use the leading uniform approximation to derive the strong wave optics suppression of off-axis caustics, which clarifies their subdominant role.

We revisit the extensively debated star-forming main sequence (SFMS)-a tight correlation between the star formation rate and stellar mass in both kiloparsec-resolved and integrated galaxies. We statistically explore the fundamental drivers of star formation at global scales, using a large volume-limited sample of 24,954 local star-forming galaxies to overcome the limitations of previous works. Based on the mid-infrared 12 micron luminosity, stellar mass, and g-r color, we estimate the molecular gas mass for the considered sample. At galaxy-wide scales, we establish global relations between the surface densities of the star formation rate, stellar mass, and molecular gas mass . These global density relations are connected with and follow similar trends as the resolved SFMS, the Kennicutt-Schmidt (KS) relation, and the molecular gas main sequence (MGMS). Taking advantage of this large catalog, we show that the scatters in the global KS and MGMS relations are smaller than that of the global relation between the star formation rate surface density and stellar mass surface density, and their Pearson correlation coefficients are higher. More importantly, multivariate regression and partial correlation analyses demonstrate that the apparent correlation between the star formation rate surface density and stellar mass surface density is entirely mediated by the molecular gas surface density, with its best-fit parameters directly derivable from those of the KS and MGMS relations. Overall, our findings suggest that the correlation between stellar mass and molecular gas, as well as that between molecular gas and star formation, are more direct and fundamental. The star-forming main sequence thus appears to be a natural by-product of these two tighter relations.

It is widely expected that the obscured accretion stage can be the initial turning-on stage of active galactic nuclei from quiescent galaxies. We present mid-infrared light curves of NGC 6447 in 3.5$\mu$m and 4.6 $\mu$m bands observed by WISE/NEOWISE, which show an almost monotonic increasing trend of 1.2 mag over 14 years. The optical light curve from ASAS-SN during the same period is consistent with a constant showing no variability. The mid-infrared color evolution shows that the galaxy transitioned into an active galactic nucleus (AGN) in 2018. The SPHEREx spectrum reveals an increasing continuum resembling warm to hot dust emission from an AGN. NuSTAR detected an X-ray source with a 2-30 keV luminosity of $8.4\times10^{41}$ ergs/s at the lower boundary of AGN X-ray emission range, and a factor of >7 variability in one year compared to the Swift upper limit. NGC 6447 was classified as a quiescent galaxy in the literature. The multi-wavelength timing and spectral properties of NGC 6447 are consistent with the expected AGN turning on event, where the obscuring material around the AGN central engine is gradually dispersed, revealing the central engine. This example shows that long-term infrared variability can be a powerful tool to find similar sources. Based on the sample selection statistics, we estimate the duration of the episodes of AGN accretion (duty cycle) signified by the turning-on event as $10^4$-$10^6$ yr.

We present independent imaging analyses of Event Horizon Telescope (EHT) observations of the active galactic nuclei in radio galaxy Centaurus~A and quasar 3C~279 using Generative Deep learning Image Reconstruction with Closure Terms (GenDIReCT), a recently developed machine-learning framework built on conditional diffusion models that uses interferometric closure invariants as primary observables. For Centaurus~A, our reconstruction reveals two prominent emission ridges ($\simeq 80\,\mu$as each) along the jet sheath with a brightness ratio of $1.4\pm 0.1$ and an opening angle of $12.3\pm 0.3$~deg. For 3C~279, we identify three distinct components in the image, with the southern jet ejecta on sub-parsec scale exhibiting a proper motion of $4.6\pm 1.0\,\mu$as over $\approx 5.39$ days away from the northern components, corresponding to an apparent superluminal velocity of $\simeq 10\pm 2$ times light speed. These measurements are consistent with those reported by the EHT Collaboration. The results are significant because we demonstrate that: (1) imaging from interferometric aperture synthesis data, especially in VLBI and most acutely in extremely sparse arrays like the EHT, remains a severely ill-posed and challenging inverse problem, yet closure invariants preserve robust morphological information that can strongly constrain structural features, and (2) more importantly, closure-invariant imaging largely avoids calibration systematics, thus providing a fundamentally independent view of spatial structure with very high angular resolution. The generative nature of GenDIReCT further allows us to sample and characterise clusters of plausible image solutions for each dataset. As a calibration-independent, generative imaging approach, GenDIReCT offers a robust and truly independent blind-imaging tool for current and future VLBI experiments.

GRB 250221A is a short gamma-ray burst (GRB) at redshift $z=0.768$, with a duration of 1.8 s and no extended emission in either Swift/BAT or Konus-Wind bands. A remarkable re-brightening feature in both optical and X-ray bands was observed at $\sim 0.6$ days after the burst trigger, but no supernova or kilonova signature was detected. The burst properties and empirical correlations or distributions (e.g., duration, spectral hardness, location in the Amati correlation, $\varepsilon-$value, $f_{\rm eff}$ parameter, and physical offset) favor a compact binary merger origin. However, a dense circumburst medium with $n\sim 80\rm~cm^{-3}$, obtained by adopting the energy injection into a jet to interpret the late-time re-brightening is inconsistent with the compact binary merger origin. In this paper, we propose a two-component jet model to explain the multiwavelength afterglow observations of GRB 250221A, in which the relativistic narrow jet ($\rm \theta_{c} \sim 3.8^\circ$) produces the prompt and the early decay afterglow emission, while the mildly relativistic wide jet ($\rm \theta_{w} \sim 4.4^\circ$) dominates at later times, resulting in the observed re-brightening feature. If this is the case, one can obtain a lower medium density with $n\sim 0.72\rm~cm^{-3}$ which is a little bit higher than that of short GRBs in merger environments, but falls into the reasonable and acceptable range. Finally, a possible kilonova emission is also discussed within the scenario of compact star merger origin of GRB 250221A.

IceCube has detected neutrino emission from the Galactic Plane (GP) at a significance of $4.5\sigma$, though its origin remains uncertain. Utilizing ten years of IceCube muon-track data, we investigate potential correlations between the GP neutrinos and $\gamma$-ray sources in the first LHAASO catalog (1LHAASO). To avoid issues caused by spectral extrapolation, this analysis focuses on sources detected by the Water Cherenkov Detector Array (WCDA). We employ an unbinned likelihood analysis to search for neutrino emission and constrain the hadronic $\gamma$-ray component of these sources. Neither single-source searches nor stacking analyses reveal significant neutrino signals. The stacking analysis indicates that the 1LHAASO WCDA population contributes at most $\sim$20\% to the diffuse GP neutrino flux measured by IceCube. The total hadronic contribution to the cumulative $\gamma$-ray emission from all WCDA sources is constrained to be at most $\sim$$60\%$, suggesting a predominantly leptonic origin for the $\gamma$-ray emission from the LHAASO source population. Even accounting for unresolved sources below the detection threshold, we estimate the total neutrino flux from all discrete sources (resolved plus unresolved) reaches at most about 40\% of the observed GP neutrino flux. These results support that the bulk of the GP neutrino emission is mainly from truly diffuse processes, i.e., cosmic-ray interactions with the interstellar medium, rather than from unresolved point sources.

High-resolution tri-axial simulations are indispensable for realistically co-modeling the dynamical signatures and the radiative fingerprints of astrophysical jets, which are becoming increasingly important in modern computational studies of jet physics. However, such simulations impose extreme computational requirements that often exceed the capabilities of conventional CPU-based codes. GPU-accelerated simulations offer a transformative solution to mitigate these limitations. In this work, we present a detailed performance benchmarking of the recently developed GPU-enabled PLUTO code (gPLUTO), demonstrating runtime speed-ups ranging from an order of magnitude to (approximately) over 30 relative to CPU-only configurations. A direct comparison between computations of extragalactic jet propagation performed at different grid resolutions confirm the physical fidelity and production readiness of the gPLUTO code, while underscoring the importance of resolving the jet radius adequately to capture the jet dynamics accurately. Leveraging GPU-PLUTO's capabilities, we finally present an application by performing high-resolution simulations of giant radio galaxy jets (GRGs $\gtrsim 1$ Mpc), representing the first such well-resolved 3D study to our knowledge (resolving scales down to 500 pc). These simulations probe a range of environmental effects on GRG jets, clarifying their formation from central galaxies within host cosmic structures, rapid peripheral expansion, and the development of asymmetric cocoon morphologies.

We investigate the gas motions in the core region of the Abell~3395 South subcluster (A3395S) using high-resolution X-ray spectroscopy with XRISM/Resolve. By analyzing the Fe~XXV He$\alpha$ emission line, we directly measure the line-of-sight bulk and turbulent velocities of the intracluster medium. We find that the one-dimensional turbulent velocity is low, at the level of $124\pm21~{\rm km\,s^{-1}}$, while a significant line-of-sight bulk velocity of $263\pm23~{\rm km\,s^{-1}}$ is detected. The coexistence of low turbulence and finite bulk motion suggests that A3395S has not yet reached a dynamically relaxed state. These results are consistent with the non-detection of a radio halo in A3395S, implying that turbulent particle reacceleration is currently inefficient in the cluster core. This study demonstrates that high-resolution X-ray spectroscopy with XRISM provides a powerful means to directly constrain intracluster medium dynamics in merging galaxy clusters, and it provides a reference for future comparative studies of A3395N and A3391 within the same large-scale structure.

We study AT2024kmq and AT2024lhc, two tidal disruption events (TDEs) with blue featureless spectra associated with high-mass black holes ($M_{\rm BH}\sim 10^8\,M_\odot$). Both events show optical precursors consistent with shock dissipation from stream self-intersection. Their X-ray emission is luminous ($L_{\rm X}\sim 10^{44}\,{\rm erg\,s^{-1}}$), highly variable (with minimum observed variability timescales of 1.3\,hr and 4.8\,hr for factor of $\sim3$ flux changes), long-lasting ($>1\,\rm yr$), emerging no later than the optical peak, and well characterized by power-laws with $1.7<\Gamma<3$ (where $f_\nu \propto \nu^{1-\Gamma}$). The X-ray properties and radio non-detections support a compact corona ($\lesssim 10 r_{\rm g}$) producing Comptonized X-ray emission. Using all published featureless TDEs, we find statistically significant bimodality in the distribution of their peak UV/optical blackbody luminosities and radii. We assemble a comparison TDE sample with early-time X-ray observations with eROSITA, in which we find different $M_{\rm BH}$ distributions in TDEs with different X-ray spectral evolution properties: low-mass black holes ($M_{\rm BH} \sim 10^6 M_\odot$) remain soft ($\Gamma>4$) within $t\lesssim 2$\,yr, intermediate masses ($\sim 10^7 M_\odot$) transition from soft to hard at $\sim$1 yr, while high masses ($\sim 10^8 M_\odot$) are hard ($1.5<\Gamma\lesssim 3$) from the outset. We interpret this result as evidence that the soft-to-hard state transition in TDEs occurs at the critical threshold of $\dot{M}_{\rm acc} \sim 0.03 \dot M_{\rm Edd}$ (similar to X-ray binaries), using the fact that the transition timescale predicted by simple disk theory scales with black hole mass as $t_{\rm tr}\propto M_{\rm BH}^{-3/4}$.

The GW190814 event, involving a black hole of mass $22.2$--$24.3 M_{\odot}$ and a compact object of mass $2.50$--$2.67 M_{\odot}$, challenges our understanding of the mass gap between the heaviest neutron stars and the lightest black holes. If the secondary is a neutron star exceeding $2.5 M_{\odot}$, hyperons are likely to appear in its core, softening the equation of state. Rapid rotation can offset some of this softening, enabling higher maximum masses, but it may simultaneously excite the Chandrasekhar--Friedman--Schutz $r$-mode instability. Bulk viscosity arising from nonleptonic weak interactions in hyperonic matter provides an efficient damping mechanism that can stabilize such configurations. In this work, we investigate the combined effects of rotation, thermal evolution, and hyperon-induced bulk viscosity on the stability of massive neutron stars. We demonstrate a direct connection between the suppression of $r$-mode instabilities and the long-term dynamical stability of hyperon-rich stars, offering a plausible interpretation of the GW190814 secondary as a rapidly rotating, hyperon-rich neutron star rather than a low-mass black hole. Our unified framework extends beyond previous studies restricted to static equations of state or extreme viscous damping assumptions, providing new insights into the stability of massive, exotic neutron star configurations.

We present results from coordinated multiwavelength observations of the SS 433, obtained with XRISM, optical telescopes, and near-infrared camera during 2024 April and 2025 March. The XRISM exposures amounted to ~200 ks in 2024 and ~100 ks in 2025. With XRISM/Resolve's high spectral resolution and large effective area, we clearly resolved numerous emission lines even in short time segments, achieving improved accuracy in Doppler-shift measurements relative to earlier observations. The simultaneously obtained X-ray and optical Doppler shifts suggest a possible tendency for the optical emission to lag slightly behind the X-rays. In the Resolve data, the Doppler shifts of the two jet components exhibited apparent asymmetries, with jet speeds fluctuating around ~0.26$\pm$0.01$c$ in 2024 and ~0.30$\pm$0.01$c$ in 2025. The velocity variations indicated modulations on a timescale of ~6.3 d, with a phase offset of about -90$^{\circ}$ relative to the nutation cycle. The observed line widths and flux of the approaching and receding jets appear consistent with the expected geometrical effects, indicating systematically larger line widths in the inner regions of the jets, as proposed by Shidatsu et al. (2025). Optical light curves show flares of ~400 s in 2024 and ~1600 s in 2025, with amplitudes up to ~15% during out-of-eclipse intervals, while the XRISM/Xtend light curves show no significant variability within the overlapping intervals and given the statistical uncertainties. Near-infrared photometry in 2024, obtained during an out-of-eclipse interval at a different epoch from the optical observations, showed no flare-like variability, and the X-ray band also remained constant within uncertainties. These coordinated observations provide a foundation for future XRISM studies aimed at probing the dynamical properties of the relativistic jets in SS 433.

On the periphery of galaxy clusters, moderately high galaxy densities and velocity dispersions favour interactions and mergers that influence galaxy evolution prior to cluster infall. Observational studies of this phase in dwarfs remain rare. We present a high-resolution study of the merger remnant VCC 693 in the outskirts of Virgo cluster, using observations from the Atomic gas in Virgo Interacting Dwarf galaxies (AVID) project. We explore the origin of VCC 693 and the consequences of the merger on its star formation and structure through a joint analysis of VLA and FAST HI emission line observations, together with complementary optical imaging and spectroscopy. We employ hydrodynamical simulations to help interpret the observations. Our analysis favours a near-major merger between two dwarfs with a stellar mass ratio of 3:1-4:1, with one likely gas-poor progenitor (i.e., a damp merger). The optical appearance of VCC 693 is dominated by complex tidal structures throughout the system, whereas the HI gas has settled to a regular rotating disk. Compared with similar-mass dwarfs, the central star formation and gas-phase metallicity are moderately enhanced. The global star formation rate, HI gas content, and HI-to-optical size ratio of VCC 693 are broadly consistent with those of typical dwarfs of similar mass, albeit somewhat lower. Decomposition of the HI rotation curve into baryonic and dark matter indicates a high halo concentration, suggesting post-merger relaxation into a more centrally peaked configuration. Together with two recent studies of AVID post-merger systems, these results support the view that even major dwarf mergers can produce remnants with overall stellar structures indistinguishable from ordinary dwarfs, and that the environmental effects in cluster outskirts can promote damp or mixed mergers, constituting an integral part of galactic pre-processing.

Current studies of stellar dynamos primarily focus on spherical stars, leaving their behavior in distorted stars largely unexplored. We utilize stars of varying distortions to examine the relation between stellar cycle periods ($P_{\rm cyc}$) and rotational periods ($P_{\rm rot}$), which are closely linked to dynamo processes. By analyzing a sample of tidally distorted M dwarfs in cataclysmic variables, we identify an anti-correlation between $P_{\rm cyc}$ and $P_{\rm rot}$, in contrast to the lack of such a relation in single M dwarfs. This means that stars with greater deformation have longer cycle periods, suggesting variations in dynamo behavior under non-spherical geometries. Our numerical simulations further reveal that, the thermal convection weakens in highly distorted stars, and subsequently, the differential rotation is also reduced. These effects may lengthen the conversion timescale between poloidal and toroidal magnetic fields, potentially explaining the newly discovered $P_{\rm cyc}$-$P_{\rm rot}$ relation in cataclysmic variables.

We report photometric observations in Johnson B and V bands of the short term variability (flickering) of Mira (omicron Ceti). The observations were performed during 7 nights in the period August - October 2025, in the course of the last minimum of the Mira pulsations. The observed peak-to-peak amplitude of the flickering is 0.11-0.28 mag in B band. For the flickering source we find luminosity in the range 0.10-0.46 L_sun. Using the amplitude-flux relation, we estimate an average luminosity of the accretion disc L_d = 0.91 +/- 0.28 L_sun. Assuming that the white dwarf accretes material through Wind Roche Lobe Overflow (WRLOF), we find that Mira B is a low mass white dwarf with M_wd = 0.24 +/- 0.04 M_sun accreting at a rate 6.8x10^{-9} M_sun/yr. This value of the mass is in the range of the extremely low mass white dwarfs. The data are available on Zenodo: this http URL

Stellar flares are sudden brightenings caused by magnetic reconnection and are frequently observed on late-type stars. High-cadence photometry of flares provides valuable insights into the mechanisms of these events, yet such observations remain scarce. We seek to explore the sub-second fine structure of stellar flares and assess the information content in high-speed photometry. New 0.3 s-cadence photometry from a six-year-long observing campaign of the active M-dwarf AD Leo is presented. We use time--frequency analysis to detect quasi-periodic pulsations in the decay phase of flares. We explore statistical measures of time series complexity of the detected flares to quantify the information gain achievable with high-cadence photometry. We detect 42 flares in 211 hours of observations. The flare frequency distribution is consistent with the previous literature. We find no quasi-periodic pulsations with periods below a few seconds, and identify two candidate signals with periods around 1 and 3 min. Using different measures of complexity on the binned flare light curves we confirm the advantages of high observing cadence. However, we also find a plateau up to a binning of ~4--5 s for a few complex flares, suggesting that an exposure time of a few seconds is usually enough to retain most of the information carried by a single-filter observation. New photometric observations of AD Leo revealed sub-structures of flare light curves on the timescale of a few seconds, but we found no features on timescales below that.

We present a sound horizon independent measurement of the Hubble parameter using a multiprobe large-scale structure analysis. Removing the dependency on the sound horizon with a rescaling procedure at the matter power spectrum level, we analyse the BOSS full-shape power spectrum and bispectrum (for the first time) using the effective field theory of large-scale structure up to one loop. We combine this analysis with the auto- and cross-angular power spectra from the DESI Legacy Imaging Survey DR9, the $3 \times 2$pt analysis from DES Y3, and the CMB gravitational lensing power spectrum from Planck PR3. Our baseline analysis, that does not rely on supernovae data, yields $h = 0.702^{+0.022}_{-0.024}$, $\Omega_m = 0.310 \pm 0.013$, and $\sigma_8 = 0.799 \pm 0.020$, corresponding to $3-4 \%$ precision measurements. When adding supernovae data from Pantheon+, we obtain a $2.6 \%$ measurement of $h$, with $h = 0.686 \pm 0.018$. We further note that our EFTBOSS analysis indicates a slight deviation of the BAO scale parameter (at $1.8 \sigma$) from its $\Lambda$CDM value, caused by the small scales of the bispectrum. We finally use the sound horizon-free EFTBOSS analysis as a diagnosis for the presence of new physics, finding that our results are consistent with the recent hints of evolving dark energy.

QU-fitting is a powerful tool for interpreting spectro-polarimetric radio continuum observations by linking them to physical models, enabling estimates of the magnetic fields in, for example, the Milky Way, galaxy clusters, and radio jets. We present a comprehensive investigation into the effectiveness and limitations of QU-fitting within the ASKAP POSSUM survey frequency ranges (800-1088 MHz and 1296-1440 MHz) with projections to other spectro-polarimetric radio observations. We simulate different physical polarization sources: Faraday simple, Burn slab, internal turbulence, external turbulence, and two-component models in the POSSUM frequencies, and assess their observational degeneracies and fit accuracies. Our results highlight the model-dependent nature of reliable fitting and identify specific regions of parameter space where model selection, and therefore characterization of the physical medium, becomes ambiguous. For QU-fitting we find the Bayes factor, computed using the marginal likelihood, outperforms more traditionally used goodness-of-fit metrics such as Bayesian Information Criterion (BIC), Akaike Information Criterion (AIC), and chi-squared for model selection. We provide empirical relationships to delineate the boundaries where model distinguishability is impossible. Finally, we evaluate how accurately QU-fitting recovers model parameters and their associated uncertainties, thereby assessing its ability to correctly characterize the Faraday-rotating medium in both point and extended sources in Faraday depth space.

We study the formation of stellar bars using 145 simulations of disc galaxies embedded in live and static dark matter haloes. We use the exponential bar growth timescale, $\tau_{\rm bar}$, to quantify how disc structure and kinematics regulate the onset and rate of secular bar formation. We extend previous work to thicker and more turbulent discs, motivated by those observed at high redshift ($z>1$). By revisiting several commonly used disc stability criteria - the Efstathiou-Lake-Negroponte parameter ($\epsilon_{\rm ELN}$), the Ostriker-Peebles ratio ($t_{\rm OP}$), and the disc stellar mass fraction within 2.2 disc scale radii ($f_{\rm disc}$) - we find that $\tau_{\rm bar}$, when expressed in terms of the disc's orbital period, follows a tight power law with each criteria. In Milky Way-like discs embedded in live haloes, bars form within a Hubble time if $f_{\rm disc} \geq 0.18$, $t_{\rm OP} \geq 0.27$, and $\epsilon_{\rm ELN} \leq 1.44$. We show discs with higher velocity dispersion experience delayed bar growth and introduce an empirical relation that correctly describes the bar formation timescales of all our live halo models. Bars in static haloes grow at roughly half the rate of those in live haloes and require substantially greater disc instability to do so.

We develop a dark-sector selective trace-coupled extension of gravity in which the matter--curvature coupling depends exclusively on the trace of the dark-matter energy--momentum tensor, $T_{\chi}$, defined from a canonical dark-matter field $\chi$. This construction provides a microphysically specified trace sector, removes the usual matter-Lagrangian ambiguity of $f(R,T)$-type models, and preserves minimal coupling of visible matter by design. We derive the full field equations, the exact dark-sector exchange structure, and the linear scalar-perturbation system in gauge-ready form. In the sub-horizon regime, we derive effective modified-gravity functions governing structure growth and light deflection, and show that the model generically produces correlated, scale- and time-dependent departures from General Relativity in growth and lensing observables. Building on this structure, we formulate a perturbation-focused multi-probe framework based on redshift-space distortions, weak lensing, and CMB lensing, explicitly targeting degeneracy breaking beyond background-expansion tests. The analysis establishes the action-level and perturbation-level foundations of the model and provides a conservative, reproducible framework for translated linear-regime constraints within a dark-sector selective modified-gravity setting.

Dust grains and polycyclic aromatic hydrocarbons (PAHs) actively contribute to the thermodynamics, chemistry, and radiative state of the interstellar medium (ISM), yet most ISM models and galaxy simulations either exclude them altogether or adopt simplified treatments. We present CALIMA, a new module for dust and PAH formation and evolution in radiation-hydrodynamics simulations for RAMSES, designed to self-consistently couple dust physics to radiative transfer and non-equilibrium thermochemistry in a multiphase ISM. The model employs a two-size, two-composition dust framework with log-normal grain populations, explicitly evolving stellar dust injection, turbulence-informed gas-phase accretion, shattering, coagulation, thermal and non-thermal sputtering, and shock destruction, while PAHs are separate components with their own evolution. The evolving dust populations and radiation field determine local, wavelength-dependent opacities, photoelectric heating efficiencies, grain-assisted recombination, dust-gas collisional heating/cooling, and H$_2$ formation on both grains and PAHs. Updated treatments of thermal sputtering and collisional cooling that include finite grain sizes and modern ion-solid physics reduce sputtering rates at high temperatures and extend the regime where dust significantly cools hot gas. One-zone ISM tests show that dust and PAH evolution modifies classical thermal phase diagrams and C-bearing chemistry, while isolated disc galaxy simulations reveal environment-dependent variations in dust-to-metal ratio, small-to-large grain ratio, PAH fraction, and interstellar radiation field intensity that drive non-trivial structure in infrared emission, UV transparency, and H$_2$ formation. CALIMA provides a physically motivated framework to interpret dust- and PAH-based observables and to assess dust-mediated feedback in galaxy formation across cosmic time.

Magnetized rotating neutron stars, or pulsars, are a possible end product of massive star evolution. Their relativistic wind successively interacts with the supernova ejecta of their defunct progenitor, then with the circumstellar medium of the progenitor, and eventually with the interstellar medium. If a massive star is static with respect to its ambient medium, then its resulting circumstellar medium is elongated along the direction of the local magnetic field, and its supernova remnant transiently appears as a rectangle. The pulsar wind nebula forming in it is, in its turn, elongated, as long as the pulsar axis of rotation matches the direction of the local magnetization. In this work, we explore how the angle between the direction of the local magnetic field of the interstellar medium and the pulsar axis of rotation influences the shaping of its pulsar wind nebula with 3D MHD simulations are carried out with the PLUTO. We use those models to perform radiative transfer calculations to derive non-thermal radio emission maps of the pulsar wind nebulae. When the polar elongation of the pulsar develop, they bend in opposite directions under the effects of the cavity carved by the stellar wind and already filled by supernova ejecta. This induces a complex distribution of magnetized supernova ejecta and pulsar wind, resulting in various observable structures, appearing as rectangles, circles, or irregular oblong shapes, in the radio waveband. The angle between the direction of the pulsar rotation axis and that of the local ambient magnetization is a governing parameter for the shaping and non-thermal radio properties of the pulsar wind nebulae of static massive stars; however, the mixing of material, once the pulsar wind nebula is old (50 to 80 kyr), is not strongly affected by that factor.

All gas-rich galaxies in cluster environments are expected to experience ram-pressure stripping from the intra-cluster medium. However, only a fraction of these develop ongoing star-formation in their stripped tail, becoming the so-called ``jellyfish'' galaxies. In this work we provide observational evidence that magnetic fields can signal differences in the extraplanar star formation and explore what are the physical conditions that lead to the formation of a jellyfish galaxy. We first focus on JO147, a jellyfish galaxy that features weak star formation activity in its tail. Using MeerKAT radio continuum observations, we discover polarized emission only in a small fraction of its tail, with an average fraction of $~10\%$, and a low Mach number $\mathcal{M}=1.3-1.6$, which suggests a possible association between magnetic field draping, shock-compression of the gas, and extraplanar star formation activity. Then, we test this scenario in a sample of 17 jellyfish galaxies from the GASP project. We combine dynamical models for their orbits within the host clusters with realistic cluster temperature profiles to infer their Mach number, and we find a positive correlation between it and the star formation activity in their tail. We conclude that supersonic motion is a necessary condition for triggering star formation in the stripped tails of jellyfish galaxies. Our findings provide empirical evidence that the critical factor preventing the stripped gas evaporation is the shock compression induced by the supersonic motion through the cluster. This process likely enhances the magnetic field surrounding the galaxy and the properties of the stripped material.

A common task in astronomical research is to estimate the physical parameters (temperature, mass, density etc.) of a gas by using observed line emission. This often requires a calculation of how the radiation propagates via emission and absorption (so-called radiative transfer). In radio and infrared astronomy, the Fortran code RADEX (van der Tak et al., 2007) is a popular tool to solve the non-LTE radiative transfer of a uniform medium in a simplified geometry. I present pythonradex, a Python re-implementation of RADEX. Written in Python, it provides an easy and intuitive user interface, improved performance as well as additional functionality not included in RADEX, such as continuum effects and overlapping lines. In addition, pythonradex provides a self-consistent computation of the total flux for all geometries, including spherical geometries.

WZ Sge-type stars are an extreme subclass of dwarf novae characterised by very rare, large-amplitude superoutbursts. Within the disc instability model (DIM), such events are explained as being triggered by enhanced mass transfer from the donor star. We present an analysis of observations of a sample of WZ Sge-type systems in deep quiescence to assess the consistency of DIM predictions with their observed properties. We find that accretion discs in quiescent WZ Sge-type systems have very low mass-accretion rates of a few $\times$$10^{-13}$ M$_\odot$ yr$^{-1}$. The discs are entirely optically thin, and their physical conditions -- such as surface density and effective temperature -- remain well below the DIM thresholds required to trigger an outburst. Observationally, no increase in disc brightness is detected prior to the superoutburst, indicating the absence of a transition to an optically thick state, in contrast to DIM predictions of a gradual disc thickening preceding the instability. We therefore find no observational evidence that superoutbursts in WZ Sge-type systems are triggered by enhanced mass transfer from the donor. Furthermore, the inferred mass-transfer rates in these objects ($\dot{M}_{\rm tr}$~5$\times$$10^{-12}$ M$_\odot$ yr$^{-1}$) are at least an order of magnitude lower than commonly assumed. We argue that the widely adopted value of $\dot{M}_{\rm tr}$ for the prototype object WZ Sge is likely overestimated. Finally, we show that in quiescence the accretion disc radius in all systems is close to the tidal truncation radius and exceeds the 3:1 resonance radius, confirming earlier results and calling into question the standard interpretation of superhump formation.

We forecast cosmological parameter constraints for a cosmic shear analysis of the Rubin Observatory Legacy Survey of Space and Time (LSST), defining an analysis framework that can accurately recover the $\Lambda$CDM model in the presence of astrophysical and data-related systematics. When accounting for our present uncertainty on the suppression of the non-linear matter power spectrum through baryon feedback, we find that the error on the composite parameter $S_8=\sigma_8\sqrt{\Omega_{\rm m}/0.3}$ almost doubles compared to an LSST analysis which neglects this astrophysical phenomenon. After the first year of observations, LSST will extend beyond the magnitude limit of existing representative spectroscopic calibration samples, requiring photometric redshifts to be calibrated using an alternative strategy. Adopting literature measurements of the reduced redshift calibration precision found from galaxy cross-correlation techniques, combined with current levels of baryon feedback uncertainty, we forecast final year LSST cosmic shear constraints that barely improve upon the first year analysis. This forecast therefore serves as encouragement to the community to develop methodology and observations to constrain models of baryon feedback and enhance photometric redshift calibration at depths where spectroscopy is unrepresentative. With tight priors on both these systematic terms, we forecast that LSST cosmic shear can deliver constraints on $S_8$ that are more than five times as constraining as existing cosmic shear surveys.

It has been 10 years since the initial discovery of Ultra-Diffuse Galaxies (UDGs) in the Coma cluster and the revelation that large, low surface brightness galaxies may constitute a greater fraction of galaxies than first thought. This left an open question: Are UDGs something special, or just an extension of the previously known dwarf galaxy population? Seeking to answer this question, in the decade following, dedicated simulations have studied and proposed a myriad of formation pathways to create UDGs. Observations have then pushed the limits of world-class observatories to perform detailed studies of these galaxies in large numbers across the full range of environments in the local Universe. These observations stress test simulations and challenge previous galaxy formation wisdom, with UDGs posing many open puzzles beyond just their unknown formation mechanism. To provide a few pertinent examples: there is observational evidence that not all UDGs follow the standard stellar mass -- halo mass relationship; there is evidence for UDGs with extraordinarily high levels of alpha enhancement; and there is evidence that some UDGs are much more globular cluster rich than other dwarfs of similar stellar mass. In this Dawes review, we undertake the task of summarising the decade of science since the discovery of UDGs. We focus on the quiescent population of UDGs and review their general properties, their proposed formation scenarios, their internal properties and their globular cluster systems. We also provide a brief conjecture on some future directions for the next decade of UDG research.

A significant fraction of stars in both the Galactic field and stellar clusters are members of binary systems. Understanding their properties is therefore essential for a comprehensive view of stellar structure, evolution, and cluster dynamics. Despite extensive studies of cluster binaries, key issues remain unresolved, particularly for photometric binaries among low-mass stars. While the binary fraction in the field strongly depends on stellar mass, cluster studies have generally suggested an approximately constant fraction over the limited mass ranges explored. In addition, the mass function (MF) of very low-mass stars is still poorly constrained in clusters older than a few hundred Myr. We use deep Hubble Space Telescope imaging of the intermediate-age open cluster NGC 2158 to investigate its binary population and derive the luminosity and MFs down to ~0.14 solar masses, enabling the first detailed analysis of binaries in this cluster. We measure a global binary fraction of 38%, consistent with other open clusters, and find a clear mass dependence: it decreases from ~52% at 1.0 solar masses to ~11% at 0.2 solar masses. This trend mirrors that of Galactic field stars, suggesting similar binary properties. The MF is characterized by three regimes: high-mass stars (alpha= -2.49 +- 0.19), low-mass stars (alpha= -1.11 +- 0.09), and very low-mass stars (alpha= -0.08 +- 0.07). The slope change near 1.0 solar mass agrees with recent surveys, though we find a deficit below ~0.3 solar masses. We also detect a main-sequence discontinuity around ~0.3 solar masses, possibly linked to the 3He-driven instability predicted by stellar models and analogous to the Jao Gap seen in nearby field stars.

Cosmological simulations suggest that various galaxy properties depend on their location within the cosmic web. Yet direct observational evidence of the dependence of star formation activity on distance to filaments remains scarce and is missing at z>1. We investigate how starburst, main-sequence (MS), and quenched galaxies are distributed with respect to cosmic web filaments, and how this distribution evolves with redshift. We first use the SIMBA cosmological simulation to predict the redshift evolution of the mean distance to the closest filament from z=3 to z=0 for different galaxy populations after removing stellar-mass dependencies. We then measure the corresponding signal in the COSMOS field, using COSMOS2020 and COSMOS-Web data, where accurate photometric redshifts enable reconstruction of the projected cosmic web from z=2 to z=0.5, and starbursts are identified through far-infrared spectral energy distribution fitting. In agreement with the results from SIMBA, starburst galaxies are found closer to filaments at z>1 and at larger distances at z<1, MS galaxies occupy intermediate environments with little evolution, and quenched galaxies show progressively shorter distances to filaments toward low redshift, with a crossing between starburst and MS populations around z~1. In COSMOS-Web, the relative evolution in the average distance to filaments between starburst and MS galaxies is detected at a significance level of at least 5{\sigma}. We show that a minimal toy model in which the only environmental ingredient is the sSFR-filament distance modulation measured in simulations is sufficient to reproduce the observed differential evolution of the average filament distance between starburst and MS galaxies. These results show that the imprint of large-scale environmental effects on the star formation activity of galaxies, predicted by simulations, is detectable from z=2 down to z=0.5.

Many small satellites in the Solar System have sub-synchronous orbits, meaning their orbits are decaying due to tidal dissipation. Unless they have substantial material strength, they will eventually tidally disrupt before reaching their planet's surface. We studied the fate of rubble-pile satellites as they migrate inward, with a particular focus on the case of Phobos. We used a combination of analytic estimates and numerical simulations to determine the failure mode and tidal disruption distance of a Phobos-like satellite, as a function of its shape and cohesive strength. Both analytically and numerically, we identify a regime for satellites with low cohesive strengths whereby their surface can be tidally stripped without undergoing internal failure. Our numerical simulations demonstrate that Phobos will be destroyed beyond 2 Mars radii if it has a bulk strength similar to those estimated for small bodies recently visited by spacecraft. Based on our results and some additional arguments, we suggest that previous studies on the fate of Phobos have overestimated its strength, and therefore underestimated its tidal disruption distance. We also speculate that if Phobos undergoes some tidal stripping, its ultimate fate may be determined by runaway collisional erosion rather than a pure tidal disruption. However, the ultimate tidal disruption distance for Phobos will depend on its unknown internal structure and bulk material properties, which will be constrained by JAXA's Martian Moons eXploration (MMX) mission and its IDEFIX rover. These results have implications for theories about the origin and evolution of the Martian moons and for tidal disruptions of other small, irregularly shaped satellites.

The dark photon is a new hypothetical gauge boson arising in extensions of the Standard Model, and constitutes a compelling dark matter candidate. As dark photon dark matter (DPDM), it can interact with electromagnetic fields via kinetic mixing, and the inelastic scattering process $\gamma \gamma' \to e^+ e^-$ becomes kinematically allowed for gamma rays above a characteristic energy threshold. This interaction imprints unique spectral attenuation features at very-high-energies (VHE), offering an observational probe of DPDM models. Using the Cherenkov Telescope Array Observatory (CTAO) Instrument Response Functions (IRFs), we simulate observations of VHE sources and forecast novel sensitivities to the kinetic mixing parameter for the photon-dark photon scattering process. Our study focuses on three key astrophysical targets: the Crab Nebula and the blazars Markarian 421 and Markarian 501. Additionally, we investigate the impact of dark matter spikes around black holes on the upper limits. Our results demonstrate that CTAO can probe the DPDM parameter space down to a mixing parameter of $\varepsilon \sim 10^{-8}$ for masses around $m_{A^{\prime}} \sim 10^{-1}~\textrm{eV}$ through high-energy spectral attenuation, at a $95\%$ confidence level.

We present the largest collisionless $N$-body cosmological simulations in a MOdified Newtonian Dynamics (MOND) cosmology to date. Our 4 simulations cover $\Lambda$CDM as a baseline, a MOND with hot dark matter model known as $\nu$HDM, and 2 unphysical models we call $\Lambda$HDM and $\nu$CDM to test the individual contributions of hot dark matter and MOND gravity, respectively. $\nu$HDM reproduces the CMB power spectrum whilst also theoretically matching cluster dynamics and preserving MOND predictions for galactic rotation curves. We test its viability on cosmological scales using simulations with $256^{3}$ particles in a box of size $800/h$ comoving Mpc. We find generically that the MOND models massively overproduce large-scale structures by $z=0$, with a most massive cluster in $\nu$HDM of $\approx 5 \times 10^{17} M_{\odot}/h$ and typical peculiar velocities of several thousand km/s. We also explore a local void solution to the Hubble tension in these models. Analogues to the observed "Local Hole'' do form in the MOND models, but values for the deceleration parameter $<-1.5$ in these regions prevent a satisfactory resolution to the Hubble tension. Whilst $\Lambda$CDM significantly underpredicts the observed bulk flow in Cosmicflows-4, the high peculiar velocities that arise in the MOND models create the opposite problem, ruling out $\nu$HDM at $>5\sigma$ confidence. Observations clearly require a much milder enhancement to the rate of structure growth in $\Lambda$CDM than is provided by the $\nu$HDM paradigm. Our results also suggest that replacing cold dark matter with hot dark matter is unlikely to provide a viable cosmological model, regardless of the gravity law.

We present scientific and technical justifications of a potential coordinated Euclid/CSST/JPCam/JUST survey of the Euclid Deep Field North (EDF-N), aimed at probing the multi-phase circumgalactic and intergalactic medium (CGM/IGM) at the cosmic noon over ~20 deg$^2$. The survey is structured around three connected goals: (1) improving photometric redshift (photo-z) accuracy through the combination of broad- and narrow-band photometry, enabling reliable identification of large-scale structures; (2) probing extended CGM emission with dedicated narrow-band imaging; and (3) mapping foreground IGM via absorption-line spectroscopy of background galaxies. Together, these components establish an integrated observational framework to investigate galactic ecosystems -- linking galaxies to their circumgalactic and intergalactic environments -- at cosmic noon. We show that the J-PAS-like narrow-band system used in JPCam substantially improves photo-z accuracies from only the Euclid/CSST broad-band data, especially for star-forming galaxies at z~1.0-1.4. This enables the identification of galaxy groups and (proto-)clusters directly from photo-z measurements. Stacked JPCam narrow-band imaging should also detect extended [O II]-emitting CGM halos. We then construct mock 3D gas distribution model and realistic galaxy catalog, and further construct mock CSST and JUST background galaxy spectra adding Lyalpha and Mg II absorptions. The reconstructed 3D H I field from CSST Lyalpha forest reliably recovers large-scale structures; however, our simulations indicate that detecting diffuse IGM Mg II absorption with JUST is infeasible, either through spectral stacking or via the two-point correlation function method. We conclude that constraining the metallicity of the diffuse IGM will require significantly deeper and higher-resolution spectroscopy expected from future facilities such as the 39 m E-ELT.

Mira stars are popular with amateur astronomers because they are generally bright and their pulsations produce large changes in brightness. This also means they are eminently suitable to study using spectroscopy. We report here on spectroscopic observations of five Mira stars: RY Cep, SU Cam, T Cep, V667 Cas and OR Cep. We describe the lifecycle of a Mira star and explain how the carbon-oxygen balance in its atmosphere determines its spectral classification. Atmospheric pulsations in these stars produce shocks which generate Balmer emission. We present graphically their different spectral behaviours and development of emission lines.

We present a spectroscopic study of low-mass galaxies (LMGs;$10^8\leq\rm M_*/M_\odot\leq10^9$) at $z\sim0.15$ in COSMOS field, and compare it to a control sample of intermediate-mass galaxies (IMGs;$10^9\leq\rm M_*/M_\odot\leq10^{10}$) at $z\sim0.35$. We examine their star formation rates (SFRs), dust attenuation properties, and the relationship between nebular and stellar reddening. For both samples, SFRs derived from H$\alpha$ are strongly correlated with SFRs from fitting simple star formation histories (SFHs) to the galaxies' spectral energy distributions. In fitting a joint SFR-$\rm M_*$ relation, we obtain a slope of $\rm {\Delta log(SFR_{H\alpha})}/{\Delta log(M_*/M_\odot)}=1.01\pm0.03$, indicating that fair ensembles of SFHs for galaxies at these stellar masses are well-described by scale-free, self-similar forms. We also examine their dust attenuation properties and the relationship between nebular and stellar reddening, exploring how these quantities vary with stellar mass and specific SFR (sSFR). Nebular attenuation increases with stellar mass for IMGs but is lower and less mass-dependent in LMGs, consistent with their reduced dust content. In all cases, stellar continuum attenuation is lower than nebular attenuation, as expected from the two-component dust model. The nebular-to-stellar color excess ratio in both samples is consistent with the canonical factor of 2.27. The ratio is mass-independent, but rises with sSFR in IMGs and remains constant in LMGs. These results suggest that in LMGs, efficient dispersal of birth clouds keeps the differential attenuation approximately constant across sSFR. Thus, although LMGs follow the same global SFR-$\rm M_*$ scaling as massive galaxies, their lower dust content and feedback-maintained ISM produce distinct attenuation behavior relative to IMGs.

Repeating partial tidal disruption events (pTDEs) provide a direct probe of stellar orbits and episodic mass loss around supermassive black holes, but robust identification requires multi-band and multi-epoch evidence. We investigate whether the late-time radio rebrightening of the nuclear transient IGR J12580+0134 in NGC 4845 can be explained as a repeating pTDE, using multi-epoch Karl G. Jansky VLA observations together with X-ray constraints from Swift/XRT and NICER. The radio light curves show two distinct episodes, with the L-band peaks separated by $\approx1513$ days. Modeling the second episode with a synchrotron afterglow framework using Markov Chain Monte Carlo fitting favors a non-relativistic outflow $v\simeq0.03c$, with an isotropic-equivalent kinetic energy of order $10^{50}$ erg propagating in an approximately constant-density circumnuclear medium. No significant contemporaneous brightening is detected by Swift/XRT during the 2016 radio flare, while faint NICER flares in 2023 suggest intermittent low-level accretion. The recurrence timescale and radio energetics therefore make IGR J12580+0134 a possible candidate for a repeating pTDE system, motivating continued sensitive radio and X-ray monitoring to test future reactivations.

We present results from a multi-epoch Very Long Baseline Array (VLBA) survey of compact radio sources in the Orion complex, conducted within both the DYNAMO-VLBA and the GOBELINS projects. Our observations detected 216 compact radio sources, of which 58 yielded reliable multi-epoch astrometric solutions. For these sources, we derived trigonometric parallaxes and proper motions with typical precisions of about 0.05 mas and 0.10 mas yr$^{-1}$, respectively. The measured parallaxes range between 2.26 and 2.65 mas, corresponding to distances of 380 - 440 pc, and delineate the depth of the Orion star-forming complex. We determine mean distances of $405\pm16$ pc for NGC 2068, $403\pm5$ pc for NGC 2024, $407\pm12$ pc for the $\sigma$ Orionis region, $388.5\pm1.7$ pc for the Orion Nebula Cluster (ONC), and $438\pm12$ pc for L1641. A comparison with Gaia DR3 astrometry for 28 common sources reveals negligible mean parallax offsets ($\Delta\varpi=-0.02\pm0.01$ mas) and small systematic differences in proper motions ($\sim$0.07 mas yr$^{-1}$), likely due to residual rotation of the Gaia reference frame. Our results demonstrate the capability of high-precision radio astrometry to map embedded stellar populations and to provide an independent calibration of the Gaia reference system in obscured regions.

We present population synthesis models of viscous protoplanetary disks subject to mild external far-ultraviolet (FUV) radiation fields ($F_{\rm UV}=1\text{-}100\,$G$_0$). Our simulations focus on gas disk evolution, exploring stellar masses drawn from an Initial Mass Function and a range of initial disk conditions. We quantify the fraction of surviving disks across $10\,\mathrm{Myr}$ of evolution, track the evolution of gas disk mass and size, and compare our results with observations of protoplanetary disks in the Upper Scorpius region, including the ten targets studied by the AGE-PRO ALMA Large Program. We find that models combining viscous evolution with external photoevaporation yield disk lifetimes of $3\text{-}7\,\mathrm{Myr}$, consistent with observed dispersal timescales, particularly for $10^{-4} \leq \alpha \leq 10^{-2}$. Low-mass stars ($0.1\,$M$_\odot$) are more susceptible to disk dispersal due to their weaker gravitational binding, with their fraction among all surviving disks dropping from $76\%$ at birth to $51\%$ by $10\,\mathrm{Myr}$. The majority of the long-lived disks are those with low viscosity $\alpha < 10^{-3.5}$ and initial characteristic radius $R_c < 125\,\mathrm{AU}$, while the initial disk-to-star mass ratio does not play an important role. The median gas disk mass and radius of the surviving disks exhibit a sharp decline in the first $0.2\,\mathrm{Myr}$ of evolution, followed by a slight increase that reflects survivorship bias. We also explore correlations between gas disk mass and size vs. stellar mass and FUV strength. Our findings highlight the critical role of external photoevaporation in shaping disk populations even at moderate levels of FUV radiation fields.

Overlapping galaxies, in which a foreground galaxy partially overlaps a background galaxy, offer a unique opportunity to measure dust attenuation, a key nuisance parameter in galaxy studies, empirically and in great detail by modelling the light of both the foreground and background galaxy and inferring the missing light in the overlapping region. However, the current catalogue of overlapping pairs is relatively limited in number compared to catalogues dedicated to individual galaxies. Expanding this catalogue is not only a necessity to facilitate further detailed dust studies beyond the few limited studies conducted thus far, but also to improve pair-to-pair variance and support automated identification through machine learning techniques. To achieve this, we utilise galaxies classified as "overlapping" from Galaxy Zoo DECaLS (GZD-1, -2, and -5), along with images from Data Release 10 (DR10) of the DESI Legacy Survey, in our individual citizen science project to classify these pairs directly using volunteers. This new catalogue will not only provide a wealth of targets for future dust studies but will also contribute to a deeper understanding of these pairs and dust as a whole.

Primary proton-air interactions at ultra-high energies leave a physically interpretable imprint on the correlated fluctuations of the depth of shower maximum and the muon content in extensive air showers. This imprint reflects the stochasticity in the partition of the primary energy among secondary particles in the first interaction. We show that these fluctuations can be accessed through a probabilistic description that isolates sensitivity to hadronic physics in the initial collision, while treating the subsequent shower development as effectively universal. The uncertainties resulting from this universality are smaller than the spread among current hadronic interaction models and comparable to current experimental uncertainties. Consequently, the joint observable space defined by these two quantities provides a new probe of hadron production in kinematic regimes far beyond the reach of human-made accelerators.

The current Fermi-LAT source catalog (4FGL-DR4: 7194 sources over 14 years) was built incrementally from the 8-year catalog. In a survey mission like Fermi, data accumulate on each source over time, so after 16 years (reached in August 2024) and twice the data for the original 4FGL sources we have more precise localization (by 24% on average). It is thus time to generate a new original catalog, which implies, beyond adding the sources newly detectable after two more years, changing the existing source names (derived from their coordinates) and reviewing the associations. We present an early 16-year list (FL16Y) of 7220 sources, which relocalizes all sources and improves a few aspects of the catalog analysis, but still uses the same model of interstellar diffuse emission as 4FGL-DR4.

We report the discovery of strongly absorbed [OI]63um in a sample of 12 DSFGs at 4.2<z<5.8 selected from the SPT survey. This is the first systematic survey of the [OI]63um fine-structure line at z>4. Using ALMA Bands 9 and 10, we obtain spatially and spectrally resolved observations that probe the interstellar medium on sub-kpc scales. Despite reaching sensitivities 10-100x deeper than most previous studies, we detect [OI]63um in emission in only 2 sources at low significance, with the remaining galaxies yielding stringent non-detections over the full velocity range covered by robust detections of other far-infrared lines, including [CII] and [NII]205um. We identify several compact (0.05-0.2") regions having [OI]63um absorption against the far-infrared dust continuum, some of which are possibly reaching below rest-frame CMB radiation level. We also detect narrow, spatially localised [OI]63um emission "escape channels" preferentially detected in regions with weak or absent dust continuum emission. We predict that similar absorption effects may appear in the [CII] line, particularly when concentrating on the regions with the densest foreground material along the line of sight. The [OI]63um line appears to be originate from a mix of compact, high optical depth [OI]63um emitting regions and sub-thermally excited, oxygen-rich molecular clouds dispersed throughout high-redshift starbursts that are capable of absorbing the ground-state line emission. Combined with a comparison to cosmological radiation hydrodynamical simulations, this supports the interpretation that regions with higher gas and dust column densities may lead to weakening an intrinsically strong [OI]63um line emission. We argue that the high [OI]63um optical depth is the dominant effect causing the strong absorption, limiting the diagnostic power of this line to trace regions of massive star formation in high-redshift DSFGs.

The large-scale dynamics of most conventional space and astrophysical plasmas are predominantly governed by Alfvén modes, which are low-frequency magnetohydrodynamic modes existing in magnetized media. At scales smaller than the ion gyroscale or frequencies exceeding the ion cyclotron frequency, the Alfvén modes transform into kinetic-Alfvén or whistler modes that significantly contribute to plasma dynamics. However, this scenario reverses in non-neutral pair plasmas, such as those found in the magnetospheres of pulsars and magnetars, around rotating black holes and in their relativistic jets, as well as in certain laboratory plasmas. In these systems, the large-scale dynamics is governed by hybrid whistler-Alfvén modes, which transform into pure Alfvén modes at smaller scales. We derive the nonlinear equations that describe the dynamics of whistler-Alfvén modes in ultrarelativistic non-neutral magnetically dominated pair plasma and discuss the spectrum of turbulence governed by these equations.

In low-luminosity active galactic nuclei like M87* and Sgr A*, the accretion disk around the central supermassive black hole is tenuous and collisionless. As a result, the usual ideal magnetohydrodynamics (MHD) approximation may not be applicable. In this Letter, we report on the first fully kinetic simulations of the accretion process where the plasma initially has finite angular momentum. The simulated accretion flow behaves remarkably similarly to the magnetically arrested disk (MAD) regime of ideal MHD, reproducing episodes of magnetic flux saturation and eruption typical of MADs. The resemblance to fluid models owes largely to kinetic instabilities, which regulate pressure anisotropy in the disk, allowing fluid terms to dominate the angular momentum transfer. In addition, by handling vacuum regions effectively, our kinetic approach probes the matter supply to the jet funnel. We observe no efficient penetration of the accreting material into this region, which suggests that a pair discharge may be required to sustain the Blandford-Znajek process.

Despite decades of observations, the physical processes governing mass loss from small bodies beyond our Solar System remain poorly constrained. These exocomets are often treated as analogs of Solar System comet, yet the stellar environments they inhabit spans a wide range in terms of luminosity, stellar winds, and evolutionary stage, leading to potentially very diverse physical behaviors. Within our Solar System, small bodies lose material through a range of mechanisms, including sublimation, desorption, impacts, and/or sputtering. Once released, the composition and dynamics of the ejecta are then altered by additional processes, such as dust sublimation, ionization, and radiation pressure. In extrasolar systems, these mechanisms unfold under vastly different radiative and plasma conditions, leading to a rich diversity of mass-loss pathways and observable signatures. This work reviews our understanding of the mechanisms driving mass loss from small bodies and the subsequent evolution of ejecta in diverse stellar environments. We compare the physical and chemical mechanisms that drive gas and dust production, and investigate how they scale with stellar luminosity, temperature, and activity. We then examine the processes that modify the composition of the ejecta (e.g., dust sublimation, dissociation, or ionisation) and its dynamics (e.g., radiation pressure or stellar winds). To illustrate how these processes vary across different stellar environments, we use four well-studied planetary systems as case studies: the Sun, $\beta$ Pictoris, AU Microscopii, and WD 1145+017. By exploring how cometary tails behave under such diverse conditions, this work provides a physical framework for interpreting exocometary activity and sheds light on why A-type stars, such as the famous $\beta$ Pictoris, are over-represented in the population of exocomet-hosting stars.

The far side of the Moon, shielded from terrestrial radio frequency interference and beyond the influence of Earth's ionosphere, should offer a uniquely quiet environment for radio astronomy and cosmological experiments. The radio sky below 30 MHz is largely unexplored and is thought to contain spectral signatures of new physics in the early, high-redshift Universe. Achieving precision measurements in this band requires accurate understanding of antenna performance and systematics. For upcoming lunar surface radio astronomy missions, this modeling will be challenging because antennas will deploy at heights that are only a small fraction of a wavelength above the lunar regolith, where strong coupling between the antenna and the surface can significantly alter impedance, radiation patterns, and efficiency. The challenge is compounded by the layered dielectric structure of the regolith and the tendency for permittivity to increase with depth, both of which are difficult to represent faithfully in numerical simulations. In this work, we review theoretical predictions for the behavior of a simple horizontal dipole above a dielectric half-space, representing the lunar regolith, and compare them with simulation results obtained using the Ansys HFSS integral equation solver. We quantify how the antenna impedance and beam pattern couple to the sky for a representative lunar surface radio astronomy experiment. The results show that surface induced effects decrease rapidly, even for modest increases in antenna height above the regolith. Conversely, a dipole antenna placed on or very near the lunar surface will exhibit complex spectral response that renders systematics control very difficult without detailed information on regolith properties.

Statistics of polarized synchrotron radiation carry information about the properties of the underlying turbulence. Different statistical measures constructed from observables probe turbulence properties in different ways. We consider a setup in which synchrotron radiation is emitted in a distant volume and then passes through a turbulent screen that induces Faraday rotation. Using both MHD simulations and synthetic turbulence spectra, we explore the spectrum of observed polarization directions measured at a single frequency as a diagnostic for recovering the statistics of turbulence in both the emitting region and the Faraday-rotation screen. We compare these results with our analytical expectations. We also compare the spectrum of polarization direction (SPD) with the wavelength-derivative diagnostic introduced and analytically explored by Lazarian \& Pogosyan. We demonstrate that the SPD exhibits greater sensitivity to turbulence in the Faraday screen. We provide an observer-friendly criterion to determine whether the SPD samples turbulence in the synchrotron-emitting region or in the Faraday screen. These results open a practical pathway for extracting turbulence statistics from existing and forthcoming single-band radio polarimetry.

Low-mass dwarf galaxies ($M_{\rm vir} \lesssim 10^9\rm\ M_\odot$) are fundamental cosmological building blocks, yet the physical processes driving their structural diversity remain poorly understood. Recent numerical simulations have suggested a diversity in the stellar-to-halo mass ratio in this halo mass range, but either the number of samples obtained from the same simulation setup or the numerical resolution was limited. We performed high-resolution cosmological zoom-in simulations for eight galaxies with a dark matter halo mass of $\sim 10^9\rm\ M_{\odot}$ up to $t=1.2$ Gyr at which most gas in the galaxies has been expelled. Our samples have a scatter of an order of magnitude in the halo mass at the reionization epoch. The stellar-to-halo mass ratio expected at $z=0$ scatters nearly two orders of magnitude with $5\times10^{-5}$ to $2\times10^{-3}$. We also observed variation in the compactness of their stellar distributions. Some of our simulated galaxies exhibit a stellar half-mass radius of $\sim30$ pc, which is as small as that of ultra-compact dwarfs. The formation condition for such a compact stellar distribution is understood as an analog of the condition for the formation of dense, massive star clusters. We found that when the central gas surface density exceeds a critical threshold ($\Sigma_{\rm gas} \gtrsim 30\rm\ M_\odot \rm\ {pc}^{-2}$), the star formation becomes highly efficient and results in dense stellar systems. These results suggest that UCDs can form in situ even in isolated dark matter halos.

We unify one of the most widely studied frameworks to explain the hierarchical structure of the flavor sector in the Standard Model, the Froggatt-Nielsen mechanism, with cosmic inflation. We propose that the complex scalar field, the so-called flavon, which breaks the Froggatt-Nielsen $U(1)$ symmetry and generates the Yukawa couplings of the Standard Model, to also drive inflation, which we dub as Inflavon. After inflation ends, the decay of the inflavon reheats the Universe, establishing a novel link between early Universe cosmology and flavor physics. As concrete examples, we present realizations where the inflavon potential is described by an $\alpha$-attractor potential. We then compute the resulting CMB observables, specifically the spectral index ($n_s$), the tensor-to-scalar ratio ($r$), and the amplitude of scalar perturbations ($A_s$) as functions of the underlying Froggatt-Nielsen model parameters. We identify the parameter space in Froggatt-Nielsen models involving the scale of Flavor symmetry breaking $\Lambda_{\rm FN}$ and FN charges which is ruled out by Planck and ACT data, as well as the region that could be probed by next-generation CMB experiments like CMB-S4, SO and LiteBIRD. We also discuss inflavon as dark matter and its isocurvature constraints.

We investigate the extent to which pseudo-Newtonian gravitational potentials can reproduce classic tests of general relativity without resorting to full general relativistic formalisms. This is useful for the researchers seeking intuitive insight into relativistic gravity. Focusing on the perihelion precession of Mercury, gravitational redshift, and gravitational light bending, we derive analytical expressions for orbital precession and demonstrate that, with suitable physically acceptable parameters, pseudo-Newtonian approaches can accurately reproduce the observed perihelion advance and gravitational redshift. However, we confirm that no single potential consistently captures all relativistic effects. In particular, while certain parameters yield agreement with general relativity for planetary motion and redshift, they fail to reproduce gravitational lensing over a broad range of impact parameters. Our results highlight both the usefulness and limitations of pseudo-Newtonian methods in modeling gravitational phenomena. Although the pseudo-Newtonian approach cannot serve as universal substitutes for general relativity, especially in strong-field regimes, it provides valuable semi-analytical insight and pedagogical simplicity. Our results indicate the usefulness of the pseudo-Newtonian approach to uncover more complicated phenomena involved with strong field gravity in possibly modification to general relativity.

We construct a pseudo-Newtonian potential (PNP) corresponding to a rotating black hole solution in a modified gravity (MGR) framework using a metric-based prescription. The motivation is to enable realistic accretion disk studies in MGR, where full relativistic MHD simulations remain computationally prohibitive. Effective potentials and the underlying Newtonian-like forces are derived for both massive and massless particles in the equatorial plane, relevant for disk dynamics. The reliability of the PNP is tested by comparing key orbital properties -- marginally stable, marginally bound, photon orbits and energies at marginally stable orbit radii -- with exact MGR predictions. The PNP reproduces the marginally stable and photon orbits exactly, while marginally bound orbits and specific energies deviate by less than about 7-10%. The influence of the MGR parameter on particle dynamics and effective potentials is analyzed, revealing non-trivial departures from simple Newtonian intuition. The study demonstrates that the proposed PNP accurately captures essential spatial properties of MGR spacetime and provides an efficient, physically consistent tool for investigating accretion phenomena and strong-gravity astrophysics beyond general relativity.

Black holes are some of the most interesting objects in the universe. While they first arise in the complicated behavior of general relativity, the physical laws ruling their behavior are surprisingly simple. For example, one of the core facts about black holes is that their area never decreases, much alike the entropy in thermodynamics. In this note directed at introductory physics students and their instructors, we use this similarity to understand properties of black hole physics using standard techniques from an undergraduate course in thermal physics. We explore the never-decreasing nature of black hole area to obtain bounds on the energy emitted in a black hole merger (a calculation originally done by Hawking). We show how this allows us to think of black holes in manners very similar to heat engines, and how these ideas have been used in modern gravitational wave observatories to test general relativity. This allows a research-level topic to be discussed in introductory physics lectures.

Next-generation gravitational-wave detectors such as Cosmic Explorer and the Einstein Telescope will operate 10-40 km Fabry-Perot arm cavities inside vacuum beam tubes. FFT-based paraxial tools treat propagation in free space and therefore do not explicitly enforce beam tube boundary conditions. We introduce a waveguide-like mode description of the optical field that incorporates an imposed beam tube boundary condition and enables an independent benchmark of free-space FFT tools We derive the associated modal-mixing matrices for mirrors and baffles, including a closed-form series for axisymmetric circular apertures. We quantify the strain-equivalent couplings from baffle miscentering and from a localized near-wall tube defect, and show that they are suppressed as baffle density increases. In the relevant regime of densely baffled cavities and small perturbations, beam tube boundary effects are subdominant, which supports the continued use of FFT-based codes to guide the design of 3G detectors.

Can we estimate the direction of arrival (DOA) of a gravitational wave (GW) signal from pulsar timing array observations? The present paper addresses the inverse problem, for which we consider quadrupole moments of pulsar timings due to GWs from a dominant isolated source such as a binary of supermassive black holes over an isotropic stochastic background. Correlations between the quadrupole moments are discussed, where the correlations between pulsar pairs over the full sky are taken into account. The correlations turn out to be in the form of a three-dimensional traceless matrix with rank 2 that can be closely related with a projection tensor for the GW. Thereby, we demonstrate that the rank-2 matrix allows to estimate the DOA of the GW. In expectation of the forthcoming Square Kilometer Array, angular resolutions as well as DOA estimation errors are also examined.

We consider a massive Kalb-Ramond field with a general non-minimal coupling to gravity. We first study the theory in flat space-time, taking into account the non-linearities. We show that the coupling with the Ricci scalar gives rise to the strong coupling of the two transverse pseudo-vector degrees of freedom, which are absent in the massless theory. We then show that if the theory is instead coupled to the Ricci tensor or the Riemann tensor, the two tensor modes become strongly coupled in addition to the transverse pseudo-vector modes. We then extend our analysis to homogeneous and isotropic space-time, with vanishing background value of the Kalb-Ramond field. We show that in this case, the couplings with the Ricci and Riemann tensor give rise to the runaway instability. Finally, we discuss the inclusion of the disformal coupling as a possible resolution to this unnatural behavior.

We present a new analytic model of force-free electromagnetic jet launched from a disk-fed rotating black hole. The jet solution is obtained through a systematic construction from previously developed methods. The resulting physical jet solution exhibits an asymptotically parabolic structure and is parametrized by the location of localized current concentration and sign reversal in the disk. We find, however, that the jet properties show negligible dependence on the disk parameter. The black hole jet captures the basic feature of the Blandford-Znajek mechanism for energy extraction and jet formation.

We consider a multidimensional polychromatic radiative transfer (RT) problem, accounting for scattering processes in a general form, i.e. anisotropic (dipole) scattering with partial frequency redistribution. Given a discrete ordinates discretization, we report the corresponding matrix structures, depending on model and discretization parameters. Despite the possibly dense nature of these matrices, the use of Krylov methods is effective (especially in the matrix-free context) and robust. We propose a theoretical analysis, using the spectral tools of the symbol theory, explaining why Krylov convergence is robust w.r.t. all the discretization parameters, even in the unpreconditioned case. In fact, the compactness of the continuous operators used in the modeling leads to zero-clustered dense matrix sequences plus identity, so that the clustering at the unity of the spectra is deduced. Numerical experiments confirm the theoretical results, which have a direct application, for example, in the simulation of radiative transfer in stellar atmospheres, a key problem in astrophysical research. In general, we demonstrate that optimal scaling with respect to RT discretization parameters is expected for Krylov solution strategies.

We investigate the modifications to inflationary observables that arise when adopting an $\alpha$-vacuum instead of the standard Bunch--Davies vacuum for quantum fluctuations during inflation. Within the Starobinsky inflationary model, we compute and compare the scalar spectral index, its running, and the running of the running arising from different choices of the initial vacuum state. We further examine the energy scales associated with $\alpha$-vacua and argue that, for any number of extra spatial dimensions, the relevant scale can be truncated at the Hubble scale, $\sim$$\mathcal{O}(10^{13})\,\mathrm{GeV}$, without conflict with current Cavendish-type experimental bounds on sub-millimeter gravity ($\sim$$250\,\mu\mathrm{m}$). Our analysis demonstrates that the $\alpha$-vacuum is subject to stringent constraints as a viable de~Sitter-invariant alternative to the Euclidean (Bunch--Davies) vacuum, with the corrections that it induces in the inflationary observables being strongly limited by the latest Planck data.

Assuming that dark matter is an ultralight pseudoscalar particle which couples to electromagnetism like an axion (an ALP), we demonstrate that the coupling of the cosmological magnetic field produced by the ALP field oscillations to the primordial dark matter fluctuations yields a spectrum of gauge field fluctuations which can produce a sufficient flux of Lyman-Werner photons to enable the Direct Collapse Black Hole formation scenario. The induced flux is consistent with the bounds on the excess flux of radio photons from ARCADE2 and EDGES measurements.