Papers for Friday, May 15 2026

We prove that any background-inert perturbative coupling $ \lambda $ in coincident $ f(Q) $ gravity exhibits a degeneracy with the clustering amplitude $ \sigma_{80} $, when using compressed CMB distance priors. This degeneracy is, in fact, a direct materialization of a more deeper $ A_s-D_0(\lambda) $ degeneracy between the primordial amplitude $ A_s $ and the present day growth factor $ D_0(\lambda) $. We outline a consistency check scheme, applicable to models even outside the $ f(Q) $ class, by computing $ A_s $ needed to reproduce the $ \sigma_{80} $ predicted by the sampler. We perform our analysis with two dataset pipelines, based on the coupled/decoupled $ f\sigma_8(z) $ data. To ensure theoretical diversity, we include $ \Lambda $CDM and the Hybrid model in the $ f(Q) $ framework. Our results illustrate that adding the $ \lambda_0\sqrt{QQ_0} $ correction to the models inflates $ \sigma_{80} $ to unphysical values, while showing moderate evidence in favor of the said models. However, this results in an increase of $ 20\%-30\% $ in $ A_s $ in $ 1.7\sigma-2.2\sigma $ tension with Planck values. We utilize the $ 1\sigma $ $ \ln(A_s) $ constraints from Planck as priors in order to fix the artificial increase in $ \sigma_{80} $ and find that all the constrained parameters return to their baseline values. Each model is penalized by around $ 2 $ units per extra parameter. Interestingly, the $ \Lambda $CDM$ +\lambda_0+\ln(A_s) $ + SDSS DR16 combination shows a weak preference over the vanilla $ \Lambda $CDM model, validated by the values of $ \log\mathcal{Z},\;AIC,\;DIC, $ and BIC.

We investigate the spherical collapse model and the abundance of galaxy clusters in a class of degenerate higher-order scalar--tensor (DHOST) theories in which gravitational waves do not decay into scalar perturbations and which are consistent with current constraints from gravitational-wave observations. We find that deviations from Einstein gravity can become significant at late times when the background universe is close to the scaling regime during the matter-dominated epoch. These deviations suppress the growth of linear matter perturbations on small scales while increasing the extrapolated linear density contrast at collapse, obtained from the spherical collapse model. Using the analytic mass function, we compute the corresponding cluster number counts. The minimum mass threshold in the mass integration for each redshift bin is determined by matching the predicted number counts in the $\Lambda$CDM model with those inferred from the eROSITA survey. We find that the cluster abundance reaches its maximum at low redshift bin, and that the number of clusters in the highest redshift bin is suppressed as the deviation from Einstein gravity becomes larger. The parameters of the theory are chosen such that the deviation from Einstein gravity at present is consistent with the local astrophysical bounds from binary pulsar observations. We find that even under such strict constraints, the upper bound on the deviation leads to lower predicted number counts compared with the $\Lambda$CDM model emulating the eROSITA survey results. However, this may be a consequence of the uncertainties in computing the number counts for the DHOST theories using the spherical collapse model and the analytical mass function.

Theoretical studies predict that dual AGNs are a critical stage of galaxy merger-driven supermassive black hole growth. Systematic searches for dual AGNs typically target late-stage mergers ($\leq10$ kpc nuclear separations) and select AGNs based on optical diagnostics. Yet, simulations predict that obscuration can occur early in the merger sequence, and that a significant fraction of dual AGNs can be found beyond $10$ kpc. Here, we report on a new sample of 157 spatially resolved mid-IR dual AGNs candidates selected based upon their mid-IR $W1-W2$ colors from the Wide-Field Infrared Survey Explorer and optically classified as galaxy merger candidates using imaging from the Dark Energy Camera Legacy Survey. Spectroscopic results are presented for approximately 2/3 of the sample. 76 candidates have been confirmed to reside in galaxy mergers; among these, 13 have been confirmed as bona fide mid-IR dual AGNs, while 63 represent strong dual AGN candidates that require further examination. 46 candidates have been rejected as non-merger contaminants (foreground-background AGNs, separations inconsistent with interacting galaxies, etc.). 35 candidates still await spectroscopic coverage. The confirmed and high confidence dual AGN candidates exhibit separations of 14.5-129 kpc; $>50$% reside at separations $>50$ kpc. Confirmed and high confidence candidates also exhibit a diversity of nuclear optical BPT classes. Seyfert-Seyferts and Seyfert-HIIs dominate the overall BPT pairs sample. 31% of confirmed mid-IR dual AGNs reside in multi-mergers involving three or more galaxies. The diversity in AGN properties and environments identified in this work highlights the importance of multiwavelength selection strategies and analyses in the quest to holistically understand dual AGNs as a population.

We present new fitting functions for the velocity divergence auto- and cross-power spectra, $P_{\theta\theta}(k)$ and $P_{\delta\theta}(k)$, calibrated on gravity-only $N$-body simulations. By applying the Evolution Mapping framework, we revise existing prescriptions to introduce a physically motivated parametrisation in terms of the clustering amplitude $\sigma_{12}$, the RMS density fluctuation smoothed at $12\,\text{Mpc}$. This approach improves robustness and extends the range of applicability beyond that of previous models. Our fits are calibrated using a suite of multi-resolution simulations, with numerical convergence carefully quantified and sampling artefacts mitigated through a conservative patching strategy. This yields converged measurements up to $k\simeq0.56\,\mathrm{Mpc}^{-1}$ and percent-level accuracy for both $P_{\theta\theta}(k)$ and $P_{\delta\theta}(k)$ over a wide range of $\sigma_{12}$. Validation against independent simulations spanning a broad range of cosmological models confirms an accuracy of $1$-$2$ per cent on scales where the measurements are robust, systematically outperforming existing prescriptions. We further assess the impact of deviations from the exact evolution mapping relation induced by differing growth histories. For most cosmologies of practical interest, we find that neglecting these effects introduces only subdominant errors. We show that expressing fitting functions in $h$-dependent units leads to spurious, unphysical dependencies on the Hubble parameter, even for models with identical linear clustering. This provides strong empirical support for parametrising non-linear evolution in terms of $\sigma_{12}$ rather than $\sigma_{8}$. Our fitting functions provide a robust description of velocity power spectra, with direct applications to redshift-space distortion modelling in galaxy redshift surveys.

Understanding the growth of supermassive black holes (SMBHs) requires observational constraints on how their angular momentum (or spin) varies with mass, since the relative importance of coherent accretion, chaotic accretion, and mergers will be reflected in SMBH spin populations. Here we present an updated compilation of reflection-based SMBH spin measurements from the literature and assemble a set of ancillary quantities of interest for each SMBH (including redshift, Eddington ratio, and X-ray luminosity). No obvious apparent correlation between the Eddington-scaled accretion rate and the black hole spin is seen, noting that formal statistical tests are beyond the scope of this review. We discuss the limitations of using this heterogeneous mass--spin sample to test predictions of SMBH growth from semi-analytic models and cosmological simulations, emphasizing the need for a more uniform sample. We then highlight the encouraging prospects enabled by the next-generation NewAthena X-ray flagship observatory. Finally, we summarize how hierarchical Bayesian population inference applied to observed SMBH mass--spin populations will constitute a powerful framework for confirming tentative mass--spin trends in future samples.

The recent success of gravitational-wave (GW) astronomy together with renewed plans for lunar geophysical instrumentation has revived interest in using the Moon as a resonant detector for mid-frequency (mHz-Hz) GWs. In realistic observational scenarios, the GW strain amplitude is expected to be constrained independently by networks of GW detectors, which motivates an inverse, \emph{tomographic} question: to what extent can measurements of the Moon's seismic response to known GWs be used to infer its internal structure? In this work, we develop a first-principles, perturbative framework that maps spherically symmetric perturbations of the elastic and density structure to measurable changes in observables, especially GW-driven modal amplitudes of the Moon. The formalism combines (i) a normal-mode representation of the elastic response, (ii) first-order perturbation theory for eigenvalues and eigenfunctions, and (iii) a linearized observation model that links frequency and amplitude observables to model parameters (bulk and shear moduli, density, and interface locations) and their perturbations. We show that the estimation errors of the Moon's elastic parameters can be reduced by about an order of magnitude with observations of calibrated GWs.

We present the Chandra 25-year Sagittarius A* (Sgr A*) X-ray flare catalog: a systematic analysis of 6.8 Ms of Sgr A* monitoring spanning the Chandra X-ray Observatory's mission lifetime. This is the most complete Chandra Sgr A* X-ray flare catalog to date, consisting of 100 flares with 2$-$10 keV unabsorbed luminosities ranging from $\sim$ 4$-$575 $\times 10^{33}$ erg s$^{-1}$. 18 flares are reported for the first time, including the second brightest Sgr A* flare observed by Chandra. The expanded dataset supports previous indications of a correlation between X-ray flare hardness and luminosity. Spectral modeling corroborates this finding, showing a change in the X-ray spectral index, from $\Gamma \sim 3$ to 2 with increasing flare brightness. Previously-established correlations between flare duration, fluence, and maximum count rate are strengthened via the greater sample size. These results likely reflect variations in the underlying particle distribution that produce weak and strong flares, and the new catalog serves as a rich archive for ongoing observational and numerical investigations into the physical mechanisms responsible for producing Sgr A*'s X-ray flares.

The stellar mass assembly of massive galaxies in the first few billion years of cosmic history remains a central challenge in galaxy formation. Galaxies with $M_\star \gtrsim 10^{10}M_\odot$ observed at $z \gtrsim 4$ must grow rapidly under conditions of intense gas accretion, feedback, and dust production. Observationally, their star-formation histories (SFHs) have been poorly constrained due to degeneracies inherent to broadband photometry. The advent of JWST enables direct spectroscopic access to detailed continuum shapes and rest-frame optical diagnostics at high redshift, providing a critical opportunity to reconstruct formation timescales of massive early galaxies. Here, we investigate massive galaxies using joint spectro-photometric SED fitting of JWST/NIRSpec prism spectroscopy from the CANDELS-Area Prism Epoch of Reionization Survey (CAPERS). Our sample comprises 148 galaxies selected photometrically with log $(M_\star/M_\odot) > 9.5$ at $z > 3.5$. We find that the most massive galaxies (log $(M_\star/M_\odot) > 10.5$) preferentially exhibit shallow, gray dust attenuation curves, consistent with higher dust optical depths and large grain sizes. We also find significant diversity in the time at which galaxies form 25% of their stellar mass. While formation timescales converge toward later cosmic times, galaxies with lower sSFR ($\lesssim -9$) at the observation epoch formed significantly earlier than systems with higher sSFRs. Across the full mass range, inferred assembly times are systematically earlier than model predictions, suggesting more rapid early growth than currently captured theoretically. These results underscore the importance of spectroscopic constraints and flexible SFH and dust models for reconstructing high-redshift massive galaxy formation histories.

Recent JWST observations indicate that star formation at $z\!\sim\!4-6$ is more stochastic than previously assumed, with rms log-SFR scatter $\sim\!0.6$ dex at $M_h\!\sim\!10^{11}M_{\odot}$, growing toward smaller halos and time-correlated on $\sim\!25$ Myr. This is significantly higher than the typical $\sim\!0.3$ dex phenomenological lognormal scatter assumed in standard line-intensity mapping (LIM) forecasts. We propagate the JWST-era burstiness through to the LIM shot-noise power spectrum and show that the result is a simple multiplicative correction: the deterministic shot noise multiplied by a line-dependent boost factor $B_\lambda$ derived in closed form by convolving the SFR correlation function with the stellar-population-synthesis kernel of each line. At $z\!\sim\!6$, we find $B_{{\rm H}\alpha}\!\simeq\!7$ and $B\!\sim\!2.5$-$3.5$ for longer-window tracers ([CII], CO, UV) - factors of $\sim\!2$-$5$ above the standard prescription, and growing further toward higher redshift. The enhancement transforms the LIM landscape: it improves auto-spectrum detectability and suppresses lower-redshift interloper contamination, but degrades cosmological applications such as BAO that rely on a clean clustering measurement. Crucially, it also opens a new use of LIM as a diagnostic of high-redshift star-formation physics beyond the regime of individually resolved galaxies: redshift tomography of a single line constrains the amplitude and mass dependence of the burstiness, while cross-line shot-noise correlations probe its time coherence.

We aim to develop a state-of-the-art tool to infer detailed star formation histories (SFHs) and age-metallicity relations from realistic observational data, while mitigating classical degeneracies and substantially reducing computational cost. In particular, we seek to exploit the complementarity of spectroscopic and photometric data to improve constraints on the spatially resolved SFH and metallicity evolution of nearby galaxies in the PHANGS collaboration. We construct and train a convolutional neural network (CNN) that combines convolutional layers, attention mechanisms, and a shared latent space to jointly predict SFHs and metallicities in 16 age bins. The network simultaneously processes integral-field spectroscopic data from PHANGS-MUSE and five-band photometric fluxes from PHANGS-HST. Training is performed on a dataset of 165\,000 synthetic spectra and photometric measurements spanning a broad range of SFH shapes, metallicity evolution, dust attenuation, and signal-to-noise ratios representative of the observations. The CNN accurately recovers SFHs and age-metallicity relations over a wide range of evolutionary scenarios. The inferred luminosity- and mass-weighted mean ages and metallicities show negligible bias, with dispersions of $\sim0.12$ dex in age and $\sim0.03$ dex in metallicity. When applied to real PHANGS-MUSE and PHANGS-HST data for NGC\,3627, the network produces smooth, spatially coherent maps of stellar age and metallicity that recover physically meaningful structures, including younger populations tracing the spiral arms and star-forming regions. The CNN is approximately $5\times10^{3}$--$2\times10^{4}$ times faster than traditional full spectral fitting codes, providing a powerful and efficient alternative for the analysis of large spectro-photometric surveys.

