Papers for Monday, Apr 27 2026

The female ratio in science field, including astronomy, is still quite low in Japan. We, the Astronomical Society of Japan, are making efforts to equalise the gender balance. In this paper, we summarise our statistics, member's thinking shown in our questionnaire, the history and accomplishments of the day-care system during annual meetings, and other activities.

The STIS team monitors the time dependent sensitivity (TDS) of each grating with one from a set of three secondary CALSPEC standard stars: GRW+70D5824, AGK+81D266, and BD+28D4211. Here, we use the three primary CALSPEC White Dwarf standard stars, dubbed the standard star "triad" (GD71, GD153, G191B2B), as an independent set of standards to verify the accuracy of STIS TDS corrections derived from the TDS monitoring stars, increasing the sample for each STIS L-mode from one up to three or four standard stars. We focus on triad star observations using the STIS L-mode gratings (e.g., G140L, G230L, etc.) with the same configuration as our standard TDS monitoring programs, and compare the triad observations to the TDS pipeline trends. Our analysis indicates the relative net count rates inferred from the triad standards agree with the TDS trends derived from the TDS monitoring stars with average residuals $<$ 2$\%$ across the full wavelength range of STIS, suggesting our current TDS L-mode trends are reliable and robust. We note that the dispersion in the residuals does vary with wavelength, with the NUV showing the lowest spread ($\pm$ 0.32$\%$ at 2400-2500 Angstroms) and the NIR the largest ($\pm$ 1.32$\%$ at 9500-9900 Angstroms); however, this scatter is also seen in our measurements of the TDS monitoring stars and is more indicative of other instrumental effects. Our findings rule out long term deviations, such as variability in our TDS monitoring stars, within measurement uncertainties.

X-rays are emitted from the corona above the orbiting matter of the accretion disk and travel either directly to us or illuminate the disk. This illumination of the inner disk is enhanced by gravitational light bending, which focuses the rays towards the black hole and therefore towards the inner radii of the disk. These rays that hit the inner radii are reflected back to us, and we observe them in the X-ray reflection spectrum. In this work, we create novel general relativistic ray tracing simulations to investigate the effects of altering the geometry of the accretion disks of black holes on the most dominant part of the reflection spectrum, the iron K$\alpha$ line. Work demonstrating the effect of disk geometry on the iron line has been performed, though many previous analyses have assumed a simplistic system, consisting of a point-source corona with a flat and infinitesimally thin accretion disk. We extend these models to more realistic accretion disk approximations. These include a constant aspect ratio disk, a radiation-pressure-dominated Shakura-Sunyaev disk, an expanded inner disk that has a non-negligible scale height in its inner regions due to radiation pressure, as well as various warped disks. Using measurement uncertainties from XRISM, we find that non-negligible thickness in accretion disks underestimates the black hole spin, corona height, and inclination angle if fitted with a flat disk model. The warped disk model could not be fit with the flat disk approximation.

The Simons Observatory Large Aperture Telescope Receiver (SO LATR) is a next generation Cosmic Microwave Background camera equipped with > 60, 000 detectors operating at 100 mK. Maintaining these detectors at the correct temperatures and locations requires stiff, cryogenically insulating struts. In this paper we report the design and performance of a novel glue joint in a strut used in the SO LATR to achieve the required performance. We use a tapped hole and set screw to create a profile on the exterior and interior wall of the glue joint, respectively, which greatly increases the strength of that joint by changing the failure mode from adhesive to cohesive. The failure mode of the resulting glue joint is cohesive with a yield strength 10% higher than a comparable smooth-walled design, and an ultimate strength 33% higher. Comparisons of the measured yield strength to the predicted axial load on the strut from simulations results in a factor of safety for the strut of 7. These struts have been installed in the SO LATR for three years and have undergone numerous thermal cycles from 300 K to 100 mK with no evidence of damage to the glue joint.

In this paper, we introduce the first implementation of magnetic field-aligned hyperbolic diffusion for standard smoothed particle (magneto-)hydrodynamics (SPH), and its linear-exact gradient extension (LESPH). Hyperbolic diffusion differs from traditional parabolic methods by incorporating the physical characteristic speed of diffusing particles and is computationally faster. This work extends it to encompass field-aligned diffusion, linear-exact gradients, and linear reconstruction to limit dissipation. Several standard test problems are presented: a diffusing slab, diffusion around a ring, a Gaussian pulse, and the magneto-thermal instability (MTI). The MTI only grows for for LESPH with reconstruction, and not for SPH. Both LESPH and SPH remain stable while fully aligning diffusion to magnetic fields. LESPH is more accurate and converges faster in the L1 error norm. SPH and LESPH both see improvements when using when also using linear reconstruction. These methods apply to other diffusive transport such as cosmic rays, viscosity, or magnetic resistivity.

We present a Python spectrophotometry extraction tool tailored for fast-moving point sources detected in the SPHEREx mission, and use it to construct a set of 0.75-5.0 {\mu}m low-resolution ({\lambda}/{\Delta}{\lambda} ~ 50) spectrophotometry data products based on the SPHEREx Quick Release 2 (QR2) for a set of 6003 L0-Y1 ultracool dwarfs: 2050 known ultracool dwarfs, 3008 known photometric ultracool dwarf candidates, and 947 newly identified ultracool dwarfs. This work more than doubles the number of ultracool dwarfs with spectroscopy, from 3449 to 7402. We provide SPHEREx templates for each spectral subtype and a set of tools to assign automated spectral types. The QR2 data release generates spectrophotometry with an average signal-to-noise per spectral channel above ~10 for most objects with WISE W2 magnitudes of 14.0 mag and brighter. The compiled data set is made available publicly at this https URL, where new spectral compilations from future data releases will also be made available as they are published. These new data provide a significant increase in the number of substellar objects for which the 2.4-5.0 {\mu}m window is now accessible, making it possible to probe important molecular chemistry of key CNOS-bearing species for the coolest brown dwarfs. We flag 2668 ultracool dwarfs as candidate young brown dwarfs, 250 as candidate subdwarfs, and 865 as possibly otherwise peculiar for future investigation. The SPIFF library presented here opens the doors to efficient confirmation of candidate substellar objects and follow-up studies of population-level atmospheric properties of cold brown dwarfs.

We improve upon the parametric model for the evolution of the density profiles of self-interacting dark matter (SIDM) halos introduced in Yang et al. (2024b), by considering the effects of mass accretion on a SIDM halo's gravothermal evolution. The original parametric model accurately predicts parameters $V_{\max}$ and $R_{\max}$, but with a tendency to overpredict $V_{\max}$ at $z=0$ for a subset of field halos. This discrepancy results from the parametric model predicting a faster rate of gravothermal evolution for these field halos compared to that measured in cosmological zoom-in simulations. We propose that the effects of mass accretion on the evolution of SIDM halos are not fully captured by the original parametric model. Our extended parametric model assumes that smooth mass accretion delays core-collapse by driving the SIDM halo back toward a Navarro-Frenk-White (NFW) profile (as it would have in the case of cold dark matter). We find that this extended model is able to substantially reduce the error in predicted $V_{\max}$ for halos compared to the original model, providing a more accurate model of SIDM halo evolution.

The Habitable Worlds Observatory (HWO) aims to image and characterize at least 25 ExoEarth candidates (EECs). Achieving this goal requires a detailed understanding of the observatory's design trade space, including the operational efficiency of the EEC survey. This study quantifies the impact of two critical parameters: the instantaneous field of regard (FoR) and the number of characterization observations required per EEC ($N_\text{char}$). We introduce a novel dynamic scheduling algorithm implemented within the EXOSIMS framework that models information gain during the mission. The scheduler models the orbital information known about each planet and forecasts detection probabilities to make scheduling decisions. We explore a multi-dimensional trade space, varying aperture size (6.5 m and 8.0 m), dedicated EEC survey time (2.5, 5.0, 7.5 years), $N_\text{char}$ (1 to 4), and FoR ($15^\circ$ to $135^\circ$). Our results demonstrate that the FoR is a major driver of the mission yield, with the yield decreasing significantly when the FoR is less than $90^\circ$. We find that increasing $N_\text{char}$ imposes a significant cost to mission yield, as each additional characterization required reduces yield by approximately 22%. The cumulative impact of requiring four characterizations instead of one lowers the yield by approximately 52%. This harsh penalty can be partially mitigated by increasing the survey duration. The relative yield loss when increasing $N_\text{char}$ from 1 to 2 is 38% for a 2.5 year survey and 14% for a 7.5 year survey. Our results highlight the complex interactions between HWO's engineering constraints and science requirements, and emphasize that the EEC survey efficiency is a critical component of HWO's design space.

