Papers for Friday, May 01 2026

HI-rich starless halos, should they exist, hold great promise for elucidating properties of dark matter halos. This Letter examines the properties of HI-rich failed halos at redshift zero across state-of-the-art cosmological simulations (FIREbox, NIVARIA-LG and Recal-EAGLE). First we compare two numerical analogs with Cloud-9, purported to be the first discovery of a starless HI-rich halo. We argue that differences may be driven by environmental factors, and/or the treatment of gas self-shielding -- which might further limit existing analytic schemes aimed at inferring dark matter halo information from 21 cm HI observations. We also find that the failed halo samples in the three simulations span different regions of the HI-gas-halo mass ($M_{\rm HI}-M_{\rm gas}-M_{\rm 200}$) plane. FIREbox objects occupy a very narrow regime, while NIVARIA-LG extends to a wider range of $M_{\rm 200}$ values - and achieves higher $M_{\rm HI}$ and $M_{\rm gas}$ values. Recal-EAGLE $M_{\rm HI}$ values are similar to FIREbox, albeit with lower gas and halo masses. Lastly, we predict that more HI-rich starless halos can be discovered by exploring the HI-poor regime in the local universe, rather than HI-rich populations at high redshift. Overall, we advocate for the allocation of resources to detect and characterize other HI-rich (and HI-poor) failed halos in the local universe, plus dedicated follow-up spectroscopic observations that scrutinize claims to the absence of a faint stellar component, and that assess their isolation status in detail.

Lopsidedness is common in disk galaxies, yet its origin and evolution remain unclear. Previous studies typically examined stellar and gas asymmetries separately, but a combined analysis offers a stronger probe of the mechanisms driving lopsidedness, recent galaxy evolution, and environment. We analyze the density and kinematics of stellar and atomic hydrogen (HI) components in nine Milky Way type galaxies from the Auriga Superstars cosmological zoom-in simulations. The high stellar mass resolution improves the visibility of disk features while reducing noise, enabling a detailed study of dynamical processes in a cosmological context. Morphological and kinematical lopsidedness are quantified using the first Fourier mode (m=1) of the face-on mass distribution and radial velocity maps, measured consistently for stars and gas between 0.5 and 1 stellar optical radius. At z=0, morphological lopsidedness in old stars (>0.5Gyr) strongly correlates with HI, tracing distortions in the global gravitational potential. In contrast, young stars (<0.5Gyr) trace asymmetric star formation along spiral arms. Stellar morphological and kinematical lopsidedness are strongly correlated, whereas HI shows a weaker correlation, with kinematic asymmetries dominating. We also find an anti-correlation between stellar lopsidedness and bar strength. Strongly barred galaxies tend to host more symmetric disks and higher central stellar mass densities. Tracing lopsidedness evolution over time, tidal interactions with massive satellites (mass ratio >1:50) induce coherent lopsidedness in both stars and HI. In contrast, smooth gas accretion mainly affects HI and young stars, leaving the total stellar component largely symmetric. Overall, these results demonstrate that lopsidedness is a powerful diagnostic of internal disk evolution, gas accretion, and environmental interactions across cosmic time.

IRAS 13224-3809 is one of the most intensively studied narrow-line Seyfert 1 galaxies, with a rich literature reporting diverse and sometimes contrasting interpretations of its complex X-ray spectra and variability. Notably, a fast and variable ultra-fast outflow (UFO) was discovered in this source, sparking debate over its nature and driving mechanisms. Motivated by these open questions, we present a systematic, time- and flux-resolved reanalysis of the full 2016 XMM-Newton (1.5 Ms) and NuSTAR (500 ks) datasets, employing careful background treatment and equal-count spectral selections. We uniformly apply three spectral models, including photo-ionized absorption, broad emission, and relativistic reflection, to all intervals. We unambiguously confirm the presence of a strong, variable outflow with velocities exceeding 0.2$c$, and find that models including absorption consistently reveal robust physical trends: a velocity-luminosity correlation of the UFO, persistently large line widths, and no compelling equivalent-width-flux anti-correlation. When emission or reflection components are included, the significance of the absorption features decreases, but significant UFO detections remain in most intervals. We also report clear evidence for rapid acceleration of the wind in response to X-ray flares, with the outflow carrying momentum and kinetic power sufficient to drive an efficient AGN feedback. The observed rapid response favors magnetic driving, analogous to coronal mass ejections, over radiative acceleration. Our results reconcile contrasting previous claims and underline the need for high-resolution spectroscopy to resolve the wind substructure. The observed UFO variability and structure are consistent with a multiphase, clumpy wind produced by thermal and hydrodynamic instabilities, with magnetic reconnection providing the rapid acceleration mechanism.

We present a detailed study of radio-detected dwarf galaxies (with stellar masses less than 3 billion solar masses) to characterize extreme star formation and search for (variable) radio AGNs. Our sample comes from Reines et al. (2020) (arXiv:1909.04670), who used the Karl G. Jansky Very Large Array (VLA) with 0.25 arcsecond resolution to observe 111 dwarf galaxies with lower-resolution (5 arcsecond) detections in the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) survey. While that work identified and focused on 13 compact radio AGN candidates in dwarf galaxies, here we focus on 16 compact radio sources consistent with star formation in dwarf galaxies. We find that these objects are dominated by thermal HII regions with ages less than 10 Myr, and the most extreme sources have ionizing luminosities requiring the equivalent of around 10,000 to 100,000 O-type stars. We also investigate the dwarf galaxies detected in FIRST but not detected in the high-resolution follow-up observations. Using the infrared-radio correlation parameter, we identify eight sources consistent with radio-excess AGNs. Five of these objects plus another 15 dwarf galaxies have no corresponding detections in the VLA Sky Survey (VLASS) indicating variability between the FIRST and VLASS observations. The FIRST radio luminosities of these sources are significantly higher than expected for supernova-related emission, suggesting the radio variability is likely associated with AGNs. Together, these results provide new context for the presence of compact star formation and massive black holes in dwarf galaxies, and highlight the utility of radio variability and multi-resolution data for identifying the dominant power sources in low-mass systems.

New observations from the Gaia spacecraft have traced an emerging demographic of low-mass associations disconnected from larger associations or GMCs. The first of these associations were recently characterized, but the star-forming environments they trace remain unknown. Using new velocities and ages alongside literature catalogs, we uncover the origins of 16 low-mass associations ($M\lesssim100$ M$_{\odot}$, $\tau\lesssim50$ Myr) using dynamical traceback. We reveal that three groups of currently disparate populations share common formation sites, comprising the Leo, CaNMoS, and AquENS associations. Twelve of 16 associations have plausible connections to larger complexes, six of which form while moving outward from well-established multi-generational star-forming events that drive known or suspected bubbles. We find that feedback from the oldest co-spatial and co-moving relatives of these associations can explain the current morphologies of the Local and Orion-Eridanus Bubbles, along with the formation of related associations like Sco-Cen and Orion OB1. Most remaining populations show evidence for triggered star formation. In the Leo Association, high vertical velocities and a deceleration signature suggest that it formed out of an intermediate velocity cloud colliding with gas in Orion, which would make it the first known case of star formation in one of these clouds. The other newly defined associations show similar asymmetric velocity signatures, such as CaNMoS, which may trace bubble-driven acceleration or a cloud collision. We conclude that the lowest-mass young associations remain undiscovered, and that these populations may have a critical role revealing the small gas overdensities that trace the processes sculpting galactic star formation.

Measuring the properties of cold molecular gas available for intense star formation in galaxy protoclusters at $z>2$ is a crucial step in understanding large scale structure formation. We present ALMA observations of CO(3$-$2) in five dusty star-forming galaxies within $\sim0.5-4$ cMpc of the core of the Spiderweb protocluster at $z=2.16$ to measure the molecular gas mass and kinematics in the most starbursting members of the protocluster. All five galaxies exhibit evidence for disturbed kinematics including non-Gaussian CO line profiles, irregular spatial morphology, and strong residuals when fitting the galaxies with a classical disk model. This could be indicative of an elevated merger rate in the outskirts of the mature Spiderweb protocluster, as all of the galaxies in our sample have multiple companions detected in H$\alpha$. Both the gas fractions and the gas depletion timescales of the galaxies are similar to field relations at cosmic noon, indicative of the fact that their prodigious star formation rates are compensated by similarly high gas masses. The most massive galaxies, as well as all of the galaxies identified as X-ray AGN in previous works, have gas fractions $<30$%, compared to the sample average of 49%, indicating declining availability of gas for star formation. Finally, we find that the gas fractions and specific star formation rates decline with distance from the Spiderweb Galaxy, supporting the reversal of the SFR density--radius relation in high-redshift protoclusters.

A high-resolution ($\sim$45000), high signal-to-noise ($>$100) K-band spectral atlas of massive stars is presented. It includes 81 stars consisting of known optical standards, spanning spectral and luminosity subclasses from O2 to O9, and supergiant luminosity and spectral subclasses from O2-B1. The telluric-corrected reduced spectra are publicly available, and are discussed here.

Solar analogs, stars that closely match the fundamental properties of the Sun, provide key benchmarks for testing stellar structure and evolution across different ages and activity levels. Their detailed characterization helps place the Sun in context within the broader population of solar-like stars. This study presents the characterization of six seismic solar analogs observed by the NASA Kepler and K2 missions. Combining asteroseismic constraints from space-based photometry with high-resolution spectroscopy and \textit{Gaia} astrometry, we derived their fundamental parameters and assessed their resemblance to the Sun. Global seismic properties and individual oscillation modes were extracted from the photometric light curves, while atmospheric parameters were obtained from data collected by the HERMES spectrograph at the Mercator telescope. Stellar modeling using seven independent stellar evolution codes yielded consistent masses, radii, and ages. These stars have masses between 0.91 and 1.04~$\mathrm{M}_\odot$, radii between 0.95 and 1.08~$\mathrm{R}_\odot$, and ages from about 1.8 to 9.1~Gyr, with typical systematic uncertainties of $\pm$ 0.02~$\mathrm{M}_\odot$, $\pm$ 0.01~$\mathrm{R}_\odot$, and $\pm$ 0.7~Gyr, respectively. One star, EPIC~206064678, exhibits properties very similar to those of the Sun, with $M = 1.016 \pm 0.033\,\mathrm{M}_\odot$, $R = 0.990 \pm 0.011\,\mathrm{R}_\odot$, and an age of $5.40 \pm 0.12$\,Gyr. It can therefore be considered a close solar twin, although it is slightly older and more metal-rich ($0.25 \pm 0.07$\,dex). Four targets display binarity signatures and all exhibit very low chromospheric activity. This work broadens the sample of well-characterized seismic solar analogs towards a larger sample of metallicities and ages, providing new references for comparative stellar studies and future asteroseismic investigations.