Observations with the JWST in lensed fields have revealed that galaxies at cosmic dawn may concentrate their star formation in highly dense, compact, star clusters. The high columns and low metallicities encountered in their birth environments suggest that Lyman-alpha (Ly$\alpha$) radiation pressure may be crucial to their formation and evolution. In this study, we address this question by post-processing snapshots from radiation hydrodynamic simulations of dense star cluster-forming clouds ($\Sigma_*\gtrsim10^3{M_\odot{pc}^{-2}}$) with a range of dust abundances ($Z_d=0-0.1Z_{d,\odot}$) using the COLT Monte Carlo code. We infer that Ly$\alpha$ is likely to have mild (~10%) effects on the gas-to-star conversion efficiencies ($\epsilon_*\gtrsim60$%) for $Z_d\gtrsim0.01Z_{d,\odot}$, and even in dust-free environments, $\epsilon_*\gtrsim25$% - much higher than the <10% values typical of star-forming regions in the local Universe. This is because the densest filaments dominating stellar mass assembly ($n\gtrsim10^4{cm}^{-3}$) remain sub-Eddington ($f_{Edd}<1$). On the other hand, the bulk of the gas volume ($n\lesssim10^3{cm}^{-3}$) has $f_{Edd}>1$, with noticeable fractions having $f_{Edd}\gtrsim10$, implying that Ly$\alpha$ can launch dynamically significant winds from these systems rapidly ($\lesssim$4Myr), with possible implications for ionizing photon escape and galactic outflows. The Ly$\alpha$ force multiplier $M_F$ is highly sensitive to $Z_d$, with $M_F\lesssim3$ ($\lesssim 500$) for $0.1Z_{d,\odot}$ (dust-free) environments respectively. Nevertheless, Ly$\alpha$ dominates over UV and IR radiation pressure at all values of $Z_d\lesssim0.1Z_{d,\odot}$, by factors of ~3-500. Our results suggest that Ly$\alpha$ radiation pressure reinforces the emerging picture of locally efficient, bursty star formation accompanied by rapid outflows in galaxies at cosmic dawn.

With the next generation of both electromagnetic and gravitational wave observatories beginning to come online, rapid analysis methods for kilonova data are becoming increasingly important in astronomy. Traditional Bayesian parameter estimation using Markov chain Monte Carlo (MCMC) is time-consuming and relies on explicit likelihood approximations that can break down when modeling uncertainties are significant. We develop a simulation-based inference (SBI) framework for kilonova parameter estimation using density-estimation likelihood-free inference. The framework uses a Gaussian process emulator trained on $\sim1300$ radiative transfer simulations generated with the POSSIS code. We demonstrate that SBI provides a rapid alternative to MCMC for inference with emulators or approximate likelihoods that is robust to emulator uncertainty and likelihood misspecification. On simulated data, the SBI method accurately recovers injected parameters and produces posterior predictive light curves consistent with the data, but the MCMC posterior recovery suffers from systematic bias caused by likelihood misspecification. When analyzing AT2017gfo, the SBI and MCMC methods yield similar light-curve predictions but different posterior distributions, with a subset of the MCMC posteriors piling up at prior boundaries. The likelihood in the MCMC fails to capture the non-Gaussian, correlated structure of the emulator uncertainty, but SBI learns the posterior directly from forward simulations that include the full predictive distribution. Once trained, the SBI framework generates $\sim2\times10^4$ posterior samples in seconds.

Sgr A*, the 4 x 10^6 solar-mass supermassive black hole at the Galactic Center, exhibits frequent flaring with X-ray luminosities of L_X ~ 10^35--10^36 erg s^-1, while its ultraviolet (UV) emission remains unconstrained due to extreme extinction (A_V ~ 30 mag). We use JWST/MIRI time-resolved spectroscopy of the central Galactic Center's 0.3 arcsec region to search for mid-infrared emission-line variability driven by Sgr A* flares, comparing the results to CLOUDY photoionization models spanning flare luminosities of L_UV = 10^32--10^39 erg s^-1 in a dense medium. We detect no statistically significant variability in any mid-infrared line, including [Fe II] 5.34 micron, [Ne II] 12.813 micron, [Fe II] 17.936 micron, and [S III] 18.713 micron. Despite expectations of a flare-driven response, we show that the lack of variability is consistent with the physical conditions in the spatially extended line-emitting gas, where light-crossing timescales of ~0.1--10 days and recombination and cooling timescales much longer than the flare timescale suppress any observable response to individual flares. We further find that the predicted emission is continuum dominated and that even the brightest lines are intrinsically weak and broadened by velocities of order 10^3 km s^-1, reducing their contrast against the continuum and limiting their detectability. Extending the analysis to higher-ionization mid-infrared and near-infrared lines does not improve sensitivity. These results demonstrate that infrared emission lines trace a steady-state radiation field rather than individual flaring events, and therefore infrared line diagnostics cannot be used to constrain the instantaneous UV flux of Sgr A*.

Sgr A* is the nearest quiescent supermassive black hole, and its proximity offers a unique opportunity to study its surrounding fuel supply. We leverage extensive spatial and spectroscopic information provided by the \jwst/MIRI MRS instrument to disentangle mid-infrared ionized gas structures in the central 0.1 parsec of the Galaxy. The Galactic Minispiral's Bar and Northern Arm are revealed by their distinct morphological and kinematic signatures. Several compact ($<1$\arcsec) gas structures including X7 also appear within $\sim 0.05$ parsec of Sgr A* in the plane of the sky, moving with blue-shifted radial velocities $\gtrsim 600$ km/s. Fine structure line measurements spanning ionization energies $\sim 7 - 55$ eV are used to constrain the incident radiation field, metal abundances (neon, argon, sulfur, nickel, and iron), and dust depletion/destruction for each identified gas structure. Overall, the Minispiral gas metallicity is $\sim 1-2.5~Z_\sun$, with a Wolf-Rayet star-driven ionizing radiation field, and significant nickel and iron dust destruction. Increased flux at energies $\gtrsim 41$ eV suggests that the compact gas structures experience an additional harder ionizing radiation source, which is most likely driven by localized fast radiative shocks from stellar winds, a hypothetical Sgr A* outflow, and/or interactions with the ambient medium.

During the hard and quiescent spectral states in X-ray binaries, a non-linear correlation is observed between radio and X-ray luminosities, providing a valuable tool to probe the connection between accretion and jet production. This relation was originally thought to define a single 'standard' correlation spanning several orders of magnitude in X-ray luminosity, and was extended to active galactic nuclei by including a mass term. However, subsequent studies revealed a more complex picture, with some sources deviating from the standard correlation and instead populating distinct tracks. To date, all large studies of the radio:X-ray plane have combined data from multiple telescopes, introducing uncertainties due to differing instrument systematics and flux conversions between observing frequencies, thereby complicating comparisons and limiting constraints. ThunderKAT was a five-year programme on the MeerKAT radio telescope that monitored X-ray binaries in outburst, and ran alongside SwiftKAT which provided quasi-simultaneous Swift/XRT X-ray coverage. We present the full set of light curves from these programmes, comprising 948 radio and 1029 X-ray data points. An important finding is the frequent detection of unresolved radio emission during the soft state, likely dominated by previously launched jet ejecta. Using these data, we construct the largest, observationally homogeneous X-ray binary radio:X-ray plane to date. We relate these results to the physical mechanisms proposed to drive inter-source diversity, and outline directions for future observational and theoretical work. This paper is accompanied by a public data release of the ThunderKAT and SwiftKAT measurements and a compiled radio:X-ray plane, available through an interactive website.

Procyon is a great system to probe stellar evolution of non-interacting binaries. We present an extensive grid of MESA (Modules for Experiments in Stellar Astrophysics) evolutionary tracks to constrain the evolution of Procyon A and B. We systematically vary the initial parameters of our grid anchored by precise dynamical masses and spectroscopic determinations of effective temperature (T$_{eff}$) and luminosity ($L$) to match the stars' positions on the H-R Diagram. Our goal is two-fold: (i) to quantify how the inferred system age and the progenitor mass of Procyon B depend on metallicity ($Z$), mixing length ($\alpha$), and core overshoot ($\beta$), and (ii) to determine the best fitting model of Procyon B within a model-based initial-to-final mass relationship (IMFR) for hydrogen-deficient white dwarfs. Our best-fit models reproduce the observed properties for both components, yielding $M_\mathrm{A}=1.487 \pm 0.095$ M$_{\odot}$, $M_\mathrm{B}=0.592\pm 0.082$ M$_{\odot}$, a system age of $2.23 \pm 0.90$ Gyr, and a white dwarf cooling age of $1.20\pm0.49$ Gyr for Procyon B, consistent with independent determinations. Our results point to higher core overshoot than the standard adopted range, with the best fits ranging from $\beta=0.5-1.0$. From our model grid, we map Procyon B to the initial-to-final mass relationship for H-deficient white dwarfs in the $1.9\!-\!2.6$ M$_{\odot}$ progenitor range. Additionally, we implement the accretion of heavy metals onto the surface of the WD and fit our isotopic abundances to spectroscopic observations. We outline the physics used in our analysis.

The Galactic Center hosts the densest known stellar environment in the Milky Way, dominated by the massive black hole Sgr A* and the surrounding nuclear star cluster. Theory predicts that this region should also contain a large population of stellar compact objects (SCOs) - black holes, neutron stars, and white dwarfs - forming a "dark cluster" whose distribution and properties remain observationally unconstrained. These unseen stellar remnants are central to questions of mass segregation, cluster dynamics, and the expected rate of extreme mass ratio inspirals (EMRIs) detectable by future gravitational-wave observatories including LISA. Current evidence for SCOs in the Galactic Center is indirect, relying on dynamical mass measurements, X-ray surveys, and a small number of transient sources. Direct detections remain elusive due to crowding, extinction, and the sensitivity limits of existing instruments. We explore how upcoming facilities, in particular the Extremely Large Telescope (ELT) with its first-light imager MICADO, can fundamentally transform this field. MICADO's combination of deep photometry, high spatial resolution, and precise astrometry will enable systematic searches for SCO-star binaries via photometric variability and orbital astrometric signatures, as well as direct detection of isolated accreting black holes interacting with the gas-rich Galactic Center environment. We outline the observational pathways, technical challenges, and expected sensitivities, showing that ELT/MICADO observations can provide the first quantitative constraints on the dark cluster population. Establishing these constraints will be pivotal for understanding the dynamical evolution of the Galactic Center, the role of compact remnants in nuclear star clusters, and the astrophysical context of gravitational-wave sources in galactic nuclei.

Binary stars undergoing mass transfer provide unique laboratories for testing stellar evolution. Here, we present a comprehensive photometric and spectroscopic analysis of the semi-detached system ET Cru. Using spectral disentangling, we independently determined the effective temperatures and chemical abundances of both components with high precision, including nine elements (eleven species). We find masses of $13.41\,M_\odot$ and $6.00\,M_\odot$ for the primary and secondary, respectively, with uncertainties of only $\sim$1.3%. The radii are $5.58\,R_\odot$ and $5.68\,R_\odot$, measured to within 0.4% and 0.5%. Surface gravities are constrained to better than 1%, while effective temperatures are determined to within 3-5%. The secondary exhibits extreme chemical anomalies, with severe carbon depletion and nitrogen enrichment far exceeding those reported in classical Algol systems. Multi-wavelength spectral energy distribution modelling yields a distance of $\sim$2.5 kpc, inconsistent with the $Gaia$ DR3 parallax, suggesting systematic astrometric uncertainties in the parallax distance. Together, these results establish ET Cru as a benchmark Algol-type binary, revealing direct spectroscopic evidence of deep CNO-cycle exposure in the donor and confirming the primary star as a rejuvenated gainer. ET Cru thus provides a chemically and dynamically illustrative case for understanding advanced binary interactions and the late evolutionary stages of massive-star evolution.

Ultralight primordial black holes (PBHs) provide a compelling window into early-Universe cosmology. Following our earlier work, we explore a mechanism for the formation of ultralight PBHs sourced by primordial inflationary magnetic fields, without invoking an ultra-slow-roll phase of inflation. We propose a magnetogenesis model in which large curvature perturbations are induced at small scales, leading to the efficient production of ultralight PBHs across a broad mass spectrum. We analyze the phenomenological implications of these ultralight PBHs for early-Universe cosmology, particularly during reheating. We compute the resulting stochastic gravitational wave (GW) background generated by both the electromagnetic spectrum and evaporating PBHs, which exhibits distinctive features tied to the underlying magnetogenesis model parameters. Our results demonstrate that inflationary magnetic fields can serve as a viable and testable origin for ultralight PBHs, opening new avenues for probing the interplay between inflation, magnetogenesis, PBHs, and primordial gravitational waves.

