Papers for Tuesday, Feb 17 2026

I consider recent observations of ultramassive black holes. These appear to confirm theoretical predictions that the relation between central black hole mass $M$ and spheroid velocity dispersion $\sigma$ has the same form $M \propto \sigma^4$ in spiral galaxies, elliptical galaxies, and cluster ellipticals, but has differing normalizations. These arise from the need for longer black hole accretion episodes to expel the gas otherwise potentially able to feed the holes in the latter two types of host. In a sample drawn from a mixture of galaxy host types the fitted power of $\sigma$ will slightly exceed the theoretically-derived value of 4 because of the differing normalizations. The observed hole masses do not currently reach the theoretical maximum values possible for disc accretion, set by the equality of the ISCO and self-gravity radii, probably because the host galaxies have insufficient gas.

Dark sectors consisting of atomic constituents (electrons, protons, and photons) offer a well-motivated extension to the Standard Model while providing multiple avenues for phenomenological study. In this work, we explore the impact of conversions between the dark and Standard Model photons in the primordial CMB spectral distortion epoch ($10^3 \lesssim z \lesssim 10^6$). These conversions are resonantly enhanced when the induced thermal masses of both photonic species are equal, thus leading to the possibility that sizeable distortions can be produced. To this end, we solve the Boltzmann equation at early times to determine the (irreducible) freeze-in or freeze-out abundance of dark photons. This procedure also allows us to update the limits on generic milli-charged dark sectors using the ACT DR6 bound on the number of effective radiative degrees of freedom ($N_{\rm eff}$). By then modeling the evolution of the thermal masses in both sectors, we compute the primordial CMB distortion using the Landau-Zener formalism. We find that when the dark electron and proton are roughly similar in mass (the positronium limit), current spectral distortion data from the COBE/FIRAS instrument is able to rule out novel regions of parameter space. We also forecast bounds from the proposed FOSSIL satellite, finding that spectral distortions can also be used to probe the ultra-low dark electric charge regions of parameter space, which are difficult to investigate by other means.

Globular cluster stellar streams probe galaxy-formation processes and can potentially reveal the distribution of dark matter in galaxies. In many theoretical studies, streams are modeled with particle-spray or direct N-body codes. But particle-spray methods abstract away the internal dynamics of the progenitor by making strong assumptions about the escape physics, while direct N-body is prohibitively expensive for realistic (N>10^5) systems. In this paper, we present the stream-modeling capabilities of KRIOS, a new basis-expansion N-body code for collisional stellar dynamics, that bridges this runtime vs. accuracy gap. We show that KRIOS reproduces NBODY6++GPU cluster models, and their associated streams, more accurately than particle spray in a fraction of the NBODY6++GPU wall-clock time. We then compare KRIOS to various particle-spray methods on 10 orbits similar to known Milky Way streams. The morphology and kinematics of these streams most disagree when the progenitor is tightly bound to the host, as these systems are often subject to stronger tidal forces. Finally, we discuss which elements of the progenitor physics are most important for modeling stellar streams and how these might be incorporated into particle-spray methods.

Dyad is a Python library for studying the dynamics of binary stars considered as gravitational two-body systems. The dynamics of a binary star are determined by its primary and secondary stars' masses along with the secondary star's six orbital elements as defined in a frame comoving with the primary star. In a population of binary stars these eight parameters vary from member to member and can each be treated as a random variable having some probability distribution. Dyad provides a class, this http URL, and a module, this http URL, for dealing with such a population of binary stars. The this http URL class represents a gravitational two-body system while the this http URL module provides a suite of classes representing the probability distributions of (1) stellar masses (including those proposed by Kroupa, 2001, and Salpeter, 1955) as well as (2) mass ratios and orbital elements (including those proposed by Duquennoy & Mayor, 1991, and Moe & Di Stefano, 2017). The software repository is at this https URL .

The Lyman $\alpha$ (Ly$\alpha$) line from high-redshift galaxies is a powerful probe of the Epoch of Reionization (EoR). Neutral hydrogen in the intergalactic medium (IGM) can significantly attenuate the emergent Ly$\alpha$ line, even in the damping wing of the cross-section. However, interpreting this damping wing imprint relies on our prior knowledge of the spectrum that escapes from the galaxy and its environs into the IGM. This emergent spectrum is highly sensitive to the composition and geometry of the interstellar and circumgalactic media, and so exhibits a large galaxy to galaxy scatter. Characterizing this scatter is further complicated by non-trivial selection effects introduced by observational surveys. Here we build a flexible, empirical model for the emergent Ly$\alpha$ spectra. Our model characterizes the emergent Ly$\alpha$ luminosity, the velocity offset of the Ly$\alpha$ line with respect to the systemic redshift, and the H$\alpha$ luminosity, with multivariate probability distributions conditioned on the UV magnitude. We constrain these distributions using $z\sim5-6$ galaxy observations with VLT MUSE and JWST NIRCam, forward-modeling observational selection functions together with galaxy parameters. Our model results in Ly$\alpha$ equivalent width distributions that are a better match to (independent) Subaru observations than previous empirical models. The extended distributions of Ly$\alpha$ equivalent widths and velocity offsets we obtain could facilitate Ly$\alpha$ transmission during the early stages of the EoR. We also illustrate how our model can be used to identify GN-z11-like outliers, potentially originating from merging systems. We publish fitting functions and make our model publicly available.

The landmark multi-messenger observations of the binary neutron star (BNS) merger GW170817 provided firm evidence that such mergers can produce short gamma-ray bursts (sGRBs). However, the limited number of BNS detections by current gravitational-wave (GW) observatories raises the question of whether BNS mergers alone can account for the full observed sGRB population. We analyze a comprehensive set of 64 BNS population synthesis models with a Monte Carlo-based framework to reproduce the properties of sGRBs detected by Fermi-GBM over the past 16 years. We consider three jet geometry scenarios: a universal structured jet calibrated to GW170817, a universal top-hat jet, and a non-universal top-hat jet with distributions of core opening angles. Our results show that models characterized by low local BNS merger rates ($R_{BNS}(0) \lesssim 50$ Gpc$^{-3}$ yr$^{-1}$) predict too few observable sGRBs to reproduce the Fermi-GBM population, effectively disfavoring them as sole progenitors. Even when relaxing assumptions on jet geometry, low-rate models remain viable only for wide jets ($\theta_c \ge 15^\circ$), in tension with the narrow jet cores ($\theta_c \approx 6^\circ$) inferred from sGRB afterglow observations. In contrast, models with local merger rates of order $R_{BNS}(0) \approx 100$ Gpc$^{-3}$ yr$^{-1}$ successfully reproduce the observed sGRB population, assuming a plausible fraction of BNS mergers launch relativistic jets and realistic jet geometries. This analysis highlights the power of combining GW observations of BNS mergers with electromagnetic observations of sGRBs to place robust constraints on the BNS merger population and to assess their role as progenitors of sGRBs.

Galaxy evolution is shaped by star formation and stellar feedback at scales unresolved by current high-resolution cosmological simulations. Precise subgrid models are thus necessary, and different approaches have been developed. However, they are degenerate and often primarily calibrated to reproduce stellar masses from observations. To explore these degeneracies, we perform three cosmological zoom-in radiation-hydrodynamics simulations of the same galaxy within a $5\times10^{11}\rm\ M_\odot$ dark matter halo at $z\sim1$, each with a different subgrid model: mechanical feedback, a combination of mechanical feedback and thermal feedback, and delayed cooling. We calibrate the simulations to match in stellar mass, isolating the effect of the models on the circumgalactic medium (CGM). Our findings demonstrate that despite producing galaxies with comparable stellar masses, the three models lead to distinct feedback modes, resulting in notable variations in the CGM properties. The delayed cooling run is dominated by ejective feedback and exhibits high burstiness, whereas mechanical and the hybrid models primarily feature preventive feedback, respectively acting at the galaxy and halo scales. Delayed cooling reduces the baryon mass to half the universal baryon fraction while mechanical feedback retains most baryons, with the hybrid model standing in between. Delayed cooling also ejects significantly more metals into the CGM than both other models. While for delayed cooling and mechanical feedback metals are almost evenly distributed in the CGM, they are concentrated around satellites in the hybrid model. These discrepancies emphasize the need to design an appropriate subgrid model to understand how stellar feedback regulates galaxy growth.

Galaxy evolution is driven by star formation and stellar feedback on scales unresolved by current high-resolution cosmological simulations, requiring robust subgrid models. However, these models remain degenerate, often calibrated primarily to match observed stellar masses. To explore these degeneracies, we conduct three state-of-the-art cosmological zoom-in simulations of the same galaxy, each incorporating different subgrid models: mechanical feedback, a combination of mechanical and thermal feedback, and delayed cooling. We compare their circumgalactic media (CGM) through quasar absorption sightlines of HI, MgII, CIV, and OVI. Our findings demonstrate that despite producing galaxies with the same stellar masses, the models lead to distinct feedback modes and CGM properties. Column densities and covering fractions serve as effective diagnostics of subgrid models, with all four ions providing strong constraints as they trace diverse gas phases, exhibit complementary spatial distributions, and originate from different mechanisms. Although all simulations bracket observed column density distributions, direct comparisons are limited by scarce detections and significant scatter in absorption strengths. Covering fractions of weak absorbers provides the most robust constraints. All models fail to reproduce HI and MgII covering fractions, and delayed cooling overproduces OVI covering fractions, while the other models underproduce them. The simulation including mechanical feedback reproduces the observed CIV covering fractions well, whereas the other models show slight offsets. We argue that this discrepancy is likely driven by unresolved thermal structures for HI and MgII, and insufficient metals for CIV and OVI, arising from missing physics such as AGNs or cosmic rays.

At low redshifts, hydroxyl megamasers (OHMs) have been shown to trace galaxy mergers, obscured starbursts, high molecular gas densities, and candidate dual supermassive black hole systems. Given this astrophysical utility, exploring these sources at larger cosmological look-back times is therefore of key interest. While previous OHM surveys have been limited to redshifts of $z \lesssim 0.25$, the ability to expand the OHM frontier is significantly enhanced with new high-sensitivity radio facilities such as MeerKAT. In this Letter, we report the discovery of an OHM in the gravitational lens system HATLAS J142935.3-002836 at $z = 1.027$, the most distant OHM source yet detected. The spectrum has blended 1667 and 1665 MHz emission and exhibits a highly complex profile, with spectral components ranging in widths of $<8$ km s$^{-1}$ to $\sim300$ km s$^{-1}$. The integrated (magnification uncorrected) luminosity of log($L_{\rm OH} / L_{\odot}$) = 5.51 $\pm$ 0.67 makes this the most apparently luminous OHM known to date. In the same wide-band dataset, we have also detected a previously unknown ${\rm H I}$ absorption line. The signal-to-noise ratio of over 150 with just a 4.7 h observation highlights the potential that MeerKAT and the future Square Kilometre Array mid-frequency array offer to explore the high-redshift OHM universe.

We examine the multi-wavelength period-luminosity-color relations for Cepheid variables in the Large and Small magellanic Clouds and the Milky Way. From first-principles stellar physics, the luminosity of a Cepheid is determined by its radius and surface temperature, yielding a fundamental PLC relation whose observational proxies are pulsation period and intrinsic color. Using Cepheids in the Magellanic Clouds, we show that the PLC relation recovers the known geometries and line-of-sight tilts of their disks, confirming its ability to detect true distance-modulus variations that are achromatic and consistent across all filters. Surprisingly, for Milky Way Cepheids with individually determined reddenings and HST and Gaia parallaxes, the residuals from multi-wavelength PL fits are also found to be achromatic, identical in sign and amplitude across all passbands, in this case indicating that parallax errors are the dominant source of scatter. Applying bandpass-averaged corrections to individual Cepheids recovers the theoretically expected wavelength-dependent narrowing of the instability strip, and results in revised parallaxes with a median improvement in precision of roughly a factor of two. In addition, they show no statistically significant correlation with metallicity over the range -0.2 < Fe/H < 0.05 dex. The final extinction- and reddening-corrected PLC relation yields an rms scatter of 0.04 mag, corresponding to 2 percent precision in distance per star. Use of a physically grounded PLC will provide a more robust foundation for the Cepheid-based extragalactic distance scale and the determination of the Hubble constant.

The concept of a system-wide measure of the sustainment of life (habitability) for space-faring interplanetary species is introduced and explored. Although largely agnostic to the details of how interplanetary life might operate (e.g., via technology or by utilizing organism traits that are, as of now, unknown to us), some assumptions must be made about energy harvesting, orbital mobility costs, radiation risks, and resource requirements. A multi-modal figure of merit is developed for evaluating an interplanetary habitable zone (IHZ). An agent-based model is also developed to simulate the dispersal of interplanetary life in a planetary system and characterize the IHZ. For the solar system, resource weightings between planetary bodies dictate many overall behaviors, including the sequence of migration from Earth to the Moon, Mars, and asteroid belt. Comparisons with the Trappist-1 exoplanetary system also point to critical sensitivities in the balance between resource availability and risk or cost factors (e.g., radiation risks and orbital Delta-v costs) that determine the structure of an IHZ. Results suggest that our solar system may have an inherent, and significant, advantage for a space-faring species over a system like Trappist-1. This modeling approach may also have application to emerging space economies in our own solar system.

Re-derivation of the sensitivities of all of the Space Telescope Imaging Spectrograph (STIS) observing modes were required after major updates were introduced to the model atmospheres of the three primary standard stars. The new predicted continuum fluxes were up to 2-3% different from the models used to originally calibrate STIS. This work focuses on the re-derivation of spectral sensitivities for the five STIS low-resolution (L-mode) gratings: G140L, G230L, G230LB, G430L, and G750L, which span wavelengths from the far-ultraviolet through the near infrared. Updated photometric throughput tables were delivered to the Calibration Reference Data System (CRDS) on April 7, 2022 and April 14, 2023, which triggered a recalibration of all historical STIS datasets taken with these modes. The sensitivities derived from each of the standard stars typically agree with one another to better than 1%, though discrepancies as large as 1.5% are found in spectral regions most impacted by hydrogen absorption.

Measurements from the Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission 2 (ARCADE-2) reveal a strong radio background in the GHz frequency range. Since the cosmological 21-cm signal is measured relative to the background radiation temperature, the presence of a radio excess can significantly alter its characteristics. Previous studies have explored the impact of an inhomogeneous radio background on the global 21-cm signal and 21-cm power spectrum. This non-uniform radio background is also expected to introduce substantial non-Gaussianity. In this work, using the bispectrum, we analyze the non-Gaussianity in the 21-cm signal in the presence of an excess galactic radio background and investigate how line-of-sight radio fluctuations from early galaxies influence its nature. We find that even a moderate enhancement in radio efficiency in early galaxies significantly affects the small-scale 21-cm bispectrum. Furthermore, the delayed heating transition caused by a galactic radio background shifts the sign change in the squeezed-limit bispectrum to lower redshifts ($z\sim11$), providing a potential observational signature for distinguishing different radio background models. These results demonstrate that the 21-cm bispectrum, particularly in the squeezed limit, is highly sensitive to radio background fluctuations, making it a powerful tool for probing high-redshift radio-loud sources and the physics of the early cosmic epoch.

Protoplanetary disc mass is one of the most fundamental properties of a planet-forming system, as it sets the total mass budget available for planet formation. However, obtaining disc mass measurements remain challenging, since it is not possible to directly detect H$_2$, and CO abundance ratios are poorly constrained. Dynamical measurements of the disc mass are now possible, but they are not suited to all discs since the measurements typically require well-behaved emission surfaces. A long-standing method is to obtain continuum flux measurements from the dust emission, and convert to a total disc mass by assumption of the dust-to-gas mass ratio, $\epsilon$. This quantity is poorly constrained in protoplanetary discs. % We investigate the impact of $\epsilon$ on the morphology of planet-containing hydrodynamical simulations of dusty protoplanetary accretion discs, and suggest that if a planet mass estimate can be obtained, then disc morphology could be used to constrain $\epsilon$ in observed systems relative to each other, improving the total disc mass estimates of protoplanetary discs.

