Papers for Tuesday, May 12 2026

We present the Microlensing Object high-Resolution Imaging Analysis pipeline, or MORIA. This is an automated procedure to reduce high-resolution HST images of microlensing targets, build empirical point-spread function models from the data, and perform simultaneous multi-star PSF fitting to blended sources, lenses, and neighbor stars. We have developed and tested this pipeline using HST observations of the microlensing event KMT-2019-BLG-0253, where we determine a host mass of $M_{host} = 0.65 \pm 0.04M_{\odot}$. We have reduced the number of possible solutions for this target by a factor of two, with the remaining solution subject to the well-known close-wide degeneracy. We determine a planet mass of $m_{p} = 7.18 \pm 0.40 M_{\oplus}$ (close) or $m_{p} = 9.48 \pm 1.13 M_{\oplus}$ (wide), and distance to the lens system of $D_L= 2.64 \pm 0.22$ kpc. This work demonstrates the importance of using an automated high resolution imaging tool to inform light curve modeling for microlensing planets found during the upcoming Nancy Grace Roman Galactic Bulge Time Domain Survey (GBTDS).

Environmental effects are a primary driver of elevated gas-phase metallicities in galaxies around massive clusters, but the underlying physical mechanisms for this satellite metallicity enhancement (SME) are still unclear. Using the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, we present the first measurement of the average SME as a function of projected cluster-centric distance. The resulting profile reveals three distinct regimes: a steep decline from the cluster center, a plateau near the cluster boundary, and an extended downturn across several cluster radii. Remarkably, the complex shape and amplitude of this observed SME profile are successfully reproduced in the EAGLE cosmological simulation. Drawing insights from EAGLE, we develop a novel satellite chemical evolution model to decompose the observed SME into physical contributions from suppressed star formation, stellar mass loss, and enriched gas inflow. Our analysis shows that continuous accretion of enriched intracluster medium dominates the SME plateau within the cluster virial radius, while mass loss and quenching jointly drive the rapid metallicity decline in the cluster core. Our method disentangles the impacts of three environmental processes on galactic chemical enrichment in the cosmic web, providing a powerful framework for understanding cluster galaxy evolution with current and future spectroscopic surveys.

Perseus is the brightest X-ray strong cool-core (SCC) cluster, with a bright central radio and $\gamma$-ray source plus low-frequency radio mini and giant halos. It is the archetype of the cooling flow (CF) problem, with X-rays implying mass cooling rates orders-of-magnitude larger than observed in other channels. Recent work suggested that ancient ($\gtrsim$\,Gyr-old) cosmic ray (CR) halos (ACRHs), injected by the central source, would produce thermal-like soft X-ray inverse-Compton (CR-IC) emission 'boosting' the CC and alleviating the CF problem. We examine Perseus and show that a simple model of CRs injected by NGC 1275 (+satellites) simultaneously accounts for the excess CF luminosity and minihalo. The models reproduce Perseus's soft X-ray surface brightness and X-ray inferred density/temperature/pressure/metallicity/cooling time/mass deposition rates; $\gamma$-ray spectra; extended hard X-rays; and radio surface brightness and spectral index data, from kpc-Mpc. These also reproduce independent constraints on magnetic field strengths and mass/potential models. The evolution of the minihalo spectral index and surface brightness are predicted by an aging population of CRs boosting the apparent SCC luminosity via CR-IC, and match well the observed hard X-ray slopes. The 'giant' low-frequency halo can be predicted by the sum of ACRHs around satellites distributed throughout the cluster, dominating diffuse synchrotron at $\gtrsim 100\,$kpc. Re-acceleration is neither needed nor important in these models, and implied CR transport speeds are consistent with buoyant advection. Previous claims of upper limits to non-thermal X-rays and CR pressure relied on strong assumptions which are not valid at the CR energies of interest, e.g. a power-law spectrum of CRs. This could resolve many historical puzzles about Perseus, and makes new predictions for future observations.

Using Extreme Value Statistics applied to the observed galaxy stellar mass and the UV luminosity functions, we empirically estimate masses and luminosities of the most extreme galaxies in cosmological surveys, including the full sky. We incorporate uncertainties in stellar mass measurements (Eddington bias) and the scatter in the stellar-halo mass relation to derive empirical limits for galaxies residing in the most massive halos. The maximum observed $M_\ast$ strongly depends on survey area and redshift, ranging from M_\ast \sim 7 \times 10^{12} \, M_\odot for full-sky surveys at $z\sim0$ to M_\ast \sim 10^{10} \, M_\odot at $z\sim16$. Massive galaxies, particularly at high redshift, approach the theoretical maximum baryonic mass available in halos, M_\ast \sim 0.16 \times M_{\mathrm{vir}}, consistent with previous claims. Accounting for measurement uncertainties significantly reduces the inferred maximum $M_\ast$ by up to $\sim1$ dex at $z\gtrsim10$, yielding stellar masses consistent with M_\ast < 0.16 at all redshifts. Assuming a perfect rank-order correspondence between the most massive halos and galaxies would guarantee this inequality at all redshifts. At 2 \lesssim z \lesssim 6, the most massive galaxies have stellar masses comparable to the total cold gas reservoir from cold and cooling flows, suggesting near-maximal star formation efficiencies, SFEs. At higher redshifts, halos are predicted to host galaxies undergoing starburst phases. When accounting for dust attenuation and adopting empirically inferred SFEs, we find good agreement between the model and the brightest observed UV galaxies at high redshifts. At lower redshifts, however, observed UV galaxies are too bright. Overall, our results indicate that current observations remain broadly consistent with $\Lambda$CDM once statistical and observational effects are properly accounted for.

We identify compact groups of galaxies (CGs) in the IllustrisTNG-300 simulation using a Friends-of-Friends (FoF) algorithm. Our approach is designed to be comparable to systematic CG searches based on spectroscopic surveys, while avoiding the conventional Hickson selection criteria, which can bias samples toward relatively low-density environments. We construct two CG catalogs: one based on a three-dimensional distance linking length of 73 kpc (i.e., $50~h^{-1}$ kpc), and another based on projected and radial linking lengths of 73 kpc and $1000~\rm km~s^{-1}$. We refer to these as the position-position-position (PPP) and position-position-velocity (PPV) CG catalogs, respectively. The PPV catalog provides a direct analog to observed CG samples. At $z = 0$ in TNG300, we identify 383 PPP CGs and 1666 PPV CGs. A large fraction ($\sim 80\%$) of PPV CGs are not physically compact systems but are contaminated by line-of-sight interlopers. We demonstrate that the scaling relation between total group stellar mass and velocity dispersion is an effective diagnostic for identifying false positives with line-of-sight interlopers. We further examine the large-scale environments of CGs and show that they reside in a wide range of densities, including the central regions of galaxy clusters. These CG catalogs provide a robust foundation for studying the formation and evolution of CGs in cosmological simulations.

Our knowledge of the space velocity of Messier 87, which is the dominant galaxy in the Virgo cluster, has been limited to the radial velocity component. Using a cadence of precision position measurements with the global geodetic very long baseline interferometry (VLBI) system over 28 years, we determined the proper motion vector of the radio-emitting core by a robust statistical method involving 1-norm optimization and bootstrapping. The proper motion vector is directed at a position angle $189.2\degr \pm 3.5\degr$ in the equatorial International Celestial Reference Frame, and its magnitude is $10.19$ $\mu$as yr$^{-1}$ with an uncertainty of $0.64$ $\mu$as yr$^{-1}$. The projected velocity of the AGN in the tangential sky plane is ($787\pm50$)~km~s$^{-1}$. The peculiar velocity of Messier 87 with respect to the preferred rest frame of the cosmic microwave background field is approximately 1037 km s$^{-1}$ (assuming a distance of 16.1 Mpc) with an angle of 65$^\circ$ to the current line of sight, which implies a tangential relative motion of M87 and the Galaxy. The peculiar velocity of M87 is directionally concordant with the reconstructed and $\Lambda$CDM-simulated motion of the Virgo filament towards the Great Attractor, but the Milky Way moves slower by 470 ~km~s$^{-1}$ in that direction.

We present 30 newly discovered repeating fast radio burst (FRB) sources from the second catalog of bursts detected by the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB). These repeaters have extragalactic dispersion measures (DMs) spanning $99.4-1446.0\ \text{pc cm}^{-3}$ and burst rates between $10^{-5.7}$ and $10^{-0.5}$ hr$^{-1}$ scaled to a fluence threshold of 5 Jy ms. We report evidence of monotonic, linear DM variations in four repeaters on years-long timescales. The newly discovered sources bring CHIME/FRB's total number of observed repeating FRBs to 80, 79 of which were discovered by CHIME/FRB, between 2018 July 25 and 2023 September 15. In the full CHIME/FRB sample, only 2.4$\pm 0.4\%$ of sources have been observed to repeat, and we do not find evidence for significant evolution of this value over the duration of the experiment. We find no substantial evidence for bimodal populations of one-off and repeating FRBs in their burst rate distributions; the distribution of upper limits on repeat rates implied from observations of as-yet one-offs is entirely contained within the observed range of repeater burst rates and the distributions do not appear inconsistent. Similarly, using the population analysis framework of C. W. James (2023), we find that our observations of repeating and yet-one-off FRBs are equally well fit assuming a power-law distribution of repeat rates with 50$-$100% of the population repeating.

The integral polarization of spiral galaxies in the radio band has been proposed as a new tracer of the intrinsic galaxy shape that augments lensing shear measurements. We revisit the method of shear estimation in this context. We introduce a new statistical model in which galaxy shape and polarization are Gaussian random variables with their covariance characterizing the quality of polarization-shape alignment. Applying the principle of likelihood maximization, we then analytically derive unbiased, minimal-variance estimators, which allow to simultaneously estimate gravitational shear, intrinsic shape alignment and line-of-sight polarization rotation, all at once and accurate to first order in these three effects. New to the literature, our estimators have the merits of being free of biases, robust in situations of few galaxies or poor polarization-shape alignment, allowing analytic reconstruction noise covariance, and minimizing uncertainties in power spectrum estimation, thus resolving conceptual issues of the existing estimation methods. This new analytic framework is generally applicable to future research that exploits the polarization-shape alignment effect of galaxies.

We present a catalog of 30 stars that are candidates for KH 15D-like binary systems, in which the observed brightness variations are caused by a circumbinary dusty disk that periodically obscures at least one of the stellar components as it moves along its orbit. Thanks to the regular observations conducted within the Optical Gravitational Lensing Experiment (OGLE) project, we provide unique light curves in the I and V bands with very long time baselines, in some cases beginning as early as 1997 and extending to the present day. Such long-term monitoring allows us to identify changes in eclipse widths, amplitudes, and light-curve shapes on timescales of many years. We highlight several circumbinary disk occultation (CBO) systems of particular interest and present spectra for three of them.

In this paper, we introduce the radio frequency interference monitor deployed at the Dominion Radio Astrophysical Observatory. It provides 2 GHz of instantaneous bandwidth, supporting channel bandwidths as fine as ~100 Hz for 1 s integrations, or integration times as low as ~50 ms for the standard 3.33 kHz channel bandwidth. After operating as a prototype instrument for several years, the monitor was commissioned to improve the calibration method, analog section temperature, and gain stability. It now operates both as a transient detector and as a long-term radio environment characterization tool. We introduce novel applications for the monitor and derive a new method for calculating the effect of gain drift on integrated data.

The identification and tracking of stellar feedback-driven galaxy bubbles is an important topic in star formation and galactic structure research. However, current observational analysis of bubbles is limited in scope; information on bubble lifetime is inaccessible. Simulation data thus provides a unique opportunity to glean some of these characteristics at high resolution. We present an investigation into the characteristics and evolution of hot, ionized bubbles in the interstellar medium of a dwarf spiral (NGC300-like) galaxy. We calculate the average radius, lifetime, temperature, density, and spatial distribution of the simulated feedback-driven bubbles using Lagrangian gas parcels, and we examine the relationship between these characteristics and the local galactic environment. We find exponential distributions of bubble lifetime and size, and we find a positive correlation between bubble lifetime and galactocentric radius. Finally, we predict how the data would appear in H$\alpha$ tracers and compare the simulated values to observations. We find an additional positive correlation between the size of the bubbles and the galactocentric radius using their H$\alpha$ tracers.

We present an analysis of specific star formation rates (sSFR) across the galaxy merger sequence using visual classifications from the Zooniverse citizen science project "Cosmic Disco: Characterizing Galaxy Collisions". Our sample comprises 4884 galaxy systems pre-selected as merger candidates from SDSS DR17 ($0.01 < z < 0.05$, $M_* > 10^{8.5}M_\odot$) using Zoobot, of which 3690 were classified as mergers spanning pre-interaction through post-coalescence stages by citizen scientist volunteers. We find a weak but statistically significant positive correlation between $\log(\mathrm{sSFR})$ and visual merger stage ($r = 0.161$, $p = 7.23 \times 10^{-23}$), with a best-fit relation $\log\left(sSFR\right)=(0.148\pm0.015)\, S_{\rm Merg}-(1.865\pm0.038)$. The large RMS scatter (0.661 dex) reflects visual merger stages capturing wide merger timescales, and our results corroborate previous findings of increasing SFR enhancement with merger progression. This work shows that citizen science is a viable complement to automated and pair-based approaches to evaluate timescales for galaxies across the merger sequence.

We present a fully Bayesian, pixel-based Doppler imaging framework that enables the simultaneous inference of surface brightness maps and geometric parameters, including the inclination $i$ and equatorial rotation velocity $v_{\mathrm{rot}}$, from high-resolution spectral time series. We treat the inference as a Bayesian linear inverse problem conditioned on nonlinear geometric parameters. The surface map is modeled as a Gaussian Process prior over pixel intensities, introducing a characteristic spatial scale that sets the map resolution. This allows analytical marginalization of the linear coefficients and efficient sampling of the nonlinear parameters with Hamiltonian Monte Carlo. {Validation with synthetic data demonstrates that our method recovers the longitudes of large-scale surface inhomogeneities and constrains $v_{\mathrm{rot}}$ and $i$ under the adopted model assumptions, while also revealing the limited latitudinal sensitivity intrinsic to Doppler imaging.} We applied this framework to high-resolution VLT/CRIRES observations of the brown dwarf Luhman 16B. Our analysis reveals a large-scale dark region at mid-latitudes, consistent with previous studies but now with spatially resolved uncertainty estimates. Furthermore, we successfully constrained the geometric parameters without fixing \(v_{\mathrm{rot}}\sin i\) or $i$ to literature values, deriving an inclination of $i = 61.0_{-12.3}^{+14.3}$ degrees and an equatorial rotation velocity of $v_{\mathrm{rot}} = 31.2_{-3.1}^{+5.3}~\mathrm{km\,s^{-1}}$. These results indicate a radius broadly consistent with evolutionary models and suggest a possible spin-axis misalignment under the assumption of comparable equatorial rotation velocities for the two components. Our code is publicly available under the MIT license.

We investigate the cosmological implications of an extended gravitational framework based on biconnection gravity, constructed from the Schr$\ddot{o}$dinger connection and its dual. In this approach, the difference between the two connections defines the mutual curvature, which encodes the non-Riemannian geometric degrees of freedom, while their symmetric combination reduces to the Levi-Civita connection and hence reproduces general relativity at the background level. Within this setting, we derive the generalized Friedmann equations for a spatially flat Friedmann-Lemaître-Robertson-Walker Universe. The resulting equations contain additional geometric contributions that may naturally encode an effective dark energy sector induced by the biconnection degrees of freedom. We explore this extra dark energy by adopting five commonly used parametrizations, namely B$\Lambda$CDM, $\omega$CDM, Chevallier-Polarski-Linder, Barboza-Alcaniz, and a logarithmic equations of state. These considerations are confronted with recent observational data, including DESI DR2, Pantheon$^+$, and CC observations. Our analysis shows that the four parameterizations enter the acceleration phase at almost the same redshifts and share the same current value of the Hubble rate. Furthermore, the statistical comparison based on the Akaike, Bayesian, and Deviance Information Criterion shows that Barboza-Alcaniz, and logarithmic parameterizations have strong evidence and are competitive with $\Lambda$CDM. To classify this biconnection gravity in the plethora theoretical models describing the current cosmic acceleration, we examine its implications through cosmographic tools, including the deceleration, jerk, and snap parameters, as well as through the Statefinder analysis and $Om(z)$ diagnostic.