Forthcoming Stage-IV dark energy optical surveys, such as LSST, have the ambitious goal of measuring cosmological parameters at sub-percent precision. Realizing their full scientific potential requires very precise measurement of the cosmic shear signal and control of corresponding systematics. In this work, we present a modern implementation of the Bayesian shear inference framework in Schneider et al. (2014), in the case that the PSF and sky background are known. This framework automatically propagates the pixel-noise measurement error from each galaxy into the final shear estimate, and thus requires no external calibration to handle noise bias. As a first application of this new implementation, we infer the cosmic shear posterior from simulated images consisting of isolated exponential galaxies with LSST-like levels of shape and pixel noise. In this simplified scenario, we estimate the absolute multiplicative bias $|m|$ of our approach to be below $0.9 \times 10^{-3} \, [3\sigma]$ when the intrinsic distribution of galaxy properties is known, and below $1.3 \times 10^{-3}\, [3\sigma]$ when these distributions are inferred alongside shear. Both results are within the LSST requirement of $|m| < 2 \times10^{-3}$. Additionally, we make progress towards the algorithm's computational feasibility in the context of modern wide-field surveys, where billions of galaxies must be processed, by leveraging differentiable forward models of galaxies, gradient-based samplers, and GPUs. Our final galaxy-fitting MCMC produces $300$ effective samples of galaxy properties in $0.45$ seconds per galaxy using a single A100 GPU. In the future, we seek to generalize our algorithm to handle selection, detection, and model shear biases so it can be applied to real survey data.

If a Large Language Model (LLM) can replicate your scientific contribution, the problem is not the LLM. What does it say about our field that so much of the anxiety about AI comes down to the fear that a machine could do what we do? Perhaps it says we should be doing something better.

Silicon-platelet feedhorn arrays are an established technology at millimeter wavelengths that, for some applications, can provide significant advantages over traditional direct-machined metal feedhorns. The Prime-Cam focal planes operating in the 350 GHz ($\sim$860 $\mathrm{\mu}$m) and 850 GHz ($\sim$350 $\mathrm{\mu}$m) bands are anticipated to carry the first silicon-platelet feedhorn arrays to operate fully at submillimeter wavelengths, representing a significant step forward in the application of this technology. In particular, the feedhorns designed for operation in the 850 GHz band represent a 3x increase in frequency compared to previously demonstrated and deployed devices of this type. Here we present a demonstration of silicon-platelet feedhorns at these submillimeter wavelengths, including in-lab performance characterization. We present fabrication metrology, room-temperature beammaps, and cryogenic optical efficiency measurements where the feedhorns are coupled to prototype CCAT Prime-Cam detectors. We show that feedhorn performance measurements are well matched to simulation and compare that performance directly to traditional, direct-machined metal feedhorns.

Ultra-high-energy (UHE) neutrinos are unique cosmic messengers that can traverse cosmological distances unattenuated, providing direct insight into the most energetic processes in the universe. Radio detection offers significant advantages for detecting highly inclined air showers induced by UHE neutrinos. This is due to a larger exposure range compared to particle detectors, which is a result of minimal atmospheric attenuation of radio signals combined with good reconstruction precision. Furthermore, this technique improves the air shower longitudinal reconstruction, which can be used to identify neutrinos with their first interaction far below the top of the atmosphere. In this work, we present a method for identifying UHE neutrinos using ground-based radio antennas. A reconstruction algorithm is introduced based on the radio emission maximum ($X^{\text{radio}}_{\text{max}}$), which demonstrates its power in distinguishing deeply developing neutrino-induced showers from background cosmic rays. Using simulations of $\nu_e$-CC-induced air showers, we evaluate the trigger efficiency, reconstruction performance, and resulting effective area and aperture prediction for a reference array. Our results show that radio detection significantly enhances the sensitivity to very inclined showers above 1 EeV, complementing traditional surface detectors. This technique is highly scalable and applicable to future radio observatories, such as GRAND. The proposed reconstruction and identification strategy provides a pathway toward achieving the sensitivity required to detect UHE neutrinos.

Fast Radio Bursts (FRBs) probe baryons permeating the cosmic web through their dispersion measures (DMs), which encode the integrated electron density along cosmological sightlines. Using 3,455 unique FRB sources from CHIME/FRB with $\sim 15$ arcmin localizations, we present an anthology of DM correlations with tracers of large-scale structure and baryonic matter at redshifts $z \lesssim 1.5$. We measure statistically significant correlations at $2.6-5\sigma$ with ten probes, including galaxies ($2.8\sigma$), weak gravitational lensing ($2.6\sigma$), cosmic infrared background ($4.0\sigma$), cosmic microwave background (CMB) lensing ($3.3\sigma$), thermal Sunyaev Zel'dovich (tSZ) effect ($3.8\sigma$), X-ray emission tracing galaxy clusters ($5.0\sigma$) and superclusters ($3.3\sigma$), soft X-ray background (SXRB, $4.1\sigma$), and radio continuum emission ($3.2\sigma$). These measurements reveal a consistent picture in which FRB sightlines intersecting overdense environments carry systematically larger DMs. Correlations with hot-gas tracers provide additional leverage on the strength of feedback, as they are strongly weighted towards the dense, bound gas. The measured amplitude of tSZ$\times$DM and SXRB$\times$DM correlations are consistent with theoretical predictions of baryon distribution from a DM-$z$ relation-inferred model with moderate feedback at $\sim 0.5\sigma$ level. Weaker feedback scenario is ruled out at $\sim 3.5\sigma$ by the SXRB$\times$DM correlation. Taken together, these measurements constitute a quantitative multi-tracer foundation for a new era in which FRBs from next generation facilities, such as BURSTT, CHORD, DSA, and SKA, in harmony with other probes, will map the baryon content of the full extent of the cosmic web.

Pulsar search is one of the main tasks for the Square Kilometre Array (SKA), implemented in the central signal processor (CSP) sub-element. As most the characteristics of undiscovered pulsars are unknown by definition, exhaustive searches over a multi-dimensional parameter space are employed. One main compute-intensive task of the pulsar search modules in the CPS is the matched filter group, which convolves the input signals with a group of large FIR filters. High-performance designs on FPGAs have been proposed that can process multiple large filters efficiently. But given that in many applications, including the here targeted pulsar search, FIR filters have many different sizes, there is further potential for optimisation. This paper investigates the optimisation of matched filtering designs. While the results are tranferable to other domains, we are motivated by the needs of the SKA pulsar search engine. The influence of changing number of filters and the difference in sizes is analysed. The generic design in time-domain (TD) is optimised by employing the longest processing time (LPT) first rule to distribute filter templates across filter processing pipelines. For the Fourier-domain (FD), the relationship between the required off-chip memory space and speedup over the generic design is investigated. To put the results into relation with with GPU design, we compared with a well-optimised design for top-end GPUs (NVIDIA Tesal P100). While a mid-range Intel Arria 10 is up to 7.5x slower than the P100, the performance per watt is slightly better on the Arria 10.

We study short-timescale 340 GHz flux-density variability of Sgr A* using ALMA Cycle 3 observations. Careful self-calibration enabled 10 s snapshot imaging with very high effective image-domain SNR, allowing high-cadence monitoring of Galactic Center sources. To reduce atmospheric and instrumental effects, we measured Sgr A* relative to multiple non-variable sources in the same field and corrected apparent variability caused by time-dependent u-v coverage and PSF changes using simulations with a static input model. We then searched for characteristic timescales over 20 s < tau < Tobs/3 using structure functions, the Lomb--Scargle method, and state-space-model autoregressive spectral analysis. No dominant narrow periodicity is found. Instead, the data show a short-timescale flat, white-noise-like regime at tau below about 2.3--6.3 min, followed by red-noise-like behavior at longer timescales. This flat regime appears in both active and quiescent phases, suggesting statistically independent fluctuations on these timescales. We interpret its upper boundary as an empirical transition timescale between decorrelated short-timescale fluctuations and longer-timescale correlated variability. The physical origin of this flat component remains uncertain, since previous theoretical and numerical studies more commonly report red-noise-like or broken-power-law variability.

Light bridge (LB) jets offer a unique window into small-scale eruptive phenomena within sunspots, the Sun's strongest magnetic environments; however, their generation mechanism remains a subject of debate. Using high-resolution observations from the New Vacuum Solar Telescope (NVST), we investigated six recurrent light bridge jets and the slipping motions of their jet base points (JBPs). Analogous to coronal jets, our observations show that these LB jets are characterized by a preceding JBP followed by a collimated jet spire. The JBP of each repeated jet along the LB displays apparent slipping motion at velocities of 0.6-1.5 km/s, which is temporally correlated with quasi-periodic enhanced photospheric horizontal motion of 1.3-6.5 km/s. Following the slipping JBPs, the resulting jet spires' fronts display similar slipping behaviors within the upper solar atmosphere. The Chinese Ha Solar Explorer (CHASE) reveals Ellerman-bomb-like spectral signatures at the JBPs, confirming that magnetic reconnection is operating at the jet base. Based on these results, we propose that repeated 3D reconnection occurring between the horizontal LB field and the ambient vertical umbral field may drive these LB jets. This process appears to be driven and/or modulated by quasi-periodic horizontal motion fueled by convective upflows and the transport of magnetic flux along the light bridge. This work suggests that some LB jets share a common reconnection-driven mechanism with coronal jets and provides direct evidence of slipping reconnection occurring along the sunspot light bridge.

To date, three samples from near-Earth asteroids have been delivered to Earth by Japan's Hayabusa (2010) and Hayabusa2 (2020) missions, and the United States OSIRIS-REx mission (2023). Free from terrestrial contamination, these pristine materials provide new opportunities to investigate planetary formation processes, the delivery of organics and water to the early Earth, and the nature of potentially hazardous asteroids. As analysis of the asteroid samples proceeds in laboratories around the world, we visit each of the missions, review the initial scientific findings, and explore the value of sample return in understanding our origins and protecting our future.