Planets around pulsars were the first discovered exoplanets, found thanks to the extremely precise pulsar timing. Here we suggest that they could also be found through the radio emission generated by the pulsar-planet magnetospheric interaction. We present the results of special relativistic numerical simulations of planets in a pulsar wind of velocity $v=0.985~c$, corresponding to a Lorentz factor $\gamma=5.795$. Planets, modeled as a perfectly conducting solid surface in an external magnetic field originating from the pulsar, produce radio emission in the extended magnetic structure on the planet's nightside. We find that the planet around a known pulsar, PSR J0636+5129 b, could be detected via its radio emission. We outline the observational prospects for such objects.

We present non-local thermodynamic equilibrium radiative-transfer modeling of SN2023ixf during the nebular phase out to 1000d, using the same ejecta that matched its photospheric evolution, namely a partially stripped red-supergiant star of initially 15Msun whose terminal explosion yielded ejecta with 7-8Msun, kinetic energy of 1.2e51erg, and 56Ni mass of 0.05Msun, augmented with a cold dense shell (CDS) of 0.2Msun at 8000km/s. Interaction with circumstellar material persists at all epochs, powering the ultraviolet (UV) flux at all times, but dominating the optical only after ~600d. Matching the V-band light curve requires invoking both enhanced gamma-ray escape and dust formation after ~200d, first in the CDS and eventually in the inner ejecta as well. Depending on where they form relative to the dust, emission lines are uniformly attenuated or skewed with a blue-red asymmetry. Our models suggest a rising dust mass (chosen as an C-rich and Si-rich mixture) in the CDS and inner ejecta, possibly reaching 1e-4Msun at 700d, while an external cold dust component is required to match the mid-infrared emission. The UV radiation, largely unaffected by dust, is influenced by the emission and absorption from Fe lines, together with strong, blueshifted emission from Lyalpha and MgII2800, both present at >~200d and with a strengthening fractional flux thereafter. Optical-depth effects play a critical role for the UV flux, and most notably on Lyalpha whose strength depends strongly on the CDS structure (mass and extent) and the treatment of power injection. The CDS is continuously slowing down from 8000km/s at 112d to ~6500km/s at 998d, suggesting a growth in mass of several 0.1Msun. SN2023ixf shares many similarities with SN1993J at 1-3yr, but it is eventually fainter due to dust extinction and cooler (i.e., weak [NII] and no [OIII] lines) likely as a result of greater CDS and ejecta masses.

In this work, we perform Markov Chain Monte Carlo (MCMC) analyses using the Planck 2018 cosmic microwave background (CMB) datasets, including temperature, polarization, and lensing, in order to compare matter bounce models with observational data. The particular model we considered contains a scalar field with an exponential potential, which behaves as dust in the asymptotic past of the contracting phase, it realizes a quantum bounce, and then behaves as a transient dark energy field at large scales in the expanding phase. The parameter $\lambda$ appearing in the exponential potential is directly related to the model's scalar spectral index, $n_s$, which is set free in the MCMC analyses, as well as the deepness of the bounce, which controls the amplitude of the power spectrum. We provide constraints on the cosmological parameters and compare the model's performance against the standard inflationary $\Lambda$CDM scenario. Our results indicate that Planck data alone cannot favor one model with respect to the other, showing that the model we investigate can be a viable alternative to inflation.

Available JWST observations TRAPPIST-1 system have suggested that several of the planets are likely airless, or possess a very tenuous atmosphere. However, the high atmospheric escape rates expected for these planets suggest that any tenuous atmosphere must be replenished by constant outgassing, and past studies on modeling potential atmospheres for the planets have not widely considered surface pressures <1 bar. Here, we show that tenuous atmospheres on the TRAPPIST-1 planets are likely possible, supported by constant plausible rates of water and/or CO$_{2}$ outgassing against assumed high escape rates (up to ~10$^{30}$ s$^{-1}$). We use a coupled photochemical-climate model and sample from a broad phase space of outgassing, surface deposition, and top-of-atmosphere escape rates to test hundreds of atmospheres per planet. Critically, our model also allows surface pressure to vary based on the balance of sources and sinks. We find that 6 different compositional archetypes are generated via H$_{2}$O and/or CO$_{2}$ outgassing across our phase space, and atmospheres commonly fall between 10$^{-4}$ -- 1 bar. We find that potentially habitable surface environments are possible for TRAPPIST-1d and e at pressures between 0.05 -- 2 bar and 0.5 -- 1 bar, respectively. Where possible, we compare our models to JWST observational data for TRAPPIST-1b, c, d, and e; all atmospheres found in this study for these planets match available transmission data to <3$\sigma$. However, emission data are consistent with atmospheric outcomes constrained to thin O$_{2}$-dominated compositions for TRAPPIST-1b ($\lesssim$0.01 bars) and c ($\lesssim$0.2 bars), which may or may not contain trace SO$_{2}$.

Magnetic switchbacks, large-amplitude, localized Alfvenic like rotations of the solar wind magnetic field, have been the subject of intensive investigation, with approximately 200 refereed papers published in the last decade. Yet, fundamental controversies persist regarding whether switchbacks can be described with Ideal MHD (magnetohydrodynamic) physics or Hall-MHD physics and whether their origin is at the solar surface or in the solar wind. To settle these controversies, we present Parker Solar Probe electric field measurements between 13 and 40 solar radii, which show that switchbacks have non-zero electric fields in the plasma frame, a finding that definitively settles the physics controversy by proving that switchbacks are Hall-MHD, not Ideal MHD, structures. Along with these electric fields, there are enhanced Poynting vectors having three components with similar magnitudes that exist only inside the switchbacks. These facts contradict the view of switchbacks as simple outward-propagating pulses. Together, they resolve one controversy by showing that switchbacks in the young solar wind are a non-MHD process. They contribute to the second (source) controversy by identifying switchbacks as sites of active, in-situ, evolution. These findings provide a new framework for understanding energy transport and dissipation in astrophysical plasmas.

Evolutionary tracks for $2-8M_\odot$ stars, with metallicities of $Z=0.014$, $0.006$, and $0.002$, including rotation, are computed with Modules for Experiments in Stellar Astrophysics (MESA). We study how rotation impacts the evolutionary properties of classical Cepheids. We examine whether rotation can offer a plausible explanation for the mass discrepancy problem when it is included in the evolutionary code using the fully diffusive approximation for rotationally induced mixing processes. We find that rotation barely influences the appearance and luminosity levels of the blue loops. While luminosity increases with increasing initial rotation rate, the increase does not exceed 0.04 dex, a fraction of the increase resulting from including the main sequence (MS) core overshooting of $0.2H_p$. As a consequence, rotation alone cannot resolve the mass discrepancy problem without simultaneously requiring significant MS core overshooting. Similar to the mass-luminosity relation, the period-radius and period-luminosity relations are barely affected by rotation, while the period-age relation predicts Cepheid ages to be only a few per cent longer compared with models without rotation. The predicted surface rotational velocities are too large compared with observations. These results are in contrast with those obtained with the Geneva code, which implements rotational mixing using the advective-diffusive scheme. In that approach, the luminosity levels of the loops are significantly higher, their luminosity extent increases, and the predicted rotation velocities are lower, compared with MESA models. The differences between the two approaches arise from significantly more efficient rotation-induced mixing during the MS evolution in models computed with the advective-diffusive scheme.

Scattered light is one of the most common sources of noise in the LIGO gravitational wave detectors. Light scattering is a highly non-linear process through which motion at low frequencies gets up-converted and creates noise in a higher frequency band in the detector data. From the beginning of the fourth observation run, many glitches appeared in the data of LIGO Livingston detector in the frequency range 10-40 Hz, and the morphology of these glitches suggested that they were produced by scattered light. From our analysis, we identified two different populations of scattered light glitches, one group having higher SNR than the other. The glitches of the high- SNR group were solely modulated by microseismic ground motion (ground motion in 0.1-1.0 Hz) and in this paper, we present models of possible coupling mechanisms for these glitches. We also present results of a statistical correlation analysis based on our models, which indicates that the microseismic ground motion at the corner station along the X direction is the one most correlated with the noise which create these high SNR glitches. After installing baffles very close to the test mass mirrors, we have noticed a significant reduction in the rate and SNR of these glitches. The low-SNR glitches were primarily modulated by high frequency (10-30 Hz) vertical ground motion at the corner station, and this motion was coupling through a specific vacuum chamber at the corner station. After installing an additional seismic isolation platform in that vacuum chamber, these glitches have disappeared.

The Chandra Interactive Analysis of Observations (CIAO) software, developed by the Chandra X-ray Center, has been the data analysis package for the Chandra X-ray Observatory since its launch in 1999. Over nearly three decades, CIAO has grown from a small software suite into a widely used system for X-ray data analysis and beyond. CIAO provides tools for calibration, spectral, imaging, and timing analysis, together with high-level scripts and the \sherpa\ modeling and fitting application. Its modular design and unified data model allow users to build flexible analysis workflows while maintaining consistency with the Chandra data processing pipeline. Visualization capabilities are provided through integration with SAOImageDS9 and Python-based tools, and simulation components such as ChaRT and MARX extend the analysis environment to include detailed modeling of instrumental effects. In this paper we describe CIAO's design, evolution, and capabilities after 25 years of Chandra operations. We also describe its core architecture, scripting environment, modeling, visualization tools, simulation components, and testing infrastructure, as well as the documentation and user support system that have contributed to its widespread use. CIAO's continued development and broad adoption highlight its important role in X-ray astronomy and its usefulness in multiwavelength astrophysical research.

We examine the blue straggler star (BSS) populations of six old ($\geq$4 Gyr) open clusters: M67, NGC 188, NGC 6791, Berkeley 32, Berkeley 39, and Trumpler 19. We find that 50% of BSSs have color-magnitude diagram (CMD) locations corresponding to single stars in the final third of their main-sequence lifetimes. This build-up of BSSs near the terminal-age main sequence (TAMS) is primarily, but not solely, driven by more massive BSSs. Eleven of the BSSs have white dwarf companions with measured cooling ages; their evolution age distributions indicate that more massive BSSs typically form far from the zero-age main sequence, whereas lower mass BSSs can form at every evolutionary age. We show that inferred core helium amounts (above primordial) of late-evolution-age BSSs correspond to the core helium fused by cluster main-sequence stars near the turnoffs. We also find that the masses of asymptotic giant branch (AGB) mass-transfer BSSs require evolved main-sequence accretors and conservative mass transfer. These findings indicate that helium enrichment of progenitor accretors leads to the prevalence of BSSs near the TAMS. We further classify the evolutionary stages of the progenitor donors in M67 and NGC 188 and find mass transfer during the AGB accounts for at least half of the BSSs. We trace how the main-sequence binary population of NGC 188 evolves, and find that only 30-40% of interacting binaries create BSSs and that progenitor orbits must change to match current BSS periods.

The hitherto neglected close binary V486 Car is studied with the aid of newly applied satellite photometry (HIPPARCOS and TESS), high dispersion spectrometry (HERCULES) and ground-based B and V photometry. While the sinusoidal light variations are suggestive of a near-contact system, the stars have only shallow eclipse, so highly confident parametrization becomes challenging. We find: $M_1 = 2.1 \pm 0.1$, $M_2 = 0.4 \pm 0.1$; $R_1 = 3.20 \pm 0.02$, $R_2 = 1.48 \pm 0.01$; (${\odot}$); $T_{e1} = 10000 \pm 500$, $T_{e2} = 6200 \pm 200$ (K); distance = 162 $\pm$ 12 (pc). New times of minima for V486 Car have been examined, including recent observations from TESS. The role of the relatively significant O'Connell effect is examined. As well as the conspicuous asymmetry from the main effect of about 0.036 mag (V), a jitter, with amplitude of about 0.005 V mag and quasi-period of order $\sim$ 10 d is noticed. There is a tendency for such photometric excursions at one maximum to precede those at the other. As well, the O -- C data indicate the presence of a low mass star $\sim$0.3 M$_{\odot}$ in an orbit separated by a few AU from the close binary. More accurate and plentiful spectroscopic data would be requisite for further investigations. A brief discussion reviews possible approaches to understanding the system in the context of near-contact binary scenarios.