Understanding the role of galactic scale bars in disk galaxy evolution requires detailed measurements of bar properties across galaxies hosting bars at many redshifts. We present measurements of bar lengths and widths in a sample of 8230 disk galaxies from Hubble Space Telescope (HST) Legacy surveys. The highest-redshift barred galaxies in the sample have $z \sim 3$; most have $z \leq 1$. Using a mass-complete sample from the COSMOS field, we examine bar properties and evolution within $0.25 < z < 1$ in galaxies with stellar mass $\log(M_{\ast}/M_{\odot}) \geq 9.5$. The lowest-mass galaxies in our sample have similar star formation rate (SFR) distributions whether or not they host bars. For galaxies with $\log(M_{\ast}/M_{\odot}) \geq 10$, barred galaxies are more likely to be quiescent or quenched, consistent with bars mainly participating in slow quenching processes. The median physical bar length increases with increasing stellar mass. Relative bar lengths and widths (as a fraction of disk radius) peak at stellar mass $\log(M_{\ast}/M_{\odot}) \sim 10.25$, and change together with mass such that the median ratio, a proxy for bar strength, does not significantly change with stellar mass. Bars in our sample tend to be slightly ($\approx 13$%) weaker at higher redshift. Quiescent and quenching galaxies have longer and wider bars than those in galaxies on or above the star-forming sequence, especially at lower redshift and higher masses; at the low-mass end of our sample, starburst galaxies host relatively longer and stronger bars. Our findings are consistent with other results from studies at both higher and lower redshift, cementing the fundamental importance of bars in disk galaxy evolution.

The early universe need not be described by an incomplete inflationary phase connected to a separate, more exotic prehistory. Recent results show that, within non-static FRW cosmology, only positive spatial curvature permits a nonsingular, geodesically complete universe with ANEC-respecting matter. We construct a geodesically complete closed $k=+1$ bounce-plus-inflation cosmology in ordinary general relativity, sourced by a single canonical scalar field with a positive vacuum offset. The bounce is supported by curvature rather than exotic stress energy: the matter content satisfies the NEC throughout and violates only the strong energy condition, as in any accelerated expansion. The solved branch remains sub-Planckian and evolves onto a curvature-diluted slow-roll phase with inflationary observables consistent with current constraints. The pivot-scale predictions are $n_s=0.9617$, $r=0.0045$ at $N_*=55$ and $n_s=0.9650$, $r=0.0037$ at $N_*=60$. Direct evolution of closed-universe infrared perturbations shows regular tensor and scalar propagation through the bounce and inflationary era, with the physical curvature perturbation freezing in the standard way. This gives a minimal explicit realization of a complete early-universe cosmology in the closed FRW branch selected by completeness and ANEC compatibility.

The evolutionary state of the 198 day eclipsing binary BM Cas is examined using spectra that cover five orbital cycles. Radial velocities measured from SiII 6347 and SiII 6371 track the motion of the primary, and a mass function is found that is similar to that obtained by Popper(1977) from MgII 4481. Absorption from a circumsystem shell complicates efforts to measure stellar velocities from FeII lines. Many of the characteristics of Halpha emission and absorption that are associated with the shell vary in sync with the motion of the primary, and it is suggested that the shell may form from material that exits the system from L2. The infrared spectral-energy distribution departs from that of an A supergiant only at wavelengths > 5um, and models are examined in which the secondary is obscured by an opaque envelope. Archived V band photometry is compared with model light curves, and it is concluded that the A supergiant is close to filling, or is filling, its Roche lobe, and that the as-yet undetected secondary may be more massive than the primary. Based on the overall properties of BM Cas and its environment, we suggest that it is an Algol or post-Algol system, in which the A supergiant was originally the more massive component. If this is the case then the stars in BM Cas had intermediate initial masses. Some of the photometric characteristics of the primary are indicative of alpha Cyg-type variability.

We investigate the survival of dust trapped in magnetically confined cool gas clouds (or {\it prominences}) around rapidly rotating M-dwarfs exhibiting the ``scallop-shell'' light-curve morphology. Using a two-dimensional magnetohydrodynamic simulation, we extend previous coronal prominence models to include a passive tracer field to allow for a single injection of collisionally charged dust grains. The tracer evolution reveals how recurrent centrifugal breakouts--the slingshot process--remove dust and gas from the prominence while chromospheric evaporation replenishes gas from below. For our simulated star, which has $R_{\ast} = 0.6 R_{\odot}$, $M_{\ast} = 0.3 M_{\odot}$, and $P_{\ast} = 0.32$ days, the resulting dust content decays exponentially with a minimum half-life of approximately 6 stellar rotations, representing a lower limit set by our assumption of fully coupled dust and gas dynamics. Synthetic velocity-phase diagnostics show a single, phase-locked feature that fades steadily, reproducing the behaviour of dips seen in TESS and K2 light curves. Comparison with observed river plots suggests a natural classification: (i) persistent, non-decaying features formed by quiescent prominences below co-rotation; (ii) gradually fading features produced by slingshot prominences near co-rotation; and (iii) abrupt disappearances linked to magnetic reconnection and flare-driven ejections. These results demonstrate that dust-bearing prominences--undergoing repeated slingshots--can persist for tens of rotations, linking the observed longevity of the scallop-shell photometric features with the dynamic cycle of prominence slingshot ejections.

Giant radio halos are the Mpc-scale extended sources associated with the merging clusters, while the mini-halos are preferentially associated with cool-core clusters. Both trace the ICM plasma physical process, and recent low-frequency observations increasingly blur the distinction between the two classes. We present the first multi-frequency spectral analysis of the galaxy cluster RXCJ0232--4420, hosting a cool core, using uGMRT (400 and 650 MHz) and MeerKAT (1283 MHz) observations. The central radio emission extends beyond $\sim 1$ Mpc at all frequencies, confirming it as a giant radio halo. One candidate relic (in the east) has also been detected, with an extent of $\sim 300$ kpc. The integrated spectral indices of halo and candidate east relic are $\alpha = -1.17 \pm 0.17$, and $\alpha = -0.85 \pm 0.17$, respectively. The resolved spectral map of the halo is mostly uniform ($-1.0$ to $-1.3$) and does not show any radial steepening. The radio surface brightness profile is well modelled by a single exponential law, with the e-folding radius constant across frequencies. The radio halo emission is morphologically well correlated with the thermal emission. Point-to-point radio-X-ray correlation analysis gives a sublinear relationship (slope $\sim 0.80$), with no frequency evolution. The presence of Mpc-scale emission in the cool-core cluster shows that such emission can arise in dynamically intermediate systems. Our results demonstrate that merger-driven turbulence, even from minor disturbances, can sustain cluster-wide particle re-acceleration without destroying the cool core.

Recent observations of protostars with the James Webb Space Telescope have revealed unprecedented chemical complexity from their ice absorption features. However, these spectra are likely influenced by radiative transfer effects, and there is little understanding of how this impacts our ability to identify, quantify, and interpret the observed ice features. We have developed a new modeling framework to investigate the radiative transfer through icy protostellar envelopes, and apply this to the IRAS 15398-3359 protostar observed by the JWST CORINOS program. The modeled H$_2$O and CO column densities are similar to previous empirical studies, but we require a high CO$_2$/H$_2$O ratio of 76% to match the optical depth of the 15 $\mu$m band. We use our modeled continuum to calculate a 6-10 $\mu$m optical depth spectrum, and see considerable differences compared to a simple polynomial continuum model, underscoring the challenges with quantifying trace ice species in this range. For this source, we find that the observed absorption predominantly originates along the viewing line of sight between 1000 - 2000 au, peaking at the transition from the outflow cavity to the envelope; the spectra are largely insensitive to absorption from ices in the outer envelope, which extends out to 20,000 au. Lastly, we show that depending on how the line of sight intersects the cavity, the apparent CO$_2$/H$_2$O and CO/H$_2$O column density ratios can be underestimated compared to the underlying ice abundance ratios. Together this provides important context for interpreting the ice constraints derived from JWST observations of protostars.

We present Northern Extended Millimeter Array (NOEMA) observations of the CO (2-1) molecular gas kinematics in the nearby Compton-thick Seyfert 2 galaxy NGC 3079, with an angular resolution of 0.5" ($\sim$40 pc). To interpret the observed CO (2-1) kinematics, we model the rotating disk using two software tools, 3D-Barolo and DysmalPy, to generate mock 3D data cubes. Both models indicate, in addition to the rotating disk, the presence of a spatially unresolved nuclear component characterized by high velocity dispersion. Analysis of the visibility data reveals that the blue-shifted, high-velocity component is spatially offset from the continuum peak by 0.17" ($\sim$ 14 pc) and exhibits line-of-sight velocities of $v$ - $v_{sys}$ = -350 to -450 km s$^{-1}$, which we interpret as a nuclear molecular outflow. We calculate a molecular gas mass outflow rate of 8.82 $M_\odot$ yr$^{-1}$, with a kinetic power ($\dot{E}_{\text{out}}$) of 3.8 $\times$ 10$^{41}$ erg s$^{-1}$ and a momentum rate ($\dot{p}_{\text{out}}$) of 2.05 $\times$ 10$^{34}$ Dyne. The momentum rate exceeds the AGN radiation momentum rate by a factor of $\sim$15, suggesting an energy-driven outflow. Furthermore, we argue that the derived kinetic power of the nuclear molecular outflow favors a jet-powered scenario that explains the slowdown and brightening of the parsec-scale radio source observed with the Very Long Baseline Array.

XRISM observations to date have shown that gas kinetic pressures in the intracluster medium (ICM) tend towards the low end of predictions from cosmological simulations. Here, we present a XRISM observation of the merging cluster Abell 2034, which exhibits the broadest emission lines yet observed in a galaxy cluster. We measure a velocity dispersion of ~470 km/s, corresponding to a kinetic pressure fraction of ~15%. This places A2034 at or above the high end of the theoretical predictions for similar-mass clusters. This large velocity dispersion may reflect Mach ~0.5 turbulence in the ICM and/or result from a core disruption driven by the ongoing head-on merger. We also detect a ~380 km/s gas bulk velocity gradient along the merger axis with an opposite sign to the galaxy velocity gradient, indicating a decoupling of the cluster galaxies (and dark matter) from the ICM. Finally, we report tentative evidence of several spectral anomalies, including a suppressed Fe He$\alpha$-z line, an enhanced Fe Ly$\alpha$-2 line, and a potential absorption feature at ~8.7 keV. The first two features may be explained by the combination of a multi-phase ICM and a non-equilibrium ionization state in the wake of a merger shock. Deeper XRISM observations of this cluster are required to confirm these features. This work highlights the importance of kinematic measurements across a large sample of merging clusters as well as the need for deep XRISM observations to unveil more exotic physics in the ICM.

