Papers for Friday, May 08 2026

Some high-mass stars likely end their lives in underluminous implosions that leave behind a black hole, known as failed supernovae (FSNe). However, neutrinos radiated during proto-neutron star formation generate a weak (Mach $\gtrsim 1$) shockwave in the outer layers of the star, which produces a unique transient as it breaks out of the dying star and signals its imminent disappearance. It was recently shown that there are two self-similar solutions that describe the propagation of this weak shockwave, and these solutions simultaneously contain outward-moving ejecta and fallback accretion onto the black hole. Here we show that the larger Mach number solutions are unstable, such that the Mach number of the shock grows with time $t$ and deviates from the self-similar prediction as $\propto t^{\alpha}$, with $\alpha \lesssim 0.1$, whereas the smaller Mach number solutions are stable. We also show that, above a critical mass loss that is readily achievable in core-collapse supernovae, the shock asymptotically strengthens and approaches the strong limit. Our results imply that it is the mass lost to neutrinos \textit{relative} to the mass enclosed by the shockwave, as well as the stellar density gradient where the shock forms, that primarily dictate its strength and the amount of material it ejects. These criteria explain why red supergiants, which have relative mass losses well in excess of the critical value at the time of shock formation, more readily eject material and create more luminous explosions compared to more compact progenitors.

The majority of discovered exoplanets have been observed orbiting field stars as opposed to within a star cluster. To determine whether the lack of observed exoplanets in star clusters is due to gravitational perturbations or observational limitations, we consider the possibility of studying exoplanets in stellar streams. We present the results of direct $N$-body simulations of planetary systems around stars that orbit within a star cluster. Our simulations demonstrate that stars with early cluster escape times tend to retain all their planets as they spend most of their time orbiting in the cluster's low-density outskirts. Alternatively, stars with later escape times can have a wide range of survival fractions as they are subjected to a range of local densities and encounter types. With respect to the stellar stream that forms as the result of the cluster's dissolution, stars near the edge of the stream are therefore more likely to have unperturbed planetary systems. Conversely, stars near the centre of the stream have a higher chance of having planets pushed to eccentric orbits, inclined orbits, or stripped from the system entirely. From our suite of simulations, we provide an estimate of the probability that a star will host a planet with a given initial semi-major axis $a_0$ based on the star's location along a stellar stream $\Delta \phi$.

We analyze the spectra of the central nuclei of extremely red quasars (ERQs) observed as part of the JWST ERS Q3D program. We focus on the complex kinematic structures of the prominent rest-frame optical emission lines. Our modeling allows us to deblend the lines and separate the emission into distinct kinematic components that imply velocity- and density-stratified gas structures on a range of physical scales within the ERQs. Supplementing the JWST data with archival data, we analyze the spectral energy distributions (SEDs) of the ERQs and find they are consistent with a significantly dust-obscured central source with a small amount of relatively-unobscured UV/optical flux that is scattered into our line-of-sight. While the kinematics of the UV and optical emission lines largely agree, the UV lines are dominated by scattered light. In contrast, the optical emission-line ratios indicate a combination of scattered and obscured emission. Our analysis focuses on one ERQ, J0834, because its distinct spectroscopic features allow the emission to be easily decomposed into separate kinematic components.

The degree of black-hole spin-orbit misalignment ("tilts") in the astrophysical population could be a powerful diagnostic to distinguish between binary formation in isolation, in dynamical environments, or in hierarchical triples. However, robust population-level spin tilt measurements are complicated by model misspecification as well as numerical and Poisson variance, ultimately owing to poor single-event constraints on tilts. Motivated by reports of a possible peak in the spin tilt distribution, we analyze the fourth LIGO-Virgo-KAGRA gravitational-wave transient catalog to test for preferred spin orientations at different black hole masses. We find that a peak in spin tilts is not statistically significant nor model independent. Since the data cannot be used to reliably identify subpopulations based on their spin tilt properties, we also consider a complementary approach: measuring the spin magnitude and tilt distributions at fixed mass scales. We find no confident correlation between mass and spin tilt, but we do confirm a confident correlation between spin magnitude and mass, corroborating recent analyses.

Active galactic nuclei (AGN) showing dramatic spectral and flux variations, either due to changes in the accretion rate (changing-state, CS-AGN) of the supermassive black hole or in the line-of-sight column density (changing-obscuration, CO-AGN), have been classified as changing-look (CL) AGN. Here we present a peculiar source, 2MASX J0413-0050, first identified as a narrow-line Seyfert 1 (NLS1s) galaxy in 2004. When re-observed twice in 2021, it showed a transition in the spectral type (towards a Seyfert 1.9) and the complete and mysterious disappearance of the Hbeta line while source was in a high accretion state. In the meantime, the X-ray flux decreased between observations taken in 2020 and 2022, and again in the most recent spectrum of 2023. Shortly after this, another optical spectrum revealed the re-emergence of both the narrow and broad Hbeta components (Seyfert 1.8). Despite the fact that it was not possible to retrieve the line-of-sight column density from the X-ray spectra, which would have helped in assessing whether this event could be attributed to a CO AGN scenario, the observational evidence does not necessarily support such an interpretation. J0413-0050 may have undergone several switch-on switch-off phases over the past 20 years, on an unknown timescale, which could have affected the accretion power and, consequently, the optical continuum and so the emission lines coming from the broad-line region (BLR). For these reasons, it is reasonable to classify this source as a CS-AGN. The case of J0413-0050 supports the hypothesis that NLS1s can indeed experience CL phenomena.

Molecular clouds form stars from the interstellar medium via gravitational collapse, following a sequence from low-density gas to high-density cores and eventually the formation of stars. In classical density wave theory, gas clouds orbiting the galaxy experience gas compression and triggered star formation, while encountering the gravitational well of spiral arms. We aim to trace these different phases of the molecular cloud life cycle via tracers of molecular gas (CO), dense molecular gas (HCN), and star formation (H$\alpha$, 24 $\mu$m) within the spiral arms of two grand-design spiral galaxies: NGC 4321 and M51 (NGC 5194). In the spiral arms of these galaxies, we investigate the relation between molecular gas, dense gas, and star formation (CO-HCN-SFR) at matched physical resolutions of 270 pc and 125 pc in NGC 4321 and M51, respectively. We employed spiral arm masks for these galaxies and investigate trends of HCN/CO and SFR/HCN (SFR/CO), which serve as proxies for the dense gas fraction and dense (molecular) gas star formation efficiency, perpendicular to the spiral arm spines. We find that HCN/CO, SFR/CO, and SFR/HCN increase from the upstream towards the downstream side of both spiral arms of NGC 4321, while their trends are less prominent in M51. Our results indicate that large-scale galactic dynamics (e.g. density waves) can induce a sequence of gas density and star formation-to-gas density variations perpendicular to the spiral arms. This sequence contributes to the increased scatter seen among spectroscopic ratios such as HCN/CO and SFR/HCN at sub-kiloparsec scales.

We report the discovery and multi-wavelength characterization of the Galactic transient EP J174942.2$-$384834, first detected by the Einstein Probe during a faint X-ray outburst in March 2025. Coordinated follow-up observations revealed two major outbursts and a rebrightening over a seven-month period. Broadband X-ray spectral modeling shows that the outburst emission was dominated by thermal Comptonization of very soft seed photons. The absence of a detected thermal disk component, together with the low inferred seed-photon temperature, is consistent with a cool and possibly truncated accretion disk. The X-ray spectrum remained consistently hard throughout the outburst activity, with a power-law photon index of $\Gamma \approx 1$-2, gradually softening as the flux declined. The optical/UV counterpart brightened in tandem with the X-ray emission and exhibited a blue continuum with broad Balmer absorption features. Together with the optical/UV - X-ray luminosity correlation, this supports a disk-dominated origin of the optical/UV outburst emission, with viscous heating likely playing a major role and irradiation possibly contributing, especially in the UV. No radio counterpart was detected, implying at most very faint jet activity. Taken together, the observed properties support the classification of EP J174942.2$-$384834 as a very faint X-ray transient black hole candidate. This study demonstrates the ability of Einstein Probe to uncover and characterize the faintest accreting compact objects in the Galaxy.

In the coming decade, the millihertz gravitational wave observatory LISA will provide the best constraints yet on the tens of thousands of close white dwarf binaries in the Milky Way, yielding unprecedented insights into the most abundant class of compact object binaries. Following inspiral via gravitational wave emission, interacting white dwarf binary pairs can lead to a multitude of outcomes, including AM Canum Venaticorum (AM CVn) binaries, R Coronae Borealis stars, young, rapidly-spinning single white dwarfs, (millisecond) magnetars, and a variety of explosive transients, most notably Type Ia supernovae. Current and future electromagnetic observations of these various outcomes coupled with the forthcoming flood of data from LISA place us on the precipice of a significant advance in our understanding of the long-term fate of white dwarf binaries. In this paper, we present a suite of mock catalogs of the Milky Way's white dwarf merger history, created using the population synthesis code $\texttt{COSMIC}$. We summarize the various merger outcomes expected (based upon varying white dwarf masses and chemical compositions) and explore ways the rates of these outcomes may vary with model uncertainties pertaining to binary evolution. We publicly release these merger catalogs as a tool for facilitating connections between gravitational wave science and white dwarf binary astrophysics.

NANOGrav and various pulsar timing array experiments recently reported evidence for a common red noise signal across millisecond pulsars. This signal exhibits Hellings-Downs inter-pulsar correlation patterns, providing compelling evidence for a stochastic gravitational wave background (SGWB) signal. In general, such a background can come from several astrophysical and cosmological phenomena. Assuming such SGWB has an inflationary origin, we use latest NANOGrav 15-year dataset to constrain the inflationary parameters e.g., tensor spectral index ($n_t$), tensor-to-scalar ratio ($r$), and explore the implications for the reheating phase through constraints on the reheating equation of state ($\omega_{\text{re}}$) and reheating temperature ($T_{\text{re}})$. We find the preference for an extremely blue-tilted tensor spectrum $n_t=2.20^{+0.36}_{-1.2}$ and the radiation-like reheating scenario $\omega_{\text{re}}=0.33^{+0.14}_{-0.36}$. Despite having no concrete evidence for the nature of the primordial vacua, the computation of gravitational wave (GW) sourced by tensor perturbations assumes the inflationary vacuum to be a Bunch-Davies vacuum. In this work, we examine modifications to the GW spectrum originating from the non-Bunch-Davies primordial vacuum. We find that NANOGrav observations favour a specific type of non-Bunch-Davies vacuum, known as the alpha-vacuum. Furthermore, our analysis demonstrates that the observations strikingly narrow down the range of the parameter $\alpha$ characterizing the vacua. On top of that, we find that a frequency-dependent parametrization of the vacuum parameter $\alpha$ beyond a threshold frequency can yield a minimal solution to alleviate the blue-titled issue. Finally, we highlight the possibility of testing such frequency dependence of $\alpha$ by probing the GW spectrum through future GW experiments.

Stellar rotation is a fundamental ingredient in shaping the evolution of stars and it can also be used to trace past stellar interactions. Yet, systematic studies of stellar rotation in large samples of stars belonging to different populations have only recently been made possible, thanks to spectroscopic surveys. We profit from the catalogue of rotational broadening and rotation periods released with Gaia DR3. We focus on open clusters to study the rotational behaviour of several interesting populations including, among others, blue stragglers and extended main sequence turnoffs (eMSTO). We use literature lists of almost a million member stars in several thousand open clusters in the Milky Way. We collect properties of stars and clusters from large surveys, including Gaia, and from various literature sources. We include a comprehensive collection of known variables and binary stars from various databases. We manually select (exotic) stellar populations from the color-magnitude diagrams of individual clusters and study their rotational properties. Our catalogue contains more than 44 000 rotationally characterised stars, almost 57 000 variables (excluding binaries) and more than 22 000 binary stars. We find several interesting results, including a few hundred new blue stragglers, several fast rotating red giants, and we increase the number of clusters with an eMSTO to 96. We discover that most clusters more massive than $10^3$ $M_{\odot}$ display an eMSTO. We present a new parametrization of the number of blue stragglers as a function of cluster mass and age. We find that the percentage of binary stars in the equal-mass binary sequence and in the main sequence are similar. We present the first large-scale statistical exploration of stellar rotation in open clusters, which already yielded new interesting results and which can be used as the basis for several detailed follow-up studies.