In this eighth paper of the SOFIA-FORCAST series on Milky Way GHII regions, we present an analysis of the massive star-forming complex W43 Main. We compared our 11 - 37 micron maps with multi-wavelength observations from the near-infrared to radio, and investigated the physical nature of compact sources and dust substructures. We applied SED fitting to constrain properties of the compact infrared objects, and examined the evolutionary states of the extended subregions. We identified 20 compact infrared objects, 16 (80%) of which we classify as massive young stellar objects (MYSOs) or candidate MYSOs. W43 Main resides at the junction of the Scutum spiral arm and the Galactic Bar, a location where enhanced turbulence is anticipated and has been proposed as a potential influence on star-formation activity. Nevertheless, our analysis shows that its Lyman continuum photon production rate, the mass of its most massive MYSO, and its MYSO density are all consistent with the survey-wide median values. We therefore conclude that, despite W43 Main's unique Galactic environment, its present star formation activity appears broadly consistent with that of an average Galactic GHII region.

With the forth observing run of the LIGO-Virgo-KAGRA gravitational-wave network, which enabled the discovery of the kilonova (KN) counterpart to GW170817, ending with no new confirmed neutron star mergers, the intrinsic rate of these events must be even lower than previously estimated. As a result, building a sample of KNe will remain challenging even with continued GW observations, motivating complementary discovery strategies that do not rely on gravitational-wave triggers. In this work, we consider how leveraging bright short gamma-ray burst afterglows can aid in the discovery on KNe with the Rubin Observatory's upcoming Legacy Survey of Space and Time (LSST), whose unprecedented depth will make such detections feasible. We find that nearly on-axis (${\theta}_{\rm view} \leq 30°$) afterglows can enhance KN detection rates in the LSST $g$-band from $29^{+51}_{-21} \ \rm yr^{-1}$ to $91^{+160}_{-65} \ \rm yr^{-1}$. We further show how the colors of the observed events can be used to distinguish between neutron star merger counterparts with and without KN emission. This study demonstrates how critical multi-wavelength and multi-survey observations are for these rare events, especially without context from gravitational waves. Fortunately, detectable events will likely be discovered near peak with LSST, allowing for rapid follow-up and confirmation. We discuss key uncertainties in our study, particularly volume rate of merger events, and the degeneracy between the empirically determined explosion energy and ambient medium density.

The superorbital period of LMC X-4 is among the most stable known in Roche-lobe overflow, high-mass X-ray binaries. We analyzed 33 years of monitoring data from the Compton Gamma Ray Observatory Burst and Transient Source Experiment (CGRO BATSE), the Rossi X-ray Timing Explorer All-Sky Monitor (RXTE ASM), the Neil Gehrels Swift Burst Alert Telescope (Swift BAT), the Monitor of All-sky X-ray Image Gas Slit Camera (MAXI GSC), and the Fermi Gamma-ray Burst Monitor (Fermi GBM). The measured phases show a smooth long-term trend with superposed systematic fluctuations. Fits with cubic, quartic, and sinusoidal models indicate that the quartic and sinusoidal forms provide significantly better descriptions, with the sinusoidal model yielding an $8900^{+210}_{-230}$-day modulation. Such a long timescale is unlikely to arise from orbital motion around a tertiary companion. The fluctuations resemble stochastic, glitch-like events on several-hundred-day timescales. Their rms period variation exceeds that of the smooth trend, yet the total rms period variation over 33 years remains only 0.55\%, demonstrating the exceptional stability of the superorbital period. During MJD 57000-60461, we detect a phase offset of 0.044$\pm$0.010 cycles between the soft and hard X-ray bands. This offset can be reproduced by including a higher-harmonic term in the azimuthal disk model, allowing a transition from antisymmetric to asymmetric structure. A contemporaneous decline in the hard X-ray flux suggests a partial obscuration of the emission region, similar to the anomalous low state in Her X-1. This evolving-disk scenario may also explain the superorbital phase shift previously reported in Her X-1.

This is a review article for The Review of Particle Physics 2026 (aka the Particle Data Book), appearing as Chapter 25. It forms a compact review of knowledge of the cosmological parameters near the end of 2025. Topics included are Parametrizing the Universe; Extensions to the standard model; Probes; Bringing observations together; Outlook for the future.

We present a systematic study of the internal mass structure of early-type galaxies (ETGs) based on 106 galaxy-scale strong gravitational lenses with background quasars, all having spectroscopic redshifts. From this parent sample, we select 24 systems with high-quality ancillary data for joint analysis of strong lensing geometry and stellar kinematics. A key contribution is the derivation of new single-aperture stellar velocity dispersions for 11 lens galaxies via an iterative spectroscopic fitting procedure that mitigates quasar contamination, providing previously unavailable data. We model the total mass-density profile as a power law, $\rho \propto r^{-\gamma}$, and parameterise its logarithmic slope as $\gamma = \gamma_0 + \gamma_z \cdot z_l + \gamma_s \cdot \log \tilde{\Sigma}$, where $z_l$ is the lens redshift and $\tilde{\Sigma}$ the surface mass density. Within a flat $\Lambda$CDM framework and using DESI BAO measurements as a prior, we constrain the parameters via Monte Carlo nested sampling to $\gamma_0 = 1.62^{+0.11}_{-0.12}$, $\gamma_z = -0.35^{+0.08}_{-0.09}$, and $\gamma_s = 0.37^{+0.08}_{-0.07}$ ($68\%$ confidence intervals). Our results robustly demonstrate that $\gamma$ increases with surface mass density ($\gamma_s > 0$) and decreases with redshift ($\gamma_z < 0$). This implies that, at fixed redshift, galaxies with denser stellar cores have steeper mass profiles, while at fixed density, profiles become shallower at higher redshifts. By successfully applying the joint lensing--dynamics method to a substantial, independently acquired sample of lensed quasars, this work provides crucial validation of structural trends previously observed in galaxy--galaxy lensing systems, reinforcing the established evolutionary picture for massive ETGs and establishing lensed quasars as a potent probe of galaxy structure.

Understanding turbulence in interplanetary coronal mass ejections is fundamental to space plasma research and critical for assessing the impact of space weather on geospace. Turbulence governs energy cascade, plasma heating, magnetic reconnection, and solar wind magnetosphere coupling, thereby influencing both ICME evolution and geoeffectiveness. While previous event-based and statistical studies have examined ICME turbulence and its radial evolution in great detail, no significant measurements of ICME magnetic turbulence at a specific vantage point have been made using multiple observatories separated azimuthally. Here, we present the first multipoint analysis of MHD turbulence across ICME plasma regions, using four spacecraft at the Sun-Earth L1 point, separated by 80 RE along the dawn-dusk direction. Previous studies reveal that ICME shocks, sheaths, and magnetic clouds are highly non-uniform, with strong azimuthal variability. Using high-resolution magnetic field observations from ISRO's Aditya-L1, NASA's Wind and ACE, and NOAA's DSCOVR, we analyze turbulence associated with the 10th October 2024 solar storm, which triggered the second-strongest geomagnetic storm of solar cycle 25. Our results reveal significant variability and differing turbulence maturity across small separations, supported by analysis of field-aligned and perpendicular magnetic field cascades, indicating strong anisotropies. Sheath turbulence is substantially modified by shock induced energy injection. Evidence of compressible turbulence and plasma energization at the flux rope interaction region indicates that internal processes, such as magnetic reconnection, strongly influence ICME plasma evolution, highlighting pronounced spatial variability in turbulence and plasma states observed by multiple L1 monitors near Earth and underscoring their potential role in space weather impacts.

The systematic variation of solar active region (AR) properties with their magnetic flux has been the subject of numerous studies but the proposed scaling laws still vary rather widely. A correct representation of these laws and the deviations from them is important for modelling the source term in surface flux transport and dynamo models of space climate variation, and it may also help constrain the subsurface origin of active regions. Here we determine active region scaling laws based on the recently constructed ARISE active region data base listing bipolar ARs for cycle 23, 24 and 25. For the area $A$, pole separation $d$ and tilt angle $\gamma$ we find scalings against magnetic flux $\Phi$ and heliographic latitude $\lambda$. Residuals from these relations are also modelled. These scaling relations are recommended for use in space climate research for the modelling of future data or missing past data, as well as for the identification of candidate rogue ARs. We confirm that the tilt angle distribution of non-Hale ARs shows a significant excess at low tilts (anti-Hale ARs). In contrast to earlier studies we show that neither the anti-Hale ARs nor non-Hales in general follow Joy's law: instead, their tilt angle distribution is best represented by vanishing mean tilt. These results are most easily reconciled with a scenario where the AR flux loops originate in the deep convective zone or below, gaining tilt during their rise under the action of the Coriolis force. A small fraction of the loops is subjected to extreme, intermittent torques resulting in either very large tilts or anti-Hale orientation. Anti-Hale ARs are suggested to be fully curled `XO-loops', and their excess is caused by a simple mechanical effect, as the contact of their legs increases resistance against further deformation by the torque.

We present photometric and spectroscopic observations of luminous red nova (LRN) AT 2025abao, the fourth discovered in M 31. The LRN, associated to the AGB star WNTR23bzdiq, was discovered during the fast rise following the minimum phase. It reached the peak at $g=15.1$ mag ($M_g=-9.5\pm0.1$ mag), and then it settled onto a long-duration plateau in the red bands, lasting 70 days, while it was slowly linearly declining in the blue bands. The object showed similarities at peak with the canonical LRNe V838 Monocerotis, V1309 Scorpii, and with the faint and fast-evolving AT 2019zhd, the third LRN in M31, though the later evolution is different. Spectroscopically, AT 2025abao evolved as a canonical LRN: the early spectra present a blue continuum with narrow Balmer lines in emission; at peak, the spectral continuum has cooled to a yellow colour, with a photospheric temperature of 6000 K. Balmer lines have weakened while absorption lines from metals (Fe I, Fe II, Sc II, Ba II, Ti II) have developed, and in particular broad (FWHM$\sim$700 km/s) from the UV Ca II H&K lines. Medium- and high-resolution spectra reveal a counter-P Cygni absorption profile in H$\alpha$. Finally, late time spectra show an orange continuum ($T\sim$4000-5000 K), a return in strength of the Balmer lines and the formation of molecular absorption bands. AT 2025abao is the rare case of a LRN with detailed archival information regarding the progenitor system. For the first time, we obtained the spectral energy distribution in the infrared of the precursor of a LRN, which is consistent with that of an M giant/AGB. We propose that the dichotomy of light curve behaviour in LRNe (two peaks vs. plateau) can be explained by the extent and H-richness of the common envelope.

The observations of protostellar systems show a transition in the radial profile of specific angular momentum (and rotational velocity), evolving from $j\sim{\rm constant}$ ($v_{\phi}\sim r^{-1}$) in the infalling-rotating envelope to $j\propto r^{1/2}$ ($v_{\phi}\sim r^{-1/2}$) in the Keplerian disk. We employ global MHD disk simulations of gravitational collapse starting from a supercritical prestellar core, that forms a disk and envelope structure in a self-consistent manner, in order to determine the physics of the Envelope-Disk Transition Zone (ENDTRANZ). Our numerical results show the transition from the infalling-rotating envelope to Keplerian disk happens through a jump in the $j-r$ profile over a finite radial range, which is characterized by the positive local gravitational torques. The outer edge of the ENDTRANZ is identified where the radial infall speed ($v_r$) begins a sharp decline in magnitude and $j$ begins a transition from $j\sim{\rm constant}$ toward $j\sim r^{1/2}$. Moving radially inward, the centrifugal radius ($r_{\rm CR}$) is defined where $v_{\phi}$ first transitions to Keplerian velocity at the disk's edge. Farther inward of $r_{\rm CR}$, model disk develops a super-Keplerian rotation due to self-gravity. The inner edge of the ENDTRANZ is defined at the centrifugal barrier ($r_{\rm CB}$) where $v_r$ drops to negligible values. Inside $r_{\rm CB}$, a net negative gravitational torque drives mass accretion onto the protostar. On observational grounds, we identify a jump in the observed $j-r$ profile in L1527 IRS for the first time using the ALMA eDisk data. Comparison with the numerical radial behavior from our MHD disk simulations suggests the observed $j-r$ jump can be used as a kinematical tracer for the existence of ENDTRANZ. Our results offer insights into the observable imprint of angular momentum redistribution mechanisms during star-disk formation.

Neutron stars in globular clusters follow complex evolutionary pathways shaped by binary interactions, mass transfer, and dynamical exchanges. Direct N-body simulations such as NBODY6++GPU successfully model stellar dynamics and compact object formation, but they usually do not track pulsar spin evolution or magnetic field decay explicitly. Building on Papers I and II of this series, we identify this gap and present a case study from an existing simulation with N = 105000 particles, showing how a neutron star forms and evolves for 200 Myr without any pulsar-physics tracking. We compare this situation with recent implementations and outline a seven-scenario framework that includes magnetic dipole spin-down, exponential magnetic field decay, environmental torques, accretion-driven spin-up, gravitational-wave emission, and merger-driven evolution. As an example, the neutron star we label Pulsar973 forms at t = 800 Myr with a post-supernova mass of 5.35 solar masses and evolves to 2.52 solar masses by t = 1000 Myr, but still lacks period P, period derivative Pdot, magnetic field B, and scenario classification. We provide mathematical formulations and specific integration points within NBODY6++GPU (Hermite scheme, Ahmad-Cohen neighbors, KS regularization, and BSE stellar evolution) to enable scenario-based pulsar evolution within direct N-body simulations.

GPM 1839-10 is an intriguing long-period radio transient (LPT), distinguished by its activity spanning at least three decades and its highly unusual emission characteristics. These features include orthogonal polarization mode (OPM) switches, down-drifting sub-structures, and distinct linear-to-circular polarization conversion behaviors. In this work, we present follow-up observations utilizing the FAST telescope at L-band, yielding a total of seven detected radio pulses. We find a consistent association between OPM switches and a decrease in polarized intensity. This feature strongly supports the hypothesis that the OPM switches are generated by the incoherent summation of OPMs. Our measured Rotation Measures (RMs) are consistent with previous observations, indicating that the magneto-ionic environment is stable. If the source is in a binary system, such stability suggests it may host a weakly magnetized companion. Crucially, we firstly observe clear evidence of a cyclotron absorption feature in one radio pulse, a signature rarely observed in radio sources. This feature allows us to infer that the magnetic field strength at the absorption site has a lower limit of tens of Gauss, which is necessary for the phenomenon to occur. This characteristic can be explained in a scenario where GPM 1839-10 possesses a weakly magnetized companion star.

The presence of dust in supernova remnants (SNRs) is confirmed by extensive infrared data from observatories such as Spitzer, Herschel, and JWST, alongside theoretical models of dust formation. This study explores the existence of dust and ices, particularly water ice via 62 {\mu}m in SNRs such as the Crab Nebula and N49, using observational data and preliminary modeling with Cloudy. Observations suggest that water ice may be present in IC 443 and possibly other remnants, though the 63 {\mu}m band could also indicate [OI] emission. Theoretical models indicate that water ice could survive under certain conditions in SNRs, with densities and temperatures analyzed. Further observations and refined simulations are needed to confirm these findings.

Water ice has been detected in several supernova remnants despite the strong heating and radiation in these environments. This challenges standard expectations for dust survival. Using two dimensional magnetohydrodynamic simulations, we study how a supernova shock interacts with a dense interstellar cloud. The simulations show that the shock naturally compresses the cloud into dense structures similar to those inferred in well known remnants. Although temperatures remain high in the adiabatic phase, simple considerations indicate that cooling would act quickly once included. These results suggest that shock cloud interactions create the physical conditions needed for water ice to form. Future work including radiative cooling and grain surface chemistry will allow direct modelling of ice growth in these compressed regions.

The Milky Way hosts astrophysical accelerators capable of producing high-energy cosmic rays. These cosmic rays can interact with the interstellar medium (ISM) across the Galaxy to produce neutrinos and gamma rays (propagation component), while their interactions with ambient material at their acceleration sites, such as supernova remnants, can give rise to the source component of the gamma-ray and neutrino flux. In this paper, we estimate the source component of the Galactic neutrino flux using simulated populations of Galactic gamma-ray sources. We compare our results with observations from neutrino experiments in the energy range of 1-30 TeV. Using simulated populations of Galactic TeV gamma-ray sources, we exploit the correlation between gamma rays and neutrinos and introduce a bracketing approach to constrain the range for the source contribution of the Galactic neutrino flux. For the upper limit, we used a simulation describing the entity of Galactic gamma-ray sources, whereas the lower limit was estimated using the hadronic component of the Galactic supernova remnant population. Our results show that the difference between this maximum and minimum is less than an order of magnitude and the flux range is comparable to the Galactic neutrino flux from the cosmic-ray interaction with the ISM. The results agree with the observed signals from IceCube and ANTARES and suggest that the propagation component, combined with the minimum source contribution predicted by the supernova-remnant model, approaches the observed neutrino flux, leaving little room for significant enhancements of the emission originating from propagating cosmic rays.