We report the discovery of six new variable stars identified through an exploratory analysis of several sky fields observed by the author using a small telescope and a CMOS camera. The search employed simultaneous photometry of hundreds of objects with AstroImageJ, supplemented by custom Python-based tools developed by the author to generate aperture lists and visualize light curves across the full field of view. The variables were further characterized using TESS data, and their properties are presented. We also describe the data analysis workflow, including the software packages LCV and VS-fit, developed by the author for periodogram analysis and light-curve modelling.

The detection of CO2 on the Jovian satellite Europa by Galileo NIMS and recent mapping of the leading side by JWST has revealed that it is most concentrated in geologically young terrains, and its v3 asymmetric stretch appears as a spectral doublet centered at 4.25 and 4.27 um. Since crystalline CO2 is unstable at Europan surface conditions, this observation implies an active source and a trapping medium, which may be separate. To this end, several hypotheses have been proposed, but no laboratory work has successfully reproduced the spectral features of CO2 on Europa so far. Radiolyzed carbonates have also been discussed as plausible precursors and host materials for CO2, though their role has not been experimentally validated in a Europa-like environment. Here, we report the first laboratory experiments investigating CO2 production from carbonate salts exposed to 10 keV electron irradiation at 50, 100, and 120 K in ultrahigh vacuum. Using diffuse reflectance FTIR spectroscopy, we observe the emergence, growth, and saturation of an absorption doublet centered near 4.25 and 4.27 um, consistent with the CO2 v3 band. Postirradiation thermal desorption studies using residual gas analysis reveal that the radiolytically formed CO2 is stable at temperatures beyond Europa's surface. This work provides the first experimental evidence that low-energy electron irradiation of carbonates in cryogenic, vacuum conditions can produce and retain CO2, and suggests that carbonates can serve as endogenous reservoirs of CO2 on irradiated icy bodies in the outer solar system.

A merger of clusters naturally drives shocks with Mach number $\mathscr{M}\lesssim 3$ in the intra-cluster medium (ICM). This process creates several distinct signatures, including sharp surface brightness "edges", temperature, and gas velocity jumps. The low density of the ICM implies that the ionization balance and electron-ion equilibration times can be long enough to produce a set of additional observable signatures. Here, we focus on two "transient" spectral signatures accessible with the high-energy-resolution telescopes such as XRISM, even for unfavorable geometry, e.g., when we are looking inside the Mach cone of the shock, precluding the appearance of sharp edges in X-ray images. In this work, we focus on (i) the $\mathtt{Z/W}$ line ratio of the Fe~XXV triplet and (ii) the contribution of ions with $T_i>T_{\rm e}$ to the line width, which might be mistakenly interpreted as the gas turbulence. We demonstrate that the $\mathtt{Z/W}$ ratio can serve as a proxy for the non-equilibrium state of the shocked ICM and facilitate interpretation of the line broadening. We conclude that these spectral signatures are within reach with missions like XRISM and can be used to constrain the heating of electrons at the collisionless cluster shocks, as well as the rate of subsequent temperature equilibration between different particle species.

The death of massive stars produces central accreting compact objects and sometimes relativistic jets. Not all jets escape the stellar envelope: unsuccessful, or choked, jets dissipate their energy into a pressurized cocoon, which expands and may break out as a mildly relativistic outflow. We investigate the plasma physics of collapsing massive stars hosting choked jets through relativistic, non-resistive magnetohydrodynamical simulations. We delineate the parameter space for jet choking and quantify the acceleration rate and efficiency of charged particles at strong shocks, which are potential sources of high-energy neutrinos and electromagnetic transients. Our study focuses on blue and red supergiant progenitors, both promising candidates for jet choking.

In this work, we aim to advance the Universal Rotation Curve (URC) paradigm by leveraging new data and extending its observational domain. Building on previous studies that established the URC using optical rotation curves reaching the galaxy optical radius $R_{opt}$, we exploit the SPARC sample's extended HI rotation curves to construct the URC out to $2R_{opt}$. This crucial extension enables us to investigate the mass distribution of spiral galaxies in a region dominated by dark matter an important step to better constrain galaxy mass models and to explore the nature of the dark matter particle. We find that the URC constructed from the SPARC sample's extended HI rotation curves maintains its universal character out to $2R_{opt}$, with the double-normalized rotation curves collapsing onto a single profile. This extended URC provides new insights into the interplay between baryonic matter and dark matter in shaping galaxy rotation curves, particularly in the outer regions where dark matter dominates. Our results not only reinforce the URC paradigm but also refine our understanding of the mass distribution in spiral galaxies, offering new constraints on galaxy mass models and implications for the nature of dark matter.

Massive stars (> 8 $M_\odot$) are known to have high degrees of multiplicity, e.g., with about 60% in triples or higher-order multiples. Such high levels of multiplicity may arise during formation (primary multiplicity) or through dynamical processing of already formed stars in dense clusters (secondary multiplicity). The level of primary multiplicity is an important metric to help distinguish between different formation scenarios, such as core accretion and competitive accretion. The level of secondary multiplicity is expected to evolve with time and be sensitive to local cluster environment. Here we analyze a suite of $N$-body simulations to study bound multiplicity and local projected stellar density, $N_*$, around massive stars within gradually forming star clusters with 50% primordial binaries in the Turbulent Clump Core Accretion (TCCA) paradigm. We find that massive stars rapidly gather triple or higher-order bound companions and enhancements in local $N_*$ via dynamical processes. We study these metrics as a function of environment in a given cluster, quantifying the increasing multiplicity that arises towards cluster centers. We find that secondary multiplicity tends to decrease in more massive clusters due to their higher velocity dispersions, but rises as the mean density of the bound cluster increases. We find our $N_*$ radial profiles are shallower compared to those in the STARFORGE simulations, which form massive stars via competitive accretion. A comparison to the AFGL 5180 system suggests it is better described by TCCA models. However, a larger number of observed systems is needed to better discriminate between these formation models.

The Solar System Notification Alert Processing System, SNAPS, is a downstream broker that ingests moving object data from ZTF and LSST and serves these data and derived properties to the public. This document describes how users can access our SNAPS data and products. This is intended to be a living document that will be updated on the arXiv when significant improvements are made to our data access schemes, and will therefore always contain the most up to date information about interacting with our databases and infrastructure. This is version 1.0.

Most galaxies have supermassive black holes (SMBH) at their centres, surrounded by stars with binary systems also present in this environment. We use two schemes - post-Newtonian (PN) and a scalar perturbation to a background metric to numerically solve the three-body problem of a binary with a SMBH. We test three different PN formulations for the PN scheme: The Einstein-Infeld-Hoffman equation, pair-wise implementation of two-body PN-terms for three bodies and the Arnowitt-Deser-Misner Hamiltonian. We compare these approaches for one million solar mass and one billion solar mass black holes, and find a statistical match between the two approximations for stellar mass binary interacting with a million solar mass black hole. We also perform a statistical study for encounters with this black hole, and find that the higher order PN formulation matches with metric-with-perturbation scheme. However, we find a decrease in separation of the binary, and eccentricity variations between different schemes around the billion solar mass black hole. This behaviour is not present if binary has a large separation or is further away from the black hole due to decreased general-relativistic effects. We find that the pair-wise PN method results in a decrease in separation at pericentre in all test cases irrespective of the distance from the black hole or mass of the black hole, making this the least reliable method for solving this problem. Our work highlights the need for caution when interpreting the results in different formulations around SMBHs. This also shows that when understanding extreme mass ratio inspirals (EMRIs) using simulations, one should beware as the binary gets closer to the black hole.

We report the indication of a long-term quasi-periodic oscillation (QPO) in the $\gamma$-ray emission of the BL Lac object 4FGL J1037.7+5711 (GB6 J1037+5711) using more than 17 years of monthly binned Fermi-LAT observations. Since blazar $\gamma$-ray variability is typically dominated by stochastic red-noise processes from turbulent jet activity and accretion fluctuations, we applied multiple independent timing techniques to test periodic modulation. These include the Lomb--Scargle periodogram (LSP), weighted wavelet $Z$-transform (WWZ), first-order autoregressive red-noise modeling (REDFIT), and epoch-folding analysis. The LSP reveals a significant periodic signal at $478.74 \pm 17.55$ days, exceeding the $99.99\%$ confidence level from Monte Carlo simulations. The WWZ analysis independently recovers a comparable period of $474.72 \pm 27.24$ days at $99.7\%$ significance, while the REDFIT analysis identifies a similar periodic feature at $481.67 \pm 35.41$ days above the $99\%$ confidence level. The epoch-folding analysis further confirms the same modulation timescale. The small differences in period estimates across methods are expected given the distinct mathematical frameworks and sensitivity functions. The long-term flux distribution of the source is better described by a lognormal profile than a Gaussian, suggesting that the underlying variability arises from multiplicative processes. The indication of a $\sim$478 day QPO, independently confirmed across all four timing techniques, may be associated with the orbital dynamics of a supermassive binary black hole system driving Newtonian jet precession or periodic Doppler-factor modulation, Lense--Thirring precession of the inner accretion disk around a rapidly spinning SMBH, or accretion-driven instabilities at the disk--jet interface amplified through relativistic beaming.

We investigate the physical origin of critical mass, a threshold where many galaxy properties and scaling relations undergo fundamental transitions, using the Horizon Run 5 simulation. Focusing on massive ($M_{\rm tot} \geq 10^{12}{\rm M_\odot}$) central galaxies, we examine the mass-dependent turnover of the stellar-to-total mass ratio (STR) and the physical processes driving it. We decompose STR into the stellar-to-baryon mass ratio ($M_*/M_{\rm bar}$) and baryon retention fraction ($M_{\rm bar}/M_{\rm tot}$) to examine galaxies' ability to retain baryons and convert them into stars. We find that STR evolution is dominated by variation in $M_*/M_{\rm bar}$, which changes by over a factor of three, peaking within a narrow range of $M_{\rm tot} \sim 10^{12.4\text{--}12.7}{\rm M_\odot}$ independent of redshift, while $M_{\rm bar}/M_{\rm tot}$ varies by at most 30%. A redshift-independent critical mass at $M_{\rm tot} \sim 10^{12.5}{\rm M_\odot}$ ($M_* \sim 10^{10.7}{\rm M_\odot}$) arises from the changing nature of gas accretion. At this scale, a dynamically stable hot gas halo develops that suppresses cool gas inflow, reducing in-situ star formation efficiency such that $M_{\rm tot}$ growth exceeds in-situ $M_{*}$ growth. Consequently, the hot gas reservoir grows while $M_{*}$ growth slows, producing upturns in $M_{\rm gas}/M_{\rm tot}$ and $M_{\rm bar}/M_{\rm tot}$ and a downturn in $M_{*}/M_{\rm bar}$ that ultimately drives the STR turnover. We also identify a secondary critical mass at $M_{\rm tot} \approx 10^{11}{\rm M_\odot}$ (or $M_{*} \approx 10^{9\text{--}9.5}{\rm M_\odot}$) where gas retention fraction peaks, above which increasing hot gas fraction gradually suppresses in-situ star formation efficiency.