Galaxy evolution is commonly described through the time evolution of observational statistics such as luminosity functions and stellar mass functions. However, these quantities are projections of an underlying multivariate galaxy state space rather than fundamental dynamical variables. We develop a unified framework in which galaxy evolution is formulated as the time evolution of a probability measure on the galaxy manifold. Representing galaxy states by latent variables $\theta\in\mathcal{M}$ and the population by a density $\rho(\theta,t)$, the evolution is governed by a general equation containing continuous transport and nonlocal jump processes. By reinterpreting manifold learning as the pushforward of measures, we distinguish observational, representation, and physical measures, and emphasize that manifold coordinates themselves need not carry direct physical meaning. In this picture, luminosity functions and stellar mass functions arise as projected observables of a single underlying dynamics, and generally do not form closed equations in observational space. The framework contains existing models as limiting cases: reduction to a single mass variable yields continuity-equation models, while additive post-merger states recover the Smoluchowski coagulation equation. We further show that luminosity-function evolution is naturally described within the Schechter family, whose apparent stability is interpreted as an effective consequence of projection. Since observables are projections of measures, inference of galaxy evolution becomes a statistical inverse problem of recovering manifold dynamics from data. This framework shifts the focus from fitting observed statistics directly to inferring the underlying state-space dynamics, thereby bridging manifold learning and physical theory.

We apply caustic skeleton theory to the Manticore-Local simulations, which are Bayesian constrained reconstructions of the Local Universe from the 2M++ galaxy catalogue, and extract the three-dimensional multi-scale caustic skeleton of two canonical weblike structures in our Local Universe, namely the Coma Cluster and the Pisces-Perseus ridge as they represent the most prominent cluster node and filamentary artery in the nearby Universe. We show that the Caustic Skeleton network of caustic singularities accurately reproduces the observed large-scale organisation of galaxies in redshift space for one of the Manticore realisations. The hierarchy of caustic features allows us to establish a multi-scale classification of the large-scale environment in which observed 2M++ galaxies reside. One of the most interesting aspects of the theory is that it predicts two topologically distinct classes of filaments (A_4 swallowtail and D_4 umbilic caustics) that form through fundamentally different folding histories yet appear morphologically similar enough, on the surface, to be overlooked by conventional structure identifiers. We find that the influence of D_4 filaments only becomes increasingly relevant towards smaller scales, and the Pisces-Perseus Supercluster in particular is revealed to be a distinctly D_4-dominated structure compared to the extended Stickman structure around the Coma Cluster. In other words, caustic skeleton theory enables a novel topological characterisation of one of the most studied filamentary complexes in the nearby Universe. [Shortened]

Studies of magnetospheric accretion and magnetic field topology in T Tauri stars have advanced over the years, but their applications to fully convective, very-low-mass T Tauri stars remain relatively unexplored. We aim to analyze the circumstellar environment of the very-low-mass dipper-like star JH 223 by investigating the accretion process and characterizing its large-scale magnetic field topology. We analyzed the photometric variability of JH 223 using observations from multiple telescopes, including K2, TESS, and LCOGT. Additionally, we used Gemini/GRACES spectroscopic and CFHT/SPIRou spectropolarimetric data to investigate the star-disk interaction and characterize the large-scale stellar magnetic field using Zeeman-Doppler imaging. JH 223 is a fully convective classical T Tauri star with an age of about 3 Myr and a mass of 0.4 M$_{\odot}$. The large-scale surface magnetic field is predominantly poloidal, with a 250 G dipolar component. The dipole field strength and mass accretion rate indicate that the disk truncation radius is near the corotation radius. The star-disk interaction, combined with the inclined dipole, generates accretion columns that warp the inner disk. As the star rotates, this warp periodically obscures the stellar surface every 3.31 days, producing dipper light curves. The same period is also detected in radial velocity and longitudinal magnetic field variability. The accretion columns, traced by redshifted absorption in H$\alpha$ and He I 1083 nm, are associated with the inner disk warp at the same rotational phase. The accretion process in JH 223 is dynamic, transitioning from an unstable to a stable regime over a few weeks, consistent with magnetohydrodynamic simulations of star-disk interaction. Results from multi-technique observations suggest that the magnetospheric accretion model remains valid for fully convective very-low-mass young stars.

Asteroseismology enabled measuring the rotation rate in the deep stellar interiors of stars across several evolutionary phases, advancing the theory of angular momentum transport in single stars from the main sequence to the white dwarf phase. However, binary stellar evolution products have not yet been studied in the context of angular momentum transport constrained by asteroseismology. Hot subdwarf B (sdB) stars can pulsate in non-radial modes, enabling probing of their internal rotation. Those in binary systems form through mass transfer, thus they can be used to probe theories of internal rotation in post-mass transfer stars. Here, we interpret observed asteroseismic core and envelope rotation rates of sdB stars in unsynchronised binary systems that formed through the common-envelope channel, using stellar evolution models of rotating sdB stars with internal magnetic fields. We find that when sdB stars form with the angular momentum content of red giant cores prior to common-envelope ejection, their predicted core rotation rates are two to ten times lower than measured asteroseismic rotation rates, and their envelope rotation rates are lower by two to five orders of magnitude. This suggests that the angular momentum content of sdB stars increases during their formation. Since sdB stars in close binary systems may host circumstellar matter from a past common-envelope ejection, we show that if they accrete a small amount of matter, the combination of internal magnetic fields with angular momentum transfer through accretion spins up both the core and envelope to match their measured asteroseismic rotation rates.

Europium plays a key role in studies of nucleosynthesis in the rapid (r-) process of neutron capture nuclear reactions and the evolution of the r-process element abundances in galaxies. We refine the method for analyzing the Eu~II lines in stellar spectra by updating the Eu~II model atom with recent data on the rate coefficients for inelastic processes in the Eu II + H I collisions. The method was tested by deriving abundances from the Eu~II 4129 A and 6645 A lines in the Sun and five reference stars with well determined atmospheric parameters and high-resolution spectra available. It was shown that abandoning the local thermodynamic equilibrium (LTE) assumption allows the abundances from the two lines to be matched within the error bars, while the difference in the LTE abundances amounts to -0.09 dex for the Sun and -0.07 dex to -0.18 dex for the metal-poor stars. Accounting for the departures from LTE (non-LTE effects) for Eu~II leads to reducing the statistical errors of the derived stellar Eu abundances. The non-LTE abundance corrections for 11 lines of Eu~II were calculated for a grid of the plane-parallel MARCS model atmospheres with effective temperatures, surface gravities, and metallicities relevant to F-G-K type stars. They are publicly available and can be applied for improving stellar Eu abundances in studies of the galactic chemical evolution.

We present a detailed X-ray study of the quiescent and flaring coronae of three active main-sequence K-type stars, V834 Tau, LQ Hya, and BY Dra, using \textit{XMM-Newton} observations. The quiescent coronae are well described by two-temperature thermal plasma models, with cool and hot components at 0.26-0.30 keV and 0.93-1.01 keV, respectively. Despite similar coronal temperatures, X-ray luminosities (10$^{29.18\mbox{--}29.75}$ erg s$^{-1}$) and overall abundances, the relative emission measures of the cool and hot components differ among the stars. High-resolution spectroscopy reveals significant iron depletion by factors of 5-10 relative to photospheric values, and an inverse first-ionisation-potential effect in all three stellar coronae. Six energetic flares are detected, with peak temperatures of 30 -- 133 MK and released energies of $0.6\mbox{--}4.2\times$10$^ {33}$ erg, classifying them as superflares. Most flares exhibit decay times roughly twice their rise times, although one event shows a decay phase nearly twenty times longer than its rise. Time-resolved spectroscopy and loop scaling laws yield flare parameters consistent with previous studies of active stars. LQ Hya displays recurrent superflares at the same rotational phase across observations separated by six months, suggesting a long-lived, complex magnetic field structure. These results provide insights into the magnetic activity and flare energetics in active stars, and their implications for stellar and exoplanetary environments.

Broad, skewed iron K$\alpha$ emission lines in the X-ray spectra of accreting black holes encode key information about the spacetime geometry of the innermost disk. While the Kerr metric is standard for spin measurements, horizonless alternatives like traversable "Kerr-like" wormholes can mimic many black hole signatures, challenging current data interpretations. We develop a relativistic reflection framework incorporating Kerr-like wormhole geometries to predict iron line distortions and assess the feasibility of distinguishing event horizons from wormhole this http URL a custom ray-tracing subroutine, we implement two \textsc{XSPEC} modules: \texttt{kwline} for $\delta$-function profiles and \texttt{kwconv} for full reflection spectra, parameterized by spin, throat radius, and shape-function coefficients. We compute a dense grid of line profiles and generate synthetic \textit{NuSTAR} spectra with realistic response matrices. By fitting these simulations with canonical Kerr models, we quantify deviations attributable to wormhole this http URL find that Kerr-like wormholes produce narrower Fe K$\alpha$ lines with suppressed red wings as the throat parameter $\lambda$ increases. In 50 ks \textit{NuSTAR} simulations ($\lambda=0.9, a_*=0.998$), simple convolutional models (\texttt{kerrconv}) can mimic the wormhole spectrum. However, self-consistent models like \texttt{relxillCp} result in statistical failure, yielding structured residuals and unphysical parameter pegging (e.g., emissivity $q_{\rm in} \to 10$). We conclude that large-throat wormholes are detectable in high-quality X-ray spectra if analyzed with fully consistent reflection models rather than post-processing approximations.

