Papers for Friday, Mar 20 2026

We investigate whether low-frequency pulsar observations can provide LoS magnetic field estimates and whether these are consistent with synthetic LoS signatures extracted from a three-dimensional CME reconstruction constrained by Solar Orbiter data. We analyze a CME occultation of the LoS to PSR J1022+1001 on 20 August 2021, observed simultaneously with LOFAR and NenuFAR. From LOFAR, we derive time-resolved dispersion measure (DM) and rotation measure (RM) and isolate the CME contributions using background estimates for interstellar, solar wind and ionospheric components. We then infer the density-weighted LoS-averaged magnetic field component <B||>_PSR from the ratio delta-RM/delta-DM. In parallel, we reconstruct the CME using a semi-empirical 3DCORE model fitted to Solar Orbiter in-situ magnetic field observations at 0.65 au. We sample the modeled magnetic field along the pulsar LoS using fixed spatial sampling points and compute synthetic LoS-averaged signatures <B||>_3D for different flux rope configurations. The derived <B||>_PSR increases from approximately -9 nT to a peak near 63 nT during the observed interval. Comparison with synthetic signatures shows that the polarity and temporal evolution of the LoS signal are strongly dependent on the flux rope configuration and only a South-West-North (SWN) configuration (confirmed by Solar Orbiter in-situ data) reproduces the observed sign and overall evolution, whereas alternative configurations are incompatible. The modeled amplitudes, however, are systematically larger than the pulsar-derived values by roughly a factor of five. We show that simultaneous low-frequency pulsar DM and RM measurements can provide LoS magnetic field estimates for a CME and can be used to test CME magnetic structure against data-constrained three-dimensional reconstructions.

We simulate Venus' evolution with a coupled one-dimensional solar-atmosphere-lithosphere-mantle-core model to predict currently unobservable features and its eruptive mass flux. We identified four distinct evolutionary pathways that simultaneously match the atmospheric abundances of water and carbon dioxide as well as the lack of a core dynamo. These scenarios are characterized by I) generally monotonic cooling, II) a low mantle melt fraction in which Venus' volcanically active phase is ending, III) a small inner core, and IV) oscillations of internal properties. Through random forest classification we determined that the key parameters that distinguish these types are the initial mantle water abundance, the mantle viscosity, the dehydration stiffening strength, the eruption efficiency, and the melting point of the core. In each of the plausible histories, Venus retains at least one Earth ocean's worth of water in its mantle and remains volcanically active today. Venus' lack of a current geodynamo allows thermal histories with an initially large inner core in our parameter sweep. In 88% of plausible histories we found that Venus possessed a past magnetic field. The results strongly disfavor recent high eruption rate estimates, but are consistent with lower estimates. Current resurfacing estimates also strongly disfavor the low melt scenario, implying that Venus is not nearly volcanically ``dead.'' These predictions are testable with anticipated data and the model can be applied to exoplanets to predict their properties.

Extreme mass-ratio inspirals (EMRIs), comprising a stellar-mass compact object (CO) orbiting a supermassive black hole (BH), are key targets for future space-based gravitational-wave (GW) observatories. Incorporating the spin of the secondary body into waveform models not only enhances measurement precision but also offers insight into the spin distribution of stellar-mass objects. In this work, we construct the flux and waveform for an EMRI with a spinning secondary in a Kerr background under the linear-spin approximation. Using the radiative prescription (half-retarded minus half-advanced field), we derive orbit-averaged evolution equations for the fundamental constants of motion, including the energy, angular momentum, Carter-like constant, and the parallel spin component. This framework provides a tractable route to generating waveforms that incorporate the secondary spin, with the potential for further simplification in future work.

We report the discovery of bimodal structure in the drift rate distribution of upward-drifting burst clusters from the hyperactive repeating fast radio burst FRB 20240114A. Using unsupervised machine learning (UMAP dimensionality reduction combined with HDBSCAN density-based clustering) applied to 233 upward-drifting burst clusters from the FAST telescope dataset, we identify a distinct subpopulation of 45 burst clusters (Cluster C1) with mean drift rates 2.5x higher than typical upward-drifting burst clusters (245.6 vs 98.1 MHz/ms). Gaussian mixture modeling reveals strong evidence for bimodality (delta-BIC = 296.6), with clearly separated modes (Ashman's D = 2.70 > 2) and a statistically significant gap in the distribution (11.3 sigma). Crucially, we demonstrate that this bimodality persists when restricting the analysis to single-component (U1) burst clusters only (delta-BIC = 19.9, Ashman's D = 2.71), confirming that the result is not an artifact of combining single- and multi-component burst clusters with different drift rate definitions. The extreme-drift subpopulation also exhibits systematically lower peak frequencies (-7%), shorter durations (-29%), and distinct clustering in multi-dimensional feature space. These findings are suggestive of two spatially separated emission regions in the magnetosphere, each producing upward-drifting burst clusters with distinct physical characteristics, although confirmation requires observations from additional epochs and sources.

In this study, the behavior of high-frequency quasi-periodic oscillations (QPOs) is investigated around a Kerr black hole immersed in a uniform Bertotti-Robinson magnetic field. The motion of the test particle is analyzed by determining the geodesic equations and evaluating the corresponding orbital, radial, and vertical epicyclic frequencies. These fundamental frequencies are used to construct the theoretical framework of QPO models based on parametric and forced resonance mechanisms. Observational data obtained from several black hole X-ray binaries (GRO J1655-40, XTE J1550-564, XTE J1859+226, GRS 1915+105, H1743-322, M82~X-1, and Sgr~A$^{*}$) are used to constrain the black hole parameters through Bayesian inference and Markov Chain Monte Carlo (MCMC) analyses. For the X-ray binaries GRO J1655-40, GRS 1915+105, H1743-322, and M82~X-1, nonzero values of the dimensionless parameter $b=Bm$ are obtained at the $68\%$ confidence level within the framework of the parametric resonance model, while only upper bounds at the $90\%$ confidence level are obtained for the remaining sources. In contrast, in the case of the forced resonance model, only an upper bound at the $90\%$ confidence interval is obtained for the magnetic field parameter for all considered X-ray binary sources. The analysis indicates that the value of the magnetic field parameter is small but not negligible, producing minor modifications to particle dynamics and epicyclic frequencies. The influence of the magnetic field is further examined through the properties of the innermost stable circular orbit and the radiative properties of the thin accretion disk, including the energy flux and temperature profiles, within the allowed parameter range inferred from the MCMC analysis.

We present a simplified and general description of the high-redshift information in acoustic scale measurements from the cosmic microwave background and large-scale structure. The transverse distance interval between photon--baryon decoupling and a late epoch in the matter era provides an analytically tractable summary statistic thereof and a general diagnostic of the current tension between the Dark Energy Spectroscopic Instrument and the CMB. We show that this "matter-era distance excess" is unlikely to be explained by modified dynamics at low redshift. We then analytically derive the matter-era distance interval's sensitivity to new physics at high redshift, including nonstandard recombination, nonminimal dark matter dynamics, and spatial curvature; in particular, we explain how this observable represents a direct geometric measurement of (and underlies the current incompatibility with) neutrino masses. Finally, we demonstrate that phenomenological models of dynamical dark energy mediate the matter-era distance excess in a manner reliant on their unphysical, extrapolated behavior at high redshift. Invoking alternative explanations of the excess removes the CMB's contribution to the evidence for these models; the residual preference of around $1.7\sigma$ mostly derives from DESI's two lowest-redshift measurements of the Alcock--Paczynski distortion, without which it drops to $0.5 \sigma$.

Mass discrepancies in galaxies are empirically known to appear only below a characteristic acceleration scale a0. Here we show that this behaviour is not limited to galaxies: it extends continuously across the full hierarchy of self-gravitating stellar systems, from gas-rich dwarfs and spirals to massive early-type galaxies, and further down to compact stellar clusters. We introduce the Milgromian dynamics (MOND) depth index DM, together with dynamical maturity index T = tcross/tH, dynamical collisionality index T1 = tcross/trelax, with tcross being the crossing time, tH the Hubble time and trelax the median two-body relaxation time, and the MOND acceleration index A = abar/a0. We uncover a well-defined two-dimensional dividing surface in dynamical space. The "dark matter phenomenon" is found only in systems that are both in the deep-MOND regime (abar < a0) and collisionless (trelax > tH), while high-acceleration, collisional systems (abar > a0, trelax << tH), including globular clusters and UCDs, show no evidence for a mass discrepancy. This clean dynamical separation defines a new, physically motivated classification scheme for stellar systems, unifying galaxies and clusters under one framework. The observed division emerges naturally within the MOND framework and provides a useful diagnostic for examining how different gravitational paradigms account for the origin of the mass discrepancy.

We present the first spectropolarimetric time-series analysis of the Maunder Minimum analog HD 166620, using 12 nights of data from CFHT/SPIRou and a single epoch from CFHT/ESPaDOnS. While individual Stokes $V$ profiles exhibit no significant polarization signatures, we leverage the rotational coverage of the SPIRou dataset to compute a grand average LSD profile. Forward modeling of the cumulative Stokes $V$ signal, assuming a purely axisymmetric dipole, yields a best-fit dipole field strength of $B_{\rm dip} = 1.10^{+0.95}_{-0.90}$G ($3\sigma$). This field strength matches simulations of the solar dipole during the Maunder Minimum phase. Our results are consistent with independent constraints on the dipole field strength from an LBT/PEPSI snapshot and exclude the presence of strong non-axisymmetric fields potentially missed by this single-epoch observation. These findings provide direct empirical evidence that the transition to weakened magnetic braking involves a weakening of the large-scale magnetic field and suggest that HD 166620 represents a state comparable to the Sun near the peak activity of a grand minimum.

A recent development in the study of magnetic cataclysmic variable stars (mCVs) has been the identification of asynchronously spinning mCVs with orbital periods <2 h that have significantly higher white dwarf spin-to-orbital period ratios than their longer-period counterparts. We report the discovery of two additional mCVs in this class. The first, Gaia21akb, is a candidate asynchronous polar at the period minimum. While TESS photometry cannot, in isolation, lead to a conclusive identification of the orbital period, the probable orbital period of 1.29 h would be the second-shortest of any known polar and would result in a spin-to-orbit ratio of 0.9879. The second system in our study, ZTF18aazmehw, is an eclipsing mCV with a 1.50 h orbital period and a spin-to-orbit ratio of 0.867. Contrary to expectations for an asynchronous polar, ZTF18aazmehw does not show discernible evidence of pole switching and might possess a disk-like structure. The increasing number of short-period asynchronous mCVs with large spin-to-orbit ratios lends credence to theoretical predictions that asynchronously rotating mCVs with sufficiently strong white dwarf magnetic fields can achieve synchronization when their orbital separations have shrunk sufficiently.

In this Letter, we establish a robust, physically motivated classification method using a Support Vector Machine (SVM) trained on a "gold-standard" sample of 24 GRBs with spectroscopically confirmed progenitors (associated SNe or KNe). By isolating the prompt main spike to excise contamination from extended emission, we derive a quantitative classification index, I_SVM = 5.01 log_10 E_p,i - 1.25 log_10 E_iso - 0.34 log_10 T_90,z - 12.90 (units: keV, 10^52 erg, s). Events with I_SVM > 0 are classified as mergers. Analysis of the standardized classification weights reveals that the discriminative power of E_p,i is approximately 5 times that of T_90,z, while E_iso contributes a weight comparable to E_p,i. This quantitatively demonstrates that spectral hardness and energetics, rather than duration, are the primary physical signatures distinguishing mergers from collapsars. The derived boundary implies a stringent hardness ceiling for collapsars, while mergers are identified as outliers with excessive hardness relative to their energy budget. The classifier successfully identifies the nature of historic test cases, including the ultra-long GRB 111209A (collapsar) and the short GRB 050709 (merger), independent of instrumental eras. This tool paves the way for cleaning archival and future high-redshift GRB samples for precision cosmology.