Strong gravitationally lensed supernovae (LSNe) are rare but extremely valuable probes of cosmology and astrophysics. Prompt identification within the alert streams of time-domain surveys such as the Rubin Legacy Survey of Space and Time (LSST) is essential for timely follow-up observations. In our previous study, Bag et al. (2026), we introduced a deep-learning framework for detecting LSNe Ia directly from multi-band, multi-epoch image cutouts. The model employs a convolutional LSTM architecture to capture spatiotemporal correlations in time-series imaging data, enabling classification updates as new observations arrive. In this work, we extend that framework by incorporating greater realism into the simulations. In particular, we present a method to construct realistic image time series from single-epoch observations by introducing epoch-to-epoch point spread function variations with corresponding variance-map corrections. The dataset is based on HSC PDR3 observations and includes simulated lensed host-galaxy arcs, SN light-curve variations, and Poisson noise. We also introduce an additional negative class consisting of SN Ia occurring in the foreground lens galaxy, representing a challenging source of false positives. Despite these additional complexities, the model retains strong performance. The receiver operating characteristic improves rapidly during the first few observations, reaching a true-positive rate of $\sim60\%$ at a false-positive rate of $\mathcal{O}(10^{-4})$ by the seventh observation and $\sim80\%$ by the tenth. We also investigate potential confusion with sibling SNe occurring in LRGs and identify the configurations that best mimic lensed systems. These results demonstrate that the image-time-series approach remains robust under more realistic observing conditions, and is well suited for real-time LSN searches in LSST and other time-domain surveys.

We present the results for the dynamics and emission profiles of axi-symmetric numerical simulations of structured gamma-ray burst afterglow jets, computed using the relativistic moving-mesh hydrodynamics code GAMMA. We find that the spreading of jets of average opening angle is moderately impacted by the initial steepness of the angular structure, although the effect disappears once the working surface of the jet substantially exceeds its initial width, and that the travel time of a sound wave across the front surface remains the best indicator of the onset of spreading also for structured jets. When computing the afterglow spectrum using a local cooling approach that traces the electron population following shock-acceleration, we observe a significant impact on the synchrotron cooling break. Similar to earlier results for top-hat jets, the cooling break is found to shift upward in frequency by well over a factor of ten relative to approaches that assume a global cooling timescale across the jet. The cooling break transition in the spectrum also becomes substantially smoother. For both local and global cooling, jet breaks become sharper with increasing frequency. Local cooling is found to initially lead to a steeper slope post jet-break. The local-cooling emission is shown to originate from a narrow frequency-dependent sized region behind the shock front, as expected, but in strong contrast to a global cooling approach.

Primordial magnetic fields (PMFs) can enhance the abundance of low-mass halos during Cosmic Dawn by sourcing additional small-scale matter fluctuations. This enhanced small-scale power can accelerate early galaxy formation, shifting the timing of Lyman-$\alpha$ coupling, X-ray heating, and reionization toward earlier times and imprinting correlated signatures on the global and fluctuating 21-cm signals. We extend the fast analytic framework {\tt\string zeus21} to include a physically motivated PMF contribution to the linear matter power spectrum, including radiative damping before recombination and magnetic-pressure suppression below the magnetic Jeans scale. The implementation preserves the speed and modularity of {\tt\string zeus21}, enabling efficient exploration of PMF parameter space. For $n_B=-2.9$, we quantify the impact of PMFs on early structure formation and 21-cm observables across a range of fiducial magnetic amplitudes, and forecast detectability with \textit{HERA} and \textit{SKA}. Combining 21-cm forecasts with external CMB priors, we find that upcoming experiments can probe PMFs through their impact on small-scale structure, providing constraints complementary to existing cosmological probes.

We present a simple and interpretable representation of TESS light curves designed for large-scale exploratory analysis. Our goal is not to optimize classification performance, but to construct a computationally efficient mapping in which proximity reflects meaningful similarity, without using labels or explicit period information as inputs. We represent each light curve using either quantile graphs or scattering transforms, reduce dimensionality with principal component analysis, and project the resulting features onto a self-organizing map (SOM). We evaluate ~1500 model configurations using a combination of standard embedding diagnostics and a light-curve-shape-based cohesion metric, and select a compact quantile-graph-based model that balances interpretability, stability, and performance. Applying the model to ~1.5 million TESS 2-minute cadence light curves, we find that the map organizes sources primarily by variability amplitude, signal-to-noise ratio, characteristic timescale, and light-curve shape. Repeat observations of the same stars show that most sources occupy stable and contiguous regions of the map, indicating that the representation captures persistent properties rather than noise and systematics. We provide an interactive web interface at this http URL that enables inspection of nodes, nearest neighbors, and individual sources across sectors. The resulting representation serves as a practical tool for exploration, anomaly detection, and dataset characterization, and illustrates how simple, deterministic encodings can yield useful structure in large astronomical time-series datasets.

We investigate the mass-metallicity relation (MZR) from z=1 to z=9 using electron-temperature-based gas-phase metallicities and examine how auroral-line selection, star-formation history, and secondary abundances affect its interpretation in the early Universe. We compile a homogeneous sample of 286 star-forming galaxies observed with JWST/NIRSpec medium resolution spectroscopy, selected through detections of the [O\,III]\,$\lambda$4363 auroral line from the public DAWN JWST Archive (DJA). We derive electron densities, temperatures, and oxygen abundances using the direct $T_e$ method, along with relative N/O and He/H abundances. Stellar masses are obtained via SED fitting and star-formation rates from reddening-corrected Balmer emission lines. To quantify auroral-line selection biases, we additionally stack galaxy spectra with and without auroral-line detections, extending the MZR into regimes inaccessible to individual measurements. The auroral-line-detected sample spans log(M*/Msun)=6.77-10.5 and 12+log(O/H)=6.9-8.4. A linear fit gives an MZR slope of $\gamma$=0.38 $\pm$ 0.09. Stacked galaxies without individual $\lambda$4363 detections define a relation with a similar slope but metallicities higher by ~0.2-0.3 dex at fixed stellar mass. Auroral-line detections also show higher SFRs, larger equivalent widths, and larger offsets from the fundamental metallicity relation, whereas non-detections appear more chemically evolved and closer to it. Several stacked bins also show enhanced N/O and He/H ratios. These results indicate that the low-mass high-redshift MZR traced by JWST is shaped by both recent star-formation history and auroral-line selection effects. Auroral lines preferentially identify high-EW, high-sSFR galaxies in the low-metallicity envelope, whereas non-detections reveal a more enriched sequence closer to the metallicity expected from the FMR.

Radio Frequency Interference (RFI) increasingly contaminates the radio astronomy spectrum, often exceeding astronomical signal amplitudes by 50-70 dB. Reliable detection and mitigation are therefore essential for studies of faint transient phenomena such as pulsars and fast radio bursts (FRBs). Existing practical methods (including Spectral Kurtosis (SK), Median Absolute Deviation (MAD), and SumThreshold) perform well in many settings but depend on assumptions about the RFI environment and data statistics, limiting their effectiveness for weak, broadband, or non stationary interference. We develop a transform based RFI detection method that requires no prior knowledge of RFI origin or type. Using Green Bank Telescope (GBT) data containing PSR J1713+0747, with 4096 channels spanning 1.1-1.9 GHz and 5.12 micro second sampling, we apply a Short Time Fourier Transform (STFT) to each channel and use an image segmentation algorithm on the STFT magnitude to generate a binary RFI mask. The masked data are inverse transformed and reassembled into a cleaned time series. Performance is assessed using the Signal to Noise Ratio (S/N) of a single pulse of PSR J1713+0747, with SK serving as the baseline. The cleaned spectrogram is dedispersed, integrated across frequency, and evaluated through the resulting S/N. Experimental results show that refining each channel's frequency content via STFT, followed by segmentation in the STFT domain, yields measurable improvements in RFI suppression.

X-ray polarimetry is now providing a new way to look at the high energy sky. The addition of two observables, polarization fraction and angle, reveals crucial new information on the structure of accretion flows and magnetic fields in astrophysical systems. Here, we review the basic physical processes that produce polarized X-rays in astrophysical contexts. Then, we briefly describe the physical processes used to measure X-ray polarization and the detectors that have been flown or are under construction.

Silver is an important light neutron-capture element whose stellar abundances help constrain the origin of the weak r-process. The Sun is an important reference point for such studies; moreover, being a moderately volatile element in CI chondrites, the solar silver abundance is interesting as a diagnostic for the debated Sun-CI abundance vs. condensation temperature trend. These studies require accurate silver abundances that go beyond the commonly used assumptions of 1D atmospheres and local thermodynamic equilibrium (LTE); however, no consistent 3D non-LTE analysis of silver has been available to date. We present a new Ag I model atom built from carefully curated radiative and collisional data, including newly computed oscillator strengths using an ab initio multi-configurational Hartree-Fock method and inelastic hydrogen collision rates based on a combined asymptotic and free-electron model approach. We assess modelling uncertainties via targeted sensitivity tests, finding the results most sensitive to hydrogen collision data. Applying the model to the solar Ag I 328 and 338 nm resonance lines, we find severe positive abundance corrections from coupled 3D and non-LTE effects. Using revised equivalent width measurements, we derive a recommended solar 3D non-LTE silver abundance of 1.15 +/- 0.08. This is an increase of 0.19 dex relative to the current reference value. Our ab initio model significantly reduces the discrepancy with the meteoritic value from 0.25 to 0.06 dex; moreover, this residual offset is consistent with recent results for other moderately volatile elements. The Sun provides the benchmark test for the first Ag I non-LTE model atom presented here. In subsequent work, this model will be applied to determine 3D non-LTE silver abundances in metal-poor dwarfs and giants, enabling improved constraints on Galactic chemical evolution and weak r-process nucleosynthesis.

We present high-resolution zoom-in simulations of molecular clouds exposed to an interstellar radiation field and cosmic ray ionisation rate up to 1000 times stronger than that of the solar neighbourhood. We detail the evolution of the accretion discs that form around the first protostar in each simulation, for a total of 7 discs, for up to 100 kyr. The use of a zoom-in procedure allows for the au-scale discs to be well resolved (with resolution < 0.25 au) whilst retaining the structure of the wider parsec-scale molecular cloud. We find that discs exposed to a stronger radiation field tend to be more massive, hotter and denser. Similarly, their host stars grow to become more massive as a result of accreting more rapidly from their surroundings. All the discs show evidence of recurrent instability during the simulations, but only some of them fragment. We investigate whether stability metrics, such as the Toomre $Q$, $\alpha$ viscosity, and $\beta$ cooling parameter, can predict fragmentation by calculating them just before the discs fragment. We find that the metrics are generally unable to do so, as the discs appear stable even up to a few hundred years before fragmenting. In solar-like environments fragments are typically of planetary mass and often migrate to the centre of the disc, whereas fragments in a high-radiation environment are massive ($\rm > 0.1 \, M_\odot$) and fully disrupt/accrete from the progenitor disc. We conclude that the evolution and properties of circumstellar discs depend on both their radiation and physical environment.