JWST has unveiled an abundant population of compact broad-line emitters largely at $z\gtrsim4$, the Little Red Dots (LRDs), which might represent a previously unprobed supermassive black hole evolution channel predominant at high redshift. However, the LRDs have remained mostly elusive at lower redshift ($z\lesssim2$) where detailed studies are possible from ground-based observatories. We searched for low-redshift LRDs in the Dark Energy Spectroscopic Instrument (DESI) survey. Our search is primarily based on emission line properties, as opposed to earlier approaches that searched for compact sources with specific photometric spectral energy distributions. We report the discovery of eight LRDs at $z=0.2-0.45$, which show spectral features akin to the high-redshift LRDs in the rest-frame optical. The sources are characterized by broad Balmer lines, steep Balmer decrements, compact morphologies, Balmer absorption features and/or strong He I emission, but weak or absent He II, [Ne V] or other high excitation lines typical of Type I AGN. For 7 out of 8 sources, we retrieve dense-cadence light curves from time-domain surveys and for most sources we find weak to no intrinsic variability ($0.0-0.1$ mag) over $4-17$ years in the rest-frame. We also highlight the identification of a quasar with similar Balmer line profiles as LRDs, but shows differences in Balmer decrement, significant variability, and high-ionisation lines. Given the effective volume $4.9{\rm Gpc^3}$ covered by DESI DR1 at $z<0.45$, our sample corresponds to a number density of $1.6\times10^{-9}$Mpc$^{-3}$, indicating a number density $\sim$10,000 times lower than in the first billion years of cosmic time. We find a dearth of luminous and red LRDs at $z<1$ compared to higher-redshift, which could suggest lower gas feeding rates of LRD activity due to higher metallicities at later cosmic epochs.

Catastrophic collisions between proto-satellites have been proposed as a possible origin of Saturn's rings. This argument relies on the concept of the equivalent circular orbit. Here, we re-examine the post-impact dynamical evolution of collision debris using analytical arguments and $N$-body simulations with fragmentation. We focus on the long-term evolution of debris distributed in a broad V-shaped region in the $a$--$e$ plane, with two arms for particles sharing a common collision radius. Because particles on the two arms possess significantly different angular momenta, inter-arm collisions dominate the evolution and drive behavior fundamentally different from the simple circularization assumed in the equivalent circular orbit approach. As a result, the classical equivalent circular orbit concept cannot predict the long-term fate of collision debris. Both our analytical framework and $N$-body simulations show that, although some debris initially passes within the Roche limit on eccentric orbits, successive collisional evolution drives the particles approximately along the original V-shaped constraint curves toward the apex of the V-shape, i.e., the original collision radius. Instead of spreading inward to form a ring, the debris converges and reaccretes near the original collision location. We therefore conclude that catastrophic proto-satellite collisions do not produce massive Saturnian rings. Rather, the debris evolves toward reaccretion into a new generation of satellite-sized bodies near the impact radius. These results fundamentally revise the dynamical interpretation of collision-generated debris and establish a more general framework applicable beyond the Saturnian system, including other planetary ring systems and debris produced during planet formation.

A recent analysis of photospheric radius expansion X-ray bursts from the low-mass X-ray binary 4U 1746-37 reported unusually small mass and radius estimates for the neutron star, suggesting it could be a quark star or quark-cluster star. Here, we propose an alternative interpretation: the star's mass and radius could be underestimated from significant blocking of the X-ray flux. Significant blocking factors ($\mathcal{B} \gtrsim 0.8$, reducing the observed flux to $\sim17\%$ of the intrinsic emission) permit neutron star parameters consistent with the canonical values: $M = 1.59 \pm 0.69 M_{\odot}$, $R = 13.0 \pm 5.45\,\mathrm{km}$, or $M = 2.12 \pm 1.08 M_{\odot}$, $R = 9.80 \pm 4.13\,\mathrm{km}$. The blocking factor, which varies with the photospheric radius, provides a natural explanation for the anomalously large peak-to-touchdown flux ratio ($\sim2.0$) and highlights the importance of accounting for geometric system configuration in neutron star mass--radius estimates.

Atom interferometers represent a promising approach for gravitational wave detection in the decihertz frequency band, complementary to existing light-based detectors. The South Pole offers unique advantages for such experiments: exceptionally low seismic noise, established infrastructure for large scientific projects, and a location that strengthens gravitational wave source localization through global triangulation. Here we discuss the scientific case and practical considerations for deploying a long-baseline atom interferometer at the South Pole, which has the potential to expand the global network of gravitational wave detectors while enabling precision tests of fundamental physics.

We present a multi-wavelength physical characterization of 14 luminous Ly$\alpha$ emitters (LAEs) at $z\approx6$, integrating deep ground-based Magellan/M2FS spectroscopy with heterogeneous JWST/NIRCam broad- and medium-band imaging. Identified via strong Ly$\alpha$ lines with extreme Ly$\alpha$ luminosities of ${>}10^{42.6}$ erg s$^{-1}$, the sample exhibits very large rest-frame equivalent widths (${\gtrsim}100$ Å) and steeply blue UV continua ($\beta_{\rm median}\simeq-2.2$, $-18.2>M_{\rm 1500}>-20.2$ mag). Crucially, the integration of NIRCam medium-band photometry (F410M) breaks the degeneracy between strong rest-optical nebular emission and Balmer breaks, resolving prior mass overestimations. The tightly constrained spectral energy distribution modeling demonstrates that these luminous LAEs tend to be unequivocally low-mass, ultra-young dwarf starbursts; half the sample is characterized by stellar masses of $M_* < 10^8 M_{\odot}$, ages $\lesssim10$ Myr, and negligible dust attenuation. We also map the production efficiency of ionizing photons and Ly$\alpha$ escape fractions ($f_{\rm esc}^{\rm Ly\alpha}$). The $f_{\rm esc}^{\rm Ly\alpha}$ values are exceptionally high, with a median of ${\gtrsim}40$%, increasing for the bluer UV continua. Finally, analyzing spatial offsets between the Ly$\alpha$ centroid and the stellar counterpart, we demonstrate empirically that internal dust content, rather than neutral hydrogen gas, dominate the suppression of Ly$\alpha$ radiative transfer. Our study reveals that strong Ly$\alpha$ emission of the luminous LAEs are generally attributed to both the vigorous starburst activities and the high $f_{\rm esc}^{\rm Ly\alpha}$. Resembling Lyman continuum leakers, these extreme dwarf systems function as highly efficient ionizing engines at the conclusion of the Epoch of Reionization.

We present a $Z_1^2$-based framework for estimating the spin frequency and frequency derivative of high-energy pulsars from Poisson-limited photon event lists. The key point is that the width of a coherent detection peak is not, by itself, the statistical uncertainty on the recovered rotational parameters. We develop and compare three computationally efficient estimators: segmented frequency regression, a coherent derivative scan, and a localized two-dimensional coherent fit. For sinusoidal signals, we derive the local form of the Z-squared response as a function of frequency and frequency derivative, and show that expressing the frequency at the midpoint of the observation removes the leading-order covariance between the two parameters. This gives simple uncertainty estimates in terms of the fitted peak amplitude and local widths, without requiring an exhaustive Monte Carlo simulation for each observation. We test these estimates with Monte Carlo simulations over a range of observing spans, signal strengths, grid resolutions, and good-time-interval structures, and show that the predicted uncertainties reproduce the run-to-run scatter of the recovered parameters in the tested regimes. We then apply the framework to AstroSat/LAXPC event lists for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4-3658. The results provide a practical and statistically motivated route to rotational-parameter estimation for targeted high-energy pulsar searches.

The origin of hot Jupiters remains a key open question. In the high-eccentricity migration scenario, traditional coreless models predict a strict tidal exclusion zone within $\sim 2.7$ tidal radii $r_\textrm{t}$, in which giant planets are either fully disrupted or ejected. We revisit this limit using three-dimensional hydrodynamic simulations of giant planets with realistic dense cores (10 - 20 $M_\oplus$). We find that even a few-percent-mass core fundamentally changes the outcome: \textbf{no total disruptions} occur within the previously suggested destruction zone ($\lesssim 2.7 \, r_\textrm{t}$). For deep encounters ($\lesssim 1.7 \, r_\textrm{t}$) planets suffer severe envelope stripping and are either progressively downsized to dense remnants or ejected after a few close encounters, possibly contributing to the free-floating planet population. In the intermediate regime ($ \sim 1.7 $--$2.0, r_\mathrm{t}$), planets experience significant partial mass loss over repeated encounters. For wider encounters ($ \gtrsim 2.0\, r_\mathrm{t} $), mass loss is minimal, allowing the planets gradually circularize into hot Jupiters. Furthermore, we show that for highly eccentric orbits ($e\gtrsim 0.9$), the change in specific orbital energy $ \Delta E_{\mathrm{orb}} $ depends primarily on periastron distance $ r_\mathrm{p} $ rather than semi-major axis $ a $. This enables us to extrapolate our fixed-$ a $ results across a broad ($a$, $e$) parameter space and identify a well-defined tidal ejection zone whose sharp boundaries converge asymptotically. Our results highlight the crucial role of planetary internal structure in high-eccentricity migration and suggest that the survival and transformation of core-bearing giant planets are far more common than previously thought.

Early-time flux excesses in Type Ia supernovae (SNe~Ia) offer a unique insight into their progenitor systems and explosion mechanisms. Although individual early-excess events and larger searches have been reported, demographic studies remain limited by sample size. We present a systematic search for early-time excess emission in a volume-limited sample ($z<0.06$) of SNe~Ia based on the Zwicky Transient Facility Data Release 2 (ZTF DR2). Using ZTF $g$- and $r$-band light curves, we identify candidates showing early-excesses shortly after the explosion time, and we apply conservative coverage and quality requirements to build reliable ``excess'' and ``no-excess'' bump and no-bump catalogs. From an initial sample of 1547 SNe~Ia, our final catalogs contain 42 early-excess and 110 no-excess events. We compare the two populations using SN and host environment parameters from ZTF DR2 and quantify the differences using two-sample statistical tests. We find the strongest differences are in SN light-curve properties: early-excess events have larger SALT2 stretch $x_1$ ($7.91\sigma$) and larger $r$-band secondary-maximum flux $\mathcal{F}_{r_2}$ ($6.25\sigma$), while differences in SALT2 color $c$ are weak ($0.57\sigma$). Early-excess events also favor bluer $(g-z)_{\rm local}$ ($3.41\sigma$) and lower $\log_{10} (M_*/M_\odot)_{\rm local}$ ($2.73\sigma$). Our results connect early excesses with SNe~Ia diversity, and motivate further analyses of upcoming larger samples.

We present an analysis of 17 years of Fermi-LAT observations of the microquasar SS~433. We detect four GeV sources in the region: a newly identified source, PS J1910+0550, located outside W50; the previously reported source J1913+0512; and two features, denoted as the East and West excesses, apparently associated with the X-ray lobes. We focus on the three sources located within W50. We do not confirm the previously reported periodic modulation from J1913+0512, as no significant periodicity is found in the full 17-year dataset. The East and West excesses exhibit distinct morphological and spectral properties, suggesting different physical origins. The East excess shows a hard spectrum with photon index $\sim1.7$, consistent with inverse Compton emission from relativistic electrons accelerated together with the particles responsible for the X-ray and TeV emission. In contrast, the West excess has a much softer spectrum with photon index $\sim2.6$ and is spatially offset from the known X-ray and TeV emission regions in the western lobe. The spectral shape and offset position of the West excess make it strikingly similar to J1913+0512. The emission from these two regions can be explained by GeV particles accelerated in SS~433, distributed throughout the source volume, and interacting with localized dense gas targets. Under reasonable assumptions regarding particle transport and energetics, both proton-proton and bremsstrahlung scenarios are viable, although the hadronic scenario is more naturally accommodated. These findings may therefore represent the first observational evidence for the acceleration of cosmic-ray protons in large-scale outflows from Galactic microquasars.

Many complex organic molecules (COMs) in star-forming regions are believed to form on dust grains. We thus expect both the reduced metallicity and dust-to-gas ratio in the outer Galaxy to have an impact on the chemical composition of these regions. We investigate if certain COMs are more sensitive than others to metallicity by measuring the chemical composition of hot cores in the outer Galaxy. We used NOEMA to perform an imaging spectral line survey of G135.27+2.79, located at a galactocentric distance of 13.1 kpc. We derived the rotational temperatures and column densities of the detected molecules while assuming local thermodynamic equilibrium and compared the chemical composition of G135.27+2.79 to other sources and to the predictions of the three-phase astrochemical code MAGICKAL. G135.27+2.79 hosts three continuum cores, labeled MM1, MM2, and MM3. Most species in MM1 trace a hot, compact region, confirming MM1 as a hot core. The chemical composition of MM1 correlates rather well with that of the inner and outer Galaxy hot cores G31.41+0.31 and WB89-789 SMM1, but its molecular abundances relative to methanol lie in between, which may reflect the influence of metallicity on COM formation. The model results agree reasonably well, though with a few notable exceptions, with the COM abundances of MM1 relative to methanol and with the abundance ratios between MM1 and G31.41+0.31. Sensitivity to the reduced metallicity and dust-to-gas ratio varies between molecules, with carbon chains and nitriles most negatively affected. The lower dust-to-gas ratio leads to slower adsorption under low-metallicity conditions so that more carbon is locked up into CO in the gas. Slow adsorption means that CO is hydrogenated more efficiently on grains, enhancing CO-related COM abundances above expectations. These results demonstrate that metallicity has a significant impact on the formation of COMs.