The ~ 200,000 stars observed by the Kepler mission have provided unprecedented constraints across astrophysics. With the advent of modern spectroscopic and photometric surveys, new limits in stellar characterizations are within reach. In this work, we report a compilation of atmospheric parameters (Teff, logg, and [M/H]) for the Kepler stars by crossmatching with several spectroscopic and spectro-photometric surveys. We use these to calculate bolometric corrections, which combined with color-magnitude diagram (CMD) information from Gaia yield self-consistent luminosities on a survey-by-survey basis. These properties will aid in future explorations of Kepler data towards new astrophysical insights. We make our catalog publicly available online in Zenodo (doi:https://doi.org/10.5281/zenodo.18620911).

We present the X-ray analysis of coronal properties in a statistically representative sample of 23 mostly radio-quiet AGN from the SUBWAYS campaign (SUpermassive Black holes Winds in XrAYs), focusing on quasars at redshifts $0.1 < z < 0.4 $ and bolometric luminosities $2 \times 10^{44} <L_{bol}(erg/s) < 2 \times 10^{46}$. The main aim of this work is to investigate the properties of the hot corona through the study of the hard X-ray band emission, including a proper treatment of the soft X-ray band. High-quality X-ray spectra from XMM-Newton, complemented by NuSTAR data extending up to 30-40 keV in the rest frame, are available for this sample. The soft X-ray band (0.3-2 keV) spectrum is best fitted by a warm corona model with a median temperature of 0.40 keV, and an optical depth in the range $\tau$=20 - 40, consistent with previous results on lower luminosity sources. The hard X-ray band is well described using a hot corona model, with a median high-energy cut-off of 87 keV, at the lower end of the distribution of typical values found in Seyfert galaxies (100 - 200 keV). The derived median value of the optical depth ($\tau$ = 1 - 5) suggests the presence of a moderately optically thick corona. Combining the SUBWAYS results with literature samples at low and high redshift, we assemble the largest sample to date of AGN with E$_{cut}$ and accretion parameter measurements, finding a significant anticorrelation of E$_{cut}$ with both $\lambda_{Edd}$ and $L_{bol}$ with the median E$_{cut}$ decreasing from 250 - 300 keV at low accretion rates and luminosities to 90 - 100 keV at high accretion rates and luminosities - consistent with enhanced coronal cooling, possibly driven by pair-production. These results favor cooler, optically thicker coronae in luminous AGN compared to those in lower-luminosity Seyfert galaxies.

The icy material within protoplanetary disks plays a central role in planet formation, yet remains poorly characterized by observations. We present 1.6-28$\mu$m spectra of five disks obtained as part of the JWST Edge-on Disk Ice (JEDIce) program, representing the largest survey of disk ices to date. The major ice species H$_2$O, CO$_2$, and CO are detected towards all disks, and exhibit a wide range of absolute optical depths and optical depth ratios across the sample. This is suggestive of a range of ice abundances and compositions, but quantitative constraints will require radiative transfer modeling. All disks exhibit ice features across the entire spatial region where the IR continuum is detected; vertically elevated ice grains therefore seem to be ubiquitous in disks. The CO ice is consistently dominated by apolar CO:CO$_2$ mixtures, implying that the disk ice compositions are neither completely reset nor pristinely inherited from the protostellar stage. The presence of these mixtures also suggests that entrapment may be important in shaping the spatial distribution of CO within the disks. Small molecules commonly seen in protostellar ices (CH$_4$, CH$_3$OH, NH$_3$) are generally not detected in our sample, though tracers of ammonium salts (OCN$^-$ and the 6.85 $\mu$m band) are common, potentially reflecting an evolution towards comet-like ice compositions. The spectra also contain a wealth of information about the micron-sized dust, atomic and molecular gas, and PAH content, which together with the ice constraints will provide a comprehensive picture of the chemical, physical, and dynamical state of these systems.

In this work, we introduce TUNeS (Temporal UNet emulator for Structure formation), a neural network framework for accelerating N-body simulations by predicting the nonlinear evolution of the matter density field from an initial particle distribution. TUNeS employs a two-stage modeling strategy, combining particle-based inference with a density-field refinement on a regular grid, enabling accurate reconstruction of both large- and small-scale structures. The model is designed to operate across redshift, taking particle snapshots at arbitrary input redshifts and predicting density fields at arbitrary target redshifts. In this work, we evaluate its performance using simulations initialized at $z=100$, with predictions generated at multiple lower redshifts. Trained on only eight N-body simulations, TUNeS reproduces reference results with good agreement in both Gaussian and non-Gaussian statistics, including two-point correlations, one-point distributions, peak counts, and three-dimensional Minkowski functionals. In particular, at $k \simeq 1\,h\,\mathrm{Mpc}^{-1}$, the power spectrum error remains at the few-percent level. End-to-end inference from $256^3$ particles to a $256^3$ density grid can be completed in $\sim25\,\mathrm{second}$ on a single GPU. Thanks to its architectural design, the model naturally scales to larger particle numbers and larger volumes through particle batching and window-based refinement.

The emergence of habitable conditions on the early Earth and on rocky exoplanets requires persistent energy sources that can drive both prebiotic chemistry and climate warming under magnetically active young G to M stars. To quantify the contribution of stellar energetic particle (StEP) events associated with superflares to the atmospheric chemistry of young planets with primitive atmospheres, we carried out a suite of laboratory proton irradiation experiments on mildly reduced gas mixtures. We present first proton irradiation experiments of $N_{\mathrm{2}}O$/$CO_{\mathrm{2}}$ rich gas mixtures that yield abundant nitrous oxide ($N_{\mathrm{2}}O$) at mixing ratios up to 1000 ppmv, together with amino acid precursors including glycine, corresponding to global production rates of order $2\times10^{10}$ kg/yr on the early Earth. Our photochemical modeling of StEP driven proton irradiation reproduces the experimentally inferred $N_{\mathrm{2}}O$ production rates and provides self-consistent atmospheric $N_{\mathrm{2}}O$ profiles. We then use these profiles of $N_{\mathrm{2}}O$ as input to a 3D global climate model to evaluate the radiative and climatic impact of StEP generated $N_{\mathrm{2}}O$ in primitive atmospheres representative of the early Earth and young rocky exoplanets. Our results show that frequent StEP events can help alleviate the faint young Sun paradox on the early Earth and can maintain temperate surface conditions on young rocky exoplanets beyond the outer edges of habitable zone, while simultaneously enhancing the buildup of prebiotic molecules. Together, these processes may constitute a robust pathway toward early planetary habitability.

The TeV neutrino signal observed by IceCube from the active galactic nucleus (AGN) NGC 1068 can probe its innermost coronal regions. If these neutrinos originate from hadrons accelerated within a magnetized turbulent corona, their intensity and spectrum depend on the turbulent magnetic field strength and turbulence coherence scale. The gamma rays accompanying neutrino production are absorbed in this optically thick environment, in a way that depends sensitively on the size of the corona. By a joint fit of the IceCube and Fermi-LAT observations, we translate the multimessenger signal from NGC 1068 and the tentative signal from NGC 7469 into quantitative constraints on coronal properties. NGC 1068, with a significant TeV neutrino excess, favors a compact, strongly magnetized corona with a large turbulence coherence length relative to the coronal size. NGC 7469, with two $\sim 100$ TeV neutrino events, points instead to a somewhat larger corona with much smaller coherence length and high magnetization, but a very small fraction of energy in non-thermal protons. We obtain the diffuse flux from a population of Seyfert galaxies identical to either NGC 1068 or NGC 7469. Finally, we consider a third scenario, motivated by the spectral break observed in the diffuse neutrino flux at tens of TeV, with coronal properties intermediate between the two point-source-inspired models. To enable detailed comparisons with the IceCube and electromagnetic observations, we release our model predictions in a GitHub repository.

Testing the $\Lambda$CDM model requires cosmological probes spanning the wide redshift interval between Type Ia Supernovae (SNe Ia, $z\lesssim2.9$) and the Cosmic Microwave Background (CMB, $z\approx1100$). Gamma-Ray Bursts (GRBs), observed up to redshift $z=9.2$, offer the opportunity to explore this regime. Here, we investigate how many GRBs are needed to become a useful cosmological probe capable of independently testing deviations from $\Lambda$CDM suggested by the recent DESI BAO observations. We develop forecasts based on the two-dimensional X-ray and optical Dainotti relations, between the luminosity at the end of the plateau phase and its rest-frame duration. Using simulated GRB samples constructed from the observed population, we evaluate the constraining power of GRBs on cosmological parameters within the $w$CDM and $w_0w_a$CDM models, both independently and in combination with CMB observations. Our results show that GRB samples containing several tens to hundreds of well-characterized plateau can already approach the precision currently achieved by CMB measurements on the Dark Energy (DE) equation-of-state parameter $w$. Particularly, a sample of $\sim66$ optical GRBs can reach a precision $\sigma_w \approx 0.47$, comparable to that obtained from Planck within the $w$CDM framework. Such sample sizes are already attainable through Machine Learning techniques that double the number of GRBs using inferred redshifts. These forecasts indicate that future GRB observations, when combined with next-generation transient missions and improved statistical techniques, will provide an independent high-redshift probe of cosmic expansion and will play an important role in testing the robustness of potential Dynamical DE signals suggested by other cosmological datasets.

A particularly intriguing subclass of rocky exoplanets are the ultra-short period (USP) worlds that orbit their host stars in less than a day. These planets are particularly rare around M dwarf stars, with so far only ten that have a constrained mass and radius. We present the validation and characterization of the ultra-short period (0.3-days), Earth-sized planet TOI-4552b orbiting a nearby (27.26-pc away) M4.5V dwarf. Complementing the TESS photometry, ground-based transit observations from LCO, ExTrA and SPECULOOS validated the planetary radius and cleared the field of any contaminants. Speckle imaging with Zorro (Gemini-S) rules out false positive scenarios caused by eclipsing binary sources. Spectroscopic observations with NIRPS and HARPS were used to obtain stellar abundances, constrain the planetary mass, and, in conjunction with the transit observations, estimate the orbital parameters. TOI-4552 is a quiet star exhibiting no short-term stellar variations seen in photometric or radial velocity data that can be associated to stellar rotation. TOI-4552b ($M_p=1.83\pm0.47\,M_e$, $R_p=1.11\pm0.04\,R_e$) lies between the Earth-like and iron-rich composition tracks on the Mass-Radius diagram. The EXOPIE interior structure model, without constraints from refractory abundance ratio, yields a core mass fraction (CMF) of 0.54 and a bulk density of 7.74g/cm$^3$. Since the CMF spans a wide range due to the large uncertainty on the mass, the definitive interior composition cannot be determined with the current dataset. TOI-4552b hints as being marginally more iron-rich compared to the Earth but confirmation of its status requires additional, precise radial velocity measurements. Combined with its high emission spectroscopic metric (ESM=19.5), negligible stellar activity and short orbital period, TOI-4552b emerges as a compelling target for atmospheric and surface composition studies with JWST.

We study the non-Gaussianity of the large-scale clustering of high-redshift halos, seeking to assess which terms of standard bias expansions are needed to understand these highly biased populations. We find that the clustering can be well modeled with only linear and quadratic bias parameters while assuming a Gaussian underlying matter field. Our analysis focuses on AbacusSummit halos at redshift $z=8$. We work with halos of mass at least $1\times10^{11}h^{-1}M_\odot$ in boxes of side length $2h^{-1}$Gpc. Measurements of bias coefficients are made by fitting bias expansions to the halo power spectrum and bispectrum. Tidal bias is not detected with only a ~$0.1\sigma$ deviation from $0$, but we see a $17\sigma$ level detection for a bias term of the form $\delta^2$. A bias term of the form $\delta^3$ is weakly detected at the $1.3\sigma$ level. Nonlinear matter is also detected at a $1.3\sigma$ level. To test how bias evolves, we run one test at $z=5$. We use a mass threshold for halos that gives the same variance in the halo field as our $z=8$ sample. Bias is smaller at $z=5$ and a tidal bias is detected at the $1\sigma$ level. Bias coefficients at $z=5$ match a linear evolution of the $z=8$ bias coefficients to within $10\%$.

Nuclear Star Clusters (NSCs) are dense stellar systems located at the centers of galaxies. Employing Enzo-Abyss, which integrates hydrodynamics with a direct N-body solver, we introduce a simulation capable of resolving the evolution of NSCs within a live galaxy. This includes live dark matter, gaseous dynamics, star formation and feedback, collisional dynamics for star clusters. The evolution of NSCs is typically shaped by two main processes: mergers of star clusters and in-situ star formation. Our simulation enables investigation of the contributions of these mechanisms to the growth of NSCs. This work focuses on the impact of stellar physics and gas content on the growth of NSCs within a dwarf galaxy. To this end, we carry out four simulations, a fiducial simulation, one without supernova feedback, one with low star formation efficiency, and one with higher galactic gas content. This study shows a likelihood that both mergers and in-situ star formation contribute to NSC evolution comparably. In addition, mergers result in disruption of dense gas clumps within star clusters, indicating that in-situ star formation is suppressed when mergers occur. However, the limitations -- such as the lack of individual star physics and limited spatial/particle mass resolution -- hinder drawing a definite conclusion. Nevertheless, with further development, our simulations will serve as a cornerstone that untangles the complex interplay between mergers and in-situ star formation in shaping the structure and mass of NSCs, thereby providing insights into their formation and evolution.

High-energy astrophysical events, particularly Gamma Ray Bursts (GRBs), have been proposed as significant contributors to mass extinction events on Earth-like planets in most of the galaxy, internal to our radius in it. This paper examines the extent to which GRBs may reset the evolutionary progress of complex life through repeated extinction-level disruptions. While resilient extremophiles may survive even the most intense GRBs, more complex surface-dwelling organisms are vulnerable to indirect atmospheric effects, primarily UV exposure following ozone depletion. By identifying evolutionary milestones and estimating how frequently GRBs would need to occur to prevent recovery between such milestones, this work proposes that GRBs could act as evolutionary filters, limiting the emergence of advanced life, but only much closer to the galactic center. We consider the implications for searches of various biosignatures versus technosignatures.

Changing-look active galactic nuclei (CL-AGN) exhibit spectroscopic and photometric changes on timescales of months to years, making them powerful laboratories for studying accretion variability onto supermassive black holes. Motivated by the growing relevance of large spectro-photometric time-domain surveys, especially the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), we compiled a master catalog of known CL-AGN from the literature and evaluated its spatial overlap with the Rubin survey footprint. Using a geometric cross-match based on sky coordinates, we identify 79 sources located in high-cadence regions of the main survey footprint (Wide-Fast-Deep, or WFD), including 5 particularly favorable targets lying within the Deep Drilling Fields (DDFs) of COSMOS and XMM-LSS. These sources represent especially promising candidates for future variability studies in the Rubin era. This Research Note presents a first proof of concept for connecting known CL-AGN with Rubin observing fields, while the full catalog and a more comprehensive analysis will be presented in a forthcoming paper.