Massive stars dominate the radiative and mechanical feedback of young stellar populations, yet their intense ultraviolet fields and strong winds are typically presumed to preclude Earth-like habitability. We quantify this expectation by mapping time dependent habitable zones (HZs) for solar-metallicity stars with initial masses of $0.8$-$120\,M_\odot$. From rotating and non-rotating \textsc{GENEC} tracks we derive bolometric ``climate'' HZ boundaries and enforce XUV energy-limited escape and wind ram-pressure retention constraints for a dipole-magnetized Earth analogue. The operational inner edge is set by the most restrictive limit, and we measure the annulus lifetime, the longest continuous residence at fixed orbit, and the maximum number of dynamically packed terrestrial planets it can host. We find a sharp main-sequence ceiling: while a $9\,M_\odot$ star sustains an operational HZ for $\sim 30$~Myr at $\sim 70$-$130$~AU, the main-sequence annulus becomes brief and extremely narrow by $12\,M_\odot$ and disappears by $15\,M_\odot$. Post main-sequence evolution can reopen HZs up to $\sim 25$-$30\,M_\odot$, but only for $\sim 0.03$-$1.5$~Myr at hundreds to $\sim 10^3$~AU, disappearing by $\sim 40\,M_\odot$. Rotation modestly increases habitable lifetimes near the upper main sequence without altering the high mass ceiling. Initial Mass Function (IMF) weighting shows that massive stars contribute only $\sim 10^{-4}$ of the habitable planet-time budget. Even so, they still add of order a few $10^{5}$ operationally habitable Earth analogues to the Milky Way at any instant. This implies that massive star systems are unlikely to dominate the Galaxy wide habitability budget, but they may still provide a set of short-lived, observationally distinct targets for biosignature searches.

Modern studies in astrophysics and cosmology increasingly rely on simulations and cross-survey analyses, yet differences in data generation, instrumentation, calibration, and unmodeled physics introduce distribution mismatches between datasets (``domain shift''). In machine-learning pipelines, this occurs when the joint distribution of inputs and labels differs between the training (source) and application (target) domains, causing source-trained models to underperform on the target. Transfer learning and domain adaptation provide principled ways to mitigate this effect. We study a concrete simulation-to-observation case: semi-supervised domain adaptation (SSDA) to transfer a four-class spectral classifier -- high-redshift quasars, low-redshift quasars, galaxies, and stars -- from J-PAS mock catalogs based on DESI spectra to real J-PAS observations. Our pipeline pretrains on abundant labeled DESI$\rightarrow$J-PAS mocks and adapts to the target domain using a small labeled J-PAS subset. We benchmark SSDA against two baselines: a J-PAS--only supervised model trained with the same target-label budget, and a mocks-only model evaluated on held-out J-PAS data. On this held-out J-PAS data, SSDA achieves a macro-F1 score (balancing precision and recall) of $0.82$ and an overall true positive rate of $0.89$, compared to $0.79/0.85$ for the J-PAS--only baseline and $0.73/0.87$ for the mocks-only model. The gains are driven primarily by improved quasar classification, especially in the high-redshift subclass ($\mathrm{F1}=0.66$ vs.\ $0.55/0.37$), yielding better-calibrated candidate lists for spectroscopic targeting (e.g., WEAVE-QSO) and AGN searches. This study shows how modest target supervision enables robust, data-efficient simulation-to-observation transfer when simulations are plentiful but target labels are scarce.

We present the absolute flux and wavelength calibration of the James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam) Wide Field Slitless Spectroscopy (WFSS) mode. Each of NIRCam's two modules (A and B) provides independent long wavelength (LW) grism spectroscopy over the 2.4-5.0 micron range, with orthogonally oriented R and C grisms. Using commissioning and calibration data from programs 01076, 01536, 01537, 01538, 01479, 01480, 04449, 04498, 06606, and 06628, we have measured the field-dependent geometry and wavelength dispersion of both first and second order spectra across the full detector area. The trace geometry was modeled using two-dimensional third-order polynomials that reproduce the observed spectral positions with an RMS accuracy better than 0.1 pixel. Wavelength calibration, derived from observations of the planetary nebula SMP LMC 58, achieves a precision of 1.6-1.9A for the +1 orders and 0.5A for the +2 orders. Absolute flux calibration, established from observations of the G-type star standard P330E, provides a consistent sensitivity function across all grisms and modules with an absolute flux accuracy of 3\%. The resulting calibration framework defines the geometric, wavelength, and photometric reference for all NIRCam WFSS observations and ensures cross-consistency between modules and grism orientations. These calibrations form the basis for accurate slitless spectroscopy with NIRCam and will support ongoing improvements to the JWST calibration pipeline and data products.

We report a newly identified nonlinearity in H4RG-10 near-infrared detectors operating under moderately elevated-temperature conditions (114 K). This component, that potentially arises from illumination-independent defect currents, introduces additional nonlinearity not captured by conventional correction models. To address this issue, we propose a data-driven nonlinearity correction (NLC) method that models the nonlinear behavior of both the classical response and the defect currents, using dual-illumination measurements and a dark exposure. Applied to H4RG-10 detectors on the PRIME telescope, the method significantly improves signal linearity, especially for pixels with high defect current, while maintaining comparable performance elsewhere. By selecting the optimal correction model per pixel, reliable NLC is achieved across the full array. This study characterizes a nonlinearity intrinsic to H4RG-10 detectors and demonstrates that data-driven post-processing can effectively restore linearity in the presence of large defect currents. Although these effects are unlikely to be significant under nominal operating temperatures, the approach may provide a practical calibration framework for future warm-operation scenarios.

An analysis of the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has produced a sample of 27 candidate dense molecular clumps with large collapse motions, as revealed by large ``blue'' asymmetrical line profiles of the optically thick \hcop\, line. %with respect to the optically thin \nthp\, line. New, more sensitive molecular line observations of this sample, conducted with the Mopra 22-m telescope, confirm the blue asymmetries in the \hcop\, line profiles, with large, positive values of the asymmetry parameter $A$ ($\bar{A}_{HCO^+} = 0.69\pm0.01$), and positive, but smaller asymmetries in the \hcn\, and \hnc\, lines: ($\bar{A}_{HCN} = 0.35\pm0.01$ and $\bar{A}_{HNC} = 0.28\pm0.01$), as expected for a less optically thick tracer in collapsing clumps. The small, positive mean asymmetry parameters for \cch\, and \htcop, $\bar{A}_{C_2H} = 0.15\pm0.02$ and $\bar{A}_{H^{13}CO^+} = 0.18\pm0.03$, likely indicate slightly optically thick emission for at least some clumps. The hyperfine ratios for \nthp\, are in their optically thin, LTE, values, but for \hcn\ they are not; the $F=1 \to 1$ hyperfine line shows abnormally weak intensities. A simple two-component model shows that self-absorption of the background $F = 1 \to 1$ hyperfine line by the main $F = 2 \to 1$ hyperfine line of a cold, foreground, redshifted cloud can reproduce the observed \hcn\, hyperfine intensities and match the \hcn\, and \hcop\, line profiles. All of these results are consistent with self-absorption of the optically thick lines on the red side of the profile, as expected for collapsing clumps. A simple two-cloud model suggests that this sample represents dense clumps with extreme collapse velocities, $V_{inf} \sim 2.4$ \kms.

Since the discovery of superhumps in 1974, these photometric modulations have provided a crucial observational window into disk instabilities in cataclysmic variable stars, particularly the tidal instability associated with the 3:1 resonance. Over the past few decades, extensive time-resolved photometry has revealed a rich diversity of superhump-related phenomena, including delayed superhump development, early superhumps in WZ Sge-type dwarf novae, systematic stage A-B-C evolution, negative superhumps, and superhumps observed in related systems such as intermediate polars and AM CVn stars. In this invited review, we summarize key observational advances since the establishment of the thermal-tidal instability framework, discuss their theoretical interpretations within the disk instability model, and highlight remaining open problems. These developments have been driven by coordinated networks of amateur observers, wide-field robotic surveys, and continuous high-precision space-based photometry from Kepler and TESS. Together, they demonstrate that superhumps remain a powerful probe of disk dynamics, binary parameters, and the interplay between thermal, tidal, and geometric effects in accretion disks.

We review the algebraic definition of the efficiency of a polarization modulation scheme, which is commonly adopted for solar and stellar spectro-polarimetry applications, and generalize it to allow distinct states of the modulation cycle to have arbitrary throughput and different photon-noise statistics for each state. Such a generalization becomes necessary to model and optimize the polarimetric efficiency of instruments implementing spatial polarization modulation schemes, where different optical paths are assigned to different polarization analysis states, which may be characterized by different throughput values. The proposed algebraic extension also proves essential for introducing a workable concept of the efficiency of a polarization calibration scheme, which can then be used to create a merit function for the optimization of calibration sequences, which take into account the specific characteristics of the polarimetric instrument and of its calibration optics.

We expand on our previous study of the impact of atmospheric seeing on polarization cross-talk, and show how the formalism that was developed in that work can be applied to treat the case of spatial modulators of polarization. Beside formally demonstrating how the problem of cross-talk is fully eliminated in such devices, we also gain insight on the meaning of polarimetric noise of temporal modulation schemes in the limit of very high modulation frequency. We also describe the problem of spectrograph instabilities, and how the spectral gradients that are naturally associated with a line spectrum feed into the problem of polarimetric errors induced by mechanical vibrations, thermal drifts, and pointing jitter. Finally, we show how this formalism can be used to estimate the contribution of polarization cross-talk to the errors on the elements of the 4$\times$4 Stokes response matrix, for the purpose of producing realistic error budgets for polarimetric instrumentation.

The past decade has been a transformative period for pulsar timing arrays (PTAs) and their search for nanohertz gravitational waves (GWs). This progress has been driven by collective advances in instrumentation for pulsar timing observations, increasingly sophisticated data-analysis techniques, and improved theoretical understanding of the origins of nanohertz GW signals. PTA sensitivity has steadily improved, leading first to progressively more stringent upper limits on the gravitational-wave background (GWB), and subsequently to the identification of a common red-noise process in pulsar timing data, the first hint of a GWB. In 2023, multiple PTA collaborations reported evidence for the Hellings-Downs correlation, widely regarded as the definitive signature of a GWB. These developments place PTAs on the threshold of a confident GW detection and the opening of a new low-frequency window on the GW Universe. In this article, we present an overview of PTA experiments, with particular emphasis on the rapid progress achieved during this pivotal period for PTA and nanohertz GW science.

Lorentz symmetry is a fundamental pillar of modern Physics, yet high-energy theories often predict its violation. One potential signature of such a violation is cosmic birefringence - rotation of the polarization plane of photons due to Chern-Simons coupling in Maxwell's electrodynamics. This rotation angle, aka birefringence angle, depends upon the distance travelled by the photon and is thus different for CMB and 21cm photons. While the rotation angle in CMB, i.e., $\alpha_\mathrm{CMB}$, has been tightly constrained by CMB experiments, the potential of the 21cm cosmological signal to constrain this parameter, as well as constrain $\alpha_\mathrm{21cm}$, remains largely unexplored. In this work, we provide constraints on both these angles by cross-correlating 21cm and CMB signals. Using the Fisher matrix formalism, we give our forecasts for 21cm experiments, including SKA, HIRAX, and PUMA, and Planck like CMB experiment. We find that best constraints $\sigma_{\alpha_\mathrm{CMB}} \sim 4.4^\circ$ and $\sigma_{\alpha_\mathrm{21cm}} \sim 100^\circ$ are found using $C_\ell^{T_{21} B_\mathrm{CMB}}$ and $C_\ell^{T_{21} B_{21}}$ respectively. Since birefringence hasn't yet been detected in 21cm, we choose the fiducial value $\alpha_\mathrm{21cm}^\mathrm{fid}=0$ assuming the null hypothesis.

We present a parametric strong-lensing model for the galaxy cluster MACS J1931.8-2635 ($z_l = 0.35$), accompanying the detection of the spectroscopically confirmed SN Eos at $z = 5.13$ (Coulter et al. 2026). We identify 10 new multiple-image systems in recent VENUS JWST/NIRCam imaging, so that the model is constrained with a total of 19 robust multiple-image systems -- nine of which also have a spectroscopic redshift. For the point-like source corresponding to SN Eos, our model predicts a total of five images, with the observed radial image pair having a similar magnification of $\mu \simeq 25 - 30$ and a small time delay of $< 5$ days, in agreement with their simultaneous observation. According to the model, the other three predicted images arrived earlier, with time delays of $3.7 \pm 0.7$, $3.5 \pm 0.7$ and $54.0 \pm 10.8$ years prior to the two observed images, and with magnifications of $12.9 \pm 2.6$, $13.0 \pm 2.9$ and $2.2 \pm 0.4$, respectively. The absence of detections at the predicted positions, where the host galaxy's images are also visible, confirms the transient nature of the source. SN Eos and its host galaxy are studied in separate articles, and we here focus on the lens model. The final model reaches a very good $r.m.s.$ distance between model and observations of $0.44''$. We present the lens-modeling results, including newly identified systems such as a triply imaged, grand-design spiral galaxy candidate at $z \simeq 3.65_{-0.09}^{+0.04}$, and briefly discuss the potential of using high-redshift lensed SNe for cosmography.

We combine narrow/medium-band filter photometry from the Southern Photometric Local Universe Survey (S-PLUS) DR4 with ultra broad-band filter photometry from Gaia EDR3 to derive fundamental stellar parameters ($T_{\rm eff}$, $\log g$, [Fe/H], ages) and elemental-abundance ratios ([C/Fe] and [$\alpha$/Fe]) for 5.4 million stars in the Galaxy (4.9 million dwarfs and 0.5 million giants), as well as for over 0.7 million red giant stars in the Large and Small Magellanic Clouds (LMC and SMC). The precisions of the abundance estimates range from 0.05-0.10 dex for metallicity in the relatively metal-rich range ([Fe/H] $> -1.0$) to 0.10-0.30 dex in the metal-poor regime ([Fe/H] $<-1.0$), 0.10-0.20\,dex for [C/Fe], and 0.05 dex for [$\alpha$/Fe]. The stellar parameters for LMC and SMC member stars are somewhat less precise than those from the S-PLUS main survey, primarily because of the effect of high reddening. The use of both metallicity- and carbon-sensitive filters provides unbiased measurements of both [Fe/H] and [C/Fe], of particular importance for very low-metallicity ([Fe/H] $< -2.0$) stars, where carbon enhancement can lead to systematically high estimates of [Fe/H] when only a single metallicity-sensitive filter is employed. Furthermore, multiple narrow-band filters enable metallicity estimates down to [Fe/H] $\sim -4.0$ with an accuracy of around 0.3 dex, exceeding the precision typically achieved by low/medium-resolution spectroscopy. This extensive photometric dataset, combined with the other three datasets in this series, will serve as a valuable legacy resource for Milky Way and Magellanic Clouds studies.

Magnetic properties of dust grains play a fundamental role in their alignment with ambient magnetic fields and magnetic dipole emission. In the radiative torque (RAT) paradigm, superparamagnetic inclusions (SPIs) embedded within dust grains are expected to significantly enhance magnetic susceptibility and alignment efficiency. Previous studies have generally assumed SPIs of a single characteristic size. In this work, we develop an effective superparamagnetism model that explicitly accounts for a power-law size distribution of SPIs. We show that the effective zero-frequency susceptibility can be described by the superparamagnetic susceptibility of uniform-sized inclusions evaluated at the critical blocking size, reduced by a factor $F_{\rm eff}\sim 0.1$. It exhibits a slight increase with dust temperature $T_{d}$, in contrast to the rapid decrease for the case of single-size SPIs. For rotating grains at angular frequency $\omega$, we identify a characteristic resonance size of SPIs that dominates the magnetic response, $N_{\rm res} = (T_{d}/T_{\rm act}) \ln (\nu_{0}/\omega)$ with $T_{\rm act}$ activation temperature and $\nu_{0}$ the characteristic attempt frequency of SPIs. The frequency-dependent effective susceptibility is well described by the maximum susceptibility $\chi_{\rm eff}^{\rm max}(\omega)$ at $N_{\rm res}$, reduced by a factor $G_{\rm eff}\sim 0.1$. Unlike models assuming uniform-sized inclusions, we find that the effective susceptibility exhibits a nearly flat spectrum for frequency below $\nu_{0}$, arising from the progressive activation of larger inclusions at lower frequencies. This effective superparamagnetism model based on the SPI size distrbution has important implications for magnetic grain alignment, dust polarization, and magnetic dipole emission across diverse environments.