We report the detection of linear polarization in the radio afterglow of GRB 260310A, representing the first centimeter-wavelength polarization detection of a gamma-ray burst (GRB) afterglow and the first measurement of Faraday rotation in a GRB environment. We detect linearly polarized emission across $11-25$ GHz, with a polarization fraction decreasing monotonically from $(3.18 \pm 0.18)\%$ at 25 GHz to $(0.69 \pm 0.22)\%$ at 11 GHz. Interpreting the radio data as emission from a reverse shock in a structured, relativistic jet, the observed depolarization toward lower frequencies is consistent with suppression by synchrotron self-absorption, while the low observed polarization at high frequencies relative to the theoretical maximum suggests a patchy magnetic field in the jet with a coherence scale, $\theta_{\rm B}\approx10^{-2}$ rad. We identify a frequency-dependent rotation of the polarization angle consistent with Faraday rotation, with a rotation measure of ${\rm RM} = -(8300 \pm 90)~\rm{rad/m^2}$ at the GRB redshift. The magnitude of the rotation measure is consistent with propagation through a dense, magnetized environment, such as a progenitor HII region. These findings demonstrate that GRB afterglows exhibit measurable linear polarization at centimeter wavelengths, and that their polarimetric properties probe both intrinsic jet magnetization and the surrounding medium. Future multi-frequency polarimetric monitoring over timescales of days to weeks will enable detailed studies of the evolution of magnetic field structure and provide new constraints on the role of magnetic fields in GRB afterglows.

We revisit the production of primordial black holes (PBHs) by a U(1) gauge field with a pseudo-scalar coupling to the inflaton. We improve upon the existing literature by working in the homogeneous backreaction regime with numerically computed gauge mode functions, adopting state-of-the-art PBH abundance calculations, and incorporating the uncertainty in the statistics of $\delta\rho$. We find that PBHs can account for all of the dark matter in the asteroidal mass range, even when the inflaton gradient energy density is highly subdominant ($10^{-4}$--$10^{-3}$ of the kinetic energy), supporting the validity of the backreaction scheme. This mechanism also unavoidably generates a stochastic gravitational wave background with an amplitude that will be measured at LISA and that will allow to indirectly discriminate between different statistics of $\delta \rho$.

We analyzed 19 years of $R$-band data of the blazar 3C 454.3 from the Whole Earth Blazar Telescope (WEBT) archive, along with new data from its members and from public archives such as those provided by the Small and Moderate Aperture Research Telescope System (SMARTS) and the Steward Observatory projects to search for quasi-periodic oscillations (QPOs). We detected a QPO of $\sim$ 433 days using Lomb-Scargle periodogram, which lasted from MJD 54980--58450 as detected by the weighted wavelet Z-transform technique, making it one of the most persistent QPOs ever detected in the optical regime. The phase dispersion minimization technique was also performed to further validate this QPO claim. We detected this signal at a global significance of $2.53\sigma$ across all methodologies. To explain the observed QPO, we have considered both models focused on the accretion disk around the super-massive black hole (SMBH), and those based purely on jet emissions. Plausible jet-based models involve a shock moving down the jet in a helical magnetic field, whereas the SMBH models could involve Lense-Thirring effect-induced jet precession or dual jets in a binary SMBH system. We introduce a novel approach to distinguish genuine QPOs from spurious signals arising from annual seasonal gaps, a common limitation of ground-based observations.

We present a model-independent, photometry-only framework for identifying strongly lensed supernovae when multiple images are unresolved and blended into a single point source. Building on the simulation-based methodology of Bag et al. (2021), we apply this approach to real Zwicky Transient Facility (ZTF) data using a validation sample of spectroscopically confirmed Type Ia supernovae. The method models the observed flux as a superposition of two time-shifted components, and Bayesian inference is used to estimate the relative scaling and time delay. Applying this framework to 445 well-converged supernovae, we find that only a single object satisfies the selection criteria when adopting a conservative threshold of $\Delta t \ge 12$ days, corresponding to a false positive fraction of $1/445 \approx 0.22\%$. A laxer threshold of $\Delta t \ge 10$ days yields fourteen objects, for a false positive fraction of $3.15\%$. The method provides a scalable and model-independent first-stage filter for identifying lens-like candidates in large time-domain surveys such as the Rubin Observatory's Legacy Survey of Space and Time (LSST).

The inner regions of protoplanetary disks are promising formation sites of rocky planetesimals. Theoretical studies have proposed a scenario in which thermal ionization activates the magnetorotational instability (MRI) in the hot inner disk, and the resulting pressure maximum at the MRI activation boundary accumulates dust and promotes planetesimal formation. However, the inner disk may be thermally unstable, and the activation boundary can vary in time, potentially preventing the maintenance of a dust trap sustained by a steady pressure maximum. We propose an alternative scenario in which planetesimals form in a thermally unstable inner disk through dust self-accumulation driven by the coevolution of dust and disk temperature. To this end, we perform simulations that simultaneously calculate the non-equilibrium thermal evolution, the gas and dust surface density evolution, dust growth, and planetesimal formation. Our results show that thermal instability triggers cyclic MRI activation and deactivation, during which planetesimals are formed. The MRI is activated in the inner disk, and driven by thermal instability, the active region expands outward and then reverts to an inactive state. Triggered by a local enhancement in the dust surface density, dust undergoes self-accumulation while migrating inward in the MRI-inactive phase, causing planetesimal formation. Once the MRI is reactivated at a smaller radius, the cycle restarts. For a typical accretion rate of $10^{-8}M_{\odot}~{\rm yr^{-1}}$, a planetesimal belt forms near 1 au. This mechanism can produce sufficient planetesimal mass to form multiple super-Earths. This work provides a framework for a self-consistent model of planetesimal formation based on the coevolution of dust and disk temperature, serving as an initial condition for subsequent planet formation simulations.

We present a comprehensive analysis of dense cores and filamentary structures in the M16 Eagle Nebula using high-resolution ($11.7^{\prime\prime}$) surface density and temperature maps derived from \textit{Herschel} observations. Using the \textit{hires} algorithm for map construction and the \textit{getsf} method for source and filament extraction, we identified 233 cores and 111 filaments in this massive star-forming region. The filaments exhibit a median width of 0.4\,pc -- and a median linear density of 61\,$M_\odot$\,pc$^{-1}$, with 76\% being supercritical for gravitational fragmentation. Our radial analysis of the $\sim$60\,pc diameter shell driven by the central NGC 6611 cluster reveals strong enhancements in structure formation: filament formation efficiency (FFE) is 2.3 times higher within the shell (peaking at 22\%), while core density shows a concurrent 1.5-fold enhancement. The moderate correlation between core density and FFE ($r=0.67$) indicates coupled formation processes. Theoretical analysis demonstrates that observed surface densities exceed the critical threshold for fragmentation by a factor of $\sim$8, with a fragmentation timescale ($\sim$1.5--2.0\,Myr) comparable to the shell's dynamical age ($\sim$1.0--1.3\,Myr), indicating we are observing fragmentation in progress. These results reveal a hierarchical fragmentation sequence -- shell compression $\rightarrow$ filament formation $\rightarrow$ core formation -- providing clear observational evidence for positive feedback where massive star formation triggers secondary structure formation in the surrounding molecular cloud.

Recent technological advances in astronomy, particularly the growing popularity of smart telescopes for the general public, make it possible to develop highly effective detection solutions that are accessible to a wide audience, rather than being reserved for major scientific observatories. Published in 2023, DeepSpaceYoloDataset is a collection of annotated images created to train YOLO-based models for detecting Deep Sky Objects, particularly suited for Electronically Assisted Astronomy. In this paper, we present an update to DeepSpaceYoloDataset with the addition of a new split, test2026, designed to evaluate detection models with a greater diversity of images.

Extended gamma-ray sources provide significant information about particle propagation. 4FGL J1626.0-4917 was labeled as an unassociated source in 4FGL catalog without known counterparts at other wavelengths. We report an analysis on 4FGL J1626.0-4917 with 17 years Fermi Large Area Telescope data and archival Chandra X-ray Observatory data. We find extended GeV emission around this source, which can be modelled by a Gaussian disk of 0.28 degree radius with a significance of the extension of 7.2 sigma. The gamma-ray spectrum of 4FGL J1626.0-4917 has a photon index of 2.73. The gas content, including molecular, neutral and ionized gas, was investigated and the potential hadronic origin is discussed. The diffuse GeV gamma-ray emission may likely originate from the interaction between accelerated protons in 4FGL J1626.0-4917 and the target proton in surrounding gas, although the leptonic process cannot be ruled out. The X-ray spectral analysis was performed, which reveal a point source inside 4FGL J1626.0-4917. We investigate potential counterparts, including the stellar cluster NGC 6134 and the supernova remnant G335.2+0.1. Our results highlight the complexity of unidentified extended gamma-ray sources and the need for further observations.

The nature of eclipsing time variations (ETVs) in post-common-envelope binaries (PCEBs) is still unknown. Circumbinary planets routinely fail the test of time and the Applegate mechanism has energetic constraints and problems in reproducing observations. Based on recent analytic models of magnetically-induced ETVs and stellar dynamo simulations, we aim at explaining ETVs via non-axisymmetric magnetic fields that drift in the azimuthal direction of the star, know as azimuthal dynamo waves (ADWs). We implement a time-varying non-axisymmetric quadrupole moment ($Q$) in a binary system. We solve for the dynamics of the system, compute the resulting eclipsing times, and construct O-C, diagrams. We perform several simulations with different amplitudes of $Q$, periods, stellar masses and binary separations. ADWs naturally give rise to characteristic shapes in the O-C diagram that resemble observations. Depending on how fast $Q$ changes, the solutions can have a sharp decrease in O-C producing amplitudes such as the one obtained in QS Vir, or sinusoidal-like shapes such as in V471 Tau or NN Ser. We also find that the amplitude of the eclipsing times varies from tens to hundreds of seconds. ADWs offer a self-consistent explanation for ETVs as they are expected in dynamo theory. They can explain a variety of features in the observed O-C diagrams. As suggested by dynamo simulations, ADWs are easily excited in rapidly rotating stars, alleviating energetic constrains required in the context of the Applegate mechanism. They produce non-axis $Q$ that in turn produce ETVs that can account for the long-term variation of the O-C, diagrams. We expect in this case that the resulting O-C diagrams are not strictly periodic, unlike explanations based on a third body that would imply a strict periodicity unless additional mechanisms are being invoked.