We present a multi-wavelength study of the massive protocluster G286.21+0.17 (G286) using \emph{JWST} near-infrared imaging and ALMA H$^{13}$CO$^{+}$(1--0) observations. The \emph{JWST} images uncover a compact ($\sim$0.5 pc) hub-filament system (HFS), comprising a dense central hub connected by at least four converging filaments seen in absorption, along with multiple H$_2$ protostellar jets/outflows. The hub hosts dense core G286c1. The H$^{13}$CO$^{+}$ emission confirms this HFS over [$-$19.2, $-$16.4]~km~s$^{-1}$. The \emph{JWST} images further trace prominent photodissociation regions around the H\,{\sc ii}~region~A, powered by a B-type star. The radial distribution of ALMAGAL 1.38 mm core properties reveals steep power-law slopes toward the hub center. Within the inner hub (r < 8'', $\sim0.1$~pc), the core number density follows $\rho~[\rm pc^{-2}] \propto r^{-2.4\pm0.5}$, the surface density scales as $\Sigma~[\rm g~cm^{-2}] \propto r^{-1.0\pm0.2}$, and the enclosed core mass varies as $M_{\rm core}~[M_{\odot}] \propto r^{-1.2\pm0.2}$, while core diameters remain approximately constant ($D_{\rm core}~[\rm AU] \propto r^{-0.1\pm0.1}$). These trends, along with filament mass accretion rates of $7\times10^{-6}$--$1.8\times10^{-4}$~$M_\odot$~yr$^{-1}$, support a competitive accretion scenario in which gravitational focusing enhances core growth toward the hub center. Filament linewidths increase from tail/outer-region to head/hub-region, consistent with gravity-driven turbulence. However, the absence of a preferred alignment between velocity gradients and gravitational force directions may indicate a dynamically evolved system. The HFS likely formed through large-scale gas layer interactions and compression by the adjacent H\,{\sc ii} region. Overall, star formation in G286 appears regulated by filamentary accretion, competitive core growth in the hub, and stellar feedback.

The NASA TESS mission is conducting a nearly full-sky survey, enabling the photometric characterisation of millions of stars. The forthcoming ESA PLATO mission will provide long-duration, high-precision photometry of tens of thousands of bright stars to be characterised through asteroseismology. The TESS Luminaries Sample is a catalogue of 196 bright naked-eye (V < 6) main-sequence (MS) and sub-giant (SG) stars exhibiting solar-like oscillations. Among them, the subset located within the PLATO Long-duration Observation Phase (LOP) fields constitutes an exceptional set of targets that will be observable by PLATO from the earliest phases of the mission, making them ideal calibrators during commissioning and the first months of science operations. This paper aims to provide an in-depth asteroseismic characterisation of 32 Luminaries stars that fall within the PLATO LOP fields of view. Individual mode parameters were extracted for the first time for 26 of them. We used three independent seismic pipelines, one of which is similar to the algorithms used in the official PLATO pipeline. Statistical criterion were applied to identify the optimal combination of data calibration, observing cadence, and fitting pipeline for each star. For all stars, we derived large and small separations, the asymptotic phase term, radial mode amplitudes, and mean linewidths per order. Comparisons reveal consistent trends in the seismic parameters, confirming the robustness of our analysis. In sub-giants, mixed-mode identification differs in the three pipelines, revealing extraction inconsistencies requiring longer datasets to improve our mode identifications. The Luminaries stars located in the PLATO LOP fields constitute a unique sample that will play a crucial role in validating, calibrating, and optimising PLATO's seismic performance.

Clean Hoag-like ring galaxies are often characterized by an old compact central component, a depleted gap, and a detached outer ring. We identify a boundary-supported radial-layering mechanism in a shell-deformed Kepler control model. A compact inner boundary supplies the core state, while a localized effective shell deformation, interpreted as the reduced imprint of externally supplied material settled near a finite circularization radius, needs to create only an internal maximum and a subsequent outer minimum. These act as the gap barrier and ring-supporting well. The onset of this structure is organized by a saddle-node threshold of the critical-point equation. In a 10^4-point Monte Carlo scan, shell-localized boundary-supported candidates occupy finite parameter volume under the adopted priors, and none of the localized candidates contains an ordered interior minimum--maximum--minimum subsequence. The same branch gives a scale-free gap-to-ring interval overlapping representative ratios for Hoag's Object, UGC 4599, and PGC 1000714, but not for the environmentally processed comparison object JO171.

Recent advances in artificial intelligence (AI) agents are pushing AI beyond tools toward autonomous scientific discovery. We discuss two complementary agentic systems for cosmology: \texttt{CMBEvolve}, which targets tasks with explicit quantitative objectives through LLM-guided code evolution and tree search, and \texttt{CosmoEvolve}, which targets open-ended scientific workflows through a virtual multi-agent research laboratory. As preliminary demonstrations, we apply \texttt{CMBEvolve} to out-of-distribution detection in weak-lensing maps, where it iteratively improves the benchmark score through code evolution, and \texttt{CosmoEvolve} to autonomous ACT DR6 data analysis, where it identifies non-trivial pair- and scale-dependent behaviour and produces analysis-grade diagnostics. These examples show how cosmology can provide both controlled benchmark tasks and realistic open-ended research problems for the development of AI scientist systems.

Circumbinary planets (CBPs) currently identified are in nearly coplanar configurations relative to their host binaries, yet the dynamical origin of this preference remains unclear. We investigate this question by simulating the secular spin-orbit evolution of CBP systems with tidal decay. A representative case shows that the system evolves through three stages (coplanarization, spin-orbit synchronization, and spin-orbit alignment) through the angular momentum exchange between stellar spin and orbital motion. The evolution of mutual inclination is strongly coupled to stellar obliquity. Phase-space analysis and examination of stellar Cassini states reveal that arbitrary initial inclinations are gradually damped to coplanarity by tides, while stellar obliquity is adiabatically captured into Cassini states with diminishing oscillation amplitudes. This study provides a self-consistent analytical and numerical framework for determining stellar Cassini states and understanding coupled spin-orbit evolution in CBP systems. It shows that tidal dissipation, combined with adiabatic capture into Cassini states, drives the observed dynamical behavior.

The spectrum and composition of ultra-high-energy cosmic rays (UHECRs) suggest that the population dominating above the ankle releases particles with an unusual hard spectrum at low rigidity, below the EV scale. In self-confinement scenarios, such an apparent hardening arises from transport: escaping UHECRs generate magnetic turbulence that delays their own release from the magnetized environments surrounding their sources. We extend the self-confinement scenario based on the non-resonant streaming instability to a mixed nuclear composition. We describe the confinement region with an effective leaky-box model including escape, photodisintegration, and secondary production. We then compare the resulting spectrum and composition with Auger measurements and compute the associated cosmogenic neutrino and gamma-ray emission. We find that self-generated turbulence can suppress the escaping flux below the EV scale for source luminosities and magnetic-field coherence lengths compatible with UHECR sources hosted in galaxy clusters and propagating through cosmic filaments. During confinement, heavy nuclei efficiently photodisintegrate, producing secondary protons that contribute below the ankle and help account for the observed composition. The predicted neutrino flux remains compatible with current limits, while the diffuse gamma-ray background provides a potentially strong constraint on the most extreme configurations.

Brown dwarfs, being transitional objects between giant planets and low-mass stars, possess dense, cool interiors that provide optimal conditions to explore non-standard physics. Capture and accumulation of dark-matter particles can alter the thermal, structural and dynamic of these substellar objects. We aim to apply a self-consistent two-fluid framework to model the internal structure of self-gravitating brown dwarfs and to quantify how the presence of a dark-matter component modifies their mass--radius relations and dynamical properties. The brown dwarf is modeled as a composite system of a baryonic fluid, described by a polytropic equation of state, and an independent dark-matter fluid. Both components are coupled through their shared gravitational potential in hydrostatic equilibrium. We solve numerically the coupled Lane-Emden equations for a range of dark-matter mass fractions. We find that dark matter accumulating in the core reshapes the baryonic density profile, modifying both the radius and the second-order Love number. Radius and dynamical anomalies in brown dwarfs can serve as diagnostic tools to constrain dark-matter properties. Future high-precision astrometric missions could identify these structural signatures, establishing brown dwarfs as possible detectors of dark matter in the Galaxy.

We present a deep DREAMS JWST/NIRSpec MSA medium-grating spectrum of MACS0416-Y1, a galaxy at $z=8.312$ with the highest-redshift ALMA dust continuum detection to date, in order to characterize its properties together with archival IFU and ALMA data. The deep NIRSpec spectrum reveals a broad H$\beta$ line with a width of $\sim1100$ km s$^{-1}$. We interpret it as a broad-line AGN whose line diagnostics are consistent with AGN activity across its clumpy structure, given the absence of little red dot signatures. MACS0416-Y1 clearly shows [OIII]4363 emission, suggesting a moderately low metallicity of $12+\log(\mathrm{O/H})=7.86^{+0.09}_{-0.08}$ ($0.15~Z_\odot$). The combination of [CII]158$\mu$m and dust continuum emission indicates low dust mass ratios of $\log (M_{\rm dust}/M_{\rm gas})=-3.60^{+0.29}_{-0.22}$ and $\log (M_{\rm dust}/M_{\rm metal})=-0.95^{+0.29}_{-0.20}$. Because the metallicity of MACS0416-Y1 is around the critical metallicity of $0.1\textrm{-}0.2~Z_\odot$, the system is expected to undergo dust growth, explaining these low dust mass ratios as well as its small dust mass, $M_{\rm dust}\sim10^6~M_\odot$. The intense UV radiation from the AGN may contribute to a high dust temperature of $T_{\rm dust}\simeq 91^{+62}_{-35}$ K, boosting the dust-continuum emission above the ALMA detection limit despite the small $M_{\rm dust}$ at $z>8$. We find a very high total flux ratio of [OIII]88$\mu$m/[OIII]5007 = $0.26 \pm 0.06$ in MACS0416-Y1, above predictions from single ionized nebular models at any electron density. This discrepancy suggests that the [OIII]88$\mu$m and [OIII]5007 trace largely distinct regions, with the optical line suppressed in dusty nebulae, and thus requires careful interpretation when combining optical and infrared emission lines in JWST+ALMA studies.

The environment plays a critical role in galaxy evolution, with galaxy clusters and their infall regions offering diverse conditions that shape galaxies before they enter the dense cluster core, a process known as ``pre-processing''. However, identifying environmental substructures, particularly galaxy groups in these transitional zones, remains challenging due to projection effects and ``fingers-of-god'' distortions. In this work, we present a supervised machine learning framework for classifying galaxies into three environmental categories: main cluster, group, and neither, using observable galaxy properties such as positions, line-of-sight velocities, and stellar mass. The model is trained on mock observations derived from cosmological simulations designed to replicate survey conditions and achieves an overall accuracy and class-size-weighted precision of 81%. The neither and main cluster populations are reliably recovered, whereas group galaxies remain the most challenging to identify, achieving 30% completeness and 76% purity. Within $1\times R_{200}$, classification performance is suppressed, but it improves beyond this radius, reaching 40% completeness and 80% purity. Resampling and thresholding strategies allow the model to be tuned toward either higher purity or higher completeness; in this study, we adopt first-past-the-post thresholding to emphasise purity. Model performance is consistent across cluster masses and dynamical states, and it outperforms both Friends-of-Friends and Gaussian Mixture Modelling. This flexibility makes it well suited to upcoming spectroscopic surveys of cluster infall regions, providing a robust statistical tool for disentangling environmental influences on galaxy evolution.

The upcoming observations from the Square Kilometer Array Observatory will provide the astronomical community with a wealth of observations of important objects at long wavelengths. Full analysis of these outputs will necessitate specialized methods and software. Using synthetic observations of protoplanetary discs as an example, we present a machine learning-based visibilities-informed reconstruction for enhanced observations (VIREO) method for denoising data. This method explicitly provides a denoising U-Net with the interferometric observation's point spread function as both an additional input and term in the model's loss function. VIREO outperforms traditional cleaning methods and PSF-ignorant denoising models by producing data that is quantitatively cleaner and more conducive to analysis of the planets within the disc. Applying VIREO to archival ALMA data creates images with significantly less background noise, while maintaining, and in some cases enhancing, the substructure. By demonstrating the general utility of visibility-informed models, our results suggest that VIREO is generally applicable across the interferometric observatories.

Radio waves propagating through the interstellar medium are influenced by variations in plasma density. For spatially localised plasma structures along the line of sight, time-delay Doppler analyses of pulsars often reveal scintillation arcs in the secondary spectrum, frequently exhibiting a parabolic morphology. In the thin-screen approximation, the arc curvature is commonly used to infer the distance to the plasma concentration, which is modelled - via Kirchhoff-Fresnel diffraction theory - as an effective phase screen imposed by the column density of a localised disturbance. Here, we identify several limitations of the thin-screen model that necessitate a fully three-dimensional treatment, without reducing the problem to a projected screen density. We show that the arc curvature can vary depending on the three-dimensional structure of the plasma, rendering it a less reliable indicator of distance. Moreover, when volume propagation is considered, asymmetries and a richer variety of features emerge in the secondary spectrum compared to those predicted by the thin-screen approximation. We conjecture that these phenomena are linked to the onset of branched flow produced by a sequence of weak but correlated scattering events.