We present Active Learning for Accelerated Bayesian Inference (\texttt{alabi}): an open-source Python package for performing Bayesian inference with computationally expensive models. Given a forward model and observational data to construct a likelihood and priors, \texttt{alabi}\ uses a Gaussian Process (GP) surrogate model trained to predict posterior probability as a function of input parameters, and employs active learning to iteratively improve GP predictive performance in high-likelihood regions where the GP is most uncertain. \texttt{alabi}\ provides a uniform interface for using Markov chain Monte Carlo (MCMC) with different packages, including the affine-invariant sampler \texttt{emcee}, and nested samplers \texttt{dynesty}, \texttt{multinest}, and \texttt{ultranest}. This approach facilitates accurate estimation of the desired posterior distribution, while reducing the number of computationally expensive model evaluations required by factors of thousands. We demonstrate the performance of \texttt{alabi}\ on a variety of test cases, including where inference is challenging due to complex posterior structure or high dimensionality. We show that \texttt{alabi}\ offers a substantial improvement for likelihood functions with evaluation times $\gtrsim 1$\,s, speeding up MCMC computations by a factor of $10-1000\times$ when tested on problems with up to 64 dimensions.

We identify different dust features in our compilation of infrared spectra for 267 planetary nebulae (PNe) from the Spitzer, ISO, and IRAS telescopes. We classify 206 objects according to their dust type: mixed dust (MD), oxygen-rich dust (ORD), carbon-rich dust (CRD), PNe with only polycyclic aromatic hydrocarbons (PAHs) in their spectra (oPAH), and featureless (F) PNe. We study statistically the distributions of surface brightness and diameter of PNe with different types of dust as well as their distributions in the this http URL find that both MD and ORD PNe are closer to the Galactic centre than CRD and oPAH PNe, and that the Galactic distributions of each pair of groups are statistically compatible, suggesting that they have similar progenitors. Since oPAH PNe have, on average, larger diameters and lower surface brightness than CRD PNe, we suggest that oPAH PNe are evolved CRD PNe. On the other hand, F PNe have the lowest surface brightness and the largest diameters, suggesting they could contain evolved PNe from any initial type of dust. Among the PNe with silicates, we find that only a few ORD PNe have just amorphous silicates in their spectra, and their distributions of Galactocentric distances and Galactic heights suggest that they had low-mass progenitors. MD PNe with both amorphous and crystalline silicates have the largest surface brightness and the smallest diameters and might be the earliest stages of PNe with the most massive and metal-rich progenitors.

We investigate the evolutionary pathways of green valley (GV) galaxies drawn from the SDSS-IV/MaNGA survey. The GV sample is divided into fast- and slow-rotating galaxies based on stellar spin, and their stellar and gas-phase metallicities are compared. Fast-rotating galaxies exhibit systematically higher metallicities than slow-rotating galaxies in both gas and stars. However, the gas-phase difference is significant only at low stellar masses, while the stellar metallicity offset persists across the full mass range. Using a simple yet physically motivated chemical evolution model, optimised to jointly fit gas-phase metallicities and integrated stellar spectra, we reconstruct the star formation and chemical enrichment histories of individual galaxies and constrain gas inflow and outflow parameters. At low stellar masses, fast- and slow-rotating galaxies show similar gas-infall and star formation timescales, but the the slower population experienced stronger outflows which reduce their chemical content in both gas and stars. At high masses, the combination of reduced pristine gas inflow and more efficient gas removal in slow-rotating galaxies produce gas-phase metallicities comparable to fast-rotating galaxies but systematically lower stellar metallicities. These differences suggest distinct evolutionary pathways for GV galaxies. Slow-rotating galaxies likely experienced more mergers, usually associated with strong gas removal processes, leading to their systematically lower metallicities. At low masses, stronger supernova-driven outflows reduce their chemical content while leaving star-formation timescales similar to fast-rotating galaxies. At high masses, merger-triggered AGN feedback may rapidly deplete and suppress gas infall, producing the shorter star-formation timescales seen in slow-rotating galaxies. Alternative environmental and assembly-driven scenarios are also discussed.

We present the first atmospheric observation and characterisation of the aligned, 1468 K hot Jupiter, NGTS-2b, with one JWST NIRSpec/G395H transit. These observations complete the GO 3838 observing campaign of the BOWIE-ALIGN program, which aims to investigate the link between hot Jupiter atmospheric composition and formation history through the atmospheric analysis of planets orbiting F stars that are aligned and misaligned with the host stellar spin axis. The 2.84-5.18 micron spectrum shows weak absorption features attributed to H$_2$O and CO$_2$ absorption, which our free chemistry retrievals fit with posteriors that converge on high mean molecular weight solutions attained through significant H$_2$O mixing ratios. By comparing our results to interior modelling, we show that some of these solutions exceed the 43.5x solar upper limit we obtained from our interior structure models. Such solutions are likely due to cloud-metallicity degeneracies and insufficient wavelength coverage to resolve them. We show that, in the case of our observations, the likelihood distribution of H$_2$O abundances is flat and uninformative, such that our retrievals are biased by the prior. Additionally, our statistically favoured atmospheric solution contains absorption from SO. The chemical abundances retrieved with this model are likely not astrophysically feasible and we demonstrate that the presence of SO is driven by only two data points. Our equilibrium chemistry retrievals hint at a subsolar C/O ratio and supersolar metallicity; however, we find wide posterior distributions that extend to solar values.

The release of heat by a planetary embryo modifies the local density perturbations, forming thermal lobes in its vicinity, and thereby altering the torque exerted by the disk on the embryo. In laminar disks, these thermal torques can dominate the disk-embryo interaction, rendering the classical Lindblad and corotation torques largely subdominant. The aim of this work is to investigate how turbulence driven by the MRI instability affects the thermal lobes formed around a planetary embryo, and to analyze the resulting torque acting on the embryo. We evaluate the thermal torques exerted on a planetary embryo of mass $M_p=0.33M_{Mars}$ and on a planetary core with mass $M_{p}=1M_{\oplus}$, each embedded in a turbulent gaseous protoplanetary disk, by means of high-resolution 3D magnetohydrodynamics simulations that include thermal diffusion and an initially toroidal magnetic field. The magnetic field strength is characterized by the $\beta$-plasma parameter with $\beta\in\{50,1000\}$. We consider two values for the luminosity of the planetary embryo: $L=0$ (cold thermal lobes) and $L=L_c$ (hot thermal lobes), where $L_c$ represents the critical luminosity. We find that, even in the presence of a weak magnetic field and irrespective of the luminosity, for both planetary masses, the development of turbulence in the disk (which takes between 1.5 to 3 orbital periods) completely disrupts the thermal lobes. As a result, the torque acting on both the planetary embryo and the Earth-mass core displays a strongly oscillatory behavior. This suggests that planets with masses in the range $0.03M_{\oplus}\lesssim M_{p}\lesssim 1M_{\oplus}$ experience stochastic migration, as expected in turbulent disks. Thermal torques become inefficient in turbulent regions of protoplanetary disks, such as outside the dead zone, in both the inner and outer disk regions where the magnetorotational instability operates.

The information recoverable from galaxy spectra depends fundamentally on spectral resolution, yet assembling large samples at high resolution remains observationally expensive. We present a deep-learning framework for spectral super-resolution that enhances low-resolution galaxy spectra by a factor of $\sim$10 in resolving power ($R\sim100$ to $R\sim1000$). The model is trained on 1,187 paired JWST/NIRSpec observations from the JADES program, where low-resolution prism spectra are matched with medium-resolution grating spectra (G140M, G235M, G395M) combined into a unified reference covering 1-5 $\mu$m. Our three-stage architecture performs an initial super-resolution, infers the redshift from the coarse reconstruction, and then applies a physics-informed residual refinement that uses attention across emission-line tokens to learn inter-line relationships and predict parametric line profiles, alongside a convolutional branch for continuum corrections. Evaluated on a 20% held-out sample, the model achieves noise-limited residuals over most of the spectral range and systematically improves the signal-to-noise ratio of key diagnostic lines including [OII], H$\beta$, [OIII], and H$\alpha$, often by factors of several. The super-resolved spectra successfully deblend features that are entirely unresolved at prism resolution, such as the [OIII] $\lambda\lambda4959,5007$ doublet and H$\beta$. As a proof of concept using JWST data, this approach is readily extensible to the low-resolution grism spectroscopy that will be delivered by Euclid and the Roman Space Telescope, potentially enabling population-level diagnostics across millions of galaxy spectra that would otherwise be inaccessible at grism resolution.

We present the weak gravitational lensing dataset from the Super-pressure Balloon-Borne Imaging Telescope (SuperBIT), which imaged 30 galaxy clusters during its 45 night flight in April to May 2023. SuperBIT is a first-of-its-kind balloon-borne imaging telescope that achieved near diffraction-limited observations in near-space conditions above 98% of the Earth's atmosphere. We use the metacalibration algorithm to obtain calibrated galaxy shapes for our target clusters and several calibration fields, enabling unbiased reconstruction of the weak-lensing signal. We employ several diagnostics throughout the pipeline, including assessments of point-spread function (PSF) modeling residuals and their impact on weak-lensing measurements, as well as tests for correlations between galaxy shapes and measured galaxy and PSF properties. To assess the multiplicative shear bias of the pipeline, we analyze a parallel set of simulated images that incorporate the real observing conditions from the flight, including measured SuperBIT PSFs, observed sky backgrounds, and detector noise, yielding a bias of $(1.1 \pm 7.8)$~per~cent.

The galaxy-halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a non-parametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past lightcones constructed with the Santa Cruz semi-analytic models. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the non-parametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modeling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

Weak gravitational lensing is a widely used probe in cosmological analysis. It allows astrophysists to understand the content and evolution of the Universe. We are entering an era where we are not limited by the data volume but by systematic uncertainties. It is in this context that we present here a simple python-based software package to help in the computation of E-/B-mode decomposition, which can be use for systematic checks or science analysis. As we demonstrate, our implementation has both the high precision and speed required to perform this kind of analysis while avoiding a scenario wherein either numerical precision or computational time is a significant limiting factor.

More than 200 moons exist in our Solar System, yet no exomoon has been confirmed to date. While the innermost two planets of the Solar System lack natural satellites and most studies favour the existence of exomoons around long-period planets, some theoretical studies that take tidal dissipation, orbital decay, and migration processes into account suggest that exomoons may survive around short-period exoplanets. We investigated the impact of exomoons on planetary thermal phase curves and assessed their detectability within a theoretical framework. We simulated the thermal phase curves of exomoon-exoplanet systems, including mutual transits and occultations, and explored their dependence on planetary orbital periods across a wide range of systems. Close-in airless exomoons maintain large day-night temperature contrasts, amplifying the thermal phase-curve signal of the system. When the exomoon transits or is occulted by the exoplanet, the transit depth varies with the planetary phase, and the occultation depth varies with the exomoon's phase. The maximum occultation depth can reach $\sim$ 20 ppm for long-period systems. For short-period planets, the signal can reach up to $\sim$100 ppm, although such configurations may not be dynamically stable over long timescales. If exomoons are not accounted for, the planetary temperature distribution retrieved from observed thermal phase curves may overestimate the planetary day-night temperature contrast and underestimate the planetary horizontal heat transport. In principle, the periodic exomoon-exoplanet mutual occultation signal could be extracted using methods such as box-fitting least squares, providing a framework for future observational studies and instrument planning.

The physical origin of fast radio bursts (FRBs) remains an unsolved mystery in astrophysics, with the magnetar central engine model as the leading framework. Systematically searching for physical associations between FRBs and the energetic astrophysical transients (ATs) that form magnetars provides a critical test of this scenario, and key clues to FRB progenitors. We perform a systematic search for FRB-AT associations using a sample of 3765 unique FRBs, combining the second CHIME/FRB catalog with 124 additional localized FRBs with measured redshifts. We develop a 3D Bayesian inference framework that jointly incorporates angular separation, positional uncertainty, and redshift constraints to quantify the association probability of candidate pairs. Through spatial cross-matching, we identify 14 FRB-optical transient and 15 FRB-gamma-ray burst (GRB) candidate pairs. Our framework recovers the previously reported high-significance association between FRB 20180916B and AT 2020hur, with an association probability of 0.9998. For the proposed candidate FRB 20190309A and short GRB 060502B, our analysis yields an association probability of 0.83, which is insufficient to claim statistically significant association. No new statistically significant FRB-AT associations are found for all remaining candidates. Our work demonstrates that small angular separation alone is insufficient to confirm FRB-AT associations, and high-precision FRB localization is essential for definitive identification.