We investigate the potential birthplace of Galactic binary neutron star (BNS) systems through a kinematic analysis. Using high-precision astrometry from Gaia DR3, updated pulsar distances, and Monte Carlo sampling of astrometric errors, we integrate the past trajectories of 11 Galactic BNSs and 167 globular clusters plus 2967 open clusters, to search for past encounters. Our results suggest that BNS origin in globular clusters is unlikely, with low encounter probabilities (e.g., $\lesssim 0.5\%$ for NGC 5139) and requiring excessive ejection velocities. Conversely, our analysis indicates that open clusters are a non-negligible formation channel. Specifically, the double pulsar J0737$-$3039 shows a $13.9\%$ ($5.4\%$) probability of originating from the young cluster OC 0450 (Theia 58). Based on encounter proximity and time, we argue that Theia 58 is its more plausible birthplace. Our work provides kinematic evidence consistent with an open-cluster origin for a subset of field BNSs.

The center of the Galaxy is a prominent source in X-rays and gamma-rays. The study of its high-energy (HE) emission is crucial in understanding the physical phenomena taking place in this dense and extreme environment, where the closest supermassive black hole (SMBH) to us, Sgr A*, is lurking nearly invisible, today, in most of the energy spectrum. These phenomena are probably common to other galactic nuclei and may explain the feedback processes between nuclear regions and galaxies, so important for the overall evolution of the Universe. The Galactic center HE emission is very complex and consists of both thermal and non thermal radiation produced by compact and extended sources, surrounded by more diffuse components. All these objects and media are interacting with each other in the narrow and dense Central Molecular Zone (CMZ). Some of them also show relevant extensions towards the Galactic poles, indicating energetic outflows that seem to link the center to the recently observed large Galactic polar structures. In spite of the fundamental advances obtained in the last twenty five years with the most sensitive X-ray and gamma-ray observatories, several questions remain open to investigations. We review here the main observational results and the open issues on the high-energy diagnostics of the Galactic nuclear activity, focusing on processes that take place in the CMZ, and in particular discussing the role of the present and past SMBH activities in powering this region and possibly the whole Galaxy.

This study demonstrated the pulsar navigation capability of the CubeSat X-ray observatory NinjaSat, which is equipped with two Gas Multiplier Counters (GMCs). The GMCs are sensitive to the 2-50 keV energy band and have an effective area of 16 cm^2 per module at 6 keV. We verified the timing accuracy by observing the Crab Pulsar and confirmed stable timing performance within 100 microseconds. To demonstrate pulsar navigation, we applied a method that optimizes orbital parameters to maximize the significance of the pulsar X-ray pulse profile, known as the Significance Enhancement of Pulse-profile with Orbit-dynamics (SEPO) method. We observed the Crab Pulsar with a total exposure of approximately 100 ks at different epochs and analyzed the data transmitted to the ground. By comparing the optimized orbit with the satellite position derived from Global Positioning System data, we quantitatively evaluated the navigation performance. The results show that the position component along the Crab line of sight was consistently constrained within approximately 40 km, and the three-dimensional position error ranged from 27 to 370 km depending on the observation epoch. These results demonstrate the feasibility of applying a CubeSat-class X-ray observatory to pulsar navigation and provide the first experimental verification that the accuracy of the SEPO method depends on the seasonal geometry between the orbital plane and the pulsar direction.

The discovered planets in apparent single-planet systems (singles) and those in systems with multiple detected planets (multis) exhibit a rich diversity of physical and orbital properties. We investigate the differences and similarities between 1730 singles and 1522 multis in a catalogue of confirmed transiting planets orbiting main-sequence stars with spectral classes ranging from late-M to late-F. After we removed the hot Jupiters, the planet types and their fractional numbers were similar for the multis and singles hosted by FGK-type stars. Furthermore, the median radii of both the singles and the multis increase with host star temperature already from late- to early-type M dwarfs and further up to F stars. Our analyses show that the singles are larger on average than the multis in our F, G, K, and M samples. However, after we excluded the hot Jupiters, the radius distributions of the singles and multis orbiting FGK stars are statistically indistinguishable, particularly at $R < 4\, R_\oplus$. In the FGK sample, we also identified an unexpected and significant overabundance of multis, compared to singles, at radii of $\approx1.4-1.6\, R_\oplus$. For the early- and late-type M samples, our work indicates that the multis are smaller on average than the singles and that the radius distributions of the multis and singles are different, except for the planets with $R < 4\, R_\oplus$ hosted by early-type M dwarfs. Nevertheless, these results for the two M samples are inconclusive because the sample sizes of 167 and 101 planets are limited. In conclusion, our analyses reveal that the singles and multis, excluding the hot Jupiters, orbiting FGK stars are overall consistent with originating from the same underlying population based on their planet types and radii. This interpretation is, however, not applicable within the region of the overabundance of multis identified in this work.

We develop a self-consistent nonlinear extension of diffusive shock acceleration that incorporates cosmic ray (CR) backreaction on the shock precursor together with a physically motivated upstream-escape mechanism that produces an exponential high energy cutoff. The CR pressure gradient decelerates the upstream flow facing the shock wave, generating an extended precursor in which higher rigidity particles sample a larger cumulative velocity gradient and thereby acquire a progressively harder spectrum. Finite-size/escape effects are modeled by a momentum-dependent loss term, which naturally terminates acceleration and steepens the spectrum near the cutoff. The precursor compression ratio is not imposed as a closure condition but is determined dynamically by requiring consistency between the injection rate inferred from thermal leakage at the subshock and the injection strength demanded by the nonlinear shock modification, with CR-driven wave heating providing stabilizing negative feedback. Applying the model to young supernova-remnant-like parameters and standard one-zone Galactic diffusion, we reproduce the main features of the latest DAMPE proton spectrum: gradual hardening from hundreds of GeV to multi-TeV energies and a subsequent exponential cutoff at tens of TeV. The resulting spectral evolution follows directly from the competition between precursor-mediated nonlinear feedback and upstream escape.

We report the Fe K relativistic reverberation feature for the first time in the Narrow-line Seyfert\,1 galaxy Ton\,S180. Using a long observation from {\it XMM-Newton} we find that the Fe K emission lag peaks at $117\pm49$ s in the lag energy spectrum computed for frequencies $(0.3-8.5) \times 10^{-4}$ Hz. The lag amplitude drops to $22.85\pm14.20$ s as the frequency increases to $(8.5-30) \times 10^{-4}$ Hz. The time-averaged spectrum of the source shows a relatively narrow Fe K line at $\sim6.4$ keV, resulting in black hole spin to be low ($a=\rm 0.43_{-0.14}^{+0.10}$) found from the reflection modelling. We perform general relativistic transfer function modelling of the lag energy spectra individually. This provides an independent timing-based measure of the spin at $a=0.30_{-0.17}^{+0.34}$, and black hole mass $M_{\rm BH} = 0.29_{-0.16}^{+0.01}\times10^8M_{\odot}$, comparable to the previous measurement, and height of the corona $h = 2.59_{-0.33}^{+5.17}r_{\rm g}$. Further, we observe that the Fe K lag and the black hole mass fit well in the linear lag-mass relation shown by other Seyfert 1 galaxies.

The SS433, a well-known binary system with an internal black hole, have shown since half a century, an inner (a few year light distances) twin precessing jets spirals. These beams are made by tidal forces while stripping mass from large stellar companion feeding an inner BH accretion disk and an orthogonal accelerating twin jet. From it, the radio, X gamma jet emission. A couple of years ago H.E.S.S telescope as well as HAWC and LHAASO array detectors, discovered also the surprising signature of an unexpected far twin separated gamma beam at tens TeV energy. At a hundred light years distances from its central source. We suggest that it is the legacy of a past rare eruption, a century ago, of tens PeV (10^16 eV) relativistic twin neutron beams. Their beta decay in flight at far distances, into proton, neutrino and in particular into tens TeV electrons, could feed the observed TeV gamma traces. They are originated by the same secondary tens TeV electrons emitting hard gamma, by Inverse Compton Scattering onto interstellar infrared photons.

Interstellar comets provide rare opportunities for probing the diversity of refractory and volatile inventory around other stars. As the second ever interstellar comet, and the third interstellar object, 3I/ATLAS has been the focus of telescopic observations since its discovery in July 2025. Following the previous observations at multi-wavelengths, we present further radio observations of the 1665/1667 MHz ground-state OH lines and millimeter observations of the CO($J$=1-0) transition at 115.271 GHz that trace the coma $\rm H_2O$ and CO abundances, respectively. We derived OH production rates of $(1.32\pm0.47)\times10^{28}\ \rm s^{-1}$ at 2.27 au and $(1.89\pm0.37)\times10^{28}\ \rm s^{-1}$ at 1.96 au as well as an average CO production rate of $\rm (5.75\pm1.91) \times 10^{27}\, s^{-1}$ between 2.33 and 1.75 au, inferring a CO/$\rm H_2O$ ratio of ($28\pm11\%$). With the mean HCN production rate of $2.5\times 10^{25}\ \rm s^{-1}$ at 2.1 au reported by \citet{2025arXiv251120845R} and \citet{2025arXiv251002817C}, we infer a CO/HCN ratio of ($230\pm76$). By synthesizing water production rates measured with instruments of different apertures, we found that the sublimation from extended source in the coma contributes significantly to 3I's pre-perihelion water measurements, accounting for up to 80\% from 3 au to 2 au.

In the I-band color-magnitude diagrams (CMD) of resolved nearby galaxies, the reddest asymptotic giant branch (AGB) stars form a previously unremarked-upon, but nevertheless distinct and easily-identified population of high-luminosity stars. Hereafter we refer to this population as being comprised of I-Band AGB (IAGB) stars. Identifying these stars in the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC) and in NGC4258 (for all three of which there are published geometric distances) we find that the marginalized luminosity functions are each well approximated by single-peaked Gaussians, having one-sigma dispersions of +/- 0.22 mag, +/- 0.25 mag and +/- 0.24 mag, respectively. The zero points for the modal I-band absolute magnitudes of IAGB stars are found to be M_I = -4.49 +/- 0.003 mag (stat) in the LMC (4204 stars), M_I = -4.67 +/- 0.008 mag (stat), for the SMC sample (916 stars), and M_I = -4.78 +/- 0.030 mag (stat) for NGC4258 (62 stars). A global average over these three independent calibrations of the IAGB zero point (weighted inversely by squares of their systematic errors) gives <M_I> = -4.65 +/- 0.119 mag (stat) +/- 0.025 (sys). In Paper II we will show the results of applying the IAGB Method to 92 galaxies additional galaxies resolved by HST, reaching out to distances just short of 10 Mpc.

The advent of new and more powerful infrared spectrographs has significantly motivated the advancement of the study of atomic and molecular line lists and stellar atmosphere models. While optical abundance determinations rely on extensively validated line lists and modeling frameworks, infrared measurements still face larger uncertainties, largely driven by the choice of atmospheric models and the quality of the available atomic data. In this work, we aim to deliver a homogeneous and reproducible set of atomic absorption lines in the Y, J, and H bands (9800 - 18000 (Angstrom)), based exclusively on laboratory atomic data. We analyse CRIRES spectra of six Gaia FGK Benchmark Stars spanning a wide range in effective temperature, surface gravity, and chemical composition. Synthetic spectra are computed using the benchmark stellar parameters, and each transition is evaluated independently in every star through a quantitative sequence that examines line depth, saturation, blending (purity), and the agreement between observed and synthetic line profiles. We identify a set of robust atomic transitions in these bands that remain consistent across the full range of stellar parameters represented in our sample. Lines of alpha-elements such as Mg I, Si I, and Ca I, together with several Fe I transitions, satisfy all robustness criteria. Among the neutron-capture species explored, only Sr II provides lines that consistently meet our requirements. Beyond the specific list of accepted transitions, this study demonstrates that a fully quantitative, multi-criteria framework provides a transparent and reproducible foundation for near-infrared line validation as laboratory data, stellar atmosphere models, and instrumentation continue to improve.

We demonstrate that magnetospheric interactions between merging neutron stars (NSs) generate dual-jetted current outflows, analogous to the Alfvén wings observed during planetary interactions in the Solar System. Using 3D relativistic MHD simulations, we model the interaction as a conducting sphere moving through a highly magnetized plasma of the companion's magnetosphere. Unusually, the interaction operates in a regime that is relativistic yet sub-Alfvénic. Electromagnetic draping amplifies magnetic fields in a narrow layer near the stellar surface, leading to the generation of electric currents along the local magnetic field. The generation of beamed outflows enhances instantaneous power of the pulsar-like radio and high energy emission, produces spin/orbital modulations, and is likely to lead to observable precursor emission preceding the main gravitational wave event.

Hubble Space Telescope (HST) observations of 92 galaxies that have a strong showing of I-band Asymptotic Giant Branch (IAGB) stars in their color-magnitude diagrams (CMDs), are used to measure the relative offset between the mean apparent I-band magnitudes of the IAGB population and the corresponding apparent I-band magnitudes of the TRGB as measured in the same frames (and CMDs) of those individual galaxies. This first exploratory, large-sample comparison is independent of any extinction (foreground or internal) that may be shared by these two populations. The marginalized luminosity functions used to determine the modal value of the {\it IAGB } population are well fit by a single, symmetric Gaussian. The difference in the two apparent magnitudes (in the sense IAGB minus TRGB) is -0.589 mag, with a combined standard deviation of +/- 0.119 mag. Adopting M_I = -4.05 mag for the TRGB stars, the modal absolute magnitude of the IAGB is then calculated to be M_I(IAGB) = -4.64 +/- 0.12 mag. The ensemble dispersion quoted above gives a standard error on the mean of +/- 0.012 mag (based on the full sample of 92 galaxies). Independently, the three geometry-based zero points for I-band AGB stars are found (in Paper I) to be M_I = -4.49 +/- 0.003~mag in the LMC (4204 stars), M_I = -4.67 +/- 0.008 mag, for the SMC (916 stars) and M_I = -4.78 +/- 0.030 mag for NGC4258 (62 stars), leading to a global zero-point (weighted) average of <M_I> = -4.64 +/- 0.15 mag (stat). The scatter found in the anchors is comparable to the scatter in the field sample discussed here, but the calibration sample is small. The application of this method to galaxies well outside of the Local Group, shows that these standard candles can readily be found and measured out to at least 9 Mpc, using already available archival data

The recent detection of a gravitational wave background in the nano-Hertz frequency range by Pulsar Timing Array (PTA) collaborations, including NANOGrav, EPTA, and PPTA, has opened a new avenue for exploring fundamental physics in the early universe. In this work, we analyze a supercooled first-order phase transition in a hidden sector with a spontaneously broken $U(1)_X$ gauge symmetry as a source for this signal. We demonstrate that the thermal history of the hidden and visible sectors plays a crucial role in the gravitational wave power spectrum analysis. Our analysis shows that supercooled phase transitions can generate gravitational waves strong enough to explain the PTA observations while satisfying cosmological constraints from Big Bang Nucleosynthesis.

Europa, Ganymede, and Callisto are key targets in the search for habitability due to the potential presence of subsurface oceans. Detecting complex organic molecules (COMs), essential for prebiotic chemistry, is crucial to assessing their potential. Though COMs remain undetected on these moons, ESA's JUICE and NASA's Europa Clipper missions aim to fill this gap with their science payloads. This study explores the formation and transport of COMs within Jupiter's circumplanetary disk (CPD), a critical environment for the formation of the Galilean moons. Using a time-dependent model that couples the evolving CPD structure with the dynamics of icy particles of varying sizes and release times, we assess two primary COM formation pathways: thermal processing of ices and UV photochemistry. The results indicate that heating, particularly of NH3:CO2 ices, occurs efficiently before substantial irradiation, making it the dominant pathway for COM formation in the Jovian CPD. However, the relative efficiencies of these two processes are governed by particle density, disk viscosity, accretion rate, and UV flux, which collectively determine drift timescales and exposure to favorable thermodynamic environments. Existing models indicate that Europa's accretion was relatively cold and prolonged, possibly allowing some COMs to survive incorporation, whereas Ganymede and Callisto likely formed under even cooler conditions conducive to preserving COM-rich material. These results highlight the potential inheritance of complex organics by the Galilean moons and offer a framework for interpreting upcoming compositional data from JUICE and Europa Clipper.