Interacting supernovae probe the twilight years of massive stars, exhibiting signatures of interaction between the supernova ejecta and surrounding material expelled from the progenitor. We present the peculiar interacting supernova, SN\,2025ngs in NGC5961 (37.8 Mpc). This transient toes the line between strongly interacting supernovae (type IIn) and type IIP supernovae. SN 2025ngs presents photometrically as a short-plateau supernova, with a plateau duration, t$_{\mathrm{PT}}^{}\approx70$ days. Interaction features subside within a week post-explosion, consistent with the growing number of flash supernovae, giving way to a short period where a typical IIP spectrum is exhibited. Towards the drop off the plateau, interaction features re-emerge, exhibiting complex H$\alpha$ profiles throughout the rest of the transient evolution. We compare with models of early spectra, finding the abundances generally consistent with a supergiant progenitor with a high mass-loss rate (10$^{-3}$ M$_\odot$ yr$^{-1}$). Early, high-resolution spectra reveal a double-horned H$\alpha$ profile, providing strong evidence for shock interaction with a proximate disk-like circumstellar medium. Spectroscopically, SN 2025ngs closely resembles the luminous SN 1998S, despite photometric differences, with SN 2025ngs having a relatively modest peak magnitude of $M_\mathrm{V}=-17.9$ mag, adding another member to the surprisingly diverse 98S-like group.

We have been conducting a search for narrowband radio signals with the L-band receiver (1.15-1.73 GHz) of the 100 m diameter Green Bank Telescope (Margot et al., 2023). So far, we have captured radio emissions from 70,000+ stars and planetary systems in the ~9 arcminute beam of the telescope. Our data-processing pipeline has a demonstrated 94%-99% efficiency for the detection of narrowband signals across the full range of frequency drift rates (+/-9 Hz/s). All 100 million candidate signals detected to date were either automatically (99.5%) or visually (0.5%) confirmed to be anthropogenic in nature. These results allow us to place stringent limits on transmitter prevalence: at the 95% confidence level, fewer than 6.3e-5 of stars within 20,000 ly host a transmitter that is detectable in our search (EIRP > 5e16 W). Our most interesting signals have been uploaded to a citizen science platform (this http URL), where 40,000+ volunteers to date have contributed insights and classifications. We are using artificial intelligence (AI) to accelerate our search, automatically excise radio frequency interference, and improve signal detection. UCLA SETI research has involved ~200 undergraduate and ~20 graduate students so far.

We present sensitive ($\sim$9 $\mu$Jy), sub-arcsecond resolution radio continuum observations at 9.0 GHz (3.3 cm) and 14.0 GHz (2.1 cm) obtained with the Karl G. Jansky Very Large Array (VLA) toward the nearby Coronet Cluster in Corona Australis (d $\approx$ 150 pc). We monitored the region from March 2012 to February 2015 using all available VLA configurations, allowing us to construct deep X- and Ku-band maps at multiple angular resolutions. We detected 20 radio sources, including 14 previously known Young Stellar Objects (YSOs), five sources possibly associated with shock emission, and one background galaxy. We resolved IRS 5, previously known to be a binary system, and identified IRS 7A and IRS 7B as multiple systems at centimetre wavelengths. The younger Class 0 and I YSOs exhibit spectral indices $\alpha_{pk}$ ranging from -0.4 to 1.7, while the more evolved Class II YSOs show flatter values between 0 and 0.8, consistent with free-free emission, with minor contributions from non-thermal emission. The Class III source is only constrained by an upper limit. Radio variability, measured as a fraction of the mean intensity peak, is found to be ubiquitous and independent of evolutionary stage. Variability structure functions computed for nine sources indicate no preferred timescales for most of them. We also investigate spectral index variability for six sources and find significant variations in only one object. Finally, we analyse the extended radio emission toward IRS 7B, where some subcomponents exhibit negative spectral indices suggestive of non-thermal processes.

We present a multi-wavelength analysis of the nearby millisecond pulsar PSR J0437--4715, combining Hubble Space Telescope (HST) far-ultraviolet, ROSAT soft X-ray, and XMM-Newton X-ray data, to model its broadband emission and energy-resolved pulse profiles, and infer key stellar parameters via Bayesian inference. The broadband emission includes cold thermal, hot thermal, and non-thermal components: cold bulk surface emission is modeled with a non-magnetized partially-ionized hydrogen atmosphere; hot-spot emission adopts the pulse profile modeling technique with a non-magnetized fully-ionized hydrogen atmosphere model; and non-thermal emission is included as a phase-invariant power-law component. By adopting an informative prior on the hot-spot geometry informed by radio polarization position angle measurements, the joint multi-instrument analysis yields a statistically viable and radio-consistent solution with a gravitational mass of 1.38$\pm$0.03~M$_\odot$ and an equatorial circumferential radius of 13.25$_{-0.35}^{+0.34}$~km (68\% confidence intervals). The hot-spot geometry consists of two spherical caps with uniform temperature distributions: the primary hot spot is situated at a colatitude of $\approx$130$^{\circ}$, and the secondary hot spot lies at a colatitude of $\approx$9$^{\circ}$, close to the north pole. It yields tighter radius constraints than HST+ROSAT fits and shifts the radius posterior distribution to larger values relative to NICER-only fits. This work demonstrates the importance of multi-wavelength data in refining neutron star mass-radius measurements and resolving geometric degeneracies.

We report on a gravitational wave search for compact binary coalescences involving at least one component with mass between $0.2 M_\odot$ to $1 M_\odot$, and ratio of component masses between 0.1 and 1. The analysis uses data collected by the LIGO detectors between May 24 2023 15:00 UTC and January 16 2024 16:00 UTC. No statistically significant sub-solar mass candidates were identified by the participating search algorithms. We report the detection sensitivity of the current searches to the target sub-solar mass black hole population, while also reporting the sensitivity of the search to low-mass neutron star binaries for the first time. With the absence of detections, we place upper limits on the merger rate of sub-solar mass black holes, ranging from $110$ $\mathrm{Gpc^{-3}yr^{-1}}$ to $10000$ $\mathrm{Gpc^{-3}yr^{-1}}$ at 90% confidence. We use the merger rate limits to constrain two illustrative dark matter scenarios that can form sub-solar mass compact objects: primordial black holes, and dark black holes forming in a dissipative dark matter model. For late-forming primordial black hole binaries, our search excludes the fraction of dark matter in primordial black holes to be 1 for masses above $0.9 M_\odot$. In the early-formation scenario, we limit this fraction to be $\leq$ 7% at $1 M_\odot$, and $\leq$ 40% at $0.35 M_\odot$. For the dissipative model, the excluded region in the parameter space of dark matter fraction in dark black holes and their minimum possible mass extends down to (1.2 to 1.3) $\times10^{-5}$ when the minimum mass is $1 M_\odot$. For binary neutron stars that include sub-solar mass components, we estimate the sensitive space-time hypervolume to be $10^{-3}$ $\mathrm{Gpc^3yr}$, and report the upper limit on their merger rate for a simple, fixed population as ~86 $\mathrm{Gpc^{-3}yr^{-1}}$.

Three decades ago, humanity entered the Exoplanet Era, with the discovery of the first planets orbiting other stars. Today, more than 6000 exoplanets are known - a tally recently bolstered by NASA's TESS spacecraft. Whilst TESS is an exceptional planet finding machine, dedicated follow-up observations from the ground are required to confirm the existence of the planets it discovers. To achieve this, we constructed the southern hemisphere's only dedicated exoplanet detection and characterisation facility, MINERVA-Australis, at the University of Southern Queensland's Mt Kent Observatory. Funded in 2015, MINERVA-Australis saw first light in 2018, in time for the launch of TESS. MINERVA-Australis has since been scouring the skies, working to confirm and characterise the incredible harvest of planets detected by TESS. To date, the facility has contributed to the discovery of 40 new exoplanets, and continued the legacy of radial velocity data from the Anglo-Australian Planet Search program.

This work focuses on the development and demonstration of tunable superconducting on-chip resonator, leveraging the intrinsic current-dependent non-linear kinetic inductance of superconducting aluminium, and investigating the effect of oxygen content. Thin films are deposited using standard metal evaporation. We present results from a comprehensive study based on a series of evaporated Al thin films. This research aims to inform and constrain optimisation strategies for the design of mm-wave on-chip spectrometers, particularly regarding yield, resolution, and efficiency. By systematically varying film stoichiometry, we use a series of DC measurements to extract fundamental film properties such as resistivity, critical current and critical temperature. Furthermore, we employ low-loss DC-coupled microwave resonators to characterise both their microwave properties and the non-linear kinetic inductance, comparing these findings to a determined non-linear kinetic model. Finally, we discuss the possibility of usage in a parametric amplifier.

We present the first direct computation of spatially averaged dynamical quantities in the local Universe, employing the Cosmicflows-4++ reconstruction and a covariant scalar averaging formalism. We extract the domain-averaged density, expansion rate, spatial curvature, and kinematical backreaction over cosmologically relevant domains around our Galaxy, extending up to a comoving radius of $300~\mathrm{Mpc}/h$. The resulting domain-averaged present-day energy budget features nontrivial variations with scale that reflect a nested structure within the cosmic neighborhood, including a large-scale void shell encompassing the local cosmic web. Remarkably, we find significant contributions to this energy budget from the average spatial curvature at the $\mathcal{O}(10\%)$ level on all probed scales. By contrast, the kinematical backreaction remains much smaller throughout the surveyed volume, reaching at most a $\mathcal{O}(1\%)$ contribution on the smallest scales considered, i.e., $30~\mathrm{Mpc}/h$. Convergence to the global $\Lambda$CDM background is not observed within this range of scales.

TRAPPIST-1 is an ultra-cool dwarf that flares frequently. These flares shape the surrounding planets' high-energy irradiation environments, with consequences for atmospheric chemistry and escape, and they can contaminate transmission spectroscopy of those planets. A quantitative flare-frequency distribution (FFD) spanning the full energy range is therefore essential for both interpreting JWST spectra and modeling the planets' irradiation histories. Here we present a unified FFD over four orders of magnitude in energy by jointly analyzing $\approx$87\,hr of JWST/NIRISS and JWST/NIRSpec time-series spectroscopy together with $\approx$74\,days of \textit{Kepler}/K2 photometry. To enable a consistent comparison across these heterogeneous datasets, we convert all events to energies in the TESS bandpass. For the Kepler-to-TESS conversion we adopt a cooler flare continuum appropriate for ultra-cool dwarfs ($T_{\rm flare}=3500$\,K). After correcting for flare-detection sensitivities, the combined JWST+K2 cumulative FFD is consistent with a single power law, $N(\ge E_\mathrm{TESS})\propto E_\mathrm{TESS}^{-\beta}$, with $\beta=0.753$ over $E_{\rm TESS}\simeq10^{29}$-$10^{33}$\,erg. The slope of the distribution indicates that the time-averaged flare energy budget is dominated by rare, high-energy events rather than by the more numerous low-energy flares. This bandpass-consistent FFD provides a practical basis for JWST transit-spectroscopy planning and for modeling the flare-driven irradiation environment of the TRAPPIST-1 planets.

Context: The Circular Restricted Three-Body Problem provides a fundamental framework for understanding resonant dynamics in binary star systems. Aims: We develop a unified Hamiltonian formulation for mean-motion resonances that encompasses both circumstellar and circumbinary planetary orbits within the CR3BP. Unlike the Solar System case where the perturbing body is a planet of negligible mass, here the perturber (a stellar companion) has a non-negligible, finite mass, a crucial difference that we fully incorporate. Methods: Starting from the full Hamiltonian in each configuration, we perform canonical transformations to resonant action angle variables and derive reduced one-degree-of-freedom Hamiltonians through systematic averaging over the fast orbital motion. Leading-order scaling laws for the Fourier coefficients of the resonant perturbation are obtained, revealing their dependence on the binary mass ratio and the planet's orbital distance. Results: The resulting effective potential is shown to exhibit bistability under the well-defined condition |epsilon2/epsilon1| > 1/4, where epsilon1 and epsilon2 are the amplitudes of the first two resonant harmonics. This bistability creates the essential dynamical setting for stochastic resonance. Scaling laws for the Fourier coefficients are derived for both S-type and P-type configurations. Estimates for known binary-planet systems show that while currently observed systems lie below the bistability threshold, the theory predicts that extreme configurations (a/a_b <~ 1.5 for P-type, almost equal mass binary) could host bistable resonances accessible to future observations. Conclusions: This work provides a natural Hamiltonian framework for studying stochastic resonance in binary planetary systems, bridging analytical celestial mechanics and the nonlinear dynamics of exoplanetary systems subject to realistic perturbations.