The long-standing assumption of symmetric radio emission beams from the two magnetic poles of pulsars is challenged by observational evidence of asymmetry and underfill. Direct testing of this symmetry remains difficult for most pulsars. As an indirect test, we collected polarization profiles of 11 interpulse pulsars observed with the Five-hundred-meter Aperture Spherical radio Telescope, MeerKAT, and Parkes. We developed a rotating vector model incorporating aberration and retardation effects to fit the position angle swings of selected pulsars, thereby determining the intrinsic emission region corresponding to the observed pulse windows. Based on both the conal and fan beam models, we compared three key parameters-beam radius, magnetic azimuth width, and emission intensity-between the intrinsic emission regions of the main pulse and interpulse. Among the eight pulsars with a confirmed double-pole geometry, none exhibits similarity in the azimuth width. Only two show potentially similar beam radii, while six demonstrate comparable emission intensities within specific parameter spaces. These results indicate that the emission beams from the two magnetic poles of a pulsar may be generally dissimilar in size, suggesting that the physical conditions governing pair production and particle acceleration differ between the two poles. The random distribution of active emission regions further implies inhomogeneity within the polar cap, which may originate from the differences in local magnetic field structure or surface properties.

The blazar 3C 279 is well known for its rapid and large-amplitude variability. On 20 December 2013, the source exhibited an orphan {\gamma}-ray flare characterized by a flux-doubling timescale of a few hours, a very hard spectrum, a time-asymmetric light curve with a slow decay, and no significant optical variability. We propose a new interpretation of this event based on a two-zone scenario in which a stationary emission region produces the quiescent emission, while a second zone accelerates within the broad-line region(BLR). We compute the time-dependent radiative output of both zones with the EMBLEM code, including synchrotron, synchrotron self-Compton, and external inverse-Compton processes, as well as bulk acceleration, adiabatic expansion, and a Fokker-Planck treatment of the electron distribution. This is the first attempt to precisely model the asymmetric {\gamma}-ray flux evolution during this flare. A model with a stationary region outside the dusty torus and an accelerating plasma blob reproduces the main features of the event: a short and intense {\gamma}-ray flare with a hard spectrum and no optical counterpart. The flare results from the variation of the external photon field in the blob frame as the blob crosses the BLR and reaches its terminal Lorentz factor not far from the inner radius of the BLR. Bulk acceleration and the propagation of a plasma blob within the jet provide a natural mechanism for producing high-energy flares and asymmetric light curves without requiring an ad hoc time-dependent particle injection. The model predictsa delayed EUV/X-ray enhancement once the blob exits the BLR. No very-high-energy data are available for this event, but if {\gamma}-rays were emitted in this band, a delay would be expected with respect to the Fermi-LAT flare.

The evolution of protoplanetary disks, especially in the early stages of planetary formation, as dust grows, is the cornerstone of the birth of planets. The mechanisms involved in the growth of sub-micrometric dust grains into planetesimals within a very short time frame are a challenging field of study, while the initial conditions remain relatively undefined. One of the main challenges is to unambiguously identify the dust properties within the disk, and our goal is to break this barrier by investigating the light scattered by dust particles lying on the protoplanetary disk surface from many recent promising observations. In this study, we used a set of 30 polarized light images composed of new VLT/SPHERE observations to examine the light scattered by dust grains. For each ring-shaped system, we used the new DRAGyS tool to estimate the disk geometry using the substructures visible on the surface and to extract the limb-brightening-corrected scattering phase function, which encodes the dust grains' physical properties. Finally, we compared our results with the AggScatVIR database of numerical scattering phase functions of nonspherical dust. We combined our measurements of disk geometry to estimate an average disk flaring of about 1.357. First, we recovered the two categories of scattering phase functions based on their shape, as determined in previous studies. Category I is monotonically decreasing and can be explained by fractal organic aggregates with small monomers of 100nm, or compact aggregates with medium porosity and big monomers of 400nm. Category II is defined by a bell-shaped scattering phase function and can be explained by sub-micrometric irregular grains or compact aggregates with low porosity. This statistical study offers general trends about dust populations, but the degeneracy is too strong to apply this method to a unique disk analysis.

Cygnus X-3 is a microquasar composed of a compact object of unknown nature closely orbiting around a Wolf-Rayet star. The particularities of this source make it a unique case among microquasars. This fact, together with its recent establishment as a PeV particle accelerator, makes Cygnus X-3 a very interesting target for the investigation of the physical processes leading to gamma-ray production. In this work, the TeV and GeV gamma-ray emission of Cygnus X-3 is studied in order to determine its origin and constrain the properties of the system. For that purpose, a point-like analysis of 130 h of data taken with the MAGIC telescopes between 2013 and 2024 was performed, which represents the largest available sample for Cygnus X-3 at $\sim$TeV energies. Additionally, contemporary data from Fermi-LAT were also analysed to better contextualize the MAGIC observations. For a more detailed investigation of the source physics, the data were divided into three subsets according to the flaring state of the source and orbital phase. No significant detection of Cygnus X-3 is found between 0.1 and 7 TeV for any of the datasets, and differential and integral flux upper limits are reported over the long-term monitoring of the source. The Fermi-LAT fluxes can be considered compatible with previous results, taking into account the different data samples used across studies. The MAGIC upper limits presented in this work represent the most constraining ones up to date at $\sim$TeV energies. An eventual detection of Cygnus X-3 at these energies would significantly constrain the source properties, and is not unreasonable to expect given that the source has already been detected in both the GeV and PeV regimes during flaring states. Further observations of Cygnus X-3 at energies above tens of GeV would be valuable for this purpose.

In recent years, radio detection of ultra-high-energy cosmic rays (UHECRs), with energies above $10^{18}$ eV, has become an established technique. The radio emissions can be simulated with high accuracy using Monte Carlo codes such as ZHAireS and CoREAS. These simulations are essential but are computationally intensive. In this work, we present a machine-learning-based emulator that reproduces radio signal simulations with high accuracy in milliseconds rather than hours. Primary particle properties can then be reconstructed by comparing measured signals to emulated traces using a Markov Chain Monte Carlo approach. Using ZHAireS simulations carried out over the GRANDProto300 experiment layout, the method achieves an 8.9\% resolution on electromagnetic energy and a 0.08° angular resolution, matching state-of-the-art reconstruction performance. Finally, we apply the method on real data, successfully reconstructing cosmic-ray candidates detected by the GP300 prototype.

We propose that the size of the universe and its rate of expansion cannot be simultaneously specified with arbitrary precision, a quantum mechanical statement encoded in a deformed commutation relation for the scale factor. The deformation modifies the Friedmann equation by adding a geometric correction to the expansion rate, and the sign and magnitude of a single free exponent determine the cosmological behavior. When the exponent is positive, the model predicts late-time dark energy with $w > -1$, testable with current and next-generation surveys. When the exponent is sufficiently negative, the same deformation produces a non-singular classical bounce that resolves the Big Bang singularity. The model introduces no new particles or fields and preserves a scale-invariant primordial power spectrum. The deformation has a natural interpretation as a horizon-scale phenomenon, with the cosmological horizon, and not the Planck length, setting its characteristic scale. The late-universe regime is then its generic application, with the expansion history as the primary observable signature. Cosmic acceleration may be the macroscopic imprint of quantum gravity at the cosmological horizon.

Unraveling the growth of supermassive black holes and their connection to host galaxies requires disentangling the Active Galactic Nuclei (AGN) emission from that of the stellar populations. When an AGN spectrum is observed at different activity phases, if the spectral decomposition properly recognizes the nuclear and stellar components, key physical properties - such as black-hole mass, stellar mass, and stellar velocity dispersion - should remain consistent. We present a novel optical spectral-fitting approach that combines pPXF and PyQSOFit to robustly decompose spectra into stellar and AGN components. We apply this technique to three SDSS samples with repeated optical spectra of the same objects at z<0.55: 32 changing-look AGN in bright and dim states, and 15 quasars and 15 galaxies with three single-epoch and one stacked spectrum each. To compare with the literature, we use SDSS spectra and photometric data from AGN in the eFEDS field, as well as Gemini and VLT observations of some of our changing-look AGN. We evaluate the reliability of stellar mass, velocity dispersion, and black-hole mass measurements, especially in relation to the AGN-to-total continuum contribution (fAGN). For host-derived properties, especially when fAGN<0.8, our method yields consistent results. For single-epoch black-hole mass estimates from Ha and Hb, 3-sigma confidence in the broad-line flux and FWHM provides effective criteria for selecting reliable measurements. After applying these quality cuts, measurements across different epochs agree within uncertainties, and their reliability is confirmed by the alignment with previously established scaling relations. Many changing-look AGN in our sample exhibit "breathing" broad-line regions, as determined from Ha analysis, while some deviate significantly, suggesting non-virialized systems across the spectral transition.

Open clusters are fundamental laboratories for investigating stellar and Galactic evolution, and serve as important benchmarks for asteroseismic analyses. Using a boutique method to analyze TESS photometry, we study red giants in two old metal-rich open clusters: NGC 188 & NGC 6791. By comparing Kepler and TESS observations for NGC 6791, similar oscillation mode frequencies are recovered, however we find a systematic offset of 2.2% with a scatter of 9% in the $\nu_{\text{max}}$ measurements. We attribute this discrepancy to the lower signal-to-noise of the TESS data for these relatively faint stars. For the brighter cluster NGC 188, we present new seismic measurements in 17 red giants. We estimate average seismic masses for the RGB of $M_{\text{RGB,NGC188}} = 1.13\pm0.04$(rand)$^{+0.12}_{-0.19}$(sys) $M_{\odot}$ and RC of $M_{\text{RC,NGC188}} = 1.11\pm0.01$(rand)$^{+0.11}_{-0.19}$(sys) $M_{\odot}$, consistent with independent mass estimates for this cluster and with similar precision to previous Kepler studies. From the difference between the average evolutionary phase masses, we estimate an integrated RGB mass loss of $\Delta M = 0.02 \pm 0.04$(rand)$\pm0.01$(sys) $M_{\odot}$, supporting the evidence for lower mass loss at higher metallicities. Using asteroseismology and chemical abundances, we identify three binary interaction candidates: two under-massive stars and one over-massive star potentially exhibiting dipole-mode suppression. Finally, we derive an average seismic cluster age of $7.0\pm0.9$ Gyrs, in good agreement with previous literature ages. Our analysis demonstrates the strong potential of TESS asteroseismology for open clusters, and motivates extending this investigation to other TESS clusters that span a wider range of ages and metallicities.

