Papers for Monday, Mar 02 2026

Type Ia supernovae (SNe Ia) are common luminous astrophysical transients. SNe Ia serve as distance indicators for measuring the expansion rate of the universe and play important roles in galactic nucleosynthesis. However, ambiguities persist regarding the nature of their stellar progenitors and explosion mechanisms. The recent discovery of \textit{Gaia} hypervelocity white dwarfs (WDs) has provided direct evidence in support of helium-ignited double degenerate SNe Ia. In this study, we investigate the outcomes of helium-ignited double-degenerate WD mergers by performing a set of 3D hydrodynamical simulations with two different codes: \texttt{AREPO} and \texttt{FLASH}. We consider two distinct binary WD systems close to helium ignition, evolving each with both codes while keeping initial conditions fixed. The first binary WD model produces a double detonation of the primary WD and the hypervelocity ejection of the surviving secondary, similar to the canonical dynamically driven double degenerate double detonation (D6) scenario. In the second model, the secondary also undergoes a core detonation, resulting in the complete disruption of both WDs. Notably, despite utilizing distinct numerical solvers, nuclear reaction networks, and mesh strategies, \texttt{AREPO} and \texttt{FLASH} produce broadly consistent outcomes for both sets of initial conditions. While the nucleosynthetic yields differ due to the different nuclear reaction networks employed, the overall agreement between the simulations demonstrates the robustness of the numerical modeling of this scenario. Our results strongly support the viability of both the D6 and quadruple detonation channels for at least some SNe Ia. We explore the prospective observational signatures of this channel, including in the X-rays using \textit{XRISM's} \textit{RESOLVE}.

We conduct simulations of $\sim$kpc-scale cool clouds in the circumgalactic medium (CGM), using initial conditions sampled from a highly resolved cosmological magneto-hydrodynamical zoom-in of a Milky Way-like galaxy. We select ten distinct cold clouds with masses of $m_{\rm{cloud}}$ $\sim$ $10^{4.5-5}$ M$_\odot$, originally resolved at a mass resolution of $m_{\rm{gas}}$ $\sim$ $200$ M$_\odot$. To further resolve small-scale features and physics, we implement a targeted refinement scheme within spherical regions co-moving with each cloud, thereby boosting the local mass resolution by a factor of 1000, reaching $m_{\rm{gas}}$ $\sim$ $0.2$ M$_\odot$ (spatial resolution, $r_{\rm{gas,cloud}}$ $\sim$ $O(\rm{pc})$). The selected clouds have diverse properties, across a broad parameter space, resulting in heterogeneous evolution. For the clouds we study, radiative cooling is the dominant physical process enabling cloud survival, while magnetic fields play a comparatively smaller role. The motion of these clouds is governed not only by drag forces that decelerate them, but also by acceleration from momentum exchange with the complex background velocity field, which can cause them to move faster than ballistic projectiles set by their initial velocities. Our results suggest that the non-trivial details of realistic cosmological initial conditions -- specifically the complex density, temperature, and velocity fields -- may play an important role in subsequent cloud evolution, and that sampling the output of an existing large-scale simulation provides a self-consistent approach to capture these effects without ad hoc assumptions.

We present new X-ray, optical, and UV observations of the tidal disruption event candidate eRASSt J234402.9-352640, (hereafter J2344). Between 50 and 60 days after peak optical brightness, J2344 exhibited large-amplitude modulations in its 0.2-2 keV emission, when the flux repeatedly dimmed and re-brightened by a factor of ~6, over a ~3-day timescale. These modulations exhibited harder-when-brighter behaviour but were not detected in high-cadence observations obtained 60-70 days and 170-200 days after peak optical brightness, when the system instead exhibited stochastic X-ray variability over timescales of hours. We discuss the different physical mechanisms responsible for such exotic X-ray variability and explore the possibility that the modulations in J2344 were caused by the Lense-Thirring precession of the inner accretion flow around the disrupting black hole.

We present the most sensitive direct imaging search for the nearest ($d = 3.2$ pc) Jupiter-analog exoplanet, $\epsilon$ Eri b, with JWST/NIRCam coronagraphy between 4-5 $\mu$m (F444W). We achieve a 5$\sigma$ contrast sensitivity $\approx3.0\times10^{-7}$ ($\Delta \approx 16.3$ mag) in the F444W filter at the expected planet separation of $\approx$1". This is the deepest 4-5 $\mu$m contrast performance achieved for any JWST/NIRCam observation to date at these separations (and $>10\times$ better than ground-based limits). Yet, the planet remains elusive to imaging. We update the star's age to $1.1\pm0.1$ Gyr, older than previous age estimates, using the latest gyrochronology relations. This significantly impacts $\epsilon$ Eri b's inferred effective temperature ($T_{\rm eff}$), which is now expected to lie between 150-200 K based on evolutionary models for a 1 $M_{\rm Jup}$ planet. Using cloud-free Sonora Flame Skimmer models and custom PICASO patchy cloud models in the above $T_{\rm eff}$ range, we find that the F444W non-detection of $\epsilon$ Eri b can be explained by a metal-enriched atmosphere and/or an atmosphere containing water ice clouds. Both possibilities suggest that $\epsilon$ Eri b's atmosphere is strikingly similar to that of Jupiter in our Solar System. Alternatively, if we do not enforce the dynamical mass ($0.98 \pm 0.09\;M_{\rm Jup}$), a solar metallicity, cloud-free, $\lesssim0.81\;M_{\rm Jup}$ planet would be consistent with the NIRCam upper limit based on the Sonora Flame Skimmer evolutionary models. Finally, we place limits on the size of a potential ring system using the NIRCam/F210M data and discuss the opportunity to directly image $\epsilon$ Eri b with additional JWST observations, the Roman Coronagraph Instrument, the ExtraSolar Coronagraph on the Lazuli Observatory, and EELT/METIS.

We investigate the gravitational wave (GW) signals emitted by planetary-mass primordial black holes (PBHs) passing nearby or traversing neutron stars (NSs). While previous studies mainly focused on the detailed waveforms of the signals, we estimate the rate of PBH-NS gravitational-wave events originating from the Galactic center and compute the probability of detecting a signal over 10 years of LIGO-Virgo-KAGRA observations. We examine in detail the case of PBHs bound to NSs, focusing on eccentric orbits that give rise to repeated GW bursts emitted in correlated series, each burst corresponding to a periastron passage. Despite the enhancement from the large number of bursts produced by a single PBH-NS pair, the total number of signals produced in this way remains subdominant to those due to random unbound encounters of PBHs with NSs. We also find that both types of signals have a very small probability $P\lesssim 10^{-8}$ to be detected in a 10 year period.

Recently, an intermediate-mass black hole (IMBH) candidate was announced in the Galactic globular cluster Omega Centauri. IMBHs at the lower end of the traditional mass range have also been detected through gravitational-wave transients, though their formation and subsequent growth linking the two mass scales remains a mystery. One way IMBHs may be produced is through the collapse of very massive stars produced by stellar collisions in dense stellar environments. However, IMBHs may be ejected from such environments by either dynamical recoil from binary-single scattering or gravitational-wave recoil following the merger of two black holes. We conduct Newtonian and post-Newtonian binary-single scattering experiments to study dynamical ejection in greater detail. We obtain fits to the probabilities for dynamical ejection, gravitational wave capture, and per-encounter hardening as a function of the binary mass ratio and hardness with respect to its environment. We borrow techniques from survival analysis (commonly used in studies of medicine, epidemiology, engineering, etc.) to develop a model to calculate the probability of IMBH binary ejection vs in-cluster merger. We confirm that the dynamical ejection probability strongly depends on both the mass ratio of the IMBH compared to other BHs in its environment and on the environment's velocity dispersion. We estimate that for a typical Milky Way globular cluster, IMBHs with mass $\lesssim10^3\,\mathrm{M}_\odot$ are unlikely to be retained until the present. Our results also suggest that IMBH mergers with $q\lesssim0.2$ may be detectable at higher redshifts with future gravitational wave instruments such as the Einstein Telescope and Cosmic Explorer.

The Lagrangian Particles (LP) module of the PLUTO code offers a powerful simulation tool to predict the non-thermal emission produced by shock accelerated particles in large-scale relativistic magnetized astrophysics flows. The LPs represent ensembles of relativistic particles with a given energy distribution which is updated by solving the relativistic cosmic ray transport equation. The approach consistently includes the effects of adiabatic expansion, synchrotron and inverse Compton emission. The large scale nature of such systems creates boundless computational demand which can only be satisfied by targeting modern computing hardware such as Graphic Processing Units (GPUs). In this work we presents the GPU-compatible C++ re-design of the LP module, that, by means of the programming model OpenACC and the Message Passing Interface library, is capable of targeting both single commercial GPUs as well as multi-node (pre-)exascale computing facilities. The code has been benchmarked up to 28672 parallel CPUs cores and 1024 parallel GPUs demonstrating $\sim(80-90)\%$ weak scaling parallel efficiency and good strong scaling capabilities. Our results demonstrated a speedup of $6$ times when solving that same benchmark test with 128 full GPU nodes (4GPUs per node) against the same amount of full high-end CPU nodes (112 cores per node). Furthermore, we conducted a code verification by comparing its prediction to corresponding analytical solutions for two test cases. We note that this work is part of broader project that aims at developing gPLUTO, the novel and revised GPU-ready implementation of its legacy.

Rapid TeV flares have been observed from the core of the active galaxy M87. These have been attributed to inverse Compton scattering of disk photons by electrons and positrons accelerated in transient reconnection layers formed in baryon-poor regions of the magnetosphere of the central black hole, M87*. It was previously shown that even a small number of protons accelerated in the same layers can lead to bright GeV proton-synchrotron flares, if protons receive $\gtrsim20\%$ of the dissipated power for reconnecting fields of $\sim$100 G. We aim to investigate the role of Bethe-Heitler pair production in the emission of reconnection-driven flares from M87* in this physical regime. We perform numerical calculations that incorporate inelastic collisions between relativistic protons and photons, as well as photon-photon pair production, and compute the non-thermal radiation from the layer. The numerical calculations are also supported by analytical estimates. We find that disk photons act as targets for Bethe-Heitler pair production. The resulting pairs emit very high-energy synchrotron photons ($\gtrsim$0.1 TeV), which are subsequently attenuated by the disk photon field, leading to further pair production. The synchrotron emission of these secondary pairs produces soft photons, as part of an electromagnetic cascade, enhancing pion production and photon-photon attenuation down to tens-of-GeV energies.

Star-forming galaxies spend most of their lifetimes on the star-forming main sequence, which establishes a tight empirical and statistical relation between stellar mass and star-formation rate. Occasional episodes of rapid star formation can push them temporarily above this sequence, turning them into starbursts. Yet some galaxies display starburst-like traits -- rapid, dense, and compact star formation -- while still remaining within the scatter of the main sequence. These "starbursts in the main sequence" (SBMSs) reveal the complexity and diversity of star formation modes, making them crucial for understanding how galaxies evolve and transition between different regimes. In this paper, we identify SBMSs in the cosmological simulation NewHorizon and follow their evolution across time to uncover their physical origins and the role of this special regime in shaping galaxy evolution. We explain the existence of SBMSs by a comparatively earlier assembly of their stellar mass, driven in particular by more frequent and repeated mergers as the other galaxies, as well as exceptionally productive starburst events triggered by these interactions. As a result, this regime appears preferentially -- though not exclusively -- in the most massive galaxies. The SBMS behavior is not continuous within individual galaxies but instead arises intermittently as a short-lived (~ 30 Myr) evolutionary mode. Nevertheless, such SBMS episodes exist throughout cosmic time across the galaxy population... [abridged]

Nearby radio galaxies are among the most promising candidates for the acceleration of ultra-high-energy cosmic rays (UHECRs). In this work, we develop a physically motivated, source-resolved framework to quantify the contribution of the three nearest FR-I radio galaxies$-$Centaurus A, Virgo A, and Fornax A to the UHECR flux measured by the Pierre Auger Observatory. Acceleration spectra derived from detailed jet-acceleration models are combined with numerical simulations of extragalactic propagation, while the more distant radio-galaxy population is treated as a continuous background. By fitting exclusively the measured UHECR energy spectrum, we determine the relative contribution of each source and constrain the fraction of jet power converted into UHECR luminosity. We find that a small number of nearby radio galaxies can account for the highest-energy UHECR flux with acceleration efficiencies of order $10^{-3}-10^{-2}$, while the background contribution remains subdominant. The resulting scenarios yield mass-composition trends broadly consistent with observations and predict distinct levels of secondary neutrino fluxes. These results demonstrate that physically grounded, source-specific modeling of nearby radio galaxies provides a viable and predictive explanation for the origin of the highest-energy cosmic rays.

Gravitational-wave (GW) astronomy has advanced our understanding of compact mergers through instruments like the Laser Interferometer Gravitational-Wave Observatory (LIGO). However, the extreme sensitivity required for these detections makes the instruments susceptible to short-duration transient noise, or glitches, which obscure GW data. Current tools such as Omega Scan and GravitySpy assist in identifying and classifying such noise, but are limited by manual inspection or dependence on large training sets. To address these challenges, we present \textit{OmegaNeuron}, a machine-learning tool that integrates GravitySpy's image similarity methods with Omega Scan's transient analysis to automate the identification of auxiliary channels that witness glitches. Applied to multiple glitch examples, OmegaNeuron consistently highlighted plausible witness channels and showed strong agreement with existing correlation tools, while providing clearer ranking through a quantitative similarity metric. Integrated into the \texttt{gwdetchar} package, OmegaNeuron enables faster analysis that improves glitch witness identification, enhancing both detector sensitivity and the reliability of gravitational-wave observations.

We present nanoCMB, a minimal but accurate calculator for the unlensed CMB temperature and polarisation angular power spectra ($C_\ell^{TT}$, $C_\ell^{EE}$, $C_\ell^{TE}$) of flat $\Lambda$CDM cosmologies. Written in $\sim$1400 lines of readable Python, the code implements the full line-of-sight integration method: RECFAST recombination, coupled Einstein--Boltzmann perturbation equations in synchronous gauge with a tight-coupling approximation, precomputed spherical Bessel function tables, and optimally constructed non-uniform grids in wavenumber and conformal time. Despite its brevity, nanoCMB achieves sub-percent agreement with CAMB across the multipole range $2 \le \ell \le 2500$, running in $\sim$10 seconds on a modern laptop. The entire calculation lives in a single, easily modifiable Python script, designed as a pedagogical bridge between textbook treatments and research-level Boltzmann solvers, with every approximation and numerical choice made explicit. We describe the physics, equations, and computational methods in detail, with code snippets illustrating each stage of the calculation. The code is publicly available at this https URL.

It has long been suggested that X-ray transients are produced at periastron of stellar-compact object binaries with eccentric orbits. Recoil of matter evaporated from the star by X-rays from matter transferred at periastron increases the orbital semi-major axis and eccentricity. After periastrons the object would be a transient X-ray source, a Galactic analogue of a tidal disruption event (TDE), but recurrent with a gradually increasing period rather than catastrophic.

High-J CO emission in high-redshift galaxies has been studied primarily on an individual-source basis, limiting our ability to draw population-level conclusions about molecular-gas excitation. To address this limitation, we present a catalog of CO spectroscopy based on archival data from the Atacama Large Millimeter/submillimeter Array (ALMA) for a sample of galaxies at z > 3, focusing on high-J transitions (Jup = 9-17). Combining ALMA archival data with published measurements, we compile CO spectral line energy distributions (SLEDs) for 38 well-studied systems spanning z ~ 3.1-6.9, including 5 hot dust-obscured galaxies (Hot DOGs), 17 submillimeter-bright galaxies (SMGs), and 16 optically selected quasars. The class-median SLEDs rise steeply to Jup = 9 and remain approximately flat through Jup ~ 11-12. SMGs show relatively stronger low- to mid-J emission relative to CO J=9-8, while Hot DOGs exhibit tentative evidence for higher excitation. Comparison with simple excitation models suggests that X-ray dominated region (XDR) heating or dense, shock-heated gas can account for the extended high-J CO SLEDs. A tentative anti-correlation between the CO(9-8)-to-infrared luminosity ratio and excitation among the dusty galaxy populations suggests that the enhanced excitation in Hot DOGs may be driven by XDR heating from obscured AGN activity rather than by shocks.

Extending the sensitivity of terrestrial gravitational-wave detectors below 20 Hz is a long-standing challenge, limited by ground motion and inertial sensing noise. In this letter, we demonstrate ultra-high-vacuum compatible inertial isolation and position sensing technologies that achieve active platform stabilization down to 10 mHz. Our laser position sensors reach a sub-pm/$\sqrt{\rm Hz}$ sensitivity above 10 mHz, independent of the input light polarization, representing a 100-fold improvement over the current LIGO position sensors. In addition, our inertial sensors provide at least a factor of 5 improvement in low-frequency sensitivity compared to state-of-the-art commercial seismometers. We integrate these technologies into a LIGO-like interferometer model and predict a low-frequency sensitivity improvement of up to an order of magnitude at 10 Hz, with enhanced linearity and calibration stability. This extension increases the detection horizon for intermediate-mass black hole binaries of mass $10^3 M_\odot$ by a factor of 3. Our results provide the first experimental demonstration of a practical pathway to sub-10 Hz operation of terrestrial gravitational-wave detectors and establish key technologies for next-generation observatories such as Cosmic Explorer and Einstein Telescope.

In our study spanning 2015-2021, we examined sink pixels (SPs) in the Advanced Camera for Surveys Wide Field Channel (ACS/WFC) using dark and SP reference files. SPs are pixels with values $\le$ $-10$ electrons below the local background of LED-flashed short (0.5 sec) darks, that collect and trap significant charge during readout. Analyzing seven years of short dark data, we assessed SP creation and persistence. In this time frame, 5,430 SPs were created in WFC1 and 5,649 SPs in WFC2, with creation rates of about 2.15 pixels/day and 2.23 pixels/day, respectively. These calculations allowed us to detect 44,068 SPs, not including SP trails, in the detector by the end of 2021, constituting approximately 0.25\% of the science frame. We found it is rare for SPs to return to a typical, non-negative pixel value. We observed more flagged SPs near the serial register than the chip gap. Skewed histograms for the $y$-position distribution, exhibiting a local peak in the distribution of SPs very near the chip gap described as the ``bounce-back" effect, were evident for both WFC1 and WFC2, while the $x$-position distribution remained uniform. Examining CTE-corrected images from 2015, 2018, and 2021 revealed consistent trends, with the gradient getting steeper over time due to CTE losses, which is also worse for pixels further from the serial register. We simulated the CTE-impacted readout of a short dark exposure with uniformly distributed SPs, to assess how CTE influences SP detectability. While the gradient effect was reproduced, the local peak near the chip gap was not. Filling in of SPs by CTE charge-release during readout appears to explain most of the gradient in the $y$-position density of SPs.

The Solar Blind Channel (SBC) typically exhibits elevated dark rate levels at temperatures exceeding 25.5°C. However, instances of rapid dark rate increases have been observed before the detector reaches this threshold. To more closely monitor these anomalies, the existing calibration program was expanded, scheduling 24 total orbits per year distributed across eight visits for Cycle 31 and beyond. With enhanced data availability, we provide further updates and analysis of the dark rate in this report. We investigated the influence of the South Atlantic Anomaly (SAA) passage on dark rates. Using orbital parameters from SPT products, we plotted the Hubble Space Telescope's (HST) locations during dark exposures relative to the SAA. Comparison of HST paths during elevated and stable dark rate visits revealed no significant correlation between proximity to the SAA boundary and dark rate increases. Furthermore, we examined the stability of dark rates in the vicinity of the SBC-LODARK aperture within the detector. Our findings confirm that this region maintains consistently low dark rates across visits, unaffected by the elevated dark current observed elsewhere in the detector. The SBC-LODARK aperture therefore continues to be recommended for small sources.

Interstellar objects provide a unique view into the formation of other star systems. Here we present spectroscopic observations of the recently discovered interstellar object 3I/ATLAS between a heliocentric distance of $3.7$ to $1.8$~au on either side of its travels through perihelion. We obtained several observations with the Keck-I/LRIS, Keck-II/NIRES, Gemini/GMOS, and UH88/SNIFS spectrographs, covering a wavelength range of $0.3 - 2.5~\mathrm{\mu m}$. We report the continued emission of both Ni and CN, along with post-perihelion detections of Fe and a weak detection of $\mathrm{C_3}$. We determine the spectral slope across optical and NIR wavelengths and find a positive spectral slope in the optical, with values ranging from $\sim 21 - 27\%$ in the blue regions ($0.4 - 0.55~\mathrm{\mu m}$) to $\sim 6 - 10\%$ in the red ($0.65 - 0.9~\mathrm{\mu m}$) regions. In contrast, the NIR showed a negative spectral slope of $\sim -0.9 \%$ between $0.9 - 1.5~\mathrm{\mu m}$ and $\sim -2.3\%$ between $1.9 - 2.5~\mathrm{\mu m}$. 3I/ATLAS shows a clear turnover in its spectral shape at $\sim 1.1~\mathrm{\mu m}$, corresponding to scattered light from the dusty coma. Finally, in the NIR, we do not find an increase in the depth of the water features identified in an earlier NIR observation of 3I/ATLAS. Our observations of 3I/ATLAS in the NIR show a similar shape to the NIR spectrum of 2I/Borisov as it approached perihelion.

When two massive stars orbit each other, their winds create a shock cone. In some cases, an evolved, carbon-rich Wolf-Rayet (WR) star's wind collides with that of an orbiting OB star, condensing into dust downstream. This dust is then seen as large spiral structures that eventually move into the interstellar medium. Among these colliding wind binaries, the archetype system WR104 has become an enigma. Aperture masking interferometry with Keck revealed an evolving face-on dust spiral with multiple rungs of dust visible from years of observations. In contrast to direct imagery, recent spectroscopic results implied that the orbit must have an inclination quite different from the face-on geometry. We examined the ASAS and ASAS-SN photometry to put further constraints on the geometry of the orbit. Through a phase-binning of the light curve, we find that the recent g-band light curve is brightest at a time when the OB star is in front of the WR star in our line of sight, with the lowest flux happening at the opposite conjunction. We fit the light curve with an illustrative model for scattering eclipses, which then allows us to infer an inclination of the system of $(41.8^{+13.0}_{-14.9})^\circ$. This inclination agrees with the recent spectroscopic orbit and presents challenges to previous interpretations of high-angular resolution images of the dust plume. We provide a qualitative geometric model for the dust plume to reconcile these results and show how WR104 can provide a means to study the properties of WR dust in detail.

The solar cycle - most notably characterized by its sunspot activity patterns - serves as a cornerstone of heliospheric physics. This research uncovers a fundamental magnetic dichotomy in the Sun's full-disk field, identifying two functionally separate populations: the Strong-Field Group (SG) and Weak-Field Group (WG). The solar cycle exhibits a dual nature, much like Janus, with the SG and WG operating in opposing phases regardless of low or high latitudes. The SG-dominated cycle represents one facet of this duality and is visually prominent at the solar surface. It is well-established that this component synchronizes with the sunspot cycle at low latitudes but operates in anti-phase at high latitudes. In contrast, the WG-driven cycle acts as its hidden counterpart, functioning in opposition to the SG at both high and low latitudes -- a behavior that had not been identified until now. Influenced by these magnetic field groups, this dual nature permeates the entire solar atmosphere, revealing that the full-disk solar activity is globally modulated by the Janus cycle.

Binary neutron star mergers are thought to produce hot, rapidly rotating neutron stars with masses that can far exceed their Tolman-Oppenheimer-Volkoff mass. The gravitational-wave emission from such remnants provides a unique opportunity to measure the nuclear equation of state at densities and temperatures not available to terrestrial experiments. Current detector design is informed by gravitational-wave signals from general relativistic hydrodynamics simulations of neutron star mergers, typically with hybrid thermal treatments for the equation of state, where a cold equation of state is modified by adding a thermal component. We use realistic equations of state based on the relativistic mean field model with consistent treatment of thermal effects to compute the distribution of expected peak gravitational-wave frequencies. Marginalising over equation of state and progenitor neutron star masses, we show the peak frequency of emission ranges from $\sim2.5$ to 4 kHz. The width of this distribution suggests the need for broadband observatories with kHz sensitivity, and calls into question some of the so-called post-merger optimised configurations. We show the proposed KAGRA high-frequency design is well-suited to measuring post-merger remnants when compared to the KAGRA broadband design.

The increasing size and heterogeneity of exoplanet catalogs have made systematic habitability assessment challenging, particularly given the extreme scarcity of potentially habitable planets and the evolving nature of their labels. In this study, we explore the use of pool-based active learning to improve the efficiency of habitability classification under realistic observational constraints. We construct a unified dataset from the Habitable World Catalog and the NASA Exoplanet Archive and formulate habitability assessment as a binary classification problem. A supervised baseline based on gradient-boosted decision trees is established and optimized for recall in order to prioritize the identification of rare potentially habitable planets. This model is then embedded within an active learning framework, where uncertainty-based margin sampling is compared against random querying across multiple runs and labeling budgets. We find that active learning substantially reduces the number of labeled instances required to approach supervised performance, demonstrating clear gains in label efficiency. To connect these results to a practical astronomical use case, we aggregate predictions from independently trained active-learning models into an ensemble and use the resulting mean probabilities and uncertainties to rank planets originally labeled as non-habitable. This procedure identifies a single robust candidate for further study, illustrating how active learning can support conservative, uncertainty-aware prioritization of follow-up targets rather than speculative reclassification. Our results indicate that active learning provides a principled framework for guiding habitability studies in data regimes characterized by label imbalance, incomplete information, and limited observational resources.

Using the Spitzer and WISE images, we discovered 42 mid-IR luminous dusty supernovae with local integral-field spectroscopy data. The observed mid-IR emission indicates the presence of newly formed dust, or pre-existing dust heated by the radiation from the supernovae or circumstellar medium interactions. We carried out a systematic analysis of the supernova host environments and their dust properties, for understanding the dust-veiled exploding stars, and whether such an intense dust production process is associated with their local environments. We find that dusty supernovae prefer the locations with higher EW(H{\alpha}), lower metallicity, and heavier host extinctions compared to typical SN types, and they show the same increasing sequence in the values of EW(H{\alpha}) and oxygen abundance from hydrogen-rich, type IIn and hydrogen-poor dusty supernovae. These differences in environmental properties of different dusty SN types indicate the diversity of their progenitors. We also found that one marginal correlation is a negative correlation between the SN dust mass and star formation rate. This means that SNe would be more mid-IR luminous and more dust-rich at the region with lower star formation rate. However, the SN dust mass show no correlation with the metallicity and the host extinction, which were thought to be key factors affecting the mass-loss history of progenitors and the CSM environment of SNe. Therefore, the dust formation process in SNe might be insensitive to metallicity and the dust condition of their host environments.

We present a comprehensive search for non-thermal high-energy $\gamma$-ray emission from the nearby merging galaxy cluster Abell~119 ($z=0.044$) using archival observations spanning over seven decades in photon energy. Our analysis combines hard X-ray data from INTEGRAL/ISGRI (30--100~keV), MeV $\gamma$-ray observations from COMPTEL (0.75--30~MeV), and GeV--TeV data from the DArk Matter Particle Explorer (DAMPE; 3~GeV--1~TeV). No statistically significant emission is detected at the cluster position in any energy band. In the hard X-ray regime, we derive a $3\sigma$ upper limit of $F_{30-100\,\mathrm{keV}} \lesssim 8.5 \times 10^{-11}$~erg~cm$^{-2}$~s$^{-1}$ from ISGRI mosaic imaging. Reanalysis of archival COMPTEL data yields 95\% confidence-level upper limits ranging from $\sim 9 \times 10^{-11}$ to $\sim 8 \times 10^{-10}$~erg~cm$^{-2}$~s$^{-1}$ across 0.75--30~MeV. In the GeV--TeV range, DAMPE constrains the differential energy flux to $\sim 10^{-12}$--$10^{-10}$~erg~cm$^{-2}$~s$^{-1}$ (95\% confidence level). These results provide independent multi-band constraints on the reported GeV excess from recent Fermi-LAT studies. While our DAMPE limits do not exclude the flux levels claimed in those analyses, the absence of confirmation across keV--MeV--GeV bands indicate that any non-thermal emission from Abell~119 remains tentative.

Astrophotonics will be central to astronomical instrumentation, enabling lightweight, compact, and environmentally stable systems for both ground-based observatories and space missions. One key application is beam combination for nulling integrated with a photonic lantern, and long baseline nulling interferometry, which suppresses starlight to reveal exoplanets and companions. Compact, broadband photonic beam combiners are essential for providing a pathway to complex circuitry on a single chip and scalable solutions for single and multi-telescope instruments, and are investigated herein. Two-waveguide photonic combiners rely on symmetric coupling to interfere light, which is chromatic and requires modification for broadband operation. A three-waveguide configuration (tri-coupler) offers the potential for deeper, broader, and stable achromatic nulls. This work compares simulations of two evanescent tri-couplers and a multimode interference coupler (MMI) across the 1.5 - 1.8 ${\mu}$m band, evaluating exoplanet throughput, starlight attenuation, sensing characteristics, and estimations on fabrication tolerance. The standard tri-coupler was outperformed by both a tapered tri-coupler and the MMI, each of which achieved exoplanet throughput >85% throughout the band. The standard design recorded a minimum exoplanet throughput of 50% at the waveband's extremes. The tapered tri-coupler was further redesigned to achieve a non-degenerate sensing state. The MMI, while limited to a starlight attenuation of 40 dB $\left(10^{-4}\right)$ by uncoupled light, showed the greatest tolerance to fabrication errors, offering strong practical potential. Future designs aim to combine high exoplanet throughput, deep starlight attenuation, and non-degenerate sensing within a single integrated architecture.

We investigate the impact of a Coronal Mass Ejection (CME) on the transport and acceleration of relativistic protons in the solar wind using a coupled 3D Magnetohydrodynamics (MHD) simulation and a test-particle approach. The CME is driven by a spheromak injected into a Parker solar wind at a heliocentric distance of 0.139 AU. The trajectories of 5 GeV protons, injected towards the CME from 3 AU, are integrated in the guiding-centre approximation limit and scattered in velocity space with a mean free path $\lambda_{\|}$. Our results show that the CME can increment the protons energy by several GeV. The acceleration occurs during the time particles stream along the portion of a magnetic field line subject to compression downstream of the quasi-perpendicular portion of the CME-driven shock. In our configuration, the maximum energy gain, which is of the order of a few percent per shock crossing, occurs when the shock approaches 0.3 AU. Large energy gains require multiple passes through the acceleration region, which is made possible by the combined action of the mirror force and pitch angle scattering. The efficiency of the acceleration on time scales of the order of hours scales as $\lambda_{\|}^{-3/2}$. Energy spectra harden for decreasing parallel mean free path $\lambda_{\|}$.

We present an investigation of the full width at half maximum (FWHM, or {\gamma}) of absorption features of Type Ia supernova (SNe Ia). We found that, the average value of FWHM can be well predicted with the rest wavelength ({\lambda}). The velocity also plays an important role, as objects with a higher velocity tend to have a larger FWHM. Temperature may be the third factor, as we found that, at the same velocity (but different phases), a normal-velocity (NV) object tends to have a larger FWHM than high-velocity (HV) object. Also, 1991T/1999aa-like objects that are believed to have relatively high temperatures show the largest FWHMs if compared at the same velocity. Generally speaking, FWHM evolves very slowly with time and shows no correlation with {\Delta}m15, but 1991T/1999aa-like objects are characterized by relatively fast decreasing FWHM. On the other hand, we found that, objects with relatively small FWHMs shows a tighter correlation between absorption depth (A) and {\Delta}m15, possibly a sign of higher degree of homogeneity. We also found that A/{\gamma} of Si II {\lambda}5972 has a strong correlation with {\Delta}m15, and more importantly, a relatively slow time evolution, making it a useful luminosity estimator even in the absence of phase information.

Pulse-profile modeling of rotation-powered millisecond pulsars targeted by NICER has enabled mass--radius constraints of several neutron star sources, with implications for the dense-matter equation of state. For the bright isolated pulsar PSR J0030+0451, the inferred mass--radius was previously found to depend strongly on the assumed hot spot model. These hot-spot models yielded different mass--radius constraints, with the statistically preferred model exhibiting some mild tension with results inferred for PSR J0437$-$4715, PSR~J0614$-$3329, and GW170817. We present an updated pulse-profile analysis of PSR J0030+0451 using new NICER observations obtained between 2017 July to 2023 January, increasing the number of X-ray counts by about 50% compared to previous analyses. We jointly analyze the NICER data with archival XMM-Newton observations to better constrain the source spectrum and background. The new analysis significantly reduces the discrepancy between the hot spot models. The inferred mass and radius are $M = 1.43^{+0.20}_{-0.17}\,M_\odot$ and $R_{\rm eq} = 12.68^{+1.31}_{-1.04}$ km (68% credible intervals), reducing the tension with the results from other sources. In addition, the inferred hot spot configurations suggest the presence of intra-spot temperature gradients.

The circumgalactic medium (CGM) surrounding galaxies is believed to be a significant reservoir of baryons, yet its total mass remains poorly constrained. We present a novel approach to probe the CGM of the Andromeda galaxy (M31) using fast radio bursts (FRBs) from the CHIME/FRB Catalog 2. By comparing the dispersion measures (DMs) of 171 FRBs whose sightlines intersect M31's halo (within $r_{\rm vir}= 302\,\mathrm{kpc}$) to a control sample of 684 FRBs, we estimate the DM contribution from M31's CGM. We find evidence for an excess DM of $\delta\mathrm{DM} = 5.9$--$59.6\,\mathrm{pc\,cm^{-3}}$ in the inner halo ($\approx 0$--$151\,\mathrm{kpc}$) and $\delta\mathrm{DM} = 25.7$--$64.6\,\mathrm{pc\,cm^{-3}}$ in the outer halo ($\approx 151$--$302\,\mathrm{kpc}$). Using a generalized halo model parameterized by a closure radius $r_{\mathrm{close}}$, we constrain the baryon distribution and infer a best-fit value for $r_{\mathrm{close}}$ to be $9.2^{+9.9}_{-4.9}\,r_{\rm vir}$ ($1\sigma$), corresponding to a total CGM mass of $M_{b,\mathrm{halo}} = 18.6^{+7.9}_{-8.4} \times 10^{10}\,M_\odot$. Our results suggest that M31 may harbor a substantial fraction of its cosmic baryon budget in diffuse, ionized gas. This work demonstrates the potential of FRBs as a powerful tool for studying the CGM of nearby galaxies, with future larger samples expected to provide tighter constraints on the baryon content of galactic halos.

Accurate spatially resolved dust corrections are critical for interpreting the structure and evolution of star-forming galaxies (SFGs). We present an empirical model for predicting spatially resolved dust attenuation ($A_V$) in SFGs using integral field spectroscopy from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. Using a sample of 5,155 galaxies over $7.20<M_\ast<11.14$ and $0.0002 < z < 0.1444$, we derive $A_V$ maps from the Balmer decrement across more than 1,898,954 star-forming spaxels. Using local star formation rate surface density ($\Sigma_{\text{SFR}}$) as a predictor, the model achieves $R^2 = 0.69$ and RMSE $=0.22$ mag, with residuals that are approximately Gaussian and centred near zero. It predicts $A_V$ within a factor of $\sim$1.3 on kpc scales. We also demonstrate that the relation can be applied iteratively to recover dust-corrected $\Sigma_{\mathrm{SFR}}$ from uncorrected values, converging by the fourth iteration with minimal residual bias ($-0.01$ mag) and low RMSE ($0.42$ mag). The model accurately reproduces $A_V$ maps across diverse morphologies and orientations, including edge-on systems. It also recovers the observed radial $A_V$ profiles, capturing their dependence on stellar mass and relative star formation activity, with more massive and more strongly star-forming galaxies showing steeper gradients.

Stellar winds interact with the partially ionized interstellar medium (ISM), forming astrospheres. A key feature of this interaction is the hydrogen wall - secondary interstellar atoms produced via charge exchange near the tangential discontinuity separating the stellar wind from the ionized ISM component. This secondary component is decelerated and heated compared to primary interstellar hydrogen, making the hydrogen wall detectable through Lyman-$\alpha$ absorption spectra toward nearby stars. Such structures have been observed by the Hubble Space Telescope around the Sun and other stars. In this paper, we propose that another feature of the stellar wind/partially ionized ISM interaction may also be detectable in Lyman-$\alpha$ spectra: the neutral stellar wind. It forms via charge exchange between supersonic stellar wind protons and interstellar atoms penetrating deep into the astrosphere due to their long mean free paths. We present a parametric numerical analysis of astrospheric structures and their synthetic Lyman-$\alpha$ absorption spectra. Using a 2D kinetic-hydrodynamic model, we vary the terminal wind velocity while maintaining constant dynamic pressure to keep the astrosphere size consistent. For winds slower than the solar wind (terminal velocities $V_0 \lesssim 200$ km/s), charge exchange efficiency in the supersonic region increases dramatically, producing a distinct and observable absorption feature from the neutral wind. This signature is negligible for solar-like winds ($V_0 \approx 400$ km/s) but emerges as a direct spectroscopic diagnostic for winds up to $\sim 200$ km/s. Detecting this neutral wind absorption offers a novel method to directly constrain stellar wind velocities.

When a dark matter halo falls into a more massive object and becomes a subhalo, it typically loses much of its mass through tidal stripping. The reverse process is also possible in principle. The subhalo may gravitationally capture material from its host. If sufficiently efficient, this process could make an initially starless subhalo visible. We use high-resolution N-body simulations to estimate the efficiency of capture. We find that after an extended period orbiting within its host, at most $\sim 10^{-4}$ of a subhalo's remaining mass has been acquired since infall. This captured material is less concentrated to subhalo centre than material retained from before infall. It is also very much less abundant than host material that is instantaneously passing through the subhalo on almost unperturbed orbits. Captured stars are not sufficiently spatially concentrated to be distinguished from the dominant background of "field" stars, and their concentration in velocity space is no greater than that of typical stellar streams in the halo. Unfortunately, stellar capture is not efficient enough to allow initially starless low-mass subhaloes to be detected.

The early evolution of unstable hydrodynamic and magnetohydrodynamic equilibria is often governed by a few dominant linear eigenmodes. We investigate whether initialising a simulation with a superposition of linear eigenmodes that contains the most unstable mode saves computation time, and how the selection of the included modes affects the non-linear evolution. Using the non-linear (magneto)hydrodynamic simulation code MPI-AMRVAC, the evolutions of a flow-sheared fluid interface, a Harris current sheet, and a flow-sheared plasma interface were simulated for various initial perturbations. The novel initial perturbations were linear eigenmodes of the equilibrium, or superpositions thereof, and calculated with the Legolas code. We benchmarked to initialisation with velocity noise and, in the case of the Harris sheet, initialisation with an analytic magnetic field perturbation. By perturbing an unstable equilibrium with a superposition of linear eigenmodes that contains the most unstable mode, significantly less computation time is spent on the linear stage of the evolution compared to traditional perturbation methods. In the best case, the simulation time needed to reach the non-linear stage is reduced by an order of magnitude. The inclusion or omission of certain modes in the initial perturbation is observed to affect the non-linear evolution to various degrees. The perturbation of equilibria with a superposition of linear eigenmodes that contains the most unstable mode allows simulations to reach a late-evolution stage faster, thus saving computation time. Additional eigenmodes can be included alongside the fastest growing mode to obtain further benefits, for example, to accelerate symmetry breaking in the non-linear stage, or to delay their effect on the non-linear evolution. Coupling spectroscopic codes with (magneto)hydrodynamic codes therefore offers significant advantages.

Observations of linear structure connecting pulsars to gamma-ray halos reveal injection of TeV electrons into the interstellar medium (ISM). In some cases, this could be attributed to nearly scattering-free electron transport along large-scale magnetic fields connected to pulsar winds with very slow diffusion across the field lines. In this work we model this process with a magnetic flux tube emerging from the pulsar and attached to the ISM magnetic field. We show that particles in this case have an anisotropic distribution of magnetic pitch angle, such that the overall emission is highly beamed. We apply this model to pulsar tails and filaments, including the extended X-ray and TeV emission associated with PSR J1740+1000 and the misaligned X-ray jet in the Guitar Nebula, to constrain their particle population and magnetic fields.

We present a timing and scintillation study of the young Galactic halo pulsar PSR J1740+1000 using observations from the Nanshan, FAST, and Parkes radio telescopes. From timing analysis, we measure the pulsar's proper motion for the first time, indicating motion away from the Galactic plane at a position angle of 16.7 +/- 4.8 degrees (Galactic coordinates), with a total proper motion of 56.9 +/- 8.0 mas/yr and a corresponding transverse velocity of 329 +/- 80 km/s. This velocity suggests that PSR J1740+1000 is a typical-velocity young pulsar born within the Galactic halo. In scintillation studies, we detect scintillation arcs, arclets, and double-layered adjacent arcs in the secondary spectra. Under isotropic and anisotropic scattering assumptions, the screen-to-pulsar distance is 370 +/- 72 pc and 1 +/- 12 pc, respectively. The latter closely matches the scale of the pulsar wind nebula associated with PSR J1740+1000 and provides a better fit, suggesting that scattering is likely dominated by the PWN. The double-layered adjacent arcs observed on MJD 60180 imply that the pulsar's scattered image consists of two dominant components (A and B) and multiple weaker components. Component A is located at the pulsar's geometric position (0 uas), while Component B is located 112 +/- 16 uas and 23 +/- 17 uas from the central component under isotropic and anisotropic scattering, respectively. The frequency-independent angular position of Component B hints at refraction by an AU-scale structure within the scattering region, possibly originating from the PWN.

The Hubble tension and the recently reported anomaly in data from the Dark Energy Spectroscopic Instrument (DESI) are considered to pose serious challenges to the standard $\Lambda$CDM model of cosmology. In this work, we show that resolving the Hubble tension with a scenario featuring dark radiation-matter decoupling (DRMD) predicts the presence of dark acoustic oscillations (DAO) similar in scale to baryon acoustic oscillations (BAO). Using an inference independent of large-scale structure data, relying only on Planck measurements of the cosmic microwave background and SH$0$ES-calibrated supernova data, we find evidence for a DAO signal with drag-horizon scale $r_{d,\mathrm{DAO}} \in[54,65]\,\mathrm{Mpc}/h$ ($68\%\,\mathrm{C.I.}$) and amplitude $A_\mathrm{DAO} \in [0.02,0.05]$ ($68\%\,\mathrm{C.I.}$). These predictions provide a concrete target for current and upcoming large-scale structure surveys, including DESI, Euclid, and the Roman Space Telescope. Remarkably, the predicted DAO properties are consistent with those required to explain the DESI anomaly, offering both an alternative to evolving dark energy and a preliminary validation of the relevance of a dark radiation-matter decoupling scenario for addressing the Hubble tension.

Over the past decade, runaway stars have been identified, believed to originate either as surviving donors of Type Ia supernovae or as partially deflagrated accretors producing Type Iax supernovae. While the former have been extensively studied recently, the origins of the latter (also called LP 40-365 type stars) remain under-explored and therefore less well understood. So far seven such objects are known. In this paper, we report the discovery of a new LP 40-365 type runaway star, notably hotter than previously studied members of this class. Spectral analysis confirms that its atmosphere is neon- and oxygen-dominated, consistent with earlier analyses of other LP 40-365 type stars. Kinematic analysis indicates that the star has a high probability of being unbound from the Galaxy and was most likely ejected from the Galactic disk approximately 2.8 Myr ago with an ejection velocity exceeding 600 km/s. This result further emphasizes the discrepancy between the abundance yields and kick velocities predicted by white dwarf deflagration models and those observed in stars of LP 40-365 type, underscoring the need for a reassessment of such systems.

The Milky Way's nuclear star cluster (NSC) is a unique laboratory to study the formation and evolution of dense stellar systems around a supermassive black hole. Previous work suggests that most stars in the NSC are old; however, the detailed age and metallicity distributions remain uncertain. We constrain the star formation history (SFH) and metallicity of a poorly explored region located $\sim$3 pc from SagittariusA*. We analyse VLT/NACO imaging in an intermediate-band filter centred at 2.24 $\mu$m, complemented by $H$-band data. We construct completeness-corrected $K$-band luminosity functions (LFs), clearly identifying the Red Clump and Red Giant Branch Bumps. The SFH is derived by fitting cumulative LFs with MIST, PARSEC, and BaSTI models spanning a wide range of ages and metallicities, using Monte Carlo sampling to estimate uncertainties. Metallicity constraints are refined using spectroscopic measurements from the literature. The stellar population is predominantly old and metal-rich: $75.6 \pm 9.5$% of the stellar mass formed $\gtrsim 10$ Gyr ago, with median [M/H] $\sim +0.35$. An intermediate-age component at 2-3 Gyr contributes $20.8 \pm 8.7$%, while minor populations are present at $\sim$400 Myr ($0.9 \pm 0.8$%) and 20 Myr ($3.6 \pm 1.4$%), the latter representing a small but non-negligible young population. Systematic uncertainties from stellar models, binning, photometric range, unresolved binaries, and filter choice are assessed. These results indicate early dominant formation, a significant 2-3 Gyr episode, and minor recent activity, consistent with spectroscopic measurements and with properties of the inner NSC and nuclear stellar disc.

We use the galaxy ultraviolet luminosity function measurements at $z=6-9$ to constrain modification to standard inflationary power spectrum. These observables are sensitive to the matter power spectrum which itself depends on inflationary initial conditions. We consider specific models where a bump or oscillatory features are introduced to the standard power law inflation spectrum. We find that the galaxy luminosity observables can probe such modifications at wavenumbers $0.5\lesssim k \lesssim 20$ Mpc$^{-1}$. We obtain upper limits on the amplitude of bump-like features at the mentioned wavenumbers. We obtain constraints which are similar to previous constraints on these models using measurements of optical depth of reionization. However, the galaxy luminosity functions are a more direct probe for these type of models and, therefore, can complement indirect constraints coming from measurements of IGM properties.

The ionization mechanisms of low-ionization nuclear emission-line regions (LINERs), which are common in the local Universe, have been debated for decades. Our nearest large neighbor, M31, is classified as a LINER based on its optical emission line properties within the central kpc. In this work, we present a detailed photoionization modeling of the circumnuclear ionized gas in M31, explicitly tailored to its well-constrained physical conditions, including the absence of ongoing star formation and a currently inactive active galactic nucleus (AGN). Using spatially resolved CFHT/SITELLE observations, we find that photoionization by hot, evolved low-mass stars distributed throughout the bulge can roughly reproduce the observed radial intensity profiles of H{\alpha}, H\b{eta}, and [NII]. However, these models fail to match the observed [OIII] emission, producing radial profiles and [O III]/H\b{eta} ratios that are significantly steeper than observed. This discrepancy indicates a deficit of high-energy ionizing photons in standard stellar photoionization models, even with extended ionizing sources. We explore whether this tension can be alleviated by invoking either a bulge-filling, low-density ionized medium surrounding a denser H{\alpha}-emitting disk, or enhanced AGN activity in the recent past. While both scenarios can partially increase the [O III] emission, neither provides a fully satisfactory explanation under physically plausible conditions. Together with our earlier results for M81, these findings underscore persistent challenges in explaining LINER-like emission solely through conventional photoionization mechanisms.

This work investigates symbolic regression (SR) as an interpretable alternative to black-box machine learning for the classification of stars, galaxies, and quasars in the Sloan Digital Sky Survey Data Release 17 (SDSS DR17). We conduct a systematic comparative study of four state-of-the-art SR frameworks: {\tt PySR}, Exhaustive Symbolic Regression ({\tt ESR}) with MDL-based selection, Physical Symbolic Optimization ({\tt PhySO}) using deep reinforcement learning, and Multi-View Symbolic Regression ({\tt MvSR}). By deriving compact analytic functions (complexity $\leq 10$) on a representative training subset and subsequently evaluating them via a 5-fold stratified cross-validation protocol on 100,000 spectroscopically confirmed objects, we map spectroscopic redshift ($z$) to continuous classification scores. Our results demonstrate that these low-complexity expressions achieve high predictive reliability, with {\tt MvSR} reaching a Cohen's Kappa of 0.8948 and {\tt PhySO} achieving exceptional parametric stability ($\sigma < 0.002$). We show that these models not only match the performance of traditional baselines but also provide a transparent, mathematically concise characterization of the astrophysical boundaries separating galactic and extragalactic populations.

Wide binaries (WBs) offer a unique opportunity to test gravity in the low-acceleration regime, where modifications such as Milgromian dynamics (MOND) predict measurable deviations from Newtonian gravity. We construct a rigorous framework for conducting the wide binary test (WBT), emphasizing high quality sample selection, filtering of poor astrometric solutions, contamination mitigation, and uncertainty propagation. We show that undetected close binaries, chance alignments, and improper treatment of projection effects can mimic MOND-like signals. We introduce a checklist of best practices to identify and avoid these pitfalls. Applying this framework to Gaia DR3 data, we compile a high-purity sample of WBs within 130 pc with projected separations of 1 - 30 kAU, spanning the transition between the Newtonian and MOND regimes. We find that the scaled relative velocity distribution of wide binaries does not exhibit the 20% enhancement expected from MOND and is consistent with Newtonian gravity across all separations. A meta-analysis of previous WBTs shows that apparent MOND signals diminish as methodological rigour improves. We conclude that when stringent quality controls are applied, there is no observational evidence for MOND-induced velocity boosts in wide binaries. Our results place strong empirical constraints on modified gravity theories operating between a0/10 and 200 a0, where a0 is the MOND acceleration scale. Across this range of internal accelerations, Newtonian gravity is up to 1500x more likely than MOND for our cleanest sample.

Recent high-cadence observations by Subaru-HSC have identified a population of ultrashort-timescale microlensing events, providing a compelling window for planet-mass primordial black holes (PBHs) to constitute the entirety of dark matter. In this Letter, we demonstrate that this PBH population and the nanohertz stochastic gravitational-wave (GW) background reported by pulsar timing arrays (PTAs) can be naturally unified by a single primordial origin: a broad, nearly-flat enhancement of the curvature power spectrum with an amplitude of $O(10^{-2})$. The resulting PBH mass function spans the planet-to-solar mass range, while remaining consistent with all current observational constraints. This unified PBH--induced-GW framework makes concrete multi-messenger predictions, which can be decisively scrutinized by forthcoming microlensing surveys, next-generation PTAs, space-borne interferometers, precision astrometry, and laser ranging experiments.

We present observations of the classical T Tauri star DO Tau collected with the near-infrared SPIRou spectropolarimeter and precision velocimeter at the Canada-France-Hawaii Telescope from early 2020 to late 2025. Circularly polarized Zeeman signatures were clearly detected at most epochs in the atomic spectral lines of DO Tau, yielding longitudinal magnetic fields of up to 280 G modulated with a period of 5.128+-0.002 d which we identified as the rotation period of DO Tau. Applying Zeeman-Doppler imaging to the SPIRou data recorded in 2021, 2024 and 2025, we found that DO Tau hosts an unusual large-scale magnetic field that is weaker, less poloidal, more inclined to the rotation axis, and varies more rapidly with time than those of previously studied T Tauri stars, possibly as a result of intense accretion between the inner disk and the stellar surface. The dipole component of this large-scale field of about 0.2-0.3 kG even flipped polarity toward the end of our observing campaign, making DO Tau the first T Tauri star for which a magnetic polarity reversal is reported. The magnetospheric gap surrounding the central star was quite compact, extending to ~1.6 Rstar (0.014 au) as a result of the strong accretion rate (log Mdot = -7.7 Msun/yr), with the inner accretion disk being warped by the tilted stellar magnetic field. Radial velocity variations suggest the presence of a close-in planet of a few Mjup or a density structure in the inner accretion disk at an orbital period of 21 d (corresponding to 0.12 au), which might be linked to the wiggle in the jet axis of DO Tau.

The origin of the X-ray plateau phase observed in a large fraction of gamma-ray burst afterglows remains debated. We present a novel analytic framework for the hydrodynamics of ultra-relativistic, radially stratified outflows interacting with an external medium. By explicitly accounting for a continuous distribution of Lorentz factors within the ejecta, we derive analytic expressions describing the evolution of a long-lived, mildly relativistic reverse shock and determine its crossing time. Then, we compute the resulting synchrotron emission from both the forward and reverse shocks. The forward shock naturally produces a shallow, long-lasting X-ray decay consistent with the observed properties of X-ray plateaus (including the Dainotti relation). We further show that reproducing the observed plateau durations requires the stratified ejecta to extend to Lorentz factors of order $\gtrsim 100$, consistent with the ultra-relativistic outflow that powers the prompt $\gamma$-ray emission. The reverse shock generates a long-lived millimeter emission component that outshines the forward shock emission at these wavelengths. Both the plateau and reverse shock emission terminate smoothly once the slowest ejecta are processed, marking a transition to the standard Blandford-McKee self-similar evolution without requiring late-time energy injection or an additional emission component. Such stratified outflows are expected on physical grounds, as the ultra-relativistic ejecta responsible for the prompt $\gamma$-ray emission are unlikely to be launched with a single Lorentz factor. This model provides a unified picture in which the same outflow powers the prompt emission, the X-ray plateau, and the subsequent afterglow evolution.

Relativistic jets from Active Galactic Nuclei (AGN) are highly energetic and emit radiation across a wide range of frequencies. Despite several observational studies, their particle composition still remains a key open question. The detection of high-energy neutrinos from blazar sources such as TXS 0506+056 has highlighted the plausibility of hadronic/lepto-hadronic models for AGN jets. To understand the origin of high-energy neutrinos from such sources, it is imperative to capture the complex interplay between the jet dynamics, their composition, and the mechanism of particle acceleration and cooling in relativistic jets. In this pilot study, we have coupled a numerical multi-zone framework for lepto-hadronic modeling, with 3D relativistic magneto-hydrodynamic simulations of AGN jets, including external photon fields. Our framework provides synthetic multi-wavelength and neutrino flux by spatially sampling the simulated jet into multiple zones. We investigate the implications of such a framework in exploring the different intrinsic and extrinsic pathways for proton-enrichment in jets. Essentially, we find that for low proton-to-electron number density ratios, producing a substantial jet neutrino flux, requires the underlying proton energy distribution to have a relatively flat spectrum with a power-law index of $\simeq 2.0$. We further find that while intrinsic shocks triggered by kink-instabilities in the jet can accelerate electrons to high energies, they may not be sufficient to produce such flat particle energy distributions for the chosen set of parsec-scale jet parameters. Finally, to produce a significant jet neutrino emission, our simulations suggest the need to consider particle acceleration mechanisms through alternative pathways, either internal or external.

The reionization process is driven by ionizing photons from dwarf galaxies in halos with virial temperature $T_{\rm vir} \gtrsim 10^4$ K, while minihalos whose $T_{\rm vir}\lesssim 10^4$ K consume ionizing photons and have negative contributions to reionization. Since ionizing sources and minihalos have different clustering characteristics, not only the reionization history, but also the morphology of the ionization field, is sensitive to the small-scale power spectrum. If the power spectrum at small scales is enhanced compared with the standard six-parameter $\Lambda$CDM model, then both the sources and sinks of ionizing photons would be boosted and the net impact depends on the competition between them. Therefore, the 21 cm signal that can probe the morphology of the ionization field will be a useful tool for detecting the small-scale power spectrum. Using the power spectrum proposed by Cielo et al. (2025) (C25) as a demonstration, we investigate the influence of small-scale power spectrum enhancement on the ionization field and the 21 cm signal. We find that for the C25 model, even under the constraints of observed UV luminosity functions for high-$z$ galaxies and reionization history, the 21~cm power spectrum and the bubble size distribution could be still significantly different from the regular $\Lambda$CDM model. The upcoming SKA-low AA* telescope, and a further imaging telescope, have the potential to detect the small-scale power spectrum more deeply.

The evolution of the global solar magnetic field directly impacts the interplanetary magnetic field (IMF). During the solar maximum of Cycle 24, the monthly averaged IMF strength doubled over five Carrington rotations in late 2014. To understand the physical origin of this increase, we investigate the temporal evolution of open magnetic flux resulting from the emergence and decay of bipolar magnetic regions (BMRs). Using surface flux transport and potential field source surface models, we simulated how BMR characteristics, spatial distributions, and interaction with background magnetic fields affect open flux evolution. Our simulation confirmed that the relative configuration of BMRs can either inhibit open flux expansion via closed loops or promote it through favorable connections. The increase in open flux is primarily driven by the equatorial dipole component, which is enhanced by differential rotation acting on tilted BMRs. These behaviors suggest that large open field structures develop from equatorial dipole components formed by these stretched BMRs. We attribute the rapid IMF increase in 2014 (Carrington rotations 2152-2157) to the combination of the following three factors: (1) a specific sunspot configuration that facilitated the expansion of the southern coronal hole; (2) the emergence of a giant sunspot group (active region 12192) with high magnetic intensity; and (3) the diffusion of these regions, which reinforced the global magnetic field. These results imply that rapid open flux variations during solar maximum are governed not only by the characteristics of emerging BMRs but also by their interaction with pre-existing large coronal holes.

Across the peak of cosmic star formation at $z\sim1-2$, inflow, processing, and feedback drive rapid changes in the spatial distribution and chemical composition of baryons in galaxies and surrounding reservoirs; this baryon cycle can be tomographically mapped by line intensity mapping (LIM) of atomic hydrogen, ionized carbon, and carbon monoxide. We present a simulation-based forecasting framework for detecting auto- and cross-power spectra between spectroscopic surveys of four such tracers at $z\sim0.5-1.7$ mapping the same deep field - TIM, EoRSpec/FYST, MeerKAT, & Euclid. We forward-model 3-D distributions for these tracers from magnetohydrodynamic simulations, directly capturing the two-halo, one-halo, and shot statistics without relying on analytical decompositions. We further detail a signal-to-noise formalism, tailored to LIM surveys with highly anisotropic geometries and Fourier-space coverage. We demonstrate that galaxy cross-correlations will be the dominant discovery channel for current-generation surveys. These instruments will detect the auto-spectra for CO and HI 21 cm and the CO $\times$ 21 cm cross-spectrum at modest S/N $\sim 1-10$, while placing upper limits on the [CII]-158$\mu$m signals. [CII], CO, and HI LIM will be $\sim3-30\times$ ($0.5-1.5$ dex) more sensitive to cross-correlation with the Euclid survey, however, than their respective auto-correlations, constraining all three models of line emission at high significance (S/N $\sim 10-40$) within this decade. Finally, we formulate a staged instrumental trajectory with planned or reasonable improvements, including the as-proposed SKA-Mid. We forecast advancing the per-$k$-mode sensitivities of each auto-, galaxy-line, and line-line spectrum by several orders of magnitude, enabling new percent- and sub-percent level constraints on cosmology and the redshift evolution of star formation and the baryon cycle.

Going beyond second order in weak lensing (WL) statistics is known to break degeneracies among cosmological parameters. We take a step further here, investigating whether higher-order statistics (HOS) in weak lensing can disentangle among General Relativity (GR) and modified gravity (MG), also taking into account the presence of massive neutrinos. To this end, we rely on mock convergence maps obtained from GR and $f(R)$ gravity N - body simulations, and we look for MG signatures in a wide set of higher-order WL probes. We rely on different metrics to quantify the discriminatory power of each probe, also varying the measurement setup. We find out that WL HOS can indeed disentangle MG and GR also in the presence of massive neutrinos.

We present a detailed study of HUDFJ0332.4-2746.6 (HUDF46), a $z \approx 1.84$ overdensity in the Hubble Ultra Deep Field, previously identified with HST as a proto-cluster. JWST/NIRISS spectroscopy expands its confirmed membership from 18 to 37 galaxies, while deep HST/ACS, JWST/NIRCam, and JWST/MIRI imaging provide a comprehensive multiwavelength view from the rest-frame UV to the mid-infrared. This dataset probes the population across three dex in stellar mass ($M_\bigstar \approx 10^{7.5\text{--}10.5}\,M_\odot$), delivering the first direct view of a young cluster down to such low-$M_\bigstar$ at $z\gtrsim1$. Assuming virialization, we derive a velocity dispersion of $\sigma \approx 670\pm 91\,\mathrm{km\,s^{-1}}$ and a halo mass of $M_{200} \approx (1.2\pm0.2) \times 10^{14}\,M_\odot$, in agreement with X-ray constraints from deep {\it Chandra} data. Despite residing in a massive halo likely in the hot-halo regime, the population is overwhelmingly star-forming, with no established red sequence and no extended X-ray emission from a hot intracluster medium. HUDF46 members have stellar and structural properties nearly indistinguishable from coeval field galaxies, and the structure hosts only one AGN candidate, found in its brightest galaxy, which lies at the cluster center. Overall, HUDF46 appears to be in a transitional phase prior to the onset of environmental quenching, making its galaxy population a key benchmark for tracing the processes that will later build a passive population and shape the assembly of massive clusters at later cosmic times.

When in situ solar energetic electron (SEE) events are closely associated with nonthermal flares, the escaping electron population is frequently observed to be much smaller than the nonthermal-radiation-emitting population near the solar surface. If a single accelerated population drives both signatures, the physical mechanism causing this severe deficit of upward-propagating electrons remains poorly this http URL on one of the 2022 November 10--12 SEE events associated with recurrent solar jets and interplanetary type III radio bursts, we present a new, combined microwave--X-ray analysis using the Expanded Owens Valley Solar Array (EOVSA) and the Spectrometer/Telescope for Imaging X-rays (STIX) aboard Solar Orbiter. This synergy enables, for the first time for such an event, spatially resolved diagnostics over a broad energy spectrum of the near-Sun energetic electrons, complemented by in situ measurements made by spacecraft at multiple heliocentric longitudes and distances. Consistent with earlier results based on in situ and X-ray data, our results show that only 0.1--1\% of energetic electrons escape into interplanetary space. Crucially, the new microwave spectral imaging analysis suggests that energetic electrons are strongly concentrated in a compact region just above a mini-flare arcade at the base of the jet spire, and that their number density decreases by at least two orders of magnitude in the direction of the jet spire away from this region. This steep gradient, revealed by the microwave diagnostics, points to efficient local acceleration and trapping in the region analogous to the above-the-looptop ``magnetic bottle'' region in major eruptive flares, allowing only a small fraction of electrons to access open magnetic field lines and enter interplanetary space.

Using the POLISH2 polarimeter at the Lick Observatory Shane 3-m and Nickel 1-m telescopes, we discover rotation phase-locked variations in linear polarization to be common among asteroids and a NEO in a clear, 383 to 720 nm bandpass. Essentially all bodies in our eight-year study harbor statistically significant, repeatable linear polarization variations at the $0.01\% = 100$ ppm level or above (1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, 6 Hebe, 7 Iris, 12 Victoria, 15 Eunomia, 16 Psyche, 132 Aethra, 216 Kleopatra, and 65803 Didymos). Since polarimetry is a differential technique, such variations cannot be due to shape changes and must be caused by heterogeneity in surface albedo and/or composition. While (4) Vesta has long been known to exhibit large, repeatable polarization variations across its surface, we discover the variations on (6) Hebe, (12) Victoria, and (65803) Didymos to be 1.5 to 3.5 times as large. We hypothesize that the \textit{DART} impact with Dimorphos blanketed Didymos with depolarizing ejecta, which suggests pristine variations across Didymos to have been even larger. As the only NEO in this study with data quality sufficient to investigate polarization variations, Didymos' huge variations suggest they may be common among NEOs. We also discover optical circular polarization to be enhanced for low-albedo, M type asteroids, which is correlated with their large radar albedos. Thus, we present optical circular polarimetry as an alternative method for the identification of metalliferous bodies.

We forecast that the Terahertz Intensity Mapper (TIM) cross-correlated with Euclid's Fornax deep field (EDF-F), TIM$\times$EDF-F, will detect the [CII]-galaxy cross-power spectrum at a median redshift of 1.1 with $\gtrsim 7 \sigma$ confidence. The Poisson component of the cross-power spectrum at $0.1 \leq k \leq 10$ hMpc$^{-1}$ (i.e. cross-shot noise) will be detected at $\gtrsim 3 \sigma$ in 4 bins spanning $0.5 < z< 1.7$. This measurement will constrain the mean [CII] specific intensity over half of cosmic history and assess the degree to which Euclid-selected galaxies account for the [CII] intensity observed by TIM. We find that TIM can detect the cross-power spectrum across a wide range of [CII] intensity models.

Shock waves generated during cluster mergers offer a powerful probe of how large-scale structure grows and evolves in the Universe. As part of the MeerKAT-South Pole Telescope (SPT) survey, we report the discovery of a single arc-like radio relic in the galaxy cluster Abell 4067 ($z=0.099$), one of the lowest-mass clusters known to host such a structure. MeerKAT UHF-band (0.58--1.09 GHz) observations reveal a relic with a largest linear size of $\sim 1.48 \pm 0.02$ Mpc, located at a projected distance of 0.95 Mpc from the cluster centre. XMM-Newton X-ray data show that the relic's position and orientation relative to the intracluster medium (ICM) elongation are consistent with a merger-driven shock-wave scenario. The relic has an estimated radio power of $3.10 \pm 0.03 \times 10^{24}$ W Hz$^{-1}$ at 150 MHz. When placed in the $P_{150\,\mathrm{MHz}}$--$M_{500}$ scaling relation, the Abell 4067 relic appears less luminous compared to relics in more massive clusters, suggesting an association with weak merger shocks. This finding supports the idea that relics in low-mass clusters may form through less energetic merger events, leading to weak merger shocks. This is further supported by the absence of a detectable central radio halo in Abell 4067, reinforcing the idea that luminous radio halos are not a universal outcome of cluster mergers and highlighting the role of cluster mass, merger energetics, and evolutionary stage in shaping diffuse radio emission in the intracluster medium.

One key piece of evidence for dark matter is the rotation-curve problem: the disagreement between measured galactic rotation curves and their luminous mass. A novel solution to this problem is presented here, in a model that predicts observed Doppler-shifted spectra based only on the luminous matter estimates and one free model parameter. This model is applied to fit the rotation curves of the SPARC sample of 175 galaxies, yielding mass-to-light ratios, goodness of fit measurements, and the free parameter. The model's average reduced chi square compares favorably with the dark matter model for the same data, and more galaxies are successfully fit by this model. The model provides a useful formulation linking luminous matter to the observed rotation curves, with the dark matter contribution to galaxies encoded in two transformation terms of the luminous mass. It also offers a lower-parameter characterization of the rotation curve problem, and a power law relationship between the model's free parameter and galactic photometric quantities is observed, potentially removing the need for the free parameter.

We present a broadband study of the WLQ SDSS J101353.45+492758.1, which displays a nearly featureless UV-optical spectrum with only a weak Mg II line alongside an exceptionally low X-ray flux. We model its spectral energy distribution using the relativistic thin-disk model kerrbb with a power law, and the multicomponent AGN model relagn, a physically motivated extension of agnsed incorporating warm and hot Comptonizing regions. Our fits constrain the black hole mass, accretion rate, X-ray loudness, and coronal energetics. Both approaches yield consistent BH masses of M_{BH} \approx 2 \times 10^{9} M_\odot and an Eddington accretion rate of \dot m \approx 0.1. The relagn fit including a warm Comptonizing region provides a significantly improved representation of the UV-soft X-ray continuum. The warm corona, characterized by kTe \simeq 0.20 keV, {\Gamma} \simeq 3.8, and an optical depth {\tau} \simeq 7.26, extends to \sim 34 R_g. The hot corona appears compact and energetically suppressed, leading to an intrinsically weak X ray output with log(L_{X}/L_{bol}) \simeq -4.29, among the lowest reported for WLQs. The {\alpha}_{ox} \sim 2.06 indicates the source to be in high/soft AGN spectral state. The combination of a luminous, standard disk and extremely weak hot corona suggests that this quasar hosts a highly inefficient inner coronal region. This explains its X-ray faintness and extreme deficit of high-ionization emission lines. The source may represent an AGN analog in "ultrasoft" accretion state, or a system in which the ionizing continuum is suppressed by a compact or quenched corona. Our study suggests that the source is not accreting at high Eddington ratio, highlighting the physical diversity of WLQs, and supports the view that geometric and radiative effects jointly shape their extreme spectral properties.

Being the most massive known young stellar cluster in the Milky Way, Westerlund 1 (Wd1) constitutes an ideal benchmark for understanding the evolution of massive stars. However, the cluster age remains highly controversial (~4-10 Myr), hindering the use of Wd1 as a reference for massive star evolution. One of the main issues is high foreground extinction, which has so far prevented the detection of the main sequence. Using infrared spectroscopy we seek to detect the cluster's main sequence for the first time, to characterise the Hertzsprung-Russell diagram, and to use the cluster's turn-off to obtain a robust age estimate. We obtained multi-epoch, near-infrared VLT/KMOS spectroscopic observations of Wd1 to map its population of massive stars. The spectra of ~110 members were analysed with CMFGEN models to derive stellar parameters, populate the cluster Hertzsprung-Russell diagram, and compare it with isochrones from evolutionary models. Our observations returned 47 new spectroscopically identified cluster members, with spectral types O9-B1 III-V. The cluster turn-off indicates an age of 5.5+/-1.0 Myr at a distance of 4.23+0.23-0.21 kpc, displaying a moderate degree of coevality. We demonstrate that our estimate of the age of Wd1 is robust against reasonable changes in the distance and extinction law, and the adopted rotational velocity and metallicity of the stellar isochrones. We further find that ~65% of the OB stars with multi-epoch coverage exhibit radial-velocity variability. Infrared observations of the unevolved stellar population support a single episode of star formation with an age of ~5.5 Myr, reinforcing its potential as a benchmark for massive star evolution and providing a reference sample for future binary population studies.

We examined the long-term behavior of the superbias calibration frames for the Advanced Camera for Surveys Wide Field Channel (ACS/WFC) aboard the Hubble Space Telescope (HST). Superbias frames are used to remove detector-level bias structure from science images and are currently generated after an anneal and delivered monthly. The primary goal of this study was to determine whether the frequency of superbias generation could be reduced without compromising calibration quality, potentially aligning with the Wide Field Camera 3 UVIS (WFC3/UVIS) approach of generating only one superbias per year. We analyzed superbias frames produced from 2007 through 2024 to investigate whether these calibration products have changed significantly over time, and whether the frequency of superbias generation and delivery could be safely reduced without loss of calibration accuracy. In addition to visual inspections and pixel-level comparisons, we employed Principal Component Analysis (PCA) to evaluate whether any long-term, global structure exists beneath the apparent noise in these frames. Our findings show that the superbias structure has remained fairly stable post-Servicing Mission 4 (SM4), a 15-year period, and no significant or unexpected global trends or systematic shifts were detected. However, due to unstable hot columns and increasing readout dark observed in ACS/WFC data, it is likely that these calibrations still benefit from more frequent superbias updates than the annual cadence adopted for WFC3/UVIS.

The RADiation-hard Electron Monitor (RADEM) aboard the Jupiter Icy Moons Explorer (JUICE), launched on 14 April 2023, measures high-energy protons and electrons during the cruise phase and will continue throughout the nominal mission. Initial in-flight observations could not be explained by pre-flight ground calibration, motivating an in-flight calibration campaign. We calibrated the RADEM sensors using galactic cosmic rays by progressively increasing detector thresholds, thereby modifying their response to high-energy particles. Threshold-dependent in-flight count rates were compared with theoretical expectations derived from the Badhwar-O'Neill 2020 galactic cosmic ray model and corresponding response functions. These results were used to derive new in-flight calibration coefficients and to develop a flux reconstruction algorithm based on the bow-tie method. In several cases, the in-flight calibration slopes differ by up to an order of magnitude from ground calibration values. Proton fluxes from solar energetic particle events reconstructed with this method agree within a factor of two with independent measurements from the Solar and Heliospheric Observatory. These results demonstrate that RADEM provides accurate proton flux measurements in interplanetary space and is well suited for both single-spacecraft analyses and coordinated multi-mission studies of solar energetic particles. While electrons were clearly detected during the JUICE Lunar-Earth gravity assist, reliable reconstruction of their fluxes requires further analysis.

Intense solar eruptions occasionally accelerate solar energetic particles (SEPs) and can trigger ground-level enhancements (GLEs). Among the 77 known GLEs, the first four GLEs, #1 -- 4 in the 1940s took place before the advent of the standard neutron monitors and were missing from the International GLE Database. This data gap challenged their quantification. To overcome this difficulty, we systematically gathered, digitised, and quantified contemporaneous cosmic-ray records pertaining to these GLEs. These data allow us to reconstruct the temporal evolution, with the 1 -- 15 min resolutions, of these GLEs, and broaden their geographical coverage to a global scale. GLEs #1 and #3 exhibited gradual increases in their rise times, measured at 45 +/- 15 and 105 +/- 15 min, respectively. In contrast, GLEs #2 and #4 both exhibited abrupt increases of 15 +/- 15 min. We also compared integral ionisation increase on the standard ionisation chambers and their local geomagnetic cutoff rigidities Pc to qualitatively compare these GLE's spectral hardness: Our result indicates that their spectra are extremely hard for GLEs #2 and #4 and mildly hard for GLEs #1 and #3. GLE #3 showed the greatest integral ionisations for polar detectors among them.

We report the first unambiguous detection of refractive substructure in an active galactic nucleus (AGN) using ground-based Very Long Baseline Interferometry (VLBI). Our analysis of TXS 2005+403 - observed at 1--5 GHz along a line of sight through the Cygnus region - reveals clear signatures of turbulence-induced substructure on long baselines that cannot be explained by the smooth scatter-broadened profile from diffractive effects alone. This signal persists across multiple observations spanning 2010-2019, demonstrating stable scattering properties along the line of sight. The combination of high flux density, compact intrinsic structure, and strong scattering establishes TXS 2005+403 as an exceptional laboratory for probing Galactic turbulence. This detection demonstrates that AGNs can serve as cosmic lighthouses illuminating interstellar plasma across the sky, complementing pulsar scintillation studies and informing scattering mitigation for millimeter-wavelength imaging of Sagittarius A*.

Hypothetical light scalar particles trigger the superradiant instability around spinning black holes (BHs), causing clouds of scalars to grow around the BH. In the presence of sufficiently strong particle self-interactions (characterized by the decay constant $f$), scalars are ejected from BH orbits, resulting in coherent, non-relativistic emissions that continuously carry away the BH's angular momentum. Parameters exist for which cloud growth is much faster, and scalar depletion is much slower, than the age of the Galaxy. This defines a distinct class of astrophysical sources of scalars, which we call \emph{BH scalar sirens} -- BHs that persistently emit scalars effectively forever. We compute the scalar background from the expected population of $N_\text{BH}\sim 10^{8}$ isolated stellar-mass BHs in the Milky Way, which are sirens for scalars in the mass range $10^{-13}$--$10^{-11}\,$eV and $f\lesssim 10^{14}$--$10^{9}\,$GeV. This provides a detection target independent of early-universe scalar production or cosmological initial conditions. The generated observable signals are up to two orders-of-magnitude larger than those expected from a misaligned cosmic scalar in this mass range. The energy spectrum of emitted scalars is distinctly broader and at higher velocities (up to $\sim 10^{-1}c$) than that of virialized dark matter, and encodes the mass and spin distributions of the BH population. While stellar-mass Milky Way BHs are our primary target, our framework extends to supermassive, intermediate-mass and light BHs. Given the difficulty of directly observing populations of isolated BHs, scalar emissions offer a novel probe of these otherwise invisible objects, highlighting the potential for joint discovery between scalars and BHs, and broadly motivating searches for scalars over many orders-of-magnitude in mass.

The Blandford-Znajek (BZ) mechanism is widely recognised as the most compelling process to extract rotational energy from an accreting black hole and power the emission of relativistic jets. We explore the universality of this process for generic black-hole spacetimes within the Konoplya-Rezzolla-Zhidenko formalism and find that the lowest-order contribution to the BZ power is invariant across different black-hole spacetimes. We also show that at the next-leading-order, different black-hole spacetimes will lead to different BZ luminosities. As a result, while slowly rotating black holes cannot be distinguished via measurements of their jet power, rapidly rotating ones have the potential of providing information on the strong-field properties of the spacetime when independent measurements of the BZ luminosity and of the black-hole angular velocity are available.

Departure from idealised plane waves gives rise to intricate geometric structures in wave fields. One such structure is the polarisation singularity, which emerges when multiple monochromatic waves interfere (such as would be the case for stochastic backgrounds), producing loci of purely circular or linear polarisation. In this work, we extend the theory of polarisation singularities to gravitational waves and higher spin fields. Building on the electromagnetic description, we formulate the gravitational analogue of polarisation singularities and show that they are generic features of gravitational waves. Their dimension, however, depends on the spin of the field. We illustrate these results with simulations of plane-wave interference and analyse the resulting singularity densities.

Monte Carlo simulations were conducted using the Particle and Heavy Ion Transport code System (PHITS) to investigate the role of secondary cosmic rays in the generation of long-duration bursts from thunderclouds and to clarify the conditions of the electric field region responsible for particle acceleration. The simulations utilized realistic secondary cosmic-ray spectra, including gamma rays, electrons, positrons, and muons, as input. The simulation results indicate that gamma rays provide the dominant supply of seed electrons for long-duration bursts, regardless of the geometry or strength of the electric field region. They also reveal the structure and strength of the electric field region required to produce gamma rays exceeding several tens of MeV, which have so far been detected only by high-altitude observations. Furthermore, the fluxes of long-duration bursts estimated from the simulation results were compared with observational data to constrain the properties of the electric field region. In particular, the comparison with measurements at Yangbajing, located at an altitude of 4.3~km, helps narrow down the possible range of electric field strengths and configurations.

In Einstein-scalar-Maxwell theories with a coupling between the scalar field $\phi$ and the electromagnetic field strength $F$ of the form $\mu(\phi) F$, we investigate the existence of exotic compact objects (ECOs) and their observational signatures in photon and massive-particle dynamics. For $\mu(\phi)$ diverging at the origin while all physical quantities remain finite, we demonstrate the existence of electrically charged ECOs with a shell-like structure whose density peaks at an intermediate radius. We compute their mass and radius, together with the scalar and vector field profiles, on a static and spherically symmetric background. We then examine the existence of photon rings and place bounds on a model parameter by requiring the absence of a linearly stable photon ring. Under this condition, photon echoes from ECOs are absent. We also compute the gravitational-lensing deflection angle $\Psi$ and show that it attains a maximum for an impact parameter of the same order as the ECO radius. Finally, we study the parameter space in which innermost stable circular orbits of massive particles exist.

Using the ZeroDegree multi-reflection time-of-flight mass spectrograph of the CRISMASS project at RIKEN Radioactive Isotope Beam Factory, we performed high-precision mass measurements of proton-rich nuclei near the doubly magic nucleus 100Sn, achieving uncertainties on the order of 10 keV. The masses of 91Rh, 92Pd, and 96Cd were determined for the first time with high precision, and the accuracy of several additional masses was substantially improved. Incorporating the new data into X-ray burst simulations significantly reduces the abundance uncertainties in the $A$ = 90-100 region, shifting the reaction flow toward $A$ = 90 production and suppressing the synthesis of heavier nuclei. Further investigation of the $\nu p$-process indicates that 99Rh plays a significant role in the reaction flow within the mass region studied. These high-precision mass measurements refine the mass surface near 100Sn and provide critical constraints on models of proton-rich nucleosynthesis.

Neutrino non-standard interactions (NSI) with electrons, predicted in many extended theoretical models of particle physics, are known to alter the picture of neutrino decoupling from the cosmic plasma. We update previous analyses of neutrino decoupling in presence of NSI with electrons, extending the parameter space in order to provide, for the first time, a full study of their effect on the production of light elements during Big Bang Nucleosynthesis (BBN). We compare the BBN bounds on non-universal and flavour-changing NSI parameters with the constraints from terrestrial experiments. Our results show that the limits from BBN are significantly less stringent than the experimental bounds, but they are complementary and can provide a test of neutrino physics at different temperature scales and epochs.

The gravitational-wave (GW) detections reported by the LIGO-Virgo-KAGRA (LVK) collaboration have so far been consistent with quasi-circular compact binary coalescences (CBCs). Nevertheless, a small fraction of binaries driven to merge through dynamical interactions in dense stellar environments or in field triples may retain measurable orbital eccentricity when entering the sensitive frequency band of LVK detectors. Confident measurement of eccentricity in the LVK band would provide strong evidence for such dynamically driven mergers; however, eccentric gravitational-waveform models are computationally expensive, and performing production-level inference on all detected signals is not an efficient use of resources when eccentric signals are expected to be rare. An intermediate step between detection and analysis, in which the signal is assessed for the potential presence of eccentricity, could provide quick recommendations for which signals should undergo full eccentric inference. We apply the wavelet scattering transform (WST) to a large set of synthetic waveforms in realistic noise and assess its discriminatory power using simple linear and shallow neural-network classifiers. We find that the WST representation enables effective discrimination between eccentric and quasi-circular binaries and provides a compact multi-scale representation of GW signals. Our approach achieves ~64% percent detection accuracy at a false alarm rate of 10%, with an AUC of 0.844 and an average precision of 0.876. We also examine the ability of our classifiers to distinguish eccentricity from spin-induced precession and find robust performance across a range of spin-precession magnitudes.

This workshop report from "Shaping the Digital Future of ErUM Research: Sustainability & Ethics" (Aachen, 2025) reviews progress on sustainability measures in data-intensive ErUM-Data research since the 2023 call-to-action on resource-aware research. It evaluates short-, medium-, and long-term actions around monitoring and reducing CO2 emissions, improving data and software FAIRness, optimizing workflows and computing infrastructures, and aligning operations with low-carbon energy availability, including concepts such as "breathing" computing centers, long-term data storage strategies, and software efficiency certification. The report stresses the need for systematic teaching, training, mentoring, and new support formats to establish sustainable coding and computing practices, particularly among students and early-career researchers, and highlights the importance of dedicated steering and funding instruments to embed sustainability in project planning. Ethical discussions focus on the transformative use of AI in ErUM-Data, addressing autonomy, bias, transparency, explainability, attribution of responsibility, and the risk of deskilling, while reaffirming that accountability for scientific outcomes remains with human researchers. Finally, the report emphasizes that sustainable transformation requires not only technical measures but also targeted awareness-building, communication strategies, incentives, and community-driven initiatives to move from awareness to action and to integrate sustainability and ethics into everyday scientific practice.

We present a class of models based on a gauged $U(1)_F$ flavor symmetry that explains the hierarchical structure of fermion masses and mixings via the Froggatt-Nielsen (FN) mechanism, while also solving the strong CP problem by the Peccei-Quinn (PQ) mechanism. A global $U(1)_{\rm PQ}$ symmetry with a nonzero QCD anomaly emerges accidentally in this setup as a byproduct of the gauged $U(1)_F$ symmetry. The resulting axion is shown to be of high quality, with the axion potential safeguarded against quantum gravity corrections by the gauge symmetry. Three models, which are generalizations of the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) axion model, are presented realizing this idea. The right-handed neutrino mass scale is identified as the Froggatt-Nielsen scale in these models. We present explicit UV completions of the FN sectors of these models and show that they preserve the high quality of the axion. In these models, the axion acts as a flavon field, leading to testable predictions in flavor-changing decays of neutral mesons. The axion also serves as the dark matter of the universe with the right amount of relic abundance without causing cosmological domain wall problems. Baryon asymmetry of the universe is realized via leptogenesis which is calculable in these models and found to be of the right order of magnitude.