We propose a non-parametric approach to construct the statistical equation of state (EOS) continuously from the nuclear crust to the asymptotic-freedom regime. Driven by the observationally required stiffening to support two-solar-mass neutron stars (NSs) with relatively small radii for low-mass NSs, this global thermodynamic constraint suggests a clear peak of squared sound speed ($c_s^2$) in massive NSs. To prevent overshooting perturbative QCD (pQCD) energy-density bounds, this early stiffening must be actively compensated by an extended density range of softening, with $c_s^2$ not approaching $1/3$ until $\sim\!30\,n_{\rm sat}$. Consistently, the trace anomaly $\Delta \equiv 1/3 - p/\epsilon$ becomes positive beyond NS densities and approaches the pQCD limit from above. This natural emergence of $\Delta > 0$ organically aligns with some anticipated microphysics, likely arising from the pressure dilution in a quark-hadron mixed phase, non-conformal pQCD corrections to quark-gluon interactions, or the symmetry-breaking effects of finite strange quark mass. By measuring the degree of this non-monotonic behavior in the posterior, we find evidence for a hadron-quark phase transition in the cores of the most massive neutron stars. This indicates that the non-perturbative quark matter is intrinsically soft, fundamentally distinguishing it from the stiff scenarios associated with the quark-star picture.

The discovery of the gravitational-wave event GW170817 from a binary neutron star merger, together with its multi-wavelength electromagnetic counterparts, marks the beginning of the era of multi-messenger gravitational wave astronomy. Observations of gravitational-wave signals from compact binary mergers enable an independent measurement of the luminosity distance to the source. This implies that gravitational-wave sources can serve as standard sirens to probe the expansion history of the Universe, providing a new approach to constrain cosmological parameters. In this paper, we review the basic principles of using gravitational-wave standard sirens to constrain cosmology. We discuss various methods for determining the source distance and redshift, as well as the capabilities of second and third generation ground-based detectors and space-based detectors in constraining cosmological parameters, especially the Hubble constant and dark energy parameters. By examining two types of standard sirens, binary neutron star mergers with electromagnetic counterparts as bright sirens and stellar-mass binary black hole mergers as dark sirens, we illustrate the methodology, challenges, and future prospects of the standard siren approach.

Recent high-precision atmospheric observations with JWST is enabling detailed characterization of sub-Neptune atmospheres and motivating efforts to understand and constrain their interiors. Theoretical studies suggest that sub-Neptunes possibly host hydrogen-dominated atmospheres that are chemically coupled with an underlying magma ocean. However, a quantitative retrieval framework directly linking atmospheric spectra to magma ocean properties has yet to be established. Here we introduce MELTYQ, a coupled magma-atmosphere retrieval framework that links transmission spectra to the oxidation state and volatile inventory of underlying magma oceans. MELTYQ combines a magma-atmosphere equilibrium model, which includes the solubility of H-/O-/C-/N-bearing species in the melt and redox reactions, with a Bayesian spectral retrieval scheme. Using simulated retrieval tests, we validate the approach and show that magma redox state and volatile content can be constrained under favorable observational conditions. As a proof of concept, we apply MELTYQ to JWST transmission spectra of the benchmark sub-Neptunes K2-18 b and TOI-270 d. We find that coupled magma-atmosphere retrievals are generally capable of reproducing the observed spectra of these planets. However, we identify several key limitations in the current framework. Specifically: more flexible free-retrieval approaches remain statistically preferred; the CO/CO$_2$ absorption feature near 4.5 $\mu$m for TOI-270 d is not fully captured; and a number of underlying model assumptions may not be strictly valid. Nevertheless, embedding coupled magma-atmosphere models directly within Bayesian retrievals enables quantitative assessment of degeneracies and sensitivities, establishing a pathway for directly connecting atmospheric spectra to magma composition in this underexplored exoplanet regime.

Blazars are promising candidates for high energy neutrino sources, yet the physical origin of their neutrino emission remains uncertain. In this work, we extend our previous study by modeling the broadband spectral energy distributions (SEDs) of 103 neutrino blazar candidates (NBCs) within a hadronic framework. To estimate the maximum possible neutrino output, we adopt an assumption in which the high energy emission is dominated by p gamma interactions and the contribution from leptonic inverse Compton scattering is strongly suppressed. From the SED modeling, we constrain nine key parameters describing the emission region and particle energy distributions. We perform a partial correlation analysis to investigate the relationship between neutrino luminosity and electromagnetic emission, and we found a weak or moderate correlation between optical R band and neutrino emission. Our model predicts prominent proton synchrotron emission peaking in the MeV band for most sources, with 99 out of 103 NBCs exhibiting proton synchrotron peaks within 0.1 to 100 MeV, highlighting the MeV band as a key window for distinguishing between leptonic and hadronic scenarios. Based on the model-predicted maximum neutrino fluxes, we find that three NBCs are potentially detectable by IceCube, while up to 22, 45, and 62 sources may be detectable by KM3NeT, NEON, and TRIDENT, respectively. These results provide testable predictions for future multi-messenger observations and offer new insights into the composition and radiation mechanisms of blazar jets.

Filamentary structures are ubiquitous in the interstellar medium, yet the extent to which magnetic fields influence the morphology of cold atomic gas remains an open question. The nearby Riegel-Crutcher cloud, composed of long and narrow H I filaments observed in self-absorption, provides a critical test case. We present the most extensive optical polarimetric survey of this region to date, comprising more than 90,000 high signal-to-noise stellar polarization measurements combined with Gaia DR3 data. Using stellar polarization, extinction estimates, and archival Na I absorption data, we locate the cloud at a distance of $150 \pm 15$ pc, consistent with that of the Pipe Nebula. The plane-of-sky magnetic field traced by optical starlight polarization closely matches that inferred independently from Planck 353 GHz dust-emission polarization, revealing a coherent large-scale magnetic field across the region. A Rolling Hough Transform analysis shows that the H I filaments are tightly aligned with this field orientation. Together, these results provide strong observational evidence that the structure of the cold neutral medium in the Riegel-Crutcher cloud is closely linked to a highly ordered magnetic field. This level of coherence supports a scenario in which magnetic fields play a dynamically important role in shaping the cloud structure, and suggests that the Riegel-Crutcher cloud is part of a larger magnetized complex influencing gas flows in the solar neighborhood.

On 2023 February 27 at 18:15:55.77 UT, a bright fireball streaked across the sky above northern Sweden. The event offered a valuable opportunity to study the phenomenon using an optical system primarily designed for auroral studies, the Auroral Large Imaging System (ALIS_4D), that captured the event. In this study we show the capability of ALIS_4D to perform observations in support of meteor event analysis. We estimated the trajectory from the recorded data and computed the orbit. In addition, we investigated the origin of the meteoroid searching for its parent body. Fitting the analytical ablation model known as $\alpha$-$\beta$ to the trajectory as well as incorporating local wind-field data in Monte-Carlo dark-flight simulations, strewn-fields were computed and physical properties of the meteoroid were estimated. Trajectory analyses delineate a strewn field along the border between Kiruna and Gällivare in northern Sweden. Our findings indicate that the meteoroid's parent body was likely an Apollo family object. We performed an orbital similarity analysis to identify candidate parent bodies of the fireball. Our simulations suggest that close approaches with Earth could have disrupted the meteoroid's orbit, placing it on a collision course.

HI intensity mapping is a promising technique to probe large-scale structure, traditionally analyzed via two-point statistics such as the angular power spectrum. This technique has proven very powerful but may miss key non-Gaussian information present in the signal. We extend the starlet l1-norm, a multi-scale higher-order statistic previously applied to weak lensing maps, to the brightness temperature fluctuations of the HI density field. The HI signal is highly non-Gaussian at late times (z < 1) due to nonlinear structure growth, motivating the use of advanced summary statistics. We simulated full-sky HI lognormal brightness temperature maps using CAMB and GLASS, generating 10,000 realizations with associated cosmological parameters. We extracted both the starlet l1-norm and angular power spectrum from these maps. Using the JaxILI framework, we performed neural density estimation for implicit likelihood inference. The analysis considered simulated maps incorporating realistic noise and telescope beam, capturing the impact of observational effects on parameter inference. In this work, we focus on the redshift range 0.4 < z < 0.45, chosen to match the interval already targeted by existing MeerKLASS observations. The starlet l1-norm significantly outperforms the angular power spectrum in constraining cosmological parameters, achieving almost a 3x improvement in the figure of merit relative to the angular power spectrum by capturing non-Gaussian features missed by two-point statistics. Moreover, our results suggest that the starlet l1-norm is robust to several of the systematic effects included in our simulations. Our findings highlight the potential of multi-scale higher-order statistics such as the starlet l1-norm to enhance cosmological inference from future HI intensity mapping surveys.

We present new constraints on gravitational vector perturbations ($\mathcal{V}$-modes) using Cosmic Microwave Background (CMB) data, including temperature and $E$-mode polarization from SPT-3G D1, ACT-DR6, and $Planck$, as well as $B$-mode data from BICEP/Keck and SPTpol, which provide the strongest constraints on $\mathcal{V}$-modes. We consider three initial conditions (ICs) that source $\mathcal{V}$-modes: neutrino isocurvature (ISO), neutrino octupole (OCT), and a sourced mode (SMD) generated by an anisotropic stress before matter-radiation equality. We also consider including tensor modes along with $\mathcal{V}$-modes for each of these ICs. Combining all datasets, we obtain 95\% confidence level upper limits of $r_\mathrm{v} < 1.3\times10^{-4}$ (ISO), $r_\mathrm{v} < 6.8$ (OCT), and $r_\mathrm{v} < 4.2$ (SMD), with slightly tighter bounds when tensors are included, at a pivot scale $k_p\ =\ 0.05$ Mpc$^{-1}$. Interestingly, for SMD without tensors, using SPTpol $B$-modes alone yields $r_\mathrm{v} = 4.7 \pm 2.1$, consistent with zero at $2.2\sigma$. Similar result is found for SMD when including tensor perturbations. No statistically significant deviation from $\Lambda$CDM is found. However, $\mathcal{V}$-modes are not fully excluded by current $B$-mode data and should be considered when interpreting primordial signals.

Over the past couple of decades, several networks for the automatic detection of fireballs have been deployed. Their primary scientific goal is to facilitate the rapid recovery of meteorites, determine their pre-atmospheric orbits, and look for possible dynamic links with parent bodies. The Bocosur network is a contribution to the global deployment of automated fireball networks and to the increase of the number of recoverable meteorite falls. It is located in Uruguay, South America (Lat: -30$^{\circ}$ to -35$^{\circ}$). Its main scientific goal is the detection of fireballs of asteroidal origin, massive enough to produce meteorites, and also to inspire secondary-level students and teachers through their involvement in this citizen-science oriented project. The deployment of this network started in 2019, and was completed in March, 2023, when we installed 20 stations separated $\sim 120$ km, covering an area of $\sim 180,000$ km$^2$. During this period of time, one major technological upgrade was made when we migrated from a well-known camera to a higher-resolution, more sensitive system. We were able to build a completely autonomous system at an affordable cost that can be replicated in all the stations. A comparison between the astrometric and photometric performance of these two detection systems is reported. Also, a photometric methodology for estimating the brightness of very bright fireballs is presented and validated against the known magnitudes of Jupiter and the full Moon. We obtain mean residuals of the astrometric reduction of $\sim$5', and the discrepancy between the obtained brightness of Jupiter and the Moon average to 0.18 and 1.2 magnitudes, respectively. Results on the processing of a very bright (M$_{peak}\sim$-9.0 mag) fireball detected in four stations are also presented.

The reionisation time field treion(r) captures the entire history of cosmic reionisation by mapping the moment where each region of the Universe became ionised. Previous work has shown that treion(r) can be inferred from 21-cm observations, using convolutional neural networks (CNNs). However, these CNN predictors are trained on specific reionisation models, raising critical concerns about their reliability when applied to observational data potentially differing from their training assumptions. This paper aims to propose and test a method to evaluate the coherence of our CNN predictors with respect to their input model, thereby enabling the validation or exclusion of underlying reionisation models based on their reconstruction behaviour. By setting the CDM model as reference input, we evaluate the coherence of treion(r) reconstructions by comparing them across different redshifts for several prediction models as the statistics of treion (r) reconstructions should be the same for every redshift of the input maps. Our study particularly investigates CNNs trained on cold and warm dark matter (WDM) models, with WDM particle masses of 2, 3, 5, and 7 keV. We find that the predictors trained on 5 and 7 keV WDM models exhibit high-level self-consistency similar to the CDM predictor, while the 2 keV predictor, and to a lesser extent the 3 keV predictor, display significant deviations across several metrics. These findings seem to demonstrate that CNN predictors retain sensitivity to differences in the underlying reionisation model and can be used to assess model compatibility with observations. Our results highlight the necessity of validating machine-learning predictors against their input models before applying them to real data. The method proposed here offers a pathway to more trustworthy applications of CNNs in the study of reionisation.

We present James Clerk Maxwell Telescope (JCMT) observations of the $^{12}$CO(J = 2-1) emission of 42 low-mass, star-forming disk galaxies of morphological type Scd or later from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The sample, which probes the M33-like stellar-mass regime, is complemented with metallicities, star formation rates, and \hi\ masses used to investigate the star formation process and to test scaling relations involving molecular gas mass in low-mass systems. We detect CO emission in 55% of the sample and derive H$_2$ masses using both a constant Galactic and a metallicity-dependent CO-to-H$_2$ conversion factor. The 12 $\mu$m luminosity, which includes polycyclic aromatic hydrocarbon features, exhibits a tight linear correlation with the CO line emission, making it a robust tracer of global molecular gas content. The molecular gas mass - star formation rate relation, i.e. the Kennicutt-Schmidt law, is the most fundamental one and it is found to remain linear down to low stellar masses. We also find that the mean molecular gas depletion time is slightly shorter in low-mass late-type galaxies than in more massive systems, consistent with their higher specific star formation rates. Finally, while the specific molecular gas mass ($M_{\rm H_2}/M_*$) shows no significant dependence on stellar mass and a large intrinsic scatter, the HI-to-stellar mass ratio ($M_{\rm HI}/M_*$) decreases with increasing stellar mass and molecular fraction ($M_{\rm H_2}/M_{\rm gas}$), highlighting the progressive transition from atomic- to molecular-dominated interstellar medium along the galaxy population.

Accurate estimates of the absorption of X-rays by interstellar gas and dust are of crucial importance for the analysis and interpretation of almost all astronomical soft X-ray observations. However, the present X-ray absorption data extensively used by the community were derived from a reduced interstellar abundance (~70% of solar) and ignoring dust scattering. Therefore, these X-ray absorption data, although highly popular, could have been substantially underestimated. Here we update the interstellar X-ray absorption and scattering by making use of updated atomic cross sections, updated interstellar abundances, and realistic X-ray dust physics, and appropriately distributing metal elements in gas and dust. The resulting X-ray absorption and scattering data are publicly available on GitHub.

Active galactic nuclei (AGNs) exhibit complex variability across multiple wavelengths, reflecting diverse physical processes near their central engines. This work investigates the temporal variability of four AGNs Mrk~509, NGC~5548, NGC~4151, and NGC~4593 using multifractal detrended moving average (MFDMA) analysis and Fisher-Shannon information plane applied to their X-ray, ultraviolet, and optical light curves. These methods quantify the scaling behavior and complexity of the variability, revealing persistent correlations and distinct variability patterns across energy bands. The Fisher-Shannon analysis further characterizes the degree of stochasticity and structural complexity in the emission processes. Our findings support the interpretation that multifractal and information-theoretic measures provide effective diagnostics of the physical mechanisms driving AGN variability. This study demonstrates the utility of advanced time series techniques as effective diagnostics of AGN variability mechanisms.

We present a machine-learning model for generating super-resolution $N$-body simulations with non-vanishing spatial curvature, conditioned on a given low-resolution field, $\Omega_k$, $\Omega_\mathrm{m}$, $\sigma_8$, $h$, and redshift. By upscaling the resolution of $N$-body simulations, such models can drastically reduce the computational cost of producing high-resolution simulations suitable for modelling current and future surveys of large-scale structure. Our model is trained as a generative adversarial network, allowing injected noise to be interpreted as stochastic structure and enabling the generation of an ensemble of plausible high-resolution realisations. We evaluate the model performance by comparing key cosmological summary statistics in the generated simulations to their high-resolution counterparts. We find that the model accurately reproduces large-scale statistics, robustly recovering most of the power that was missing from the low-resolution input, but exhibits a residual suppression of power on small scales of up to $\sim 10\%$ at $k \sim 1\,h\,\mathrm{Mpc}^{-1}$. The abundance of halos around $10^{14}\,M_\odot$ is affected at a similar level, and we find that the profiles of these halos have a lower central density. Although the overall performance is decent, we anticipate that the fidelity of the generative model can be further increased with more and better training data, as well as through improvements in the model architecture and training process. To show a production-scale use case, we apply our model to upscale the resolution of a light cone from a large-volume $N$-body simulation with spatial curvature, producing a first-of-its-kind catalogue that simultaneously captures geometric effects at large scales and accurate nonlinear structure at small scales.