International student mobility plays a critical role in shaping future research careers, particularly in highly globalized fields such as astrophysics. The Leiden/ESA Astrophysics Program for Summer Students (LEAPS) offers a 10-week, fully funded research program at Leiden Observatory and the European Space Agency's ESTEC centre for undergraduate and master's students. Designed to foster early research involvement, LEAPS supports students from diverse academic and cultural backgrounds. Since its inception in 2013, LEAPS has hosted 194 students from over 40 countries. Data collected for 165 participants reveal that over 50% have progressed to Ph.D. studies, with some members of earlier cohorts already securing competitive international fellowships in astronomy. LEAPS participants have collectively contributed to at least 25 peer-reviewed publications and 13 international conference presentations. LEAPS has contributed successfully in preparing undergraduates for research careers in astrophysics through hands-on experience, mentorship, and scientific exposure. By addressing barriers related to financial means and promoting diversity, the program not only enhances individual career trajectories but also contributes to the broader goal of inclusive academic mobility. Continued efforts are needed to further increase global representation and assess long-term impacts on participants' scientific careers.

We conducted an extensive identification and analysis of various morphological classes and subclasses of radio galaxies using the latest high-resolution data from the second data release of the LOFAR Two-Metre Sky Survey (LoTSS DR2). This paper presents the first results of our large-scale investigation: a new catalog of ``winged" radio galaxies (WRGs). These objects represent a fascinating class of irregular radio galaxies, characterized by a pair of secondary radio lobes (``wings") in addition to the primary active lobes. We identified and cataloged 621 new WRGs and 403 additional candidates. Among the confirmed winged sources, 382 are classified as ``X"-shaped radio galaxies (XRGs), while the remaining 239 are ``Z"-shaped radio galaxies (ZRGs). We also estimated several basic parameters for these winged sources and performed a Fanaroff-Riley (FR) classification. Our results show that the majority of the sources ($\sim$88\%) exhibit edge-brightened radio lobes and high average radio power ($\rm log_{10}[P_{144MHz} / W Hz^{-1}]$ = 26.25), consistent with an FR-II classification. The average spectral index between 144 MHz and 1.4 GHz is --0.84, which is steeper than that found for previously identified winged sources based on higher-frequency data from the VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. This indicates that our study is capable of detecting fainter sources. The median linear size of the winged sources, 498 kpc, confirms that these are large-scale structures, with approximately 16\% having sizes exceeding 0.7 Mpc, making them potential candidates for giant radio galaxies.

Mid-infrared astronomy from the ground faces critical challenges in accurately detecting and quantifying sources due to the dominant spatially and time-variable background noise. Moreover, chopping and nodding, the traditional methods for dealing with these background issues, will not be technically feasible on the next generation of extremely large telescopes. This limitation requires the development of novel computational methods for a robust background reduction. We present and evaluate a novel method named LOw-RAnk Background ELimination (LORABEL) to improve the sensitivity of mid-infrared astronomical observations, without the need for classical telescope nodding, source masking, or other overheads in observing time. We applied a low-rank background-reduction strategy to (1) data taken on the ground with the VISIR with synthetically injected sources, and (2) airborne data from SOFIA. We compared the performance of our new method to classical chopping and nodding techniques, and analysed the effect on source photometry and detection precision for different observational scenarios. In regimes with a low signal-to-noise ratio (S/N $<5$) in the ground-based VISIR data, LORABEL reduces variation in the photometric error with respect to chopping differences alone and even the classical chop-nod sequence, at the cost of introducing a bias. Secondly, we demonstrate that LORABEL increases detection precision in comparison to traditional background-reduction methods. For the SOFIA dataset, we achieve a $20-100$ fold decrease in mean background flux with respect to the traditional chop-nod method while preserving most of the source flux. Our findings suggest that LORABEL is applicable to a wider range of instrumental observation, that is, both ground-based and airborne, and it is a suitable tool in the context of faint-source detection.

A sample of 139 young open star clusters closely associated with the Radcliffe wave is considered. Modeling their spatial distribution and kinematics over a time interval of 30 Myrs ago and 30 Myrs into the future revealed that they exhibit the main properties characteristic of a Radcliffe wave over the past 10-15 Myr. They are distributed on the galactic XY plane as a long and narrow chain inclined to the Y axis, and exhibit a wave-like behavior of their vertical coordinates up to 15 Myr in the past. This behavior of their vertical coordinates will persist over the interval of 15-20 Myr in the future. A new finding is the presence of vertical perturbations with an amplitude of deviation from the galactic symmetry plane of up to 200 pc over the entire time interval considered in the past, up to -30 Myr. This result calls into question the possibility of using a scenario in which the initial disturbance of the interstellar medium is assumed to be the Parker instability of the galactic magnetic field.

In modern cosmology, the rapid growth of high-precision observational data, along with significant theoretical advances, has intensified the challenge of identifying a robust, model-independent framework to probe the expansion history of the Universe. In this work, we propose a novel artificial neural network (ANN)-based framework for the non-parametric reconstruction of the late-time cosmic expansion. The framework is trained and validated through a three-stage screening pipeline prior to its application to real observational data. As a demonstration of its effectiveness, we reconstruct the Hubble parameter $H(z)$ using the latest cosmic chronometer measurements. Our results show that the reconstructed expansion history aligns with the predictions of the $\Lambda$CDM model within observational uncertainties, thereby supporting the robustness and reliability of the proposed approach.

Cold accretion and quenching are closely related aspects of galaxy evolution, as sustained gas supply is required to maintain star formation. High-redshift galaxy groups therefore provide a valuable laboratory for testing how the thermal state of accreting gas relates to the emergence of quiescence. We measure quiescent fractions in a sample of 16 spectroscopically confirmed galaxy groups at $1.6<z<3.6$, spanning halo masses from $10^{12.8},{\rm M_\odot}$ to $10^{13.9},{\rm M_\odot}$, by fitting the SEDs of candidate member galaxies selected from the COSMOS2020 catalog and using a membership-probability approach to estimate group quiescent fractions. We compare these quiescent fractions to the expected cold or hot accretion state of each halo and find evidence for a correlation: quiescent fractions reach about 50 percent in groups in the hot-accretion regime and are consistent with zero in groups in the cold-accretion regime. In mature hot-accreting groups, massive quiescent galaxies are preferentially found in the inner regions ($R<0.5R_{\rm vir}$), with a 4.4-sigma excess relative to the outskirts. Most groups lack a clearly established brightest group galaxy and instead show small stellar-mass gaps, typically $M_{*,1}/M_{*,2}<3$, indicating that they remain in an active assembly phase rather than being dynamically evolved systems. Consistently, the stellar-mass excess of the dominant galaxy, measured relative to the SHMR expectation, does not predict the group quiescent fraction. Taken together, our results support a picture in which the cold-to-hot transition in gas accretion contributes to the onset of quiescence, possibly through inside-out starvation associated with filament disruption in shock-heated intra-group gas, and suggest that environment plays a greater role than internal processes in shaping the quiescent galaxy population in these structures.

Gravitational-wave (GW) observations of compact binary coalescences (CBCs) are traditionally interpreted under the assumption that the binary evolves in isolation. However, in realistic astrophysical environments, brief three-body encounters may perturb the binary's orbital evolution and imprint deviations on the emitted GWs. We develop a physically motivated model for such interactions, retaining Newtonian three-body dynamics supplemented by leading-order ($2.5$PN) radiation-reaction within the binary. We show that such encounters produce a distinctive morphology of dephasing and amplitude modulation in GWs. We search for this kind of distortion from the LIGO--Virgo--KAGRA (LVK) GW catalog GWTC-4 on three events: GW170817, GW190814, and GW230627\_015337, chosen based on high SNR and in-band duration $\gtrsim 10~\mathrm{s}$. We find no statistically significant deviation in the data, which translates into constraints on the absence of any intermediate-mass black hole in the mass range above $\sim 10^2$ M$_\odot$ in the vicinity of these binaries of radius approximately $10^{-1}~\mathrm{AU}$. This arises from robust exclusions arising from fly-by interactions that would dynamically disrupt the binary and are directly ruled out independent of waveform modelling, placing the first upper bound on intermediate-mass black holes near these GW events. In future, with the availability of long-duration GW signals, this new avenue can probe encounters of the binary GW sources with compact objects of lighter masses at distances farther away than 1 AU and hence opens a new window to probe the population of individual compact objects of both astrophysical and primordial origin in astrophysical systems of dense environments ranging from galactic centers to dense globular clusters.