We construct a concept-object knowledge graph from the full astro-ph corpus through July 2025. Using an automated pipeline, we extract named astrophysical objects from OCR-processed papers, resolve them to SIMBAD identifiers, and link them to scientific concepts annotated in the source corpus. We then test whether historical graph structure can forecast new concept-object associations before they appear in print. Because the concepts are derived from clustering and therefore overlap semantically, we apply an inference-time concept-similarity smoothing step uniformly to all methods. Across four temporal cutoffs on a physically meaningful subset of concepts, an implicit-feedback matrix factorization model (alternating least squares, ALS) with smoothing outperforms the strongest neighborhood baseline (KNN using text-embedding concept similarity) by 16.8% on NDCG@100 (0.144 vs 0.123) and 19.8% on Recall@100 (0.175 vs 0.146), and exceeds the best recency heuristic by 96% and 88%, respectively. These results indicate that historical literature encodes predictive structure not captured by global heuristics or local neighborhood voting, suggesting a path toward tools that could help triage follow-up targets for scarce telescope time.

In this study, we aim to constrain the atmospheric composition of GJ 1214 b using all available transits observed with the upgraded CRIRES+ spectrograph at the VLT by searching for the signatures of water vapour, methane, and carbon dioxide. We analysed eight CRIRES+ transit datasets covering the K band (1.90-2.45 microns) at a resolving power of R ~ 100,000. We used the SysRem algorithm to correct for telluric and stellar contributions and employed the cross-correlation technique with templates from petitRADTRANS to search for H2O, CH4, and CO2. Injection-recovery tests across a grid of metallicities (Z) and cloud-deck pressures (pc) were performed to quantify detection limits. We also generated predictions for ANDES observations using end-to-end simulated datasets with EXoPLORE. We detect no significant H2O, CH4, or CO2 signatures. Injection-recovery tests show that such non-detections exclude atmospheres with low-altitude clouds and moderate or low metallicities. CH4 yields the tightest empirical limits, with CO2 unexpectedly ruling out intermediate metallicities (~ 100xsolar) with clouds deeper due to its rapidly rising opacity in compressed, high-Z atmospheres. Our constraints are in line with either a high-Z or a high-altitude aerosol layer, in agreement with recent JWST inferences. The combined analysis of eight CRIRES+ datasets provides the most stringent high-resolution constraints on the atmospheric properties of GJ 1214 b to date. Simulations of a single transit observed with ANDES on the ELT predict modest improvements for H2O, a substantially expanded detectable region for CH4, and the strongest gains for CO2, making the latter a particularly effective tracer for characterising high-metallicity atmospheres in sub-Neptunes.

By identifying quasi-periodic oscillations (QPOs) reported in FRB 20240114A (from the Five-hundred-meter Aperture Spherical Telescope) with neutron star crustal torsional oscillations, together with experimental constraints on the incompressibility $K_0$ of symmetric nuclear matter at saturation density, we constrain the mass and radius of an extragalactic neutron star at redshift $z\approx0.13$. Identifying the low-order QPO frequencies as fundamental oscillations, and frequencies of $567.7\,\mathrm{Hz}$ or $655.5\,\mathrm{Hz}$ (rest frame) as first overtone candidates, implies neutron star mass ranges of $1.00$--$1.55\,M_\odot$ or $1.17$--$1.76\,M_\odot$, respectively. The radius is also constrained, with a self-consistent value around $13$~km, further supported by the calculation of the NS structure within the low-mass/low-central density regime. Simultaneously, we also constrain another nuclear saturation parameter, namely the density dependence of the nuclear symmetry energy at saturation density (i.e., the slope parameter), $L$, and determine it to be $L=59.5-96.8$ MeV with $\sim 10\%$ systematic uncertainty, which is broadly consistent with previous constraints on $L$ obtained from experiments and astronomical observations. Thus, a mapping of FRB QPOs to crustal torsional modes seems reasonable. This can be confirmed with upcoming FRB surveys over a broad range of redshifts and more elaborate data analyses.

Combinations of recent cosmological observations, including Dark Energy Spectroscopic Instrument (DESI), show hints of a dynamical nature for dark energy. While the data suggest the possibility of the phantom crossing, it is worth thoroughly exploring quintessence models. Given that phenomenological parametrisations of the equation-of-state parameter $w(a)$ with a sharp transitional feature fit the data well, we study the realisation of such models in quintessence. In the late Universe, the quintessence field begins to oscillate abruptly, changing the behaviour of $w$. Naturally, such a model entails tachyonic instability, and particle production modifies $w$. We perform numerical lattice simulations to study the time dependence of $w$. In addition, the violent particle production produces significant density perturbations and the stochastic gravitational-wave background, whose characteristic scale depends on the mass scale of the quintessence around the minimum of the potential. We discuss the observability of these late-time cosmological signals through cosmic microwave background, quasar astrometry, pulsar timing arrays, and other observational probes.

PS1-12sk is a type Ibn supernova (SN) found at the host environment showing no obvious ongoing star formation, which challenges the massive star explosion scenario. We use the ejecta-circumstellar medium (CSM) interaction (CSI) and the CSI plus $^{56}$Ni models in the context of double white dwarf (WD) merger to fit the bolometric light curve (LC) of PS1-12sk, since the He emission lines at the photospheric phases indicated the interaction between the SN ejecta and He-rich CSM. We find that the CSI model failed to explain the LC, but the CSI plus $^{56}$Ni model can account for the bolometric LC. The derived masses of the two WDs and $^{56}$Ni are $\sim 0.70 M_\odot$, $\sim 0.40 M_\odot$, and $\sim 0.09\,M_\odot$, respectively. The facts that the ejecta mass ($\sim 0.984 M_\odot$) is well below the Chandrasekhar limit ($\sim 1.4 M_\odot$) and that the $^{56}$Ni mass is comparable to the $^{56}$Ni yields of the explosions of some sub-Chandrasekhar explosion models support the scenario that PS1-12sk might be from a sub-Chandrasekhar explosion induced by the merger of two low-mass WDs. The derived innermost radius ($\sim 13.81 \times 10^{12}$ cm) and the mass of the CSM ($\sim 0.116 M_\odot$) disfavor the possibility that the CSM was formed in the merger phase. We suggest that the flybys before the merger can account for the position and mass of the CSM.

We present a systematic investigation of the radiative properties and optical appearance of rotating SV black holes surrounded by a thin accretion disks, and mainly analyze the influence of the regularization parameter $g$ on related observables. The results show that although the kinetic quantities and the location of the innermost stable circular orbit (ISCO) depend on the regularization parameter $g$, the radiative efficiency of the rotating SV black hole is the same as its Kerr counterpart. Within the Novikov-Thorne thin-disk model, the radiative flux, effective temperature, and spectral luminosity are studied, and by adopting observational parameters relevant to SgrA* and M87*, concrete examples of the rotating SV black holes are calculated and compared with that of the Kerr black holes. The results show that the parameter $g$ suppresses the maximum values of these quantities. In addition, using a backward ray-tracing technique, we numerically simulate the optical appearance of rotating SV black holes and analyze the corresponding intensity images, redshift and observed flux distributions. Our results show that these quantities are affected by $g$. In particular, as $g$ increases, the observed intensity is significantly suppressed and the photon ring region has remarkable increase in its width. Our findings suggest that accretion-disk-related observables may provide important avenues to distinguish rotating SV black holes and Kerr black holes, and offer theoretical guidance for future high-resolution observations.

The James Webb Space Telescope (JWST) has enabled the discovery of hundreds of supermassive black holes (SMBHs) at redshifts $z\gtrsim 4-7$. A non-negligible fraction of these SMBHs are hosted in galaxies with BH-to-galaxy mass ratios ($M_{\rm BH}/M_\star$) being excessively larger than that for local SMBHs by $\sim 1-2$ dex. The origin of these ``overmassive'' BHs remains elusive, demanding either a heavy seed formation scenario or rapid growth of seed BHs. Their deviation from local scaling relations challenges our understanding of how SMBHs and their host galaxies coevolve across cosmic time. In this paper, we apply phenomenological modelings for BHs and galaxies to dark matter halo merger histories from N-body simulations to investigate the subsequent evolution of JWST-discovered ``overmassive'' SMBHs. We find that early evolution of ``overmassive'' SMBHs is dominated by stunted accretion leading to gradual decreases in $M_{\rm BH}/M_\star$ ratios. In contrast, less massive SMBHs experience super-Eddington accretion during their early evolution, resulting in a slow increase of mass ratios toward $M_{\rm BH}/M_\star \sim 0.01$. Convergence occurs at $M_{\rm BH}\sim 10^8~M_\odot$ with $M_{\rm BH}/M_\star \sim 0.01$. At lower redshift, nearly all SMBHs evolve onto local relations, as expected given that our models adopt empirical relations derived from low-redshift observations. This suggests that the global feedback mechanisms regulating the coevolution of $M_{\rm BH}/M_\star$ ratios are implicitly encoded in local relations in terms of star-formation rate distribution, black hole accretion rate distribution and their active (quiescent) fractions.

A model-agnostic description of the baryon acoustic oscillation (BAO) feature in redshift space requires a number of ingredients. Physically, one must describe the impact of cosmological bulk flows which progressively and anisotropically smear out the feature over time. One must also model the effects of the scale dependence of tracer bias and the mode coupling between short and long scales. All of these can be incorporated using the Zel'dovich approximation alone, without reference to any particular cosmological model. On the technical front, one needs a robust, complete and cosmology-independent basis to describe the shape of the real space BAO feature in linear theory, which can then be propagated to the nonlinearly evolved, measured feature in redshift space. In this work, we describe how these ingredients -- which we have systematically constructed in recent work -- come together in an accurate framework capable of describing the BAO-scale pairwise measurements of state-of-the-art galaxy surveys. Using mock observations and $N$-body simulations, we show that our template-free framework can produce unbiased and precise cosmological constraints for samples with realistic levels of nonlinearity. This work represents one of the final steps in constructing a data-ready analysis framework for model-agnostic cosmological inference from the BAO feature.

Foreground masking and incomplete sky coverage complicate CMB polarization analyses by inducing mode coupling and imperfect E/B separation, with particularly strong impact on searches for primordial $B$-modes. We present SkyReconNet-P, a convolutional neural network for inpainting CMB polarization maps that extends the SkyReconNet framework to jointly reconstruct the polarization $(Q,U)$ maps from partial-sky observations. The method combines regional processing with a hybrid design, utilizing standard convolution and dilated convolution to do a multiscale feature integration. We evaluate performance at both the map and power spectrum level using two masking scenarios: a generated random mask and the Planck 2018 common polarization inpainting mask. For both masking scenarios, SkyReconNet-P reproduces the large-scale morphology of the target maps. In power-spectrum space, we find that the reconstructed $E$-mode spectrum closely tracks the target at low multipoles, while small biases emerge at higher $\ell$. For $B$-mode, the raw reconstructed spectra exhibit a larger multipole-dependent bias, which we mitigate using a simulation-based linear calibration. We show that the calibrated $B$-mode spectrum preserve more information by comparing it with spectrum estimation using pseudo-$C_\ell$. Finally, we demonstrate cosmological parameter inference from calibrated reconstructed spectra by fitting $(r, A_{\rm lens})$ with a Gaussian bandpower likelihood, recovering posteriors consistent with injected parameters across three test ensembles down to $r \sim 10^{-3}$. These results support inpainting as a complementary route to cut-sky approaches when downstream pipelines can greatly benefit from statistically well-characterized, gap-filled polarization maps.

In a recent arXiv posting [1] I reported the result of an experiment: asking this http URL to compare three items concerning (ordinary) Gamma Ray Burts (GRBs): the data, the standard paradigm(s) and the "Cannonball" (CB) model. Here I ask the same URL to extend this comparison to long--lasting GRBs, binary Neutron-Star mergers and their associated short--hard GRBs, low--luminosity GRBs, X--ray flashes, X--ray transients, and non--solar cosmic rays. The results of this experiment are enlightening but worrisome. Except for this abstract, two footnotes and two other references to standard [2] and CB-model [3] articles and talks, all of what follows is, verbatim, what the cited AI "opines".

Magnetic fields are known to efficiently redistribute angular momentum in stars. They have been recently measured in the cores of red giant stars using asteroseismology. It was shown that core magnetic fields, if unaccounted for, can bias the measurements of the gravity offset $\epsilon_g$, which is otherwise well characterised for red giants. Exploiting this bias as a way to detect magnetic fields in the cores of red giants, we wish to increase the number of magnetic field detections in red giants, but also to establish a method that could be widely applicable to all red giants. We selected 218 Kepler red giants showing abnormal measured values of $\epsilon_g$. We used robust statistical criteria based on the expected lifetime of mixed modes to identify significant modes. We then adjusted an asymptotic expression for mixed modes to the observed frequencies, taking magnetic field and rotation into account, using Bayesian inference. We then assessed the probability of magnetic field detection and measured the magnetic field intensity using stellar models for the favourable cases. We found new magnetic red giant stars with fields ranging from 34 to 260 kG. For these stars, we measured values for $\epsilon_g$ now in agreement with the expected value. Adding these new detections to those of previous studies, we showed that the mass distribution of magnetic giants is similar to that of the complete catalogue of red giants, but different from the mass distribution of red giants with suppressed dipole modes. We also found that the core rotation of magnetic red giants follows a similar distribution as red giants in general. This could either mean that the detected fields do not have a predominant impact on the redistribution of angular momentum, or that other red giants also harbour an internal field that is currently non-detectable or has not yet been detected.

Context. Solid methane (CH$_4$) is an important molecule in interstellar and planetary environments, serving as a precursor to complex organic compounds and a potential biosignature in exoplanetary studies. Despite its significance, laboratory data on low-temperature phase of methane below 10 K remain limited. Aims. We aim to obtain spectra of methane in binary mixtures at 10 K and compare it to the spectra obtained at 6.7 K. These temperatures correspond to phases II and II* of pure methane and are representative of dark molecular clouds and protostars at early stages. We also aim to test the obtained data applicability to JWST data interpretation. Methods. Laboratory reference spectra were obtained on the ISEAge setup via FTIR spectroscopy in transmission mode. A weighted $\chi^2$ minimization is used for the fitting. Results. We present infrared spectra with corresponding band strengths of pure methane and binary mixtures with methane: CH$_4$:H$_2$O,CH$_4$:CO$_2$, CH$_4$:CH$_3$OH, CH$_4$:NH$_3$ at 6.7 K and 10 K showing a 20\% increase in mixtures compared to commonly used 10 K band strength value of pure methane. We also test the usability of the spectra on open JWST data by probing the spatial distribution of methane in B335. We also present additional experiments concerning the phase transition of methane between phase II* and phase II. Conclusions. Our results reveal distinct spectral features for methane in non-H$_2$O environments, enabling more accurate interpretation of JWST observations. The dataset of spectra, publicly available on Zenodo, can be used for fitting JWST data.

To explore the environmental mechanisms causing quenching in nearby star-forming galaxies, we study the variation with local and large-scale environments of a star formation concentration index, C-index $\equiv\log{(r_{50,{\rm H}\alpha}/r_{50,\rm cont}})$, that traces the spatially-resolved distribution of H$\alpha$ emission. Our analysis combines (i) GAMA spectroscopic redshift survey data to optically select galaxy groups and reconstruct the cosmic web, (ii) eROSITA data to identify X-ray-emitting groups, and (iii) SAMI Galaxy Survey data to characterise spatially-resolved star formation. We find that galaxies in X-ray+optical groups exhibit the lowest median C-index and the highest fraction of centrally-concentrated star-forming galaxies relative to optical groups and the field (independently of group or stellar mass). Star-forming galaxies in more X-ray luminous groups at fixed dynamical mass show more concentrated star formation. At large scales, nodes show the lowest median C-index and the highest fraction of centrally-concentrated star-forming galaxies relative to filaments and voids, which have similar C-index distributions. C-index correlates most strongly with the distance to the closest node, leaving no significant role for other local or large-scale environment metrics. Finally, regular star-forming galaxies tend to have spins aligned parallel to filaments, consistent with smooth gas accretion, while centrally-concentrated galaxies tend have spins aligned perpendicular to filaments, likely driven by mergers and associated with bulge growth. These results suggest that multi-scale environmental processes, i.e. locally and at large-scale, act to concentrate star formation toward galaxy centres, via gas-related mechanisms in nodes and ram-pressure stripping in X-ray+optical groups.