Determining whether temperate rocky exoplanets orbiting M stars retain atmospheres is currently a central goal of exoplanet astronomy. To this end, the James Webb Space Telescope has begun searching for atmospheres on these worlds with MIRI secondary eclipse spectroscopy and photometry. Here, we develop a novel climate-constrained Bayesian inference framework that yields atmospheric pressure and composition constraints from these datasets, while accounting for planetary, stellar, and model uncertainties. Our approach fits observations with model spectra derived from self-consistent pressure-temperature profiles at radiative-convective equilibrium, thus maximizing the information extracted from the data and providing more robust inferences than retrievals that use parameterized pressure-temperature profiles. We demonstrate the framework on the existing MIRI LRS eclipse spectrum of LTT 1445A b (1.34 $R_\oplus$ and $T_{\mathrm{eq}} \approx 431$ K). An atmosphere does not need to be invoked to explain the data, meaning a bare rock model produces an adequate fit. If the planet has an atmosphere, the $2\sigma$ upper limits on surface partial pressures are $\lesssim 1$ bar for an optically thin gas like O$_2$, N$_2$ or CO, $\lesssim0.1$ bar for CO$_2$, $\lesssim 10^{-3}$ bar for H$_2$O, and $\lesssim 10^{-4}$ bar for SO$_2$. Scheduled MIRI F1500W observations could detect one of the thicker atmospheres permitted by the existing data (1 bar O$_2$ and 0.01 bar CO$_2$), if a precision of 20 ppm or better is achieved. This case study demonstrates that climate-constrained Bayesian inversion can turn rocky-planet eclipse spectra into the quantitative constraints necessary to test population-level atmospheric retention hypothesis, like the cosmic shoreline.

The search for an atmosphere on the closest rocky M dwarf planet, LTT 1445 A b, has been the subject of intense investigation from both the ground and space. Here, we present the first JWST transmission spectrum of LTT 1445 A b using a single visit spanning 3-5~$\mu$m using NIRSpec/G395H. We conduct two independent reductions of the data using both the Eureka! and ExoTiC-JEDI pipelines. Overall, we measure the NRS1 transit depths to a median precision of $\sim23$~ppm in 41 spectroscopic channels with uniform widths of 30 pixels ($\sim$ 0.02 $\mu$m), and the NRS2 transit depths to $\sim36$~ppm precision in 65 spectroscopic channels, also with uniform widths of 30 pixels. We rule out any statistically significant spectral features at this precision and place limits on atmospheric metallicity using a grid of chemical equilibrium models with grey opaque clouds. Using NIRSpec/G395H alone, we can place limits on the atmospheric metallicity of $\gtrsim350~\times$ Solar when the opaque pressure level is greater than 0.01~bars. We also conduct a combined analysis of JWST/NIRSpec and HST/WFC3 transmission data and find our atmospheric limits can be extended $\gtrsim500~\times$ Solar when considering both datasets. Future analyses both in transit and emission will uncover whether there are detectable atmospheric features.

Context. The spatially coherent multipath propagation of pulsar radiation leads to a temporal and spectral interference patterns called scintillation. It is caused by density variations in the ionized interstellar medium, which often take the form of thin scattering screens filled with multiple subimages of the pulsar. PSR B1508+55 is known to be scattered by one or two such screens. Aims. We investigate appropriate methods to achieve precise astrometry for a scattering screen from simultaneous observations of only two telescopes on a very long baseline without forming visibilities. Methods. Two simultaneous observations of PSR B1508+55 were performed with the 100-m telescope at Effelsberg and the Five-hundred-meter Aperture Spherical Telescope (FAST). Using and improving existing scintillometry techniques, we leveraged the evolving, very long baseline to precisely measure the screen orientation, effective velocity, and scintillation arc curvature. We inferred the one-screen and two-screen model parameters and we imaged the closer screen. Results. Each single epoch leads to much tighter angular constraints than long-term monitoring of scintillation arcs, revealing an ongoing evolution of the orientation of the closer screen. Images of the scattered pulsar were obtained with a resolution on the order of 0.1 mas. These results confirm the highly anisotropic alignment of the scattered images, while also revealing small-scale deviations from a large-scale straight line. Conclusions. We demonstrate that simultaneous observations of scintillation can be used as a powerful substitute for very long baseline inferometry.

Millimeter continuum emission and self-scattering polarization from protoplanetary disks are widely used to constrain dust properties. Interpreting these observations requires practical prescriptions for the disk emission. However, only approximate formulae are available for the continuum emission, and no widely applicable formula has yet been established for the polarized emission. We aim (i) to assess the validity of commonly used analytic approximations for the (sub)millimeter continuum emission from protoplanetary disks, and (ii) to derive realistic prescriptions for the disk emission for both the continuum and the polarization. We numerically solve the radiative transfer equation in an isothermal, constant-density plane-parallel slab, including dust absorption, emission, and self-scattering with full Stokes parameters. We find that commonly used analytic approximations for the continuum emission are systematically about 10 to 15% lower than our numerical solutions. Consequently, SED analyses of (sub)millimeter observations that adopt these formulae are likely to overestimate the optical depth (and thus the disk mass) and the dust temperature, and underestimate the albedo (and thus altering the inferred constraints on grain size). We also provide empirical fitting formulae that reproduce our numerical results for the continuum emission and polarization fraction. These formulae will enable observational data analyses to be carried out more accurately and efficiently than with the conventional approaches. For the analysis of (sub)millimeter observations, we recommend using our new empirical formulae or interpolation of our numerical results, rather than commonly used approximations.

Aims. We present a homogeneous observational study of electron temperature ($T_{\rm e}$) relations between ionic species: $T_{\rm e}$([N II]), $T_{\rm e}$([O II]), $T_{\rm e}$([O III]), $T_{\rm e}$([S II]), $T_{\rm e}$([S III]) and $T_{\rm e}$([Ar III]), using 699 spectra of Galactic and extragalactic H II regions and local star-forming galaxies (SFGs). Methods. We use the DEep Spectra of Ionised REgions Database Extended (DESIRED-E), comprising more than 3000 spectra with direct $T_{\rm e}$ determinations, selecting those with at least two $T_{\rm e}$ diagnostics. We recompute electron density ($n_{\rm e}$) and $T_{\rm e}$ using updated atomic data and a consistent methodology. The resulting $T_{\rm e}$--$T_{\rm e}$ relations are analysed using orthogonal distance regression, quantifying total and intrinsic dispersions and comparing slopes with previous works and photoionisation models. Results. Relations involving low-ionisation $T_{\rm e}$ diagnostics show large intrinsic dispersions, especially for $T_{\rm e}$([O II]) and $T_{\rm e}$([S II]), likely due to sensitivity to $n_{\rm e}$ inhomogeneities, recombination contributions, and uncertainties. In contrast, relations using $T_{\rm e}$([N II]) show lower dispersions, indicating that this diagnostic provides a more reliable estimate of the low-ionisation zone temperature when only higher-ionisation $T_{\rm e}$ diagnostics are available, despite observational difficulties at low metallicity. Overall, slopes agree with model predictions, particularly for relations with low intrinsic dispersion, such as those involving $T_{\rm e}$([N II]) and $T_{\rm e}$([S III]). These results provide a robust empirical basis for estimating $T_{\rm e}$ when limited diagnostics are available.

During near-Sun crossings of the heliospheric current sheet (HCS), Parker Solar Probe (PSP) observed populations of high-energy protons and heavier ions indicating possible energization by magnetic reconnection up to 10s -- 100s keV nucleon$^{-1}$. Here we study ion acceleration by magnetic reconnection at the HCS. To estimate ion energization, we solve the Parker transport equation coupled to a large-scale 2D MHD reconnection simulation. We find that multiple ion species develop power-law distributions with both spectral index and high-energy cutoff $E_{\text{max}}$ consistent with in-situ data. By accounting for the injection physics determined by kinetic simulations, we confirm that the charge-to-mass ratio scales as $E_{\text{max}} \propto (Q/M)^{\alpha}$ with $\alpha \sim 0.8-1.1$, approximately consistent with PSP measurements in the broader range $\alpha \sim 0.6-1.7$. In the limit where ions are injected at the same energy per nucleon, $\alpha$ can be as low as $\sim 0.3$. These findings further support the role of magnetic reconnection in producing high-energy heavy ions at the HCS.

The dependence of the stellar initial mass function (IMF) on star-formation environment, particularly at low metallicities and high redshifts, remains poorly constrained. Ultra-faint dwarf galaxies (UFDs) are local fossils of high-redshift galaxies hosting old, metal-poor populations, and their resolved stellar populations provide unique pathways to constrain the sub-solar IMF. We investigate the low-mass IMF in the Boötes I (Boo I) UFD with JWST/NIRCam, leveraging its capability to resolve over 10,000 stars reaching $\lesssim0.15 M_{\odot}$, obtaining one of the largest, deepest resolved stellar samples for UFDs. We explore three different functional forms of the IMF with machine learning and statistical techniques, combining forward modeling of synthetic color-magnitude diagrams with simulation-based inference. We find that a single power-law IMF fails to reproduce the observed luminosity function and also deviates from the canonical Salpeter IMF. Our best-fit broken power-law and lognormal IMF parameters are consistent with the Milky Way within 68% confidence level, providing evidence that star formation at metallicities as low as [Fe/H]$\approx-2.4$ follows a similar IMF as in the Milky Way. By treating Boo I as a local relic analogous to a high-redshift galaxy with a stellar mass of $\lesssim10^5 M_{\odot}$ at $z\gtrsim6$, our results provide evidence for the universality of the IMF across both local and high-redshift environments.

Tomographic redshift distributions from photometric data are crucial ingredients in cosmic shear analysis, since they are required for the theoretical calculation of the signal based on the redshift distribution of the galaxies where the shear field is sampled. In this paper, we develop as a proof of concept Leistedt et al.'s template-based Bayesian Hierarchical Model framework into an application to weak lensing data by sampling the redshift distributions of the galaxies in the KiDS+VIKING-450 survey. We also use a principal component analysis to provide a set of representative templates drawn from a large superset. For computational tractability, subsets of $10^5$ galaxies are chosen to determine the redshift distributions, and we test the sensitivity of the cosmological inference to the subset chosen, finding it to be subdominant compared to the statistical error. We marginalise over the inferred redshift distributions and find that the Bayesian method increases the clustering parameter compared with previous studies, alleviating the $S_8$ tension with Planck, where $S_{8}\equiv\sigma_{8}\sqrt{\Omega_{\tm{m}}/0.3}=0.756\pm 0.039$, assuming flat $\Lambda$CDM. The tension with Planck for this survey is reduced from $2.3\sigma$ to $1.9\sigma$. We also infer a value for the matter density, $\Omega_{\tm{m}}=0.31\pm 0.10$.

Accurate estimates of fundamental physical properties of galaxies, such as star formation rates (SFRs) or stellar masses, are essential for testing and constraining models of galaxy formation and evolution. Spectral energy distribution (SED) modeling has become the standard method for deriving these quantities. However, the influence of the underlying stellar population synthesis (SPS) models on the inferred parameters remains poorly quantified. This work investigates how the choice of SPS models affects the estimation of SFRs and stellar masses derived from SED modeling. Four widely used SPS models are applied to a sample of 17 230 galaxies with spectroscopic redshifts, selected from recently published Hubble Space Telescope and James Webb Space Telescope photometric catalogs. SEDs are modeled using the Code for Investigating GALaxy Emission. The analysis is performed in two steps: (i) estimating galaxy properties with each SPS model, and (ii) employing synthetic catalogs to assess the relative impact of model choice on the recovered parameters. Systematic differences are found among the models, with stellar mass estimates varying by up to ~ 0.6 dex and SFRs by up to ~ 0.4 dex between certain model pairs. The choice of stellar population model introduces significant systematic uncertainties in derived galaxy properties. This dependence should be accounted for when interpreting SED-based measurements and comparing results across different studies of galaxy evolution.

Solar pores possess irregular and evolving boundaries that are often far from the ideal circular flux tubes assumed in many magnetohydrodynamic (MHD) oscillation models. To tackle this irregularity in a consistent way for wave analysis, we introduce a framework that employs the convex hull of the pore boundary - derived from intensity minimum error thresholding - as the domain to perform further analysis. Using the modal assurance criterion, we find the range of pore shapes for which this approximation is valid. We demonstrate the usefulness of this framework by applying it to multi-height, high-cadence observations (4170 angstrom continuum, G-band, Na~\textsc{i}, and Ca~\textsc{ii}~K) of a solar pore, and apply Proper Orthogonal Decomposition of the convex hull to determine wave modes. The fundamental sausage ($m=1$) and kink ($m=2$) mode is found to remain reliable, while higher-order fluting modes ($m\ge3$) are strongly degraded by small-scale boundary irregularity. As expected, sausage-like modes dominate the variance at all heights and exhibit a systematic upward shift in frequency, consistent with freely propagating compressive waves. In contrast, the kink-like motions appear weak, confined to a persistent low-frequency peak, and most plausibly interpreted as a forced response to granular buffeting rather than a propagating mode. Together, these results establish a practical methodology for boundary-mode analysis in real, highly structured pores and provide new constraints on the nature and height evolution of MHD waves in the lower solar atmosphere.

