Papers for Tuesday, Apr 07 2026

Here we report null results in the search for radio emission below $200$~MHz from six known fast radio burst sources. The observations reported here were taken using the Irish LOFAR station's high-band antennas over the course of 2020, 2021 and 2022; the cumulative observing time was $218$~h.

The PLanetary Transits and Oscillations of stars (PLATO) mission will begin its four-year nominal mission in early 2027 by monitoring its first Long-duration Observation Phase field at South (LOPS2), for at least two years continuously. The primary aim of PLATO is a very ambitious and challenging one: the discovery of Earth-like planets in the habitable zone of nearby and bright solar analogues. To this purpose, the PLATO Mission Consortium, through its Ground-based Observing Program (GOP), will perform the follow-up needed to confirm part of the candidate planets photometrically detected by PLATO and measure their masses through radial velocity curves. For the LOPS2, the GOP is committed (as part of the PLATO mission) to follow-up the candidate exoplanets discovered orbiting the 15,000 high-quality target subset of the PLATO Input Catalog (PIC) known as the Prime Sample (PS). The PS will be made public nine months before launch in the context of the first Guest Observer (GO) call for proposals to be issued by the European Space Agency (ESA). Here, we present the quantitative metrics and thresholds defined to select and prioritize the PS. Our method is perfectly general and suitable to rank any list of stars surveyed for transiting planets. We also describe the astrophysical properties of the LOPS2 PS, both in a statistical sense and for some specific targets of interest.

This chapter provides an almost comprehensive guide to the Imaging X-ray Polarimetry Explorer (IXPE) data analysis. The chapter briefly introduces the IXPE spacecraft and the instrument onboard; subsequently, the strategies adopted to extract polarimetric information and to optimize the response are given. Moreover, the data formats and processing steps to avoid potential systematic errors and to achieve the best results from the IXPE data are reported. Both the model-independent and the model-dependent analyses are summarized, as the instrument response functions and their different available flavors. The idea behind this chapter is to collect suggestions and answers to questions that are typically raised by users, aiming to enable researchers to maximize the scientific return from IXPE observations.

The ESA PLAnetary Transits and Oscillations of stars (PLATO) mission is designed to detect terrestrial planets in the habitable zones of solar-type stars. Owing to telemetry constraints, the selection of PLATO targets must be performed in advance. In this paper, we present the first public release of the PLATO Input Catalogue of targets (tPIC2.2), which provides the list of stars that will be observed during the PLATO first Long-duration Observation Phase field at South (LOPS2) as part of its core science program. We exploit astrometric and photometric data from Gaia Data Release 3 (DR3), together with three-dimensional maps of the local interstellar medium, to identify stars belonging to the PLATO stellar samples as from mission requirements. The tPIC comprises 217,741 stars, including 202,315 FGK dwarfs and subgiants, 15,037 M dwarfs and 789 known planet host stars. The median distances of the samples are 512 pc for FGK stars and 133 pc for M dwarfs. We estimate interstellar reddening for almost all targets and develop an algorithm to infer fundamental stellar parameters (effective temperature, radius, and mass) in an homogeneous way from astrometric and photometric observables. The tPIC fulfills all the science requirements of the PLATO mission. The tPIC also includes a list of stars that host known exoplanets (confirmed or still candidate), located within the LOPS2 field.

It has been suggested that "Little Red Dots" (LRDs) might be accreting black holes enshrouded by dense gas in a nearly closed geometry, which completely covers the central black hole, leading to an atmosphere-like structure known as the "black-hole star" ($\rm BH^\star$). We test this scenario by analysing new JWST spectroscopy in the far ultraviolet (FUV, rest-frame) of the prototypical LRD Abell2744-QSO1, at $z=7.04$. We found the presence of broad Ly$\alpha$ emission with an FWHM of $\sim 1000$ km/s, and detections of OI, CIV, and/or FeII emission lines. The NIRCam imaging and NIRSpec slit images indicate that the low-velocity component ($v\lesssim 200$ km/s) of Ly$\alpha$ is likely spatially extended, but the high-velocity component ($v\gtrsim 200$ km/s) of Ly$\alpha$ remains unresolved. Based on the multi-component kinematics and flux of Ly$\alpha$ relative to Balmer lines, we conclude that the observed line profile is unlikely to be broadened by subsequent resonant scattering through the interstellar medium. This suggests that the high-velocity component of Ly$\alpha$ originates in the broad-line region, although resonant scattering in the dense gas likely makes Ly$\alpha$ broader than H$\alpha$ as observed. The nebular features of this LRD indicate that there is at least one relatively optically thin direction where Ly$\alpha$ can escape from the broad-line region (BLR). We also found indications that photons from the BLR are powering fluorescence of FeII and OI on a larger physical scale. The FUV features thus challenge the fully-covered geometry interpretation and suggest that there are "holes" in the $\rm BH^\star$, or the absorbing medium is simply clumpy.

Contemporary all-sky surveys have observed thousands of extragalactic transients in the nearby universe, and upcoming surveys will discover exponentially more at higher redshifts. With these large samples, population-level analysis of the photometric behavior of different transient classes is now possible, allowing for photometric classification and physical parameter inference from relatively sparse individual light curves. To enable such studies, we introduce Gaussian process Optimized Photometric Regression of Extragalactic Archival Ultraviolet-infrared eXplosions, a.k.a GOPREAUX--a Python package for Gaussian Process Regression of multi-wavelength transient photometry. Our modeling is unique in that it interpolates transient emission across phase and wavelength in a non-parametric, data-driven way. This allows for predictions of light curves and spectra at higher redshifts, where the rest-frame ultraviolet (UV) emission is redshifted into the observer-frame optical or infrared. To this end, we aggregate a sample of almost 1,300 transients observed in the UV and optical with the Neil Gehrels Swift Telescope, complemented with additional optical and infrared coverage from surveys such as ZTF and open-source data releases. Our sample includes 275 Type II SNe, 172 stripped-envelope SNe, 72 superluminous SNe, and 58 tidal disruption events, among other classes. Our code and reduced photometry--comprising over 146,000 photometric observations--are available as open-source software and data products. Here we discuss our sample criteria, data reduction and modeling methodologies, the multi-wavelength light curves and spectral templates produced by our models, and the future directions in photometric classification and physical parameter inference this code and data repository enables.

High-resolution X-ray spectroscopy provides a powerful probe of the interstellar medium (ISM), giving direct access to the composition and physical state of dust grains and atomic species in dense environments. We present a study of the gas and dust along the line of sight to the bright low-mass X-ray binary GX 340+0, which samples higher-density gases in the ISM. Using deep Chandra/HETG spectra, we characterize X-ray absorption fine structure from dust, gas absorption lines, and the optical depths of the Si, S, and Fe K-edges. By fitting these three edges simultaneously, we reduce degeneracies in the dust composition and find that amorphous olivine dominates the fractional contribution among the dust columns ($\sim$65%), followed by metallic iron ($\sim$19%), iron sulfides (pyrrhotite and troilite; $\sim$10%), and fayalite ($\sim$5%), with the remaining species contributing only a few percent in total. From the inferred stoichiometry, we estimate that $\sim$74% of Fe is associated with silicates, $\sim$8% with sulfides, and $\sim$18% with metallic iron, suggesting that Fe is predominantly incorporated in iron rich silicate grains along this sightline. We detect S ii absorption and infer a sulfur dust fraction of $\sim$35%. We also detect absorption structure near the Ca and Ar K edges, highlighting the need for improved atomic photoabsorption data. The Chandra/HETG spectral resolution remains essential to disentangle dust and gas contributions at the Si and S K edges, providing a benchmark for dust characterization in high-density regions in the ISM.

We investigate whether Big Bang nucleosynthesis (BBN) remains compatible with the Covarying Coupling Constants plus Tired Light (CCC+TL) cosmology. In this framework, only quantities with explicit length dimensionality covary through a universal scaling function $f \left( z \right)$, while dimensionless constants and dimensionless ratios remain invariant. At the redshifts $z$ relevant to BBN, $f \left( z \right )$ approaches a constant plateau $f_{\text{max}} \left( z \right) \simeq 3$, and the tired-light contribution is negligible, so the early-time dynamics reduce to a global rescaling of dimensioned quantities. In particular, the Hubble expansion rate $H$ at fixed temperature $T$ satisfies $H_{\text{CTL}} \left( T \right) = f^{-1}_{\text{max}} H_{\Lambda\text{CDM}}\left( T\right)$, implying a longer cooling time $\Delta t$ between weak freeze-out and the onset of nucleosynthesis by the same factor (CCC+TL labeled as $\textit{CTL}$). We find that BBN predictions are preserved provided the relevant interaction rates $\Gamma$ and decay rates governing the neutron lifetime ${\tau}_n$ share the same plateau scaling as $H$, so that governing combinations such as $\Gamma\text{/}H$ and $\text{exp} \left( -\Delta t \text{/} {\tau}_n \right)$ remain invariant. Implementing these plateau rescalings in the Kawano/NUC123 network (via a single control parameter $\texttt{fctl} \equiv f_{\text{max}}$) yields identical light-element abundances for $\texttt{fctl}= 1$ ($\Lambda$CDM) and $\texttt{fctl} = 3\left( \text{CCC+TL} \right)$ to within $10^{-3} - 10^{-4}$ level, consistent with numerical rounding. We also illustrate that adopting the lower late-time CCC+TL baryon density from the Pantheon+ data fit can reduce the ${}^7\text{Li}$ discrepancy but simultaneously increases D/H, implying that BBN alone does not select between the late-time baryon-density inferences considered here.

PSR J1928+1815 is a 10.55 ms millisecond pulsar in a 3.6 hr orbit with a massive ($1.0$-$1.6\,M_{\odot}$) companion that produces extended radio eclipses. The companion, proposed to be a stripped helium star, is undetected in optical and infrared surveys. We present deep near-infrared imaging using Keck/NIRC2 with laser guide star adaptive optics. No source is detected at the pulsar position down to a $5\sigma$ limit of $K_s \approx 21.3$. Using stripped-star atmosphere models and conservative extinction estimates, we show that any plausible helium star companion would have been detected, ruling out this interpretation. A massive white dwarf (WD) companion remains consistent with the non-detection. We consider two possible origins for the eclipses: (1) absorption in a wind driven by a young, hot WD, and (2) material ablated from the WD by the pulsar. The former can naturally arise following Case BB mass transfer, which produces $\sim 1.2\,M_\odot$ WDs capable of sustaining winds of $\dot{M} \gtrsim 10^{-12}$-$10^{-13}\,M_\odot\,{\rm yr}^{-1}$ for $\sim 10^4$-$10^5$ yr, sufficient to obscure the pulsar at GHz frequencies. The latter requires efficient coupling of the pulsar's spin-down luminosity to the companion to drive the needed mass loss, which may be difficult to achieve. If the eclipse is powered by a WD wind, the system is likely observed in a short-lived phase; alternatively, if the companion is an older WD, the origin of the eclipsing material remains unclear. The apparent uniqueness of PSR J1928+1815 is consistent with a short detectability lifetime, though formation rate estimates remain uncertain.