We identify a systematic distortion of the gain-vs.-frequency function of radio telescopes caused by digital flattening ("whitening") of the signal's spectrum followed by re-quantization, a common pair of processes in the signal processing of modern telescopes. Wide-bandwidth telescopes often have a large variation of signal power over frequency. Flattening of the spectrum allows samples of the channelized signal to be represented in a small number of bits, allowing efficient downstream processing. However, we show that this produces subtle systematic error in the measured spectra. We explore this effect in data from the Owens Valley Radio Observatory's Long Wavelength Array (OVRO-LWA) and through detailed semi-analytic simulations. Although the effect can be small so that it has heretofore been unrecognized, we demonstrate that it produces distortion of the spectrum at a level that is problematic for some science, in particular 21 cm cosmology. Finally, we explore mitigation strategies, showing that the effect can be substantially reduced by careful choice of the gain distribution along the signal path or by incorporating dithering in the re-quantization step.

Recent observations of metal-poor, star-forming dwarf galaxies reveal He III regions, traced by nebular He II 4686 emission that require a strong source of extreme-ultraviolet (EUV) radiation. The origin of this hard ionizing radiation remains poorly understood, as standard stellar populations fail to account for it, posing key implications for the understanding of early galaxy formation. We present a systematic Chandra X-ray study of 21 nearby star-forming galaxies with He II emission but lacking Wolf-Rayet spectral signatures. Using 7 new and 36 archival Chandra X-ray observations combined with optical stellar population synthesis modelling, we constrain the ionizing continuum required to sustain the observed He II line, the ionizing continuum available from X-ray objects, and the properties of the host H II regions. We find that the inferred EUV output from accreting X-ray sources in our sample is systematically lower than what is required to produce the observed He II emission. Our sample is consistent with established empirical scaling relations for X-ray luminosity, indicating that this discrepancy cannot be attributed to an anomalously low number or luminosity of X-ray sources. These results indicate that accreting X-ray sources alone cannot account for the observed He II-ionizing photon budget, pointing to additional or alternative sources of hard EUV radiation in metal-poor star-forming environments. Potential alternative or additional contributors are discussed.

We investigate the population properties of binary black holes (BBHs) from the LIGO-Virgo-KAGRA collaboration, focusing especially on those in the high-mass range, using the newly released GWTC-4 catalog. For the first time, we search for a subpopulation of low-spin intermediate-mass black holes (IMBHs) that would indicate formation via stellar core collapse. With the currently available catalog, we find no evidence for such a subpopulation, and set a 90\% upper limit on the merger rate of collapse-formed IMBHs at $0.077~\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$. The mass distribution of low-spin (stellar-origin) black holes truncates at $65^{+23}_{-22}\,M_\odot$, consistent with the lower edge of the pair-instability mass gap (PIMG), although we cannot directly determine its upper boundary from current data. Informed by stellar evolution theory, we estimate the upper edge of the PIMG to be $150\pm24\,M_\odot$. We find that the observed IMBHs belong to a high-spin subpopulation, consistent with formation through successive hierarchical mergers.

We analyzed the distribution of spin parity in spiral galaxies using the HSC DR2 data. The spiral winding parity of disk galaxies, observed as S-spiral or Z-spiral projected onto the sky plane, provides robust information on the sign of the line-of-sight component of their spin vectors, specifically whether the spin vector points toward or away from us. The distribution of 49,494 S/Z annotated spirals with spectroscopic redshift (0.05 $\le z$) was analyzed for 46,247 fiducial cubic search volumes of various sizes, 20--200 Mpc, deployed in the 3D supergalactic coordinates. We counted the number of S-spirals and Z-spirals in each cube, evaluated the binomial probability of the observed S/Z imbalance, and identified statistically anomalous cube candidates. The observed cumulative distribution functions for the 256 sets of cubes are in good agreement with the theoretical binomial distribution and with those obtained from 1000 Monte Carlo realizations assuming random S/Z spin assignments. The number of statistically anomalous cubes is also comparable to that expected from the random assignments. These results indicate that the spin-vector distribution of spiral galaxies is consistent with statistical randomness expected from the standard cosmological model of structure formation.

Modern astronomical observatories generate a massive volume of multimodal data, creating a critical bottleneck for expert human review. While multimodal large language models (LLMs) have shown promise in interpreting complex visual and textual inputs, their ability to perform specialized scientific classification while providing interpretable reasoning remains understudied. We introduce AstroAlertBench, a comprehensive multimodal benchmark designed to evaluate LLM performance in astronomical event review along a three-stage logical chain: metadata grounding, scientific reasoning, and hierarchical classification over five categories. We use a pilot sample of 1,500 real-world alerts from the Zwicky Transient Facility (ZTF), a wide-field survey that scans the northern sky to detect transient astronomical events. On this dataset, we benchmark 13 frontier closed-source and open-weight LLMs that support visual input. Our results reveal that high accuracy does not always align with model ``honesty,'' defined as the ability to self-evaluate its reasoning, which affects its reliability as a real-world assistant. We further initialize a human-in-the-loop evaluation protocol as a precursor to future community-scale participation. Together, AstroAlertBench provides a framework for developing calibrated and interpretable astronomical assistants.

Dark matter may play an important role in galaxy formation through its non-trivial properties. For example, self-interacting dark matter may contribute to the formation of the widely observed core structures in galaxies. However, galaxy formation is a complex process, and such core structures can also arise from baryonic effects within the cold dark matter framework. To clarify the role of dark matter self-interactions, it is necessary to study systems that evolve without significant baryonic disturbances. Low-surface-brightness galaxies in the field, which are gravitationally isolated and have evolved with minimal external influence, are suitable candidates for this purpose. Since these galaxies typically contain only a small amount of baryonic matter, strong baryonic effects are not expected in their evolutionary history. In this study, we assume that these galaxies decoupled from proto-clusters at high redshift. Based on this assumption, we set initial conditions and estimate the time required for core formation, which we compare with the time corresponding to the redshift of proto-clusters. We examine five low-surface-brightness galaxies in the field and three observed proto-clusters at redshifts z=2.45, 7.66 and 7.88. Our analyses, based on order-of-magnitude estimates without numerical simulations, excludes a self-interaction cross section of sigma/m = 1 cm^2/g, while sigma/m = 0.1 cm^2/g is favored. This result is consistent with constraints derived from the shapes of present-day cluster cores.

The extremely metal-poor nature of the C-19 stream indicates that its progenitor was a primordial stellar system born in the very early Universe. Current observations show that it has a small metallicity dispersion (0.18 at the 95% confidence level), which is the signature of a globular cluster origin, while at the same time displaying an unusually large velocity dispersion ($\sim10$ km/s) typical of dwarf galaxies. To reconcile this conflicting observational evidence, previous simulations have focused on potential interactions with dark matter subhalos, which can efficiently make a cluster stream dynamically hot. In this work, we explore internal dynamical processes in star cluster formation, focusing on initial conditions shaped by gas expulsion and a top-heavy initial mass function. We find that the large observed velocity dispersion and broad stream morphology can be reproduced by a cluster that underwent severe gas expulsion and expansion during its birth phase, which is potentially a typical formation scenario of extremely metal-poor star clusters. A top-heavy IMF and binaries can also increase the velocity dispersion. The formation of C-19 may involve a combination of these effects.

We present spatially resolved Keck/LRIS spectroscopy of three clumpy star-forming galaxies at $z\sim0.1$, comparing outflow properties traced by H$\alpha$ and Mg II emission with those probed by Mg II and Na I D absorption. Outflow velocities measured using Mg II absorption ($\langle v_{\rm out} \rangle = -560 \pm 30$~\kms) are consistently higher than those traced by H$\alpha$ emission ($\langle v_{\rm out} \rangle = -124 \pm 3$~\kms) across $\sim$5 kpc$^{2}$ regions. Despite this offset, the correlation between $v_{\rm out}$ and galaxy properties, such as SFR and $\Sigma_{\rm SFR}$, show similar slopes for both tracers, with Mg II absorption systematically offset by $\sim 0.4$ dex. In two galaxies, Mg II emission is also detected, yielding velocities consistent with H$\alpha$. In one galaxy we also detect outflows in Na I D absorption and find similar velocities as Mg II in absorption, which leads to a $\sim$0.4 dex higher Na I D outflow velocities compared to those measured in emission. Our spatially resolved results are consistent with those found for galactic-scale measurements, implying the outflow relationships are similar from the sales of $\sim$1-2 kpc to global measurements. Combined with literature measurements, these results suggest that the offset in velocities is driven not by ionisation state, but rather by the systematics associated to how absorption and emission measures trace the gas density.

Balmer-dominated shocks in Type Ia supernova remnants offer powerful probes into collisionless shock physics and hints towards supernova progenitor environments. Prior studies focused on the hydrogen Balmer lines, which manifest as a superposition of broad and narrow emission lines. Using integral-field spectroscopy with MUSE, we discovered broad and narrow helium emission lines from Balmer-dominated filaments of three Type Ia supernovae remnants in the Large Magellanic Cloud: SNR 0509-67.5, SNR 0519-69.0 and N103B. We detect broad and narrow He~\textsc{i} 5876~Å~,7065~Å emission in SNR 0519 and N103B and He \textsc{ii} 8236~Å in SNR 0519. In SNR 0509 we detect narrow He~\textsc{i} 5015~Å, 6678~Å, 7065~Å and 7281~Å, with only 7065~Å~ exhibiting a broad component. The detection of narrow He\,\textsc{ii} challenges existing shock models, where such emission is not expected, and may indicate either incomplete ion-ion equilibration behind the shock or an origin in shock precursors. For SNR 0509 and N103B, the neutral He/H line ratios indicate enhanced helium abundances, whereas SNR 0519 is consistent with the primordial He/H value. We therefore propose helium emission in Balmer-dominated shocks as a new diagnostic of shock physics and Type Ia supernova circumstellar environments. Although our modeling is primarily a proof of concept, it demonstrates the possibility to infer the total He-to-H abundance ratio, with dominant uncertainties arising from the assumed initial ionization fractions. Despite the uncertainties, we demonstrate that narrow helium lines can serve as effective probes of circumstellar conditions and progenitor evolution when analysed alongside reliable constraints on the preshock neutral H/He abundance ratio.

Polar brightening (PB) observed at microwave frequencies serves as an important probe to study the thermal and magnetic properties in the Sun's polar regions. Building on earlier studies that linked microwave PB to polar faculae, small-scale loops, and the polar coronal holes (PCHs), we present a comprehensive analysis of the long-term behaviour of 17 GHz microwave PB and its relation to polar magnetic field and coronal hole evolution. Using daily Nobeyama Radioheliograph observations spanning 1992 to 2018, we quantify microwave PB peak temperature variations and compare them with the temporal evolution of PCH area extracted from SDO/AIA-based SPoCA coronal hole catalogues during the period 2010-2018. We also examine the correspondence between microwave PB and the polar magnetic field to assess the nature of their association. Our results show a strong correlation between microwave PB peak temperature and PCH area, as well as with the polar magnetic-field strength. In addition, we found that regions of enhanced microwave emission are frequently associated with small-scale loop structures, consistent with Coronal Bright Points (CBPs), which are often associated with the eruption of jets. Overall, this study aims to investigate the impact of coronal holes, polar magnetic fields, and small-scale polar activity on polar brightening observed at 17 GHz and its long-term evolution.