We present detailed asteroseismic modelling of 95 main-sequence solar-like stars and Kepler exoplanet host stars using the FICO procedure, a three-step method that combines forward and inverse techniques that enables precise inference of fundamental stellar parameters such as mass, radius, age, and mean density. We applied the FICO procedure to a catalogue of stars with high-quality asteroseismic and classical observations, and compared its results against literature values. We also compared its performance with direct frequency fitting using semi-empirical surface corrections. The FICO procedure achieved statistical precisions of 2.3%, 0.82%, 6.9%, and 0.49% in mass, radius, age, and mean density, respectively on average, well within PLATO quality requirements. We reconfirmed that surface-independent methods more effectively mitigate biases inherent to semi-empirical surface corrections, particularly for stars more massive than 1.15 Msun or above 6050 K. Two regimes were identified: near-solar conditions, where both approaches perform similarly, and higher-mass stars, where surface-independent methods consistently outperform direct fitting methods. While our results are consistent with literature values, we observed age biases (~11.5% on average for the Kepler LEGACY sample) that are comparable to the PLATO accuracy requirement of 10% for a Sun-like star, and therefore not negligible in that context. The FICO procedure provides a robust framework for high-precision stellar characterisation in the PLATO era. Its hybrid architecture effectively addresses surface effects, making it a promising tool for the accurate determination of exoplanet host-star properties. Our findings also highlight the importance of carefully selecting and validating the physical assumptions embedded in stellar models, particularly in the context of next-generation space missions such as PLATO.

Using the SIMBA, EAGLE, and IllustrisTNG-100 galaxy formation simulations, we examine the anisotropy of the satellite distribution and its dependencies on central galaxies, host halos, and cosmic filaments. We find that in all simulations the satellite anisotropy is robustly aligned with the halo/central galaxy major axis. This correlation is both redshift- and halo-mass-dependent and also extends to filamentary structures outside the halo to several virial radii. The alignment persists up to $z=1.5$ at high redshifts, and the mass dependence remains down to $M_\mathrm{200c} \approx 10^{11}M_{\odot}$. We identify a clear $3\sigma$ scale-dependent transition in the structural tracers of satellite anisotropy: satellite distributions correlate with central galaxy morphology at small scales ($<0.3R_{\rm 200c}$), are governed by host halo triaxiality at halo scales ($0.3$-$2R_{\rm 200c}$), and align with cosmic filaments beyond $2R_{\rm 200c}$. By tracing satellite trajectories in SIMBA, we uncover the kinematic origin of this transition, demonstrating that satellites prefer halo major-axis aligned regions because their trajectories intersect this axis far more frequently and stay in it for a longer time under the host's gravitational potential. This dynamical processing effectively erases primordial filament-related signals upon accretion ($<2R_{\rm 200c}$), explaining the shift in dominant structural tracers across scales.

In this work, we investigate a classical cosmological test - the dependence of galaxy surface brightness on redshift z (the Tolman test). We analyzed 6 860 galaxies with reliably determined spectroscopic redshifts from the ASTRODEEP-JWST photometric catalogue. We find that (a) the mean surface brightness of galaxies indeed decreases with increasing distance, and (b) the observed trend shows a significant departure from the prediction of the standard cosmological model, which expects the mean surface brightness to decline as ~ (1 + z)^-4.

Star-forming galaxies (SFGs) are a subclass of $\gamma$-ray emitters and a correlation between their $\gamma$-ray luminosity ($L_{\rm \gamma}$) and the total infrared (IR) luminosity ($L_{\rm IR}$) has been established based on the Fermi Large Area Telescope (LAT) data. NGC 2403 and NGC 3424 have been reported as outliers in the $L_{\rm \gamma}$-$L_{\rm IR}$ correlation with light curves showing significant variability, which contrasts with the temporally stable $\gamma$-ray emission in other SFGs, originating primarily from cosmic rays interacting with interstellar medium. In this study, we reanalyze the $\gamma$-ray emission in the directions of NGC 2403 and NGC 3424 using more than 16.5 yr Fermi-LAT data. NGC 3424 is found to be spatially coincident with the detected $\gamma$-ray source, while NGC 2403 is significantly offset from the nearest $\gamma$-ray source, suggesting an implausible association. We confirm the previously reported variability of both $\gamma$-ray sources and the significant deviation from the $L_{\rm \gamma}$-$L_{\rm IR}$ correlation when assuming an association of both $\gamma$-ray sources with the two galaxies. Our findings lend further support to the interpretation that their $\gamma$-ray emission is driven primarily by alternative radiative processes-rather than by star formation activity-such as the ejecta of the Type IIP supernova SN 2004dj in NGC 2403 interacting with a surrounding high-density shell and an obscured active galactic nucleus in NGC 3424.

Galaxy physical properties are influenced by their environments, but the processes responsible for mass and environmental quenching and structural transformations remain debated. Galaxy clusters are ideal laboratories for investigating galaxy formation and evolution, offering a full range of galaxy properties and environments. Observations of large-scale structures, particularly filaments in cluster outskirts ($r \sim5r_{200}$), are currently constrained to the low-redshift Universe. To explore galaxy evolution at intermediate redshifts, deep photometric data, ideally combined with spectroscopic redshifts, are essential. Abell S1063 cluster ($z$ = 0.346) is observed within the Galaxy Assembly as a function of the Mass and Environment program with the VLT Survey Telescope (VST-GAME) combined VISTA Public Survey program Galaxy Cluster At Vircam. We investigate galaxy evolution across a wide range of stellar masses and environments. We release a multiwavelength photometric catalog with photometric redshifts for 64394 sources in $1x1 deg^2$. The analysis of overdensity regions provides insights for future studies on galaxy properties in cluster outskirts. The dataset is obtained through deep ($r<$24.65 mag) and wide optical ($u$, $g$, $r$, $i$, VST) and near-infrared ($Y$, $J$, $K_s$, VISTA) observations. The photometric catalog includes all detected sources, excluding nearby or overlapping objects, saturated stars, and image artifacts. The multiwavelength catalog enabled photometric redshift estimates and identification of cluster members. The density field allowed comparison of galaxy properties, colors, and masses across environments. We detect a very dense structure near the cluster center, and with such a large field of view, we find another dense region to the north-west, in the opposite direction to the cluster elongation. Filaments connecting the regions are also visible.

Archival XMM and ROSAT X-ray data are used to investigate the structure of the Abell 548 - Abell 3367 region. Based on previous optical studies, this is a region likely to be rich in structure though studies are in disagreement regarding the connection between Abell 3367 and Abell 548. We use the available archival X-ray data together with kinematic data of counterpart galaxies to address this question and to determine the structure in this region. The region is particularly rich in X-ray structure elongated along a SW-NE axis consisting of numerous extended X-ray sources. In general, the structure consists of many galaxy groups and clusters which appear segregated in X-ray luminosity with the least luminous $\sim$ 30% toward the outer region of the clusters, possibly tracing a filament. We find evidence to suggest a supercluster of 3 clusters at redshifts: $\sim$ 0.04, 0.045, and 0.06. Some of the X-ray sources coincident with Abell 3367 have counterpart galaxy redshifts consistent with Abell 548 and others are significantly higher. This supports that Abell 548 and Abell 3667, form a supercluser and the higher redshift X-ray source is a background object. They are part of a larger structure consisting of a previously identified cluster at redshift 0.04, and two groups at redshift $\sim$ 0.06. In addition, there is a filamentary structure at z $\sim$ 0.103. The ubiquity of groups in the large scale structure suggests that they provide an environment where galaxies are in close proximity and evolution via interaction can proceed well before the galaxies make their way into the dense central region of a cluster.

Context. The Large Array Survey Telescope (LAST) is a wide-field visual-band survey designed to explore the variable and transient sky with high cadence. Its raw data stream is automatically processed in near real time at the observatory site, producing science-quality images, catalogs, and transient alerts. Transient alerts are then reported to the Transient Name Server (TNS). Aims. The LAST pipeline comprises two major components: (i) processing and calibration of single images followed by coaddition of $20\times20$s exposures, producing single-image and coadded-image catalogs; and (ii) subtraction of coadded images from calibrated reference images followed by transient detection. In this work we present a detailed description and validation of the second component of the pipeline. Methods. Transient detection is based on the algorithm for proper image subtraction (ZOGY). We combine ZOGY subtraction with the Translient statistic for sub-pixel motion discrimination, together with a sequence of deterministic filtering steps, to produce a clean stream of transient candidates without the use of machine learning. Results. Using commissioning data, the pipeline achieves a preliminary $5\sigma$ limiting magnitude of $20.3$-$20.7$mag, a single-epoch transient detection efficiency of $\sim80$%, and a purity of $\gtrsim90$% at signal-to-noise ratio of $\geq7.5\sigma$.

We present a highly simplified nonlinear hydrodynamical model to emulate the main observed features of amplitude modulation (commonly known as Blazhko effect) in RR Lyrae stars. The model is based on the assumption that the periodic flow generated by the pulsation carries surplus helium in the ionization zones He I and II. Once this extra helium reaches a critical amount, a Rayleigh-Taylor-type instability leads to a back-flow of the surplus helium and the process starts over again, due to the continuing effect of pumping helium upward by the pulsation. This periodic variation of helium leads to various efficiency of radiation flux blocking in the helium ionization zone that shows up as a long-term periodic variation of the pulsation amplitude.

In this work, we extend the standard pre-inflationary misalignment mechanism for axion-like particles (ALPs) by introducing a pre-oscillatory phase with constant equation of state $\omega_\phi\in[-1,1]$, generated by a tracking potential. During the radiation-dominated era, the potential undergoes a rapid transition to the conventional cosine potential. The resulting change in the potential energy across the transition can drive the ALP into a kinetic misalignment phase ($\omega_\phi=1$) prior to the onset of oscillations. Motivated by persistent cosmological tensions, such as those in $H_0$ and $S_8$, we also investigate an ALP coupling to a dark radiation sector (DR), allowing for its decay. Using a Bayesian analysis, we constrain the ALP parameter space with current cosmological data. Our analysis shows that ALP-induced DR does not resolve the existing tensions. Instead, the data place robust constraints on the model, favoring negative values of $\omega_\phi$ and constraining the symmetry-breaking scale to $f_\phi\in[80,1.5\times10^{10}]~\mathrm{TeV}$, corresponding to ALP masses in the range $m_\phi\in[10^{-20},10^{-2}]~\mathrm{eV}$.