We analyze a clean sample of 1,118 late-type, face-on galaxies without AGN contamination from the MaNGA survey. Their photometric structures are quantified via two-component (bulge+disk) decompositions on deep $g$-band images from the DESI Legacy Survey. Using a disk central surface brightness of $\mu_{\rm 0,d,cor}$(g) = 22 $\pm$ 0.3 mag arcsec$^{-2}$ (corrected for inclination and cosmic dimming) as the classification threshold, we identify 159 low surface brightness (LSB) galaxies, 388 LSB candidates, and 571 high surface brightness (HSB) galaxies. LSB galaxies are predominantly low-mass ($M_\ast < 3 \times 10^{10}$ M$_\odot$), exhibiting 29\% larger effective radii, 15\% lower star formation rates (SFRs), and 12\% reduced gas-phase metallicities than HSB counterparts at comparable masses. These differences cause systematic offsets from standard scaling relations. Despite comparable gas content, LSB galaxies host older stellar populations, longer gas depletion times, and less efficient star formation. Spatially resolved analyses further reveal that LSB galaxies display centrally suppressed $\Sigma_{\rm SFR}$, flatter SFR gradients, and rising specific SFR profiles toward their outskirts. Together with steeper negative metallicity gradients, these trends suggest ongoing gas accretion fueling outer-disk star formation. Consistently, the outer regions of LSB galaxies exhibit stronger H$\delta_A$ absorption and lower D$_n$4000 indices, indicating fading A-star populations. Moreover, LSB galaxies show lower $\Sigma_{\ast}$ across all $R/R_e$ and more centrally depleted stellar mass profiles on an absolute radial scale, compared with HSB and large-size star-forming galaxies. Collectively, LSB galaxies represent a distinct population with slow evolution, inefficient star formation, and continued susceptibility to late-time gas accretion and peripheral star formation.

We present high spatial resolution ($\lesssim$1.0''), multi-wavelength observations of UGC 2369S, a nearby luminous infrared galaxy showing three distinct cores separated on kpc scales in near-infrared (NIR) imaging with significant X-ray emission. Utilizing optical/NIR adaptive optics (AO), radio, \chandra X-ray, as well as archival HST imaging, we perform a comprehensive study of AGN activity, obscuration, and host properties. As one of the clearest cases of a triple-nucleus merger at $\simeq$3 kpc separations, UGC 2369S is the first to be studied with high-resolution observations at multiple wavelength. We find that the northern core, having possibly the most massive black hole in the system ($\rm M_{BH}\simeq10^{8}\,M_{\odot}$) is consistent with a heavily obscured AGN. However, its high dust extinction ($\rm A_v>5$), hydrogen column density ($N_\mathrm{H}\gtrsim 10^{25}\,\rm cm^{-2}$) and non-detection of optical coronal lines and coronal X-ray emission leave the identification inconclusive. The other two cores show no evidence for black-hole activity and instead exhibit signatures of tidal disruption. From stellar mass surface density and stellar velocity dispersion maps, we infer that the strongly varying gravitational potential in this three-body system may have cannibalized the stellar bulge of the southwestern core, leaving a metal enriched remnant. An ongoing survey focusing on similar triple systems could help us understand how they evolve and help benchmark numerical simulations, providing insight into gravitational wave predictions and the formation of the most massive black holes.

We report our optical observations of the Crab pulsar using the Imager of MPPC-based Optical photoN counter from Yamagata (IMONY), a high-time-resolution photon-counting imager with 100 ns timing resolution, mounted on the 3.8 m Seimei telescope in Japan (f/D~6). The detector format was upgraded from a $4\times4$ to an $8\times8$ GAPD array with larger pixels ($100$ to $200~{\mu m}$), resulting in a 14".5 field of view on the Seimei telescope. We conducted nightly optical observations for one week, including two nights of simultaneous optical and radio observations with the 64 m Usuda radio telescope. Thanks to the large diameter of the Seimei telescope and the high time resolution of IMONY, we successfully detected optical Single Pulses (SPs) emitted in each rotation. Moreover, we found an optical peak timing drift of $30\pm7.9~\mathrm{\mu s}$ over three days, with a significance of $3.9{\sigma}$. The corresponding emission region size is 9.1 km, which is equivalent to 0.006 times the light cylinder radius of the Crab pulsar. We ruled out the possibility of a pulsar glitch and suggested that the optical pulsed emission region of the Crab pulsar may fluctuate due to the spatial drift and variations in the magnetospheric caustics.

We present the first statistical census of emission-line variable active galactic nuclei (EVA) at cosmic noon by combining untargeted and deep HETDEX spectroscopy with multi-epoch spectra from SDSS, DESI, and LAMOST. Anchoring all candidates to a HETDEX spectroscopic epoch and requiring AGN classification in either the HETDEX or the external epoch(s), we identify a homogeneous sample of 100 EVA at z~1.5, including 98 newly identified. Emission-line variability is selected primarily through statistically significant line-flux changes, supplemented by extensive visual inspections using contemporaneous photometric light curves. The resulting incidence fraction is $f_{\rm EVA} \approx 0.9\%$. The rest-frame intervals between spectroscopic epochs span $\sim$1--10 yr, with brightening and dimming events exhibiting statistically indistinguishable characteristic timescales ($\Delta T\sim2.2$ and $\sim2.6$ yr, respectively). A key result is the characterization of the Baldwin effect in the time domain: while many EVA follow the ensemble Baldwin effect (eBeff) between two epochs, a substantial fraction exhibit apparent anti-eBeff responses. Time-resolved spectroscopy of an individual source reveals that the intrinsic EW--luminosity relation is non-stationary, with the line-to-continuum responsivity systematically evolving from stronger to weaker across successive variability cycles; sparse two-epoch sampling of this evolving intrinsic Baldwin evolution (iBeff) naturally produces both eBeff-like and anti-eBeff behaviors. Finally, EVA show no strong preference for extreme Eddington ratios but exhibit a mild tendency toward lower $\lambda_{\rm Edd}$ values relative to matched control samples, driven primarily by sources observed in their dim states. Together, these results establish a coherent framework for interpreting emission-line variability in AGN at the peak epoch of cosmic black hole growth.

As stars evolve, they undergo significant changes in their physical properties, which can have a profound impact on the planets orbiting them. In particular, the mass lost through stellar wind may be partially accreted by orbiting planets. We present the results of 18 simulations of one-planet systems with planetary masses of 0.5, 1, 2.5, 5, 10, and 13~$\mathrm{M}_\mathrm{J}$, each at initial orbital distances of 5, 10, and 20~AU, orbiting a 2~M$_\odot$ star through its red giant branch and thermally pulsating asymptotic giant branch phases. Our results show that planets with smaller orbits and higher masses accrete more stellar wind material than their wider-orbit and lower-mass counterparts, although the total mass accreted across all simulations remains small compared to their initial planetary mass. Even for the most massive planet, 13 $\mathrm{M}_\mathrm{J}$ at 5 AU, the total mass accreted was $\sim0.56$\% of the planet's initial mass; nevertheless, we find that the accretion luminosities of the simulated planets, with the exception of one planet, exceed their expected equilibrium luminosities, suggesting that such emission could be potentially detected. This result is key for the detection of planets around AGB stars, which have no confirmed detections as of yet. We also estimated the accretion and luminosities of two detected two-planet systems over a few orbits, obtaining results consistent with the one-planet simulated systems. Additional tests without wind accretion and with stellar wind drag force showed that, while both have a negligible effect on the orbital evolution, wind accretion remains relevant for the planetary luminosity.

Using a sample of 116 wide binary systems as coeval and chemically homogeneous stellar pairs, we investigate the factors governing lithium depletion in main-sequence stars. We recover the well-established morphology of the lithium--effective temperature ($T_{\mathrm{eff}}$) relation, including the Li dip (6200--6600\,K), the Li plateau (6000--6200\,K), and a linear trend for cooler stars ($T_{\mathrm{eff}}$ $<$ 6000\,K), where lithium abundance increases by $\sim$0.15\,dex per 100\,K. We demonstrate that the apparent correlation between projected rotational velocity ($v\sin i$) and lithium abundance is secondary to the underlying $T_{\mathrm{eff}}$ dependence; $v\sin i$ is not an independent driver of lithium depletion in our sample. Notably, we identify an anomalous system within the Li dip where the primary star exhibits a $\sim$1.4\,dex lithium excess compared to its secondary companion at nearly identical $T_{\mathrm{eff}}$. We discuss two plausible origins for this anomaly: external enrichment via planetesimal accretion or planetary engulfment, and binary interactions with an unresolved tertiary companion. Our results confirm $T_{\mathrm{eff}}$ as the dominant parameter controlling lithium depletion, while highlighting that additional, non-standard processes can occasionally produce significant lithium enrichment.

Comet 12P/Pons-Brooks is a periodic comet with an orbital period of approximately 71 years. Because of the period duration, aphelion of 35.3 au, and highly inclined orbit (74.2 deg), it is classified as a Halley-type comet. In this work, we present the results of the analysis of photometric and spectral observations of comet 12P/Pons-Brooks while it was at a heliocentric distance beyond 1.1 au, before perihelion passage. Quasisynchronous observations were carried out from July 2023 to March 2024 using 0.61-m and 1.3-m telescopes at the Skalnaté Pleso Observatory, 1-m telescope Zeiss-1000 Sanglokh International Astronomical Observatory, and 0.7-m telescope AZT8 at the observation station of Taras Shevchenko National University of Kyiv. Photometric observations were conducted using B, V, and R filters of the Johnson-Cousins and Bessel photometric systems. We completed our data collection with broadband photometric and spectral observations on the 2-m HCT telescope at the Indian Astronomical Observatory. During this period, the dust coma morphology and photometric data measured in the R filter, such as the apparent, absolute magnitudes and dust activity level, indicate several repeated outbursts in brightness and dust activity occurred around July 22, September 25, October 17-23, and a significant outburst was recorded on November 1-7. We detected strong gas emission features in the cometary spectrum belonging to CN, C2, and C3 molecules. Furthermore, gas production rates of these molecules were estimated using Haser's model. A comparison of dust activity level and CN gas production rate shows that dust activity increased significantly as the comet approached the Sun, while gas activity exhibited more moderate changes.

We present a feasibility study of the Antarctic 12m Terahertz Telescope (ATT12), a next-generation facility to be constructed at New Dome Fuji in Antarctica, designed to open up the FIR and THz windows for extragalactic astronomy. While ATT12 will enable a wide range of Galactic and extragalactic science, this paper focuses on its potential for studies of dusty star-forming galaxies (DSFGs) across cosmic time. Using realistic atmospheric transmission models and the planned instrumental specifications of heterodyne spectrometers and wide-field multi-color continuum cameras, we assess the expected sensitivity and scientific capabilities. We show that spectroscopic observations will enable detections of [CII]158um from galaxies with log(LIR/Lsun)>12 out to z~7, while [OIII]88um will remain observable for HyLIRG-class systems up to z~10. Line ratios including [OIII]52/88um, [NII]122/205um, and [OIII]/[NIII] will provide unique diagnostics of electron density and O/N abundance at z~4-8. Wide-field continuum surveys with the continuum cameras (KIDS-1/2; 300-850 GHz) will reach confusion-limited depths of ~1-2 mJy over ~10,000 deg$^2$, detecting of order $10^{6}$-$10^{7}$ DSFGs with log(LIR/Lsun)>12 at z<5 and $\lesssim10^{3}$--$10^{4}$ HyLIRGs up to z~7 or higher. Higher-frequency cameras (KIDS-3/4; >850 GHz) are designed for targeted follow-up observations and to extend coverage toward the THz regime. Taken together, ATT12 will provide the first statistically representative samples of DSFGs across cosmic time and, through synergy with ALMA, JWST, and the proposed FIR Probe PRIMA, will establish a multi-wavelength framework in which ATT12 discovers large samples through wide-area surveys, ALMA provides high-resolution follow-up of gas and ISM structure, JWST probes stellar populations and metallicity in the rest-frame optical/NIR, and PRIMA delivers ultra-sensitive FIR spectroscopy.

We leveraged reliable age and distance estimates from LAMOST-DR3 and APOGEE-DR17+AstroNN combined with \gaia\ data to perform a detailed analysis of the stellar age distribution in the Milky Way's (MW) outer disc using giant stars. Selecting stars near the midplane ($|z|<0.3$ kpc) on near-circular orbits ($\lambda_c > 0.9$), we analysed these independent datasets that employed different age-estimation methods. Our stringent kinematic selection criteria effectively exclude halo stars, ensuring that the observed age trends reflect genuine disc properties rather than contamination from older halo populations. Our results reveal a 'U-shaped' stellar age profile, where a negative gradient in the inner disc transitions to a positive gradient in the outer disc region. We identify the minimum in the stellar age profile at $R_{\rm min}=11.28 \pm 0.58$ kpc and $R_{\rm min}=12.15\pm 0.62$ kpc for the APOGEE-DR17 and LAMOST-DR3 samples, respectively. Using N-body+SPH simulations, we demonstrate that $R_{\rm min}$ corresponds to the break radius in the stellar density profile ($R_{\rm br}$), marking the edge of the Galaxy's star-forming disc. This break arises from a sharp decline in the star formation rate, with the outer positive age gradient produced by the radial migration of stars born inside $R_{\rm br}$. The cessation of star formation in the outer disc might be due to several mechanisms, including the dynamical influence of the bar's outer Lindblad resonance, the onset of the Galactic warp, or thermally regulated star formation. Overall, our results support the picture that the MW has a Type II (down-bending) stellar disc with a break at $R_{\rm br} \approx 11.28-12.15$ kpc, where the combination of star-formation cut-off and radial migration produces the observed U-shaped age profile.