Studying the interstellar medium (ISM) in merging high-redshift galaxies is crucial for understanding early galaxy assembly, star formation, and black hole growth, predicted by hierarchical $\Lambda$CDM models. Deep imaging and spatially resolved spectroscopy with JWST enable unprecedented insight into these processes, even for galaxies in the Epoch of Reionization. We present NIRSpec and MIRI integral field spectroscopy and MIRI imaging of the merging galaxy Gz9p3 at z=9.3 of the UV and optical rest-frame showing a clumpy morphology in the continuum as well as line emission covering the entire galaxy over a range of 5 kpc from the central clump to the tail region. We analyze the integrated spectrum as well as different apertures in the galaxy allowing a spatially resolved characterization of the ionized ISM of this galaxy. We compare our measurements with archival NIRCam imaging and ALMA data. We measure a total star formation rate of 13.4 $\pm$ 1.8 Msun yr$^{-1}$, a metallicity of 12+log(O/H) = 7.84 $\pm$ 0.05 and $\xi_{ion}$= 25.4 $\pm$ 0.1 erg$^{-1}$ Hz and a burstiness parameter of 0.9 $\pm$ 0.1 for the integrated spectrum. We find large spatial differences in these parameters between the central clump and the tail region. The optical [OIII] emission peaks in the main galaxy, the far-infrared [OIII] emission peaks towards the tail, indicating different physical conditions in the ISM of the tail and main galaxy. This study presents the spatially resolved ISM analyses of a galaxy at z>9, revealing nebular line emission and strong spatial variations in star formation, metallicity, physical conditions, and ionizing efficiency. The results indicate a recent, metal-poor starburst in a tail alongside a more evolved, enriched central clump with evidence for extreme excitation. This demonstrates the power of spatially resolved JWST spectroscopy of galaxies in the Epoch of Reionization.

Parametric decay instability (PDI) of Alfvén wave is thought to play an important role in the dissipation of the large-amplitude Alfvén waves and in the heating of magnetized plasmas. Temperature anisotropy is frequently observed by spacecraft, including Parker Solar Probe (PSP), in the near-Sun solar wind, yet its impact on PDI in the near-Sun solar wind has been understudied. We calculate the maximum growth rates of PDI, $\gamma_{\max}/\omega_{0}$, where $\omega_0$ is the frequency of the parent wave, by solving the linear dispersion relation of Chew-Goldberger-Low (CGL) equations under several expanding background models. To assess the effect of temperature anisotropy, the growth rate is compared with that derived from ideal magnetohydrodynamics (MHD). From $R_0$ ($ \sim 1.02R_\odot$) to $30R_0$, we consider three expansion cases: (i) spherically symmetric adiabatic expansion with constant wind speed, (ii) Multi-source observation- and model-constrained expansion, and (iii) a PSP-constrained profile of $(\beta_{\parallel},\xi)$, where $\beta_\parallel=8\pi p_{\parallel0}/B_0^2$ is the parallel plasma beta and $\xi=T_{\perp0} / T_{\parallel0}$ is the temperature anisotropy, that includes Parker-spiral effects. We find that temperature anisotropy increases $\gamma_{\max}/\omega_{0}$ for $\beta\lesssim 0.1$ in the near-Sun solar wind: in the case of (iii), temperature anisotropy with $T_{\perp0} > T_{\parallel0}$ increases $\gamma_{\max}/\omega_{0}$ by factors of $\sim 1.5$ over $R\simeq 1$--$15\,R_0$, whereas temperature anisotropy with $T_{\parallel0}>T_{\perp0}$ decreases $\gamma_{\max}/\omega_{0}$ at larger $R$. Our results suggest that the temperature anisotropy plays an important role in the onset of PDI even in low-$\beta$ regimes, such as the near-Sun solar wind.

PSO J083.8371+11.8482, a quasar at $z = 6.34$ with a nearby companion galaxy, provides an opportunity to study the impact of active galactic nucleus (AGN) activity on the surrounding environment during the epoch of reionization. We analyze ALMA observations of the [C\,\textsc{ii}] 158~$\mu$m emission line and the far-infrared (FIR) continuum, which trace cold interstellar gas and dust-reprocessed radiation from star formation and AGN heating. The quasar host shows star formation rates (SFRs) of $544$--$3764~\mathrm{M_{\odot}~yr^{-1}}$ from [C\,\textsc{ii}] and $1861$--$2932~\mathrm{M_{\odot}~yr^{-1}}$ from FIR emission, while the companion galaxy exhibits lower SFRs of $21$--$145$ and $76$--$211~\mathrm{M_{\odot}~yr^{-1}}$ from the same diagnostics. Both galaxies follow typical $L_{\mathrm{[C\,II]}}/L_{\mathrm{FIR}}$ ratios observed in star-forming galaxies and show no evidence for a [C\,\textsc{ii}] deficit, indicating that stellar heating dominates the interstellar medium energetics. The [C\,\textsc{ii}] moment maps reveal compact emission with centrally peaked intensity and ordered rotational kinematics in both systems. Velocity dispersions remain well below values associated with powerful AGN-driven outflows, and no significant morphological asymmetries or disturbed velocity fields indicative of AGN feedback or major mergers are detected, although marginal kinematic substructure in the quasar's high-velocity channels warrants further investigation. Although the companion lies at a projected distance of $18.248 \pm 0.277$~kpc within the quasar proximity zone, neither morphological nor kinematic signatures indicate AGN-driven outflows affecting the circumgalactic medium. We therefore interpret this system as being observed in a pre-outflow accretion phase, where rapid supermassive black hole growth precedes the development of large-scale AGN feedback.

Understanding the large-scale dynamics of molecular clouds (MCs) is crucial for constraining the processes that govern star formation and the structure and evolution of the Galaxy. While gas tracers have traditionally been used to map MC kinematics, stellar tracers such as young stellar objects (YSOs) and open clusters (OCs) provide a complementary approach that enables direct comparisons between the stellar and gaseous components. We aim to validate OCs as complementary tracers by testing whether they retain the same bulk kinematic imprint as YSOs, and to reconstruct the three-dimensional (3D) motions of the main MC complexes within 2.5 kpc of the Sun using YSOs and young OCs as tracers. Using Gaia DR3 astrometry together with complementary spectroscopic surveys for radial velocities, we compiled a unified sample of 24,732 stellar tracers. We applied robust clustering in proper motion space to identify co-moving YSOs and derived cloud-averaged motions via Monte Carlo sampling. These were compared with the kinematics of OCs younger than 30 Myr. Finally, we performed orbital integrations in a realistic Galactic potential to trace the past evolution of the clouds and quantify their expansion and rotation. We derive homogeneous 3D kinematics for 15 MC complexes within 2.5 kpc. YSOs and OCs exhibit strongly consistent kinematics, with a median spatial velocity offset of $\simeq 2$ km s$^{-1}$, confirming that both populations trace the bulk motion of their parent clouds. The resulting cloud kinematics show a median peculiar velocity of $\simeq 8.7$ km s$^{-1}$ with respect to Galactic rotation. We trace back the Solar System's voyage through the Orion cloud and the common origin of Lupus, Ophiuchus, and Corona Australis in Sco-Cen. Internally, we detect significant expansion in Orion and Ophiuchus ($5\sigma$) and coherent rotation in at least seven complexes.

In gamma-ray bursts (GRBs), the electron pitch angle ($\alpha$) is usually assumed to be isotropically distributed. However, recent numerical simulations indicate that only the high-energy electrons (with Lorentz factors $\gamma>\gamma_{iso}$) are distributed isotropically, whereas the low-energy electrons (with $\gamma<\gamma_{iso}$) follow an energy-dependent anisotropic distribution during magnetic reconnection. The mean value of $\sin^2 \alpha$ approximately follows the relation $\langle \sin^2 \alpha \rangle \propto \gamma^{m}$ for $\gamma<\gamma_{iso}$. In principle, polarization measurements may help us constrain the pitch-angle distribution of electrons in GRBs, since different pitch-angle distributions produce distinct synchrotron polarization signatures. The polarization of GRBs produced by isotropically distributed electrons has been extensively studied. In this paper, we investigate synchrotron polarization produced by anisotropically distributed electrons within a globally toroidal magnetic field in GRB prompt emission. Our results show that the synchrotron PDs in the $\gamma$-ray and X-ray bands produced by anisotropically distributed electrons are systematically lower than those produced by isotropically distributed electrons, while the PD in the optical band could be either lower or higher than that of isotropically distributed electrons, depending primarily on the value of the energy slope $m$. In addition, we compared our numerical results with observational data, and the comparison suggests that an anisotropic distribution of electrons may offer a potential explanation for the PD and spectral data of some GRBs.

Prior to core collapse, the neutrino emission from red supergiants (RSGs) is so large that a nearby ($\lesssim1$kpc) RSG will become visible in current and near-future neutrino detectors. The rate of emission and the spectra of the pre-supernova (pre-SN) neutrinos from RSGs are sensitive to the temperature, density, and detailed isotopic composition of the core. During the last year of the star's life, these properties change considerably. Several factors of stellar evolution modeling - such as the treatment of mass loss and convective overshooting - alter the thermal conditions and composition of the RSG core as it approaches collapse. In this paper we present the first study of how varying the treatment of mass loss and convective overshooting together affects the pre-collapse core properties and neutrino emission of RSGs. We use the stellar evolution instrument MESA and construct a grid of 32 models with zero-age main sequence masses of $\{ 12, 15, 18, 20\}$ $M_\odot$, use the so-called 'Dutch' mass-loss scheme with wind efficiencies of $\{0.2, 0.4, 0.8, 1.0\}$, and consider two convective overshooting schemes. Our models use a large 206-isotope nuclear network in order to accurately compute the structure and composition of the star. We find that, in the last few days of the star's life, the general trend of the conditions and composition in the core of the star is one of contraction, heating, and deleptonization, but that during this phase, this general trend will be interrupted by the initiation of core silicon burning and shell burning episodes that cause the core to expand and undergo convective mixing with material of a higher proton fraction that temporarily reverses the deleptonization. The pre-SN neutrino emission reflects these changes with a gradual shift to higher energies and larger flux that becomes dominated by beta processes a few hours prior to the collapse.