Simulations predict that bars in galaxies should slow down over time. This is often attributed to the exchange of angular momentum between the bar and other regions of the galaxy, such as the outer disc and dark matter halo, which implies that galaxies with a more massive halo or disc should be able to slow down the bar more efficiently. However, observational evidence for this process has been limited. In this work, we provide observational support for the slowing down of bars as predicted by simulations. We combine bar kinematics measurements obtained with the Tremaine-Weinberg method and host galaxy mass estimates derived from Jeans anisotropic models for a sample of 30 galaxies from the MaNGA survey. We find a statistically significant anti-correlation (>4sigma) between the bar pattern speed and both the stellar and total dynamical mass, which suggests that the slowest bars reside in the most massive galaxies. However, while the slope of the best-fit line between the pattern speed and dark matter mass is negative, it is not statistically significant (2.43sigma). We also find that bars with lower pattern speeds have more extended NFW dark matter profiles with lower central densities. Additionally, we find statistically significant correlations (>3sigma) between the corotation radius and the stellar mass, dark matter mass, and total dynamical mass. Finally, we find no significant correlations that involve the dark matter fraction or R, likely due to the inherent challenges associated with measuring these specific parameters accurately.

Population inference in gravitational-wave astronomy allows us to connect individual detections to the astrophysics of compact objects and their environments. Current approaches employed for population inference with LIGO-Virgo-KAGRA data approximate evaluation of the hierarchical population likelihood via post-processing of individual-event posteriors. However, the case of the Laser Interferometer Space Antenna (LISA) will be more complex for two main reasons: the transdimensional "global fit" approach to LISA data analysis which models all signals and noise simultaneously, and the presence of both individually-resolved signals and the unresolved stochastic ``Galactic foreground" arising from the Galactic binary population, which induces a circular dependence between the resolved and unresolved systems and our ability to detect the former. These challenges are not without opportunity; LISA's data will contain every mHz compact binary in the Milky Way -- either individually or within the Galactic foreground -- with great potential for Galactic and stellar astrophysics. We therefore propose an alternative approach: direct evaluation of the full hierarchical population likelihood within the LISA global fit. We develop a statistical formalism for joint inference of individually-resolved gravitational-wave sources, an unresolved stochastic foreground, and a shared, underlying astrophysical population, present PELARGIR, a prototype GPU-accelerated population inference module for the LISA global fit, demonstrate the formalism and PELARGIR via a toy model analysis, and lay out a roadmap towards an astrophysically-motivated LISA global fit with embedded population inference. While we apply the formalism here to the population of LISA Galactic binaries, it is applicable across the gravitational-wave spectrum with use cases in pulsar timing and next-generation terrestrial observatories.

We aim to confirm the compact nature and constrain the radio spectra of candidate persistent radio sources (PRSs) associated with repeating fast radio bursts (FRBs). We performed European VLBI Network (EVN) observations at 5 and 8 GHz targeting two candidates identified in a recent VLA survey. We measured flux densities and upper limits at milliarcsecond resolution and combined them with published VLBI data at lower frequencies to derive spectral constraints. We detect a compact source associated with FRB 20190417A at 5 GHz with a flux density of $150\pm45$ uJy, while no detection is obtained at 8 GHz. The source is unresolved and has a brightness temperature $T_{\rm b} \gtrsim 10^{6-7}$ K, confirming its non-thermal nature. Combining our measurement with VLBI data at 1.4 GHz, we derive a spectral index $\alpha = -0.19 \pm 0.29$, consistent with a nearly flat spectrum. This makes FRB 20190417A only the second PRS with a spectral index constrained using VLBI data. The inferred luminosity places the source on the proposed $L_{\nu}$-|RM| relation. Including this source yields a scatter of $\sigma_\Delta = 0.65$, corresponding to $\hat{\alpha}|\epsilon| = 1.5 \pm 0.7$, consistent with forward shocks in the free-expansion phase or young pulsar wind nebulae. For the candidate PRS associated with FRB 20181030A, we report upper limits of 80 uJy at 5 GHz and 150 uJy at 8 GHz, corresponding to $L_{5\,\mathrm{GHz}} \lesssim 3.8 \times 10^{25}\ {\rm erg\ s^{-1}\ Hz^{-1}}$, and implying a steep spectral index ($\alpha \lesssim -1.2$) if the VLA emission arises from a compact component. Our results highlight the importance of VLBI in isolating compact emission from FRB engines and provide one of the few spectral constraints for PRSs at milliarcsecond resolution. The consistency of FRB 20190417A with the $L_{\nu}$-|RM| relation supports a nebular origin for the persistent emission.

The diversity of galaxy morphologies and their relations with galaxy and halo properties is fundamental to understanding galaxy formation. Cosmological simulations of representative volumes can help disentangle the origin of observed correlations, but most suffer from two main limitations that affect morphologies: an over-pressurised interstellar medium and spurious interactions between stellar and dark matter particles. We present an overview of galaxy morphologies in the COLIBRE simulations, which address these limitations and reproduce many observed galaxy scaling relations. To quantify galaxy morphology, we use four (strongly-correlated) theory-space metrics, three kinematic and one spatial. We explore how different choices and limitations affect these indicators, including luminosity- versus mass-weighting, aperture size and shot noise. Overall, we find good convergence in present-day morphologies across two orders of magnitude in mass resolution. COLIBRE predicts that kinematic morphology correlates strongly with stellar mass and colour, and that galaxies with stellar masses of $\approx(1-2)\times 10^{10}\,\mathrm{M}_{\odot}$ tend to be the most rotationally-dominated. At fixed stellar mass, the morphology of central galaxies correlates weakly with the properties of their host halo. Morphology correlates more strongly with internal galaxy properties, with more disky galaxies being more gas-rich, having higher star formation rates and exhibiting younger and more extended stellar populations. Other properties, like the mass of the most massive black hole, the fraction of stars that are accreted and stellar metallicity, also correlate with morphology, but with correlation strengths sensitive to the stellar mass of the galaxy and whether it is a central or satellite.

Rest-frame ultraviolet (UV) spectra of Little Red Dots (LRDs) often show Ly$\alpha$ emission. Along with broad Balmer emission, LRDs are expected to produce broad Ly$\alpha$ emission. However, the large column density of neutral gas invoked to explain the Balmer break should significantly redshift and further broaden the Ly$\alpha$ line, making it challenging to detect without sensitive, moderate-resolution spectra. We present ultra-deep (29 hours) G140M JWST/NIRSpec observations covering the rest-UV of two LRDs in Abell2744 from the SPURS Cycle 4 Large Program. One of our targets is Abell2744-QSO1, a gravitationally-lensed LRD at $z=7.04$ with faint UV emission (M$_{\rm UV}=-16.9$), and the other source (UNCOVER-2476) is newly-confirmed at $z=4.02$ with a very bright UV continuum (M$_{\rm UV}=-19.6$). We find that Abell2744-QSO1 has a broad Ly$\alpha$ profile, along with narrow CIV, FeII$\lambda1786$, and OI$\lambda1302$ emission. The Ly$\alpha$ profile suggests an origin similar to the broad H$\alpha$, but the line is considerably less redshifted than expected from existing dense gas models. We show that the line profile can be explained if the dense neutral gas is clumpy, allowing Ly$\alpha$ to escape by scattering off of the clump surfaces. We find that UNCOVER-2476 has narrow [NeIV] emission, indicating either a hard radiation field or shocks. We confirm two close neighbors with Ly$\alpha$ emission around Abell2744-QSO1, indicating it traces a dense environment that may have ionized its surrounding IGM. We suggest that LRDs may preferentially trace bubbles carved by their dense environments, contributing to the prevalence of Ly$\alpha$ in the population.

Context. Terrestrial exoplanets are expected to host secondary, high-metallicity atmospheres derived from outgassing of volatiles such as N2, CO2, H2O, CH4, and CO. Photochemical organic hazes are likely to form in such environments, significantly affecting atmospheric observations and planetary habitability. Aims. We investigate haze formation in representative terrestrial exoplanet atmospheres and assess how CH4 versus CO as the primary carbon source affects haze production rates, particle properties, and chemical complexity. Methods. We performed six laboratory simulations by exposing gas mixtures at a few mbar to glow discharge at 300 K. Each atmosphere contained 75% N2, CO2, or H2O, 10% of each of the other two gases, and 5% CH4 or CO. Gas-phase products were analyzed with a residual gas analyzer, and solid products were characterized by production rate, particle density, atomic force microscopy, Fourier-transform infrared spectroscopy, and very high-resolution mass spectrometry. Results. CH4 experiments produced more diverse gas-phase species and much higher haze yields than the corresponding CO experiments. CO-derived hazes showed a narrow particle size range of 10-80 nm, whereas CH4-derived hazes were denser and chemically more complex. The identified molecular formulas suggest growth pathways linked to gaseous precursors such as HCN, CH2O, and C2H4. Conclusions. The atmospheric redox state critically controls haze formation in simulated terrestrial exoplanet atmospheres. CH4 is significantly more effective than CO in initiating organic growth, leading to higher haze production rates and greater chemical complexity. These results provide useful constraints for exoplanet atmospheric modeling and spectral interpretation, and further support the possibility that reducing atmospheres may facilitate prebiotic organic chemistry relevant to the emergence of life.

Moreton waves are rare chromospheric signatures of large-scale coronal disturbances, often associated with big flares and coronal mass ejections (CMEs). Using high-cadence, full-disk H$\alpha$ spectroscopic observations from CHASE, together with the EUV data from SDO/AIA and magnetograms from SDO/HMI, we analyzed a coronal EUV wave and an H$\alpha$ Moreton wave event associated with a filament eruption on 2024 July 29. The Moreton wave fronts are roughly cospatial with the fast-mode coronal EUV wave fronts, which propagate with a speed of $\sim$600 km s$^{-1}$. By comparing the Moreton wave fronts with photospheric features, we found that they preferentially appear along photospheric supergranule boundaries characterized by 1600 Å bright ridges, concentrated magnetic fields, and convective downflows. It is shown that the H$\alpha$ line profiles at the Moreton wave fronts are systematically redshifted. Gaussian fit yields a systematic downward Doppler velocity of 1.73 km s$^{-1}$. Using the bisector method, we further derived height-dependent Doppler velocities in the chromosphere. While there is an expected tendency for the downward velocity to decrease from 4.12 km s$^{-1}$ in the upper chromosphere to 1.60 km s$^{-1}$ in the lower chromosphere, it is intriguing to see an unexpected velocity enhancement in the lower chromosphere. We conjecture that coronal fast-mode MHD waves experience mode-conversion to slow-mode waves, which propagate along magnetic field lines of the magnetic canopy, forming preferential appearance of Moreton waves at magnetic networks, where the convective downflow contributes to the velocity enhancement in the lower chromosphere.