The search for supermassive black hole binaries (SMBHBs) has, in recent years, seen the dawn of exploration with several hundred candidates claimed from photometric and spectroscopic surveys monitoring AGNs. While only a handful persist to date, the advent of upcoming high-precision wide-field photometric missions motivates continuing the pursuit of confirming SMBHBs in the optical. We explore the possibility of using the ESA Plato space mission to detect photometric signatures of SMBHBs. Motivated by the Kepler observation of Spikey, the best known self-lensing flare (SLF) candidate to date, this work aims to benchmark the scientific outcome if Plato were to observe Spikey-like objects via its Guest Observer programme. Starting from the Gaia database, we assemble a catalogue of 12,226 bright ($G < 19$) high-probability Quasars for the two pointing fields of Plato's nominal mission. This Plato Quasar catalogue will be pivotal for future follow-up observations of larger photometric searches such as the Vera Rubin LSST survey. We use the Plato camera simulator, PlatoSim, to realistically explore the noise budget in Plato's faint limit, while generating mock light curves to benchmark Plato's ability to recover signatures of SMBHBs. We show that, although not at all designed for the purpose, Plato is capable of detecting Spikey-like SMBHB candidates through their relativistic photometric signatures using Bayesian inference and evidence. Plato will in particular be able to confirm or rule out Spikey and Spikey-like objects with a limiting magnitude of $G\leq18$. With a minimum 2-yr baseline per pointing field, we show that Plato not only could play an essential role in future SMBHB research, but may be an integrated part of the observational fleet of continuous high-precision facilities monitoring SMBHB candidates in the near future.

Electromagnetic (EM) sounding can constrain the electrical structure of Enceladus and, in turn, the salinity of its ocean and the porosity, fluid content, and thermal state of its hydrothermally active core. Here, we assess the feasibility of EM sounding at Enceladus using both global (orbiter) and local (lander) EM induction transfer functions. We provide a physical framework for modeling EM induction for 1-D and 3-D subsurface conductivity models and discuss how transfer functions can be estimated from global or local measurements of the magnetic and electric fields. We simulate 3-D induction effects arising from variations in ice-shell thickness. The magnitude of these effects in the magnetic field correlates with the ice-shell thickness at the surface and is strongly dependent on the ocean's conductivity. These magnetic variations, if observed, would favor a moderately to highly conductive ocean, providing lower bounds on salinity and volatile content. The absence of these effects indicates a thicker, more homogeneous ice shell and/or a lower-conductivity ocean. Given plausible magnitudes, a polar-orbiting mission with low-altitude measurements will be required to detect these effects. In summary, an orbiter will constrain global ocean conductivity using long-period induction and possibly map the ice thickness variations. The detailed EM sounding of both the hydrosphere and the core can be achieved by a lander-based broadband EM sounding at periods $\approx 10^1-10^5$ s to probe ocean salinity and thickness, as well as core properties including porosity, fluid content, and temperature.

We present new JWST/MIRI-MRS data of the z=4.601 extremely luminous obscured quasar WISEA J224607.56-052634.9 (W2246-0526). Our fits of its spectral energy distribution (SED) with the SED fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models) predict an active galactic nucleus (AGN) fraction in the range 72-81 per cent, an intrinsic AGN luminosity of 4.2-7.2 x 10^14 Lo, a polar dust luminosity of 1.6-1.7 x 10^14 Lo, a black hole mass of 1.3-2.3 x 10^10 Mo (assuming the quasar is accreting at the Eddington limit), a star formation rate (SFR) of 360-2900 Mo/yr and a stellar mass of 4.8-5 x 10^11 Mo. The stellar and black hole masses of W2246-0526 are typical of a giant elliptical galaxy at z=0. We find statistically significant evidence for the presence of a hot dust component, which we interpret as polar dust in the context of a torus geometry, based on recent results obtained for nearby AGN. We explore two smooth and two two-phase models for the AGN torus, to put constraints on the AGN fraction of the galaxy, the black hole mass and its SFR. We show that the presence of polar dust affects the estimate of the AGN luminosity and we recommend to take into account this component in SED fits of other high-redshift obscured AGN/quasars. Despite the large difference in luminosity, we discuss possible links between the presence of this hot dust component in W2246-0526 and in some local AGN, suggesting that they may have a different origin.

Obtaining robust galaxy number counts is crucial for understanding galaxy evolution, and submillimetre counts in particular have proven valuable for revising subgrid physics models in cosmological simulations. In confusion-limited surveys, which are common at these wavelengths, statistical methods such as $P(D)$ fluctuation analysis are required to recover counts of faint, unresolved galaxies. However, the standard $P(D)$ framework assumes that galaxies are Poisson-distributed, whereas in reality galaxies are clustered. Using simulations, we demonstrate that this clustering systematically biases $P(D)$-derived number counts, and present an empirical method that simultaneously measures and corrects for this bias by combining the 1- and 2-point statistics in the map, thereby maximising the information extracted from the data. Applying this method to deep Herschel-SPIRE observations of the GOODS-N field, we provide revised galaxy number counts at 250, 350 and 500$\mu$m. Our results indicate that at 500$\mu$m clustering inflates the apparent counts by a factor of 1.6 around 10mJy and slightly suppresses the faintest sub-mJy counts, with milder effects at 350$\mu$m and 250$\mu$m owing to the smaller beam sizes. This methodology is broadly applicable to other confusion-limited data sets with well-characterised beam and noise properties, including SCUBA-2 and CCAT, enabling unbiased exploitation of the full statistical information in current and future far-infrared and submillimetre surveys.

Radiative cooling can strongly influence the structure and dynamics of black hole accretion disks. Here, we perform general relativistic magnetohydrodynamic (GR-MHD) simulations of magnetically arrested disks (MADs) around a non-spinning black hole. Radiative cooling is consistently included in the simulations and its intensity is scaled by the mass accretion rate ranging from $10^{-7}$ to $10^{-4} \dot{M}_{\mathrm{Edd}}$. Considering synchrotron and bremsstrahlung emission, we quantify how radiative losses modify the disk structure and the accretion dynamics. In the inner MAD disk regions, accumulation of magnetic field regulates gas accretion, enforcing the gas into a discrete interchange-driven filamentary structure. We identify, both analytically and numerically, a transition mass accretion rate above which radiative cooling becomes faster than the heating, which is assumed to occur via local coupling to the magnetic field. Above this mass accretion rate, cooling substantially reduces the gas thermal pressure, leading to considerably thinner and denser accretion filaments, and a substantial increase in radiative efficiency, relative to lower accretion rates. We show that under these conditions, conventional measures of the disk scale height become misleading in MAD flows. We therefore introduce an alternative definition based on the polar position of the density maximum, which more robustly characterizes the filamentary structure of the disks in the presence of strong magnetic fields and cooling.

Nova-like variables are high-accretion-rate cataclysmic variables (CVs) that, in contrast to dwarf novae, do not undergo outbursts caused by thermal-viscous instability. However, a small group of nova-likes, classified as AH Pic-type stars, show recurrent small-amplitude outbursts, which are unexpected by the classical disk instability model. The physical mechanisms underlying these outbursts are not clear. In this study, we present a comprehensive time-domain analysis of the CV candidate GSC 08227-00723 using photometric data from ASAS-SN and TESS. The long-term light curves reveal a sequence of low-amplitude, recurrent stunted outbursts with recurrence times ranging from 30 to 50 days. Notably, these outbursts are frequently preceded by precursor brightenings, a feature reminiscent of super-outbursts in SU UMa stars driven by tidal instability. Period analysis of high-cadence TESS data identifies a coherent periodic modulation at 0.297 d, likely the orbital period, and a persistent positive superhump signal at 0.352 d. The latter corresponds to an exceptionally large superhump excess of $\epsilon^+ \approx 0.19$, surpassing typical values seen in CVs. Additionally, we detect short-timescale variability resembling quasi-periodic oscillations in the TESS light curves. Based on the outburst properties and photometric behavior, we classify GSC 08227-00723 as a new member of the AH Pic subclass of nova-like stars. We discuss how tidal effects may be involved in the observed behavior, although the exact mechanism is still unclear.

The formation environments of merging binary black holes remain uncertain. While hierarchical assembly in dense stellar clusters has been widely explored as an explanation for black holes exceeding the stellar-mass limit, growth through gas accretion in active galactic nucleus (AGN) disks is an alternative that has received less observational scrutiny. Here we search for an accretion-origin subpopulation using only spin magnitudes, fitting a three-component mixture model to 166 binary black hole mergers from LIGO--Virgo--KAGRA with component shapes fixed from theoretical predictions and only the mixing fractions inferred from the data. We find strong evidence ($ln B = 5.7$) that $\sim 10\%$ (90% credible interval $[1\%, 14\%]$) of detected mergers belong to a subpopulation with primary spins clustered near $a_1 \approx 0.9$, consistent with the theoretical prediction for accretion spin-up. The hierarchical-merger prediction of $a_1 \approx 0.7$ is decisively disfavored as the location of the high-spin subpopulation ($ln B = 5.7$). Post hoc validation reveals that the accretion candidates have systematically higher masses (median $m_1 = 58\,M_\odot$) and aligned spins (median $\chi_{\rm eff} = 0.33$, vs. $0.04$ for standard-dominated events). The accretion subpopulation is not limited to systems above the pair-instability mass gap: GW190517 ($m_1 \approx 39 M_\odot$) is among the top candidates, demonstrating that accretion spin-up operates across a range of masses. GW190521, previously interpreted as a hierarchical merger, shows comparable support for an accretion origin. These results provide the first population-level observational evidence for an accretion-origin subpopulation in black hole mergers.

We quantify systematic biases in rotation curves inferred from steady, axisymmetric Jeans modelling when the underlying stellar velocity field is non-axisymmetric. Using a high-resolution $N$-body/hydrodynamic simulation of an isolated Milky Way-like disk galaxy, we construct mock stellar-kinematic measurements for two observer azimuths relative to the bar. One observer is placed at a Solar-like viewing angle of $25^\circ$ from the bar major axis, and the other at $115^\circ$. For each configuration, we analyse multiple snapshots and compare the Jeans-inferred circular-speed curve, $V_{\rm c,Jeans}(R)$, with a reference axisymmetric circular-speed curve, $V_{\rm c,axi}(R)$, defined from the azimuthally averaged ($m=0$) component of the gravitational field. The Jeans analysis is performed in a wedge-shaped mock observational volume that mimics limited sky coverage. For the $25^\circ$ configuration, the mean azimuthal streaming is typically higher than the azimuthally averaged expectation by $\approx 10$--$15~\mathrm{km\,s^{-1}}$, which leads to an average overestimate of the axisymmetrically defined circular speed by $\approx 4\%$ ($\approx 10~\mathrm{km\,s^{-1}}$) in the inner disk. Across snapshots, the mean offset corresponds to a $\sim 1.5$--$2\sigma$ systematic deviation of $V_{\rm c,Jeans}$ from $V_{\rm c,axi}$. For the $115^\circ$ configuration, the bias reverses sign and $V_{\rm c,Jeans}$ tends to underestimate $V_{\rm c,axi}$. Under the usual spherical approximation, a $\approx 4\%$ bias in $V_{\rm c}$ corresponds to an $\approx 8\%$ bias in the enclosed dynamical mass at fixed radius. These results imply that steady, axisymmetric Jeans modelling of Milky Way stellar kinematics can overestimate the axisymmetrically defined circular-speed curve at the percent level unless non-axisymmetric streaming is modelled explicitly or included in the error budget.

Hypercompact HII regions (HC) are regions of ionized gas associated with the early stages of high-mass star formation. With the aim of better understanding their characteristics, we studied five candidate HC HII regions. Here, we present observations with the Jansky Very Large Array (VLA) at 2 and 6 cm, with angular resolutions in the range of $\sim$1 -- 3\arcsec and report the images of the detected sources and the measured parameters. In addition, we explore several possible scenarios, considering the regions as both uniform and non-uniform spheres, and as winds, both spherical and collimated. In most cases, the sources were unresolved, but by applying the models, we estimate that their sizes vary in a range of 0.3 to 3.7 mpc while their electron densities are in the range of $1.3 \times 10^{5}$ to $2.4 \times 10^{6}$ cm$^{-3}$, indicating that most sources are consistent with small, weak UC HII regions, although a few remain viable candidates for HC HII regions, with G40.28$-$0.22 as the strongest case. We do not rule out the possibility that some sources are jets or stellar winds.

We present a three-point inverse solution for reconstructing meteoroid deceleration and mass-loss histories from sparse observations constrained only by the entry, peak-brightness, and terminal points. The method combines the $\alpha$-$\beta$ analytical formalism with a derivative-free global optimizer and a numerical inversion of the height-velocity relation, enabling the retrieval of physically consistent solutions even when full velocity profiles are unavailable. Applied to the 2017-2018 European Fireball Network (EN) catalog, the approach achieves an 88% convergence rate when fitting only height-velocity pairs, and 63% when terminal and initial masses are also imposed. 52% of mass-constrained solutions (34% overall) yield bulk densities consistent with their $PE$ classes, with higher strength emerging as the primary discriminator among events retaining coherent classifications when only 3 points are used as input data. Rapidly evolving high-energy, high-mass events show the largest incompatibility with the $\alpha$-$\beta$ model. The inversion produces a continuous bulk-density distribution spanning $\sim$300-4000 kg$\,$m$^{-3}$, in contrast to the discrete densities fixed by $PE$-based categories. The EN fireball dataset is now supplemented with self-consistent $\alpha$ and $\beta$ estimates.

In this Paper we consider twisted accretion disks in supermassive binary black hole by analytical and numerical means. It is assumed that the disk orbiting around the more massive rotating component and that the disk rings are inclined with respect to the orbital plane. We use orbital parameters of the binary often employed in the precessing massive (PM) model of the well-known blazar OJ 287. Unlike our previous investigation of a similar problem, here we consider disks with both small and relatively large relative thicknesses $\delta=h/r$, where $h$ is the disk's height at a typical radius $r$, as well as a range of values of the viscosity parameter, $\alpha$, including the cases when $\alpha \lesssim \delta$. Similar to our previous results, we find that the twisted disk relaxes to a quasi-stationary state in the frame precessing with the Lense-Thirring frequency of the orbit. However, its shape is qualitatively different from that corresponding to the case of $\delta=10^{-3}$ and $\alpha=0.1$ considered in our previous work. In a disk with $\delta=10^{-3}$ but $\alpha \le 2\cdot 10^{-2}$, we find the new effect of generation of a twisting spiral wave near the resonance between a forcing frequency associated with the presence of the secondary and the Lense-Thirring frequency of a particular disk ring defined in the precessing frame. We propose an analytic theory of it, which is in a good agreement with our numerical results. This effect leads to multiple crossings of the orbit with the disk per one orbital period, which contradicts the PM model. When $\delta \gtrsim 0.1$, a typical disk's inclination within the orbit of the binary is smaller than that of the orbit which results in only two crossings of the orbit with the disk per one orbital period. We suggest that the additional heating of the disk gas by the secondary-disk collisions may result in $\delta \sim 0.1$.

Assuming that the compact source at the Galactic center, Sagittarius A*, is a massive boson star, we fit the near-infrared flare astrometry data. We consider 12 discrete boson star configurations and model the flare as a hotspot on a circular equatorial orbit. The analysis is performed in a Bayesian framework using nested sampling, yielding the marginal posterior distributions of all parameters as well as the Bayesian evidence for each model. For comparison, the same procedure is applied to a Schwarzschild black hole. The resulting Bayesian evidence values differ only marginally between the boson star and black hole cases, and the well-determined mass of Sgr~A* (${\sim}4.296\times 10^6\,M_\odot$) falls within the 68\% highest density interval in every configuration. We conclude that, under current near-infrared astrometric constraints and within the considered parameter ranges, a massive boson star and a Schwarzschild black hole remain statistically indistinguishable as the compact object at the Galactic center.

Spectropolarimetric interpretation of chromospheric lines requires solving the radiative transfer problem under non-local thermodynamic equilibrium (non-LTE) conditions. This means computing atomic-level populations self-consistently with the radiation field. While traditional inversion codes employ 1.5D approximations, they neglect horizontal radiative transfer, which can be significant near magnetic structures and in the chromosphere. We present a method to solve 3D atomic-level populations using Graph Neural Networks (GNNs), extending prior 1.5D work to the full 3D domain. By discretizing the solar atmosphere as a directed graph, in which nodes encode physical properties and edges encode spatial distances, an Encode-Process-Decode GNN propagates information to efficiently capture radiative coupling. The network is trained on a Bifrost simulation using Ca II populations from Multi3D as ground truth. The trained GNN accurately predicts populations of the five-level Ca II atom plus continuum. Correlations exceed 0.99 in the photosphere and mid-chromosphere; errors in the upper chromosphere remain unbiased. Inference is $\sim 10^6$ times faster than traditional iterative solvers. Spectral synthesis of the Ca II 8542 Å line yields intensity profiles with $< 2 \%$ mean residuals relative to the full 3D solution. This framework bypasses the computational bottleneck of iterative solvers while preserving essential non-LTE physics, including horizontal transfer, paving the way toward routine 3D non-LTE inversions.