There is now strong evidence that the Milky Way (MW) hosts a nuclear stellar disc (NSD). However, whether the NSD is purely axisymmetric or contains a nuclear bar remains unresolved. Since approximately $50\%$ of barred galaxies with MW-like mass in the local Universe host a nuclear bar, investigating whether the MW hosts one is of interest. We conduct a systematic analysis to identify robust kinematic diagnostics capable of determining whether the MW hosts a nuclear bar. Using N-body simulations, we explore the kinematic signatures indicative of a nuclear bar. Using the phase-space coordinates longitude $(\ell)$, latitude $(b)$, proper motions ($\mu_\ell$ and $\mu_{\rm b})$ and line-of-sight velocity $(v_{\rm los})$, we test various diagnostics assuming different nuclear bar orientations. We also evaluate how sample size, dust extinction and bar amplitude influence the efficacy of the diagnostics. We identify two independent kinematic diagnostics capable of revealing a nuclear bar in the MW: (1) the vertex deviation, $l_{\rm v}$, of the ($v_{\ell}-v_{\rm los}$) velocity ellipse; and (2) The asymmetry in the $\mu_{\ell}$ vs $\ell$ distribution. While both are impacted by the sample size and extinction, the vertex deviation proves more robust, especially when combining stars from multiple observational fields. We also assess the correlation between the line-of-sight velocity and the $h_3$ Gauss-Hermite moment ("skewness") of the line-of-sight velocity but find no clear distinction between an NSD and a nuclear bar based on this metric. Our results suggest that data from the current KMOS survey may allow a marginal detection of a nuclear bar using the vertex deviation method. A companion paper provides further validation and detailed analysis of this approach. Nonetheless, future surveys will provide the high quality data necessary to fully exploit the diagnostics outlined in this study.

Context. The High Resolution Telescope of the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI-HRT) operates in an extreme observational environment, observing the Sun as close as $0.28$ au. The high thermal load and large illuminating field puts high demands on the instrument in terms of both imaging performance and false light control. Aims. To characterise the amount of stray light (false light) within SO/PHI-HRT, apply a correction, and re-compare the data products with the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory (SDO/HMI). Methods. We analyse solar limb profiles and a Mercury transit to quantify the amount of stray light and add a correction term when partially reconstructing the SO/PHI-HRT images. For the comparison with SDO/HMI we use data from the 2023 March Solar Orbiter inferior conjunction and compare the magnetic fields on a pixel-by-pixel basis. Results. Increased continuum intensity contrast in the quiet Sun, and darker intensity levels are found in strong magnetic features. Consequently, much stronger fields are inferred in these features. Comparing the stray light corrected data with that from the standard SDO/HMI data products results in a much closer agreement across all vector magnetic field components, particularly when the cadence and noise levels are identical. In most solar features, SO/PHI-HRT infers stronger fields than the SDO/HMI line-of-sight magnetograms. Compared to the vector magnetic field from SDO/HMI the two are very well aligned, with only slight differences in the strongest field regions (where $|\mathbf{B}|>1600$ G or $|\mathbf{B_{\text{LOS}}}|>1300$ G).

Hydrogen cyanide (HCN) is crucial for the RNA World hypothesis, forming biomolecules essential for early life. Life likely emerged around 4 billion years ago during the early Archean Eon, a period on Earth with a fainter sun, frequent impacts, and a weakly reducing atmosphere. Warm little ponds (WLPs) are hypothetical protective aqueous environments that help explain the emergence and evolution of fragile prebiotic chemistry in such a hostile environment. WLPs need to undergo cycles of evaporation and rehydration, concentrating prebiotic molecules that increase the likelihood of (de-)polymerisation and forming early RNA molecules. We use a 1-D model of atmospheric chemistry to compare atmospheric HCN delivery to WLPs with exogenous sources. Using early Archean Earth as our baseline, we examine the sensitivity of atmospheric HCN delivery to the atmospheric C/O ratio, semi-major axis, assumed stellar host type, and methane budget, exploring conditions across rocky exoplanets. We find that atmospheric HCN delivery is sensitive to these parameters but its values generally exceed that of meteoritic delivery and our baseline Archean Earth. Planetary atmospheres with higher C/O ratios within the habitable zones of G stars and those closely orbiting M-dwarfs deliver the most atmospheric HCN. We find that atmospheric HCN delivery is remarkably robust, so this molecule is likely not the rate limiting step for the emergence of prebiotic chemistry on rocky exoplanets. This finding, with important caveats, potentially increases the probability of life emerging on other worlds.

$\require{mediawiki-texvc}$Galaxy assembly was already well underway in the first 400 Myr of cosmic time, as recently revealed by JWST. However, the contribution of these early galaxies to cosmic reionisation remains uncertain. Here we present new JWST/NIRSpec observations of GS-z13-1-LA obtained as part of the OASIS and JADES programmes, whose combined deep (56 h) NIRSpec/PRISM spectrum confirms the Lyman-$\mathrm{\alpha}$ line detection and blue UV continuum at redshift $z = 13.1$ presented in a previous work. The measured Lyman-$\mathrm{\alpha}$ emission (rest-frame equivalent width of $66_{-9}^{+10}\,Å$) and steep continuum slope ($\beta_\text{UV} \approx -3$) point towards GS-z13-1-LA hosting a remarkably hot and powerful ionising source, and allow at most a modest contribution from the nebular continuum. The steep turnover of the continuum is still present, but less pronounced in the new OASIS spectrum. Combined, this implies that ionising photons may escape GS-z13-1-LA at a sufficient rate to weaken the other, still undetected UV lines, and to lead the formation of a small ionised bubble ($R_\text{ion} \approx 0.2\,\mathrm{pMpc}$). A yet larger bubble could alleviate the required ionising production efficiency of GS-z13-1-LA from $\xi_\mathrm{ion} \approx 10^{26.4}\,\mathrm{Hz\,erg^{-1}}$ down to $\approx 10^{25.9}\,\mathrm{Hz\,erg^{-1}}$, still extremely high but more readily reconcilable with stellar models. In turn, this would require a notable overdensity of galaxies with highly efficient ionising capabilities, a scenario for which tentative evidence is found in the form of 16 nearby photometric candidates and one spectroscopically confirmed source, JADES-GS-z13-0. The new OASIS observations therefore confirm the overall picture of GS-z13-1-LA as an early beacon of reionisation, providing compelling evidence for its start only 330 Myr after the Big Bang.

Using $JWST$/NIRCam data over a 0.28 deg$^{2}$ area from COSMOS-Web survey, together with $HST$/ACS data, we investigate early-type fraction of massive galaxies with $M_{star}>10^{10.5}M_{\odot}$ at $0.2<z<2.0$, and explore the formation of their early-type morphology. We measure concentration index $C$ ($=R_{80}/R_{20}$) and asymmetry index $A$, and select early-type galaxies with $C>C_{n=2.5}$ and $A_{cor}<0.2$. Here $C_{n=2.5}$ is the concentration expected for a Sersic profile with $n=2.5$ under the spatial resolution and depth of the data, and $A_{cor}$ is the asymmetry corrected for resolution effects. The fraction of early-type galaxies with $M_{star}>10^{11}M_{\odot}$ ($=10^{10.5}$-$10^{11}M_{\odot}$) decreases with increasing redshift from ~70% (~40-60%) at $z$ ~ 0.3 to ~20-25% (~15-25%) at $z$ ~ 1.8. We also examine the evolution of their $R_{20}$ and $R_{80}$, which enclose 20% and 80% of the total flux of the galaxy, respectively. The median $R_{80}$ shows strong mass dependence and significant redshift evolution, whereas the median $R_{20}$ shows little dependence on either stellar mass or redshift. In contrast, morphological differences are more pronounced in $R_{20}$ than in $R_{80}$: the median $R_{20}$ of early-type galaxies is smaller than that of late-type and irregular galaxies by 0.25-0.45 and 0.3-0.6 dex, respectively. The median SSFR of sample galaxies strongly correlates with $R_{20}$, and early-type galaxies have lower SSFRs by ~1 dex. We further find that early-type galaxies at $z>1.3$ have younger mass-weighted stellar ages of $t_{mw}<2$ Gyr than late-type and irregular ones. Their SSFRs, $t_{mw}$, and morphological properties suggest that these high-$z$ early-type galaxies experienced rapid formation of a dense stellar core through starburst, followed by quenching of star formation, and subsequently resumed star formation ~1-2 Gyr later.

Proto-stellar clusters, likely progenitors of globular clusters, are compact with typical mass $\sim 10^6\,{\rm M}_\odot$ and size $\sim 1\,{\rm pc}$, as revealed recently by JWST observations at $z\sim 10$. Sufficiently high compactness can provide a time window for early-formed stellar black holes (BHs) to accrete primordial gas. We develop a model to determine the final spin distribution of stellar BHs which grow in mass via gas accretion within compact gaseous proto-stellar clusters. The velocity shear within a BH's sphere of influence induces the formation of an accretion disk which is repeatedly disrupted by stochastic perturbations to the BH motion. We assume low initial BH spins $a_{*,{\rm ini}} = 0.01$, and restrict initial BH masses below the upper BH mass gap, $m_{\rm BH,ini} < 55\,{\rm M}_\odot$. Our analysis shows a strong BH spin-mass correlation, obtained within $\sim 10 \,{\rm Myr}$ when gas is depleted. Low-spin BHs, $a_{*} \leq 0.3$, are predominantly low-mass, $m_{\rm BH} \lesssim 25\,{\rm M}_\odot$, in contrast to high-spin black holes, $a_{*} \geq 0.7$, which are predominantly high-mass, $m_{\rm BH} \gtrsim 65\,{\rm M}_\odot$. Notably, there exist also low-spin, high-mass outliers with $\sim 1$ mass-gap BH per cluster expected to have $a_{*} \sim 0.1$. The general trend, however, expressed by the median spin as a function of final BH mass is well fit by a high-spin saturating exponential with transition mass $\sim 50\,{\rm M}_{\odot}$. For $m_{\rm BH} \geq 100\,{\rm M}_\odot$ the median spin is $\bar{a}_{*} \sim 0.90$ with the central $68\%$ of the distribution spanning $a_{*} \sim 0.70 - 0.96$, in striking agreement with the estimated spins of the gravitational-wave signal GW231123. These spin values persist up to the highest masses generated by our mechanism, $m_{\rm BH} \sim 10^3\,{\rm M}_\odot$.

Large high-quality samples of HII regions and their parent Giant Molecular Clouds (GMC) are now available for local galaxies. It is therefore possible to investigate links between the CO and H$\alpha$ luminosity functions and whether massive stars form in GMCs of all masses. The CO luminosity functions (LF), representing the distribution of GMC masses, are consistently steeper than the H$\alpha$ luminosity functions. The CO LF invariably steepens in the outer disk where fewer massive GMCs are present beyond the median cloud galactocentric distance. The H$\alpha$ LF also steepens in the outer disk for most of the galaxies examined. Using Salpeter, Kroupa, and Chabrier Initial Mass Functions (IMF) along with stellar mass-luminosity-radius relations, we compute numerically the bolometric luminosity and H$\alpha$ emission from young star clusters. The cluster masses are linked to the GMC mass by assuming that the cluster mass is a constant fraction (3\%) of the parent cloud mass. In particular, results for a fully stochastic IMF are compared to suggestions that very massive stars only form in massive clusters or clouds. Within the limits of the observations -- no small molecular clouds or low-luminosity HII regions can be detected at the typical $\sim 10$~Mpc distance of the sample galaxies -- we find no evidence for a maximum stellar mass which varies with cloud or cluster mass.

Astrophysical luminous objects such as the first stars have not yet formed in the Dark Ages. However, primordial black holes (PBHs) always exist throughout cosmic history since the inflation epoch. During the Dark Ages, PBHs may accrete the ambient gas and release radiation like astrophysical luminous objects, change the cosmic radiation field, the thermal status of the intergalactic medium (IGM), and the hydrogen spin temperature. The accretion rate is modulated by the relic supersonic relative streaming velocities between dark matter (DM) and baryons, imprinting Velocity Acoustic Oscillations (VAOs) features in the 21 cm power spectrum. Such VAOs features could be a promising probe for detecting the PBHs in Dark Ages. We find that even if PBHs comprise only a small fraction of DM, they can generate VAOs wiggles with a relative amplitude up to about 30% in Dark Ages. For example, for PBHs with a mass at recombination of 200 solar masses and mass fraction in the total DM f_PBH,rec around 1e-13 at the recombination era, VAOs features appear at redshift around 20; if f_PBH,rec is around 3e-10, then VAOs features could appear as early as redshift around 40. Moreover, the redshift evolution of the VAOs features exhibits clearly separated stages dominated by inhomogeneous Ly-alpha scattering, and inhomogeneous X-ray heating, respectively. It reflects the characteristics of PBHs (mass and fraction in total DM) and their interactions with the IGM. We also estimate that, the VAOs wiggles at redshift around 20 are detectable for the upcoming SKA-low AA*, while wiggles at redshift around 40 are detectable for an hypothetic lunar surface-based interferometer array in the future.