Scalar-induced gravitational waves (SIGWs) provide a powerful probe of inflationary dynamics on scales far smaller than those accessible to the cosmic microwave background and large-scale structure. In scenarios with a transient ultra-slow-roll (USR) phase, the curvature power spectrum can be strongly enhanced on small scales, potentially generating an observable stochastic GW background. In this regime, scalar dynamics during inflation can become nonlinear, challenging the validity of standard perturbative predictions. Existing semi-analytical calculations of SIGWs rely on the linear evolution of inflation fluctuations. In this work, we compute SIGWs from USR inflation using lattice simulations. We evolve the inflaton field fully nonlinearly during inflation and extract the curvature perturbation nonperturbatively, then simulate its post-reheating horizon re-entry by evolving the Newtonian potential linearly while retaining the full non-Gaussian structure of the initial conditions for the primordial fluctuations in the tensor source. For moderate non-Gaussianity, the semi-analytical prediction captures the correct order of magnitude of the GW signal but receives important corrections. When inflationary non-Gaussianities are large, it can fail dramatically in both amplitude and spectral shape, independently of the overall size of the tensor power spectrum. Our results show that reliable predictions of SIGWs in such scenarios require nonperturbative control of the inflationary scalar dynamics. The code used for this work is available at this https URL.

Probing dark matter distribution on sub-galactic scales is essential for testing the Cold Dark Matter ($\Lambda$CDM) paradigm. Strong gravitational lensing, as one of the most powerful approach by far, provides a direct, purely gravitational probe of these substructures. However, extracting cosmological constraints is severely bottlenecked by the mass-sheet degeneracy (MSD) and the unscalable nature of manual and neural-network modeling. Here, we introduce LensAgent, a pioneering training-free, large language model (LLM)-driven agentic framework for the autonomous physical inference of mass distributions. Operating as an autonomous scientific agent, LensAgent couples high-level logical reasoning with deterministic physical modeling tools, demonstarting successful reconstruction of mass distribution in SLACS Grade A strong lensing systems. This self-evolving architecture enables the robust extraction of sub-galactic substructures at scale, unlocking the cosmological potential of upcoming wide-field surveys such as the Rubin Observatory (LSST) and Euclid.

The extended main sequence (eMS) and extended main sequence turnoff (eMSTO) phenomena have been observed in some young and intermediate-age star clusters in the Milky Way and in the Magellanic Clouds. In this study, we conduct a survey of 53 galactic open clusters (OCs) to investigate the roles of stellar rotation, differential extinction, and cluster properties in the emergence of eMS and eMSTO. The projected rotational velocities are taken from the Gaia ESO spectroscopic survey and the Gaia DR3 line-broadening velocities. Stellar members of each OC are identified using the ML-MOC algorithm with Gaia DR3 astrometry. We divide clusters into four classes based on the color-rotation distribution, extinction, and MSTO morphology and report 14 clusters ($A_{\rm v} \lesssim 0.15$ mag, Class I) that exhibit split MS with fast and slow rotators populating the redder and bluer parts of MSTO. For the remaining clusters, differential extinction hampers the color-rotation distinction and also inflates MSTO width and therefore introduces a systematic offset in the MSTO-age relation. We also quantify the fraction of slow rotators among MSTO stars, finding a median value of $f_{\rm slow\, rot}^{v \sin i<100} \approx 0.41$ and the fraction reaching the spin-down limit, $f_{\rm slow\, rot}^{v \sin i<30}$, is $ \approx 0.08$. We find no statistically significant correlation between $f_{\rm slow\, rot}$ and either the binary fraction or cluster age.

We study the performance of the flat $\Lambda$CDM model and the dynamical dark energy parameterizations $w_0$CDM and $w_0w_a$CDM, in which the dark energy (DE) equation of state is either constant ($w=w_0$) or redshift-dependent [$w(z)=w_0+w_a z/(1+z)$], without and with a varying CMB lensing consistency parameter $A_L$, using combinations of Planck PR4 CMB data (PR4 and lensing), and a compilation of non-CMB data composed of baryon acoustic oscillation (BAO) data that do not include DESI BAO data, Pantheon+ type Ia supernova observations, Hubble parameter measurements $H(z)$, and growth rate $f\sigma_8$ data. We also compare results from earlier Planck PR3 data with those obtained using PR4 data in order to assess the stability of cosmological constraints. For the largest data combinations, PR3/PR4+lensing+non-CMB, the cosmological parameters inferred from PR3 and PR4 data are consistent, almost all differing by $1\sigma$ or less. For the $\Lambda$CDM$+A_L$ model, we have $A_L=1.087 \pm 0.035$ for PR3 and $A_L=1.053 \pm 0.034$ ($1.6\sigma$ above unity) for PR4, which indicates that the CMB lensing anomaly is reduced when PR4 data are used. For the $w_0 w_a$CDM parameterization, we find $w_0 = -0.863\pm0.060$ (quintessence-like) and $w_0+w_a=-1.37^{+0.19}_{-0.17}$ (phantom-like), suggesting that the current observations favor dynamical DE over a cosmological constant at about $1.8\sigma$. For the $w_0w_a$CDM$+A_L$ parameterization, we find $w_0=-0.877\pm 0.060$ and $w_0 + w_a =-1.29_{-0.17}^{+0.20}$, corresponding to a preference for dynamical DE over a cosmological constant of about $1.5\sigma$ and with $A_L = 1.042 \pm 0.037$ exceeding unity at $1.1\sigma$. These results indicate that while the PR4 data mildly favor a time-evolving DE, part of this preference may be associated with possible residual excess smoothing present in the Planck PR4 CMB anisotropy spectra (abridged).

The initial distance to PEARLSDG estimated from the Tip of the Red Giant Branch suggested it was an exotic isolated quiescent dwarf galaxy. We combine recent and archival Hectospec spectroscopy to place it at $z = 0.02843 \pm 0.00012$ ($D \approx 124$\,Mpc) within a galaxy group, revising the distance from 30\,Mpc to $\sim$124\,Mpc. We then carry out {\sc Prospector} SED fitting using parametric and non-parametric star-formation histories sampled with \texttt{dynesty}, \texttt{nautilus}, and \texttt{emcee}, recovering metallicity $\log(Z/Z_\odot) = -0.44^{+0.35}_{-0.06}$, stellar mass $\log_{10}(M_*/M_\odot) = 9.25^{+0.02}_{-3.73}$, and dust attenuation $\hat{\tau}_V = 0.67^{+0.02}_{-0.05}$. The updated metallicity places PEARLSDG squarely on the standard mass--metallicity relation, resolving its former outlier status, with its quenched star-formation history consistent with environmental quenching in a group setting.

We present the discovery of a low-mass companion located at $\rho$ $\sim$ 0\farcs{}85 ($r_{\rm proj} \approx 62~au$) from the early-type 1.2 Gyr-old star HIP 53005 using direct imaging data from the Subaru and Keck Telescopes and astrometry from the Hipparcos-Gaia Catalog of Accelerations. The companion, HIP 53005 C, is a component of a multiple system also including a $\approx$ 12\farcs{}4-separation M dwarf companion inducing a negligible proper motion acceleration. HIP~53005 C's position on color-magnitude diagrams, the fit of its spectral energy distribution to atmosphere models, and its location on an empirical mass-magnitude diagram all suggest that it lies at the M/L transition and near the hydrogen-burning limit ($\sim80~M_{\rm Jup}$). However, our orbital fitting combining direct-imaging relative astrometry with proper motion acceleration favors a much higher dynamical mass of $\sim185\ M_{\rm Jup}$. An additional unseen, more closely-orbiting companion below the detection limit (at $\rho\lesssim0\farcs2$)) may explain this discrepancy. Alternatively, HIP~53005C could be a low-mass binary like Gliese~229Bab, making this system an intriguing laboratory for studying multiple star formation.

We present constraints on the Hubble constant ($H_0$) derived from the observed dispersion measure (DM) distribution of unlocalized fast radio bursts (FRBs). While localized FRBs with redshift measurements have been used to investigate the Hubble tension, their sample remains limited. Here we demonstrate that unlocalized FRBs -- which are far more numerous -- can independently constrain $H_0$ without requiring redshift information, as cosmic expansion imprints itself on their DM distribution. Analyzing a selected sample of 2124 unlocalized FRBs from the CHIME Catalog II, we obtain $H_0 = 73.8^{+14.0}_{-12.3}~\mathrm{km\,s^{-1}\,Mpc^{-1}}$ at the $1\sigma$ confidence level, corresponding to an uncertainty of about 18%. Breaking the degeneracy between $H_0$ and the characteristic cutoff energy $E_*$ of the FRB isotropic energy distribution would reduce this uncertainty to 9%. This work constitutes the first $H_0$ measurement derived solely from the DM distribution of unlocalized FRBs, highlighting their potential as a new cosmological probe. Future joint analyses with localized FRBs promise even tighter constraints.

Reconstructing the assembly history of the Milky Way requires precise constraints on the dynamical age of its merger remnants -- the time elapsed since a progenitor satellite was disrupted by the Galactic tidal force. We present a new framework to derive this dynamical age for disrupted stellar systems by extending the Fourier analysis of the orbital frequency distribution proposed by Gomez and Helmi. To overcome the smearing of frequency-space structures caused by observational noise, we introduce the Greedy Optimistic Clustering algorithm. This method allows for an optimistic exploration of the density contrasts in the orbital frequency space by taking into account the observational uncertainty in the data, effectively sharpening the signal required for age estimation. By applying this method to the Helmi stream, we derive a dynamical age of $6.8 \pm 0.8$ Gyr. Our derived accretion epoch is consistent with the observed kinematic properties of the Helmi stream. In particular, the marked asymmetry in the vertical velocity distribution -- where approximately two-thirds of the stars have negative $v_z$ in the solar neighborhood -- supports a relatively recent arrival. This suggests that the progenitor of the Helmi stream was accreted during an epoch of Galactic growth distinct from the much earlier Gaia-Sausage-Enceladus merger ($\sim 10$ Gyr ago). We validate our methodology using error-added mock simulations, demonstrating the reliability of our approach. Our results establish the Greedy Optimistic Clustering framework as a powerful chronometric tool for reconstructing the hierarchical assembly of the Milky Way using current and future high-precision astrometric datasets.