In the last year, several pieces of evidence have pointed to a possible deviation from the standard cosmological model, $\Lambda$CDM. The recent work by the Dark Energy Survey (DES) collaboration reports a preference in the ballpark of $3\sigma$ in favor of dynamical dark energy against the standard cosmological model. For that, it used its final analyses of Baryonic Acoustic Oscillations and type Ia Supernovae, both sensitive to the expansion history of the Universe, in combination with the Cosmic Microwave Background (CMB) from Planck. This adds to the growing debate about the nature of dark energy. \textit{Published as a Perspective in Nature Astronomy in August 2025}.

The origin of the highly ionized iron emission (Fe XXV at $6.7\,\mathrm{keV}$) characterizing the Galactic ridge X-ray emission (GRXE) remains a fundamental puzzle in high-energy astrophysics. Although the GRXE continuum is largely resolved into discrete populations of cataclysmic variables and coronally active stars, these sources exhibit Fe XXV equivalent widths significantly lower than that of the total GRXE, leaving the intense iron line emission unexplained. In this work, we cross-correlated the XMM-Newton survey of the inner Galactic disk with Gaia DR3 astrometry to identify and characterize hard X-ray sources ($>2\,\mathrm{keV}$) with reliable stellar counterparts. We selected 107 X-ray sources located within the red giant branch of the color-magnitude diagram, many of which are verified long-period variables. These sources exhibit high X-ray luminosities ($L_{\mathrm{X}} \approx 10^{31}$--$10^{33}\,\mathrm{erg~s^{-1}}$), significantly exceeding the typical coronal saturation levels of single giants. Their X-ray spectra are notably harder than those of quiescent stellar coronae, with plasma temperatures reaching up to $kT \approx 6\,\mathrm{keV}$ and a prominent emission feature at $\sim 6.7\,\mathrm{keV}$. The combination of high $L_{\mathrm{X}}$, hard spectra, and intense Fe XXV emission identifies this population as accretion-powered binaries associated with late-type giants. Our analysis demonstrates that this population contributes $\sim 20\%$ of the total GRXE continuum and $\sim 40\%$ of its iron line emission, providing a key component to resolving the Galactic X-ray background puzzle.

We study the postinflationary dynamics of an Einstein-Cartan-Holst gravity-motivated inflationary scenario, known as Einstein-Cartan pseudoscalaron inflation, coupled to a type-I seesaw extension of the Standard Model with three heavy right-handed Majorana neutrinos. In particular, we show that nonthermal leptogenesis emerges as a necessary and self-consistent mechanism for generating the observed baryon asymmetry of the Universe, mainly because of the universal coupling of the inflaton to the additional heavy Majorana fermions. The resulting framework provides theoretical predictions that are fully compatible with the latest cosmological constraints from the Cosmic Microwave Background, Baryon Acoustic Oscillations, and Big Bang Nucleosynthesis, as well as with neutrino oscillation experiments, for a wide range of the fundamental Barbero-Immirzi model parameter $\gamma$, which controls the inflationary and postinflationary phases. In particular, for $\gamma \sim -1/100$ and a lightest Majorana-neutrino mass of order $10^{13}$ GeV, we find a scalar spectral index $n_s \sim 0.970$, a tensor-to-scalar ratio $r \sim 0.004$, for a number of e-folds before the end of inflation $N_e \lesssim 60$, and a baryon-to-entropy ratio $n_B/s \sim 8.7 \times 10^{-11}$.

We investigate the optical appearance of a Kerr-Bertotti-Robinson (Kerr-BR) black hole illuminated by a geometrically and optically thin accretion disk. Instead of using a phenomenological power-law emissivity, we adopt a magnetically driven synchrotron emissivity proxy coupled to the local electromagnetic environment. With a backward ray-tracing framework, we examine the effects of the spin $a$, magnetic parameter $B$, and observer inclination $\theta_O$ on the ray-classification maps, redshift distributions, and specific-intensity images. We show that the ISCO position is modified by both $a$ and $B$, and that rapidly rotating prograde configurations can develop an additional model-dependent inner cutoff when the magnetically dominated approximation underlying the emissivity prescription ceases to be applicable. High-resolution one-dimensional intensity profiles further separate the direct image, the $n=1$ lensing-ring contribution, and the higher-order $n\geq 2$ photon-ring subimages, while quantifying the Doppler-induced brightness asymmetry. Retrograde disks exhibit a wider emission-depleted central region because of the outwardly shifted ISCO, making the higher-order lensed components more clearly distinguishable from the direct emission. These results show that the disk inner boundary and the magnetic-field-dependent emissivity can substantially influence the observable appearance of Kerr-BR black holes.

Eclipse shadow bands are transient intensity patterns that can appear on the ground near solar eclipse totality. This study presents a reusable image-processing framework for analyzing shadow-band video recordings collected with consumer-grade cameras. The framework quantifies band orientation, band prominence, and band power spectral density from video recordings. Applied to two eclipse datasets, the method detected statistically significant shadow-band activity during eclipse windows that align with the scintillation theory for shadow bands. The results also highlight simultaneous superimposed eclipse shadow band modes with orthogonal orientations. This demonstrates that consumer grade cameras can support quantitative analysis of shadow bands and may support future observational and atmospheric studies.

Neutron star (NS) low-mass X-ray binary (LMXB) systems depict a variety of X-ray spectral and timing features, which can be useful to probe the accretion-ejection mechanism in the strong gravity regime. Here, we study the relatively unexplored and faint NS LMXB 1A 1246-588, which is also an ultra-compact X-ray binary (UCXB) with a white dwarf donor. We investigate its temporal and spectral behavior using pointed NICER and AstroSat observations, supported by long-term MAXI/GSC monitoring. The MAXI light curve shows modest, recurrent outburst-like enhancements, providing the long-term flux context for interpreting the pointed observations. During the AstroSat observations in 2017, the source exhibits an absorbed 0.4-20 keV flux of $(1.18 \pm 0.02)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$, while during the NICER observations in 2019, it spans an absorbed 0.5-10 keV flux range of $(0.7-3.7)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$ and traces an atoll-like pattern in the hardness-intensity diagram. Broadband spectral modeling shows that the emission is well described by a soft blackbody and a hard Comptonized component, with no statistically required multicolor disk contribution. The blackbody temperature increases from 0.28 to 0.39 keV, with an emitting radius consistent within 6.9-13.8 km, while the Comptonization photon index varies from 1.8 to 2.3. We find that the observed spectral-state evolution is driven by a redistribution of accretion power between thermal emission from the NS boundary layer and Comptonized emission, consistent with atoll-type behavior. These results provide the first quantitative, multi-epoch view of accretion-state evolution in 1A 1246-588, revealing systematic changes in the thermal boundary-layer emission and the Comptonizing region in this UCXB system.

We investigate the thermal evolution of isolated neutron stars containing hyperon--mixed kaon--condensed matter, focusing on the role of proton superconductivity. The equation of state utilized for cooling calculation is based upon the minimal relativistic mean--field framework supplemented by chiral SU(3) dynamics for kaon condensation with an additional component on the three-baryon force, which ensures stiffness at high densities enough to meet astrophysical constraints on neutron-star masses and radii. We show that the nucleonic direct Urca processes operate at relatively low stellar masses ($M \gtrsim 1.3\,M_\odot$), erasing any observable signature of strangeness in the absence of superfluidity. However, if the proton $^1{\rm S}_0$ superconductivity works, because of suppression of fast neutrino cooling processes, the cooling scenario could become relevant with the strangeness, depending on the density regions of the pairing gap. In particular, if the proton superconductivity is so strong in high-density regions ($T_{c,p}\sim10^{10}~{\rm K}$), the nucleon and hyperon direct Urca processes shut down, which makes the kaon-induced Urca processes dominant in massive neutron stars. This scenario is in good agreement with several cold isolated neutron stars identified recently. Hence, we suggest that strong proton superconductivity can render kaon condensation observationally visible through cold neutron-star observations, providing a potential signature of strangeness in dense matter.

We investigate the morphology and structural properties of 246 massive compact quiescent galaxies (MCGs; $\log M_{\star} \sim 10$-$11$, $\sigma_{\mathrm{e}} \sim 150$-$350\,$km\,s$^{-1}$, $R_{\mathrm{e}} \sim 0.7$-$2.5\,$kpc) at $z \sim 0$, selected as outliers in the stellar mass-velocity dispersion and velocity dispersion-size relations, using $g$-, $r$-, and $i$-band Hyper Suprime-Cam images. We compare them to a control sample of average-sized quiescent galaxies (CSGs) matched in stellar mass, star formation rate, redshift, and $g-i$ color. Both samples are dominated by S0 galaxies, comprising $93\%$ of MCGs and $71\%$ of CSGs, while ellipticals account for $4\%$ and $11\%$, respectively. The fraction of interacting or morphologically disturbed systems is low in both samples ($13\%$ for MCGs and $16\%$ for CSGs). Multi-component decompositions of the $g$- and $r$-band images show that $75\%$ of MCGs require a three-component model (bulge, disk, and envelope), while $21\%$ are best fit by two components and $4\%$ by a single Sérsic profile. Two-component MCGs are preferentially low-inclination systems, suggesting that the three-component fraction represents a lower limit. In contrast, only $7\%$ of CSGs exhibit a comparable three-component structure. Bars are present in $29\%$ of CSGs but are absent in MCGs. For three-component systems, MCGs and CSGs have similar bulge ($R_\mathrm{e}=0.39$ vs.\ $0.45$\,kpc) and envelope ($R_\mathrm{e}=6.4$ vs.\ $5.8$\,kpc) sizes, while MCG disks are significantly more compact ($R_\mathrm{e}=1.9$ vs.\ $3.3$\,kpc). The envelope component shows a broad ellipticity distribution ($\epsilon_\mathrm{Envelope} \sim 0.0$-$0.6$), which we interpret as corresponding to either a stellar halo or a thick disk.

Recent JWST observations suggest that galaxies at z > 9 may be dominated by low-opacity SNe-produced dust before efficient ISM grain growth is established. This transition in dust origin and opacity could explain both the prevalence of galaxies with extremely low dust attenuation and the excess of UV-bright galaxies relative to most pre-JWST predictions. We investigate whether this transition, combined with evolving star-formation efficiency, can reproduce these observed properties. We develop a physically motivated attenuation framework combining (i) extinction laws for reverse-shock-processed SNe dust, (ii) metallicity- and dust-to-metal-dependent opacity scalings, and (iii) porous radiative-transfer geometries allowing partial UV-photon leakage. Unlike outflow-driven scenarios requiring large-scale gas evacuation, our approach preserves gas reservoirs while reducing effective UV opacity through dust composition and geometry. We introduce extinction-based, gas-based, and hybrid attenuation prescriptions linking SNe-dominated and ISM grain-growth dust regimes. We find that the observed A_FUV-M_star relation at z > 9 is best reproduced for an intrinsic FUV dust opacity kappa_UV(dust)=1000 cm2/g, characteristic of low-opacity SNe dust, naturally producing very low attenuation even in gas-rich galaxies. This regime reproduces galaxies with extremely low dust attenuation (GELDAs), which dominate observed samples at z > 9. Applied to intrinsic UV luminosity function models, our SNe-dominated and hybrid prescriptions mainly suppress the brightest galaxies, bringing predictions into agreement with JWST measurements without requiring extreme star-formation efficiencies or dust-free interstellar media. Our results suggest that the UV-bright galaxy excess at z > 9 reflects a transition in dust origin and opacity during the earliest phases of galaxy evolution.

We simulate star formation and star cluster assembly inside a molecular cloud with parameters we derive directly from observations of the Aquila Rift. We model the evolution of stars and gas together while resolving close encounters between stars, the formation of new stars, and stellar feedback to follow cluster formation up to the expulsion of the surrounding gas. We find that star formation takes place in clumps spaced unevenly along Serpens South and that these clumps accrete surrounding gas to grow and form new stars. Gas flows along the filament promote the merger of these clumps into a star cluster inside the Serpens South filament. The imprints of these mergers are seen in the dynamics of the Serpens South cluster in the form of velocity space anisotropies, cluster rotation, and cluster expansion. Before gas is removed from the simulation, the Serpens South cluster merges with the nearby cluster W40 non-monolithically resulting in a fractal cluster at the end of the simulation. The dynamics inherited from the mergers throughout the simulation are still seen in the final bound stellar system after the gas has been removed. We compare these results with recent observations of Milky Way clusters to comment on their formation histories. We also study how our results change when lowering the mass resolution of our simulation and removing observations of dense gas tracers from our initial condition setup. Each of the three simulations result in different final cluster configurations pointing towards the importance of gas in cluster assembly.

Resonant populations of trans-Neptunian objects serve as crucial dynamical archives for unraveling the early migratory history of the Solar System. A quantitative assessment of the capture efficiency into various mean motion resonances (MMRs) during migration is essential for understanding the origins of these populations, constraining migration parameters, and reconstructing of the primordial planetesimal disk. Using numerical simulations, this study systematically investigates the capture capability of exterior MMRs during Neptune's outward migration in a planar model. For a specific p:q MMR, the small bodies can be captured only when their eccentricities surpass a certain threshold, which increases with faster migration rates, greater distances of MMRs, and higher resonance orders. On the other hand, as long as a particle's eccentricity is suitable, its capture efficiency shows little dependence on the migration rate; instead, it mainly depends on the p value and heliocentric distance, decaying exponentially as either parameter increases. Based on our simulation results, we derive for the first time a simple empirical expression to calculate eccentricity threshold and the capture efficiency. This research provides a systematic quantitative framework for understanding capture into Neptunian MMRs during migration. Future integrations of more comprehensive observational data will facilitate a more precise reconstruction of the Solar System's early dynamical evolution.

We present RAYTHEIA, a high-performance reverse ray-tracing algorithm designed to efficiently solve three-dimensional direction-dependent equations in astronomical simulations. The algorithm uses a dual-grid framework in which the native simulation mesh -- serving as the source grid for ray emission -- and an adaptive mesh refinement (AMR) Cartesian contribution grid are constructed for efficient ray-walking and contribution accumulation. The core of the algorithm integrates a leaf-only linear-octree data structure to reduce memory overhead, the digital differential analyzer (DDA) traversal method to efficiently determine the ray-walking path, Morton Code indexing to fast leaf cell lookup during traversal, and the slab method to analytically compute the path length. Furthermore, RAYTHEIA employs a hybrid (MPI/OpenMP) distributed parallel framework with a chunk-to-chunk communication strategy, achieving exceptional, near-ideal linear speed-up ratio and delivering high-end performance. We integrate RAYTHEIA with the 3D-PDR code to solve the complex chemistry and radiation transfer in photodissociation regions (PDRs). This allowed the modelling of three-dimensional PDR chemistry in a turbulent, star-forming cloud at an unprecedented resolution of $512^3$ grid cells. The algorithm demonstrates accuracy and convergence even at low angular resolutions. We further showcase the capabilities of RAYTHEIA by producing high-resolution synthetic emission maps of key diagnostic lines of a star-forming region capturing physical effects such as [O I] $63\mu$m self-absorption, measuring the [C I]-bright but CO-dark molecular gas, and deriving a CO-to-H$_2$ conversion factor in agreement with observations.

The leading-trailing asymmetry in lunar crater distribution provides a critical record of inner solar system dynamics, yet the long-standing discrepancy between the observed higher asymmetry and lower theoretical predictions indicates a gap in our understanding of the impactor population. This paper hypothesizes that lunar impact ejecta, which can enter Earth-like orbits and return, constitute a previously unaccounted-for component. Through numerical simulations, we find that ~25% of escaped ejecta will re-impact the Earth-Moon system within 3 Myr, with about 1.2% striking the Moon. Crucially, these lunar impacts exhibit an extreme leading-trailing asymmetry with a ratio of 5.9. Our results indicate that lunar ejecta, if comprising ~15% of total impactors, can fully explain the observed asymmetry, leading to their recognition as active agents shaping the lunar impact record. This work provides new constraints for understanding the impact environment of the Earth-Moon system, with direct relevance to the interpretation of lunar geology, the transport of lunar material to Earth, and ongoing space exploration missions.