The radial acceleration relation (RAR) is a tight empirical correlation between the observed radial acceleration (a_tot) and the baryonic radial acceleration (a_bar) measured across galaxy radii: these two accelerations start to deviate significantly from each other below a characteristic acceleration scale, a0. So far, observational studies of the RAR have predominantly focused on galaxies in the local Universe, leaving its evolution with cosmic time largely unexplored. Using high signal-to-noise data from the MUSE Hubble Ultra Deep Field survey, we investigate the RAR with a sample of 79 star-forming galaxies (complete above M* >10^8.8 Msun) at intermediate redshifts (0.33 <z <1.44). We estimate the observed intrinsic acceleration and the baryonic acceleration from a disk-halo decomposition that incorporates stellar, gas, and dark matter components, with corrections for pressure support, using 3D forward modelling. We find a RAR in our intermediate-z sample offset from the local relation, with a higher characteristic acceleration scale, a0(z~1) = 2.38+/-0.1* 10^-10 m/s^2, and a larger intrinsic scatter (~0.17 dex). Dividing the sample into redshift bins and refitting the RAR in each bin, we find a characteristic acceleration scale that systematically increases with z. Parametrizing the z-dependence as a0(z)= a0(0) + a1 * z, we obtain a1 = 1.59+/-0.1 *10^-10 m/s^2, providing evidence for a z-evolution. We find similar results using various dark matter halo profiles as well as the Modified Newtonian Dynamics framework in our 3D forward modelling. Our results show that the RAR persists at intermediate redshift, with statistically significant redshift evolution of the characteristic acceleration, pointing to a possible evolution of the baryon-missing mass connection over cosmic time.

Primordial black holes (PBHs) provide a compelling interpretation for the binary black holes (BBHs) observed by ground-based gravitational-wave (GW) detectors, especially for those BBHs in the theoretical mass gap. In the early Universe, the scalar perturbations required to produce such PBHs inevitably generate scalar-induced GWs (SIGWs). These SIGWs peak in the sub-nanohertz band, and manifest secularly as measurable jerk-like drifts in the second derivative of pulsar spin periods. In this Letter, we perform the first search for SIGWs using pulsar parameter drifts, and place a 95\% confidence-level upper limit on the PBH abundance of $f_{\mathrm{PBH}} < 10^{-10}$ over the mass range $[3 \times 10^{-1}, 4 \times 10^{4}] M_{\odot}$. Our results strongly disfavor a PBH origin for the BBHs currently detected by the LIGO-Virgo-KAGRA (LVK) Collaborations.

The supermassive black hole Sagittarius A* (Sgr A*) exhibits temporal and spectral variability across the electromagnetic spectrum. However, variability at radio frequencies below ~ 5 GHz for timescales shorter than a day remains largely unexplored. We investigate the variability of Sgr A* at 2.79 GHz on short timescales (1 min), to probe an under-explored regime of its emission process. Through point-source model fitting in the uv-domain, we analyse the flux density variation of Sgr A* over an 8 h observation. We detect flux variation on a few tens of minute timescale with a modulation index of 6.11 %, a mean flux density of ($827 \pm 0.1_{\mathrm{stat}} \pm 33_{\mathrm{sys}}) \, \mathrm{mJy}$, and a mean spectral slope of $0.08\pm0.03$. Furthermore, we measure the slope of the structure function of the observed light curve as $0.81 \pm 0.05$ with a characteristic timescale of about 120 min. Our study at low radio frequencies is a critical step toward constraining the physical mechanisms that drive Sgr A*'s variable emission and its spectral energy distribution. Our study suggests that variability at centimetre and millimetre wavelengths is likely more closely related than previously thought.

Black-hole X-ray binaries (BHXRBs) in the hard and hard-intermediate spectral states commonly exhibit prominent type-C quasi-periodic oscillations (QPOs) in their X-ray power spectra. Despite extensive observational and theoretical efforts, the physical mechanism responsible for these oscillations has not yet been firmly established. The disk-corona system in BHXRBs is radiatively coupled, as hard X-ray emission from the corona can be reprocessed by the accretion disk and re-emitted as soft photons that contribute to cooling the coronal electrons. Aim of the present study is to examine whether this feedback can give rise to limit cycles having the spectro-temporal properties of QPOs. We model the coronal emission using a one-zone radiation framework and solve the time-dependent kinetic equations for electrons and photons. Electrons are energized by some unspecified process and cool via inverse Compton scattering of soft photons originating from (i) the accretion disk and (ii) disk reprocessing of the hard radiation produced in the corona. When electron cooling is dominated by soft photons reprocessed in the accretion disk, the disk-corona system undergoes limit-cycle oscillations. For a subset of the model parameters, these oscillations reproduce key properties of type-C QPOs observed in BHXRBs. The oscillation frequency depends on the coronal radius and on the energization timescale, while the resulting X-ray spectra are well described by power laws extending up to energies of ~ 10-100 keV. These calculations confirm and extend earlier semi-analytical results obtained with simplified treatments. Owing to the scale-invariant nature of the model, the results can be readily extrapolated to other accreting systems, such as Active Galactic Nuclei.

On the Hertzsprung-Russell diagram, F-type solar pulsators connect the Sun to intermediate mass stars located on the instability strip. With respect to lower mass stars, they are structurally peculiar in the sense that they are constituted of three distinct dynamical layers: a small convective core, a deep radiative interior, and a shallow convective envelope. Current asteroseismic techniques only provide limited information on the interior dynamics of these stars. Indeed only gravity modes (g modes), for which unambiguous characterisation is lacking, are able to probe the deep stellar layers. A better understanding of the excitation and behaviour in F-type solar pulsators is therefore necessary in order to consider their detection. In this work, we simulate for the first time the global stellar structure of an F-type star (core, radiative interior, envelope). We show that the contribution of the core strongly affects the spectrum of excited g modes, with low-order high-degree modes unable to form due to their interaction with the turbulent convection of the core. Finally, by computing the disc-integrated signature of the modes, we are able to demonstrate that they preserve their integrity up to the top of the convective envelope, which is a strong argument in favour of their detectability with spaceborne photometry.

Dark matter and baryons acquire a relative velocity after decoupling in the early Universe. Baryons are gravitationally unstable only above their Jeans scale, while cold dark matter (CDM) is unstable on all scales. We show for the first time that their relative drift triggers a resonant gravitational instability that drives sound waves in baryons. When the projected DM drift is subsonic, the stable oscillatory branch of baryons resonates with the Doppler-shifted DM mode, producing exponentially growing perturbations whose growth rates exceed the intrinsic CDM growth rate. The instability peaks below the baryon Jeans scale and, in baryon-dominated environments, opens a window of complete stability between the Jeans scale and the resonance. Supersonic drift suppresses growth, as previously noted. The resonant coupling also transfers momentum between the species, creating a non-viscous, collisionless drag. We derive an accurate analytical approximation for the growth rate at resonance and show that the associated timescales range from years to tens of millions of years across diverse environments -- planets, protoplanetary disks, stars, molecular clouds, galaxies, and galaxy clusters -- typically much shorter than their ages. In an expanding FLRW universe, the instability enhances baryon density perturbations at different redshifts for appropriately oriented modes while suppressing the growth of those aligned with the DM stream. The universe thus sings across all scales, and this resonant mechanism provides the means to listen: it offers a novel probe of dark matter through its seismic imprint on astrophysical objects and may explain long-standing puzzles such as the persistence of spiral arms and the heating of the intracluster medium in galaxy clusters.

We present 3D kilonova radiative transfer simulations for a series of binary neutron star merger models. The masses of the neutron stars are varied as well as the total mass of the system and two different equations of state were used (SFHO and DD2), producing a range in dynamical ejecta masses and elemental abundance patterns. In this paper, we focus on the bolometric light curves and spectra in the polar direction for comparison with observations of the kilonova AT2017gfo. We calculate line-by-line opacities and include new calibrated lanthanide atomic data. All of the simulated spectra show strong features from Sr II, La III, Gd III and Ce III, which appear to correspond to features identified in AT2017gfo, although the simulated features are generally more blueshifted. The models with the lowest lanthanide fraction in the polar direction also show a Y II feature. Ce III, Ce II, Nd III and Nd II play an important role in shaping the spectral continuum. While the bolometric luminosities in the polar direction vary with the ejecta mass of each model, we find only little sensitivity of the spectral properties to the merger configuration. Our study demonstrates that dynamical ejecta alone can reproduce (although at earlier times) many spectral properties of AT2017gfo, suggesting dynamical ejecta may have a strong impact on the early spectral evolution. However, future simulations are needed to also elucidate the role of other ejecta components for shaping the kilonova spectrum.