T Tauri stars, in more advanced stages of evolution, during the final accretion phase of stellar formation, exhibit intense stellar winds and surface magnetic fields with intensities around a kilogauss. With the growing interest in the search for rocky exoplanets with Earth-like dimensions, it is essential to deepen our understanding of the interaction between stellar winds and planetary magnetospheres. We investigated the interaction between stellar winds from 46 weak-lined T Tauri stars (WTTSs) and the magnetospheric protection of Earth-like planets located within their habitable zones. We employ two distinct stellar wind models, nonmagnetized and magnetized with both constant and resonant Alfven wave damping, to evaluate the pressure balance between the stellar wind and the planetary magnetic field. Our results show that the strong wind dynamic and magnetic pressures characteristic of WTTSs lead to systematically compressed planetary magnetospheres, significantly smaller than that of the present-day Earth. The analysis further indicates that planetary magnetospheric sizes increase with stellar age, following the decay of stellar magnetic activity, in agreement with previous findings for solar-type stars.

Gauss--Bonnet inflation typically affects the dynamics over an extended portion of the trajectory, making it difficult to isolate a controlled imprint at CMB scales. We consider a trajectory-selective coupling \(\mu(\phi,X)\) that gates the Gauss--Bonnet sector in phase space, enabling the higher-curvature contribution to be localized within a finite e-fold window while remaining negligible elsewhere. We identify stable inflationary solutions consistent with this localization and enforce standard ghost and gradient stability conditions for both scalar and tensor perturbations. For these viable backgrounds we compute pivot-scale observables and examine their dependence on the overall Gauss--Bonnet strength and on the kinetic gating. The framework offers a controlled route for realizing localized higher-curvature effects with predictable consequences for CMB-scale measurements.

Flaring activity from the jets of active galactic nuclei has been studied for several decades, closely related to the loading and evolution of the jet. In this work, we focus on the sub-hundred parsec jet region, well traced by non-thermal radio and gamma-ray emission. Only in recent years have light curves capturing the decade-long behavior of such sources become available for a large ensemble of objects. While previous studies have focused on a direct correlation or few-month lag between gamma-ray and radio activity, recent neutrino-bright blazars observed by the IceCube Neutrino Observatory present multi-year delays between initial gamma-ray activity and subsequent radio flares. In this work, we search for similar-timescale correlations between Fermi-LAT gamma-ray data and RATAN-600 radio data from ~100 blazars. We consider two gamma-ray bands, 100 MeV-1 GeV and 1 GeV-500 GeV, as well as the integral band, to compare correlations between potential opaque and unabsorbed regions of the jet. Gaussian process modeling is used for smooth light curve function prediction. We also analyze morphological AGN core data from the MOJAVE survey, forming a sub-selection to better illustrate potential dependence on location. In the broader selection, several sources exhibit delayed flares on the order of 1-3 years. In the stacked analysis, we find the highest correlation for a radio delay on the order of 180 days. The stacked correlation resulting from the MOJAVE sub-selection corresponds to a slightly smaller time lag. Delayed radio flares or extended radio emission appear to be notable features within the general blazar population.

We assess whether intracluster light (ICL) can serve as an observational discriminator of dark matter physics. The self-interacting dark matter (SIDM) model has gained increasing attention as a possible resolution to small-scale discrepancies between collisionless cold dark matter (CDM) simulations and observations, predicting distinct tidal interaction histories within galaxy clusters. We analyze Cluster-EAGLE zoom-in galaxy clusters re-simulated from identical initial conditions in both CDM and SIDM frameworks. The morphological similarity between dark matter and multiple baryonic tracers -- gas, all stars, galaxies, and the combined brightest cluster galaxy plus ICL (BCG+ICL) -- is quantified using the Weighted Overlap Coefficient, a contour-overlap statistic. We find that dark matter is traced most accurately by BCG+ICL, followed by gas, all stars, and galaxies. The BCG+ICL component remains a robust tracer even at high redshift, while gas initially traces dark matter poorly but improves over time, eventually approaching the performance of BCG+ICL. Notably, in the SIDM case the gas distribution more closely resembles dark matter than in CDM. This reflects the underlying physics: in CDM, collisionless dark matter behaves similarly to the collisionless BCG+ICL, whereas in SIDM, self-interactions introduce an effective collisionality, making dark matter evolve more like the gas component. We also find that dwarf and satellite galaxies are more sensitive to the underlying dark matter model, despite their poorer overall tracing performance. Our results demonstrate the potential of ICL as a novel observational probe of dark matter physics and provide a first step toward using diffuse cluster light to constrain the nature of dark matter.

In our previous analysis of the eclipse timing variation patterns of eclipsing binaries located in and near the Northern Continuous Viewing Zone (NCVZ) of the TESS space telescope, 135 hierarchical triple star candidates were found. Now, two additional years of TESS observations are available and, hence, we have extended the former analysis with the use of the new observational data. We now detect 168 triple star candidates in the updated and reanalyzed sample. The majority ($\sim74\%$) of them are identical to the former triples candidates. For many of them, our new solutions are more certain than the original ones. Therefore, we can now conclude that we have identified at least 66 short-period hierarchical triple stellar systems in the NCVZ with full confidence. In the case of the majority of the remaining systems in our sample, the presence of a close third stellar component appears to be very likely. We also identify additional, longer timescale period variations in 34 systems ($20\%$ of the total sample) and conclude that in at least three systems the presence of a fourth stellar component is quite plausible. Finally, we report the complete disappearance of the eclipses in two former EBs and detect eclipse depth variations in seven other EBs as well. We interpret this effect as the consequence of changing orbital inclination caused by a non-coplanar third body.

We introduce a numerical optimization method to extract the X-ray reverberation response functions from the multi-band light curves of the active galactic nuclei. This approach does not require prior assumptions about the accretion disc and corona geometry, provided that the light curves result from the superposition of direct and singly-convolved signals, consistently across all bands. By reformulating the light curve equations into the matrix form, the optimal response matrix is derived by minimizing the squared difference between the observed and reconstructed light curves using a gradient-based optimization algorithm. We demonstrate that the method can robustly accommodate up to two convolution processes, such as the reverberation and the propagation, simultaneously. When tested on the synthesized light curves, the method demonstrates robustness of the solutions to variations in the relative contributions of each light curve component as the recovered response kernel remains acceptably close to the ground truth, as evaluated by both the response geometry and the reconstructed light curves. The method's tolerance to random noise was also assessed. With appropriate denoising, the response kernel can be reliably recovered when the signal-to-noise ratio is at least $100$. We show, as a proof of concept, that the proposed method is geometrical-model independent and has the potential to offer a flexible complement to traditional approaches.

The evolution of solar surface magnetic fields is essential for understanding solar activity and the underlying dynamic process. The surface flux transport (SFT) model is a widely used and effective tool for simulating this evolution. Active regions are incorporated as magnetic flux sources of the SFT model, but their configurations are usually simplified as symmetric or asymmetric bipolar magnetic regions (BMRs). Here, we aim to quantitatively and systematically assess how such flux source approximations affect SFT results and explore improved approximation methods using our recently developed SFT code. By comparing simulations that incorporate realistic active region configurations from solar cycle 23 through the ongoing cycle 25, we show that approximating active regions as symmetric BMRs leads to a systematic overestimation of the axial dipole strength at solar minimum. This result is independently confirmed using an algebraic quantification that evaluates the axial dipole contribution of individual active regions. The overestimation can be partially reduced by monotonically decreasing the size of the approximated BMRs, but it cannot be fully eliminated. When active regions are instead represented by morphologically asymmetric BMRs, the simulated axial dipole strength exhibits a strong and nearly linear negative dependence on the size ratio between the following and leading polarities. Based on these results, we propose a combination of BMR size and polarity size ratio that yields an axial dipole evolution comparable to that obtained with fully incorporated realistic active region configurations. This study provides a new quantitative constraint for improving future simulations with approximated BMRs.

In order to meet the science goals of the PLATO space mission, an extensive science calibration and validation plan has been designed. This paper describes this plan, as well as the methodology adopted to select the science calibration and validation stars that have entered its input catalogue. This is the so-called {\tt scvPIC}, which is part of the general PLATO Input Catalogue (PIC) for the first selected long pointing field in the Southern Hemisphere known as LOPS2. While many of PLATO's science requirements needed dedicated stars as calibrators as discussed here, its most stringent requirement is the delivery of the age of the host stars of exoplanetary systems with an accuracy better than 10\% for a G0V star of {\it V} = 10 mag, i.e. a nearby Sun-like star. This is presently not within reach for large populations of dwarfs and subgiants in the Milky Way as it requires the models of their stellar interiors to be improved. We discuss how this ambitious age requirement led to the selection of tens of thousands of red giants, and of thousands of main-sequence early F-type gravity-mode pulsators in order to deduce their internal rotation profile across stellar evolution. This asteroseismic observable will then be imported as key information into improved models of dwarfs and subgiants in the Milky Way as optimal modelling tools for ever better age-dating of the exoplanet hosts as the PLATO mission moves along. Additional calibrators and validators included in the {\tt scvPIC} are a few thousands of binaries, a few hundreds of legacy and benchmark stars, a few hundred photometrically stable stars, and six transiting brown dwarfs.

Recent cosmological observations, including the latest Dark Energy Spectroscopic Instrument (DESI) data releases DR1 and DR2, have renewed interest in the possibility that dark energy may exhibit dynamical behavior rather than being a strict cosmological constant. In this work, we perform a fully model-independent reconstruction of the quintessence scalar field potential using Gaussian Process regression and current Hubble measurements. Instead of assuming a specific functional form for the scalar field potential, we reconstruct the quintessence potential and the corresponding kinetic energy directly from observational data. Our analysis is based on Hubble parameter measurements obtained from cosmic chronometers and the latest high-precision DESI DR2 baryon acoustic oscillation (BAO) data, together with Type Ia supernova data from the Pantheon+ compilation. Gaussian Processes provide a nonparametric and model-independent framework that allows the data to guide the reconstruction. We employ two covariance functions, namely the squared exponential and the Matern ($\nu = 9/2$) kernels, in order to assess the sensitivity of the reconstruction to the kernel choice. We further explore the impact of background cosmological assumptions by considering different priors on the matter density and spatial curvature. Finally, we compare the reconstructed scalar field potential with two theoretically motivated benchmark models: a power law potential and an exponential potential. We find that both models remain consistent with the reconstructed potential within the inferred confidence intervals.

We derive explicit, algebraic expressions for the steady-state number density of cosmic ray electrons as a function of position and energy using Green's function of the diffusion equation with energy losses for an axisymmetric distributions of the particle sources in the galactocentric radius $r$ and distance to the mid-plane $z$. The solution is obtained for a Gaussian distribution of the particle sources in $r$ and $z$ but we show that it can be used for an arbitrary spatial distribution of the sources. The accuracy of our results is about 10% or better in a wide ranges of $r$ and $z$ and particle energies. These solutions can be used in the interpretation of radio astronomical observations of galaxies, particularly in the studies of the radio luminosities for large galaxy samples, and represent a physically justifiable and efficient alternative to the assumption of the energy equipartition between cosmic rays and interstellar magnetic fields.