Polycyclic aromatic hydrocarbons (PAHs) are fundamental to understanding the interstellar medium (ISM) of several astrophysical objects. Normally present in Starburst (SB) galaxies, they have also been more frequently detected in active galaxy nuclei (AGNs), suggesting an inner dusty torus that can shield the radiation from the central black role. In this work, we analyze the 6.2 $\mu$m PAH band of SB-, AGN- and mixed-dominated spectra from 175 IDEOS database galaxies. After fitting of the band, the sources were distributed into the Peeters' A, B and C classes according to their profile peak positions. Class A objects are predominant in 80% of the entire sample, which could indicate the presence of PAHs with nitrogen incorporation. The water ice absorption at 6.0 $\mu$m was also studied in eleven objects, and it affected the PAH band poorly. A prominent second spectral feature after 6.3 $\mu$m is present in ten galaxies. Fitting both PAH profiles at 6.2 $\mu$m changes all the fit results: the first profile is consistently blue-shifted and classified as class A due to the presence of the second component. Further studies are needed to better comprehend these PAH trends in galactic environments.

The organosulfur biosignature gases dimethylsulfide (DMS) and dimethlydisulfide (DMDS) have recently been claimed to be present in the atmosphere of sub-Neptune exoplanet K2-18b, leading to the suggestion of possible extraterrestrial life. Abiotic formation pathways for DMS and DMDS in reducing atmospheres have also been proposed, raising concern over the use of DMS and DMDS as biosignature gases more generally. In this paper we independently test and contrast the proposed abiotic formation pathways for DMS and DMDS using K2-18b as a case study, and explore the wider implications for the atmospheric carbon and sulfur chemistry of hydrogen-rich sub-Neptunes. We demonstrate that one proposed formation pathway is capable of producing observable abundances of abiotic DMS and DMDS, however it depends sensitively on the energy barrier of the limiting step, which remains unmeasured experimentally. In contrast, hydrocarbons including C2H6 are formed abundantly in such atmospheres and offer a plausible alternative explanation to the reported suggestions of organosulfur compounds on K2-18b, having previously been shown to be degenerate observationally. Finally, we demonstrate that sulfur hazes form via the photochemistry of H2S and condense in the atmosphere of K2-18b even at trace abundances. We propose that variation in atmospheric sulfur abundance can explain the diversity of haziness observed across the sub-Neptune population so far with JWST.

We present a statistical analysis for the interaction-induced central concentration of HI gas distributions and its connection with interaction-induced central star formation enhancement, using a large sample of $\sim 10^4$ galaxies from the ALFALFA and SDSS surveys. By adopting the HI profile parameter $K$, an indicator of gas concentration inferred from the integrated 21 cm emission line, we find that galaxies with more centrally concentrated HI (higher $K$ values) or enhanced specific star foramtion rate (sSFR) exhibit significantly stronger clustering and higher probability of hosting a nearby neighbor on scales below $100h^{-1}\mathrm{kpc}$, which is more pronounced in low-mass galaxies. Furthermore, by utilizing the enhancement functions for a sample of galaxy pairs, we directly trace the evolution of HI concentration and sSFR enhancement as a function of projected separation. Our findings indicate that tidal interactions drive a statistical synchrony between the central concentration of atomic gas and the enhancement of central star formation. Gas concentration appears to be a necessary condition for central star formation enhancement in interacting systems at all but the smallest separations. Compared to satellite galaxies, central galaxies exhibit stronger enhancement of gas fraction, gas concentration and sSFR, suggesting the role of environmental regulation.

We present HECATEv2, the second release of the Heraklion Extragalactic Catalogue (HECATE), an all-sky, value-added galaxy catalogue comprising 204733 galaxies from the HyperLEDA database with recession velocity <14000 km/s (D~200 Mpc). This release focuses on qualitative upgrades of the provided information while maintaining the same parent galaxy sample as HECATEv1. Improvements include a new cosmology-based distance framework, expanded and homogenised optical and mid-infrared photometry from SDSS-DR17/NSA, PS1-DR2, and AllWISE, and new quality-control flags for stellar contamination, incorrect photometry, and coordinate inconsistencies. We also extend the galaxy-size coverage and derive stellar population parameters for a substantially larger fraction of the sample. Star-formation rates (SFR) and stellar masses (Mstar) are now available for >70% of galaxies using updated mid-IR/optical calibrations that account for stellar population age and dust attenuation, while gas-phase metallicities are derived for ~90%. Activity classifications are provided for >50% of galaxies based on spectroscopic and/or photometric diagnostics, and supermassive black hole masses for ~86%. In terms of L$_{B}$,L$_{Ks}$,SFR, and Mstar, HECATEv2 is among the most complete local-Universe catalogues with spectroscopic redshifts. We also provide spatial completeness maps as a function of distance and luminosity, highlighting variations across the sky. Compared to other catalogues (e.g. GLADE+, NED-LVS), HECATEv2 offers broader (optical, near- and far-IR photometry, metallicity, activity classifications) or comparable (mid-IR photometry, SFR, Mstar) coverage, making it a robust reference for studies of SMBH-host galaxy connections, gravitational-wave and high-energy transient hosts, population analyses, and rare galaxy subpopulations.

We present results from our ongoing campaign to follow up the satellite candidates from the Identifying Dwarfs of MC Analog GalaxiEs (ID-MAGE) survey. Previously, we published a list of 355 unresolved satellite candidates identified around 36~nearby LMC- and SMC-mass hosts (D$=$4$-$10~Mpc). We present the velocities of 83 satellite candidates from new Green Bank Telescope \hi\ observations, optical long-slit spectra, and the Dark Energy Survey Instrument Data Release 1. Based on their velocities, we identify six candidates as probable satellite galaxies ($6.5\times10^5\leq M_\star/M_\odot\leq1.5\times10^7$) and 77 as background galaxies. Our results underscore the ability of spectroscopic follow-up to effectively separate satellites from background galaxies. Using the refined sample, we update our previously derived estimates for the average satellite population per host and find 1.7$\pm$0.7 (1.0$\pm$0.3) satellites per LMC-mass (SMC-mass) host. Our current satellite sample includes 25 galaxies confirmed by distances or velocities. This set includes the complete satellite populations of three hosts (UGC~04422: zero satellites, UGC~08201: zero satellites, NGC~3432: four satellites), which we compare to simulations and known satellite systems from the literature. Our sample is nearly complete for the most massive satellites (M$_\star > 10^7~M_\odot$). We find these massive satellites have a quenched fraction of 10--25\%, placing them between the $<$5\% quenched fraction of isolated galaxies and the 40--70\% quenched fraction of MW-analog satellites with $10^7~M_\odot < $ M$_\star < 10^8~M_\odot$. This demonstrates the impact that low-mass galaxies have on the evolution of their satellites.

Baryon acoustic oscillation (BAO) analyses usually report the anisotropic shift parameters $\alpha_\perp(z)$ and $\alpha_\parallel(z)$ relative to a fiducial cosmology, and these quantities are primarily used for cosmological parameter inference. Here we show that they can also be used to construct a direct internal consistency test of the background geometry. In particular, we derive a new null test of flat Friedmann-Lemaître-Robertson-Walker (FLRW) geometry written entirely in terms of the reported BAO shift parameters. The test is calibration free: the sound-horizon ratio $r_{\rm d}/r^{\rm fid}_{\rm d}$ cancels identically, so the relation is independent of the absolute BAO scale. We also derive a calibration-free reconstruction of the deceleration parameter $q(z)$ from the radial BAO sector. Applying these results to anisotropic DESI DR2 BAO measurements, we find no evidence for a breakdown of the flat-FLRW distance relation within current uncertainties. Our results show that anisotropic BAO measurements already provide a nontrivial internal geometric consistency test before performing any model fit.

We provide the community with a homogeneous catalogue of small, significant substructures (henceforth NESTs) extracted from the spatial distribution of Young Stellar Objects (YSOs) in a large, consistent sample of star-forming regions. The catalog allows us to explore the relevance of small scale spatial substructure and discuss the interpretation of NESTs as tracers of star formation activity and remnants of the star formation process. We apply our procedure to consistent catalogues of YSOs to obtain NESTs in a sample of star-forming regions. We apply a photometric classification scheme to obtain the evolutionary stage of YSOs and statistically explore the distribution of class 0/I objects as a proxy of recent star formation activity. The region sample is diverse (in distance, size, structure, and global evolutionary stage), and we consequently find different structural properties and star formation histories. Most NESTs in regions with high recent star formation activity show even higher levels of activity. Moreover, the proportion of NESTs with higher activity than the region average increases with the global level of activity of the region. In approximately half of the regions we also find significant spans in the evolutionary stages of the NESTs, consistent with gradients and episodes of star formation. The combination of NESTs with a statistical exploration of the star formation history within each region provides robust and powerful insights into the star formation process. Our results support the role of NESTs as pristine remnants of star formation in highly active regions,stressing the role of fragmentation. The combination of small structures with large scale spatio-evolutionary patterns suggests hyerarchical, prolonged, dynamic, and complex star formation scenarios.

3I/ATLAS is the third confirmed interstellar object passing through the Solar System. In this work, we conduct narrowband radio technosignature search toward 3I/ATLAS using the Five-hundred-meter Aperture Spherical Telescope (FAST) L-band multibeam receiver from October 2025 to January 2026 on 4 separate dates (i.e. Mars closest, perihelion, Earth closest and flew away from Earth, respectively). We carry out frequency-drifting signal searching with signal-to-noise ratio (SNR) over 10 within 1.05-1.45 GHz via \texttt{bliss} pipeline. These signal hits are grouping into event by beam, frequency and drift rate matching, the events are then filtered by cluster analysis and drift rate cut-off. We also characterized the events by their significance in SNR, structure tensor as well as principal component analysis (PCA). No credible narrowband radio technosignature are detected from 3I/ATLAS after visual inspections. The null results place constraints on the presence of transmitters above $2.862\times 10^{-3}$ W. We further introduce a Bayesian inference framework to constraints on the existence probability and characteristic power of hypothetical transmitters using physically motivated priors to bracket plausible transmitter scenarios.

Neptune's present axial tilt of approximately 28 deg. with respect to its orbital plane can be explained by collisions that its primordial core may have experienced with surrounding planetary embryos during the final stages of its formation. Alternatively, Neptune could have attained its present mass solely through pebble accretion, without the formation of nearby planetary embryos. The embryo-collision scenario has the advantage of naturally explaining the large axial tilts observed in the ice giants. To account for these tilts without invoking late-stage catastrophic collisions, an alternative mechanism must be considered. In this work, I propose that Neptune's current axial tilt could result from the interaction between Triton's tidally evolving orbit and Neptune's spin axis, causing it to resonate with the solar system eigenfrequency s_8. Starting from a Triton-mass satellite captured via the binary planetesimal disruption mechanism, I show that orbital evolutions bringing the satellite near Triton's present orbit can induce a spin-s_8 resonance capable of producing a significant axial tilt of Neptune's spin axis. I develop a model for planetary spin-axis evolution using Euler's equations for a rigid body, which is incorporated into classical numerical integrations of the Newtonian equations of motion. I also include a tidal model to account for the satellite's semimajor-axis decay and orbital circularization. Several numerical simulations are performed with this model, including Neptune as the central body, the newly captured satellite, the Sun and the three other giant planets. Increases in Neptune's obliquity are observed, exceeding 50 deg. in some cases. An obliquity above 20 deg. is obtained for roughly 1/3 of the cases. If Neptune initially had a near-zero obliquity, its current value could therefore have been naturally engendered by the tidal evolution of Triton.

Cosmic-ray pitch-angle transport in magnetohydrodynamic (MHD) turbulence is governed by the interplay between magnetic mirroring and gyroresonant scattering. We develop a guiding-center (GC) Langevin model with explicit mirror drift and gyroresonant diffusion to describe the pitch angle evolution. This model accurately captures our test-particle simulation results in three-dimensional MHD turbulence, driven both solenoidally and compressively. We find that magnetic mirroring can drive anomalous pitch-angle diffusion at large pitch angles (including $90^\circ$) with non-Markovian memory effects, which arises from trapping of particles in magnetic wells. Gyroresonant scattering controls the escape rate from these wells. Across $M_{\rm A}$, large-pitch-angle particles are jointly regulated by mirror trapping and gyroresonant escape, exhibiting a transition from anomalous to normal diffusive pitch-angle transport as scattering strengthens, whereas small-pitch-angle particles remain gyroresonance-dominated and diffusive throughout. The pitch angle transport is found to be dominated by the compressible perturbations with marginal influence from Alfvén modes. In compressible turbulence with realistic damping accounted for, transit time damping (TTD) treatment fully recovers mirror interactions.

Traditional weak gravitational lensing shear estimators are carefully calibrated but struggle to fully capture realistic galaxy morphologies, point-spread-function (PSF) effects, blending, and noise in deep surveys, while blindly trained machine learning (ML) models can introduce significant calibration biases. Here we construct a fully D$_4$-equivariant deep neural network for galaxy shape measurement whose architecture enforces symmetry under 90$^{\circ}$ rotations and mirror transformations, and adopt the Analytical Calibration framework (AnaCal) to calibrate the model using its backpropagated gradients. For isolated galaxies in LSST-like single-band simulations, we demonstrate that our approach achieves $\sim$10% lower shape noise than the traditional moment-based Fourier Power Function Shapelets estimator in the high-noise regime, equivalent to a $\sim$20% gain in effective galaxy number density, while simultaneously achieving multiplicative biases consistent with zero across a wide range of noise levels, PSF sizes and ellipticities, and magnitude selection cuts, with all measurements satisfying $|m| {<} 10^{-3}$ (i.e., within the 0.2% LSST requirement) and most at the ${\sim}10^{-4}$ level. We demonstrate this framework on isolated single-band galaxy images with Gaussian noise and known PSF, establishing a rigorous, physics-informed foundation for future extensions of ML-based shear estimation to blended sources and multi-band observations in Stage-IV surveys. All codes and data products will be made publicly available upon acceptance.