Compact stars serve as natural systems where matter exists at densities far beyond those achievable in laboratory experiments. Among them, magnetars are expected to possess interior magnetic fields that may reach values of the order of $10^{17}-10^{18}$ G. These extreme conditions are expected to alter the microscopic and macroscopic properties of dense matter. In this review, we examine how strong magnetic fields affect fermionic matter through mechanisms such as Landau quantization and anomalous magnetic moment interactions. We further discuss the behaviour of magnetized hadronic matter within relativistic mean-field approaches and consider the possible emergence of additional degrees of freedom, including hyperons, $\Delta$ resonances, meson condensates and quark matter. The consequences of these effects for neutron-star structure and observational constraints are also briefly outlined.

The May 2024 solar superstorm highlighted the vulnerability of rapidly expanding low Earth orbit (LEO) satellite networks to severe space weather events. To systematically evaluate LEO network resilience, we introduce an open-source tool, CosmicDancePro. It enables a comprehensive analysis of the effects of solar storms in the LEO satellite network. It integrates real-world multimodal datasets, including space weather measurements from several satellites, upper-atmospheric density conditions from data-driven and high-fidelity physics-based models, and LEO satellite trajectory and LEO network measurement traces to quantify orbital decay driven by enhanced atmospheric density and network connectivity degradation. We utilize CosmicDancePro to analyze the Starlink constellation's behavior during two recent major solar storms. First, we identify the specific fleet management strategies Starlink adopts during the May 2024 solar superstorm and how they differ from its regular orbit-correction strategy. Second, we identify the mechanisms driving the previously unexplained 'W'-shaped altitude variation pattern across orbital planes of LEO constellations. Finally, our network-layer analysis quantifies the connectivity degradation during these storms, revealing transient disruptions that include repetitive short-lived outages, reconfiguration latency spikes above 500 ms, up to 60% increase in uplink loss, distorted diurnal latency patterns, and a 10+ Mbps drop in end-user data rates during storm peaks.

This paper introduces the Wide-band Asp-Clean (\texttt{WAsp}) algorithm, a novel scale-sensitive image reconstruction method tailored for wide-band imaging applications. This algorithm is particularly beneficial for thermal noise-limited imaging with aperture synthesis telescopes, where joint spatio-frequency modeling of the sky brightness distribution is critical. The \texttt{WAsp} algorithm replaces the use of the MS-Clean algorithm in the MS-MFS algorithm with the {\tt Asp} algorithm \citep{Asp_Clean}, which itself has been improved for both imaging and runtime performance. With the high sensitivity of current and next-generation telescopes, spatio-frequency modeling in a scale-sensitive basis becomes crucial for ensuring that residuals align with the noise model across the frequency band. Although existing wide-band scale-sensitive algorithms have demonstrated superior performance over scale-insensitive counterparts, they often suffer from well-documented deficiencies, leading to significant wide-scale residuals in Stokes-I at low levels and consequently significant relative errors in spectral index maps. The \texttt{WAsp} algorithm addresses these limitations while maintaining computational efficiency. The implementation can be configured to support narrow-band and wide-band scale-sensitive imaging, spectral-cube imaging applications and joint single-dish and interferometer imaging. To demonstrate improved imaging performance, we show comparison with existing algorithms via carefully developed simulations for stress-testing the algorithms. We also present results from its application to real-world wide-band data, underscoring its effectiveness in practical imaging scenarios.

The Square Kilometre Array (SKA) is expected to start science operations in 2030 and by that time there could be up to 10$^5$ artificial satellites in Earth's orbit, comprising an increase of an order of magnitude compared to 2024. Most of these new satellites will belong to satellite megaconstellations aimed at providing communication services all over Earth. These satellites create radio frequency interference (RFI) that can impact the observations of modern radio telescopes. In this Letter, we forecast the amount of observing time for which the SKA interferometers will be exposed to satellites, risking RFI contamination. We employed an analytical model and considered two cases of exposure to satellites; (1) satellites that only lie in the main beam and (2) satellites that lie in the main beam or the first sidelobe. We show that for SKA-Low, the exposure is high, with satellites in the beam for 30% of the observation time across half of the frequency range, rising up to 100% below 100 MHz. For SKA-Mid, high frequencies are mostly spared, but observations below 1 GHz could also end up seeing satellites for at least 30% of the time. We conclude that satellites will be unavoidable during SKA observing conditions, risking a strong impact on the RFI environment. This will necessitate a concerted effort to obtain accurate measurements of satellite RFI and to improve our understanding of the impact on various science cases. Finally, new mitigation techniques that are less data-destructive than simple flagging must be introduced.

Highly magnified stars at cosmological distances ($z \gtrsim 0.7$) become detectable thanks to microlensing by intracluster stars near the critical curves of galaxy clusters. Multi-epoch photometric campaigns targeting caustic crossing galaxies magnified by massive galaxy clusters enable the detection of these objects as transient events. Such stars provide unique opportunities to study stellar populations at early cosmic times, probe the nature of dark matter, reveal small-scale structure in the cluster, and improve lens models. To date, only a few dozen high-redshift stars have been reported, with a single lensed galaxy, the Dragon, holding the current record of 44 detections. These numbers, however, remain insufficient to exploit their full potential. In this paper, owing to the inclusion of new observations, we report the identification of more than 100 magnified stellar events in the Dragon, behind the massive galaxy cluster Abell 370. The relatively low redshift of the Dragon ($z\approx0.725$) facilitates the detection of its most massive stars. Using imaging data from three different cycles (2022--2024) with the James Webb Space Telescope, we apply a time-domain technique to identify flux variations associated with caustic-crossing events. From the spatial distribution of stellar events we constrain the high-end slope of the stellar luminosity function, finding $\beta=2.18^{+0.20}_{-0.30}$. Alternatively, assuming a fixed slope, we constrain the microlens surface mass density. In addition, we examine the parity asymmetry of the detected caustic-crossing events, a proposed probe of wave dark matter, and find that it remains present. We also use the events to trace the regions of highest magnification, offering an alternative way to map the system critical curves.

Magnetic fields with field strengths between $10^{-17}\,$G and a few Nanogauss are expected to exist today in the intergalactic medium (IGM). Their origin is unknown, but may be of primordial nature, in which case they would have influenced the thermal and ionization history of the IGM as well as the growth of small-scale matter perturbations. In this work, we revisit constraints on Primordial Magnetic fields (PMFs) by consistently accounting for their energy losses through ambipolar diffusion and decaying turbulences from recombination through the epoch of reionization, which progressively reduces the magnetic field strength over time. We implement these effects in ${\tt HyRec}$ and ${\tt exo21cmFAST}$ to model the interplay between PMFs and astrophysical processes up to reionization. Using a neural-network emulator (${\tt NNERO}$), we perform a MCMC analysis that combines late-time probes of the reionization history and galaxy UV luminosity functions. We find that including PMF energy losses significantly relaxes previous bounds, as the reduced field strength suppresses their imprint on observables. Employing a Fisher matrix analysis, we estimate the sensitivity of the 21$\,$cm signal experiment HERA to the PMFs' imprint on intergalactic medium perturbations and show that 21$\,$cm cosmology could significantly improve on current bounds. Our results highlight the importance of modeling PMF evolution self-consistently with the IGM evolution to extract current bounds and future sensitivities.

The consequences of a protoplanetary disk collision with a gas stream are being studied using three-dimensional numerical gas-dynamic simulation. The influence of orbital parameters and the stream mass on the accretion activity of the star is examined. It is shown that the orbital inclination and the initial mass of the infalling material are the most influential parameters in determining the accretion rate. The obtained accretion rate dependencies are compared with actual observational data for two FU~Ori type stars. It turns out that not only is the maximum accretion rate consistent with observational estimates, but the behavior of the accretion rate over time is very similar to available long-term light curves.

We investigate how multi-band gravitational wave (GW) observations can constrain the uncertainties in the Hubble parameter ($H_0$) using primordial black holes (PBHs) as possible sources. Our framework combines scalar-induced and merger-induced GWs from PBHs, and forecasts on a combination of two future detectors Square Kilometre Array (SKA) and the Einstein Telescope (ET), enabling a multi-band analysis. We perform a statistical forecast of the PBH parameters, $M_{\rm PBH}$ and $f_{\rm PBH}$, using signal-to-noise ratio (SNR) estimates and Fisher matrix analysis. Imposing $\mathrm{SNR} \geq 1$, we identify the accessible PBH parameter space and propagate these uncertainties to estimate the corresponding uncertainties in $H_0$. For $\delta \theta_i/\theta_i \leq 0.1$, with $\theta_i \equiv M_{\rm PBH}(f_{\rm PBH})$, we find $\delta H_0 \lesssim 2~{\rm km\,s^{-1}\,Mpc^{-1}}$ in a conservative approach, improving to $\delta H_0 \lesssim \mathcal{O}(0.1)~{\rm km\,s^{-1}\,Mpc^{-1}}$ for $\delta \theta_i/\theta_i \leq 0.01$ for an optimistic approach of precision measurement. The results are further found to be largely insensitive to the fiducial choice of the $H_0$, with only moderate dependence on the PBH collapse efficiency. These findings demonstrate that multi-band GW observations provide an independent and complementary approach to constraining the uncertainties in $H_0$, with the potential to provide a novel, cosmic distance ladder-independent measure of the Hubble parameter.