X-ray variability provides a powerful probe of the innermost regions of active galactic nuclei (AGNs), offering valuable insights into the accretion process and the structure of the corona. Previous studies have established a correlation between the X-ray variability timescale, black hole mass, and luminosity, forming the AGN X-ray variability plane. A possible link between the X-ray spectral photon index and X-ray variability was noted in early studies but has rarely been incorporated into subsequent analyses of the variability plane. Moreover, the limited sample sizes in earlier works have limited the robustness and universality of the X-ray variability plane. In this work, we compile a sample of 112 AGNs with 399 exposures from the 4XMM-DR14 catalog and constrain the correlations between X-ray variability timescale, black hole mass, luminosity, and photon index using the recently developed fitting method, BADDAT {(Baseline-Aware Dependence fitting for DAmping Timescales)}, which enables a robust exploration of an extended parameter space. Our analysis confirms the dependence of the rest-frame variability timescale ($\tau_{\rm rest}$) on black hole mass ($M_{\rm BH}$) and further incorporates the photon index ($\Gamma$) into the variability plane, yielding a best-fit relation of $\log (\tau_{\rm rest}/{\rm s}) = 1.22\log (M_{\rm BH}/M_\odot) - 0.24\Gamma - 3.53$, which is strongly favored over the model with $M_{\rm BH}$ alone. In contrast, the inclusion of luminosity does not produce a comparable improvement. The correlation with $\Gamma$ likely reflects the effects of Comptonization and the geometry of the corona.

We present a Gaussian-process (GP) emulator for the monopole of the redshift-space halo power spectrum in $\Lambda$CDM cosmologies with massive neutrinos. The emulator is trained on 1000 COLA simulations distributed in a Latin-hypercube design over the six-dimensional cosmological parameter space $\{\Omega_m h^2,\Omega_b h^2,\Omega_\nu h^2,\sigma_8,h,n_s\}$, with outputs at 11 snapshots spanning $0.5 \le z \le 2.0$. From redshift-space halo catalogues we measure shot-noise-subtracted monopole spectra over $0.01 \le k \le 0.50\,h\,\mathrm{Mpc}^{-1}$. We also generate 1000 fixed-cosmology realizations to estimate the covariance matrix and to construct synthetic data vectors for likelihood tests. On held-out cosmologies, the emulator reproduces the simulated spectra to typically better than $2\%$ across the scales and redshifts considered. Combined with its GP-based estimate of interpolation uncertainty, this speed and accuracy make the emulator well suited to repeated likelihood evaluations in Markov Chain Monte Carlo analyses. The resulting framework provides an efficient route toward neutrino-mass inference from DESI-motivated redshift-space clustering measurements.

We investigate the origin of the broad exponential wings observed in a significant fraction of the Halpha profiles of JWST-discovered little red dots (LRDs) and little blue dots (LBDs). Recent studies have shown that exponential broad-line profiles are not a prerogative of LRDs, are often also present in LBDs, and need not imply that electron scattering is the dominant broadening mechanism in every source. Motivated by our unification picture in which LRDs are the dust-reddened, high-inclination counterparts of compact blue broad-line AGNs, we model the broad Balmer emission with a virialized, radially stratified broad-line region (BLR). In this framework, the observed profile is the luminosity-weighted superposition of clouds spanning a range of radii and therefore a range of characteristic virial velocities. We show that such a stratified BLR can reproduce the extended exponential-like wings observed in three representative LRDs, without requiring electron scattering to be the primary origin of the broad wings. Our results support a picture in which the broad wings and the line cores encode different physics: the wings arise primarily from virial BLR stratification, whereas the cores retain additional imprints of absorption and radiative transfer in dense gas. The successful fits further suggest that the cloud radial distribution peaks near the dust sublimation radius, while the exponential wings are shaped by the line-emitting inner BLR shells where the higher virial velocities produce the high-velocity tails. This offers a simple physical explanation for the exponential wings of little dots, without invoking exotic new components or scenarios.

We report the detection of weak high-velocity Na$\;$I absorption at V$_{\odot}$ = $-61.0\pm0.2$ km s$^{-1}$ in the spectrum of the stellar remnant at the center of the Galactic supernova remnant of 1181 AD. This velocity is not unlike that seen in old, more evolved SN remnants, but is much less than the remnant's $\simeq10^{3}$ km s$^{-1}$ expanding optical nebula. We briefly discuss its possible nature and origin.

We present a search for strong H$\beta$+[O III] emitters at $z=9.4-11.3$ in the HUDF using ultra-deep JWST/MIRI F560W imaging from the MIDIS survey. Three galaxies are identified via pronounced F560W flux excesses, consistent with strong rest-frame optical line emission. SED modelling yields rest-frame H$\beta$+[O III] equivalent widths of $\sim 600-1300$AA (median $\simeq 1260$AA), placing these sources among the most extreme known at these epochs. Combining these with a literature sample of 16 spectroscopically confirmed galaxies at $z\geq 9$, we find a median ${\rm EW}^{\rm H\beta+[O III]}_{\rm rest}\simeq 1300$AA, similar to values at $z\sim6-9$. We find no evidence for either a strong increase or decline in EW beyond $z\sim9$. A tentative trend of higher EW with increasing UV luminosity is observed, while no statistically significant anti-correlation with stellar mass is found. We place a first constraint on the H$\beta$+[O III] luminosity function at $z\simeq9-11$ ($\Phi\sim10^{-3.4}\,{\rm Mpc^{-3}\,dex^{-1}}$ at $\log( L_{\rm H\beta+[OIII]}/{\rm erg\,s^{-1}})=42.5$), consistent with a decline relative to $z\sim7-8$. The MIDIS sources have $\log(\xi_{\rm ion}/{\rm Hz\,erg^{-1}})=25.1-25.4$. We find significant correlations between $\xi_{\rm ion}$ and EW and $\beta$, but not with UV luminosity, consistent with trends at lower redshift. These results suggest that the physical conditions governing nebular emission and ionising efficiency are already in place at $z\sim9-11$, extending trends established at $z\sim6-9$.

Bayesian inference of nanohertz gravitational-wave background models in pulsar timing array analyses often relies on Gaussian-process interpolators to avoid repeated, computationally expensive strain-spectrum calculations. However, Gaussian-process training becomes a bottleneck for large training sets. We test whether probabilistic neural networks can replace Gaussian processes in this role for both a self-interacting dark matter model and a phenomenological environmental model. We find that neural networks recover consistent posteriors while significantly reducing both training and Markov chain Monte Carlo runtime, with the largest gains for the more computationally demanding model.

We examine the impact of input neodymium (Nd) atomic data on the light curves and spectra of kilonovae, probing the sensitivity of kilonova observables to the atomic physics of this important lanthanide element. We use the SuperNu Monte Carlo radiative transfer code, simulating a simple semi-analytic 1D kilonova with a pure Nd atmosphere, fixing the radiative transfer method while using input atomic data generated by three different codes: the LANL suite of atomic physics codes, HULLAC, and Autostructure. We see that the choice of atomic data significantly shapes the resulting light curves and spectra. Peak bolometric luminosities differ by a ratio of nearly 1.5 between HULLAC/Autostructure and LANL data sets. Moreover, we observe significant near- to mid-IR differences in the structure of the spectra. We specifically attribute these differences to the choice of atomic data for neutral Nd I. Many of the results here have been adapted from a presentation at "Radiative Transfer and Atomic Physics of Kilonovae" in Stockholm, 2023. We additionally present a LANL data set with energies calibrated to available values in the NIST Atomic Spectra Database, and demonstrate that this calibration also significantly affects IR spectral structure at late time. The substantial differences in kilonova observables that arise from tuning the atomic data of just one lanthanide element highlight the special attention that must be paid to atomic physics uncertainties when modeling kilonovae, from AT2017gfo to beyond.

Gaia two-parameter (G2P) stars have cumulative errors in parallax and proper motion so great that only their mean positions were reported in DR3. One potential cause of these high errors is another star as indicated by two intensity peaks in the scans. Speckle interferometry astrometric measurements of 25 G2P stars with high multi-peak percentages were obtained with the 1.5m telescope at Mt. Wilson Observatory. Of the 25 observed G2P stars, seven had no reported Gaia companions within 5.0". We found nearby companions for all seven. The 18 other G2P stars had known Gaia companions within 2.0". Of these, 13 had separations that agreed closely with the speckle measurements but with some discrepancy in position angles, three stars did not agree in either separation or position angle and no companion was detected for the remaining two. Although some of these issues may be resolved in DR4 or DR5, others may be inherent limitations of Gaia capabilities that speckle interferometry observations may be able to fill.

Using $\Lambda$CDM and Padé-(2,1) cosmography, we study directional variations in the Hubble constant, $H_0$, using galaxy cluster and Type Ia Supernovae (from Pantheon Plus) by the hemisphere decomposition method. Since there is a degeneracy between $H_0$ and absolute magnitude $M_B$ for Supernovae, Cepheid host calibration is usually required to constrain $H_0$. Hence, in this work in order to complement the Cepheid host calibration in Supernovae, we also use calibrations based on galaxy cluster scaling relations. We find that there is a $\lesssim 1\sigma$ difference in $H_0$ variations when using galaxy clusters as calibrators compared to Cepheids highlighting that the variations in $H_0$ are robust across different calibration methods. Across all combinations of models and data sets used, we obtain a consistent deviation $\sim 2\sigma$ from isotropy. In nearly all cases, we notice that the maximum $\Delta H_0$ aligns with the CMB dipole direction.

The Hubble tension is shaped not only by shifts between early- and late-time parameter estimates, but also by the stiffness of the constraints that define them. In this work, we analyze this geometric structure in the wCDM model by separating the discrepancy into two components: a parameter displacement and a directional Fisher curvature. Within the local Gaussian approximation, the quadratic tension along a given direction factorizes into the squared shift and the combined directional curvature contributed by the datasets. Applying this framework to Planck, DESI DR2, and SH0ES, we show that extending \LambdaCDM to wCDM primarily reshapes the Fisher geometry of the CMB constraint rather than opening a genuinely new route to concordance. Allowing the dark-energy equation-of-state parameter w to vary suppresses the leading Planck Fisher eigenvalue to only \sim 2.7 % of its \LambdaCDM value, while producing only a modest rotation of the dominant acoustic-scale eigenmode. The net effect is a strong softening of the effective acoustic rigidity. At the same time, high-precision late-time data, especially from DESI DR2, inject substantial curvature along the expansion-rate direction. This added stiffness acts as a geometric wall, closing off phantom-like escape routes and sharply limiting tension relief within the extended parameter space. Our results indicate that changes in the inferred H_0 tension under model extension are best understood as a reconfiguration of the constraint manifold rather than as evidence for new physical agreement. The shift-curvature decomposition thus offers a simple, fast, and physically transparent way to diagnose cosmological tensions.