We present the first results from a program searching for extended Ly$\alpha$ halos around high redshift ($ z \gtrsim 6.5$) quasars using the red channel of the Keck Cosmic Web Imager (KCWI). Our observations reveal a Ly$\alpha$ halo extending to $\simeq11$ pkpc around the $z=6.64$ broad absorption line quasar J0910$-$0414. The Ly$\alpha$ velocity field displays a rotation-like gradient, and the gas velocity dispersion is consistent with gravitationally dominated motion ($\sigma_{\mathrm{Ly\alpha}}<300$ km s$^{-1}$). Comparison with the $[\mathrm{C\;II}]$ kinematics of the host galaxy core from ALMA observations shows that the Ly$\alpha$-emitting gas extends over a much larger region, shows distinct kinematics, and has a smaller velocity dispersion ($\sigma_{\mathrm{Ly\alpha}} \simeq 0.6\sigma_{\mathrm{[C\;II]}}$). The Ly$\alpha$ spectral region of the quasar is largely obscured by a deep $\mathrm{N\;V}$ absorption trough, and as a result, roughly $55\%$ of the total Ly$\alpha$ flux is from the extended halo. These observations demonstrate the potential of KCWI for probing the cool gas reservoir that fuels the growth of quasars and their hosts in the epoch of reionization.

Previous work showed that ultralight-dark-matter solitons can provide dynamical friction for supermassive black-hole binaries, suppressing low-frequency power in the pulsar-timing-array gravitational-wave background and constraining the particle mass and effective ultralight-dark-matter fraction. Here we extend that analysis by comparing the predictive performance of four models: simplified and realistic ultralight-dark-matter implementations, a phenomenological environmental-hardening model, and a gravitational-wave-only model. We use Bayesian leave-one-out cross-validation on the five lowest pulsar-timing-array frequency bins. The phenomenological model gives the largest expected log predictive density, but its advantage over the other models is not large compared with the estimated standard errors. The current data therefore do not decisively prefer one model overall. The clearest pairwise result is within the ultralight-dark-matter framework: the simplified model outperforms the realistic implementation in all five frequency bins. Current pulsar-timing-array data are therefore compatible with ultralight-dark-matter-induced low-frequency suppression, but do not yet distinguish ultralight-dark-matter significantly from more generic environmental descriptions of supermassive-black-hole-binary evolution.

The Hubble tension is usually expressed as a discrepancy between the low H_0 inferred from Planck CMB data within base \LambdaCDM and the higher value obtained from late-time distance-ladder measurements. This scalar comparison compresses distinct inference problems into one derived parameter: Planck CMB, DESI DR2 BAO, and Pantheon+SH0ES constrain physical densities and acoustic scales, ruler-normalized distances, and calibrated luminosity-distance relations, respectively. We reformulate the comparison in terms of the dimensionless expansion history E(z)=H(z)/H_0. This does not remove the absolute-scale discrepancy, but separates the normalization encoded in $H_0$ from the redshift-dependent shape of the expansion history. Within a common flat-\LambdaCDM framework, each probe posterior is mapped onto posterior-implied E(z) histories. Since the reconstructed values E(z_k) are strongly correlated across redshift, we quantify the global mismatch with a covariance-subspace history displacement S_{hist}, alongside pointwise redshift differences. The histories are not identical, but the discrepancies are moderate: the pointwise significance is typically 1-2\sigma, while S_{hist} simeq 1.65 for DESI DR2 and S_{hist} \simeq 2.55 for Pantheon+SH0ES relative to Planck. With two retained covariance eigenmodes, these correspond to two-sided one-dimensional Gaussian equivalents of approximately 1.1\sigma and 2.1\sigma, both below the conventional \simeq 4.9\sigma Planck-SH0ES scalar-H_0 discrepancy.

IY Lyr, historically misclassified as an eclipsing binary, is now established as a first-overtone RR Lyrae star (RRc star). Using multi-band photometry (ASAS-SN, ZTF, TESS, and our BVRI data), LAMOST spectroscopy, and Gaia astrometry, we investigate its pulsation, binarity, and Galactic population. From O-C analysis, we detect a long-term period decrease and a light-travel time effect with an orbital period of 3.94 years, eccentricity of 0.46, and a mass function of 0.65 M$_{\odot}$. The companion is independently confirmed by radial velocity residuals and Gaia proper motions. Combined constraints yield an orbital inclination of 94.8$^{\circ}$ and a companion mass of 1.37 M$_{\odot}$. Chemical abundances ([Fe/H] $\simeq$ -1.0, [$\alpha$/Fe] $\simeq$ +0.27, Xiang et al. 2019) and dynamics ($L_{\rm z}$ $\simeq$ 1250 kpc km s$^{-1}$, $Z_{\rm max}$ $\simeq$ 1.31 kpc) identify IY Lyr as an old, high-$\alpha$, thick-disk star. The companion mass lies at the peak of the neutron star mass distribution, and the system's age excludes a main-sequence star; we conclude the companion is most likely a typical neutron star, although a massive white dwarf near the Chandrasekhar limit cannot be ruled out. IY Lyr is among the few RRc binaries with a compact companion verified by multiple methods, and it has important implications for thick-disk binary evolution and neutron star formation.

We present a detailed analysis of the evolution of type-C quasiperiodic oscillations (QPOs) observed during the flaring state of the recently discovered black hole X-ray binary Swift J1727.8-1613, utilizing data from the Insight Hard X-ray Modulation Telescope. By examining the relation between the QPO fractional rms amplitude and QPO frequency across various energy bands, we discover that the behavior significantly differs between these energy bands. Below 10 keV, the QPO fractional rms generally decreases with increasing QPO frequency, whereas above 10 keV, the QPO fractional rms remains relatively stable with frequency. Additionally, we report, for the first time, the detection of a common break at around 4 Hz in the relation between QPO fractional rms and frequency in both the 2-4 and 50-100 keV energy bands. We also find that the evolution of all the spectral parameters alters its behavior at around 4 Hz, with the changes in all parameters becoming flatter. This suggests a significant change in the geometry of the accretion flow. We attribute the observed break to the overall changes in the spectrum.

Silicon monoxide (SiO) traces the physical conditions and dynamics in the circumstellar envelopes (CSEs) of AGB stars. We present high-resolution ALMA Band 6 observations of highly excited SiO emission in 14 oxygen-rich AGB stars. We cover transitions from v = 0 to v = 8, including first detections of 28SiO v = 3, 4, 8, J = 6-5, 29SiO v = 6, J = 6-5, and 30SiO v = 4, 5, J = 6-5, some of which are masers. The v = 8 transition is the highest v-state observed in an AGB star yet. Masers in v = 0 are detected clearly in V PsA and IRC+10011 and tentatively in T Mic. R Hya exhibits the richest SiO spectrum. SiO J = 6-5 absorption is seen in R Aql, R Hya, S Pav, and T Mic, with features indicative of both infalls and outflows, and tentative detection of 28SiO v = 8, J = 6-5 absorption is found towards S Pav and R Aql. Highly excited SiO emission is often distributed in arcs or clumps with velocity gradients; components in R Hya and U Her align with predicted shock fronts. Detection rates show no significant difference between low and high mass-loss rate stars, although line overlap may affect some intensities. Maser detections appear uncorrelated with pulsation period or phase. The radius enclosing 90 per cent of compact SiO emission shows a tentative correlation with mass-loss rate. These results highlight the role of mass loss and CSE geometry in shaping high-excitation SiO emission.

Coronal hole boundaries are the interfaces between closed and open magnetic field regions in the solar atmosphere. Many fundamental processes take place at these regions, including magnetic reconnection that is responsible for solar wind release and restructuring of the solar magnetic field. In this paper, we present a case study in which we investigate the physical properties of the boundary of a large low-latitude coronal hole. Differential Emission Measure analysis is used to derive the plasma properties of these regions. We also apply correlation dimension mapping analysis to measure the irregularities of the coronal hole boundary. We find that the leading boundary of this coronal hole has a slightly higher average plasma temperature, is associated with a stronger and more unipolar magnetic field, and has a smoother boundary line than the trailing counterpart. These differences are hypothesised to be direct consequences of the local magnetic field configurations of the coronal hole boundary: the leading boundary corresponds to large, well-organised coronal loops, and the trailing boundary corresponds to more dispersed, randomly orientated small magnetic bipoles. Hence, we suggest that the surrounding magnetic field structure and the nature of magnetic reconnection influence the properties of coronal hole boundaries.

Understanding Jupiter's zonal winds is crucial to unraveling the dynamics of its atmosphere. Over the last decades, multiple data sources and techniques have been used to study zonal winds in Jupiter. Here, we develop a correlation-based method for the near-infrared data from the Cassini spacecraft to investigate zonal winds at different altitudes. The new method uses Jupiter's rotation to scan the planet as it rotates, allowing retrieval of winds from the analysis of light-curves of specific pixels over the Jovian disc. The method allows the retrieval of winds at multiple wavelengths from the Cassini/VIMS spectral data despite the low spatial resolution and the non-uniform cadence of the data. By applying this method to two VIMS data cubes acquired on 15 January 2001 at 09:42 UT and 16 January 2001 at 03:22 UT, we reveal the zonal winds at five main latitudes using information from three different wavebands, as well as the wind vertical structure at the equator, showing significant vertical wind shear in the troposphere. The vertical wind shear we derived is weaker than as reported in previous studies, highlighting the intricate interactions among multiple dynamical processes in Jupiter's atmosphere and reflecting the complexity of its atmospheric circulation. Despite the uncertainty due to the low spatial/temporal resolution and non-uniform cadence of the Cassini/VIMS-IR spectral data, the new method established in this study maximizes the value of the Cassini/VIMS in understanding Jupiter's zonal winds. Further observations are essential to explore the underlying mechanisms in Jupiter's atmosphere.

Simulating coronal mass ejections (CMEs) from their origin in active regions (ARs) to their propagation to Earth remains challenging, particularly when aiming to resolve AR scales and employ realistic magnetic field strengths without compromising computational efficiency. Here we present a methodology for end-to-end CME modeling that addresses these challenges. Three nested magnetohydrodynamic simulations are coupled to jointly cover the heliosphere from solar surface to beyond $1.5$ au. A block-structured adaptive mesh refinement scheme is employed to achieve $\sim 700$ km resolution in the low corona, allowing AR scales to be resolved while maintaining the total grid count below $10^8$ across the entire computational domain. A semi-relativistic Boris correction combined with a relativistic mass-density factor is used to handle magnetic field strengths up to $10^3$ G without prohibitively small time steps. Using this model, we simulate the emergence of a bipolar AR into the corona, the initiation of a CME by shearing of the AR core field and the subsequent evolution. Our simulation captures its pre-eruption energy buildup, triggering by magnetic reconnection, rapid acceleration, and propagation to 1 au and beyond. The simulated CME exhibits a three-part structure in synthetic coronagraph images and a torus-shaped flux rope in the heliosphere, with synthetic in-situ observations showing shock formation, density compression, and a prolonged southward $B_z$ component at 1 au. The entire simulation requires about one day on a moderately sized cluster (e.g., $600$ processors), while the simulated CME takes three days to arrive at $1$ au, offering a lead time of two days if used for forecasting.

The particle nature of dark matter (DM) remains one of the most significant enigmas in modern cosmology. Axion-like particles (ALPs), as well-motivated candidates for cold dark matter, can undergo radiative decay into photon pairs, a process that is significantly enhanced in the presence of ambient radiation fields. In this work, we propose a novel probe of $\mu{\rm eV}$-scale ALP DM by cross-correlating radio intensity mapping (IM) with the large-scale galaxy distribution from the 2MASS Redshift Survey (2MRS) in the local universe ($z\leq 0.1$). We develop a comprehensive theoretical framework that incorporates stimulated decay effects driven by both the Cosmic Microwave Background (CMB) and a bottom-up modeled extragalactic radio background (ERB). By forecasting the sensitivity of the Square Kilometre Array (SKA) Phase 2, we demonstrate that this cross-correlation technique provides a promising and complementary approach to searching for ALP DM signals. This study establishes a new proof-of-concept for utilizing next-generation radio telescopes to probe ALP dark matter on cosmic scales.