Tidal fragmentation at perihelion and nontidal fragmentation elsewhere cause the orbital distribution of Kreutz sungrazers of all sizes to be extremely complicated and highly nonuniform. Among the features are (largely fortuitous) clusters of bright (naked-eye) objects and clumps of dwarf objects (often closely genetically related, as their detection primarily by the SOHO coronagraphs suggests) on the one hand; and both spectacular and less brilliant sibling sungrazers, whose perihelion times are scattered over centuries, on the other hand. Investigation of four fragment nuclei of the Great September Comet of 1882, the products of a perihelion breakup of the comet's original nucleus, showed that their orbital periods followed a distinct pattern, which likewise applied to other tidally split sungrazers and was characterized by a specific value of the second difference of parameter u_frg of neighboring fragments' centers of mass. The algorithm has a potential for the prognostication of bright Kreutz sungrazers over the rest of the 21st century and beyond. However, because of its as yet unverified empirical character, the utmost caution should be exercised when applying the procedure.

Modified Newtonian Dynamics (MOND) is a paradigm that can do away with dark matter at galaxy scales, but displays a residual missing mass discrepancy in galaxy clusters. Prompted by the updated JWST-based gravitational lens model of the Bullet Cluster, I confirm here that this cluster exhibits the same residual missing mass discrepancy as other clusters of similar mass in the MOND context. Moreover, this missing mass should be mostly collisionless, since it is centred on the galaxies of the Bullet Cluster.

Wind is a key mechanism for supermassive black hole (SMBH) feedback to their host galaxies. In tidal disruption events (TDEs), black holes spend most of their time accreting at highly sub-Eddington rates, implying that feedback from persistent sub-Eddington winds could be significant. We investigate the effects of black hole mass, viscosity parameter and stellar debris temperature on the properties of winds from hot accretion flows in TDEs. We find that more massive black holes yield a higher accreted fraction and launch faster winds, while the debris temperature has a negligible influence on the accretion flow. For $\alpha=0.1$, the mildly-relativistic unbound winds ($\sim 0.1c$) are launched predominantly from the outside of the accretion flows along the equatorial plane, with a kinetic energy of $\sim10^{-4}L_\mathrm{Edd}$. In contrast, convective bound outflows dominate for $\alpha=0.01$, which differs from the true winds typically seen in active galactic nuclei and X-ray binaries. Potential applications for explaining delayed radio brightening in TDEs at $\sim10^3$ days and for searching for intermediate-mass black holes through radio and X-ray surveys are also discussed.

Aims. We constrain the dynamical perturbations of the spiral arms in the Milky Way disk, based on the non-axisymmetric streaming motions of RGB stars revealed by APOGEE and \textit{Gaia}. Methods. We develop a revised steady-state radial-velocity response model that incorporates both the \(V_{R,\sin}\) and the dynamically important \(V_{R,\cos}\) components for a two-armed logarithmic spiral potential. The model is validated using orbit integrations with \texttt{AGAMA} and Bayesian parameter recovery with \texttt{dynesty}, and is applied to the smoothed two-dimensional radial-velocity field of RGB stars while accounting for Lindblad and corotation resonances. Results. The revised model reproduces the phase and amplitude of the mock radial-velocity field to the \(\sim2\%\) level, substantially improving upon earlier \(V_{R,\sin}\)-only formulations. Applied to the observational data, it yields a robust pitch angle of \(p \simeq 10^\circ\) and a local surface density contrast of \(\xi \simeq 5\)--\(18\%\) at the solar radius. The radial scale length is less well-constrained (\(h_{R,1} \simeq 40\)--\(50\,\mathrm{kpc}\)) due to intrinsic parameter covariance. Resonance effects strongly shape the velocity field, thus affecting the fitting: the radial velocity becomes extremely large near the Lindblad resonances, whereas it vanishes close to the corotation resonance. Conclusions. Our results demonstrate that including both the \(V_{R,\sin}\) and \(V_{R,\cos}\) terms is essential for a physically consistent interpretation of stellar streaming motions induced by a spiral potential. The observed kinematics constrain the spiral pattern speed to \(\Omega_{p} \approx 10\)--\(20\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}\).

We present a comprehensive analysis of Vela X-1 using two new NuSTAR observations, placed in the context of four earlier datasets obtained between 2012 and 2020. The energy-resolved pulse profiles demonstrate a significant transformation from an asymmetric structure at low energies to distinct double peaks above 12 keV, whereas the pulse fraction escalates with photon energy but decreases with flux. Broadband spectra validate the Fe K alpha emission line and disclose both fundamental and harmonic cyclotron resonant scattering characteristics (CRSF). We observe no substantial link between CRSF energies and luminosity, contrary to previous findings; rather, the photon index and folding energy demonstrate distinct anti-correlations with flux, aligning with sub-critical accretion and increased Comptonization in the accretion column. Our results provide the first clear evidence that the harmonic CRSF in Vela X-1 does not follow the long-term decay previously claimed. The fundamental line energy also displays an irregular evolution, without a clear monotonic trend. Notably, the harmonic-to-fundamental energy ratio departs from the canonical value of two, suggesting that the line-forming regions are located at different heights within the accretion column. These results provide new constraints on the accretion geometry and magnetic field topology of Vela X-1, highlighting the importance of continued monitoring with current and future X-ray observatories.

This chapter reviews how machine learning (ML) can be used to extract astrophysical and cosmological information from redshifted 21 cm observations of the cosmic dawn and the Epoch of Reionization, with an emphasis on SKA-Low science. We first summarize the basic physics of the global signal and spatial fluctuations, highlighting why the signal is intrinsically non-Gaussian and highly sensitive to poorly constrained properties of early galaxies and radiation backgrounds. We then discuss the main analysis bottlenecks that dominate current and future observations: bright foreground contamination, radio-frequency interference, ionospheric distortions, calibration errors, and the computational burden of repeated forward modeling in high-dimensional parameter spaces. Building on this context, we organize the ML literature by its role in the pipeline: observation-domain methods that operate on contaminated measurements and image products, theory-domain methods that accelerate or compress forward modeling, and inference-domain methods that map complex observables to astrophysical and cosmological constraints. The central message is that ML is most useful in 21 cm cosmology when it preserves physically relevant structure and propagates uncertainty explicitly, rather than acting as an opaque replacement for the underlying forward model.

Recent baryon acoustic oscillation (BAO) distance measurements, when combined with Cosmic Microwave Background (CMB) observations in the $\Lambda$CDM framework, lead to a preference for negative neutrino masses. We investigate whether this neutrino mass anomaly can be alleviated by a class of astrophysically motivated reionization histories. Using a frequentist analysis, we find that some reionization histories can move the best-fit value of $\sum m_\nu$ to a positive value and bring $\sum m_\nu\simeq0.06~{\rm eV}$ into the 95\% confidence interval. To separate the effect of the total optical depth from that of the details of the reionization history, we compare a high-$\tau$ history with a two-step tanh-like reionization history of the same $\tau$. The resulting $\Delta\chi^2(\sum m_\nu)$ profiles are nearly identical. This indicates that the effect is mainly driven by the total optical depth, while the details of the reionization history play only a minor role.

Identifying remnant radio-loud active galactic nuclei (AGNs) is challenging due to their diverse morphological and spectral characteristics. Using three-dimensional hydrodynamic simulations of 15 radio galaxies, we investigate how the spectral evolution of remnants depends on progenitor power, active lifetime, environment, and underlying dynamics. The simulations span low-density group and high-density cluster environments re-gridded from smooth-particle-hydrodynamic cosmological simulations. The resulting remnants exhibit a wide range of morphologies, from amorphous structures to double-lobed forms. We find that jet power correlates with the spectral slope. As the remnant lobes evolve, we find surface brightness depends strongly on environment: group remnants are systematically dimmer and more amorphous than cluster remnants, highlighting a potential observational bias against these low-surface-brightness sources. In our models, we estimate that the peak surface brightness of a low-redshift, 50 Myr-old remnant from a low-power progenitor in a 10^{13} M_sun group environment should be routinely detectable at the 3{\sigma} level with LOFAR, although 20-30% of the emission would remain undetectable within a reasonable integration time. We find young remnants exhibit low-frequency (150-1400 MHz) spectral indices that overlap with active sources, and follow a consistent and established spectral-evolution sequence: significant curvature ({\alpha}_{1400}^{6000} - {\alpha}_{150}^{1400} > 0.5) develops before an ultra-steep low-frequency index ({\alpha}_{150}^{1400} > 1.2). The results presented in this work are intended as a reference point for current and upcoming low-frequency studies of radio remnants.

We present a study that combines remote-sensing and in-situ observations of coronal mass ejections (CMEs) interacting with the nearby heliospheric current sheet (HCS). The sequence of eruptive events under study culminates in the largest directly observed flare of solar cycle 25 on 3 October 2024, producing a fast halo CME. Their source region can be linked to a so-called nested active region (or active longitude) that persisted over several solar rotations. Such long-lived regions reflect deep-seated magnetic structures that shape the global magnetic field configuration. By applying the drag-based CME propagation model, we connect the near-Sun observations from several CMEs during that activity period with in-situ measurements. While one of the CMEs propagated on the opposite side of the HCS from Earth, and therefore did not produce in-situ signatures near Earth, we detect, over the period October 5-10, 2024, a complex of HCS and CME structures propagating together with a shock ahead of them. The HCS seems to be locally replaced by the CME signatures, leading to a long-duration sector reversal of more than 48 hours. This event highlights the intrinsic connection between solar surface structures, the global magnetic field, and the evolution of complex eruptive events.

The search for extraterrestrial intelligence (SETI) has traditionally focused on the detection of narrowband electromagnetic beacons. However, terrestrial technology is increasingly evolving toward distributed, low-power, wideband digital infrastructure. The strict adherence to narrowband filtering that characterises most SETI surveys, therefore, risks discarding the aggregate leakage signatures of advanced civilisations by systematically misclassifying them as unstructured noise. We investigate the feasibility of detecting such planetary-scale broadband radio technosignatures (BRaTs) using a hierarchical observational framework. In this tiered approach, wide-field radio surveys conducted by next-generation arrays (such as the SKA and its precursors) perform the initial deep-field observations, with targeted Very Long Baseline Interferometry (VLBI) providing the definitive, high-resolution follow-up. Because broadband continuum emission is largely insensitive to Doppler drift, long-duration "SETI Deep Fields" are observationally viable, extending the accessible detection volume for Kardashev Type I leakage to 100 pc. To distinguish these signals from other astrophysical confounders, a multi-parameter diagnostic framework is proposed. Candidate technosignatures are identified through a convergence of high brightness temperatures, negligible circular polarisation, spectral non-uniformity, interstellar scintillation, and sub-milliarcsecond astrometric co-motion with nearby Galactic stars/exoplanets.

We study test-particle dynamics in the equatorial plane of a Kerr-Bertotti-Robinson black hole (BH) immersed in a cloud of strings (CS), with mass M , rotation a, magnetic parameter B, and string parameter {\alpha}. Using the Hamilton formalism we recover the effective potential Ueff and the conditions for circular motion, and we compute the specific energy E and specific angular momentum L together with the radial, vertical, and azimuthal epicyclic frequencies {\nu}r , {\nu}{\theta} , {\nu}{\phi}. Going beyond the analytic setup, we provide the first numerical mapping of the innermost stable circular orbit (ISCO) for this background and tabulate rISCO, EISCO, LISCO, and the accretion efficiency {\eta} = 1 - EISCO for both co- and counter-rotating motion across a wide (a, B, {\alpha}) grid. The CS parameter pushes the ISCO outward and raises {\eta} from 0.057 in Schwarzschild to above 0.25 for {\alpha} = 0.30 at a = 0.9. We then connect the model with observed twin-peak high-frequency quasi-periodic oscillations (QPOs) in three microquasars (GRO J1655-40, XTE J1550-564, GRS 1915+105) using the relativistic-precession (RP) model and find \{chi}^2-minimum fits with {\alpha} < 0.13. A general-relativistic hydrodynamical (GRH) study of Bondi-Hoyle-Lyttleton (BHL) accretion completes the picture: the CS contribution sustains shock-cone instabilities, redistributes power-spectral-density (PSD) peaks, and produces low-frequency QPO-like components that distinguish KBR+CS from pure Kerr or KBR.

The formation timescale of the Milky Way thick disk is one of the central debates in Galactic archaeology. The age-metallicity relation (AMR), formation timescale, and chemical evolution gradients are frequently used to infer a rapid assembly, short-timescale enrichment, and bursty formation history of the thick disk. However, stellar ages are not directly observable, introducing the potential risk that inferred ages may harbor a systematic compression tied to observational quality. In this paper, we use the same stellar sample and identical physical covariate matching conditions, but two independent age scales--spectroscopic inferred ages (astroNN) and asteroseismic ages (APOKASC-3)--to compare the observable signatures of the thick-disk formation history. We find that several key observables previously supporting a rapid thick-disk formation are systematically weakened under seismic anchoring: the AMR slope flattens from -3.29 to -1.86 Gyr dex-1 (Delta a = +1.43), the formation timescale widens from 3.04 to 3.55 Gyr, and the peak formation age shifts from 9.1 to 6.0 Gyr. Through transport inversion experiments, we further show that additive noise can only broaden the age distribution and cannot reproduce the above pattern, whereas a compressive transport map (lambda < 1) simultaneously reproduces a narrower age distribution, a steeper AMR, and rapid-formation-like observables. This result indicates that the compression transformation itself is sufficient to generate rapid-formation-friendly observables without requiring an intrinsically bursty formation history. Our findings reveal that statistical interpretations of the Milky Way formation history may depend sensitively on the stellar age definition itself.

Test-mass thermal noise is a limiting noise source for current and next-generation ground-based gravitational-wave observatories. Uniform-intensity higher-order laser beams, including Laguerre-Gaussian (LG) and Hermite-Gaussian (HG) modes, have been proposed as alternatives to the fundamental Gaussian beam due to their thermal-noise advantages. As interferometer power increases toward the megawatt regime, thermal aberrations from absorption in the test-mass coatings become increasingly significant. In this work, we quantify the robustness of higher-order modes against absorption-induced thermal deformation. We show that, under identical operating conditions, higher-order modes produce substantially more uniform thermal distortions than the fundamental mode, requiring significantly less thermal compensation power. The optimal curvature correction is reduced to 33% for the LG$_{2,2}$ mode and 24% for the HG$_{3,3}$ mode relative to the fundamental mode. We further show that the residual thermal deformation of higher-order modes results in lower optical loss, larger cavity power buildup, and improved modal purity in an aLIGO-like cavity. In addition, astigmatism compensation further enhances the intracavity purity of HG modes under self-heating-induced deformation. These results demonstrate that higher-order modes not only mitigate thermal noise but also intrinsically reduce beam self-heating effects, making them promising candidates for future high-power gravitational-wave interferometers.

The Microchannel X-ray Telescope (MXT) is one of four instruments on the Space-based multi-band astronomical Variable Objects Monitor (SVOM) satellite mission, launched on the 22nd June 2024. The MXT is a narrow-field-optimised, lobster eye X-ray focusing telescope, consisting of an array of 25 square MPOs, with a focal length of 1.14 m and working in the energy band 0.2 - 10 keV. The design of the MXT optic (MOP) is optimised to give a 1 degree FoV to match the detector size, but the optic has the unique characteristics of a lobster eye design, with a wide FoV of 6 degree diameter, and a PSF, which is constant over the entire FoV. The MPOs on the Flight Module (FM) MOP have a pore size of 40 um giving the optimum thicknesses across the aperture of 2.4 mm in the centre and 1.2 mm at the edges. Using specific target sources, the in-orbit calibration of the optic is here described, and compared to the extensive on-ground calibration, which was carried out at the PANTER test facility, MPE, Germany. The design and limitations of the electron diverter, situated directly behind the optic, are also discussed.

The transition region is a thin inhomogeneous region where Alfven waves' energy fluxes generated elsewhere in the solar atmosphere are effectively converted. Large-scale kinetic Alfven wave propagations, transmission, and reflection processes across the transition region are examined. The two-fluid model is adopted, and a study is conducted of how the kinetic Alfven wave electric displacement field changes across the transition region. The analysis outcomes are: the kinetic Alfven wave and electrostatic ion sound waves are coupled by the transverse wavenumber; wave electric field components (normal to the transition region) become enhanced up to two orders; the energy fluxes of transmitted kinetic Alfven waves are re-directed almost horizontally along the transition region, an evanescent electric field zone of enhanced intensity is induced tightly beyond the transition region; the ponderomotive force that emerges in that zone due to the reflected kinetic Alfven waves accelerates plasma particle upwards compared to their initial energy in the upper chromosphere.