Here we present broadband radio modeling of GRB 240205B, using observations with the Australia Telescope Compact Array (ATCA) and the South African MeerKAT radio telescope. Our observations include an automatically triggered early-time ATCA observation that began approximately 13 minutes after the gamma-ray signal and continued for 12 hours, resulting in the earliest detected GRB radio afterglow to date at about 35 minutes post-burst. Following this initial detection, we conducted an extensive radio follow-up campaign for more than 5 months. Although the observations beyond one day post-burst are well described by a standard forward shock model, the observation before one day post-bust reveals an additional synchrotron component, which can be explained as the reverse shock. This component would have been missed without the automated ATCA rapid-response trigger. We find that a combined reverse and forward shock model in a stellar wind medium best describes the radio afterglow. We constrain the spectral breaks due to synchrotron self-absorption and the minimum electron energy, and we use the light-curve peaks to constrain the microphysical parameters. We put GRB 240205B in the context of the growing sample of GRBs with radio detections in the first hours after the gamma-ray trigger. Using our rapid response observation, we estimate the highest model independent constraint on a GRB minimum bulk Lorentz factor of around 100 at about 35 minutes post burst. We also discuss future prospects of detecting similar long GRBs at centimeter wavelengths, as well as potential improvements to future strategies for targeting their radio afterglows.

The nature of the so-called G objects orbiting the Galactic Center remains unresolved. These sources exhibit compact Br$\gamma$ emission, extreme infrared colors, and remarkable dynamical stability through close passages to the central supermassive black hole, challenging conventional interpretations as stars or unbound gas clouds. We investigate the hypothesis that G objects are the remnants of neutron stars that have been converted into low-mass black holes through the capture of primordial black holes, a viable dark-matter candidate. We construct a population-level framework linking the abundance and spatial distribution of these remnants to the neutron-star population, the inner dark-matter density profile, and the primordial black-hole mass and abundance. Within this framework, the observed G-object population and the long-standing deficit of ordinary radio pulsars in the Galactic Center emerge as complementary consequences of the same conversion process. We further identify a suite of observational signatures-across infrared, radio, X-ray, and microlensing channels-that render this scenario empirically testable and distinguishable from stellar-envelope models. Our results show that G objects can act as sensitive probes of compact-object capture physics and of dark matter on sub-galactic scales.

Io's internal heat is generated by Jupiter-driven tidal dissipation and Laplace resonance. This energy partially melts the mantle, but the melt fraction, depth, and spatial distribution of dissipation remain poorly constrained. Tidal deformation is linked to the mantle's physical state via a parametric approach accounting for melting onset depth and latent heat of fusion. Io is modeled as a three-layer body comprising a fluid core, a viscoelastic mantle, and an elastic lithosphere. The degree-2 potential Love number k2 is computed by solving spheroidal oscillation equations with an adapted California Planetary Geophysics Code (CPGC). Mantle properties (viscosity, shear modulus, Andrade parameter beta) are iteratively updated based on the local melt fraction, and mantle compressibility is explicitly incorporated. Results show a self-consistent shallow-mantle enhancement of tidal heating. Melt presence decreases effective viscosity and increases anelasticity, amplifying upper-mantle dissipation. Reproducing the observed real part of k2 in 1D models requires melt fractions below the rheologically critical threshold. Compressibility yields higher values for the real part of k2, reinforcing this upper bound. Mass flux analysis confirms melt percolation capacity exceeds thermodynamic production, implying efficient drainage. The reference Andrade parameter beta strongly influences the imaginary Love numbers (k2, h2, l2) and the predicted libration amplitude. These constraints support a heterogeneous mantle characterized by a magmatic sponge structure rather than a global magma ocean, establishing a quantitative framework linking Io's interior, partial melting, and tidal dissipation to Juno observations.

Spin-orbit obliquity measurements of hot-Jupiter systems constrain giant planet migration and tidal evolution. In binary systems, combining stellar obliquities with the orbit-orbit angle ($\gamma$) between the planetary and stellar companion orbits provides further insight into the dynamical influence of stellar companions. Here we aim to determine the projected obliquity ($\lambda$) of the hot Jupiter TOI-1333Ab ($P\approx4.72$ d, $M_{\rm p}\approx2.4$ M$_{\rm J}$) and place the system in the context of hot-Jupiter migration and tidal realignment in binary systems. We analysed spectroscopic observations obtained during planetary transit to model the Rossiter-McLaughlin effect and derive the projected obliquity. We combined this measurement with published system parameters and constraints on the wide stellar companion orbit to assess plausible migration scenarios. We measure a projected obliquity of $\lambda=-5 \pm 10^\circ$, showing that TOI-1333Ab is well aligned with the stellar spin axis of its F-type host star. The low obliquity and its modest eccentricity ($e=0.073^{+0.092}_{-0.052}$) are consistent with either disc-driven migration or high-eccentricity migration followed by efficient tidal circularisation and realignment. With an effective temperature of $6274\pm94$ K, the host star lies above the canonical Kraft break where the systems are frequently misaligned. Despite this, we find the system to be well aligned. In comparison with other planetary systems in binaries, TOI-1333 occupies a relatively isolated region in projected obliquity-orbit-orbit angle ($\gamma=81.5\pm1.1^\circ$) space, making it a valuable system for studying the interplay between migration, tides, and stellar companions.

We present a high-efficiency, high resolution on-chip filterbank spectrometer designed for line intensity mapping and broadband wave-like dark matter searches. This spectrometer maximizes sensitivity to the faint, aggregate cosmic signal for line intensity mapping while providing redshift precision necessary to resolve the clustering of large-scale structure. For broadband dark matter searches, it partitions broadband signals into narrow frequency bins that directly correspond to dark matter masses, allowing signal frequency to be mapped with high precision to the particle's mass. Existing superconducting filterbank architectures used by the mm-wave community are limited by a 50\% inherent efficiency limit and are highly sensitive to thin-film dielectric loss. The design presented in this paper addresses these bottlenecks by eliminating the termination resistor and employing a niobium-on-silicon coplanar waveguide resonant structures for the filterbanks. Sonnet electromagnetic simulations of a 10-channel device around 90 GHz demonstrates >74\% per channel efficiency and a resolving power of $R=1211\pm105$. Sensitivity analyses confirm that the design is robust against typical fabrication uncertainties with the exception of dielectric thickness, providing a scalable technology solution for the next generation of millimeter-wave spectroscopic experiments.

We assess the impact of non-ideal, continuously rotating half-wave plates (HWPs) on cosmic microwave background (CMB) polarization measurements targeting large angular scale signal. Such hardware solutions are used in or planned for multiple modern CMB efforts, both ground-based, for instance, small aperture telescopes of Simons Observatory or satellite borne, such as LiteBIRD. Using a frequency-dependent parametric model based on the Mueller matrix formalism, we characterize the induced mixing of Stokes parameters. Through end-to-end simulations, we propagate these effects from time-ordered data to cosmology via map-making and component-separation stages, quantifying their impact on the $B$-modes power spectrum and the tensor-to-scalar ratio, $r$. Our analysis shows that neglecting the frequency dependence of a three-layer HWP gives rise to significant polarization leakage, biases foreground spectral parameters, and leads to residual contamination in the recovered CMB maps. To mitigate these effects, we investigate multiple analysis strategies progressively incorporating a more complete description of the instrumental response. At the map-making level, this requires generalizing the standard pointing matrix to account for the full time- and frequency-dependent instrumental response. We find that standard HWP models, reduce the biases only down to $r \sim 10^{-2}$, while a more advanced approach based on a generalization of both map-making and component separation, implemented using JAX, can suppress it down to $r \sim 7 \times 10^{-4}$. Finally, we extend this approach to a time-domain component-separation, enabling a statistically consistent treatment of instrumental response in the presence of time-domain features. We demonstrate its feasibility and validate it by performing a full end-to-end analysis, recovering results in good agreement with the map-based ones.

We present the first semi-empirical strong-line calibrations to determine metallicity in Active Galactic Nuclei (AGNs) that use the directly observable X-ray luminosity ($Ł_{\rm X}$) instead of the dimensionless ionization parameter ($U$). The calibrations are derived from an extensive grid of photoionization models computed with the {\sc Cloudy} code, which are compared with observational data of Seyfert nuclei from the Burst Alert Telescope (BAT) AGN Spectroscopic Survey (BASS). In this first paper, we develop new calibrations for two key optical metallicity diagnostics based on the $N2$ and $O3N2$ indices, which are valid in a metallicity range of $8.0 \lesssim \logOH \lesssim 9.1\, {\rm or}\, 0.2 \lesssim (Z/Z_{\odot}) \lesssim 2.6$, with precision of $1\sigma \approx 0.22$ dex ($N2$) and $\approx 0.20$ dex ($O3N2$). We systematically investigate the influence of the AGN spectral index $(\aox)$, narrow-line region (NLR) gas density (\Ne), the characteristic peak temperature of the Big Blue Bump $(T_{\rm BB})$, and $Ł_{\rm X}$. We find a strong, opposing secondary dependence on $Ł_{\rm X}$ for both indices. We demonstrate that neglecting this parameter overlooks systematic offsets intrinsic to the diagnostics, leading to metallicity errors of up to $\sim 1.0$ dex, particularly for the least and most luminous sources. This framework offers a more precise characterization of chemical enrichment in the NLRs of AGNs by leveraging their intrinsic X-ray emission to mitigate these systematic biases.

We present a comprehensive chemical and age orbital cartography of the Galactic high-$\alpha$ disc using subgiant stars with precise ages, element abundances, and full phase-space information from the \textsl{LAMOST--Gaia} data set. Specifically, we map how average [Fe/H], [$\alpha$/Fe], and age vary across present-day kinematic and orbital coordinates. We analyse the data in full and across mono-abundance populations to measure element abundance-orbital and age-orbital gradients across orbital actions and angular-momenta. Our results show that the high-$\alpha$ disc exhibits clear and coherent gradients in [Fe/H], [$\alpha$/Fe], and age with orbits; these gradients are much stronger and sharper in orbital space than in present-day kinematics, showing that orbital diagnostics recover the intrinsic disc structure of old disc populations more effectively than instantaneous kinematic coordinates. We find that older high-$\alpha$ populations display qualitatively similar element abundance--orbital and age--orbital trends to stars in the low-$\alpha$ disc, although the high-$\alpha$ gradients are generally shallower. The presence of these ordered correlations indicates that the old high-$\alpha$ disc is structured, and preserved a strong fossil record of its early assembly despite the Milky Way's subsequent accretion history. This result implies that later mergers did not fully erase the chemical-orbital and age-orbital structure imprinted during the high-$\alpha$ disc's earliest formation epoch. All together, these findings indicate that the Galactic high-$\alpha$ disc formed mainly through inside-out and upside-down growth.

Thermal misalignment provides an alternative to the standard misalignment mechanism for the cosmological production of scalar dark matter. In this framework, feeble couplings to particles in the thermal bath generate a finite-temperature potential that drives the scalar towards large field values early in the radiation era, dynamically inducing the misalignment before the onset of scalar oscillations. As a result, the relic abundance is controlled primarily by particle masses and couplings rather than the initial field value. As a light spectator field, the scalar acquires inflationary fluctuations that are uncorrelated with the adiabatic curvature mode, generically sourcing isocurvature perturbations. We show that, unlike standard misalignment, where light scalars are strongly constrained by cosmic microwave background bounds on dark matter isocurvature for high-scale inflation, thermal misalignment can naturally suppress the isocurvature signal. This occurs through a novel late-time phase offset between the background zero mode and the superhorizon perturbations, which reduces the final dark matter density contrast. Thermal misalignment therefore provides a new and generic route to isocurvature-safe scalar dark matter.

We propose the idea that our Universe is a realization among different possible branches, which can be observationally tested through the modified dispersion relation of the gravitational waves. We achieve this through a framework of spatially constrained vector fields. We show that the simplest realizations of such theories in flat and cosmological spacetimes do not introduce new propagating modes, but they give rise to tensor perturbations that differ from those of standard general relativity. We further show that such theories admit stealth black hole solutions, and we recover weak gravitational potentials, thus passing the solar system experiments. Finally, we discuss the implications of such theories and propose further generalizations.

Tentative evidence suggests that the cores of massive neutron stars consist of deconfined quark matter. We argue that the formation of such a quark matter core during a galactic supernova could be accompanied by the emission of gravitational waves in the MHz band. These signals constitute a new target for high-frequency gravitational wave detectors, demonstrating that such detectors may offer unique opportunities for testing quantum chromodynamics in an otherwise inaccessible regime.