Although the spin parameter of dark matter halos is well known to follow a log-normal distribution at fixed epoch, its quantitative redshift evolution - encompassing both the mean and the dispersion - remains only partially explored. Prior studies either lack the mass resolution required to establish reliable evolutionary trends or do not provide analytical relations that enable forward modelling. Using a suite of LCDM N-body simulations with controlled resolution across the redshift range 0 < z < 5, we characterise the evolution of the mean and dispersion of the Peebles (lambda) and Bullock (lambda') definitions of spin. We find a mild but statistically robust linear evolution for ln(lambda) and a non-monotonic trend with a turnover at z ~ 1 - 2 for ln lambda', which we verify are unaffected by mass resolution of choice of halo definition. We provide closed-form fitting functions for these trends that allow modellers to draw physically motivated spin values at any redshift within our range of validity. This is a practical, redshift-dependent alternative to the common assumption of a constant spin distribution, and provides a useful input to semi-empirical and semi-analytic models of galaxy formation.

swdatatoolkit is a Python-based scientific software library designed to support the acquisition, preprocessing, and analysis of solar and space weather data. The toolkit consolidates functionality across multiple domains, including data downloading from established heliophysics sources, image preprocessing, edge detection, image texture quantification, magnetic field analysis, and the derivation of higher-level parameters commonly used in solar physics research. Its modular structure reflects the heterogeneous nature of space weather data and enables reproducible, extensible workflows for both exploratory analysis and machine-learning-driven studies. This paper presents an overview of the library's available capabilities, its scientific motivations, and its role in the broader space weather research ecosystem.

Recent work on the Loeb Scale has provided astronomy a structured framework for assessing anomalous interstellar objects, including a quantitative mapping of a classification ranking, its evolution with the addition of data, and a broader observational strategy for firming its verdict. What remains unclear is the epistemic and methodological meaning of the threshold built into that framework. Here we argue that the central philosophical issue is no longer whether astronomy can define such a threshold, but how a threshold already in place should regulate scientific inquiry under uncertainty. We suggest that candidate technosignature status, such as Level 4 on the Loeb Scale, should be understood as an intermediate epistemic status: stronger than permissive openness, weaker than confirmation, yet sufficient to justify methodological escalation. The argument proceeds in three steps. First, it reconstructs the recent philosophical debate through the work of Lomas, Lane, and Cowie. Second, it turns to historical cases discussed by Kaplan (2026) to show that important discoveries are often delayed not only by weak evidence, but also by paradigms, prestige, and institutional filtering. Third, it interprets candidate status as a form of structured scientific commitment under uncertainty, one that justifies intensified observation, broader hypothesis management, and more deliberate allocation of attention and resources without licensing belief in artificial origin. The paper concludes by arguing that AI should not be the arbitrator in deducing an extraterrestrial origin, but can support the detection, comparison, and prioritization of anomalies once a candidate status has been formally recognized.

The lack of discovery of particle dark matter candidates within the favored mass-window range brings in the motivation for the study of new options brought by Planck-mass dark matter models. Extended supergravity theories predict the existence of non-relativistic gravitinos that could at least in part contribute to the missing mass-energy density of the Universe. The feasibility study for the discovery with DEAP-3600 experiment of Planck mass charged gravitino dark matter is presented. Additionally the expected signal events topology within the detector is discussed.

Sufficiently strong electric fields can produce charged-particle pairs via the Schwinger effect. We argue that steep matter-density gradients, as can arise in neutron star interiors, would analogously produce neutrino-antineutrino pairs. We then discuss observational signatures of these gradient-produced (anti)neutrinos and how they could provide new probes of neutron-star structure and baryon-dense QCD.

We present a framework in which three classes of dark matter number-changing processes can affect both the relic abundance via thermal freeze-out in the early universe and the generation of indirect cosmic-ray signals today. These processes are: (i) direct annihilations into Standard Model final states; (ii) annihilations into metastable on-shell mediators that subsequently decay into Standard Model particles; (iii) semi-annihilation processes featuring a dark matter particle in the final state, accompanied by a metastable mediator. A central element of our analysis is the systematic inclusion of semi-annihilation alongside the more commonly considered channels. This setup is largely model-independent, as we only assume the presence of one or more of these processes with unsuppressed $s$-wave contributions. We analyze representative benchmarks for the dominant decay modes of the mediator and show how the resulting injection spectra for $\gamma$ rays, neutrinos, and cosmic-ray antimatter vary with the relative importance of the three classes of processes. As an application, we evaluate the observable $\gamma$-ray fluxes from dwarf spheroidal galaxies in the GeV-TeV window. Finally, we provide explicit model realizations in which multiple processes coexist, and discuss how their interplay shapes indirect detection signatures. Our results provide a consistent connection between early-universe dynamics and present-day observables, revealing distinctive features that arise when multiple dark matter processes contribute simultaneously.

Solar-reflected dark matter (SRDM) consists of dark-matter particles up-scattered and accelerated by energetic electrons in the solar interior, producing a high-velocity tail that can enhance signals in direct-detection experiments, especially for MeV-scale masses. We consider an inelastic dark matter (iDM) model, in which solar scattering populates the excited state; subsequent de-excitation in terrestrial detectors releases the mass-splitting energy, substantially helping the energy release of the collision to be larger than the detector threshold. Using detailed Monte Carlo simulations, we generate the velocity and energy distributions of solar-reflected iDM over a range of dark-matter masses $m_\chi$ and mass splittings $\Delta$. We then compute event rates and energy depositions for current xenon and semiconductor experiments. Our results show that these experiments can place new constraints on the parameter space of MeV-scale iDM.

Reducing orbital eccentricity in numerical relativity simulations of binary black holes is essential for producing astrophysically relevant gravitational wave models, as many of these systems are expected to be near-circular in nature. Standard eccentricity reduction procedures rely on iterative schemes, often requiring four or more trial simulations to achieve desired thresholds. This approach is computationally expensive because each trial simulation adds ~10% to the total simulation run time of multiple weeks to months. We introduce a data-driven approach that accelerates this process by learning the values of the initial orbital frequency, Omega_0, and radial velocity, adot_0, that yield an evolution with small eccentricity. This is done using a Gaussian Process Regression model trained on an archive of previously eccentricity-reduced numerical relativity simulations. For all configurations tested, using the trained model consistently reduces the number of required eccentricity reduction iterations to just zero or one, significantly lowering computational costs relative to post-Newtonian initial guesses. These results demonstrate the power of data-driven methods in accelerating expensive numerical relativity simulations.

Gravitational-wave parameter estimation for binary neutron star (BNS) systems poses severe computational challenges due to the extended signal duration, which can reach several minutes in current detectors. Neural posterior estimation (NPE), a simulation-based inference approach, offers dramatic speedups but requires effective dimensionality reduction of the high-dimensional input data. We present a novel compression strategy based on likelihood-oriented summary statistics derived from the relative binning formalism of Zackay et al. (2018), which compresses raw frequency-domain data into the summary data. The summary data is based on a polynomial approximation of the waveform ratio using frequency banding grounded in post-Newtonian approximation, and directly evaluated with only $O(1000)$ sample points of the waveform. As a result, both the training and storage cost become more efficient than previously reported networks for BNS inference. We train a set of NPE networks on these summary statistics and validate a network against traditional nested sampling over 1024 BNS injections. The network produces well-calibrated posteriors across all source parameters we consider, with Jensen-Shannon divergences (JSD) consistent with numerical noise for most parameters. Although we find that the median JSD for the most inconsistent parameter exceeds $10^{-2}$ bits with current configurations, our results show potential for rapid parameter estimation of the BNS signal.

The multi-messenger exploration of dark matter and physics beyond the Standard Model has emerged as a central direction in modern astro-particle physics, particularly following the discovery of gravitational waves. In this work, we present a comprehensive review and forward-looking perspective on machine-learning-enhanced multi-messenger approaches, combining information from gravitational waves, cosmic rays, gamma rays, neutrinos, and collider experiments. We summarize the current state of the field, discuss recent methodological developments, and outline a coherent research program aimed at integrating heterogeneous datasets within a unified inference framework. Our collaboration proposes here a plan for forthcoming analyses aiming at extracting information on the properties and interactions of dark matter, and finally on its genesis, combining multi-messenger astronomy techniques and inputs from laboratory physics. The main objectives planned in this line of research comprise: i) the multi-messenger analysis of new physics in cosmology, including mainly, but not only, several different models of dark matter; ii) the phenomenology of new physics signatures in ground-based cosmic rays experiments, with cross-correlation to the corresponding physical, astrophysical and cosmological observations; iii) the development of machine learning methods for data analysis in ground-based cosmic rays experiments, in light of the new physics signatures. We note that several groups have explored the use of multi-messenger observations, including gravitational waves, to probe alternative dark matter candidates. The present work builds on these developments by focusing on the role of machine learning in integrating heterogeneous datasets. We foresee that such a cross-fertilizing approach will represent the right path to extract information about the main questions left in fundamental physics.