Type Ia supernova (SN Ia) standardisation parameters exhibit evidence for systematic variation across the host galaxy star-formation rate - stellar mass (SFR$-M_\star$) plane, motivating the incorporation of galaxy SFR information in cosmological inference. SFRs are commonly estimated via spectral energy distribution (SED) fitting with far-infrared (FIR) measurements to account for dust-obscured star formation. Such FIR coverage will, however, be limited for upcoming time-domain surveys such as the Rubin Observatory Legacy Survey of Space and Time (LSST), necessitating the use of alternative SFR tracers. Here, we reconstruct the SFR - $M_\star$ plane using 1.4 GHz diagnostics, to test the consistency of host classifications against FIR-constrained SED-based estimates. Within this plane, SN Ia host galaxies are divided into three regions: Region 1 (low-mass), Region 2 (high-mass star-forming) and Region 3 (high-mass passive). We find that ${\sim}84$ per cent of SN hosts retain identical region assignments when using radio versus FIR-constrained SED-derived SFRs. Measuring SN Ia nuisance parameters ($\alpha,\beta, M$) within each subregion, we find consistent values between the two SFR - $M_\star$ plane reconstructions, indicating limited sensitivity to SFR estimator choice, with the largest deviations in Region 3 at ${\sim}1.1\sigma$. Across the three 1.4 GHz SFR - $M_\star$ subregions, we confirm the region-dependent variation in SN Ia standardisation parameters - particularly $\beta$ - reported in our earlier SED-based analysis. With near-complete radio coverage of the LSST footprint anticipated from current and forthcoming radio continuum surveys (e.g., Square Kilometre Array), radio SFR calibrations will become an increasingly useful and scalable approach to host galaxy classification, supporting the construction of robust SN Ia subsamples for precision cosmology.

Despite the diversity and energy dependence of $\gamma$-ray pulse morphologies in Crab, Vela and Dragonfly, the phaseograms of these three canonical young pulsars can be organised within a single four-beam geometric template. Using \textit{Fermi} Large Area Telescope data, we fit the 60~MeV--3~GeV phaseograms with a mechanism-agnostic, geometry-first parametric model that incorporates phase-dependent Doppler shifts and constrains the three-dimensional locations and bulk motions of four emission sites. In each pulsar, the phaseogram admits a decomposition into two altitude-separated beam pairs. The lower-altitude pair is produced by plasma with bulk motion close to azimuthal corotation, sharpening the main peaks. The higher-altitude pair shows a radially outward bulk-motion component, suggestive of inertial effects in a toroidally dominated magnetic field, and contributes bridge/shoulder emission and ripple-like modulations overlapping the main peaks. As a posteriori, the lower-altitude pair is consistent with curvature-dominated outer-magnetospheric emission, while the higher-altitude pair is consistent with synchrotron-dominated emission from a current-sheet-like outflow. Higher-altitude site heights vary from $\simeq 0.7$ (Crab, $\approx 1$~kyr) to $\simeq 1.1$--$1.4$ light-cylinder radii (Vela and Dragonfly, $\approx 10$~kyr). This unified four-beam, observation-driven geometry maps an altitude-dependent azimuthal tilt of pulsed $\gamma$-ray emission, providing an observationally anchored framework amenable to systematic tests and readily extensible to other young pulsars.

Fast radio bursts (FRBs), thought to originate from magnetars, exhibit diverse polarization properties that constrain their emission physics and local magneto-ionic environments. The polarization position angle (PPA) is particularly sensitive to magnetic-field geometry in the emitting region and propagation effects in the magnetosphere and beyond. In hyper-active repeaters, PPAs are typically stable within bursts and over timescales of hours to days. Here, we present observations of the repeating source FRB~20240114A, which show significant burst-to-burst PPA variations. Using full-Stokes, high-time-resolution observations from the Nançay Radio Telescope (1.1--1.8\,GHz) and the Effelsberg 100-m telescope (1.3--1.5\,GHz) over $\sim1$~year, we measure rotation measures (RMs), polarization fractions, and time-resolved PPAs across 12 epochs. The RMs remain stable, and the emission is predominantly highly linearly polarized, with $\sim81\%$ of bursts showing $L/I > 0.8$, while circular polarization is weaker ($\sim16\%$ with $|V/I| > 0.1$). We find no evidence for Faraday conversion. The PPA exhibits rapid, stochastic variations from milliseconds to hours, spanning $\pm90^\circ$ during two active periods and $\pm50^\circ$ in a third. The distribution of PPA jumps shows that (1) there is no difference in the distribution of jumps on timescales shorter or longer than 1\,s; (2) positive and negative jumps are equally likely; and (3) a jump of $\pm90^\circ$, as expected from, e.g., orthogonal mode jumps, is not more common than any other value. This combination of stable RM, high linear polarization, and extreme PPA variability is not seen in other hyper-active repeaters. These results disfavor emission from a single fixed region and instead suggest multiple emission regions and/or strong magnetospheric and foreground propagation effects, such as plasma lensing.

The long-standing tension in the Hubble constant $H_0$ has motivated extensive explorations of both new physics and observational systematics, for example, the late-time systematics in measuring the B-band absolute magnitude $M_B$ of type Ia supernovae, which is degenerated with $H_0$ via an intercept $-5a_B=M_B+5\lg (c/H_0/\mathrm{Mpc})+25$ in the linear relation $m_B=5\lg d_L(z)-5a_B$ between the apparent magnitude $m_B$ and logarithmic dimensionless luminosity distance $\lg d_L(z)$. Therefore, this intercept can be evaluated directly from pure observational quantities ($m_B$ and the redshift $z$) for a given model of $d_L(z)$ without knowing underlying systematics in $M_B$-$H_0$ degeneracy. Hence, the constancy of this intercept across different supernova datasets and different redshift bins within the same dataset for a given late-time model serves as a powerful diagnostic for disentangling late-time new physics from local supernova systematics. In this mini-review, we will show that: (1) there is a local $a_B$ tension in PantheonPlus around $z\sim0.01$, and the elimination of it leads to a $H_0$ measurement consistent with both SH0ES typical three-rung and first two-rung measurements; (2) there is a late-time $a_B$ tension in DES-Y5 around $z\sim0.1$, and the elimination of it largely reduces the preference for dynamical dark energy. We also update the late-time $a_B$-tension analysis for both DES-Y5 and DES-Dovekie supernovae, and find that this $a_B$ tension around $z\sim0.1$ is mainly driven by the inter-data tension between DES supernovae and DESI+Planck constraint, and the dynamical dark energy is preferred as a compromise of this tension. Finally, we briefly mention an interacting dark energy model that resolves this tension among DES, DESI, and Planck, and point out a crucial difference between the effective and apparent equations of state of dark energy.

The WISE source, J141046.00+740511.2, has been recently observed from radio to $\gamma$ rays. Although the optical spectrum is consistent with a BL Lacertae (BL Lac) object, the source displays unusually weak radio emission, which challenges standard interpretations. Our aim is to understand the origin of the broadband emission from J141046.00+740511.2, using a leptonic model of an extended jet. To obtain the distribution of electrons along the conical jet, we solved a steady-state convective transport equation. Emissivities were computed along the jet and integrated over the cone volume to obtain the observed flux. Our model successfully reproduces the observed multiwavelength spectral energy distribution from radio to $\gamma$ rays and naturally accounts for the source's low radio flux without invoking extra emission zones. We also reproduce the mid-IR emission within the same framework. These results demonstrate that extended jet leptonic models can robustly describe the broadband physics of radio-weak BL Lacs.

Finding the host galaxies of stellar-mass compact binary mergers will open a new window for studying their formation histories and measuring key cosmological parameters, such as the Hubble constant. To date, only one merger, GW170817, has had its host galaxy confidently identified through electromagnetic counterpart observations. The large localization volumes from the LIGO-Virgo-KAGRA (LVK) network, combined with the lack of electromagnetic emission for most events, make host identification challenging. However, as the sensitivity of the gravitational-wave (GW) detector network improves, events are becoming increasingly well localized. Furthermore, galaxy luminosity traces mass or star formation rate, and thus correlates with the probability of hosting a merger. Focusing on the most luminous galaxies within the localization volumes of the best-localized GW events, we estimate the corresponding Hubble constant for each galaxy by combining its redshift with the luminosity distance inferred from LVK observations. For the well-localized LVK events \texttt{S250207bg}, \texttt{GW190814}, and \texttt{S250830bp}, we find only $1$, $1$, and $4$ galaxies, respectively, when restricting the analysis to the most luminous $1\%$ of galaxies above $L_{\rm th} \sim 10^{11} h^{-2} L_{\odot}$ in each event's localization volume and adopting a broad $H_0$ prior. The probability of these galaxies being random, and not associated with the GW events, is $29$-$36\%$ across the three events. We encourage further follow-up observations of these candidate host galaxies. We expect this approach to become increasingly powerful in future LVK observing runs, enabling constraints on merger formation histories and measurements of the Hubble constant.

Models based on variable miscibility among hydrogen, molten silicate, and molten iron, coupled with atmospheric escape, can reproduce the observed occurrence density structure of sub-Neptunes and super-Earths in mass-radius space. The models are also consistent with the radius gap and the observed radius-period relationship exhibited by these planets. The degree of overlap between predicted and observed planetary occurrences suggests that hydrogen-silicate-iron miscibility may serve as a unifying concept for the formation and evolution of these planet classes. The well-defined equilibrium conditions at the boundary between supercritical magma oceans and the overlying hydrogen-rich envelopes are important features of the models. Planets formed with less than ~1 % hydrogen by mass develop discrete, terrestrial-like metallic cores, while those accreting greater hydrogen concentrations are predicted to have fully miscible interiors and no discrete metal cores. Hydrogen-silicate-iron miscibility provides an overarching explanation for the full range of sub-Neptune and super-Earth architectures based on the accreted hydrogen mass fraction and the phase equilibria governing silicate, iron metal, and H$_2$ miscibility.