Context: Recent gamma-ray observations of young star clusters revealed that stellar wind termination shocks accelerate particles, with the energy reservoir provided by the mechanical power of massive-star winds. Aims: Our goal is to identify promising targets for future gamma-ray studies of stellar clusters and associations powered by massive stars. As the wind power of a single Wolf-Rayet (WR) star can rival the cumulative wind power of the most massive clusters, we also investigate isolated WR stars, many of which are indeed isolated. Methods: We ranked a large sample of stellar clusters and associations according to the number of member WR stars divided by the distance squared, a quantity proportional to the expected gamma-ray signal, and searched for spatial correlations with known gamma-ray sources. We repeated the same procedure for individual WR stars with known wind mechanical powers and distances. Results: We found a hint ($\lesssim 3 \sigma$ confidence) for a correlation between WR-hosting clusters and unidentified GeV gamma-ray sources, and identified new spatial associations for 11 clusters. We also found spatial coincidences between 4 isolated WR stars (WR110, WR114, WR111, and WR14) and unidentified gamma-ray sources. Although no significant correlation is found for isolated WR stars as a population, these 4 objects exhibit particularly large wind-power-to-distance-squared ratios, a necessary condition for detectability with current instruments. Assuming the gamma-ray emission is powered by WR winds, it can be interpreted as arising from interactions between particles accelerated at the wind termination shock and ambient matter or radiation fields. Conclusions: Since the wind power of an individual WR star can rival that of an entire stellar cluster, we provide a ranking of stellar clusters and isolated WR stars that may constitute potential gamma-ray emitters.

Conventional cosmological initial condition generators are designed exclusively for fully periodic cubic domains and cannot produce the non-periodic, observer-centric configurations required by stereographically projected N-body codes such as StePS. We present STEPSIC, an open-source initial condition generator that extends Lagrangian perturbation theory-based initial conditions to the spherical and cylindrical geometries used by StePS, while also supporting cuboid domains with arbitrary aspect ratios. The code constructs Gaussian random density fields on anisotropy-free Fourier grids with cubic voxels, applies first- and second-order LPT to obtain displacement and velocity fields, and interpolates these onto particles via B-spline mass-assignment kernels with Fourier-space deconvolution. For stereographic geometries, a multiresolution scheme maps displacement fields across the radially varying particle mass resolution intrinsic to the projection. Both standard and paired-and-fixed variance-reduced realizations are supported. In periodic cubic boxes, the recovered matter power spectrum agrees with the input linear theory prediction to better than 0.5% up to half the Nyquist wavenumber, independent of box aspect ratio (tested up to 10:1). Cross-validation against monofonic using identical white noise fields yields sub-percent power spectrum agreement, with a small residual offset consistent with differences between two independent implementations. Full N-body evolution of matched cylindrical StePS runs confirms that second-order LPT correctly suppresses the 2-3% transient power excess present in first-order initial conditions.

Pulse profiles probe the emission geometry of accreting X-ray pulsars, but their observed shapes may depend on instrumental response and observational setup. The pulsed fraction spectrum provides a compact spectro-timing observable that can both trace localized spectral features and serve as a quantitative cross-calibration diagnostic. We assess the consistency of energy-resolved pulse profiles obtained with simultaneous XMM-Newton/EPIC-pn and NuSTAR/FPM observations of Vela X-1, and investigate the broadband pulsed fraction spectrum as a diagnostic of spectral features from 1 to 70 keV. We construct energy-phase matrices for both instruments and derive pulsed fraction spectra after carefully accounting for instrumental and observational effects. We quantify the residual systematics in the overlapping 3-10 keV band. We then model the broadband pulsed fraction spectra phenomenologically and search for timing signatures of spectral features. After correcting for instrumental effects, the pulsed fraction spectra derived strictly over the common exposure intervals of the two instruments agree within 5% in their overlapping 3-10 keV range. Remaining discrepancies larger than 5% are confined to the iron-line region and can be attributed to the different energy resolutions of the two instruments. The broadband pulsed fraction spectrum reveals significant localized features corresponding to known emission lines in the soft band and to cyclotron resonant scattering features. Orbital-phase-resolved modeling of the EPIC-pn pulsed fraction spectrum shows that the soft-band features depend strongly on the equivalent absorption column, with emission-line signatures becoming progressively suppressed during highly absorbed intervals. The pulsed fraction spectrum serves both as a quantitative cross-calibration diagnostic and as a powerful spectro-timing diagnostic.

This work spotlights the Statistical Adiabatic Channel Model as an efficient and accurate method for deriving low temperature (de)-excitation rate coefficients for collisions induced by heavy projectiles. For such systems, fully quantum treatments become intractable, while quasi-classical methods fail at low temperature. Here, we demonstrate that the Statistical Adiabatic Channel Model overcomes these limitations by combining statistical sampling with an adiabatic channel representation. Its application to the HCN and HNC isomers colliding with CO yields rate coefficients in quantitative agreement with full quantum results benchmarked for the lowest total angular momentum. These systems are relevant for modeling cometary comae, where reliable molecular data remain scarce. Remarkably, this approach also reproduces near-resonant energy transfer and isomeric effects, demonstrating that essential quantum and structural features can be captured within a statistical framework.

Compact protoplanetary discs are becoming increasingly prominent in observations. Their dispersal pathways may differ substantially from those of extended discs. We aim to quantify the role of the disc outer radius in internal photoevaporation, provide a simple scaling relation for compact discs, and test whether the resulting evolutionary tracks reproduce the observed inside-out clearing of young stellar populations. We performed radiation-hydrodynamic simulations of X-ray-driven photoevaporation for discs with different outer radii, and derived the dependence of the total mass-loss rate on the cut-off radius. We find that the surface mass-loss profiles are nearly independent of disc size, but their integrated wind rates are reduced according to the cumulative mass-loss rate distribution. We incorporated this scaling into disc population synthesis models. When the internal photoevaporation is applied only up to the cut-off radius compact discs evolve via inside-out clearing consistent with observational diagnostics, while when the cut-off radius is not considered, the disc spreading is hindered and the disc dispersal proceeds from the outside-in. The introduction of mild external photoevaporation present in nearby star forming regions cannot prevent the disc spreading when the cut-off radius prescription is included, but it can much better explain the evolution of disc radii as a function of time. Disc dispersal prescriptions must include the dependence on disc cut-off radius to capture the evolution of compact discs. The proposed scaling provides a simple, physically motivated correction that better predicts the growing observational evidence for compact discs and inside-out dispersal.

Silicon photomultipliers are increasingly replacing photomultiplier tubes in Cherenkov telescope cameras, but achieving single-photoelectron resolution with nanosecond timing in a low-noise, scalable detector system remains challenging. We present a co-designed SiPM sensor and front-end application specific integrated circuit (ASIC) that meets these requirements. The custom hexagonal sensor, developed with Hamamatsu Photonics, incorporates an integrated optical filter and fourfold pixel segmentation. The readout is performed by a second prototype of the FANSIC ASIC, optimized for this application and fabricated in 65~nm standard CMOS technology, it provides eight channels with on-chip analog summing of sub-channels on a $3.5\times 3.5~\mathrm{mm}^2$ die, while consuming only 24~mW per channel. We demonstrate clear single-photoelectron peak separation with a gain of $2.7 \times 10^{-12}~ \mathrm{V \cdot s}$ , and an impulse response below 4~ns full width at half maximum with a 1.7 ns rise time, preserving the nanosecond-scale structure of Cherenkov pulses. The system responds linearly from 1 to 130 photoelectrons, and 55 distinct photoelectron peaks are resolved by varying the source intensity. These results demonstrate that the integrated sensor-electronics architecture delivers the speed, resolution, and dynamic range required for imaging atmospheric Cherenkov telescopes, and provides a scalable path toward large-area camera modules.

We present the Swift Analysis Pipeline for Lightcurve Extraction (SAPLE), a semi-automated pipeline to extract the Swift-UVOT and Swift-XRT data products and spectral information (magnitudes, photon indices, and fluxes) for a set of observations of any point source of interest. This pipeline is not meant to substitute, but to complement the tools the Swift team has already set up. Specifically, SAPLE provides a Swift-UVOT semi-automated pipeline that also returns the absorption corrected specific fluxes for any observation and filter of interest, a tool which to our knowledge is not publicly available to the community yet. Moreover, for Swift-XRT, SAPLE enables the user to extract a lightcurve of both flux and photon index (with associated uncertainties), assuming a redshifted powerlaw spectrum. The main codes are available through a GitHub repository (L. Marcotulli & N. Torres-Albà 2026), and the following paper summarizes the main steps of the analysis.

We present a comparison of the two most recent and comprehensive Spitzer phase curve studies - Dang et al. (2025) and Swain et al. (2025) - which report analyses of the Spitzer 4.5 $\mu$m phase curves. The studies employ different approaches for correcting instrument systematics and they also use different approaches for selecting the optimal exoplanet system parameters. To evaluate the level of consistency between the two studies, we compared the constraints on the ratio of planet-to-star radii ($R_P/R_\star$), eclipse depth ($F_P/F_\star$), phase curve amplitude ($A$), and phase curve offset ($\phi$). We find that the two studies produce similar results at the population level although results for individual planets can vary, especially for phase curve offset values. We examined the difference of planet system parameters to see if inconsistencies in individual planet results were due to data reduction methods or system parameter choices. We also examined whether the system parameters used by both studies were consistent with Kepler's third law. During this comparison, we identified one case where stellar mass, planet semi-major axis, and orbital period did not follow Kepler's law even though the values were all compiled from the same publication. To assess whether this kind of discrepancy was recurrent, we recalculated the orbital periods using Kepler's third law and compared them with the values listed in the NASA Exoplanet Archive. Our detailed analysis of archival system parameters strongly suggests that testing reported/selected parameters for consistency with Kepler's third law is worthwhile.

Neutrino flavor evolution in dense astrophysical environments such as core-collapse supernova (CCSN) is influenced by collective effects. While the Fast Flavor Instability (FFI) and the Collisional Flavor Instability (CFI) are recognized as key drivers of rapid flavor conversion, their non-linear competition with collisional damping in spatially varying environments remains poorly understood. Motivated by recent findings that FFI and resonance-like CFI co-occur in the post-bounce phase in CCSN, we scrutinize their dynamic competitions and asymptotic states. To this end, we perform numerical simulations of the quantum kinetic neutrino transport, incorporating both spatial advection and the collision terms. We demonstrate that the interplay between these coexisting neutrino flavor instabilities and collisions leads to rich dynamics. Rather than merely inducing simple decoherence, collisional damping can substantially alter the overall dynamics of collective flavor oscillations, driving the system through complex evolutionary pathways. In all cases where flavor instability develops, we find that the system converges to the same flavor-equilibrated asymptotic state, despite the diversity of intermediate dynamics. Our results suggest that this dynamic competition could alter the widely accepted picture of collisionless FFI, highlighting the need to incorporate realistic collisional effects into studies of flavor conversions in CCSN models.

The persistent challenge of creating stable de Sitter vacua within string theory undermines the observational validity of the $\Lambda$ cold dark matter (CDM) model. This difficulty suggests that the concordance model of cosmology, characterized by a constant dark energy $\Lambda$, may reside in the swampland of inconsistent quantum gravity theories rather than the string landscape of consistent ones. Recent observational data, particularly from the Dark Energy Spectroscopic Instrument (DESI), have significantly challenged $\Lambda$CDM cosmology. Specifically, the combination of DESI baryon acoustic oscillation measurements with cosmological surveys seem to indicate a preference for a dynamic, time-evolving dark energy rather than a constant, with roughly 10\% reduction in density over the last several billion years. This review summarizes significant advancements made over the past two years in linking DESI findings to string-inspired scenarios.

Combined HI emission-absorption studies constrain the spin temperature and phase structure of the neutral atomic hydrogen interstellar medium (ISM), but have largely been limited to the Milky Way and the Local Group. We extend this technique to galaxies at distances of 7-22 Mpc using deep data from the MeerKAT HI Observations of Nearby Galactic Objects - Observing Southern Emitters (MHONGOOSE) survey, and quantify the detection fraction and Cold Neutral Medium (CNM) properties at these distances. We search for HI absorption toward 56 background continuum sources in 21 out of the 30 MHONGOOSE galaxies (with nine galaxies lacking suitable background sources), and detect absorption associated with the galaxies' HI disks in three cases: one sight line in NGC 289 and two in NGC 7424. This corresponds to detection rates of 3/56 (5 percent) for the full sample and 3/31 (10 percent) for a clean sub-sample of sight lines, considering only unresolved background sources behind 14 low-inclination galaxies. Detections occur only where both the continuum flux and the foreground \HI column density are high, with optical-depth sensitivity as the primary limiting factor. For the detected sight lines, we model the absorption and emission spectra to derive spin temperatures and CNM fractions using the standard combined emission-absorption method. The CNM spin temperatures and line widths are comparable to Local Group measurements, but the inferred CNM fractions are systematically lower. We argue that this difference is primarily a resolution effect: at the distances of our galaxies, the emission spectra average over several hundred parsecs, diluting structured CNM relative to the smoother Warm Neutral Medium (WNM). This demonstrates that emission-absorption analyses can be extended beyond the Local Group, provided that care is taken in constructing representative emission spectra.

Relativistic jets launched from stellar-mass compact objects embedded in the accretion disk of an active galactic nucleus (AGN) can produce nonthermal emission upon successfully breaking out of the disk. In this paper, we present a comprehensive study of the long-term propagation dynamics and broadband nonthermal radiation signatures of such jets in a realistic AGN environment, explicitly modeled as wind outflows. Our modeling reveals two distinct features imprinted by the high-density AGN medium: rapid deceleration of the jet ejecta, accompanied by a prompt downshift of the emission spectral energy distribution, and persistently strong synchrotron self-absorption, giving rise to a prominent quasi-thermal hump in the emission spectrum. Crucially, both gamma-ray burst jets and jets powered by accreting binary black hole merger remnants can produce detectable multi-wavelength emissions that substantially outshine the AGN background. Moreover, the short time delays between gravitational wave triggers and these electromagnetic counterparts--typically less than 106 s--greatly facilitate secure multi-messenger associations. Besides, our findings highlight that interaction-induced radiation from AGN-embedded jet systems offers a powerful diagnostic probe of the spatial distribution,density structure, and physical properties of the AGN medium.

We investigate whether reconnection of small-scale current sheets in transrelativistic supermassive black hole (SMBH) coronae can supply the nonthermal protons needed for high-energy neutrino emission, using NGC 1068 as a test case. We model the corona as a strongly turbulent, low-$\beta$, collisionless hydrogen plasma with characteristic size $r_{\rm co}$, magnetic field strength $B$, proton density $n_p$, and radiation energy density $u_{\rm rad}$. Combining the observed IceCube-band neutrino luminosity with the X-ray luminosity and Thomson optical depth reduces these coronal quantities to a one-parameter family. Across this family, the proton magnetization $\sigma_p \equiv B^2/(4\pi n_p m_p c^2)$ is transrelativistic with $\sigma_p \sim 0.3$. In this regime, we show that repeated encounters with intermittent reconnecting current sheets can energize suprathermal protons up to tens of PeV before photomeson cooling limits further acceleration. These injected particles may then be further processed by stochastic interactions with the turbulent cascade. Motivated by PIC simulations of strong turbulence at comparable magnetization, we adopt a nonthermal proton spectrum with an independently specified index and find that the predicted TeV spectral shape is broadly consistent with NGC~1068 without fitting the proton spectral slope.

With the detection of gravitational waves (GWs), multi-messenger astronomy has opened a new window for advancing our understanding of astrophysics, dense matter, gravitation, and cosmology. The GW sources detected to date are from mergers of compact object binaries, which possess the potential to generate detectable electromagnetic (EM) counterparts. Searching for associations between GW signals and their EM counterparts is an essential step toward enabling subsequent multi-messenger studies. In the era of next-generation GW and EM detectors, the rapid increase in the number of events brings not only unprecedented scientific opportunities, but also substantial challenges to the existing data analysis paradigm. To help address these challenges, we develop GW-Eyes, an agentic framework powered by large language models (LLMs). For the first time, GW-Eyes integrates domain-specific tools and autonomously performs counterpart association tasks between GW and candidate EM events. It supports natural language interaction to assist human experts with auxiliary tasks such as catalog management, skymap visualization, and rapid verification. Our framework leverages the complex decision-making capabilities of LLMs and their traceable reasoning processes, offering a new perspective to the multi-messenger astronomy.