We present a model-independent test of anisotropy in the cosmic distance duality relation (CDDR), $D_L=(1+z)^2 D_A$, using the Pantheon+ type Ia supernova sample and baryon acoustic oscillation (BAO) data. The angular diameter distance is reconstructed via Gaussian Processes, enabling an estimate of $\eta(z)=D_L/[D_A(1+z)^2]$ without assuming a background cosmology. We also allow for a possible isotropic evolution, parameterized as $\eta(z)=1+\eta_1 z$, and find a redshift-dependent deviation whose significance depends on the assumed supernova calibration. Anisotropy is modeled through a dipole modulation and constrained using a full covariance-based likelihood. To assess statistical significance, we construct null realizations that preserve both the redshift distribution and the survey selection function. We find that the observed dipole amplitude is consistent with isotropic expectations and lies below the levels induced by statistical fluctuations and survey geometry. We obtain a robust 95\% upper bound $A_{95}=0.025$, stable across different supernova calibration choices. We find no evidence for intrinsic anisotropy in the CDDR. Our results highlight the importance of accounting for survey selection effects in anisotropy searches and provide a viable framework for testing directional deviations in cosmological relations.

We aim to investigate the feasibility of accurately determining the helium-to-metal enrichment ratio $\Delta Y/\Delta Z$ for open clusters using Gaia DR3 photometry. To test the reliability of this calibration, we performed a theoretical investigation using mock open clusters. We generated synthetic photometric data from isochrones calculated by five different stellar evolution codes (FRANEC, PARSEC 1.2s, PARSEC 2.0, BASTI, and MIST), for which the true $\Delta Y/\Delta Z$ is known. We then fitted these mock clusters with two sets of isochrones calculated with the FRANEC code, differing only in the implementation of bolometric corrections (BCs). The analysis focused on the G-band absolute magnitude range (4.3 to 6.5 mag) to minimise the impact of poorly constrained physics. Synthetic clusters were generated at [Fe/H] values from 0.0 to 0.15 dex, for different numbers of populating stars and different levels of photometric uncertainties. The Monte Carlo experiments revealed significant and code-dependent biases. Unbiased results were achieved only when the stellar models used for synthetic-cluster generation and fitting were identical. Using identical FRANEC stellar models but different BCs introduced a significant bias of up to 0.6. Furthermore, using different stellar models for synthetic cluster generations resulted in even larger biases: $\Delta Y/\Delta Z$ was underestimated by up to 0.8 for PARSEC target isochrones, while it was overestimated for BASTI and MIST isochrones by up to 0.6 and 1.5, respectively. The magnitude and the inconsistency of these biases strongly suggest that the photometric calibration of $\Delta Y/\Delta Z$ using open clusters is not reliably robust.

We infer recoil (kick) velocities for all binary black hole merger events reported up to the GWTC--4 catalog, together with candidate intermediate-mass black hole events. We obtain informative kick constraints for GW231028\_153006 ($839^{+1018}_{-681}\,\mathrm{km\,s^{-1}}$) and GW231123\_135430 ($974^{+944}_{-760}\,\mathrm{km\,s^{-1}}$). Additionally, we compute recoil velocities for recently reported events from the ongoing fourth observing run: GW241011\_233834, GW241110\_124123, and GW250114\_082203, obtaining $v_{\rm kick} = 974^{+555}_{-466}\,\mathrm{km\,s^{-1}}$, $394^{+582}_{-207}\,\mathrm{km\,s^{-1}}$, and $115^{+301}_{-95}\,\mathrm{km\,s^{-1}}$, respectively. The remnant of GW241011\_233834 is therefore inferred to have one of the largest recoil velocities among currently known events. We find that present recoil kick constraints are driven primarily by measurements of the mass ratio and spin magnitudes, while the contribution from spin orientation angles remains subdominant in most cases. We estimate typical retention probabilities of the remnant black holes in GWTC catalogs to be $\sim 1$--$5\%$ for globular clusters, $\sim 15$--$30\%$ for nuclear star clusters, $\sim 5$--$40\%$ for dwarf galaxies, and $\sim 70$--$100\%$ for elliptical galaxies. We further show that, even for remnants retained in globular clusters, recoil-induced spatial displacements from the cluster core are often significant, which can substantially suppress the chances of hierarchical mergers. We find that the probability for a GWTC merger remnant to participate in hierarchical mergers is $\sim 0.1$--$1\%$ in globular clusters and $\sim 1$--$15\%$ in nuclear star clusters.

Compact, steep-spectrum radio sources are key tracers of exotic astrophysical objects such as pulsars and high-redshift radio galaxies. All-sky radio surveys at different frequencies, like the TIFR-GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS), have been usually exploited to identify such tracers. The more recent imaging survey, Rapid ASKAP Continuum Survey (RACS), with higher angular resolution and better sensitivity offers an avenue for a far better identification and characterization of compact, steep-spectrum sources. In this work, using publicly available RACS images at 887 MHz and 1.4 GHz, we present an image-domain characterization of 171 compact source candidates between declinations -40 degrees and +41 degrees, that were detected and appeared compact at 147 MHz in TGSS but not detected at 1.4 GHz in NVSS. Our detailed characterization resulted in the identification of 66 compact sources, 87 non-compact, diffuse or resolved sources, and 18 sources that are not detected in either of the RACS or NVSS images, implying spectral indices steeper than -2.0. Out of the 66 compact sources, 34 have spectral indices steeper than -1.5. We demonstrate that a large fraction of the sources in our sample were earlier not detected and resulted in incorrect spectral index limits due to poor imaging quality of NVSS in the Galactic plane. We present the spectral indices and morphological classification of all the sources in our sample and discuss their usefulness in identifying and studying interesting sources such as radio pulsars, high-redshift radio galaxies, and other extragalactic sources.

Power-law (PL) distribution functions (DF) are prevalent in highly diverse systems. The systems range in size from nanometer to mega light years, in complexity from dust grains to living organisms, and characterize the distribution of various events in nature and in various human activities. To gain some insight on why PL DF are so prevalent, we explore the conditions leading to the formation of a PL DF in a simple system of colliding hard sphere. We follow the time evolution of the energy DF through direct Monte Carlo simulations. In statistical equilibrium, the DF evolves into the Maxwell-Boltzmann (MB) DF. A transition to a PL DF occurs when: 1. The system is initially far from equilibrium. For example, a mix of light and heavy particles with the same velocity. 2. The system dynamics is scale-free, which holds in the intermediate asymptotic regime, far from the initial and the final equilibrium states. The scale-free dynamics leads to a DF which evolves in a self-similar form. 3. The system is open with a scale-free boundary condition. For example, a constant injection of particles far from equilibrium. The DF PL index is set by the time dependence of the self-similar DF and by the boundary condition. The PL index is independent of the self-similar DF form. Conditions 1-3 are common in a great variety of systems, which may explain why PL DF are so prevalent in nature.

China has systematically collected nighttime astronomical plates since 1900, creating a large historical dataset that has been digitized with optical scanners. For astrometric registration of these digitized plates, sources were first extracted using SExtractor, and then matched astrometrically with this http URL and the Gaia catalog. However, suboptimal early storage conditions and subsequent environmental deterioration have impeded accurate source matching, resulting in processing failures for several thousand digitized plates. In this work, we introduce a Transformer-based classification model that takes cutouts of SExtractor-detected sources as input and leverages multi-scale feature fusion to identify trustworthy stellar sources on the plates. Trained on plates with successful astrometric calibration, our AI-based classifier was then applied to SExtractor detected sources of 1883 digitized plates, enabling us to complete the astrometric registration for 1353 of them. This AI-augmented pipeline streamlines the processing of historical plate archives and enhances their scientific value for long-term time-domain astronomical studies.

Cygnus X-3 has recently been established as a variable ultra-high-energy(UHE) gamma-ray source with photons detected up to 3.7~PeV. The temporal correlation between its PeV activity and GeV flares, together with the possible orbital modulation, suggests that the emission is produced within or close to the binary system. In this work, we test whether the contemporaneous GeV emission zone can also host the acceleration of the parent protons responsible for the multi-PeV photons. We jointly model the contemporaneous \textit{Fermi}-LAT spectrum and orbital light curve with a one-zone leptonic scenario dominated by anisotropic external inverse-Compton scattering. The fit places the GeV emission region at $H\sim2.8\times10^{11}\,$cm and constrains the magnetic field--size product to $BH\lesssim10^{13.3}\,$G\,cm at the $3\sigma$ level. This implies a maximum proton energy of only $\sim0.3$~PeV from the Hillas criterion, far below that required by the observed PeV emission. We therefore conclude that the GeV zone cannot be the main PeV acceleration site. Instead, the PeV emission should originate from a more compact inner region, and the jet magnetic field must dissipate rapidly between the PeV and GeV emitting zones.

Using DESI SV3 spectroscopic group centrals together with deep HSC photometric data, we measure the conditional luminosity functions (CLFs) of central and satellite galaxies, separately for red and blue populations, in dark matter halos spanning $M_h\sim10^{12}- 10^{15}M_{\odot}$ and $0<z<1$. The depth of the HSC imaging enables CLFs to be measured to unprecedentedly faint limits, reaching $M_r \approx -15$ at $0.2 \leqslant z < 0.5$ and $M_r \approx -17$ at $0.5 \leqslant z < 1.0$. We find a remarkably weak evolution in the CLF of satellite galaxies over $0<z<1$. The Blue satellite CLFs are well described by a single Schechter function across all halo masses and redshifts, with a nearly constant slope of $-1.25\lesssim \alpha\lesssim -1.2$. In contrast, red satellite CLFs exhibit a pronounced and ubiquitous faint-end upturn in all halo mass and redshift bins, with little evolution in the faint-end slope ($-1.8\lesssim \alpha_f\lesssim -1.7$). These results indicate that the low-mass end of the red sequence in clusters/groups was already established by $z\sim1$. Both satellite characteristic magnitudes and central galaxy luminosities fade with time. Red central galaxies are consistent with passive evolution, whereas the luminosity evolution of blue centrals is dominated by ongoing star formation. Satellite galaxies evolve more rapidly than predicted by simple stellar population models, highlighting the importance of environmental effects. The quenched fraction of satellite galaxies as a function of stellar mass exhibits a universal minimum at $M_{\ast} \sim 10^9M_{\odot}$, independent of halo mass and redshift. We discuss possible interpretations of these results and their implications for galaxy formation and evolution.

This Research Note is in response to the recent paper by S. M. Perriyil et al. (2026). They provide measurements of the time delay (delta t) between the hard X-ray and soft X-ray peak times for 96 flares observed with RHESSI and GOES. These delays are found to be dependent on the length of 9 the magnetic loop(s) joining the HXR footpoints seen in RHESSI images. We offer a possible explanation for this coincidence in terms of the duration of the electron beam heating, commonly inferred from the duration of the HXR emission, and the time taken for heated plasma to rise to the loop top as inferred in this paper from delta t. We suggest that the particle acceleration occurs at or near the top of the loop(s) and that it is quenched by the increase in density as the heated plasma reaches the acceleration site.