$\omega$~Centauri, the most massive globular cluster in the Milky Way, exhibits a level of stellar population complexity that has long resisted a unified chemical characterisation. We exploit high-resolution near-infrared spectroscopy from the Milky Way Mapper survey (MWM DR19) to construct one of the largest homogeneously analysed samples of $\omega$~Cen members to date. Applying Ward-linkage hierarchical clustering in a seven-dimensional chemical abundance space, without prior assumptions on population number or boundaries, we identify ten chemically distinct stellar populations. Their nucleosynthetic signatures trace four enrichment channels: iron-peak, $\alpha$-element, CNO-cycle, and high-temperature proton-capture processes. The populations organise into two dominant groups separated by a large light-element spread at a modest iron baseline, consistent with AGB-driven self-enrichment. This dichotomy reflects distinct enrichment pathways: core-collapse supernovae establish the iron baseline, while AGB stars dominate light-element and $s$-process enrichment. A decoupled rise in $s$-process abundances relative to iron-peak elements, together with sub-dominant Type~Ia contributions across all metallicities, indicates evolution on timescales shorter than the characteristic Type~Ia delay time. One intermediate-metallicity population retains a primordial composition, providing evidence for spatially segregated enrichment within the progenitor. The most metal-rich component may trace star formation continuing after accretion into the Milky Way halo. All populations lie in the accreted regime of the $[\mathrm{Al/Fe}]$--$[\mathrm{Mg/Mn}]$ plane, supporting an ex-situ origin. These results reinforce the interpretation of $\omega$~Cen as the remnant nucleus of an accreted dwarf galaxy and provide a framework for future chemo-dynamical studies.

In this thesis, we study the collective modes induced by the chirality of elementary particles in magnetized media and their implications for neutron star phenomenology. We theoretically predict that the chiral magnetic wave can arise in quark matter inside neutron stars, resulting in the emergence of novel types of seismic oscillations and associated gravitational waves in the context of asteroseismology. We also investigate the response to dynamical electromagnetic fields and find that a dynamical screening, distinct from the conventional Landau damping, occurs due to the chiral anomaly.

Gravitational Wave (GW) data bring an exceptional avenue to test the underlying models of coalescing compact objects. In the regime of strong gravity and high curvature, they allow the exploration of minute deviations from the best-fit models, which are difficult to uncover with other observational modalities. These deviations can stem from departures from General Relativity (GR) or unaccounted astrophysical effects. They may not be explainable within the current description of GW strain data, or may simply be difficult to model. However, they are expected to be correlated between detectors and across the population of observed events. The recently developed SCoRe analysis pipeline leverages these properties by focusing on the correlated power between detectors and combining results from multiple events. In this paper, we apply the framework on the Third Gravitational-Wave Transient Catalog to search for source-dependent deviations. In particular, we explore whether there is evidence for a mass-scale in the observed events, which can act like a line of demarcation in their physical properties by exhibiting a deviation that is different above and below this mass-scale. This mass scale dependency naturally arises in gravitational theories described through effective field theories, due to environmental effects or in scenarios involving exotic compact objects, where the GW signature can differ from the standard binary black holes in GR. Using the 30 highest Signal-to-Noise Ratio events in the catalog, we find Bayes factors ranging from 0.16--0.5 (depending on where the threshold mass is set), thus disfavoring the hypothesis of existence of any mass-scale between $\sim 2.5$ M$_\odot$ and $60$ M$_\odot$. We also compute the distribution of excess cross-correlated power across events and find a Bayes factor of $0.07$, which agrees with expected noise statistics.

Statistical inference in observational science typically relies on a fundamental assumption: as sample size increases and uncertainties decrease, the inferred results should converge to the true physical quantities. This assumption underpins the notion that big data lead to more reliable conclusions. In Galactic archaeology, stellar ages inferred from spectroscopic surveys are widely used to reconstruct the formation history of the Milky Way disk. The age metallicity relation (AMR) and its derived formation timescale are often regarded as key physical diagnostics of early disk evolution. This interpretation carries an implicit premise: that observational quality does not introduce systematic bias into age inference. Here we show that this premise may fail. Using a large sample of subgiant stars, we identify a region in the observational quality parameter space (signal-to-noise ratio and parallax precision) where the inferred formation timescale exhibits a systematic offset of 0.5-1 Gyr relative to an independent asteroseismic reference, while the statistical uncertainties remain small, thus producing a stable-but-wrong inference state.

The Paneitz operator is a dimension-4 conformally invariant fourth-order differential operator that has recently attracted attention for possible cancellations of the vacuum energy. We show that, in four dimensions, the Paneitz operator acting on a scalar field falls within the class of extended mimetic gravity theories. Thus, it exhibits the usual instabilities of mimetic gravity. Assuming such instabilities are cured by higher derivative terms, we derive constraints on the Paneitz operator from a modified propagation speed of gravitational waves, after including the Einstein-Hilbert action in the mimetic gravity formulation.

We study the dimension-six SMEFT four-lepton operators in the $\mu$--$\tau$ sector. These operators control both charged-lepton scattering and neutrino self-interactions, the latter being weakly constrained by direct laboratory probes despite their importance for cosmological tensions. We compare three classes of constraints on the Warsaw-basis coefficients $[C_{\ell\ell}]_{2222}$, $[C_{\ell\ell}]_{2233}$, and $[C_{\ell\ell}]_{2332}$. We use perturbative unitarity from $2\!\to\!2$ partial-wave analysis, spin-summing positivity sum rules, and the experimental bounds from NA64$\mu$ and the global fit~\cite{Falkowski:2017pss}. We find that the global fit dominates for $[C_{\ell\ell}]_{2222}$ and $[C_{\ell\ell}]_{2332}$, while NA64$\mu$ provides the leading bound on $[C_{\ell\ell}]_{2233}$, with the unitarity line for this direction entering the range of collider energies near $200~\mathrm{GeV}$. Renormalization-group running between $1~\mathrm{GeV}$ and $30~\mathrm{TeV}$ modifies these coefficients by up to $10\%$. Translating these bounds onto the effective four-neutrino coupling $G_\mathrm{eff}$, we find values many orders of magnitude smaller than the strongly interacting regime motivated by the Hubble tension; this excludes heavy-mediator UV completions of strong $\nu_{\mu}$--$\nu_{\mu}$ and $\nu_\mu$--$\nu_\tau$ self-interactions within the validity of the dimension-six SMEFT and in the absence of tuned cancellations between operators, while leaving the cosmologically motivated light-mediator scenarios unconstrained by this analysis. Finally, we comment on the bounds these coefficients place on a leptophilic $L_\mu - L_\tau$ $Z'$ UV completion. Our SMEFT-based current and projected NA64$\mu$ bounds reproduce the dedicated $Z'$ analyses already available in the literature.

We propose a novel and generalized mechanism, dubbed stochastic axion mixing. In a multi-axion framework, this mixing occurs naturally provided that the masses of all ultra-light axion-like particles (ALPs) are distinct and lighter than the zero-temperature mass of the QCD axion. Crucially, this mechanism is independent of the relative magnitudes of the axion decay constants. In contrast to the conventional maximal mixing scenario -- which strictly relies on specific decay constant hierarchies -- stochastic mixing represents a significantly broader formalism. Notably, maximal mixing emerges as a specific subset of stochastic mixing under restrictive conditions. This new mechanism offers profound implications for axion cosmology.

We introduce a fast, one-dimensional Newtonian {reduced model} to capture uniform rotation in cold white dwarfs, encoding centrifugal support as an effective pressure anisotropy. Using $\Delta_{\rm rot}(r)=\frac{1}{3}\rho(r)\Omega^2 r^2$ derived from the stationary Euler equation with $\langle\sin^2\theta\rangle=2/3$, the model incorporates rotation into hydrostatic balance without a two-dimensional solver. Applying the Chandrasekhar degenerate-electron equation of state, we compute interior structures and global sequences for $ \rho_c \in [10^6, 10^{11}]~{\rm g\,cm^{-3}} $ with rotation proxies $f \le 0.35$, finding monotonic increases in limiting mass and radius, with a percent-level mass gain at $f = 0.35$. We quantify applicability using sub-Keplerian diagnostics evaluated on the rotating configurations, $\max(\Omega/\Omega_K)$ and $\max(\epsilon)$, together with a bulk-interior smallness measure $A_{10^{-2}}\equiv \max_{p_r/p_c\ge 10^{-2}}(\Delta_{\rm rot}/p_r)$. Within the scanned domain these diagnostics remain below unity. The model is therefore best viewed as a reduced Newtonian benchmark for slow-to-moderate rotation, not as a replacement for fully axisymmetric calculations of rotating stars.

Core-collapse supernovae provide natural laboratories for the production of new light particles. In particular, axion-like particles (ALPs) can be constrained via SN1987A cooling arguments. However, significant astrophysical and nuclear uncertainties imply that such bounds may vary strongly depending on modeling choices, even when expensive simulations are employed. In this context, semi-analytic methods offer a simple and fast alternative for deriving new-physics constraints. Building on a previous semi-analytic framework, in which proto-neutron star (PNS) observables are expressed in terms of six global PNS parameters, we include a finite ALP mass in the calculation and derive bounds in the axion-nucleon coupling versus mass plane. The obtained bounds are in good agreement with previous results from numerical simulations, demonstrating the robustness of the method. We also illustrate the sensitivity of the bounds to different PNS parameter calibrations, nuclear effects and cooling exclusion criteria.

Metastable cosmic strings are a common prediction of grand unified theories and act as a source of a gravitational-wave background (GWB) that can explain the 2023 pulsar timing array (PTA) signal. In this paper, we revisit the GWB signal from metastable strings, emphasizing the need to carefully distinguish between two different time scales: (i) t_LB, the time scale of loop breaking because of spontaneous monopole nucleation on closed string loops, and (ii) t_NC, the time scale of network collapse when string segments attached to monopoles begin to enter the Hubble horizon. We discuss under which conditions these two time scales are similar or far apart from each other and illustrate the resulting consequences for the GWB signal. In doing so, we generalize the description of the GWB signal from metastable strings to a three-parameter model in terms of the string tension G\mu and the time scales t_LB and t_NC, which allows us to unify the modeling of standard metastable strings with what is known as quasi-stable strings. In the limit of a large t_LB/t_NC ratio, we, moreover, derive a compact analytical expression for the predicted GWB spectrum in excellent agreement with numerical results in the literature. We thus conclude that our new templates for the GWB spectrum from metastable strings can be readily used in the analysis of future PTA data sets.

We study geodesic motion of a test particle in Schwarzschild spacetime. Bound and scattering geodesics are commonly described using Darwin variables, which provide a convenient parametrization of the radial motion. However, this description breaks down at the separatrix and does not extend straightforwardly to plunging trajectories. We construct an analytic continuation of Darwin variables that yields a real parametrization of bound, scattering, and plunging Schwarzschild geodesics, thereby providing a unified kinematical description of all types of test-mass motion. As a proof of concept, we then apply these variables to a simple non-geodesic evolution in which the energy and angular momentum are driven by a constant external force. This toy model is not intended to represent a physical radiation-reaction model, but rather to illustrate how the extended variables can be used to follow an orbit through a transition to plunge using a single orbital phase variable across the separatrix.