Cosmic voids contain higher-order cosmological information and are of interest for astroparticle physics. Finding genuine matter underdensities in sparse galaxy surveys is, however, an underconstrained problem. Traditional void finding algorithms produce deterministic void catalogs, neglecting the probabilistic nature of the problem. We present a method to sample from the stochastic mapping from galaxy catalogs to arbitrary void definitions. Our algorithm uses a deep graph neural network to evolve "test particles" according to a flow-matching objective. We demonstrate the method in a simplified example setting but outline steps to generalize it towards practically usable void finders. Trained on a deterministic teacher, the model performs well but has considerable stochasticity which we interpret as regularization. Cosmological information in the predicted void catalogs outperforms the teacher. On the one hand, our method can cheaply emulate existing void finders with apparently useful regularization. More importantly, it also allows us to find the Bayes-optimal mapping between observed galaxies and any void definition. This includes definitions operating at the level of simulated matter density and velocity fields.

HESS J0852-463 is a TeV {\gamma}-ray source located in the Galactic plane. The region consists of a supernova remnant (SNR, RX J0852.0-4622) with a shell-like morphology, commonly referred to as Vela Junior, and a pulsar denoted PSR J0855-4644. Pulsars are among the most efficient leptonic accelerators in our Galaxy, making this region particularly interesting to study. We utilise the most recent data taken by the High Energy Stereoscopic System (H.E.S.S.), to investigate any {\gamma}-ray emission associated with the pulsar in this region, PSR J0855-4644. We applied a full forward folding method on the H.E.S.S. data. Utilising 3D modelling techniques, we evaluated the TeV {\gamma}-ray emission towards the various components of this complex system. The distinct energy-dependent morphology observed in our data motivates further investigation of this source. We resolved the emission in the Vela Junior region into various components, several of which correspond to the SNR itself. In particular, we find a new extended component which is coincident with the position of PSR J0855-4644. The spectrum follows a power-law distribution with a best-fit index of {\Gamma}E = 1.81 \pm 0.07stat which differs from the properties of the surrounding {\gamma}-ray emission of the Vela Junior SNR. A one-zone leptonic joint fit between the X-rays (from XMM-Newton) and {\gamma}-rays (from H.E.S.S.) leads to a lower limit on the magnetic field of 1.6{\mu}G and a spectral index of {\alpha} = 1.88 \pm 0.01, in line with expectations of pulsar wind nebulae (PWNe). In this paper, we report the first detection of the PWN of PSR J0855-4644 at TeV energies with the H.E.S.S. experiment at a significance of 12.2{\sigma}. This is attributed to the advanced techniques of the 3D analysis. Based on the pulsar's characteristics, its PWN is consistent with the known TeV PWNe population in the Galaxy.

4U 1556-60 is a low-mass X-ray binary that was discovered more than 50 years ago as a persistent X-ray source; however, very little was known about it. Recently, Gaia obtained a parallax for the optical counterpart that places 4U 1556-60 at a distance of only about 700 pc, making it one of the closest X-ray binaries known to date. This close distance drastically alters what was previously assumed about the source. We revisit 4U 1556-60 in light of the newly determined distance of 700 pc, reinterpreting its literature and presenting new X-ray and radio observations to better understand various characteristics of the system. We conclude that a scenario in which 4U 1556-60 is a candidate ultracompact neutron star X-ray binary at a distance of ~700 pc is able to explain the observed properties of the source. It resides at a persistent X-ray luminosity of ~2x10^34 erg/s, an unusual value for a typical X-ray binary, but similar to several ultracompact systems. The ratio of the X-ray to optical luminosity is very high, also suggesting a physically small accretion disk. The radio jet is undetected with a very deep upper limit of 3x10^25 erg/s, which is about 10^3 times fainter than the expected black hole jet correlation, strongly indicating a neutron star accretor. The X-ray spectrum is dominated by a power law, and the X-ray timing properties are also consistent with observations of other very low accretion rate X-ray binaries. No spin or orbital periodicity are found in the X-ray data. Future observations, especially to determine its orbital period, will further aid in understanding 4U 1556-60.

Quantum gravity theories often modify spacetime symmetries. In particular, Lorentz invariance may be violated when approaching the Planck scale. Although the scales at which interactions occur in extensive air showers induced by ultra-high-energy cosmic rays in the atmosphere are many orders of magnitude below the Planck scale, these violations might still be observable. In this work, the fluctuations in the number of muons in the extensive air showers measured at the Pierre Auger Observatory are exploited, for the first time, to constrain Lorentz invariance violations. The bounds derived in the hadronic sector are the strongest ever obtained, and do not rely on assumptions about the mass composition of ultra-high-energy cosmic rays. The fluctuations in the number of muons constitute a new and powerful observable to further explore Lorentz invariance in a region of the parameter space not accessible to other observables.

We investigate whether a non-minimal coupling between dark matter and gravity can influence the internal dynamics of ultra-diffuse galaxies. Within this framework, the gravitational potential is modified by an additional term that captures the interaction between spacetime curvature and the dark matter with a coupling constant determined by a length scale $L$. Using spherical Jeans modeling, we analyze the kinematic data of three ultra-diffuse galaxies namely: NGC\;1052-DF2, NGC\;1052-DF4, and Dragonfly\;44, which span the observational extremes from dark matter deficient to dark matter dominated systems. For each galaxy we explore several dark matter halo profiles, two orbital anisotropy models, and both with and without Stellar-to-Halo Mass Relation scenarios, and we perform a Bayesian parameter inference. Our results show that across all the considered configurations, the constrained astrophysical parameters are consistent with standard ones from General Relativity. The posterior distributions of $L$ show no preference for non-zero values and result only in upper limits, suggesting that any non-minimal coupling contribution must be small and perturbative on this scales. Future high precision velocity measurements will be essential to determine whether non-minimal coupling effects can become observationally distinguishable in low-acceleration systems.

Sensitivity of searches for extended emission around gamma-ray sources is naturally limited by the precision of the knowledge of the Point Spread Function (PSF) of gamma-ray telescopes. Inaccuracies in the PSF models of the Fermi Large Area Telescope (LAT) can potentially lead to false positive detections of source extension. We explore uncertainties in the Fermi/LAT PSF by comparing the PSF models provided by the Fermi/LAT Instrument Response Functions (IRFs) with signals of bright pulsars. We compare the analytical PSF models of Fermi/LAT IRFs with pulsar data and fit the pulsar data with the same analytical model as in the Fermi/LAT IRFs to derive an improved set of PSF parameters. We then apply this revised PSF parameterisation to the search of extended emission around a blazar, Mrk 501. We find that the parameters of the analytical PSF models of Fermi/LAT IRFs are inconsistent with the pulsar data. We obtain an improved set of PSF parameters from the fits to pulsar data that is consistent with observations. We find no evidence of the previously reported extended signal around Mrk 501 if the revised PSF consistent with pulsar data is used in data analysis.

Miras and other long-period variable (LPV) stars on the AGB follow period-luminosity (PL) relations. These relations have been difficult to study for Galactic LPVs because their distances were poorly known in the past. We aim to establish the PL relations of solar-neighbourhood Miras for several near-IR photometric bands. To this end, we used multi-epoch photometry from the DIRBE and unTimely/WISE catalogues, Gaia parallax distances, and contemporary pulsation periods obtained from optical observations of a well-selected sample of solar-neighbourhood Miras. We show that clearly defined PL relations in the nine investigated near-IR bands emerge from our data, and we report the slopes and zero-point magnitudes. We find that Galactic Miras are fainter in the near-IR than their Large Magellanic Cloud siblings. We derive average period-temperature, period-bolometric-luminosity, and period-radius relations from fits to synthetic SEDs constructed from the PL relations. By applying AGB evolutionary models, the scatter of stars around the PL sequences can also be used to test whether the parallax uncertainties quoted in the Gaia catalogue are realistic. Furthermore, we performed such tests based on a comparison with parallaxes obtained with the VLBI and with a sample of LPVs in the globular cluster 47 Tuc. We conclude that, for Galactic Miras with a fractional parallax uncertainty of <0.1 in the Gaia catalogue, the parallax uncertainty is underestimated by factors between 1.0 and 1.7, and most likely by $\sim1.3$. For more uncertain parallaxes, we find evidence that the distances (parallaxes) are generally overestimated (underestimated). Nevertheless, we find strong evidence that the large error-inflation factors reported for AGB stars in the literature are unrealistic. Our results lend confidence to the parallax measurements of these highly extended, variable stars.

Following the detections of the first extragalactic ``Luminous Red Nova'' (LRN) M31 RV in 1989, and its first Galactic counterpart V4332~Sgr in 1994, there have been many discoveries of similar, or closely related, objects. They are important because they bridge the luminosity gap between the brightest novae and supernovae, a largely unexplored parameter space. The cause of eruptions in LRNe is still unclear, a stellar merger being the most favored mechanism. However, barring V1309~Sco, there has been no direct evidence for a merger in the other objects. Here we present strong evidence that V4332~Sgr was a merger event. High resolution infrared observations of the CO fundamental band show an unusually small $^{12}$C/$^{13}$C ratio of $3.5\pm1$. This indicates that a violent event had occurred, whose effects penetrated deep enough to allow CNO cycle processed $^{13}$C in the inner H burning shell to be brought to the surface. We rule out planetary ingestion, and propose that the eruption was due to a merger between V4332~Sgr and a companion star. It is shown that V4332~Sgr was likely surrounded by an edge-on disk before its eruption. If this disk was a flattened common envelope containing V4332~Sgr and a companion star, then a merger scenario would not be inconsistent. Furthermore, V4332~Sgr had multiple outbursts, previously unreported but an important piece of information, since multiple outbursts are a trait shared by many LRNe.

Astrochemical observations have revealed a surprisingly high level of chemical complexity, including long carbon chains, in the earliest stages of Sun-like star formation. The origin of these species and whether they undergo further growth, possibly contributing to the molecular complexity of planetary systems, remain open questions. We present recent observations performed using the 100-m Green Bank Telescope of the prestellar core L1544 and the protostellar system IRAS 16293-2422. In L1544, we detected several complex carbon-bearing species, including $\mathrm{C_2S}$, $\mathrm{C_3S}$, $\mathrm{C_3N}$, $\mathrm{c\text{-}C_3H}$, $\mathrm{C_4H}$, and $\mathrm{C_6H}$, complementing previously reported emission of cyanopolyynes. In IRAS 16293-2422, we detected $\mathrm{c\text{-}C_3H}$ and, for the first time, $\mathrm{HC_7N}$. Thanks to the high spectral resolution, we refine the rest frequencies of several $\mathrm{c\text{-}C_3H}$ and $\mathrm{C_6H}$ transitions. We perform radiative transfer analysis, highlighting a chemical difference between the two sources: IRAS 16293-2422 shows column densities 10-100 times lower than L1544. We perform astrochemical modeling, employing an up-to-date chemical network with revised reaction rates. The models reproduce the general trends, with cyanopolyyne and polyynyl radical abundances decreasing as molecular size increases, but they underestimate the abundances of cyanopolyynes longer than $\mathrm{HC_5N}$ by up to two orders of magnitude. Current models, which include the dominant neutral-neutral formation routes, cannot account for this discrepancy, suggesting that the chemical network is incomplete. We propose that additional ion-molecule reactions are crucial for the formation of these species. Developing a more comprehensive chemical network for long carbon chains is essential for accurately interpreting present and future observations.

We summarize the past four decades of astrophysics and exoplanet direct imaging mission concept studies, technology developments, and scientific progress that have led to the initiation of the Habitable Worlds Observatory project by NASA.

The angular, harmonic cross-correlation between gravitational wave (GW) events and galaxy catalogues contains rich information on the large-scale structure and the origin of compact binary mergers. In this work, we study how uncertainties in the binary formation channel affect the predicted cross-correlation signal for both current-generation and next-generation networks of detectors. We generate five mock GW catalogues for which we vary the progenitor-to-remnant mass-transfer function and the time-delay probability distribution between progenitor and remnant. We then cross-correlate these catalogues with galaxy samples modelled on the 2MASS Photometric Redshift catalogue (2MPZ) and the Gaia-unWISE quasar catalogue (Quaia). We find that the mass-transfer function has negligible effect on the cross-correlation signal, with differences remaining within redshift uncertainties. In contrast, the time-delay distribution dramatically affects the redshift distribution of the GW events and, with it, the cross-correlation signal, particularly for shallow galaxy catalogues. In particular, current-generation facilities can achieve significant detections only for the longest time delays when cross-correlated with 2MPZ, whilst all cross-correlations with the deeper Quaia catalogue are marginally detectable or consistent with zero. Our exploratory results thus demonstrate that forecasts on cosmological or astrophysical parameters derived from GW-galaxy cross-correlations are, as expected, strongly sensitive to the assumed binary formation history.

It has long been suspected that black hole accretion-outflow coupling is invariant from the stellar to supermassive scales. Stellar mass black hole accretion flows are known to launch jets and outflows as they transition through critical accretion rate thresholds, with values well constrained observationally owing to their short evolutionary timescales. In contrast, accretion flows in typical supermassive black hole (SMBH) systems (those in active galactic nuclei) evolve over thousands of years, making the critical transitions at which jets are launched impossible to constrain in individual systems. Tidal disruption events (TDEs) provide the unique opportunity to witness the birth and evolution of an accretion flow onto a SMBH which evolves on timescales of years. Here we show that TDEs launch outflows during a super-Eddington accretion phase and a second, physically distinct outflow, at a critical accretion rate of $L_{\rm crit} \approx0.02$ $L_{\rm Edd}$, the same as the critical accretion rate for state transitions observed in accreting stellar mass black holes. This work naturally explains the mechanism, observed properties, and detection rate for prompt and delayed outflows observed in TDEs, which until now have been open problems. More broadly, we demonstrate that SMBHs exhibit the same accretion-outflow coupling as stellar mass black holes and that the critical low accretion rate threshold for jet formation in black holes is scale invariant.

We present a comprehensive validation analysis of the exoplanet candidate TOI 7475.01 (TIC 376866659), detected by the TESS mission. Using a custom pipeline combining natural flux preservation with robust BLS detection, we identified a transit signal with a period of $3.2538$ days and a depth of $\sim 4600$ ppm. To rule out false positives, we performed centroid analysis, spatial contamination checks using Gaia DR3, and a statistical validation using triceratops. Our results show a Signal-to-Noise Ratio (SNR) of 294.13 and a False Positive Probability (FPP) of $\approx 0$. Based on the clean spatial environment, stable centroids, and high statistical probability, we validate TOI 7475.01 as a planetary companion.

We investigate the bar fraction in galaxy pairs from the SDSS to assess how galaxy interactions affect bar structures. Compared to isolated galaxies, close pairs exhibit a significantly reduced bar fraction at projected separations within 25 kpc. This reduction is driven almost entirely by systems showing clear merger or disturbance signatures, indicating that tidal interactions suppress bars. The decline is dominated by a decrease in weak bars, while the fraction of strong bars remains largely unchanged. Bar suppression is primarily associated with major mergers and is strongest in massive host galaxies. A weaker but statistically significant suppression is detected in minor mergers only for massive galaxies with small bulges. In contrast, no significant dependence of bar suppression on the relative orientation between pair members is found. These findings provide observational evidence that tidal perturbations in major mergers play a key role in regulating bar evolution.

NGC 4631 X-8 is an ultraluminous X-ray pulsar (ULXP) having a spin period of about 9.7 s, discovered using XMM-Newton observations in 2025. The pulsar is known to show one of the largest spin-up rates ($\sim 9.6 \times 10^{-8}$ s s$^{-1}$) among the ULXP population. We explore the surface magnetic field of the neutron star in this source using different models, and find that the inferred magnetic field lies in the range of about $0.3-2 \times 10^{14}$G. We study the long-term magnetic field and spin period evolution of the pulsar assuming steady accretion using prevalent theoretical mechanisms and find that the pulsar will evolve to become a millisecond pulsar having decayed magnetic field of about $\sim 10^{9}$G in about a million years. The scenario of the formation of a millisecond pulsar is also probed using an estimate of the super-Eddington duty cycle of about 14% from the literature, which suggests that the neutron star would accrete sufficient matter to become a recycled millisecond pulsar. Exploring the magnetic field as well as the spin period evolution properties of ULXPs may enable us to understand the poorly understood evolutionary features of ULXPs, shed light on one of the pathways of millisecond pulsar formation and also help us to understand transient super-Eddington accretion phases in newborn magnetars, which are believed to power energetic events such as long gamma-ray bursts and Type I superluminous supernovae.