We report the discovery of a new candidate pulsating ultraluminous X-ray source (PULX) in NGC 5128 (Centaurus A). The candidate, 4XMM J132542.2-425943, is a transient source, identifiable as a clear X-ray point source for $\sim 8$ months in 2014, during its only major recorded outburst. The source flux exceeded $10^{-12}$ erg cm$^{-2}$ s$^{-1}$ at the peak of the outburst. The long-term light curve of 4XMM J132542.2-425943 shows two further, less luminous detections in 2017 and 2024, but was otherwise in quiescence. This behaviour is similar to the class of pulsating transients with outbursts that reach the ultraluminous regime, which includes the well-studied Galactic PULX, Swift J0243.6+6124. However, 4XMM J132542.2-425943 displays a soft X-ray spectrum, making this source distinct from the existing population of PULXs, which typically show hard spectra below $10$ keV. We searched the 2014 XMM-Newton observations for X-ray pulsations, revealing coherent, sinusoidal X-ray pulsations at a frequency of $1.27$ Hz in one XMM-Newton observation (ObsID 0724060801), at a marginal significance. For this signal we measure a pulsed fraction, PF$\approx~15 - 17~\%$ and $\dot{f}~\sim~4~\times~10^{-9}$ Hz s$^{-1}$. However, we find that the intrinsic randomisation employed by XMM-Newton's Science Analysis Software, XMM-SAS, during the data reduction procedure introduces considerable uncertainty in the strength of our marginal pulsations, which varies significantly between consecutive data reduction iterations. We explore the impact of this randomisation and demonstrate that it can generate widespread false positives and false negatives, which, in the context of PULX searches, may cause viable candidates to be unnecessarily discarded or vice versa.

It is of a general interest to look for signatures of stellar bodies orbiting supermassive black holes (SMBHs) in galactic nuclei other than the Galactic center. Previously stellar transits were analyzed in UV, optical, and X-ray domains as well as potential microlensing signatures due to more compact bodies orbiting SMBH accretion disks. Here we complement previous studies by considering nearby ($z=0.001$) radio-loud active galactic nuclei targeted by different facilities in the millimeter domain. At these wavelengths the radio core is sufficiently small so that it can be occulted by large evolved stars in dense nuclear star clusters. We find that in the millimeter domain evolved stars with stellar radii of $\gtrsim 500\,R_{\odot}$ can cause eclipses with the relative depth of $\sim 10\%$. Typical recurrence timescales are at least 10 years and the eclipse durations are $\sim 10$ days. Towards lower frequencies the eclipse temporal profiles become shallower and broader while towards higher frequencies they are deeper and narrower. Although expected to be rare due to selection effects and evolved stars being prone to tidal disruption, recurrent eclipses of mm radio cores can be applied to infer SMBH masses and constrain the composition of the Nuclear Star Cluster of the host nucleus.

This chapter tells the still-unfolding story of extracting polarization signatures from general relativistic magnetohydrodynamics simulations of accretion disks. In some sense, this effort is premature as there are still very few results of this kind. Much more abundant are phenomenological models. Nevertheless, we feel now is the time to rally the community to this cause. Since the focus of this book is on X-ray polarimetry, we focus exclusively on simulations of accretion onto compact objects. Most of the relevant work so far has been on black hole accretion disks, though neutron stars are also viable targets for X-ray polarimetry. The focus of our chapter is on how X-ray polarimetry coupled with accretion simulations might help us better understand properties of the disks, coronae, and jets that are the dominant components of accreting compact sources. We briefly illustrate the promise of this technique by demonstrating how it has already been used in the case of the Event Horizon Telescope (using radio polarimetry). We also speculate about where this field may be heading in the near future.

Atmospheric compositions of sub-Neptunes and super-Earths are often interpreted as tracers of formation location relative to volatile ice lines. However, prolonged magma oceans can chemically equilibrate with primordial atmospheres and modify accreted volatile signatures. In this study, we couple a synthetic planet population from the Bern Generation III formation model to an extended global chemical equilibrium framework including sulfur and nitrogen chemistry, and compare accreted and equilibrated compositions for $\sim$ 1200 young planets shortly after formation ($\sim$ 40 Myr) formed inside and outside the water ice line. We find that interior-atmosphere equilibration systematically alters elemental ratios and molecular abundances. The atmospheric C/O ratio shifts relative to the accreted state and remains systematically higher for planets formed outside the ice line. Nitrogen-bearing species NH$_3$, N$_2$ are strongly depleted through dissolution into the silicate melt, while minor amounts of HCN are produced, leading to low atmospheric nitrogen abundances. Sulfur-bearing species remain more abundant than nitrogen-bearing species; during equilibration, accreted H$_2$S partitions into the interior and small amounts of SO$_2$ form, but overall sulfur abundances depend only weakly on formation location. Silicon-bearing gases (SiH$_4$, SiO) are generated in substantial amounts, with narrower distributions for planets formed outside the ice line. We identify atmospheric C/O, SiH$_4$, and H$_2$O as potential indicators of formation location, while nitrogen depletion emerges as a generic outcome of magma ocean equilibration. Comparison with characterized sub-Neptunes such as TOI-270 d, K2-18 b, and GJ 3470 b shows broad consistency with oxygen-dominated, metal-rich atmospheres shaped by interior-atmosphere exchange.

The potential to detect axion dark matter through astrophysical processes has shown high promise in recent years. We therefore expand on previous work studying the axion-to-photon conversion efficacy of neutron stars and the interstellar medium (ISM) in this endeavor, respectively. For neutron stars (NS), we examine the possibility of a background signal emanating from all NS magnetospheres in the galaxy. Using a heuristic Galactic model, we find a significant background signal emanating from such magnetospheres in the Milky Way. This signal, while weak in absolute power ($\gtrsim 1$ mJy sr$^{-1}$ from the Galactic Center, at 2 GHz), can be detected through new statistical techniques with current instrumentation like the Atacama Large Millimeter Array (ALMA) at high radio frequencies (200 - 950 GHz). These techniques make use of higher order statistics like spectrally-limited ($\sim 300$ km s$^{-1}$) increases in confusion noise levels and kurtoses of survey images, and also show promise for general population estimation techniques. For the ISM, we consider Primakoff processes between free electrons and axions, and derive typical signal strengths of $10^{-15}$ Jy sr$^{-1}$ $\cdot$ $m_a$/eV, with a local, cosmological upper bound of $10^{-8}$ Jy sr$^{-1}$ $\cdot$ $m_a$/eV. Hence, we find that any diffuse axion signal from the ISM and other, large-scale, astrophysical plasmas to be too weak to be detected with modern technologies. We therefore find that the best avenue towards detecting a potential quantum chromodynamics (QCD) axion dark-matter particle is through the radio imaging of large swaths of the Galactic Center and other regions where we expect large numbers of pulsars and neutron stars.

The theory of General Relativity predicts that, since massive bodies curve spacetime, light from a distant source would be deflected by a foreground massive object -- a phenomenon known as \emph{Gravitational Lensing}. Historically, the strength of deflection of light from background stars by the sun, during the 1919 solar eclipse, supplied one of the first proofs for the theory of General Relativity. However, it is only in the last few decades, with the advent of the Hubble Space Telescope and other large, ground-based facilities, that lensing has become a principal tool in modern astronomy. Lensing allows us to study both the matter content of the lensing bodies such as galaxies or clusters of galaxies, mainly dominated by the otherwise-invisible \emph{dark matter}, and the distant background sources that are being lensed by them. Strong gravitational lensing, where sources are substantially magnified and multiply imaged, is particularly useful to that end. The substantial magnification allows for a high-resolution view of the sources and to detect fainter and farther objects than would otherwise be possible; and image multiplicity helps in verifying the distance to them, and for studying variable or transient sources. Paired with the unprecedented capabilities of the James Webb Space Telescope (JWST), lensing now allows us to observe, detect, and study distant sources like never before. I summarise recent advances in strong-lensing applications and near-future prospects with JWST.

A direct coupling between the inflaton and Standard Model gluons can dynamically raise the QCD confinement scale during inflation, making the axion temporarily heavy and suppressing axion isocurvature perturbations. As inflation proceeds, the confinement scale relaxes, the axion becomes light, and late-time de Sitter fluctuations can generate the observed dark matter abundance. We analyze this mechanism without specifying an inflationary potential, instead parametrizing the background by $\epsilon(N) \propto 1/N^p$, where $N$ is the number of $e$-folds before the end of inflation. The single parameter $p$ distinguishes monomial models ($p=1$), standard plateau models ($p=2$), and ultra-flat plateau or hilltop-like models ($p\ge 3$). We (analytically) show that the mechanism selects plateau-like ($p\ge 2$) inflation: monomial models generically cause the confinement scale to grow too rapidly, while plateau models keep the QCD sector under perturbative control. In the minimal scenario, reheating occurs through the same inflaton-gluon coupling, and viable axion dark matter production is obtained when deconfinement occurs after the CMB window. The early-confinement sector also shifts the scalar spectral index to larger, bluer values, opening viable parameter space for models otherwise disfavored by CMB data.

We present a novel \textit{gauge-invariant and minimal} formation mechanism of primordial black holes (PBHs) in first-order phase transition (FOPT) and domain walls (DW) separately. This is based on the first-order tensor perturbations, generated during FOPT from bubble collisions \& sound waves, and from DW annihilation, sourcing curvature, at second-order in perturbation theory. We show that the PBH formation implies \textit{model-independent constraints} on FOPT parameters ($\beta/H, \alpha, T_{\star}$ ) and on DW parameters, ($\alpha_{\rm ann}, V_{\rm bias}, \sigma$), from existing PBH constraints. We find that asteroid mass PBHs can become the entire dark matter (DM) of the Universe, for $T_{\star} \in (4 \times 10^{2}, 10^{4})$ GeV, for $\beta/H \simeq 6$, involving $\alpha >\mathcal{O}(1)$ values. The corresponding FOPT Gravitational Waves (GW) amplitude will have its characteristic peak at $\Omega_{\rm GW}^{\rm p} h^2$ $\sim \mathcal{O}(10^{-8})$ between frequencies $f_{\rm p} \in ({10^{-5},10^{-2}})$ Hz which is within the reach in LISA and SKA detectors. PBH as entire DM is possible for $\sigma^{1/3} \in [10^{6}, 10^{8}]$ TeV, for $V_{\rm bias}^{1/4} \in [10^7, 10^{10}]$ MeV with the corresponding GW amplitude peak from DW annihilation $\Omega_{\rm GW}^{\rm p} h^2$ $\sim \mathcal{O}(10^{-9})$ (for $\alpha_{\rm ann} \sim 10^{-2}$) and peak frequencies between $f_{\rm p} \in (4 \times {10^{-4},10^{-1}})$ Hz with ($T_{\rm ann} \in 4.5 \times [10^3, 10^6] $) GeV within the reach in LISA and ET detectors. We also provide semi-analytical formulae for the tensor-induced density spectrum, $P_{\delta^{(2)}}$, $M_{\rm PBH}$ and $f_{\rm PBH}$, relating them in terms of FOPT and DW parameters which in turn, are related to viable particle physics origin of such FOPT and DW, and therefore, constrain such microphysics, either in the visible, or in dark sector models.

In this paper, we discuss the spherically symmetric gravitational collapse of matter fields in the de Sitter universe. The energy-momentum tensor of the matter field is assumed to admit a wide variety including dust, perfect fluids with equations of state, fluids with tangential and radial pressure, and with bulk and shear viscosity. Under different initial conditions imposed on the velocity and the density profiles, and by combining the results from exact analytical methods with those obtained from numerical techniques, we track the formation and evolution of spherical marginally trapped spheres as the matter suffers continual gravitational collapse. We show that the quasilocal formalism of trapped surfaces provides an ideal framework to study the evolution of horizons. More precisely, black hole and cosmological horizons may be viewed as the time development of marginally trapped surfaces.

Negative masses have long been explored, but their observational viability remains unclear. In this work, we develop a unified, observationally testable framework to constrain negative masses using both coupling level and dynamical probes. We establish that while dipole radiation bounds require universality of gravitational charge, the intrinsic dynamics of negative mass binaries generically lead to anomalous behaviors such as anti-chirps, dispersal and runaway motion. These signatures are absent in current gravitational wave observations, providing a robust exclusion channel independent of modified gravity assumptions.

Geomagnetic field reversal sequences exhibit persistence times spanning a broad range, from a few $10^4$ years to superchrons lasting more than $10^7$ years. Despite extensive observational and theoretical work, the physical mechanisms governing how such reversals occur and how their broad temporal variability is organized are still not fully understood. Here we investigate the temporal variability of geomagnetic polarity in a thermally driven low-dimensional geodynamo model subject to a slow periodic modulation of the control parameter governing the large-scale induction, namely the $\alpha$-effect parameter. We find that the modulation generates a multipeaked probability density function of magnetic persistence times, with local maxima occurring at approximately integer multiples of the modulation timescale, as expected in a stochastic-resonance-like regime. The peak positions follow an approximately linear dependence on their index, showing that the characteristic timescales selected by the system are set by the imposed modulation period. These results provide a physically motivated numerical framework in which slow modulation of a geodynamo control parameter can organize reversal statistics through stochastic-resonance-like dynamics.