We summarize a set of constraints on a proposed population of primordial, ultra-massive, electrically charged compact objects, which we call astrons. The analysis combines charge generation, charge saturation, persistence of the charge in an ionized medium, screening by the intergalactic plasma, the Reissner--Nordström geometry of highly charged compact objects, lensing, and the cosmological implications of a sparse charged population. We also discuss the possible relation to the early structures revealed by the James Webb Space Telescope: if astrons are relevant there, they would be primordial dark seeds rather than luminous objects directly observed at high redshift. The resulting scenario is sharply constrained. Ordinary accretion saturation gives charges far below the large-charge phenomenological benchmark, screening is a serious plasma-physics issue, and a large charge can place the exterior geometry deep in the super-extremal regime. As expected at the level of a homogeneous Friedmann--Lema\^ıtre--Robertson--Walker (FLRW) description, the interaction energy of a population of charged objects scales as \(a^{-4}\), so the simplest perfect-fluid reduction does not generate asymptotically late-time acceleration; any acceleration era tied to that homogeneous component can only be transitory. The astron scenario should be regarded as a constrained framework whose viability depends on plasma physics and on a cosmological treatment beyond the homogeneous approximation.

Observational data have revealed a clear dichotomy in the [{\alpha}/Fe] vs. [Fe/H] diagram of the Milky Way thick and thin disc stars. Many recent studies have shown evidences of a co-evolution phase between the high- and low-{\alpha} disc sequences as well as the presence of very old low-{\alpha} stars. We aim to revise the parallel chemical evolution model that assumes two parallel histories of star formation for the two discs, by considering a pre-enriched delayed second infall episode in our revised scenario. By means of our chemical evolution models, we aim to explore the effects of a phase of co-evolution and the presence of old low-{\alpha} stars, as recently observed. We consider a new version of the parallel scenario for the Milky Way thick and thin disc formation, which consists into two distinct infall episodes of slightly pre-enriched gas. The gas is considered to be extragalactic but possibly contaminated by chemically enriched gas of a massive dwarf galaxy as Gaia-Enceladus, which merged with the Milky Way at least 10 Gyrs ago. Moreover, we test in our model observationally derived star formation histories of kinematically selected thick and thin discs, suggesting that the star formation is triggered by the passages of the Sagittarius galaxy. Our models can well explain the [{\alpha}/Fe] vs. [Fe/H] diagram from APOGEE DR17. Our revised chemical evolution model with a pre-enriched and delayed (roughly 1 Gyr) second infall episode, explains not only the abundance patterns of high- and low-{\alpha} stars but also stellar age distributions for the selected observational sample. We predict a short co-evolution period in between the two phases and we can explain the observed old low-{\alpha} stars, but still further data for precise stellar ages would be needed to put more stringent constraints on their physical nature.

$\epsilon$ Ori (HD 37128) is an $\alpha$ Cygni variable characterized by irregular and small amplitude variations. From TESS observations, we find the presence of stochastic low-frequency variability in this star. We have constructed a sequence of models for this star in the mass range of 30 to 70 M$_{\odot}$, using recently derived values of luminosity (log $(L/L_{\odot})$ = 5.92) and effective temperature. In these considered models, both radial and non-radial linear stability analyses have been performed. Low-order radial modes are excited in models having mass below 62 M$_{\odot}$. These radially excited modes have periods ranging from 6.8 days for the fundamental mode to a few hours for higher-order modes. Similar to the case of radial modes, several non-radial modes are found to be unstable in models having higher luminosity-to-mass ratios. Linear stability analysis for the case of $l$ = 2 and $l$ = 4 reveals the presence of a strongly unstable mode in models having a mass below 40 M$_{\odot}$. This mode is found to be unstable in all the considered models and the strength of the instability varies as a function of harmonic degree. The non-adiabatic reversible approximation reveals that the origin of instabilities associated with the low-order modes is indeed linked with strange modes. To find out the consequence of radial instabilities, non-linear numerical simulations have been performed in selected models of $\epsilon$ Ori. In the non-linear regime, these instabilities lead to the envelope inflation, finite amplitude regular and irregular pulsations consistent with an $\alpha$ Cygni variable.

In the search for life in the cosmos, biopigments on exoplanet surfaces are a critical target. Such pigments have been detected in Earth's spectrum (by the Galileo spacecraft and in Earthshine) via the "vegetation" or "photosynthesis red edge" (VRE or PRE), a sharp, step-like increase in reflectance with increasing wavelength at ~700 nm. Future space telescopes like the Habitable Worlds Observatory (HWO) are designed to obtain disk-integrated spectra of Earth-like exoplanets in the visible-to-near-infrared to identify such features. However, there has been no systematic analysis of the occurrence of similar reflectance edges among minerals of non-biological origin. Here, we use existing databases of mineral reflectance spectra to explore the risk that minerals may present false positives in the search for biopigments on exoplanets. We find that several sulfide and tectosilicate minerals, as well as the prebiotically important cyanide salt, potassium ferrocyanide, have PRE-like features. We characterize these features in order to assess how they may be distinguished from biopigments. We conclude that the future evaluation of the biogenicity of PRE-like features in exoplanet reflectance spectra can be informed by the atmospheric context, but may require an assessment of the prior probability of non-biological and biological hypotheses about the surface materials of exoplanets.

Context. Solar wavefront sensing has been a challenge for astrophysical instrumentalists, due to the low contrast between the Sun and the sky background compared to night-time observations, which limits the performance of adaptive optics systems. Aims. Wavefront correction in solar physics requires the analysis of extended images; meanwhile, at night the displacement of a punctual object is analysed. This technique limits the spatial resolution, and therefore the accuracy in the wavefront reconstruction. Methods. To solve this problem, a new method of direct wavefront sensing without the need for image formation was explored for this work. A novel and promising technology called integrated photonics was used to accomplish this task. It allows the direct measurement of phase differences across the wavefront without the need to form images, using the principle of interferometry. This technology offers a low-consumption, miniaturised solution to astrophysical problems. Results. For this work a mathematical model was derived to characterise the behaviour of the proposed wavefront sensor. The proposed system was verified and simulated using a Python-based adaptive optics simulator. These simulations demonstrate the physical behaviour of the proposed wavefront sensor and highlight the factors that must be taken into account for its correct functioning.

In multifield inflation driven by $d$ scalar fields, $O (d)$ symmetry renders the number of fields irrelevant at classical level. This ceases to be the case once stochastic effects are accommodated. The statistical quantities such as the mean number and the variance of $e$-folds as well as the primordial power spectrum and its scale dependence are perturbatively calculated in a small-noise regime. In particular, a general formula is derived for arbitrary higher-order statistical moments of the stochastic number of $e$-folds at all perturbative orders, keeping the dependence on the number of fields fully analytical. It is also discussed that the requirement for inflation to be successfully terminated puts a theoretical bound on the number of fields from above. Those general results are demonstrated for several $O (d)$-symmetric models.

The neutral hydrogen (Hi) 21-cm line serves as a powerful tracer of the neutral interstellar medium (ISM). Thermal stability analysis suggests that the neutral ISM is bistable in nature, consisting of the cold neutral medium (CNM) embedded within the warm neutral medium (WNM), both in approximate thermal pressure equilibrium. When turbulence is incorporated into the numerical simulations, a third thermally unstable medium (UNM) emerges between the CNM and the WNM. Although observational studies support the existence of this intermediate phase, a clear empirical correlation between the fraction of the UNM gas and the strength of the turbulence remains elusive. In this study, we investigate the various phases of neutral ISM using Hi 21-cm emission-absorption spectra from the publicly available GWA and LAB surveys and compare it with TIGRESS-NCR and TIGRESS-CLASSIC numerical simulations. From our observational modeling, we find that 19.8% of the gas reside in the CNM phase, 32.5% in the UNM phase, and 47.8% in the WNM phase, assuming phase boundaries defined by spin temperature: T_s < 250 K for the CNM, 250 K < T_s < 5000 K for the UNM, and T_s > 5000 K for the WNM. These results are entirely in agreement with the TIGRESS-NCR numerical simulation. We further expect that deep, sensitive absorption studies with the Square Kilometre Array (SKA) or the Next Generation Very Large Array (ngVLA) or existing Upgraded Giant Metrewave Radio Telescope (uGMRT) capable of robustly detecting WNM clouds in absorption will place more tighter observational constraints on the fraction of the gas in three different phases of the neutral ISM.

Simulation-based inference (SBI) enables parameter inference by training neural networks on forward simulations. It is being applied both for intractable likelihoods as well as under time constraints on the posterior sampling. After motivating situations in which SBI is useful, we give a pedagogical description of the basic techniques. These are posterior, likelihood, and ratio estimation. Alternatives, sequential versions, and learned summaries are discussed briefly. We provide a brief guide to choosing among the techniques in practical scenarios. SBI needs to be verified through diagnostics since failures can be subtle but would invalidate the inference result. We explain the most common diagnostic techniques. We briefly list some recent SBI applications in the cosmology and astrophysics literature. Before concluding, we discuss current methodological challenges. We identify training with limited simulation budgets as the critical problem for applications to cosmology and astrophysics.

In this work, we perform a numerical study of three Starobinsky--type inflationary scenarios: the $\alpha$--Starobinsky inflationary model, the power--law Starobinsky inflationary model, and the power--law $\alpha$--Starobinsky inflationary model. For an appropriate choice of parameters, each scenario reproduces the standard Starobinsky limit. For each case, we derive the relevant slow--roll expressions in order to compute numerically the scalar and tensor power spectra over the corresponding parameter space and evaluate the associated inflationary observables. Finally, we provide a comparative analysis in the $(n_\sca,A_\sca)$ and $(r,n_\sca)$ parameter spaces using contour plots. Our results indicate that, for certain choices of parameters, the $\alpha$--Starobinsky model and the power--law $\alpha$--Starobinsky model are favored by \textit{Planck} 2018 observations.

Wide-field imaging surveys now provide photometry for billions of sources, while spectroscopic observations remain limited, motivating methods that can extract spectroscopic information from photometric data. We present a generative framework for the joint probabilistic inference of galaxy redshifts and rest-frame spectra from broadband photometric fluxes. The model provides a sampling-based estimate of the photometric-redshift probability density function (PDF) for each galaxy, from which accurate point estimates are derived, and reconstructs rest-frame spectra that preserve key spectral properties. We pre-train a spectral autoencoder, SPENDER, on 5 million DESI DR1 spectra to learn a low-dimensional latent space that represents rest-frame spectra. Conditioned on galaxy broadband photometric fluxes, a diffusion model jointly infers the corresponding spectral latent representation and photometric redshift. The inferred latent representation is decoded into a high-resolution rest-frame spectrum, which can be transformed to the observed frame by redshifting and resampling. Sampling from the conditional diffusion model yields a full photometric-redshift PDF for each galaxy, with the resulting point estimates showing a precision comparable to that of a gradient-boosted decision tree model. In most cases, the reconstructed rest-frame spectra reproduce the overall continuum shape and capture the presence of prominent spectral features. For galaxies with sufficiently high signal-to-noise ratios in their observed spectra, the Dn4000 index shows good agreement between the reconstructed spectra and the observed spectra. On average, the spectral reconstruction residuals are close to the noise level of the observed spectra. Latent-diffusion generative modeling enables joint inference of galaxy photometric-redshift PDFs and rest-frame spectra from photometric fluxes.

Although non-greenhouse gases can vary substantially in abundance in Earth-like atmospheres, their climatic influences remain insufficiently understood. To investigate how such gases regulate climate, we vary the abundance of N$_2$ as a representative non-greenhouse component in one-dimensional N$_2$--CO$_2$--H$_2$O model atmospheres. Beyond pressure broadening of absorption lines and Rayleigh scattering emphasized in previous studies, our results show that changes in background N$_2$ pressure influence climate by modifying the amount of atmospheric H$_2$O, producing two effects: altering the thermodynamic lapse rate (H$_2$O-dilute warming) and changing the radiative contribution of H$_2$O to the greenhouse effect (H$_2$O-load warming). The resulting climate response to increasing N$_2$ depends on the CO$_2$ abundance. Under low CO$_2$ conditions, dilution of atmospheric H$_2$O leads to warming, whereas under high CO$_2$ conditions, increased H$_2$O loading also produces warming. At sufficiently high N$_2$ abundances, Rayleigh scattering induces cooling, an effect further amplified by the accompanying decrease in atmospheric H$_2$O. Under high CO$_2$ conditions, however, enhanced H$_2$O loading increases the absorption of stellar radiation and overwhelms the contribution of Rayleigh scattering, causing the cooling response to disappear. These results reveal multiple physical pathways through which non-greenhouse gases influence climate and provide a framework for understanding climate responses and habitability in diverse Earth-like atmospheres.

At the center of our Galaxy lies the bright {\gamma}-ray point-like source HESS J1745-290, which is compatible in position with Sgr A star, although an association between the two remains uncertain. Using data obtained between 2004 and 2019 with the High Energy Stereoscopic System (H.E.S.S.) on the Galactic center region, we studied the variability of HESS J1745-290 over 353 hours of observations collected over 16 years, representing the largest dataset gathered yet on this region at TeV energies. We performed a 3D maximum-likelihood analysis of the central source and the diffuse {\gamma}-ray emission in the Galactic center region. This analysis allowed us to extract the spectral and morphological intrinsic behavior of the two components. By performing this analysis on an annual basis, we derived the light curve of HESS J1745-290 and the diffuse emission over the past 16 years. The 3D maximum-likelihood analysis method allowed us to separate the central source from the overlapping diffuse emission, enabling a recalibration of the former by the latter and alleviating some of the systematic effects. We find no long-term or yearly variability. We also provide an estimate of the sensitivity of H.E.S.S. to variation of this specific source over 16 years. We rule out any yearly gamma-ray flux variation of this source larger than 30 percent, as well as any linear flux variation exceeding 30% over this time period.

Modified gravity theories such as Modified Newtonian Dynamics (MOND) and Scalar-Tensor-Vector Gravity (STVG) have been proposed as alternatives to dark matter, but decisive tests have been hindered by degeneracies between baryonic structure and gravitational laws. Here we break this degeneracy using independent, high-precision constraints: the Milky Way radial rotation curve, vertical phase-space spirals from Gaia, and a broken-exponential stellar disk. A joint reconstruction of the radial and vertical gravitational fields reveals a structural inconsistency in modified gravity -- no model can simultaneously reproduce both observations. Our results strongly disfavor MOND at $>13\sigma$ and STVG at $>4\sigma$. In contrast, dark matter halo models naturally explain the observations, providing a self-consistent test of gravity on galactic scales.

Distance measurements are fundamental to radio pulsars' use as astrophysical probes of General Relativity and the interstellar medium. One of the primary methods for determining pulsar distances is HI kinematics, which leverages the radial velocities of HI absorption and emission features detected along pulsar lines-of-sight. This method necessarily assumes a model for Galactic rotation, our knowledge of which continues to evolve in both accuracy and precision. In this research note, we derive kinematic distances for 66 pulsars with archival HI radial velocity measurements using a state-of-the-art Galactic rotation curve. The results and software are provided in an online repository. Our kinematic distances differ by $<1\sigma$ from published parallaxes for nearly all pulsars in the sample that have both types of distance measurement available. Comparison to the NE2025 Galactic electron density model shows general consistency between measured and predicted distances.

Solar filament/prominence eruptions can significantly perturb geospace when originating from favorable source locations and directions. While stellar analogs have been recently reported, the disk locations and magnetic environments of their source regions remain spatially unresolved on other stars. To bridge this gap, we investigate the typical Sun-as-a-star H$\alpha$ temporal spectral characteristics of solar filament/prominence eruptions with different source region locations (on-disk vs. limb, active region vs. quiet-Sun region). It is revealed that limb eruptions are characterized by blueshifted/redshifted emission caused by the bright off-limb erupting structures, whereas on-disk eruptions may show blueshifted absorptions due to the dark erupting filaments. Among the limb eruptions, front-side limb eruptions usually display line center emission before the blueshifted/redshifted emission, while far-side limb eruptions show the opposite sequence. Moreover, the magnetic environment at source also shapes the spectral characteristics. On-disk filament eruptions from active region exhibit much more intense flare-ribbon-dominated line center emission features compared with those from quiet-Sun region. Limb active region eruptions often show single-wing emissions, whereas large-scale quiet-Sun region (quiescent) prominence eruptions frequently display expansion-induced emission in both wings followed by line center absorption due to the disappearance of bright prominence. These distinct Sun-as-a-star H$\alpha$ spectral characteristics, dependent on eruption location, provide a diagnostic basis for inferring source regions of stellar filament/prominence eruptions from spatially unresolved H$\alpha$ spectra.