We present imaging and spectroscopic analyses of Chandra and XMM-Newton observations of ACT-CL J0123.5$-$0428, one of the most massive, highest redshift galaxy clusters detected within the survey fields of the Atacama Cosmology Telescope. The Chandra data are sufficient to characterize the morphology of this cluster and constrain the geometrically deprojected temperature in 2 spatial bins out to $r_{2500}$, revealing a dynamically relaxed system whose temperature drops to $kT = 1.8\pm0.6$ keV in the inner $\sim40$kpc. Within this same inner radius, the surface brightness and density of the ICM is sharply peaked, and the cooling time falls to $t_\mathrm{cool}=280^{+150}_{-120}$ Myr. A novel forward-modeling analysis of the XMM data extends imaging and spectroscopic measurements of this system out to $r_{500}$, constraining the redshift to $z=1.50\pm0.03$, with a mean temperature of $kT = 7.3\pm1.1$ keV and an emission-weighted mean metallicity of $Z/Z_\odot = 0.43^{+0.46}_{-0.25}$. We also utilize the limited optical/IR photometric coverage of the cluster to characterize the properties of the brightest cluster galaxy (BCG), which is coincident with the X-ray peak. Despite the high redshift and strong cool core, the BCG exhibits no signs of recent or ongoing star formation, suggesting AGN feedback has been acting persistently to stem star formation since $z\sim 2.5$. These measurements identify ACT-CL J0123.5$-$0428 as the highest redshift, dynamically relaxed, cool core galaxy cluster discovered to date, making it a premier target for future astrophysical and cosmological studies.

The sensitivities of STIS first-order, medium resolution modes were redetermined from on-orbit observations and CALSPEC models (version 11) of the primary white-dwarf spectrophotometric standard stars G191-B2B, GD 71, and GD 153. The sensitivity of an additional configuration was updated by comparing observations of the secondary standard BD +75 325 with the STIS low-resolution spectrum that has been calibrated consistently with the version 11 models. The procedures used to derive the sensitivities and verify the PHOTTAB reference files prior to their activation in CRDS (on May 1, 2025) are described. Results are presented in graphical form in an extensive appendix. Issues and uncertainties are discussed briefly, along with recommendations for future work.

TiO plays an important role in determining the atmospheric structure of M-type stars and in shaping the visual part of their spectra. We compute synthetic spectra for late-type M giants and identify systematic discrepancies in the wavelength range 5810-5850 A. To investigate the origin of these discrepancies, we analyse experimental TiO absorption cross-sections. We report the detection of a molecular band of the singlet system $^{1}\Pi$-a$^{1}\Delta$ of the TiO, with an R-head located at 5814.8 A, overlapping the 1-3 band of the C$^3\Delta$-X$^3\Delta$ system. The lower state of the band is identified as the a$^1\Delta$ v"=0 state, while the upper state is most likely the $^1\Pi$ in its ground vibrational state. The empirical band intensity is derived by comparing the relative strengths of neighbouring bands from the C$^3\Delta$-X$^3\Delta$ and B$^3\Pi$-X$^3\Delta$ systems in the experimental cross-section. The band intensity is further validated by synthetic spectrum calculations for the late-type giant 30 Her (M6 III) and comparison with its observed spectrum from the MELCHIORS library. The newly identified band is sufficiently strong to affect the flux distribution in the spectra of cool stars.

The origin of the red colouration in Jupiter's Great Red Spot (GRS) is a long-standing question in planetary science. While several candidate chromophores have been proposed, no clear conclusions have been reached regarding its nature, evolution, or relationship to atmospheric dynamics. In this work, we perform microphysical simulations of the reddish haze over the GRS and quantify the production rates and timescales required to sustain it. Matching the previously reported chromophore column mass and effective radius in the GRS requires column-integrated injection fluxes in the range $1\times10^{-12}$-$7\times10^{-12}$ kg m$^{-2}$ s$^{-1}$, under low upwelling velocities in the upper troposphere ($v_{\mathrm{trop}}\lesssim 1.5\times10^{-4}$ m s$^{-1}$) and particle charges of at least 20 electrons per $\mu$m. Such rates exceed the mass flux that standard photochemical models of Jupiter currently supply via NH$_3$-C$_2$H$_2$ photochemistry at 0.1-0.2 bar, the most popular chromophore pathway in recent literature. We find a lower limit of 7 years on the haze formation time. We also assess commonly used size and vertical distribution parameterisations for the chromophore haze, finding that eddy diffusion prevents the long-term confinement of a thin layer and that the extinction is dominated by particles that can be represented by a single log-normal size distribution.

The [C/N]-age relation has become a powerful tool for reconstructing the formation history of the Milky Way (MW), providing the largest age sample for field giant stars. However, at metallicities below [Fe/H] $< -0.5$, stellar surfaces are altered by a poorly understood process known as extra mixing, which modifies [C/N] in a mass- and metallicity-dependent manner. This effect complicates the application of the traditional [C/N]-age relation in metal-poor regimes. Within the MW, constraining the mass dependence of extra mixing is particularly challenging because stars at [Fe/H] $< -0.5$ are predominantly old and therefore low-mass, leading to strong degeneracies between mass and metallicity. In this work, we explore the potential of the Magellanic Clouds (MCs) to disentangle these effects and constrain extra mixing as a function of age and metallicity. By comparing empirical corrections calibrated in the MW with predictions from thermohaline mixing models, we isolate the mass dependence of extra mixing in the MCs down to [Fe/H] $\sim-0.7$. We find that the empirical calibration performs well for lower-mass stars ($< 1.25$ $M_{\odot}$), while theoretical models successfully reproduce the observed mass dependence down to $\sim$ 1.25 $M_{\odot}$. We further present the first observational evidence that extra mixing becomes ineffective above $\sim$ 1.8 $M_{\odot}$ at [Fe/H] $\sim -0.7$. Our results demonstrate the feasibility of deriving [C/N]-based ages for individual stars in external galaxies. Future observations targeting higher-$\log g$ or fainter stars in the MCs will provide stronger constraints on extra-mixing processes and enable the calibration of [C/N]-age relation that can be applied to low-metallicity individual stars in the MW or external galaxies.

The results of interferometric observations of the star-forming region S255IR in the frequency range 210--250 GHz are presented. The observations were carried out with the antenna array SMA (Hawaii, USA). Fifty-three molecules were detected, including complex organic molecules (COMs) such as CH$_3$CHO, CH$_3$CN, CH$_3$CH$_2$CN, and many others. Typical rotational temperatures in the hot core SMA1 fall in the range 100--200 K. Optical depths in the lines of methanol and some other molecules in the cores SMA1 and SMA2 were estimated. In SMA1, the optical depth of one of the strongest methanol lines, $5_{-1}-4_{-1}E$, proved to be $23.8 \pm 1.5$. Based on this value, one can assume that the lines of other oxygen-containing COMs, such as CH$_3$OCHO, CH$_3$OCH$_3$, CH$_3$CH$_2$OH, which are typically much less abundant in hot cores than methanol, are optically thin in SMA1. Most of the detected molecules can be roughly divided into two groups. The molecules of the first group emit exclusively toward the hot core SMA1, while some or all lines of the molecules of the second group, in addition to SMA1, can be seen toward a ring-like structure to the west of SMA1. This structure is most likely associated with the walls of a cavity formed by high-velocity outflows driven by young stellar objects (YSOs) in molecular cores SMA1, SMA2, and possibly SMA3. The gas temperature and density in the cavity walls were estimated using methanol lines. The temperature was found to be about 50--60 K, and the density about $10^7-10^8$ cm$^{-3}$. The column density of methanol near the brightness peaks in the lines of this molecule is about $5\times 10^{15}$~cm$^{-2}$. The column densities of other COMs in the ring-like structure will be determined in future studies with increased sensitivity achieved by spectral line stacking.

Spectral distortion anisotropies of the cosmic microwave background (CMB) provide a new probe of the early Universe that can be accessed using traditional CMB imaging techniques. It is possible to compute the creation and evolution of anisotropic signals for various scenarios using the frequency hierarchy method recently developed for CosmoTherm. However, the current treatment is not perfect and some approximations had to be made. Here, we carefully construct a modified form for the evolution equations that has the full equilibrium solutions built into the formulation. We improve the formalism to account for i) additional stimulated scattering effects, ii) kinematic corrections to the thermalization terms, iii) corrections to the standard perturbation variables and iv) direct photon sources. These effect could not be captured with the original formulation of the frequency hierarchy method but are indeed important for cleanly separating real distortions from temperature signals. However, we show that previous results are not altered significantly when compared to the improved formulation presented here. As a new worked example, which could indeed not be treated before, we also illustrate how possible changes in the temperature-redshift relation would create spectral distortion anisotropies in the pre-recombination era. The theoretical methods presented here are also an important step towards being able to consistently predict the CMB spectral distortion anisotropies in photon-dark photon and photon-axion conversion scenarios.

More than 150 measurements of the HeII 4686A emission line in spectra of SS433 were obtained during 388 nights in 2020-2025 with the Transient Double-beam Spectrograph on the 2.5 m telescope of Caucasian Mountain Observatory of Sternberg Astronomical Institute. We found that the HeII emission line formation region is not eclipsed and is significantly larger than both the donor star and the photosphere of the supercritical accretion disk. The HeII radial velocity curve was found to be independent of the precessional phase and inconsistent with the photometric curve. These findings suggest that the HeII line does not reflect the orbital motion of the compact object. Therefore, spectroscopic estimates of the masses of the components in SS433 based on the HeII emission line can be unrealistic.

ISRO's terrain characterization and hazard mapping from Chandrayaan-2 Orbiter High Resolution Camera (OHRC) stereo imagery were central to the safe landing of Chandrayaan-3 - the first successful landing in the lunar south polar region. However, these elevation products were generated with a proprietary pipeline and have not been publicly released. We present a 0.30 m/pixel digital elevation model (DEM) of the Chandrayaan-3 Vikram landing site using a fully open workflow based on ISIS, the Ames Stereo Pipeline, and ALE, achieving sub-meter resolution comparable to mission-reported products. The reconstruction covers 2.18 x 2.24 km with 91.2% valid pixel coverage, 8.1 cm median triangulation error, and 40-50 cm relative vertical precision. The Vikram lander and Pragyan rover are individually resolved. Geodetic alignment to an LROC NAC stereo DEM achieves approximately 30 m horizontal accuracy; pixel-wise validation at 3 m resolution confirms negligible vertical bias (median dz = +0.28 m) and robust dispersion (NMAD = 2.88 m). Stable OHRC stereo convergence requires Community Sensor Model (CSM) camera models; the legacy ISIS camera model failed across two independent sites. At 0.30 m, these DEMs complement LROC NAC DTMs (approximately 1 m), resolving sub-meter hazards below the NAC detection threshold. Applied to the extensive OHRC south polar archive, this methodology provides independent capability for hazard mapping and landing site analysis for upcoming missions including Chandrayaan-4, LUPEX, and Artemis.

Cross-survey generalization is a critical challenge in stellar spectral analysis, particularly in cases such as transferring from low- to moderate-resolution surveys. We investigate this problem using pre-trained models, focusing on simple neural networks such as multilayer perceptrons (MLPs), with a case study transferring from LAMOST low-resolution spectra (LRS) to DESI medium-resolution spectra (MRS). Specifically, we pre-train MLPs on either LRS or their embeddings and fine-tune them for application to DESI stellar spectra. We compare MLPs trained directly on spectra with those trained on embeddings derived from transformer-based models (self-supervised foundation models pre-trained for multiple downstream tasks). We also evaluate different fine-tuning strategies, including residual-head adapters, LoRA, and full fine-tuning. We find that MLPs pre-trained on LAMOST LRS achieve strong performance, even without fine-tuning, and that modest fine-tuning with DESI spectra further improves the results. For iron abundance, embeddings from a transformer-based model yield advantages in the metal-rich ([Fe/H] > -1.0) regime, but underperform in the metal-poor regime compared to MLPs trained directly on LRS. We also show that the optimal fine-tuning strategy depends on the specific stellar parameter under consideration. These results highlight that simple pre-trained MLPs can provide competitive cross-survey generalization, while the role of spectral foundation models for cross-survey stellar parameter estimation requires further exploration.

Gravitational waves (GWs) produced during first-order phase transitions (FOPTs) in the early universe provide a powerful probe of nonstandard cosmological histories. We study GW production from a FOPT ending a kination-dominated epoch in the Kination-Induced Big Bang scenario, in which a period of kination domination terminates through a phase transition that reheats the universe into radiation domination. A rolling scalar field drives the kination epoch. In the specific model we consider, its derivative coupling to a second scalar (tunneling field) dynamically traps the latter in a false vacuum, with the phase transition triggered as the kination field slows due to Hubble friction. We compute the resulting stochastic GW background from bubble nucleation and collisions, presenting analytic estimates and numerical results for the peak amplitude and frequency. In all cases we find an upper bound $\Omega_{\rm GW} h^2\lesssim 2\times10^{-7}$ from the bubble percolation condition. In the case where the false vacuum energy dominates at the transition (yet the kination field drives the FOPT), we find $\Omega_{\rm GW}h^2\gtrsim 10^{-12}$. We further find that the Hubble scale during the phase transition across a broad set of model parameters is bounded by $\mathscr{O}(10^{-13})M^2/M_{\rm Pl}\lesssim H_* \lesssim \mathscr{O}(0.1)M^2/M_{\rm Pl}$, where $M$ is the mass-scale controlling the strength of the interaction between the kination and tunneling fields. The predicted signal spans frequencies from nHz to MHz, allowing the model to explain the signal reported by Pulsar Timing Array experiments and to be constrained or probed by interferometers such as LISA, Advanced LIGO, Cosmic Explorer, and BBO. Interestingly, a FOPT can occur even if the bare tunneling potential has a single minimum, as metastability is generated dynamically by the coupling between the tunneling and the kination field.

Long-period radio transients (LPTs) are a new class of radio sources characterized by long spin periods ($P_{\text{spin}}>10^3$ s) and highly variable radio emission. While known magnetars are relatively young ($\tau<10^5$ yrs) with spin periods clustered between $1-10$ sec, it has been proposed that LPTs may be linked to a missing population of older magnetars. In this paper, we present an extensive parametric analysis of isolated magnetar spin evolution using various propeller spin-down models. In general, at higher initial magnetar B-fields ($B_0>\sim10^{15}$ G) and larger ambient densities ($n_0>\sim10^2$ cm$^{-3}$), magnetars will transition to the propeller phase earlier, and they start accreting gas from the ISM or molecular clouds after $\tau\sim10^8$ yrs. We found that a transition from the pulsar to the propeller phase is required to reach the observed LPT period range of $P>10^3$ s. More specifically, our population synthesis study based on Monte-Carlo simulations shows that two propeller models can account for most of the observed LPT periods ($P\sim1-400$ [min]) and their period derivative constraints ($\dot{P}<10^{-9}$ s s$^{-1}$). Our spin-down models predict that (1) nearby radio-quiet neutron stars with the estimated dipole $B$-field range of $B\sim(1-5)\times10^{13}$ G will transition to the propeller phase eventually after $\tau>\sim10^7$ yrs; (2) thermal X-ray emission from accretion-phase magnetars becomes too faint for detection after traveling ($d>\sim10$ kpc) from their birth places; (3) sporadic radio outbursts observed from LPTs may not be explained by regular radio pulsar and magnetar emission mechanisms that operate during the propeller phase.

Primordial black holes (PBHs), produced in the early Universe, can source a stochastic background of induced gravitational waves (GWs) and provide a non-thermal origin for dark matter (DM). We investigate DM production in a PBH-dominated cosmological framework, including contributions from PBH evaporation, gravitational production, and thermal freeze-in and freeze-out mechanisms, and determine the regions consistent with the observed DM relic abundance. We find that thermal freeze-in can compensate for the underabundance of PBH-sourced DM, while indirect detection remains largely insensitive due to the feeble interaction strength, making future GW observatories such as LISA and the Einstein Telescope (ET) unique probes of this scenario. For freeze-out DM, indirect detection experiments constrain regions with relatively large annihilation cross-sections, whereas GW observations probe complementary regions with heavier DM masses and smaller interaction strengths. Consequently, the same DM parameter space cannot be simultaneously probed by both indirect detection searches and GW missions. These results establish GW observations as a powerful and independent probe of DM production in PBH-dominated cosmologies, opening a new observational window into DM properties and the thermal history of the pre-BBN Universe.