We present an independent pipeline for detecting candidate vanished sources on digitized first-epoch Palomar Observatory Sky Survey (POSS-I) photographic plates. The pipeline detects and PSF-filters sources on POSS-I Red DSS cutouts, applies local astrometric registration refinement, and identifies candidates by cross-epoch matching against POSS-I Blue and POSS-II Red with Pan-STARRS DR1 rejection. On a 20-case benchmark harness, the pipeline recovers 8/9 sources in the April 1950 field and 3/3 in the July 1952 field, with a false positive rate of about 0.2 per 10 arcmin field on random non-crowded controls. A full-footprint sweep over the POSS-I coverage using 30 arcmin patches yields a filtered catalog of 2.85 million candidate vanished sources after post-processing PSF cuts, deduplication, and Pan-STARRS DR1 rejection. Cross-matching against the published 5,399-source Solano et al. (2022) catalog yields 3,450 matches (63.9%) with median separation 0.94 arcsec; among unrecovered catalog entries within our footprint, we find no Pan-STARRS DR1 counterpart within 3 arcsec. Applying Bruehl and Villarroel (2025)-style temporal windows to this catalog over the 368 POSS-I observation nights in the 1949-1957 interval gives a post-test calendar-day relative risk of 1.35 for the +1 day window, but the effect is not statistically significant (95% CI 0.91-2.00; two-sided Fisher p = 0.17) and is sensitive to coding unobserved days as zero-transient days. A negative binomial model of nightly candidate counts with nightly patch coverage as exposure is likewise null (IRR = 1.03, 95% CI 0.89-1.18, p = 0.71). The catalog-level replication is strong; the temporal association remains inconclusive.