We investigate the physical origin of kinematically persistent planes (KPPs) of satellite galaxies in a sample of 190 Milky Way (MW)/M31-like host-satellite systems drawn from the TNG50 simulation. Building on the identification of 46 early KPPs in a previous work, we analyse their formation in the context of the high-redshift evolution of the local Cosmic Web by tracking the deformation of the so-called Lagrangian Volumes (LVs) surrounding each system. Using a reduced tensor-of-inertia analysis, we characterise the time evolution of the principal directions of collapse and relate them to the clustering of satellite orbital poles. We find that in approximately 67\% of KPPs satellite orbital poles align with the LV direction of strongest collapse, $\vec{e}_3$, while a smaller fraction ($\sim20\%$) align with the intermediate axis, $\vec{e}_2$; alignments with the major axis are rare. These alignments are statistically distinct from random expectations and reflect the confinement of satellites to planar configurations normal to the corresponding LV principal directions. We perform a kinematic analysis of satellite motion within KPPs, finding that vertical and radial motions relative to these KPPs decay early, leading to rotation-dominated, ``disky'' configurations. The characteristic timescales for satellites to settle onto a common orbital plane, for satellite orbital pole clustering, and for LV shape evolution are found to be quasi-coeval, peaking at a Universe age T$_{\rm uni}\sim4$~Gyr, during the fast mass assembly phase of the host halo. These results support a scenario in which early KPPs are fossil remnants of high-redshift, anisotropic mass collapse driven by the local Cosmic Web formation process in $\Lambda$CDM.

We develop a model-independent framework to quantify phase-space "crystallization", the degree of ordered radial and kinematic substructure, in 79 Galactic globular clusters using the Gaia EDR3-based membership catalogue of E. Vasiliev & H. Baumgardt (2021a). We construct a scalar crystallization index, C_index, by combining a radial inhomogeneity metric (z_rad) and a local, cluster-centric tangential-velocity metric (z_vel) standardized against empirical nulls. The population distribution is strongly non-Gaussian: most clusters are consistent with smooth, equilibrium expectations, while a small high-C tail (C_index >= 2) identifies dynamically complex systems, including NGC 5139 (\omega Cen) and NGC 104 (47 Tuc). Correlation and fixed-N tests show that sample size affects detectability, but does not by itself explain all high-rank objects. Through synthetic injection tests in dynamically "quiet" control clusters, we demonstrate sensitivity to ultra-cold, shell-confined kinematic components, ruling out single-shell structures comprising more than a few to ~ 10-20% of core stars in the best-sampled control clusters. We find no evidence, within the sensitivity of the adopted diagnostics, for phase-space structures that require explanations beyond known dynamical processes. However, C_index provides a useful tool for ranking clusters by dynamical extremeness, serving both as a diagnostic for internal complexity and as a quantitative metric for prioritizing follow-up dynamical or technosignature-oriented observations.

Astrophysical current layers, e.g., in pulsar winds, can be electrically charged, while the plasma is charge-symmetric, $e^\pm$. Using PIC simulations, we investigate dynamics and plasmoid formation (tearing instability) in charged Harris-type and rotational current layers. Electrically charged current layers, initially in global force-balance, are electrostatically unstable: the resulting dynamics is an intricate interplay between electrostatic Bernstein waves (BWs) and the current tearing mode. Besides overall density and magnetic field, plasma temperature is an important factor. In the charged Harris sheet set-up, the quickly generated BW are trapped within the layers (internally reflected at the upper hybrid resonance). BWs quickly redistribute the charge modifying the initial stage of tearing, but without strongly affecting overall plasmoid growth; resulting plasmoids are mildly charged. In rotational current layers: (i) even initially overall uncharged configurations develop large fluctuations of charge density; (ii) overall dynamics depends on the initial overall temperature; (iii) for certain combination of parameters tearing rate is greatly increased in the charged case.

Ultraluminous X-ray sources (ULXs) are mostly extragalactic non-nuclear point sources having X-ray luminosity exceeding the Eddington luminosity of 10 $M_\odot$ black hole i.e., $L_X \geq $ 10$^{39}$ erg ~s$^{-1}$. They are observed in all types of galaxies; spirals, ellipticals and dwarf irregulars. But the rate of occurrence of ULXs per galaxy varies, some might host a single ULX, whereas some host a large number. In this work we attempt to identify possible differences in ULX properties between two extreme categories in spirals and ellipticals, i.e. ULXs occurring at a rate of one per galaxy ($N=1$) and those occurring at larger rate. We adopt an effective scheme to generate flux limited, credible samples corresponding to the two groups in spirals and ellipticals. From this study, we infer the presence of a separate population of ULXs in the $N=1$ spiral group which contains a reasonable fraction of both soft and hard sources, while the remaining categories contain mostly harder sources. We also find six ULXs in $N=1$ ellipticals with globular cluster association. In addition, we identify few luminous candidates likely hosting massive accretors. This study provides crucial hints of a potential link between ULX types and their occurrence rates and host morphology, a finding that warrants validation via targeted observations and detailed spectral analysis of these sources.

We investigate the spectrum and polarization of radiation emerging from a subcritical X-ray pulsar using self-consistent radiation-hydrodynamic simulations of an accretion channel in a strong magnetic field. The polarized radiative transfer in the channel above the hot spot is simulated for the two normal modes, taking into account resonant Compton scattering in a strongly magnetized plasma and the effects of vacuum polarization. We show that the deceleration of the accreting matter in the subcritical regime is mainly governed by resonant scattering. Our simulations provide the velocity profiles of the plasma flow and demonstrate that vacuum polarization dominates over plasma birefringence, enhancing both the cyclotron spectral feature and the radiative deceleration of the plasma. The linear polarization degree changes sign at photon energies above the cyclotron resonance when vacuum polarization is included. We also find that the centroid energy of the cyclotron scattering feature increases with accretion luminosity, indicating a positive correlation consistent with previous observational results and theoretical interpretation.

In recent years, improvements in galaxy cluster observations have enabled a variety of tests of fundamental physics using these systems. In this work, we test the constancy of the speed of light, $c$, by combining X-ray gas mass fraction measurements from galaxy clusters with SNe Ia luminosity distance measurements from Pantheon+. We adopt the SH0ES prior on $H_0$ and the $\Omega_b/\Omega_m$ ratio from galaxy clustering observations, thereby minimizing the dependence of our analysis on any specific cosmological model. We explore different assumptions for the cluster mass calibration (mass bias), including \textsc{CLASH}, \textsc{CCCP}, and Planck-based estimates. We find no deviation from a constant $c$ when adopting \textsc{CLASH} or \textsc{CCCP} priors, while Planck-based calibration yields a mild tension, with the hypothesis of constant $c$ being only marginally consistent at the $2\sigma$ level, indicating a non-negligible sensitivity of the results to the adopted calibration scheme.

We present a photometric and spectroscopic analysis of the fast-declining Type II SN 2020aze, observed in optical bands from 2.2 to 137.4 days after explosion. The V-band light curve reaches a peak absolute magnitude of about minus 16.97$\pm$0.20 mag by 15 days, followed by a recombination phase with a decline rate of 2.04$\pm$0.13 mag per 100 days, lasting about 120 days. Early spectra (younger than 6 days) show a transient weak narrow emission line at 4687 Angstroms and a feature spanning 4400-4800 Angstroms, attributed to narrow and broad blue-shifted He II 4686, indicating interaction between the ejecta and dense circumstellar material. Comparison with spectral models suggests a red supergiant progenitor with a weak wind and a mass-loss rate of about 1e-3 solar masses per year. Semi-analytical light-curve modeling gives an initial radius of about 1100 solar radii, an ejecta mass of about 12 solar masses, an explosion energy of about 1.5e51 erg, and a progenitor mass of about 14 solar masses. These early interaction signatures, the steep decline, and the extended photospheric phase highlight the role of pre-supernova mass loss and circumstellar interaction in shaping the diversity of Type II supernovae.

Gravitational waves (GWs) from distant sources such as inspiralling and merging stellar-mass compact binaries, intermediate-mass and supermassive-binary-black-hole can be gravitationally lensed by intervening objects, ranging from stars and primordial black holes to galaxies and clusters. Depending on the GW wavelength relative to the lens scale, lensing occurs in two regimes: geometric optics, producing multiple images with time delays and magnifications, and wave optics, resulting in frequency-dependent waveform modulations. Lensed signals are identified via parameter overlap between event pairs or characteristic frequency-dependent modulations that distinguish them from unlensed signals. Detection rates depend on the redshift and mass distributions of sources and lenses, with promising prospects for future observatories. Once confirmed, lensed GWs will be powerful probes of fundamental physics and cosmology: they can constrain dark matter, lensing structures, the Hubble constant, and other cosmological parameters. In this review, we provide a concise overview of GW lensing, covering the theoretical framework, predicted detection rates, search strategies, and applications. We conclude with prospects and future directions for observing and exploiting lensed astrophysical GW events.

The strong correlation between energy injection and mid-infrared (mid-IR) emission observed during solar flares can be used to probe energy deposition throughout the chromosphere, since the IR tracks prompt flare-induced changes in electron density. Despite its diagnostic value, solar mid-IR observations are relatively recent, with sporadic campaigns over the last decade resulting in only a few recorded flares. Earlier studies found time lags between mid-IR emissions from spatially resolved footpoints, offering clues about flare energy transport. Building on this, we analyse the time lags between emissions at two wavelengths (5.2 micrometers and 8.2 micrometers) for each footpoint. Using a local cross-correlation function, we show for the first time that the 8.2 micrometers emission channel peaks 0.3 s-0.45 s before the 5.2 micrometer channel. We investigate the origin of this lag, obtaining infrared emission estimates using results from the RADYN radiation hydrodynamics code. The theoretical lag values fall within the range of the observed ones. Variations in opacity-primarily due to flare-induced ionization-explain the wavelength-dependent temporal shift between emission maxima. In particular, longer wavelengths exhibit a smaller lag between the peak of energy injection and peak of intensity. These results contribute to a better understanding of how energy deposition during a flare affects the chromospheric layers of the atmosphere. Future observations with higher temporal resolution could exploit measurements of these time lags to more fully characterize the dynamics of energy deposition during solar flares, opening a new avenue for studying heating and energy transport processes in the solar atmosphere.

The third data release of the LOFAR Two-metre Sky Survey provides an unprecedented view of the northern sky at 144 MHz. While compact sources can be efficiently identified with automated software packages, the detection of diffuse radio emission associated with galaxy clusters still requires dedicated processing and visual inspection. Given the scale of current and forthcoming radio surveys, automated approaches based on artificial intelligence are becoming essential to the identification of the most interesting targets. We aim to develop an automated pipeline to construct a catalogue of galaxy clusters hosting diffuse radio emission from LoTSS-DR3 20arcsec images. The pipeline is designed to provide both the probability that a cluster hosts diffuse radio emission and an interpretable image of its shape and morphology. We employed Radio U-Net, a convolutional neural network optimised for image segmentation (i.e. pixel-level identification) of diffuse radio emission. To associate detected emission with individual clusters, we combined the network output with positional, mass, and redshift information from four X-ray- and Sunyaev-Zeldovich-selected cluster catalogues, resulting in a merged sample of 3822 clusters covered by the LoTSS-DR3. We produced a pixel-level segmentation map of the full LoTSS-DR3 and a quantitative indicator for the presence of diffuse emission in each cluster. This enables the selection of sub-samples with specific properties for targeted follow-up or statistical studies. As a demonstration of the first application, we identified a sub-sample of 357 clusters selected at the highest network accuracy (76%), and we showed some examples of newly detected systems. For the second, using a larger statistical sample, we verified that the detection fraction of diffuse radio sources in the four catalogues increases with the mass and redshift of the clusters. [Abridged]