External matter distributions can substantially reshape the strong field environment of compact objects, yet this effect is usually neglected in idealized isolated models. In this work, we investigate geometrically thin, optically thick relativistic accretion onto a static axisymmetric space-time that describes a slightly deformed compact object immersed in an external quadrupolar field as an exact solution of vacuum Einstein field equations. Our aim is to determine whether such locally geometries can produce distinctive accretion signatures and, more broadly, to identify the physically meaningful radial domain over which the local solution remains self-consistent. We show that the external quadrupolar distortion leaves a clear imprint on both orbital dynamics and accretion structure. We further find that the outer edge of the radiating region is closely tied to the transition between radiation pressure and gas pressure dominance, which may link the geometry to the thermodynamic properties of the flow. Therefore, the local nature of the distorted spacetime is not merely a formal geometric feature, but has observable consequences for the morphology and emission properties of accretion flows.

General relativity has achieved remarkable experimental and observational success. Critically, recent data from the LIGO-Virgo-KAGRA, Event Horizon Telescope, and GRAVITY collaborations are often credited with \textit{demonstrating} the existence of black holes, but in fact they only provide evidence for objects that should be regarded as black hole candidates. While current data are in striking agreement with the predictions for Kerr black holes, they can only rule out specific alternative models of compact objects rather than provide conclusive proof of black holes. More fundamentally, and independent of whether or not black holes exist, general relativity itself imposes limits on what can be observationally established. Essentially, no observational data is sufficient to confirm the existence of black holes.

We study the sensitivity of the diffuse high-energy neutrino flux observed in IceCube to new-physics effects resulting in an exponential flux attenuation along the trajectory, such as invisible neutrino decay or new interactions with the background encountered during propagation. We argue that, even though the sources and production redshifts of these astrophysical neutrinos are unknown, conservative energy-conservation arguments allow to severely constrain neutrino loss in most scenarios beyond the strongest existing bounds. By performing a fit to the High-Energy Starting Events from IceCube, we quantify the bounds and study their variation with the energy dependence of the attenuation, the assumed redshift distribution of the neutrino sources, and whether the attenuation affects neutrinos exclusively or no. We also show that including an energy-dependent attenuation at the level allowed in the fit may impact the determination of the spectral index of the diffuse flux.

We propose an explanation for the dark matter-baryon coincidence based on collapsing $\mathbb{Z}_N$ domain walls, which form a novel compact baryonic state: the baryoid. A baryoid has an asteroid-scale mass and up-to-nuclear-scale energy density, and can serve as a dark matter candidate. Starting from equal baryon numbers in the domains formed in the early universe, the collapse of the domain walls after the QCD phase transition leads to a baryon-number ratio of $(N-1):1$ between the false- and true-vacuum domains. Since baryons are slightly lighter in the false-vacuum domains than in the true-vacuum domain, the resulting dark matter-to-baryon energy-density ratio is naturally close to, but slightly smaller than, $(N-1):1$, or $6:1$ for $N=7$. We calculate the domain-wall dynamics and the efficiency of baryon-number trapping, derive the resulting baryoid properties, and discuss a broad set of phenomenological probes.

The interaction of electrically charged particles with magnetic fields is a fundamental problem in several areas of physics. An example is the motion of energetic particles through a magnetized plasma. The most accurate and reliable way to explore theoretically the interactions between particles and fields is via test-particle simulations. In such simulations one creates the turbulent magnetic field and solves the Newton-Lorentz equation numerically by employing an integration scheme. In the current article we discuss exponential integrators and derive systematically from this the Rodrigues scheme as well as the famous Boris integrator. For an approach where one creates the magnetic field anew at each time step, both integrators are overall comparable. In theory the Rodrigues approach should be more accurate due to the fact that the occurring matrix exponential is evaluated without further approximations. Practically, both methods provide very similar results. It is argued in the current article that a Rodrigues based integrator is a very strong alternative because for the specific problem discussed here, it does not require longer computing times.

Dark matter experiments are rare event search experiments that require zero background environment over very long exposures. To achieve this condition, a detailed simulation of detector geometry and experimental setup is required before the experiment is executed. Simulation plays a significant role in detector design and also provides a cost-effective and risk-free approach for predicting outcomes before real world experimentation. The present simulation work is focused on neutron background reduction for a dark matter direct detection experiment in India, the Indian Dark matter search Experiment (InDEx). The FLUKA and FLAIR simulation tools have been used throughout the simulation process. The experimental and simulation results available in the literature are being reproduced using FLUKA for validation purposes. The calibration and InDEx experiment are simulated, and the results are compared against the experimental results. For neutron background reduction in future experiments, the use of high density polyethylene (HDPE) is suggested and a shielding design using HDPE is presented. The results show that shielding reduces detector event rates by two orders of magnitude compared to the prior InDEx experiment without shielding.

We investigate the imprints on the angular power spectra of cosmological perturbations of a pre-inflationary bounce phase, as described by the hybrid and dressed metric approaches to loop quantum cosmology. For this purpose, we derive a new parametrization of the primordial power spectrum at the end of the inflationary regime. Apart from slow-roll coefficients and cosmological parameters that are present in the standard cosmological scenario without quantum modifications, this parametrization additionally depends only on pre-inflationary physics. More specifically, we find a dependence on the number of e-folds during the bounce epoch and on a characteristic suppression scale which, given the e-folds accumulated during cosmic evolution, is determined by the energy density at the bounce. Recall that this density depends on the Immirzi parameter and the area gap known from LQG. This leads to a robust and accurate parametrization of the primordial power spectrum. Since in pre-inflationary scenarios there is no preferred vacuum state, we adopt the NO-AHD proposal, which selects a vacuum that is optimally adapted to the background dynamics and yields a non-oscillatory primordial power spectrum. With this choice, we show that the tensor-to-scalar ratio in both quantization approaches coincides with its expression in the standard $\Lambda$CDM model when the observed scales are not much smaller than the power-suppressed region. Computing also the angular power spectrum, we find that, for a total cosmic expansion of about 140 e-folds, both the hybrid and the dressed metric approaches exhibit excellent agreement with Planck data at high multipoles, while apparently improving the fit with respect to $\Lambda$CDM for low multipole numbers.

Quantum Gravity remains elusive, largely because its observable effects are suppressed by powers of the Planck scale. Direct detection of single gravitons is widely believed to be impossible. Here we propose a concrete astrophysical mechanism that may overcome this suppression. We show that superradiant axion-like-particle clouds surrounding rotating black holes can generate multimode squeezed states of gravitons containing up to $10^6$ - $10^7$ correlated quanta. Such states exhibit distinctive polarization correlations and quantum-noise signatures that could be detectable in future gravitational-wave interferometers. Observation of these signatures would constitute direct evidence for the quantum nature of gravitational radiation. Conversely, their absence can place constraints on axion-cloud lifetimes. Our approach also provides a test of General Relativity as an effective field theory.

When considering how self-interaction affects an object's motion, it can be convenient to decompose the self-force into conservative and dissipative pieces. As a toy model for understanding such decompositions of the gravitational self-force, we consider objects that do not affect the spacetime, but are instead coupled to a nonlinear scalar field. There is then a standard splitting of the first-order scalar self-force into conservative and dissipative components. Multiple criteria can be used to obtain this splitting, all of which imply the same result. However, the implications of these criteria generically differ at higher orders. Demanding that any reasonable conservative sector be Hamiltonian, we identify multiple possible definitions of the conservative second-order self-force. Motivations for these possibilities and their properties are discussed and relevant Hamiltonians are obtained. We assume the existence of a three-point function with certain properties that is a generalization of the Detweiler-Whiting two-point function. These results apply to the two-body problem but are restricted to unbound scattering trajectories, due to infrared divergences that arise for bound orbits.

In recent years, there have been increasing challenges to the cosmological principle, based on new observations of e.g. supernovae and the cosmic bulk flow. As a result, the cosmological community is speaking their concern for the cosmological principle, and from which scales onwards it should apply. In this context, there is a desire to understand more fully the properties and signatures of cosmologies not obeying the cosmological principle. In this article, we let go of the demand of cosmic isotropy, and instead assume only spatial homogeneity in our cosmological models. We follow the results of our previous works [see citations in the list of references], and here bring these together into one unified picture, with the goal of describing the signature(s) of anisotropy in anisotropic cosmological models. We first introduce the Bianchi models -- a particular instance of spatially homogeneous cosmologies -- and show that a metric can be constructed for them when an appropriate collection of desired Killing vector fields is supplied. Then, we give the perturbations of the Friedmann equations in such Bianchi models, in the Newtonian gauge, derived using much the same methodology as applicable to the FLRW models. We show these can be combined into one characteristic partial differential equation. Finally, we use this equation in order to simulate the CMB of a toy Bianchi V example and produce its power spectrum. We close with a discussion, and suggestions for further research.

We investigate a scalar-vector-tensor theory in which matter is minimally coupled to a Jordan-frame metric, while a massive vector sector interacts with the baryonic current. We show that the conformal scalar coupling modifies the physical expansion rate measured by matter observers, leading to a late-time enhancement of the effective Hubble constant. By constructing a phenomenological scalar evolution that becomes relevant only at low redshifts, the model provides a purely late-time mechanism for alleviating the Hubble tension without significantly affecting early-universe cosmology. The scalar potential naturally acts as a dynamical dark-energy sector, while the vector contribution behaves effectively as a pressureless component at cosmological scales through a density-dependent vector mass. Hence, the framework connects late-time scalar dynamics, effective dark-energy evolution, and Hubble-tension alleviation within a unified setup. Finally, local gravitational constraints can be suppressed through a chameleon-type screening mechanism, allowing the theory to remain compatible with Solar-System tests while retaining nontrivial cosmological effects.

We point out that a simple late-time cosmological model where our Universe can absorb "baby universes" explains the exponential expansion of our universe without the need of a cosmological constant and leads to a z-dependence of the parameter w(z) in the equation of state. In this model w(z) is less than -1 for z sufficiently large.

Neutrino oscillation experiments present anomalous results across a vast range of baselines and energies. Here we show that a 3+1 scenario in which sterile neutrinos feel a novel matter potential $V_s$ proportional to background density of ordinary or (asymmetric) dark matter is able to explain several anomalies. At low-energies ($E\lesssim$ 1 TeV) the model behaves as an effective 3-flavor NSI-like scheme among active flavors and eliminates the tension between the two LBL experiments NOvA and T2K provided that the potential is negative and the two sterile mixing angles $\theta_{14}$ and $\theta_{24}$ are non-zero. A further indication in favor of a negative non-zero potential comes from the anomalous excess of $\nu_e$-like events observed in Super-Kamiokande atmospheric neutrinos, which, in the new scenario is explained by a modification of the 3-flavor resonance at few GeV. A high energies ($E\gtrsim $ 1 TeV) the new framework reveals its 4-flavor nature and produces a resonant behavior at $E \simeq$ 10 TeV as hinted at by IceCube. We identify an irreducible 3-level dynamics generating a new resonance in the $(\nu_e, \nu_\mu)$ sector intertwined with two conventional resonances in the $(\nu_e, \nu_s$) and $(\nu_\mu, \nu_s)$ systems. The novel amplification mechanism manifests with the emergence of effective mixing angles in matter ($\theta_{12}^m$ or $\theta_{13}^m$) involving active neutrinos. The scenario requires values of $f = V_s/|V_{NC}| \sim -20 $, $\Delta m^2_{41} \sim 60 $ eV$^2$, $|U_{e4}|^2\simeq \sin^2\theta_{14} \simeq 0.01-0.03$ and $|U_{\mu4}|^2 \simeq \sin^2\theta_{24}\simeq 10^{-4}-10^{-3}$. Such a very small size of $|U_{\mu4}|^2$ eliminates the tension between IceCube and the other $\nu_\mu$ disappearance searches. The model can be directly probed by KATRIN, which is very sensitive to the electron-sterile neutrino admixture in the region of high $\Delta m^2_{41}$.

We formulate axion magnetohydrodynamics beyond the ideal limit, retaining axion inertia and the essential physics of non-ideal plasmas from first principles. In this framework, regions where magnetic flux freezing breaks down acquire a new physical role: whenever $\mathbf{E} \cdot\ \mathbf{B} \neq 0$, magnetic dissipation acts as a localized source of axion radiation. We show that magnetic reconnection naturally excites mixed Alfvén-axion modes, enabling coherent energy exchange between magnetic fields and axions in magnetically dominated environments. In neutron stars and magnetars, this mechanism leads generically to transient axion bursts powered by reconnection--driven Alfvénic dissipation. We connect this production process to observational prospects and derive a characteristic sensitivity to the axion--photon coupling, complementary to searches based on static magnetic fields.