In this paper, we revisit the linear perturbation (including the comoving curvature perturbation and field perturbation) dynamics in the SR-USR-SR inflation with instantaneous transitions. Using the junction method and asymptotic expansions of Hankel functions, we derive accurate asymptotic expressions for the time evolution of mode functions and the resulting power spectrum, based on three systematic rules for identifying the dominant terms across transitions. Our results reveal that a finite dip of the final power spectrum arises from the cancellation between two growing modes within the linear perturbation theory, rather than between constant and growing terms as previously suggested. We also provide analytical descriptions of the amplitude enhancement and oscillatory features in the linear power spectrum, in agreement with numerical computations. These simple, tractable formulas not only facilitate theoretical calculations but also yield testable predictions for future CMB observations.

Artificial Intelligence (AI) has transformed robotics, healthcare, industry, and scientific discovery, yet a major frontier may lie beyond Earth. Space exploration and settlement offer vast environments and resources, but impose constraints unmatched on Earth: delayed/intermittent communications, extreme resource scarcity, heterogeneous expertise, and strict safety, accountability, and command authority. The key challenge is auditable coordination among specialised humans, robots, and digital services in a safety-critical system-of-systems. We introduce Agent Mars, an open, end-to-end multi-agent simulation framework for Mars base operations. Agent Mars formalises a realistic organisation with a 93-agent roster across seven layers of command and execution (human roles and physical assets), enabling base-scale studies beyond toy settings. It implements hierarchical and cross-layer coordination that preserves chain-of-command while allowing vetted cross-layer exchanges with audit trails; supports dynamic role handover with automatic failover under outages; and enables phase-dependent leadership for routine operations, emergencies, and science campaigns. Agent Mars further models mission-critical mechanisms-scenario-aware short/long-horizon memory, configurable propose-vote consensus, and translator-mediated heterogeneous protocols-to capture how teams align under stress. To quantify behaviour, we propose the Agent Mars Performance Index (AMPI), an interpretable composite score with diagnostic sub-metrics. Across 13 reproducible Mars-relevant operational scripts, Agent Mars reveals coordination trade-offs and identifies regimes where curated cross-layer collaboration and functional leadership reduce overhead without sacrificing reliability. Agent Mars provides a benchmarkable, auditable foundation for Space AI.

Cold dark matter may be a fluid (or plasma) residing in a strongly-interacting hidden sector, rather than a population of weakly-coupled particles. Such a scenario admits a holographic description in terms of a cosmological braneworld embedded in the linear dilaton five-dimensional (5D) spacetime. In this framework, dark matter originates from the linear dilaton bulk black hole, whose phase we show to be thermodynamically favored at all temperatures. We present a natural freeze-in mechanism for the production of holographic dark matter, in which the bulk black hole is fed by energy leaking from the brane after inflation. Our model is characterized by two free parameters, one of which, the position of the black hole horizon, is fixed by the observed dark matter abundance. The remaining parameter, the 5D Planck scale $M_5$, is consistent with all current experimental bounds provided that $M_5\gtrsim 3\times 10^5$ TeV.

We investigate magnetohydrodynamic (MHD) turbulence in plane shear flows with a streamwise background magnetic field in the super-Alfvénic regime. We show that the large-scale velocity shear suppresses turbulence imbalance, driving the system toward a balanced state -- the energies of counter-propagating Alfvén waves become essentially equal, even at initially perfectly imbalanced Alfvénic turbulence. This balancing is due to the shear-induced linear non-modal dynamics of Alfvén waves, including their transient growth and over-reflection. This linear route to balancing turbulence is new -- fundamentally different from nonlinear ones operative in shearless MHD turbulence -- and have direct implications for understanding balanced/imbalanced MHD turbulence in the solar wind, which is modeled as a shear flow in a thermodynamically complex plasma.

The four-spacecraft formation is essential for measurements of various physical fields. The use of this formation on substantially elliptical heliocentric Kepler orbits allows measuring gradients of gravitation field in Solar system. The accuracy of the measurements will be sufficient to confirm or to refute modified theories of gravity. In this paper a new approach for the description of this formation is presented. The analytical solutions of the linearized motion equations are obtained. The distinctive feature of the solutions is that they use Cartesian coordinates of one of the spacecraft, termed the chief. These solutions have a clear physical meaning. It is shown, that the volume of a tetrahedron formed by spacecraft is a polynomial of 3-rd degree of Cartesian coordinates of the chief. The polynomial's coefficients are functions of initial spacecraft coordinates and velocities and linearly depend on time. If all spacecraft have the same periods of rotation around the Sun, the volume is a polynomial of 2-nd degree of the chief coordinates with time-independent coefficients. In this case the volume can be zeroed from 0 to 4 times per the period. Suggested approach can significantly simplify planning missions for measurements of various interplanetary fields.

The Hellings-Downs (HD) correlation serves as the fundamental benchmark for detecting the gravitational-wave background (GWB) in pulsar timing arrays (PTAs) within General Relativity (GR). However, this canonical signature relies on the idealization of a continuum of sources without interference. In realistic astrophysical scenarios dominated by supermassive black hole binaries (SMBHBs), interference between discrete sources induces intrinsic deviations in the spatial correlation, which may mimic or obscure signatures of modified gravity. In this work, we derive the closed-form spatial correlation functions for a GWB with arbitrary polarization and generic GW dispersion relations, in the presence of source interference. Through a rigorous harmonic analysis, we demonstrate that source interference modifies the correlation shape but strictly preserves the lowest non-vanishing multipole moment characteristic of each polarization, specifically the quadrupole for tensor, dipole for vector, and monopole for scalar modes. The truncation at higher-order multipoles is governed by the interplay between pulsar distances and dispersion effects. Furthermore, we quantify the statistical degeneracy between interference-induced variation and modified gravity signatures. We conclude that access to only a single realization of the Universe imposes a fundamental theoretical limit on distinguishing modified gravity from GR using spatial correlations alone.

Data-driven approaches like deep learning are rapidly advancing planetary science, particularly in Mars exploration. Despite recent progress, most existing benchmarks remain confined to closed-set supervised visual tasks and do not support text-guided retrieval for geospatial discovery. We introduce MarsRetrieval, a retrieval benchmark for evaluating vision-language models for Martian geospatial discovery. MarsRetrieval includes three tasks: (1) paired image-text retrieval, (2) landform retrieval, and (3) global geo-localization, covering multiple spatial scales and diverse geomorphic origins. We propose a unified retrieval-centric protocol to benchmark multimodal embedding architectures, including contrastive dual-tower encoders and generative vision-language models. Our evaluation shows MarsRetrieval is challenging: even strong foundation models often fail to capture domain-specific geomorphic distinctions. We further show that domain-specific fine-tuning is critical for generalizable geospatial discovery in planetary settings. Our code is available at this https URL

Spectral and spatial fringes in polarized light are produced by the interference of transmitted and reflected waves at the interface between materials with different indexes of refraction. These instrumental artifacts can affect the accuracy of optical designs conceived for high-sensitivity spectroscopy and polarimetry. We review the fundamental mechanisms that are responsible for these artifacts and the possible design pathways that allow us to mitigate them. In order to do so, we also present an approximate treatment of the problem of the transmission and reflection of light through (possibly absorptive) birefringent layers, relying on a few fundamental results that can be found in the already extensive literature on the subject. Unfortunately, many of these results remain the domain of a niche of investigators working in the field of thin films and optical coatings, and are often overlooked even by experienced designers of spectro-polarimetric instrumentation. The treatment presented in this work is limited to isotropic materials and uniaxial crystals, which are the most common types of optics employed in polarimetric instrumentation, and it fundamentally relies on the approximation of small birefringence for its implementation. An extensive set of modeling examples is provided to highlight the salient characteristics of polarization fringes, as well as to assess how approximate treatments such as this compare to exact but more computational expensive formulations of the problem such as Berreman's calculus.

The dense matter equation of state is of great interest due to the recent development of astrophysical observations for neutron stars. A rapid increase in pressure indicates a continuous crossover from a hadron phase to a quark phase without any phase transitions, yet its microscopic mechanism remains elusive. Recently, a peak in the speed of sound and a baryon momentum-shell structure, which are predicted from a quarkyonic matter picture, have been regarded as key features of the hadron-quark crossover. In this work, we explore a field-theoretical framework to describe the hadron-quark crossover, drawing an analogy with the Bose-Einstein condensate to Bardeen-Cooper-Schrieffer (BEC-BCS) crossover established in ultracold atomic experiments. Strikingly, a peak in the speed of sound and the baryon momentum-shell structure can simultaneously be explained by the tripling fluctuation effect arising from a different context of quantum many-body physics. We demonstrate these properties in a simplified model and provide a microscopic derivation of the quarkyonic matter model within our field-theoretical framework.

We investigate the shadow cast by a rotating traversable wormhole in the Teo class endowed with a general redshift function, with particular emphasis on the emergence of cuspy structures. The shadow boundary is the common envelope of two critical orbit families: unstable circular orbits outside the throat and orbits at the throat itself. The formation of cusps, marking the transition between smooth and cuspy shadow boundaries, only becomes possible when the redshift parameter $\lambda$ is allowed to vary. Moreover, we uncover a universal critical value $\lambda_c$ that signals the onset of the cusp. A phase diagram characterized by the spin and redshift parameters reveals four distinct morphologies: smooth, cuspy, ears touching, and throat drowning. The morphology of the wormhole shadow may provide observational diagnostics for the different compact objects in future high-resolution imaging observations.

After distinguishing the role of classical and quantum inhomogeneities in cosmological backgrounds, we constrain the initial states of the relic gravitons as soon as the different wavelengths of the spectrum cross the comoving Hubble radius, without any reference to earlier timescales. According to this pragmatic perspective the quantum states with finite energy density at the crossing time consistently affect the two-point functions and the related power spectra. An initial state different from the vacuum turns out to be marginally permitted in the low frequency range (associated with the largest observable wavelengths that crossed the comoving Hubble radius) while the intermediate and high frequency domains of the spectrum are populated by the gravitons produced from the vacuum. The non classical correlations are expected to dominate between the kHz and the THz since in this region gravitons are produced quantum mechanically with a negligible contribution from the initial state.

We assess the potential for measuring neutrino mixing parameters at energies above 1~TeV, for the first time, using the flavor composition of TeV--PeV astrophysical neutrinos, i.e., the proportion of $\nu_e$, $\nu_\mu$, and $\nu_\tau$. Today, flavor measurements inferred from the 11.4-year IceCube Medium Energy Starting Events sample are insufficient to constrain the mixing parameters due to limited statistics, challenges in flavor identification, and uncertainty in neutrino production. Yet, upcoming multi-neutrino-telescope observations -- even using only existing telescopes -- may achieve sensitivity to $\theta_{23}$ and $\theta_{13}$ when combined with traditional oscillation experiments. We establish the current status and future prospects for testing the three-flavor mixing framework in the previously unexplored TeV--PeV regime and quantify the minimum detectable size of flavor-modifying beyond-Standard-Model effects, providing a roadmap for high-energy neutrino mixing measurements.

This paper proposes a novel quantum effect wherein particle spectra show extra broadening due to local short-time acceleration (as well as in a deSitter spacetime background). This effect provides a simple interpretation for the acceleration interpolation relation required to modify the kinematics of a test particle in the Modified Inertial interpretation of Modified Newtonian Dynamics (MOND). This effect can be regarded as a generalization of the thermal blackbody spectrum generated by the Unruh effect (which arises from long-time uniform acceleration in a flat background) to the scenario of local short-time non-uniform acceleration. This effect offers a unified framework for understanding the accelerated expansion of the universe and the anomalies in galactic rotation curves or radial acceleration. It is worth emphasizing that this modified kinematic interpretation of MOND necessitates a quantum equivalence principle as its physical foundation, that is, extending the equivalence at the level of expectation values (first-order moments) in the classical equivalence principle to the equivalence at the level of second-order moment quantum fluctuations. Therefore, the role of this effect within the frameworks of quantum reference frames and quantum gravity is also discussed to some extent.

Posterior inference from pulsar observations in the form of light curves is commonly performed using Markov chain Monte Carlo methods, which are accurate but computationally expensive. We introduce a framework that accelerates posterior inference while maintaining accuracy by combining learned latent representations with local simulator-guided optimization. A masked U-Net is first pretrained to reconstruct complete light curves from partial observations and to produce informative latent embeddings. Given a query light curve, we identify similar simulated light curves from the simulation bank by measuring similarity in the learned embedding space produced by pretrained U-Net encoder, yielding an initial empirical approximation to the posterior over parameters. This initialization is then refined using a local optimization procedure using hill-climbing updates, guided by a forward simulator, progressively shifting the empirical posterior toward higher-likelihood parameter regions. Experiments on the observed light curve of PSR J0030+0451, captured by NASA's Neutron Star Interior Composition Explorer (NICER), show that our method closely matches posterior estimates obtained using traditional MCMC methods while achieving 120 times reduction in inference time (from 24 hours to 12 minutes), demonstrating the effectiveness of learned representations and simulator-guided optimization for accelerated posterior inference.

Gravitational-wave observations of binary neutron star systems can shed light on the currently unknown dense matter equation of state. The equation of state determines a large number of neutron star properties, such as tidal deformability, radius, and quadrupole moment, several of which directly affect the emitted gravitational-wave signals. To reduce the dimensionality when computing gravitational-waves and when interpreting observational data, quasi-universal relations are commonly employed to connect different neutron star properties. However, quasi-universal relations are not exact and their use may introduce uncertainty and bias. We explore the potential biases arising from different quasi-universal relations in the third generation era: (i) the Love-Q relation connecting the spin-induced quadrupole moment and the tidal deformability, (ii) the relation between the fundamental mode frequency and the tidal deformability, and (iii) the binary Love relation. We find that for the quadrupole relation biases are only present for rapidly rotating systems, for the binary-Love relation induces moderate biases only in the next-to-leading-order tidal parameters, which can however propagate into the inferred equation of state at low masses. Regarding fundamental mode frequencies, we find that the employed relation introduces only negligible biases, while waveform systematic effects can become comparatively large. Our results highlight that while quasi-universal relations remain a useful tool within gravitational-wave analyses, careful treatment is needed to avoid biases in equation of state measurements with next-generation detectors.

The excavation of the Hyper-Kamiokande cavern, 600 m underground, is complete. Measuring 69 m in diameter and 94 m in height, it is among the world's largest rock caverns. A vertically oriented, dome-capped cylindrical design was chosen to optimize cost and performance. Combined with substantial overburden, the geometry posed major engineering challenges. This paper outlines the underground works, main cavern design, excavation plan, and the evolution of support design and construction methods during excavation, namely the information-based (observational) design and construction approach.

Searches for stochastic gravitational wave backgrounds from first-order phase transitions offer a powerful probe of hidden sectors, but quantitative predictions in gauge theories are obstructed by the gauge dependence of the finite-temperature effective potential and the associated tunneling action. We study a minimal gauged $U(1)$ dark sector containing a dark Higgs and a dark photon, optionally supplemented by a vectorlike dark fermion, coupled to the Standard Model through Higgs portal or kinetic mixing. Using the Nielsen identity together with a controlled derivative expansion and power counting, we construct a gauge-independent effective action in the high- and low-temperature limits, enabling model-intrinsic nucleation dynamics and robust gravitational wave predictions. We perform dedicated Monte Carlo scans in both limits and map viable microscopic parameters to detector-facing peak frequencies and amplitudes, spanning bands relevant to pulsar timing arrays and planned space-based interferometers. In our scans, supercooled transitions typically produce the strongest signals, whereas parametrically high-temperature transitions are comparatively rare and tend to be weak. We further connect the phase transition phenomenology to viable dark matter candidates within the same minimal field content, providing benchmark targets for dark photon dark matter and dark fermion dark matter, and highlighting their complementarity with gravitational wave observables. Overall, our results provide an end-to-end, gauge-independent pipeline from a minimal hidden sector Lagrangian to gravitational wave spectra and cosmologically viable dark matter benchmarks, yielding the most reliable and concrete predictions to date for a minimal gauged $U(1)$ dark sector.