Galactic cosmic rays (CRs) are a fundamental non-thermal component of the interstellar medium (ISM). Understanding the transport of super-high-energy particles is essential for interpreting observations of Galactic PeVatrons. Classical diffusion models assuming a homogeneous and isothermal medium oversimplify the multiphase ISM. We utilize high-resolution 3D MHD simulations to self-consistently generate a multiphase ISM, comprising the warm (WNM), unstable (UNM), and cold neutral medium (CNM), and investigate 1.5-15 PeV particle transport using a test-particle approach. We find that thermal phase transitions induce steep magnetic field strength gradients at phase boundaries, creating localized magnetic fluctuations that act as efficient sites for adiabatic mirror reflections and non-adiabatic pitch-angle scattering, strongly enhancing cross-field transport at these interfaces. However, because phase boundaries occupy only a small volume fraction and particles spend most of their trajectory in the weakly scattering WNM and UNM, the global pitch-angle scattering coefficient in the multiphase ISM is smaller than in an equivalent isothermal medium. This locally strong scattering nevertheless drives both parallel and perpendicular spatial diffusion coefficients to $\sim10^{30} {\rm cm^2 s^{-1}$ at 1.5~PeV, with the perpendicular component exceeding its isothermal counterpart ($\sim 10^{28}{\rm cm^2 s^{-1}$) by two orders of magnitude. Using a phase--phase diffusion matrix decomposition, we show that global CR transport is governed by the volume-filling, trans-Alfvénic WNM and UNM, where particles stream along stochastically wandering field lines. Cross-phase displacement correlations are universally positive, indicating cooperative transport between thermal phases. In contrast, the super-Alfvénic CNM acts as an efficient confinement that substantially suppresses local diffusion.

The PLATO mission is scheduled for launch early 2027. In this paper we present an overview of the performance drivers for the mission at the time where all flight models of the cameras have been tested and integrated on the optical bench. The PLATO consortium needs an estimate of the planet detection yield to dimension the ground-based radial velocity follow-up resources. We provide updated estimates on the yield of planet detections that can be expected from the mission under certain assumptions. As of today, large uncertainties remain on the planet occurrence rates, especially for small planets in long-period orbits, and on our ability to detect these planets in the presence of stellar variability and instrumental noise. To partially overcome these limitations, we compare results using different planet occurrence rates, detectability rates, and we include an estimate on the expected contribution of stellar variability to the noise budget. The final detection yield of PLATO will provide constraints to planet occurrence rates which in turn will help constraining planet formation models.

Using the hydrodynamical (fluid) approximation, we present a self-consistent theoretical framework that couples the origin, evolution and decay of spiral structures to the secular dynamical evolution of their host galactic disks. Our approach highlights non-resonant spiral excitation through azimuthal forcing that leverages mild, pervasive gradients in the disk's mass and angular momentum distributions, structural features we term cavernae. These cavernae are weaker but more widespread than the sharp features behind groove mode excitation and commonplace in exponential disks. We discuss how non-resonant features combine with other responses -- resonant dressing, steady waves, groove modes -- to produce a global, evolving spiral nexum that transports angular momentum and reshapes the disk. Using expressions for torques, angular momentum transport and heating, we demonstrate that global spirals are intrinsically self-limiting; the angular momentum changes and heating they generate quenches their own growth, dictating a finite lifetime for any single spiral episode. A succession of transient episodes, each with properties adjusted to the changed disk conditions, lays the pathway to long-lived spiral activity. This behavior suggests that the character of secular evolution shifts over time. We find that the short-lived, high-multiplicity (high-m) spirals that dominate in dynamically cold disks induce widespread, impulse-like non-resonant heating, yet with a low ratio of heating to radial migration. As the disk warms, high-m features are suppressed, leading to steadier, lower-m spirals that heat progressively more efficiently near resonances. In this light, the dynamical coldness of disk galaxies today requires a past dominated by high-m transient perturbations, whereas warmer, more compact systems reflect an advanced stage of evolution regulated by transient, low-m spirals.

Tidal disruption events (TDEs) can be observed when stars get too close to supermassive black holes and are torn apart and accreted. The delay time distribution of TDEs, or rate of TDEs as a function of time since a burst of star formation, can be used to determine what mechanisms influence the TDE rate. We compile a catalog of 41 TDE host galaxies with optical spectra, model the stellar populations with Bagpipes, and retrieve the age of the most recent burst of star formation to construct the delay time distribution of TDEs. TDEs occur more frequently in post-starburst galaxies than in other types of galaxies, though the mechanism causing this rate enhancement is unknown. We find that the TDE rate increases with post-burst age to reach a peak at ~1 Gyr relative to a control sample. We compare the observational TDE delay time distribution to theoretical models, which propose overdense stellar nuclei, radial anisotropies in stellar orbits, supermassive black hole binaries, and AGN disks as potential mechanisms that may enhance the TDE rate in post-starburst galaxies. Most models predict a TDE rate that declines with post-burst age, in contrast to our observational results, though some models are still feasible at certain ages (e.g., the black hole binary model matches at old burst ages and the stellar overdensity model matches at intermediate burst ages).

Joint analyses of high-resolution CMB temperature maps with galaxy surveys provide a unique way to reconstruct the radial velocity field of the underlying matter distribution via the kinematic Sunyaev-Zeldovich (kSZ) effect. Using data from the Atacama Cosmology Telescope (ACT) DR6 and the Dark Energy Spectroscopic Instrument (DESI) DR2, we present radial velocity reconstructions for luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs). Leveraging the spectroscopic data, we are able to reliably model the foreground contamination and report a negligible impact on our main observables. We detect the velocity-galaxy cross-correlation at $17.0\sigma$ for LRGs, and for the first time, at $8.3\sigma$ for ELGs and $6.8\sigma$ for QSOs. We further report the first detection of the velocity-velocity correlation using LRGs at $3.1\sigma$, as well as the highest cumulative detection of the kSZ effect to date at $20.8 \sigma$. Similarly to previous results, we find a lower amplitude of the kSZ signal compared to our fiducial halo model prediction and electron profile assuming a Battaglia profile. Combining these new observables, we obtain constraints on local-type primordial non-Gaussianity (PNG): $f_{\rm NL}^{\rm loc} = 15.9_{-34.4}^{+34.6}$ at 68\% confidence, which represents the tightest constraint to date derived from the velocity field. The measurements presented here already exhibit lower noise on a per-mode basis than the galaxy auto-correlation on the largest scales, $k<0.004~\rm{Mpc^{-1}}$, highlighting the key role these observables will play in the context of future CMB experiments such as the Simons Observatory.

At least one large-body (diameter > 1.1 km) hypervelocity cratering event occurred during ~ 0.8-0.90 Ma (Zhamanshin, Kazakhstan) in the Middle Pleistocene Transition period. Analysis designed to reduce uncertainty in the dimensions of the Zhamanshin structure employing high resolution topography demonstrated that it likely generated a ~ 26.5 km diameter multi-ring crater. This is at least two times larger than the current best estimates. Using a range of accepted impactor sizes, velocities, compositions, and angles of impact, such impacts typically yield kinetic energies of impact over 240,000 Megatons (TNT). Explosive energetic events of this magnitude (e.g., Yellowstone Caldera) at other times (K-Pg) have created global environmental effects. The factor of two discrepancy in the dimensions of Zhamanshin increases the kinetic energy yield by factors of 7-10, with significantly larger environmental consequences. This justifies examination of rapid climate transitions linked to biological consequences, including those related to environmental perturbations, at ~0.9 Ma.

Fast Radio Bursts (FRBs) carry fossil information about non-Gaussianity generated during inflation. This primordial signal is most accessible on the largest scales, where the scale-dependent bias correction $\propto f_\mathrm{NL}\,H_0^2/k^2$ dominates, but where systematic effects are also strongest. A central challenge is the degeneracy between the intergalactic-medium electron bias $b_e$ and the primordial non-Gaussianity (PNG) signal, which can degrade $\sigma(f_\mathrm{NL})$ by orders of magnitude when $b_e$ is marginalised. We show this degeneracy can be broken internally by exploiting the cross-power spectrum $C_\ell^{D\Delta t}$ between the FRB dispersion measure (DM) field and Shapiro timing delays along multiple interferometric sightlines. The DM field traces the biased electron density, while the Shapiro timing signal probes the Newtonian gravitational potential independently of astrophysical bias. Their cross-correlation is directly proportional to $b_e$, independently of the matter power spectrum, providing a self-calibration of the electron bias. We derive $C_\ell^{D\Delta t}$ analytically in the Limber approximation and find a correlation coefficient $|\rho(\ell)|\approx 0.51$--$0.79$ across $\ell = 2$--$100$. A joint Fisher matrix analysis over $\{f_\mathrm{NL},\,b_e^0,\,z_\mathrm{fb}\}$ shows that including the cross-spectrum reduces $\sigma(b_e^0)$ by a factor of $2.1$--$5.1$ relative to a DM-only analysis. After full marginalisation, the joint analysis recovers $\sigma(f_\mathrm{NL})$ within a factor of $1.0$--$1.9$ of the fixed-bias benchmark, compared with $1.7$--$3.3$ degradation without the cross-spectrum. For a shallow survey with a 500\,AU baseline and $10^4$ FRBs, the joint constraint achieves $\sigma(f_\mathrm{NL})\approx 790$, within 4\% of the fixed-bias result and a factor $3.3$ better than the marginalised DM-only case.

We discuss flat Friedmann-Lemaitre-Robertson-Walker (FLRW) metric-affine cosmology where the metric and connection as well as the matter energy-momentum and hypermomentum all obey the symmetry of spatial homogeneity and isotropy. In particular, we outline a scenario where a dark dust fluid carries spin hypermomentum which makes its effective equation of state dynamical and might relate to the DESI DR2 data.

We perform fully relativistic GRMHD simulations of magnetized, self-gravitating black hole-disk (BHD) systems in which the black hole spin is misaligned with the disk angular momentum. Massive disks (disk to BH mass ratios of $16-28\%$) around rapidly rotating black holes ($\chi\lesssim 0.97$) develop a nonaxisymmetric instability for tilt angles from $0^\circ$ to $180^\circ$. Magnetic stresses damp, but do not completely suppress, the nonaxisymmetric instability, and corresponding gravitational wave (GW) emission, in aligned systems, while they enhance it in antialigned BHDs: MRI-driven turbulence enhances angular momentum transport and accelerates nonlinear instability evolution in misaligned configurations. All models launch magnetically driven jets consistent with the Blandford-Znajek (BZ) mechanism, with collimation depending on spin orientation. The GWs reflect strong nonaxisymmetric structure from a persistent $m=1$ mode. The coupling between fast MRI and the slower nonaxisymmetric instability growth governs the outcome, with tilt controlling how MRI modifies the global mode. These simulations provide the first self-consistent GRMHD treatment of tilted, self-gravitating BHD systems and support their role as multimessenger sources.

This paper reports on the design, fabrication, and demonstration of a silicon photonics based heterodyne interferometric imaging system. The photonic integrated circuit (PIC) can perform one-dimensional spectroscopy for unique input spectrums using a single baseline within its 91 available baselines. The PIC uses polarization diversifying gratings to separate incoming light into two distinct polarizations, an on-chip 2x4 optical hybrid, and a strong local oscillator (LO) to perform the heterodyne measurements. The optical hybrids combine the input signals with the LO and splitting them into 2 components pairs for phase sensitive measurements. Furthermore, the PIC can perform 2-D image reconstruction by combining many baseline pairs to measure the visibility of a simple target. These demonstrations show the PIC's capabilities for 1-D spectroscopy and 2-D imaging applications.

The inflaton equation of motion including one loop radiative corrections from spectator fields is obtained. We consider a massless scalar conformally coupled to gravity and a massless fermion Yukawa coupled to the inflaton as models for spectators that \emph{do not feature} gravitational particle production, their production during slow roll is solely a consequence of their coupling to the inflaton. The one-loop self energy and the fully renormalized equation of motion of the inflaton are obtained and solved explicitly for an inflaton potential $m^2\varphi^2/2$. The solution features Sudakov-type logarithmic secular terms, which are resumed via the dynamical renormalization group and compared to the solutions with a phenomenological friction term. During $N_e$ e-folds of slow roll inflation the inflaton evolves as $\varphi^{(0)}_{Isr}(t)\,e^{\frac{m^2\Gamma}{9H^3}\,N_e(t)}$ for the phenomenological friction term $\Gamma$ and $\varphi^{(0)}_{Isr}(t)\,e^{\Upsilon N^2_e}$ with $\Upsilon = -\frac{\lambda^2}{24\pi^2 H^2} ; \frac{y^2_R}{12\pi^2}$ for the radiative corrections from bosonic and fermionic spectators respectively where $\varphi^{(0)}_{Isr}(t)$ is the slow roll solution in absence of interactions, showing that a phenomenological friction term is not reliable. A generalization of the optical theorem to a finite time domain and cosmological expansion is introduced to obtain the distribution function $f(k,t)$ and total number of spectators produced \emph{during slow roll}. $f(k,t)$ is peaked at superhorizon scales and the total number of particles grows $\propto e^{3N_e}$. A non-perturbative mean field theory is introduced to describe the self-consistent evolution of the inflaton coupled to spectators, its linearized version reproduces the self-energy, the inflaton equation of motion and the results on particle production.

Temporal offsets between Gamma-Ray Bursts (GRBs) and high-energy neutrinos provide a useful probe of propagation effects in extreme astrophysical environments. We investigate whether such offsets can be generated by photon propagation through dense axion clouds gravitationally bound to magnetars. Working within the Euler-Heisenberg effective theory extended by the axion sector, we derive the modified photon dispersion relations in the presence of a strong magnetic background and an oscillating axion field. We show that axion-photon mixing turns the magnetized vacuum into an anisotropic birefringent medium, leading to geometry-dependent deviations from luminal propagation and kinematic time delays that reach $\Delta t_{\perp}\simeq1.33\times10^{-12}\,\mathrm{s}$ for orthogonal propagation. Although this effect is many orders of magnitude larger than the delays expected in diffuse astrophysical backgrounds, it remains far too small to account for the macroscopic offsets discussed in current multimessenger candidates. We further show that the same birefringent medium constrains the survival of the intrinsic linear polarization of prompt GRB emission, yielding the environmental bound $g_{a\gamma\gamma}\lesssim6.02\times10^{-14}\,\mathrm{GeV}^{-1}$ for benchmark magnetar-scale parameters and axion masses near $m_a\sim10^{-4}\,\mathrm{eV}$. Magnetar-hosted axion clouds thus emerge as complementary environments in which dispersive transport and polarimetric observables jointly probe axion electrodynamics.

We present a theoretical framework demonstrating a deterministic initialization mechanism for Warm Inflation via classical conformal boundary conditions. A persistent challenge in dissipative inflationary models is the "cold start" paradox: initializing the requisite thermal bath to generate the dissipative friction that subsequently sustains radiation production. Postulating an idealized, asymptotically scale invariant pre-inflationary phase, we mathematically prove that a conformal Weyl mapping to the emergent metric furnishes a finite, analytically derived initial radiation density. Implementing a spontaneous conformal symmetry-breaking ansatz, an emergent inflaton field is subjected to this inherited thermal bath. We analytically derive the initial kinematics of this framework, demonstrating that for strict sub-Planckian temperatures, the universe naturally initializes in the weak dissipative regime (Q << 1). The initial Hubble friction provided by the boundary radiation enables a smooth, deterministic kinematic handoff to the warm slow-roll steady-state attractor. As a mathematical proof-of-concept, this mechanism provides a fully realized framework to bypass the bootstrap problem of warm inflation.

We present fully kinetic simulations of driven 2D turbulence in a relativistic plasma, designed for the first time to induce a fast magnetosonic cascade. As the driving strength increases, turbulence transitions from a weak wave-dominated regime to strong shock-driven dynamics. Using spatiotemporal Fourier analysis, we identify fast modes, finding that the weak turbulence regime exhibits spectral properties that are in excellent agreement with theoretical expectations. Our results are relevant for the modeling of turbulence in high-energy astrophysical plasmas.

We study oscillon/I-ball solutions in a real scalar version of the Friedberg-Lee-Sirlin (FLS) model. Using the two-timing analysis, we derive the conditions for oscillon solutions and explore multi-field oscillon configurations. In these configurations, the two fields form co-located oscillons that oscillate with frequencies set by their respective masses. These multi-field oscillons can be viewed as a bound state of two oscillons due to attractive interactions between the fields. We confirm these analytical predictions through numerical lattice calculations. This work extends the standard picture of single-field oscillons and may be relevant for cosmological scenarios involving multiple interacting real scalar fields.

We perform a comparative phase space analysis of the light mass Galileon model and standard Quintessence in the context of late--time cosmic acceleration. Focusing on a spatially flat FLRW background, we consider a cubic Galileon interaction supplemented by a scalar potential and examine three representative choices of the potential: a generalized cosh potential, a simple cosh potential, and a linear potential. By introducing suitable dimensionless variables, the cosmological field equations are reformulated as an autonomous dynamical system, allowing a systematic investigation of the stationary points and their stability properties. For the light mass Galileon scenario, we find that although the phase space admits scalar field dominated solutions, all critical points are of saddle type for the potentials considered. In particular, no stable late-time accelerating attractor emerges, even in the presence of de-Sitter like configurations. In contrast, the Quintessence limit admits stable de-Sitter attractors for cosh potentials, providing a viable description of the observed late--time acceleration. Our results highlight a key qualitative distinction between Galileon and Quintessence cosmologies and indicate that, within the light mass Galileon framework, the higher-order Galileon interactions may be required to realize a stable accelerating Universe.

In this article, we report the results of comparing the effect of using trace of stress-energy tensor versus real-valued scalar field in Nonminimal Derivative Coupling gravitation model, respectively denoted as NMDC-T and NMDC-phi. We employ the model into an incompressible star and see the effect of both models NMDC-T and NMDC-phi on the compactness and mass-radius relation. We find that coupling parameters of NMDC-T is less sensitive than NMDC-phi.

We extend a recently developed Big Bang Nucleosynthesis (BBN) code, {\tt PRyMordial}, to constrain a broad class of Grand Unified Theories to which BBN is sensitive, since these lead to varying fundamental couplings. A previously developed self-consistent perturbative analysis of the effects of these variations has been implemented in {\tt PRyMordial}, leading to robust constraints of the value of the fine-structure constant, $\alpha$, at the BBN epoch using current observations of Helium-4 and Deuterium abundances. We explored two different viable scenarios, relying on alternative assumptions on the gravitational sector: the variation of the gravitational coupling can be implemented by varying either particle masses, or Newton's gravitational constant. For the variation of masses, we obtained at $68\%$ confidence level a constraint on the relative variation of $\alpha$, between the BBN epoch and the present-day laboratory value, of $\Delta\alpha/\alpha=2\pm51$ ppm (parts per million), while for the variation of Newton's constant the analogous constraint is $\Delta\alpha/\alpha=2\pm22$ ppm. We also show that, given these constraints, these models do not provide a solution to the cosmological Lithium problem.

Considering $Z_4$ symmetry in Type I seesaw scenario, one could obtain mass-squared differences of light neutrinos, mixings and $CP$ violating phase within $3 \sigma$ confidence level based on neutrino oscillation data. This is possible with only three independent complex parameters for allowed Yukawa couplings and one real mass parameter for heavy right handed neutrino fields around electroweak scale. After considering only three more real parameters as coming from small soft-symmetry breaking terms, the lightest right handed neutrino could be considered as dark matter candidate via freeze-in mechanism and the other two heavier right handed neutrinos through their decays, could generate the baryonic asymmetry of the universe naturally via resonant leptogenesis.