Shearing motion and magnetic flux cancellation around the polarity inversion line (PIL) play significant roles in the build-up of free magnetic energy and magnetic flux rope (MFR) in source region of major solar flares. Here we investigate the magnetic evolution of a highly-sheared PIL in active region (AR) 13842, hosting the largest X9.0 flare of Solar Cycle 25. Since 2024 September 29, a positive-polarity pore persistently drifted northward along the western side of the AR's main negative-polarity sunspot. The main sunspot remained stationary until negative-polarity patches successively emerged to its east and approached. Rear-ended by these same-polarity patches, the sunspot then began moving westward toward the opposite-polarity pore around October 1, forming a collisional PIL. Meanwhile, on the PIL's other side, the pore was also rear-ended by same-polarity patches sequentially emerging behind it, accelerating the shearing motion around the PIL, where frequent flux cancellations were also observed. Synchronous rapid accumulation of free magnetic energy and formation of MFR were then observed in the PIL, where multiple major flares successively occurred within two days. Before these large flares, the area and total free energy of the high-free-energy-density PIL region gradually decreased in the photosphere, which could be caused by the initial ascent of MFR before eruption and serve as a precursor of solar eruptions. These results suggest that persistent flux emergences with cross separation directions facilitates rapid formation of collisional shearing PIL and frequent flux cancellations, leading to repeated MFR formations and multiple large flares in a relatively short time.

The redshifted 21\,cm line is an emerging tool in observational cosmology that can serve as a direct probe of the intergalactic medium throughout the cosmic timeline. However, the observation of the cosmological 21\,cm signal from early epochs is extremely challenging in practice, regardless of the scale of interest and redshift. The presence of bright astrophysical foregrounds and residual systematic errors along the line of sight poses challenges for its detection. Machine-learning-based Gaussian process regression\,(ML-GPR) has proven to be the most effective strategy for signal separation in LOFAR and NenuFAR observations to measure the 21\,cm signal power spectrum from the Cosmic Dawn\,(CD) and Epoch of Reionization\,(EoR). In this work, we extend this framework to synthetic CD/EoR SKA1-Low observations to assess its robustness in mitigating residual foregrounds against instrumental and environmental systematic effects. We use our developed end-to-end realistic simulation pipeline (\textsc{21cmE2E}) for SKA-Low observations. Our 4-hour tracking simulation includes extragalactic point sources, the AA* telescope configuration, primary beam response, and error models. The modelled errors incorporate residual antenna-based gain calibration errors, residual ionospheric phase errors, partial de-mixing of the out-of-field sources, and instrumental noise for 1000\,hours of deep integration time. We compare different Bayesian GPR frameworks to assess their ability to suppress residual foreground contamination while minimizing signal loss and providing reliable uncertainty estimates. Our analysis demonstrates that the 21\,cm signal can robustly recover within the $2\sigma$ credible interval for almost all k-modes over the range of $0.06 \leq k \leq 1.0$~h Mpc$^{-1}$.

The Spherical Fourier-Bessel (SFB) basis, in separating the angular and radial modes of the power spectrum, permits a targeted identification and mitigation of systematics in clustering surveys while retaining more cosmological signal than traditional bases. We demonstrate this principle on the eBOSS DR16 LRG and QSO samples, identifying modes which may be contaminated by systematics. Our initial inference on the LRG sample yields an fNL value consistent with zero, while the QSO value is in slight tension with zero. Using the SFB basis, we vary the selection of angular and radial modes to search for inconsistencies in the inferred value of fNL, an indicator of underlying systematics. In the QSO sample, we find evidence (p < 0.005 compared to the same cuts on EZMocks) of a systematic afflicting large physical scales, which is consistent with residual stellar contamination; we also find evidence (p < 0.05) for an unknown systematic in the QSO and LRG samples at the approximate angular plate and imaging scale of eBOSS.

On November 12th, 2025 the LIGO--Virgo--KAGRA (LVK) collaboration reported gravitational waves (GWs) from a compact object merger candidate (S251112cm) with at least one sub-solar mass component. Using the Dark Energy Camera (DECam), the Fraunhofer Telescope at Wendelstein Observatory (FTW), and the Zwicky Transient Facility (ZTF), we surveyed $56\%$ of the GW localization region beginning $2.4$~hours after the GW alert. We find no kilonova (KN) counterpart, and use radiative-transfer models to rule out $42\%$ (ZTF), $68\%$ (DECam), and $92\%$ (FTW) of the KN models as possible emission from this GW candidate. Within the recently proposed disk-fragmentation (``superkilonova'') model for generating sub-solar mass neutron star mergers from stellar core-collapse, the delay between the supernova explosion time and the GW merger time is estimated to be less than a few days. Searching this time window prior to the GW event, we identify and spectroscopically classify a IIb supernova (SN~2025adtq), with a spatial association odds ratio of $\log_{10}\mathcal{I} \approx 4.8$, a chance coincidence probability of ${\sim}2$--$9\%$, and an estimated explosion time ${\sim}2$ days prior to S251112cm. SN~2025adtq is the second Type~IIb supernova found in spatial and temporal coincidence with a sub-solar mass GW candidate, following the previously reported S250818k/SN~2025ulz association; jointly, we measure an odds ratio that favors the association hypothesis over the null, however, when conditioned on finding a coincident supernova by chance, the odds ratio disfavors association. Together, these results provide suggestive but inconclusive evidence for the superkilonova formation channel.

The design of new low-mechanical-loss, high reflectivity mirrors is crucial in the development of the next generation of gravitational-wave observatories. Currently, the state-of-the-art amorphous multilayer reflective coatings which are deposited at the surface of the test masses in interferometric gravitational-wave detectors present the limiting factor in detector sensitivity due to their thermal noise. These coatings require a large number of thin layers to achieve ultra-high reflectivity. However, the thermal noise generated by this type of stack increases with the number of layers used. These dielectric mirrors represent a very mature technology, with current research producing only incremental improvements, highlighting the need for new technical solutions that can address this specific issue. Here, we provide insights into the expected performance of mirrors that combine a resonant metasurface with a multilayer stack. The suggested mirror design ensures the high reflectivity required for interferometric gravitational wave detectors, while using fewer layers of properly selected materials. It allows to reduce the total thickness of the material with the poorest thermal-noise performance, namely TiO2:Ta2O5, by a factor of more than 3, making it a promising option for reducing thermal noise as well

In this work, we propose a novel mechanism for generating gamma rays from the Galactic Center via scattering of cosmic-ray protons off dark matter in the Milky Way halo, in contrast to conventional explanations based on dark matter annihilation. We present two examples of this framework that produce an observable photon signal. In the inelastic dark matter model, cosmic rays up-scatter a lighter dark matter particle, with the subsequent decay of the heavier particle yielding two photons. In the elastic dark matter model, an energetic photon is directly produced in the final state of a 2-to-3 scattering process. We show that, for a range of viable model parameters, this framework provides a fit to the observed Galactic Center gamma-ray excess spectrum comparable to those obtained from dark matter annihilation and millisecond pulsar models. Our results open a new avenue for interpreting gamma-ray observations of the Galactic Center.

Radar sounders are electromagnetic instruments that can probe deep into the subsurface of Earth and other planetary bodies by processing the echo of transmitted radar waves. Conventional approaches for analyzing such data rely on approximate assumptions and often produce point estimates that ignore parameter correlations as well as galactic and measurement noise. We propose a simulation-based inference approach to terrain parameter inversion from radar sounder data, where synthetic observations from a GPU-based simulator are used to train a neural network-based density estimator for neural posterior estimation (NPE). By explicitly conditioning on reference surface assumptions, the proposed framework allows systematic evaluation of posterior robustness to reference surface variability. We demonstrate that our NPE model is well calibrated on simulated data and transferable to real Mars radar profiles, where we analyze terrain parameters using literature-informed reference values.

The data on cosmic antiprotons have reached an outstanding precision on energies spanning from GeV to hundreds of TeV, thanks to the space-based AMS-02 experiment. The balloon-borne GAPS experiment, which just completed its first Antarctic flight, will address antiproton and antideuteron fluxes well below GeV energies. Antinuclei in cosmic rays, as well as being produced by spallation reactions between cosmic-ray nuclei and the atoms of the interstellar medium, may hide contributions from exotic sources, such as particle dark matter annihilation in the Galaxy. In this paper, we present predictions for cosmic antiproton, antideuteron and antihelium fluxes both from secondary and dark matter origin. We use state-of-the-art production spectra, nuclear coalescence for antinuclei, and Galactic propagation models to derive upper limits on the dark matter annihilation cross-section from AMS-02 antiproton data in different propagation scenarios (BIG and QUAINT). We quantify the impact of future GAPS data, showing that its sensitivity to sub-GV antiprotons could improve the $\langle\sigma v\rangle$ constraints by up to an order of magnitude for light DM ($m_{\chi} \lesssim 50$ GeV). For heavier antinuclei, the detection perspective with existing and upcoming experiments are derived for those scenarios consistent with AMS-02 antiproton flux. The detectability of such signals strongly depends on the experiment, the propagation model, and the hadronization tuning. Our analysis underscores the complementarity of antinuclei channels for indirect DM searches and the critical role of low-energy windows in constraining light DM candidates.

Time-dependent tidal interactions during the late inspiral of binary neutron stars encode valuable information about neutron-star structure, but systematically extending the familiar Newtonian mode-sum picture into full general relativity is nontrivial. In this paper, we develop a practical relativistic implementation of mode-sum tidal response for non-rotating neutron stars in Regge-Wheeler gauge. Using near-zone boundary conditions, we systematically define the interior tidal field, the relativistic overlap integrals, and the corresponding mode amplitudes. The good is that the dominant f-mode contribution is remarkably robust, reproducing the direct matching calculation to within $\sim 3$\% across the equations of state we consider. The bad is that the operator governing mode inner product is not positive definite on the full Regge-Wheeler-gauge function space, so the relativistic mode sum truncated at ${\cal{O}}(\omega^2)$ is not expected to strictly converge to the direct matching solution. The subtle is that the tidal field inside the star is not unique, although this ambiguity has only a limited impact on the dominant f-mode response for the classes of extensions studied here. Our results establish the practical utility of relativistic mode-sum approximations, while making clear that their predictive power comes from a controlled low-mode description, rather than from a formally convergent strong-field expansion.

In this work we derive novel analytical expressions for the mass and distance of a Schwarzschild black hole (BH), as well as for the orbital radius of test particles orbiting it, it terms of astrophysical observables measured throughout the entire orbit of the revolving particle. We use a general relativistic method to describe the frequency shifts of photons emitted in the vivinity of a BH by considering two emitters (or two positions of the same emitter) located symmetrically opposite to each other with respect to the observer's line of sight (LOS) when performing measurements along the orbit. Furthermore, the introduction of the redshift rapidity allows us to write independent expressions for the BH mass and its distance to Earth. We also extend our study to the case when astrophysical systems have a peculiar motion and derive the corresponding closed formulas.

The performance of superconducting resonators underpins a wide range of modern quantum technologies, yet their quality factor often deviates at low temperatures from standard Mattis-Bardeen predictions. This discrepancy is often attributed to nonthermal quasiparticles generated by microwave readout power, which limits the sensitivity of superconducting devices. We present a macroscopic model based on modified Rothwarf-Taylor equations that incorporates a power-dependent phonon generation term, providing an explicit relationship between quality factor, bath temperature and readout power. The model shows excellent agreement with temperature sweep measurements of NbN microstrip resonators with \b{eta}-Ta terminations over a wide dynamic range of readout power levels, accurately capturing the transition between thermally-dominated and microwave-induced loss regimes. This framework provides a predictive tool for optimizing superconducting resonators and advancing the design of high-Q devices for quantum sensing and quantum information processing.

To effectively describe the plasma with strong magnetic field, the force-free electrodynamics was introduced, within which the Grad-Shafranov equation plays the key role. The Grad-Shafranov equation governs the global structure of a electromagnetic field in equilibrium with symmetries. It is widely applicable in an amount of scenarios, such as the tokamak, the solar corona, the magnetosphere of Earth, neutron star and black hole, etc. However, in different situations, the Grad-Shafranov equation is expressed differently, and the derivations might be complicated. In this work, via the language of differential form, we provide a general expression of Grad-Shafranov equation, from which the expression in any specific situation can be quickly obtained. Additionally, we present a Lagrangian density for a scalar field whose on-shell condition is precisely the Grad-Shafranov equation.

We investigate a conservative formulation of $f(R,T)$ gravity motivated by a key limitation of several existing approaches: the gravitational function is often reconstructed from a chosen equation of state, making the gravity sector EoS-dependent and compromising universality. To avoid this problem, we reformulate the theory in terms of an effective energy-momentum tensor, so that the conservation law follows from the field equations and Bianchi identities while the gravitational action remains independent of the microphysical EoS. We derive the modified stellar structure equations, establish theoretical consistency conditions including coupling bounds and crust-singularity avoidance, and present the tidal perturbation sector in terms of effective thermodynamic variables and an effective sound speed. We then compute neutron star observables using realistic tabulated EoSs, including mass-radius relations and tidal deformabilities, and compare the model with current astrophysical constraints from massive pulsars, NICER radius measurements, and GW170817.

We study in detail a scenario in which the inflaton scalar field couples to both a visible sector (VS) and a hidden sector (HS). The VS is assumed to contain the Standard Model (SM), while the HS contains a dark matter (DM) candidate. We are in particular interested in a scenario in which the inflaton decays dominantly into the HS degrees of freedom. The DM candidate is taken to be a dark Dirac fermion $\chi$, coupled to a massive dark photon $\gamma'$, a popular model for a HS also known as Dark QED. The inflaton decays into particles of both sectors generate an initial asymmetry between the SM and HS fermion abundances, which we model as being proportional to the ratio of effective Yukawa couplings, $y$ and $y'$. We pay particular attention to the process of thermalisation of the HS, with temperature $T'$, as a function of $y'$ and $\alpha'$, the HS fine structure constant. We investigate the several possible ways of producing the observed DM relic abundance, and their interplay with the reheating of the HS and the transfer of energy between the HS and the VS. Key results, beyond the systematic character of our analysis, include: a new mechanism for DM production, which occurs when DM particles annihilate while still being produced by the inflaton decay; a study of the temperature ratio $\xi = T'/T$ and its relation with the initial energy asymmetry between the HS and VS, as parameterized by $\xi_i = \sqrt{y'/y}$; a reassessment of the domain of viable DM candidates, taking into account the constraints set by unitarity and the thermalisation of the HS, accounting for the LPM effect; and, in cases where the HS does not reach thermal equilibrium, an analysis of how non-thermal DM production fits within the domain of thermal DM candidates.

Since the early 1990s, there has been significant interest in the physics of dusty plasmas, which has now become a new discipline in plasma science. Dusty plasma exhibits new and unusual behaviour, and provides a possibility for modified or entirely new collective modes of oscillations, instabilities as well as coherent nonlinear structures.\\ First, a review of the important recurring terms -- The Cosmic Waves (CRs), the Alfven Waves (AWs), and the associated charged dust grains is presented. Starting from the basic composition of the CRs to their scattering mechanism, along with the different modes of scattering, is presented, along with the modes of confinement and a precise definition of each term. The paper also includes some useful diagrams and brief notes from the references. \\ Gravitation plays a significant role in the collapse of matter and the formation of cosmological structures. Unlike a static universe, this paper investigates the Jeans instability in a radiation-pressure-dominated expanding universe using the Einstein-de Sitter model for Euclidean geometry with zero curvature ($\kappa=0$). The fluid model for an expanding universe is constructed, and by taking small perturbations, the perturbed fluid equations are obtained. The dispersion relation of gravitational instability is derived using plane-wave solutions. In the static case (Newtonian cosmology), the classical Jeans instability criterion is revisited and modified in an expanding universe. The critical Jeans wave number of perturbations to excite Jeans instability depends upon the time-dependent expansion factor $S(t)$. It is found that short-wavelength perturbations are expected during the inflationary period of big bang cosmology, which are responsible for gravitational collapse and the formation of galaxies.

Long monodromic valleys arise in the penumbra of complex-structure moduli space. We show that their local branch data already determine whether they support controlled inflation, and thereby isolate the first controlled penumbral inflationary window. In the axion--saxion effective theory given in Eq.4, a branch-displacing odd term generates a plateau when $\Delta\equiv p+2\nu-q>0$, while covariant single-clock control further requires $p<2$, or $p=2$ with $12A_pm^2/V_0\gg1$ over the observational window. This splits penumbral valleys into no plateau, uncontrolled plateau, and controlled plateau before global completion is attempted. We identify a minimal analytic family with a closed-form valley and an invariant attractor equation for the full two-field dynamics, providing the first exactly solvable penumbral realization that remains predictive under the next penumbral order. The penumbra is thus promoted from a geometric suggestion to a predictive search principle.

Reliable uncertainty quantification is essential for the use of machine learning in physics, where scientific discoveries depend on validated probabilistic statements. We provide a structured overview of uncertainty quantification in ML for physics, introducing a unified taxonomy of uncertainty and clarifying the interpretation of predictive and inference uncertainties across frequentist and Bayesian frameworks. We discuss principled validation tools, including coverage, calibration, bias tests, and proper scoring rules, and illustrate them with simple regression and classification examples.