We present Very Large Telescope/Xshooter spectrophotometric observations of eleven low-redshift (z<0.085) compact star-forming galaxies (`high O32 sample'). These galaxies are characterised by extremely high emission-line ratios [OIII]$\lambda$5007/[OII]3727, ranging from 11 to 42. Galaxies with such high ratios are thought to be promising candidates for leaking large amounts of Lyman continuum radiation. They are characterized by low oxygen abundances 12+log(O/H)\,=7.5-8.0 and low stellar masses M*~10^6-10^8 Msun. Strong emission lines of various ions in all spectra are used to derive helium and oxygen abundances, and N/O, Ne/O, S/O, Cl/O, Ar/O and Fe/O abundance ratios. We also derived macroscopic velocity dispersions sigma(lambda) from various emission lines of different ions. We find that sigma(4861) of the Hbeta emission line is increased with increasing stellar mass and decreasing O32 ratio. On the other hand, sigma(lambda)/sigma(4861) ratios for various lines are close to 1. Exceptions are sigma(lambda)/sigma(4861) of two lines, HeII 4686 and HeI 10830, which are considerably higher than unity and of four lines, [OII] 3726,3729, [SII] 6717,6731, with sigma(lambda)/sigma(4861) lower than unity. The two former lines are likely produced in the inner parts of HII regions and are broadened by dynamical processes generated by massive stars, and by radiative scattering in the case of the HeI 10830 emission line. Emission in the four latter lines is produced mainly in the outer and likely more quiet parts of HII regions.

The Baikal-GVD neutrino telescope is the largest neutrino detector of its kind in the Northern Hemisphere. Muons produced in neutrino interaction in the vicinity of the detector leave track-like response in the detector allowing to reconstruct the neutrino arrival direction with the precision up to 0.2 degree. The Baikal-GVD collaboration has developed a variety of methods for the track-like event analysis. Methods for track-like event direction and energy reconstruction and neutrino cadidate event selection are discussed in this report. Preliminary results of application of analysis pipeline to the data-taking seasons from 2019 to 2021 are shown.

Accurate characterization of atmospheric optical turbulence is essential for evaluating astronomical sites and optimizing adaptive optics systems. The Multistar Turbulence Monitor (MTM) infers the vertical distribution of the refractive-index structure constant Cn2(z) from differential image motion measured between multiple stellar pairs in short-exposure frames. We present a comprehensive investigation of the MTM method, combining theoretical analysis, instrument-performance assessment, numerical simulations, and on-sky observations obtained at the Daocheng Astronomical Site. Simulations based on a standard HV turbulence model demonstrate that the inversion pipeline robustly recovers both the integrated seeing and the vertical turbulence profile under realistic centroiding noise and varying pixel scales. The Markov Chain Monte Carlo (MCMC) inversion achieves stable results with thirteen discrete height nodes and provides reliable uncertainties. Three nights of MTM measurements at the Daocheng Astronomical Site show that MTM-derived seeing closely tracks simultaneous Differential Image Motion Monitor (DIMM) results, accurately reproducing both short-term fluctuations and nightly averages. These results confirm that MTM provides a simple, portable, and versatile solution for atmospheric turbulence profiling and routine seeing monitoring.

47 Tucanae is one of the largest, brightest, and closest globular clusters to Earth. It hosts an exotic stellar population with stellar dynamics that indicate a complex evolution history. The cluster contains a large number of X-ray binaries and millisecond pulsars. However, given its large distance relative to the known pulsar population, previous surveys have found only the very brightest sources. Therefore, surveys with increased sensitivity should find many additional pulsars. Increasing the number of pulsars is crucial to investigate the dynamics of this globular cluster and could also lead to the discovery of unusual types of system. With a significantly increased sensitivity compared to earlier telescopes, MeerKAT is the natural choice to perform new surveys. We carried out two campaigns with different observational cadences to account for the high scintillation along the line of sight to this cluster. Here we report the discovery of fifteen new pulsars in 47 Tucanae with MeerKAT. These discoveries bring the total number of known pulsars in this globular cluster to 42, and the MeerKAT discoveries in this cluster to 17. We discuss some of their characteristics, which include preliminary localisations and estimates of orbits for most systems. Highlights include the discovery of 47 Tuc af, a 'black widow' pulsar with a short orbital period that was identified optically in 2002 as a candidate binary pulsar, and 47 Tuc ai, an eccentric binary pulsar with a massive companion, a unique system in 47 Tuc to date. Apart from the new systems, we also re-detect and localise 47 Tuc P and V, two elusive, seldom-detected systems that had no precise localisation from a phase-connected timing solution. The localisation of 47 Tuc V places it in a position consistent with a continuum source detected earlier in MeerKAT imaging data.

This work presents the first dedicated optical--mid-infrared time-domain variability and dust-reverberation analysis of the isolated NLSy1 galaxies Mrk 42 and Mrk 493. We combine ZTF optical light curves, WISE mid-infrared monitoring, archival Swift and SDSS data, and high-cadence IAC80 optical observations. Using colour--magnitude relations, flux--flux analysis, and interpolated cross-correlation functions, we trace variable optical continuum and delayed dust response. Both galaxies show positive optical--MIR lags consistent with dust reverberation. For Mrk 493, we measure an observed-frame g--W1 lag of $\tau_{\rm obs}=79.4\pm2.2$ d, corresponding to $R_{\rm dust}(W1)\simeq0.0648$ pc. For Mrk 42, the corresponding lag is $\tau_{\rm obs}=39.1\pm2.6$ d, giving $R_{\rm dust}(W1)\simeq0.0320$ pc. These lags provide optical--MIR dust-reverberation radii and BLR--dust scale comparisons for both objects; the resulting $R_{\rm dust}/R_{\rm BLR}$ ratios are $\simeq6.8$ for Mrk 493 and $\sim6$--7 for Mrk 42. For Mrk 42, we derive the first host-subtracted AGN continuum luminosity at 5100 Å from SDSS spectral decomposition, giving a self-consistent BLR--dust comparison on an AGN-only luminosity basis. Both galaxies have similar radial hierarchies but different colour behaviour: Mrk 493 shows significant optical and MIR bluer-when-brighter trends, whereas Mrk 42 shows strong optical but weak MIR colour variability. We also identify and analyse a major optical flare in Mrk 42 with four internal maxima spaced by 45--47 d. We interpret this signal as quasi-periodic substructure within a broader accretion-driven flare, rather than as a strictly coherent periodic process. These results indicate that, even in dynamically isolated environments, the variability of low-mass, high-accretion-rate AGNs is governed mainly by the intrinsic state of the accretion flow and its coupling to circumnuclear dust.

Many astronomical surveys prompt follow-up observations, but the decision process through which candidates are selected for follow-up can be difficult to model. This poses a challenge when inferring properties of the intrinsic population of astrophysical sources, rather than those of the set of objects detected by the survey and often-incomplete follow-up observations. We alleviate this problem by demonstrating that explicitly modeling of the follow-up selection process is not required for self-consistent inference of the intrinsic population. Using the framework of hierarchical Bayesian inference, we show that the intrinsic population can be accurately inferred even when the decision to follow up candidates strongly correlates with latent parameters of interest. We provide several worked examples, showing that the precision of posterior constraints can depend on the follow-up process and that one may have to model a population of contaminants if the initial selection is imperfect. Our result could dramatically simplify population inference that incorporates uncoordinated follow-up from multiple observers triggered by the deluge of candidates from surveys like LSST, Gaia, and next-generation gravitational-wave interferometers.

By aiming to find and characterise dozens of habitable exoplanets through the technique of nulling interferometry, the LIFE space mission will produce transformational science. One of the key parameters for such an interferometric mission is the nulling baseline length - the distance between nulled apertures, which past studies have assumed to be 10-100m. Advances in planet occurrence statistics and simulation tools allow us now to revisit this key assumption with significantly more detail, particularly with the intention to reduce the range of baselines considered due to mission implementation concerns. We utilise the LIFEsim mission simulator along with revised mathematical tools to identify whether the range of baselines could be reduced without significantly affecting planet yield and fringe tracking performance. Along the way, we also determine a new astrophysically motivated technique for choosing which baselines are optimal for a given science target. We find that indeed, LIFE could utilise a considerably shorter range of baselines, such as 25-80m, or even discrete baselines without much (<10%) loss of performance. Nevertheless, careful trade-offs between performance and implementation simplification must be made, especially considering any spectral weighting that may be required by the scientific goals, and the potential loss of target-specific baseline optimisation.

The radial acceleration relation (RAR) links observed and baryonic accelerations, and is best established in rotation curves of late-type galaxies. Pulsar timing, which measures line-of-sight (LOS) differential accelerations between the Sun and pulsars, provides a novel probe of this relation, including along directions outside the Galactic disc. By combining these pulsar differential accelerations with the acceleration at the Sun, we test whether current pulsar timing data carry information on a vector generalisation of the RAR, ${g}_{\rm obs}=\nu(|{g}_{\rm bar}|){g}_{\rm bar}$. Comparing the measured SPARC RAR (generalised to 3D) to 26 binary-system pulsars with literature accelerations, we find a reduced $\chi^2$ of 3.58, compared with 10.86 for Newtonian baryonic gravity alone. However, setting all accelerations to that of the Sun gives a reduced $\chi^2$ of 3.75, showing that this vector RAR test is dominated by the Solar acceleration with current data.

We review a new natural inflationary mechanism operated by repulsive-like primordial black holes (PBHs). In particular, working within the ``Swiss - Cheese" cosmological framework, we find that a Universe filled with PBHs, whose spacetime metric presents a repulsive-like behaviour, is characterised by an early quasi-de-Sitter cosmic expansion phase. Notably, for light PBHs with $m < 5 \times 10^8 \mathrm{g}$, evaporating before Big Bang Nucleosynthesis (BBN), one is met with an exponential inflationary phase with graceful exit and reheating proceeding through PBH evaporation. Furthermore, one finds as well that PBHs with $m \sim 10^{12}\mathrm{g}$ and abundances $0.107 < \Omega^\mathrm{eq}_\mathrm{PBH} < 0.5$ near matter-radiation equality can act as an early dark energy component, easing in this way naturally the Hubble tension.

We present an exact black hole solution in anti-de Sitter (AdS) spacetime with a Kalb-Ramond field coupled to nonlinear electrodynamics (NLED), characterized by mass, magnetic monopole charge, and Lorentz-violating parameters. The geometry admits two horizons (inner and outer) that coalesce into a degenerate horizon at a critical monopole charge. Beyond this critical point, no black hole solutions exist. In the limit of vanishing Lorentz-violating parameters, the solution reduces to the modified Kalb-Ramond and Bardeen black holes, while suitable parameter choices reproduce the Reissner-Nordström-AdS and Schwarzschild-AdS geometries. We analyze the thermodynamics of the solution by computing the Hawking temperature, entropy, specific heat, and Gibbs free energy. The NLED source introduces nontrivial modifications: the Hawking temperature displays nonmonotonic behavior with possible local extrema, the entropy deviates from the standard area law, and the specific heat may assume negative values, signaling thermodynamic instabilities. The Gibbs free energy exhibits swallow-tail structures, indicative of first-order phase transitions. Furthermore, we derive the first law of black hole thermodynamics in the extended phase space, together with the Smarr relation, and confirm their validity for the Kalb-Ramond black holes with NLED sources. Our findings highlight the rich thermodynamic structure induced by Lorentz-violating effects and nonlinear electrodynamics in AdS black hole backgrounds.

We explore the coincidence of scales between cosmic inflation and right-handed neutrinos in seesaw models. We show that inflation models, which will be tested by next-generation CMB experiments, can produce right-handed neutrinos in sufficient abundance to explain the observed baryon asymmetry of the universe. The model can be tested by gravitational wave signatures from cosmic inflation and particle production.

The annihilation of self-interacting dark matter with long-range interactions can be significantly enhanced at low velocities through the Sommerfeld effect. At special points in parameter space, where near-zero-energy resonances exist in the spectrum of the theory, the standard calculation of this enhancement appears to violate unitarity. Recently, several approaches have been proposed to regulate this behavior, some introducing explicit ultraviolet (UV) scales and others not, raising the question of whether these prescriptions are consistent. In this paper, we compare these approaches and show that even in nominally cutoff-dependent methods, the regulated Sommerfeld-enhanced cross sections are independent of the UV regulator to a good approximation, and that when the unitarity-preserving corrections are large, the regulation schemes coincide to leading order. We use these insights to write down a regulator-independent prescription for unitarization applicable to multi-state systems, where the modified enhancement can be written solely in terms of the standard enhancement factor, the hard annihilation amplitude, and the $S$-matrix for scattering in the long-range potential.

The scotogenic model defines a framework for radiative neutrino masses and provides a viable dark matter candidate. Since the scotogenic dark matter is leptophilic, indirect searches appear as an especially interesting possibility. Here we propose a simple variation of the model with a very distinct phenomenology. The scotogenic fermion singlets are naturally grouped into pseudo-Dirac pairs of mass of 0.1-1 TeV. We show that the lightest one constitutes a dark matter candidate that near threshold annihilates with a 90% branching ratio into neutrino pairs. The model gives the observed relic abundance consistently with the bounds from direct searches and with all neutrino and charged lepton data. We also show that, for a sub-MeV dark matter particle, the model suggests a scenario that could address the lithium problem.

We study electromagnetic perturbations and the associated quasinormal modes (QNMs) of parametrized static, spherically symmetric wormhole spacetimes, focusing on Damour-Solodukhin and braneworld geometries as well as their galactic extensions. Using the Bronnikov-Konoplya-Pappas parametrization, we express the metric functions in terms of a compactified radial coordinate and characterize the spacetime through far-field and near-throat parameters. The far-field coefficients govern the asymptotic structure and post-Newtonian behaviour, while the near-throat continued-fraction expansion captures the strong-field geometry near the throat. We first apply the parametrization to isolated wormholes and identify its range of validity, showing that non-polynomial metric functions can limit the convergence of the near-throat expansion and hence the accuracy of a truncated representation. We then extend the framework to a galactic Damour-Solodukhin wormhole embedded in a Hernquist dark matter halo. Imposing observational bounds from the shadow of Sgr A$^*$, we constrain the galactic compactness and deformation parameters and obtain an observationally viable parametrized metric. Within the allowed parameter space, we compute the fundamental QNM frequencies using the transfer matrix method and analyze the corresponding time-domain ringdown signals. We find that the damping rate is more sensitive to galactic compactness, whereas the oscillation frequency remains comparatively stable. Although the spectral shifts are small within the shadow-allowed region, the framework provides a systematic link between geometric parametrization, shadow constraints, and dynamical response. Our results establish an observationally consistent parametrized description of wormhole perturbations for strong-field tests of horizonless compact objects.

We propose a simplified model for the cosmological constant in string theory flux vacua motivated by type IIB and F-theory compactifications. Relative to the Bousso-Polchinski model, small vacuum energy spacing occurs in thin wafers rather than thin shells. The model is applied to the entire Schöller-Skarke database of Calabi-Yau fourfolds, which exhibit $532,600,483$ distinct sets of Hodge numbers. The overwhelming majority of those ($99.95\%$ percent for some choices of parameters) exhibit a vacuum energy spacing of~$10^{-120}$ in Planck units or smaller. Brown-Teitelboim membrane nucleation transitions can populate this landscape of flux vacua. In the thin-wall approximation, and ignoring gravitational corrections, we find that the bubble transitions are always dominated by giant leaps in flux space. The age of the universe places a bound on Calabi-Yau topology that is satisfied for the entire Schöller-Skarke database.

We demonstrate that the conversion gain of a superconducting hot-electron bolometer (HEB) mixer can be increased by biasing the device within the negative differential resistance (NDR) region of its current-voltage characteristic. Although NDR biasing has historically been avoided due to MHz-range resistive oscillations, we show that these oscillations arise from an LC resonance formed by the bias-T inductance and the effective thermal capacitance of the HEB. By applying stability criteria analogous to those developed for tunnel diodes, we redesigned the embedding circuit to suppress this resonance and achieve stable NDR operation. Direct measurements using two monochromatic 2.5-THz sources confirm the predicted gain enhancement. These results establish NDR biasing as a viable method for improving HEB mixer performance and motivate further studies of noise behavior and circuit optimization.

We develop a unified framework for testing gravity beyond General Relativity (GR) with continuous gravitational waves (CWs) from individual supermassive black hole binaries (SMBHBs). These long-lived, nearly monochromatic nanohertz signals offer unique strengths for precision tests of gravity, since their coherent phase evolution and inter-pulsar correlations in pulsar timing arrays (PTAs) retain detailed information about departures from GR over cosmological propagation distances. We consider three representative classes of deviations from GR: additional polarization states, modified dispersion relations, and parity-violating birefringence. For each, we derive the inter-pulsar cross correlation, the modified antenna response, and the propagation-induced pulsar-term phase delay. For non-tensorial polarizations, the CW cross correlation scales linearly in the alternative-polarization amplitude, compared to the quadratic scaling of the gravitational-wave background (GWB), provided the beyond-GR modes are sub-dominant. PTAs are also competitive for modified dispersion relations, where low frequencies enhance both the antenna-pattern modification and the pulsar-term phase delay. Birefringence, by contrast, is suppressed at nanohertz frequencies for most parity-violating theories. We validate the framework with injection-and-recovery simulations for breathing-mode and massive-graviton signals at current observational limits, recovering the injected beyond-GR parameters and distinguishing the CW signal from both correlated and uncorrelated background models. We further show that a pure-GR CW template recovers source parameters without significant bias when beyond-GR physics is present in the data, supporting a two-stage analysis strategy: identify candidates under GR, then test for deviations.

We investigate the optical appearance of rotating black holes embedded in dark matter environments through a phenomenological ray-tracing framework. Instead of focusing on a single dark-matter-dressed geometry, we use a comparative setting in which a regular Einasto-type dark-matter-sourced black hole and a cored-NFW halo black hole are promoted to effective rotating backgrounds. This allows us to separate profile-dependent effects from generic Kerr-like rotational features. We construct observer-screen images by numerical backward ray tracing and analyse the shadow boundary, lensing bands, transfer maps, and synthetic intensity distributions produced by semi-analytic accretion prescriptions. Particular attention is paid to the possible degeneracy between spin and the dark matter parameters, and to whether higher-order lensed images can provide stronger signatures than the primary shadow size alone. For the representative parameters considered here, the Einasto-supported geometry remains close to Kerr, while the cored-NFW case produces a larger apparent image scale and more visible deviations in the lensed structure. The framework developed in this work is lightweight and extensible, and can be used as a first step toward future radiative-transfer and low-resolution accretion-flow studies of rotating dark-matter-dressed black holes.

We discuss the properties of the extended three-flavor quark-meson diquark (EQMD) model as a renormalizable low-energy effective model for QCD. The effective degrees of freedom are quarks, scalar- and pseudoscalar mesons, diquarks, vector- and axial-vector mesons. We calculate the equation of state (EoS) in the mean-field approximation at $T=0$ imposing charge neutrality for electric and color charges. We match the EoS with a low-density nuclear equation of state. We discuss how the choice of parameters in the model affects the EoS and thereby the mass-radius for hybrid stars. We show that it is possible to construct hybrid stars whose masses and radii are in agreement with recent astrophysical observations and perturbative QCD (pQCD). The addition of vector and axial vector mesons to the quark-meson diquark is essential, since it makes the EoS sufficiently stiff for intermediate densities. Our results suggest that stars with a mass larger than $M\sim2M_{\odot}$ have a quark core with a central density $n_B\geq 3.9n_{\rm sat}$, where $n_{\rm sat}\approx0.165$fm$^{-3}$ is the saturation density. The speed of sound has a double-peak structure and relaxes to the conformal limit from above for large baryon chemical potentials $\mu_B$. This structure is caused by the decrease in the mass of the $s$ quark as $\mu_B$ increases.

We investigate freeze-in production of light dark matter through the quark flavor-changing operator $(\bar{s}\gamma_\mu d)(\bar{\chi}\gamma^\mu\chi)$ in a low-reheating cosmology. For reheating temperatures below the QCD crossover, kaon decays and scatterings generate the dark matter abundance through $K\to\pi\chi\bar{\chi}$ and $K\pi\to\chi\bar{\chi}$. The same interaction induces the rare kaon decays $K^+\to\pi^+\chi\bar{\chi}$ and $K_L\to\pi^0\chi\bar{\chi}$. This links the freeze-in relic abundance to searches at NA62, KOTO, and KOTO II. We find that lower reheating temperatures require larger couplings to compensate for the Boltzmann-suppressed kaon abundance, making kaon-driven freeze-in dark matter testable at rare kaon decay experiments.

The interpretation of the SN1987A neutrino data continues to be hindered by significant absolute timing uncertainties and a well-known tension in the angular distributions. We perform a quantitative, high-precision alignment of the Kamiokande-II and Baksan data with the IMB clock, using for the first time the relative time offset (RTO) and reducing the temporal uncertainty by two orders of magnitude to the sub-second level. The chi-square analysis shows that Baksan's absolute timestamps require an advancement of 30.4 s, while those of Kamiokande-II require a delay of about 6.4 s. The knowledge of the unified timeline provides a necessary and rigorous basis for testing whether the first Kamiokande-II event arose from the neutronization burst - a possibility motivated by its angular properties. Our analysis favors an accretion phase electron anti-neutrino origin over a neutronization-burst origin, with a likelihood ratio of 3-6, depending on the specific MSW oscillation scenario. This result corroborates the standard interpretation - that only inverse beta decay events were detected. Our framework yields the most stringent constraints to date on the SN1987A chronology and establishes a precision benchmark for future Galactic supernova observations.

Non-radial oscillation modes of neutron stars serve as diagnostics of their internal composition and relativistic structure. In this work, we investigate the perturbations of static and spherically symmetric neutron stars characterized by an anisotropic pressure. Given the background symmetry, perturbations decouple into polar and axial modes. To date, axial modes have remained less explored, primarily because matter and metric perturbations decouple in the isotropic limit. In this work, we provide a consistent treatment of axial modes and demonstrate that pressure anisotropy induces a direct coupling between matter and metric perturbations. We employ parameterized anisotropy models that ensure consistency with the treatment of matter perturbations. We numerically integrate the linearized Einstein field equations for the axial modes, employing a diverse set of realistic equations of state. Our results indicate that as the stellar mass grows, the frequency of the lower $w$-mode generally decreases, while its damping time increases. Softer equation of states typically yield slightly higher oscillation frequencies. Furthermore, larger anisotropy (i.e., when the tangential pressure exceeds the radial pressure) allows for more massive equilibrium configurations, which correspondingly leads to lower oscillation frequencies and prolonged damping times. Finally, we demonstrate that the frequency and damping time, both scaled by the stellar mass, exhibit a nearly universal quadratic dependence on the stellar compactness, remaining largely insensitive to both the underlying equation of state and the specific anisotropy model.

The IceCube Upgrade, an extension designed to enhance the IceCube Neutrino Observatory's detection of neutrinos with energies between 1 GeV and 500 GeV, will markedly improve IceCube's sensitivity to low-mass dark matter scenarios. In this study, we present sensitivity projections for the IceCube Upgrade to neutrino fluxes arising from dark matter annihilation. In particular, we consider dark matter with masses between 3 GeV to 500 GeV from both the core of the Sun and the Galactic Center. These projections indicate that the IceCube Upgrade will enable stringent limits on dark matter in this parameter space, achieving leading sensitivities to some dark matter models with only three years of data taking.