The meV mass range has emerged as a focal point in axion physics, where advances in theory, cosmology, astrophysics, and experimental techniques converge. Axions in this mass range are theoretically well motivated, can arise in ultraviolet-complete models, and can have significant cosmological impacts as dark matter or dark radiation. In parallel, their efficient production in stellar and supernova environments provides powerful astrophysical probes. Here, we provide a comprehensive overview of meV axions across these domains, highlighting both established results and open questions. We discuss the theoretical underpinnings of meV axions, their cosmological and astrophysical signatures, and the diverse experimental strategies -- ranging from helioscopes and haloscopes to quasiparticle systems and large-volume Cherenkov detectors -- that aim to explore this regime. The convergence of these approaches emphasizes the pivotal role of the meV mass range for axion discovery in the coming years, identifying meV axions as a key probe for testing beyond-Standard-Model physics. This review document is the direct outcome of the discussions at the dedicated workshop "The meV Mass Axion Frontier: Challenges and Opportunities", held at Laboratori Nazionali di Frascati (IT) on 27--28 October 2025, and organized by the EU funded COST Action "Cosmic WISPers in the Dark Universe: Theory, astrophysics, and experiments" (CA21106, this https URL). Its aim is to provide an overview of current efforts in meV axion research, their motivations, and the research goals that animate the community involved in this search.

In this paper, we pursue an investigation of the consequences of a mixing between supersymmetric partners - the photon and photino - analogous to the so-called Primakoff effect, but induced by a Lorentz-symmetry violating (LSV) fermionic-condensate background. In our framework, the LSV parameters are introduced as members of a non-dynamical superfield. As a consequence, we show that naturally there appears a mixing term between the gauge boson and the gaugino, which can be readily seen in the superspace/superfield approach. We inspect the kinetic photon-photino mixing matrix in the scenario of stellar physics which we apply our results to. Bounds on the strength of the fermionic LSV background are can be set by invoking the energy loss argument and the solar data we adopt.

Spacetime is foamy due to quantum fluctuations. Various gedanken experiments show that distances fluctuate by amounts consistent with the holographic principle, hence the name "holographic quantum foam" (HQF). One important prediction of HQF is that necessarily there exists a dark sector in the universe. The resulting cosmology is found (at least qualitatively) to be consistent with observations. Interestingly the quanta of the dark sector are found not to obey the familiar (fermionic or bosonic) statistics, but the exotic statistics known as infinite statistics (or quantum Boltzmann statistics). The most important challenge now is to check if HQF is consistent with experiments/observations. One way is to look for observational evidence of blurred distant point-sources due to physics at the Planck scale. For over two decades it has been debated whether those tiny inherent uncertainties in time and path-length can accumulate in transiting electromagnetic wavefronts from quasars and Gamma-Ray Bursts (GRBs). But a recent event is special: GRB221009A was extremely bright and energetic. That allowed follow-up across the whole spectrum from the optical/near-infrared through to X-rays, and including the highest-ever-recorded energy gamma-rays; all consistent with blurring by HQF. Those data, and a calculation of the HQF-widened point-spread function (PSF) for real telescopes viewing a GRB are presented.

We study the cosmological implications of the minimal non-linear realisation of scale invariance within the Standard Model (SM). This framework provides a technically natural explanation for the hierarchy between the Planck scale and the electroweak scale and introduces only a light, feebly coupled dilaton field beyond the SM particles. Although the model is almost indistinguishable from the minimal SM at low energies, its cosmological consequences differ dramatically. In particular, the electroweak Higgs field remains trapped in the symmetric phase until the Universe cools to very low temperatures, $T_c^{(\chi)}\sim 28$ MeV, where the first-order QCD chiral symmetry-breaking phase transition triggers the electroweak phase transition. This scenario offers intriguing possibilities for the production of primordial black holes, low-frequency gravitational waves, and multi-quark and lepton nuggets, which we explore in some detail using simplified approximations.

We review two types of fractional Kaehler potentials $K$ which reduce, along the inflationary path, to the form $N/(1-\phi^{q_{\rm M}})^{p}$ with $q_{\rm M}=1$ or $2$ and $0.1\leq p\leq10$. Their coexistence, within a non-linear sigma model, with chaotic inflationary potentials of the form $\phi^n$ (where $n=2$ or $4$) determines, independently from $q_{\rm M}$ and $n$, a class of $(p,N)$-inflationary attractors which leads to observables compatible with the ACT DR6. An implementation of these models in the context of supergravity can be also achieved by introducing two chiral superfields and a monomial superpotential, linear with respect to the inflaton-accompanying field, and supplementing the $K$'s above with a shift symmetry. Although inflation is attained for subplanckian inflaton values, the tensor-to-scalar ratio obtained for certain $N$ values can be possibly observable in the near future.

We investigate turbulent spin-down dynamics in self-gravitating Bose-Einstein condensates, comparing purely bosonic and axionic (higher-order interacting) systems. Through simulations of the Gross-Pitaevskii-Poisson system, we study condensates pinned to a crust potential undergoing rapid rotation slowdown. We find that axionic condensates exhibit more uniform density profiles and smaller sizes compared to their bosonic counterparts for similar interaction strengths, which facilitates earlier vortex entry. The sudden spin-down triggers vortex depinning and a turbulent cascade. For comparable sizes, both systems exhibit a short-lived Kolmogorov energy cascade ($k^{-5/3}$ scaling) followed by a transition to Vinen turbulence ($k^{-1}$ scaling). Crucially, their responses diverge with increasing interaction strength (and thus condensate size): the axionic system increasingly deviates from Kolmogorov scaling because of enhanced vortex retention, a trend quantitatively confirmed by analyzing the vortex fraction and its dependence on the final rotation frequency. Spectral analysis reveals that the growth of incompressible energy is primarily driven by quantum pressure during vortex detachment, rather than by compressible flows. The compressible spectrum shows thermalization ($k$ scaling). Our results demonstrate how distinct nonlinearities govern vortex dynamics and turbulent dissipation in self-gravitating quantum fluids.

The worldwide LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors continues to increase in sensitivity, thus increasing the quantity and quality of the detected GW signals from compact binary coalescences. These signals allow us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. This paper is the first of three, where we present the results of a suite of tests of GR using the binary signals included in the fourth GW Transient Catalog (GWTC-4.0), i.e., up to and including the first part of the fourth observing run of the detectors (O4a). We restrict our analysis to the 91 confident signals, henceforth called events, that were measured by at least two detectors, and have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. These include 42 events from O4a. This first paper presents an overview of the methods, selection of events and GR tests, and serves as a guidemap for all three papers. Here we focus on the four general tests of consistency, where we find no evidence for deviations from our models. Specifically, for all the events considered, we find consistency of the residuals with noise. The final mass and final spin as inferred from the low- and high-frequency parts of the waveform are consistent with each other. We also find no evidence for deviations from the GR predictions for the amplitudes of subdominant GW multipole moments, or for non-GR modes of polarization. We thus find that GR, without new physics beyond it, is still consistent with these GW events. The results of the two additional papers in this trio also find overall consistency with vacuum GR, with more than 90% of the events being consistent with GR at the 90% credible level. While one of the ringdown analyses finds the GR value in the tails for its combined results, this may be due in part to catalog variance.

In this second of three papers on tests of general relativity (GR) applied to the compact binary coalescence signals in the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), we present the results of the parameterized tests of GR and constraints on line-of-sight acceleration. We include events up to and including the first part of the fourth observing run (O4a) of the LIGO Virgo KAGRA detectors. As in the other two papers in this series, we restrict our analysis to the 42 confident signals, measured by at least two detectors, that have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$ from O4a, in addition to the 49 such events from previous observing runs. This paper focuses on the eight tests that constrain parameterized deviations from the expected GR (or unaccelerated) values. These include modifications of post-Newtonian (PN) parameters, spin-induced quadrupole moments different from those of a binary black hole, and possible dispersive or birefringent propagation effects. Overall, we find no evidence for physics beyond GR, for spin-induced quadrupole moments different from those of a Kerr black hole in GR, or for line of sight acceleration, with more than 90% of the events including the null result (no deviation) within their 90% credible intervals. We discuss possible systematics affecting the other events and tests, even though they are statistically not surprising, given noise. We improve the bounds on deviations from the GR PN coefficients by factors of 1.2-5.5 and provide illustrative translations to constraints on some modified theories. Also, we update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.92\times 10^{-23} \mathrm{eV}/c^2$. Thus, we see that GR holds, and many of the bounds on possible deviations derived from our events are the best to date.

This is the third paper of the set recording the results of the suite of tests of general relativity (GR) performed on the signals from the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), where we focus on the remnants of the binary mergers. We examine for the first time 42 events from the first part of the fourth observing run of the LIGO-Virgo-KAGRA detectors, alongside events from the previous observation runs, restricting our analysis to the confident signals, which were measured in at least two detectors and that have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. This paper focuses on seven tests of the coalescence remnants. Three of these are tests of the ringdown and its consistency with the expected quasinormal mode spectrum of a Kerr black hole. Specifically, two tests analyze just the ringdown in the time domain, and the third test analyzes the entire signal in the frequency domain. Four tests allow for the existence of possible echoes arriving after the end of the ringdown, which are not expected in GR. We find overall consistency of the remnants with GR. When combining events by multiplying likelihoods (hierarchically), one analysis finds that the GR prediction lies at the boundary of the $98.6^{+1.4}_{-9.4}\%$ ($99.3^{+0.7}_{-4.5}\%$) credible region, an increase from $93.8^{+6.1}_{-20.0}\%$ ($94.9^{+4.4}_{-18.2}\%$) for GWTC-3.0. Here the ranges of values comes from bootstrapping to account for the finite number of events analyzed and suggest that some of the apparently significant deviation could be attributed to variance due to the finite catalog. Since the significance also decreases to 92.2% (96.2%) when including the more recent very loud event GW250114, there is no strong evidence for a GR deviation. We find no evidence for post-merger echoes in the events that were analyzed. (Abridged)

Binary neutron star mergers and proto-neutron stars provide unique environments where dense matter is hot, lepton rich, and potentially undergoes a transition from hadronic to deconfined quark matter. We investigate the thermodynamics and stellar properties of hybrid matter under such conditions. The hadronic phase is described within a covariant density functional framework, while the quark phase is modeled using a Nambu-Jona-Lasinio (NJL) model that includes repulsive vector interactions, the axial $U_A(1)$-breaking 't Hooft determinant interaction, and two-flavor color-superconducting (2SC) pairing. The phase transition between hadronic and quark matter is constructed using a mixed-phase prescription that enforces baryon and lepton number conservation, allowing us to follow thermodynamic trajectories at fixed entropy per baryon and fixed lepton fraction. We analyze the phase structure of dense matter at finite temperature and study the composition of the hadronic, mixed, and quark phases in both neutrino-trapped and neutrino-free regimes. Our results show that neutrino trapping significantly modifies the particle composition and shifts the onset of deconfinement to higher densities. Using the resulting equations of state, we compute static stellar configurations and examine the influence of temperature and lepton content on the mass-radius relation of hybrid stars. Hot, neutrino-rich configurations are found to have larger radii and slightly higher maximum masses than their cold counterparts.

The production of very-high-energy (VHE, $E_{\gamma} \gtrsim 100$ GeV) gamma rays resulting from the scattering of high-energy cosmic-ray protons off axion-like particles (ALPs) populating the dark matter halo of the Milky Way is investigated. By employing the latest instrument response functions for current and future facilities, we demonstrate that ground-based VHE gamma-ray observatories, such as H.E.S.S., CTAO, and SWGO, provide a promising and complementary avenue to probe the yet uncharted ALP-proton coupling $g_{ap}$. Our results show that these experiments can reach sensitivity to couplings above $10^{-2}$ in the $1 - 10^{8}$ eV ALP mass range, a region that remains largely unexplored by supernova and neutron star cooling observations. Interestingly, we demonstrate that this search channel is capable of probing QCD axion dark matter models, assuming two benchmark models for it: the Kim-Shifman-Vainshtein-Zakharov (KSVZ) Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) models, specifically within the MeV mass range. These findings highlight the potential of VHE gamma-ray astronomy to provide unique constraints on the interaction between ALPs and the baryonic sector.

We perform a search for an X-ray monochromatic line arising from dark matter (DM) decay in the halo of the Large Magellanic Cloud. An emission line can be expected from two well-motivated DM candidates: sterile neturinos and axion-like particles (ALPs). We analyze the eROSITA-DE DR1 datasets in the energy range between 1 and 9 keV. No evidence for a DM line is found, and we set lower limits on the DM lifetime. We then recast these bounds into upper limits on the active-sterile neutrino mixing angle $\sin^2(2\theta)$ and on the ALP to photon coupling $g_{a\gamma}$, for DM masses between 2 and 18 keV. These results set new strong constraints for masses below 5 keV.

We investigate the phenomenon of black hole superradiance in the presence of torsion within the framework of Poincaré gauge theory. In particular, in contrast to the classical approach of General Relativity, we show that the inclusion of torsion in the space-time geometry enables the energy extraction from rotating black holes by Dirac fermions via chiral asymmetry, while preserving the Pauli exclusion principle.

We review recent progress in the understanding of the preheating stage of Higgs inflation formulated within the Einstein-Cartan framework of gravity. This setup smoothly interpolates between the metric and Palatini formulations of the theory, leading to a distinctive phenomenology in an intermediate regime. Following the end of inflation, the Higgs field undergoes a non-trivial out-of-equilibrium evolution driven by tachyonic instabilities and nonlinear self-interactions, which fragment the inflaton condensate and give rise to well-localized oscillon configurations. While early studies suggested the formation of long-lived oscillons and the possibility of an extended matter-dominated phase, more recent analyses show that self-interactions at small field values render these objects transient, eventually triggering their decay and the onset of radiation domination. We discuss the implications of this dynamics for the thermal history of the Universe, the inflationary observables, and the generation of stochastic gravitational waves.