Papers for Friday, Dec 19 2025

In this work we continue our investigations of the possibility of explanation of the positron anomaly (PA) in cosmic rays with the help of annihilating or decaying dark matter (DM) component by varying its space distribution. In the contrast of our previous studies, where we first assumed some specific spatial distribution of DM component and looked at how it agrees with data, here we solve, in some sense, the inverse problem: we search for distribution, in a mathematical way, which satisfies observational data. A unique algorithm has been implemented which, using linear algebra and adaptive grid methods, adjusts distribution to the data. It allows telling in principle whether or not is possible to solve PA problem by variation of spatial distribution of DM sources. A positive result has been formally obtained. A class of solutions can be identified. Though the distributions obtained at the chosen injection spectra may seem slightly realistic, nonetheless it demonstrates a quite powerful possibility in explaining PA that could be realized in more realistic models.

The search for extraterrestrial intelligence (SETI) commensal surveys aim to scan the sky to detect technosignatures from extraterrestrial life. A major challenge in SETI is the effective mitigation of radio frequency interference (RFI), a critical step that is particularly vital for the highly sensitive Five-hundred-meter Aperture Spherical radio Telescope (FAST). While initial RFI mitigation (e.g., removal of persistent and drifting narrowband RFI) are essential, residual RFI often persists, posing significant challenges due to its complex and various nature. In this paper, we propose and apply an improved machine learning approach, the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm, to identify and mitigate residual RFI in FAST-SETI commensal survey archival data from July 2019. After initial RFI mitigation, we successfully identify and remove 36977 residual RFIs (accounting for $\sim$ 77.87\%) within approximately 1.678 seconds using the DBSCAN algorithm. This result shows that we have achieved a 7.44\% higher removal rate than previous machine learning methods, along with a 24.85\% reduction in execution time. We finally find interesting candidate signals consistent with previous studies, and retain one candidate signal following further analysis. Therefore, DBSCAN algorithm can mitigate more residual RFI with higher computational efficiency while preserving the candidate signals that we are interested in.

The formation and evolution of the Milky Way's disc, bar, and bulge remain fundamentally limited by the lack of a contiguous, Galaxy-wide, high-precision chemo-dynamical map. Key open questions - including the survival or destruction of the primitive discs, the origin of the bulge's multi-component structure, the role of mergers and secular processes, and the coupling between stellar chemistry, dynamics, and the Galactic potential - cannot be fully resolved with current or planned facilities. Existing spectroscopic surveys provide either high resolution for small samples or wide coverage at insufficient resolution and depth, and none can obtain homogeneous abundances, precise 3D kinematics, and reliable ages for the millions of stars required, particularly in the obscured midplane, the far side of the bar, or the outer, low-density disc. A new wide-field, massively multiplexed, large-aperture spectroscopic facility, capable of both high- and low-resolution spectroscopy over tens of thousands of square degrees, is therefore essential. Such a facility would deliver the statistical power, sensitivity, and completeness needed to reconstruct the Galaxy's assembly history, constrain its gravitational potential, and establish the Milky Way as the definitive benchmark for galaxy evolution.

White dwarfs, the final evolutionary stage of the vast majority of stars, serve as critical tools for cosmochronology, studies of planetary system evolution, and laboratories for non-standard physics, including exotic cooling channels and weakly interacting particles, as well as crystallization processes. Beyond surface properties accessible via spectroscopy and model atmospheres, global pulsations exhibited by white dwarfs during various evolutionary phases provide a direct window into their deep interiors. Asteroseismology, the comparison of observed pulsation periods with theoretical models, enables us to infer internal chemical stratification, total mass, rotation profiles, and magnetic field strengths. Despite major advances from space missions providing uninterrupted, high-precision photometry, key challenges remain: many predicted pulsators remain quiet, while others oscillate outside theoretical instability strips, highlighting gaps in our understanding of mode excitation, diffusion, and convective mixing. Determining the masses of white dwarfs, particularly for massive and hydrogen-deficient stars, remains uncertain, with discrepancies between spectroscopic, asteroseismic, astrometric, and photometric methods. In the coming decades, large-scale surveys combining high-precision space-based photometry with coordinated ground-based spectroscopic follow-up will dramatically increase both the number and quality of pulsating white dwarf observations.

Very metal-poor massive stars in the Local Group are our best proxies for the Universe's first stars, making them essential for modeling reionization and early galactic chemical evolution. Studying such stars in our Local Universe is key to extrapolating our knowledge to more distant regions, where individual massive stars cannot be resolved but are dynamically and chemically shaping their environments. The MUSE integral field spectrograph has transformed massive star studies in the Milky Way and Magellanic Clouds, but resolving star-forming galaxies containing very metal-poor stars is at the limit of the current field of view and sensitivity. Therefore, only small dedicated efforts of selected regions are studied, providing us with snapshots of low-metallicity massive stars rather than a comprehensive picture. This scarcity is a major bottleneck for understanding and sufficiently modelling the evolution and feedback of massive stars across cosmic time. We therefore envision a new generation of panoramic integral-field spectrographs and high multiplex multi-object spectrographs mounted on dedicated large optical telescopes. Such facilities will not only allow to resolve very-metal-pool galaxies, but further enable the systematic exploration of the massive stellar content across the entire Local Group, thereby reaching a new era in massive star studies and understanding.

Novae are thermonuclear explosions on the surface of accreting white dwarfs and are key laboratories for studying explosive nucleosynthesis, particle acceleration, shock physics, and binary evolution. Despite major progress driven by wide-field time-domain surveys and multi-wavelength facilities, our understanding of nova explosions remains limited by incomplete temporal coverage, heterogeneous spectroscopic follow-up, and poorly constrained ejecta properties. In this white paper we outline the open scientific questions that will define nova research in the 2040s, focusing on the mass, composition, geometry, and dynamics of the ejecta, the role of the underlying binary system, and the connection between nuclear burning, shocks, and emission across the electromagnetic spectrum. We argue that decisive progress requires rapid-response, high-cadence, multi-wavelength observations, anchored by systematic high-resolution optical and near-infrared spectroscopy from eruption to quiescence. Finally, we identify key technological requirements needed to enable transformative advances in the physics of nova explosions over the coming decades.

Accretion onto compact objects is one of the most fundamental phenomena in the astrophysics, powering some of the most luminous objects in the sky. Along with this, accretion has also a key impact on the evolution of the Universe, through the launch of powerful outflows that affect the surrounding medium. In the last years sub-second optical-infrared observations of accreting X-ray binaries have opened a new window in the study of inflow-outflow connection, discovering a wide range of previously unkown phenomena. Here we review the key open questions in accretion and ejection physics and discuss how a dedicated facility, equipped with photon-counting detectors and high spectral resolution from the UV to the mid-infrared, can enable transformative advances in our understanding of accretion processes.

The first few billion years of cosmic history witnessed the rapid emergence of the most massive galaxies, yet their true space density, baryon assembly pathways, and early quenching mechanisms remain poorly constrained. Current surveys lack the wide-field, rest-frame FIR sensitivity needed to obtain a complete census of massive systems and to trace their cold gas, dust, and diffuse emission on galactic and circumgalactic scales. A next-generation facility with a very large aperture, wide field of view, and high mapping speed is essential to carry out deep, degree-scale rest-frame FIR surveys. Such capabilities are required to determine how common massive galaxies are, how they assemble their baryons, and what physical processes drive their early transformation and quenching.

Each galaxy is observed only once along its life, making galaxy evolution fundamentally an inverse statistical problem: time-dependent physics must be inferred from ensembles of single-epoch snapshots. To move beyond descriptive scaling relations toward physical regulation mechanisms of star formation, quenching, chemical enrichment and black hole growth, galaxies must be treated as realizations of a stochastic process whose hyper-parameters (e.g., correlation timescales, burstiness, duty cycles) are inferred hierarchically. This demands both depth and scale: continuum S/N sufficient for absorption-line ages and chemistry, and samples far larger than those in SDSS, DESI, 4MOST or MOONS, which provide either depth or size but not both across $0<z<3$. Once the relevant axes of mass, redshift, environment, structure and evolutionary phase are populated, the requirement naturally rises from $10^7$ to $\sim10^8$ galaxies. This is the regime where stochastic hyper-parameters can be well constrained and where comparisons to simulations and cosmological forward models become limited by theory rather than observations. We outline the science enabled by such a programme and the corresponding requirements for a future ESO wide-field spectroscopic facility capable of delivering tens to hundreds of millions of rest-UV-optical spectra over $0\lesssim z\lesssim3$.

In recent years, JWST has facilitated detections of carbon-bearing molecules in the atmospheres of temperate sub-Neptunes orbiting M dwarfs, ushering in a new era in the characterization of this intriguing planetary regime. We report the transmission spectrum of the temperate sub-Neptune TOI-732 c, observed with JWST NIRISS, NIRSpec G395H and MIRI LRS between 0.9-12 $\mu$m. The observations provide evidence for methane (CH$_4$) in a H$_2$-rich atmosphere, at a volume mixing ratio of $\sim$1\%, and non-detection of NH$_3$ and HCN, along with nominal constraints on other prominent molecules H$_2$O, CO and CO$_2$, which are typically expected in H$_2$-rich atmospheres. We conduct a comprehensive survey of 250 chemical species and find moderate to strong evidence (up to $\ln B\sim 5.9$, $3.9\sigma$) for additional absorption due to one or more complex molecules including higher-order hydrocarbons and/or sulfur-bearing molecules. The spectral features are strongly degenerate among these molecules and with methane, which we find at $\ln B=3.2-8.8$ (up to $3.0-4.6$$\sigma$) significance. Two complex molecules are preferred with at least moderate evidence ($\ln B \gtrsim 2.5$) in both the near- and mid-infrared, while several others show such evidence in at least one of the two wavelength ranges. The preferred molecules are found in trace quantities on Earth, with no significant sources identified in other planetary atmospheres, requiring future work to assess their physical plausibility in this planet. Future observations are required to resolve the degeneracies and place more robust constraints on these species. We highlight the need for further theoretical and experimental work to robustly characterize the atmospheric and internal composition of TOI-732 c and similar sub-Neptunes.

One of JWST's most unexpected discoveries is the emergence of "Little Red Dots'' (LRDs): compact sources at $z \gtrsim 3$ with blue rest-frame UV continua, red optical slopes, and broad Balmer emission lines that challenge standard models and suggest a population of early, unusual active galactic nuclei (AGNs). Using a comprehensive photometric selection and public NIRSpec/PRISM spectroscopy across six JWST deep fields, we identify a large sample of 118 LRDs with high-S/N spectra, enabling a population-wide analysis of their UV--optical continuum and emission lines. We find clear correlations between rest-frame color ([0.3-0.9\,$\mu$m]) and slopes: bluer LRDs have blue UV slopes ($\beta_{\nu,\mathrm{UV}} \sim 0.3$) and red optical slopes, while redder LRDs exhibit redder UV slopes ($\beta_{\nu,\mathrm{UV}} \sim 1.1$). The continuum shape shows a similar trend: redder LRDs display prominent Balmer breaks and curvature, while bluer LRDs follow power-law-like optical SEDs. From literature compilations, $\sim$60% of known broad-line AGNs satisfy our LRD criteria, and up to 90% of LRDs show broad Balmer lines. Emission-line diagnostics reveal a shift from high H$_{\alpha}$/H$_{\beta}$ and low [OIII]$\lambda5007$/H$_{\beta}$ in redder LRDs to the opposite in bluer ones, along with stronger narrow-line equivalent widths, suggesting a transition from AGN- to host-dominated emission. We fit the spectra with a two-component model combining a gas-enshrouded black hole (BH) and a galaxy host. Redder LRDs require higher-luminosity, unreddened BHs and modestly reddened hosts; bluer LRDs require lower-luminosity, reddened BHs and dust-free galaxies. This framework reproduces the diversity in colors and spectral shape by varying BH luminosity, obscuration, and host-to-BH luminosity ratio.

Context. The shapes of dark matter halos can be used to constrain the fundamental properties of dark matter. In standard Cold Dark Matter (CDM) cosmologies, halos are typically triaxial, with a preference for prolate configurations, particularly at low masses and high redshift. Aims. We focus on the characterization of total matter 3D shape in alternative dark matter models, such as Self-Interacting Dark Matter (SIDM) and Warm Dark Matter (WDM). These scenarios predict different structural properties due to collisional effects or the suppression of small-scale power. Methods. We measure the different halo component shapes - dark matter, stars and gas - at various radii from the center in the AIDA-TNG (Alternative Interacting Dark Matter and Astrophysics - TNG), which is a suite of high-resolution cosmological simulations built upon the IllustrisTNG framework. The intent is to systematically study how different dark matter models - specifically, SIDM and WDM - affect galaxy formation and the structure of dark matter halos, when realistic baryonic physics is also included. Results. SIDM models tend to produce rounder and more isotropic halos, especially in the inner regions, as a result of momentum exchange between dark matter particles. WDM halos are also slightly more spherical than their CDM counterparts, and are typically less concentrated. In all cases, the inclusion of self-consistent baryonic physics makes the central regions of all halos rounder, while still revealing clear distinctions among the various dark matter models. Conclusions. The general framework presented in this work, based on the 3D halo shape, can be useful to interpret multi-wavelength data analyses of galaxies and clusters.

Despite numerous search campaigns based on a diverse set of observational techniques, exomoons - prospective satellites of extrasolar planets - remain an elusive and hard-to-pin-down class of objects. Yet, the case for intensifying this search is compelling: as in the Solar System, moons can act as proxies for studying planet formation and evolution, provide direct clues as to the migration history of the planetary hosts and, in favourable cases, offer potentially habitable environments. Here, we present an investigation into how the search for exomoons would benefit from a new interferometric facility operating in the optical wavelength domain and leveraging baselines substantially longer than the ones the VLTI is currently equipped with. We find that an interferometer providing an astrometric precision of 1$\,\mu$as would be able to robustly detect Earth-mass and sub-Earth-mass exomoons on dynamically stable orbits around Jupiter-like planets at distances between 50 and 200 pc.

The nearby star Fomalhaut is orbited by a compact source, Fomalhaut b, which has previously been interpreted as either a dust-enshrouded exoplanet or a dust cloud generated by the collision of two planetesimals. Such collisions are rarely observed but their debris can appear in direct imaging. We report Hubble Space Telescope observations that show the appearance in 2023 of a second point source around Fomalhaut, resembling the appearance of Fomalhaut b twenty years earlier. We interpret this additional source as a dust cloud produced by a recent impact between two planetesimals. The positions and motion of two impact-generated dust clouds over twenty years provide constraints on the collisional dynamics in the debris belt.

In the standard Cold Dark Matter (CDM) scenario, the density profiles of dark matter haloes are well described by analytical models linking their concentration to halo mass. Alternative scenarios, such as warm dark matter (WDM) and self-interacting dark matter (SIDM), modify the inner structure of haloes and predict different profile shapes and central slopes. We employ the AIDA-TNG simulations to investigate how alternative dark matter physics and baryonic processes jointly shape the internal structure of haloes. Using dark-matter-only and full-physics runs, we measure the dark matter density profiles of haloes spanning six orders of magnitude in mass, from 10^9.5 Msun to 10^14.5 Msub, and characterise them with multiple analytical models. We provide the distribution of the best-fitting parameters, as well as the concentration-mass relation in WDM and SIDM. The Einasto profile well reproduces the inner flattening produced in WDM models, both in the collisionless and in the full-physics runs. In SIDM dark-matter-only runs, haloes are better described by explicitly cored profiles, with core sizes that depend on mass and on the self-interaction model. When baryons are included, the differences between CDM and SIDM decrease, and such large dark-matter cores no longer form because adiabatic contraction in the baryon-dominated region counteracts self-interactions. Nevertheless, the coupling between baryons and self-interactions induces a broader range of inner slopes, including cases that are steeper than CDM at Milky Way masses. Alternative dark matter physics thus leaves clear signatures in the inner halo structure, even if baryons significantly reshape these differences. Our results are useful for future studies that need to predict the properties of haloes in multiple dark matter models.

DESI DR2 data have been widely interpreted as evidence for late-time evolving dark energy (DE) with an apparent phantom crossing. Here we investigate an alternative explanation, based on early-Universe physics. If dark acoustic oscillations (DAO) are close in scale to baryon acoustic oscillations (BAO), they can bias the extraction of the BAO scale from the peak in the galaxy correlation function. This leads to an apparent shift in the inferred distance if the superposition of BAO and DAO features is misinterpreted as being due to BAO only. Taking this shift into account, we find that a DAO with percent-level amplitude can reconcile DESI DR2 with Planck 2018 as well as Pantheon+ supernovae data, with fit improvement at a similar level as compared to evolving DE. Notably, a DAO feature with the required properties has been predicted in a previously proposed scenario that resolves the Hubble tension via a pre-recombination decoupling of dark matter and dark radiation (DRMD). The presence of a DAO feature close to the BAO peak can be scrutinized with future full-shape galaxy clustering data from DESI and Euclid.

Long secondary period (LSP) variable stars are a subclass of long-period variables (LPV) that exhibit additional long-term variability alongside pulsations. Despite being observed in over 30% of LPVs, the reason behind the LSP phenomenon is still debated. The most favoured explanation, supported by recent growing evidence, is binarity, where the pulsating giant star has a substellar-mass companion. To further test this hypothesis, it is important to identify bright LSP variables, for which high-quality spectroscopic and interferometric observations can be obtained more easily. Motivated by the absence of a catalog of bright nearby LSPs, we searched the All Sky Automated Survey (ASAS) data in the $V$-band magnitude range 5.5-14 mag, and for declinations $< +28^\circ$. The resulting catalog contains 23 LSPs, 13 of which are new discoveries. We compare our catalog with the LSP lists available in the literature.

Luminous quasars at the redshift frontier z>7 serve as stringent probes of super-massive black hole formation and they are thought to undergo much of their growth obscured by dense gas and dust in their host galaxies. Fully characterizing the symbiotic evolution of SMBHs and hosts requires rest-frame optical observations that span spatial scales from the broad-line region to the ISM and CGM. JWST now provides the necessary spatially resolved spectroscopy to do so. But the physical conditions that regulate the interplay between SMBHs and their hosts at the highest redshifts, especially the nature of early feedback phases, remain unclear. We present JWST/NIRSpec IFU observations of J0313$-$1806 at z=7.64, the most distant luminous quasar known. From the restframe optical spectrum of the unresolved quasar, we derive a black hole mass of $M_\mathrm{BH}=(1.63 \pm 0.10)\times10^9 M_\odot$ based on H$\beta$ and an Eddington rate of $\lambda=L/L_\mathrm{Edd}=0.80\pm 0.05$, consistent with previous MgII-based estimates. J0313-1806 exhibits no detectable [O III] emission on nuclear scales. Most remarkably, we detect an ionized gas shell extending out to $\sim 1.8$ kpc traced by H$\beta$ emission that also lacks any significant [O III], with a $3\sigma$ upper limit on the [O III]$ \lambda$5007 to H$\beta$ flux ratio of $\log_{10} \left( F(\mathrm{[OIII]})/F(\mathrm{H}\beta)\right)=-1.15$. Through photoionization modelling, we demonstrate that the extended emission is consistent with a thin, clumpy outflowing shell where [OIII] is collisionally de-excited by dense gas. We interpret this structure as a fossil remnant of a recent blowout phase, providing evidence for episodic feedback cycles in one of the earliest quasars. These findings suggest that dense ISM phases may play a crucial role in shaping the spectral properties of quasars accross cosmic time.

Warm and self-interactive dark matter cosmologies have been proposed as non-baryonic solutions to the tensions between the $\Lambda$ Cold Dark Matter model and observations at the kpc scale. In this paper, we use the dark matter-only runs of the \textsc{aida-tng} project, a set of cosmological simulations of different sizes and resolutions, to analyze the macroscopic impact of alternative dark matter models on the abundance, the radial distribution and the clustering properties of halos. We adopt the halo occupation distribution formalism to characterize the evolution of its parameters $M_1$ and $\alpha$ with the mass and redshift selection of our sample. By dividing the halo population into central and satellites, we are able to study their spatial density profile, finding that a Navarro-Frenk-White model is not accurate enough to describe the radial distribution of subhalos, and that a generalized Navarro-Frenk-White model is required instead. Warm dark matter models, in particular, present a cuspier distribution of satellites, whereas self-interacting dark matter exhibits a shallower density profile. Moreover, we find that the small-scale clustering of dark matter halos provides a powerful tool to discriminate between alternative dark matter scenarios, in preparation for a more detailed study that fully incorporates baryonic effects, and for a comparison with observational data from galaxy clustering.

We investigate constraints on the inner stellar density profile from photometric data of dwarf spheroidal and ultra-faint dwarf galaxies. Our aim is to clarify under what conditions cored stellar profiles require dark matter halos that are also cored, deviating from the cuspy profiles expected for cold dark matter halos. We consider a variety of spherically symmetric stellar profiles, which we classify as "strong" or "weak" cores and cusps according to the behavior of the slope ($b_0$) and logarithmic slope ($\gamma_0$) at their centers. We explore which profiles lead to unphysical negative distribution functions when embedded in a cuspy halo, treating isotropic and anisotropic kinematics separately. We find that weakly-cored stellar profiles in 3D (i.e., $b_0 \neq 0$, $\gamma_0=0$) can be consistent with cuspy dark matter profiles, but strong 3D cores ($b_0=\gamma_0=0$) are not. However, both weak and strong 3D cores yield nearly indistinguishable inner profiles in projection, which implies that ruling out a dark matter cusp from photometric data alone is highly challenging. As an example, we study the profiles of ultra-faint dwarf galaxies and find that they are consistent with both weak and strong 3D cores. This is not just a result of the limited numbers of stars in these systems, since we reach the same conclusion even for Fornax, one of the most luminous and best-studied dwarf spheroidal companions of the Milky Way. We conclude that, based on current data and analysis techniques, cored surface density profiles in nearby dwarf galaxies cannot be taken as strong evidence against the presence of cuspy dark matter halos.

We present 4MOST-HR resolution Non-Local Thermal Equilibrium (NLTE) Payne artificial neural network (ANN), trained on 404,793 new FGK spectra with 16 elements computed in NLTE. This network will be part of the Stellar Abundances and atmospheric Parameters Pipeline (SAPP), which will analyse 4 million stars during the five year long 4MOST consortium 4: MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey. A fitting algorithm using this ANN is also presented that is able to fully-automatically and self-consistently derive both stellar parameters and elemental abundances. The ANN is validated by fitting 121 observed spectra of low-mass FGKM type stars, including main-sequence dwarf, subgiant and giant stars down to [Fe/H] $\approx -3.4$ degraded to 4MOST-HR resolution, and comparing the derived abundances with the output of the classical radiative transfer code TSFitPy. We are able to recover all 18 elemental abundances with a bias <0.13 and spread <0.16 dex, although the typical values are <0.09 dex for most elements. These abundances are compared to the OMEGA+ Galactic Chemical Evolution model, showcasing for the first time, the expected performance and results obtained from high-resolution spectra of the quality expected to be obtained with 4MOST. The expected Galactic trends are recovered, and we highlight the potential of using many chemical elements to constrain the formation history of the Galaxy.

Wide-field, multi-band surveys now detect millions of unresolved sources in nearby galaxy clusters, yet separating globular clusters (GCs) from foreground stars and background galaxies remains challenging. Scalable, automated classification is therefore essential to convert the forthcoming data from facilities such as the Vera C. Rubin/LSST, the Roman and Euclid into robust constraints on galaxy assembly. We introduce a supervised classification method to separate GCs, stars, and galaxies based on their locations in color-color diagrams. The main objective is to recover a clean GC sample for future scientific analysis. The method exploits broad spectral energy distribution coverage, deep photometry, and is optimized for next-generation survey volumes. We use the central 3deg2 of the Next Generation Fornax Survey (NGFS), which images the Fornax cluster in u'g'i'JKs. We build a Support Vector Machine (SVM; this http URL, scikit-learn) using 15 features: all color combinations and basic morphological parameters. Spectroscopically confirmed sources define the training classes. Color pairs connecting near-UV/optical/near-IR. The full 15 feature model achieves 97.3% accuracy and a pruned 7 feature model built from the most informative, least correlated features achieves 96.6% accuracy. Misclassifications amount 8.4% and 10.4%, respectively. Omitting the u' or/and near-IR bands degrades performance. Emulating LSST filters with NGFS u'g'i' and DES r'z'Y shows that u' and Y bands are crucial, but models lacking NIR remain suboptimal. Combining broad SED coverage with simple morphological parameters enables precise, scalable separation of unresolved sources. Including NIR bands significantly improves GC classification, and joining LSST with forthcoming Euclid and Roman data will further enhance machine-learning frameworks.

For many decades, dust has been recognised as an important ingredient in galaxy formation and evolution. This paper presents a novel self-consistent implementation of dust formation by stars, destruction by supernova shocks and hot gas, and growth within the dense interstellar medium (ISM) in the GAEA state-of-the-art galaxy formation model. Our new model, DUSTY-GAEA, reproduces well the dust buildup as a function of stellar mass out to z $\sim$ 6, the scaling relations between the dust-to-gas/dust-to-metal ratios and stellar mass/metallicty in the local Universe, and the dust mass function in the local Universe and out to z $\sim$ 1. In the framework of our model, dust growth dominates the cosmic dust budget out to z $\sim$ 8, and we find that observational constraints beyond the local Universe can be reproduced only assuming such efficient dust growth in the dense ISM. Yet, reproducing the estimated number densities of dust-rich galaxies at higher redshifts remains challenging, as found also in independent theoretical work. We discuss our model predictions in comparison with both observational data and independent theoretical efforts, and highlight how further observational constraints at high redshifts would help constrain dust models.

Stellar multiplicity is a fundamental ingredient of stellar astrophysics, yet binary statistics across the Galaxy remain poorly constrained. The \emph{Gaia} mission has revolutionised binary star astrophysics by delivering high-precision astrometry, photometry and global radial velocities, and by providing hundreds of thousands of non-single-star solutions in DR3. However, the RVS magnitude limit, mission time span and scanning law impose strong selection effects in period, mass ratio, inclination and semi-amplitude, leaving large regions of the binary parameter space either sparsely sampled or effectively inaccessible. In this white paper we outline the case for a dedicated, wide-field, multi-epoch spectroscopic survey explicitly optimised for binary science: deeper than the \emph{Gaia} RVS limit, with flexible cadence from hours to years, and with moderate to high spectral resolution. Using a simplified forward model of \emph{Gaia} DR5-like performance, we highlight the populations for which robust orbital solutions will be rare (ultra short period, very long period, low-amplitude and compact-object binaries), and show how a ``Binarity Beyond \emph{Gaia}'' survey would fill these gaps. Such a programme would deliver a bias correctable census of stellar multiplicity across the Milky Way and provide the spectroscopic backbone needed to exploit binary samples from \emph{Rubin}/LSST, \emph{Roman} and \emph{LISA}.

Here we summarize discussions and conclusions from the conference ``The Era of Binary Supermassive Black Holes: Coordination of Nanohertz-Frequency Gravitational-Wave Follow-up,'' held at the Aspen Center for Physics from February 2-7, 2025. The meeting facilitated a crucial knowledge exchange between electromagnetic and gravitational-wave theorists, observers, and cyber-infrastructure experts. The central goal was to guide the development of multi-messenger follow-up strategies for binary supermassive black hole detections by pulsar timing arrays. To build a common basis of understanding for the broader scientific community, this summary outlines the main considerations and recommendations from the meeting, summarizes the knowledge gaps identified, and ends with a potential roadmap to catalyze discussion about the search for electromagnetic counterparts to massive black hole binaries detected by pulsar timing arrays.

Fast reconnection in magnetically dominated plasmas is widely invoked in models of dissipation in pulsar winds, gamma-ray flares in the Crab nebula, and to explain the radio nanoshots of pulsars. When current sheets evolve reaching a critical inverse aspect ratio, scaling as $S^{-1/3}$ with the plasma Lundquist number, the so-called \textit{ideal} tearing instability sets in, with modes growing, independently of $S$, extremely rapidly on timescales of only a few light-crossing times of the sheet length. We present the first set of fully 3D simulations of current-sheet disruption triggered by the ideal tearing instability within the resistive relativistic MHD approximation, as appropriate in situations where the Alfvén velocity approaches the speed of light. We compare 3D setups with different initial conditions with their 2D counterparts, and we assess the impact of dimensionality and of the magnetic field topology on the onset, evolution, and efficiency of reconnection. In force-free configurations, 3D runs develop ideal tearing, secondary instabilities, and a thick, turbulent current layer, sustaining dissipation of magnetic energy longer than in 2D. In pressure-balanced current sheets with a null guide field, 2D reference runs show the familiar reconnection dynamics, whereas in 3D tearing dynamics is quenched after the linear phase, as pressure-driven modes growing on forming plasmoids outcompete plasmoid coalescence and suppress fast dissipation of magnetic energy. Taken together, these results suggest that the evolution and efficiency of reconnection depend sensitively on the local plasma conditions and current-sheet configuration, and can be properly captured only in fully 3D simulations.

We present a joint cosmological analysis of projected galaxy clustering observations from the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1), and overlapping weak gravitational lensing observations from three datasets: the Kilo-Degree Survey (KiDS-1000), the Dark Energy Survey (DES-Y3), and the Hyper-Suprime-Cam Survey (HSC-Y3). This combination of large-scale structure probes allows us to measure a set of $3 \times 2$-pt correlation functions, breaking the degeneracies between parameters in cosmological fits to individual observables. We obtain mutually-consistent constraints on the parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.786^{+0.022}_{-0.019}$ from the combination of DESI-DR1 and DES-Y3, $S_8 = 0.760^{+0.020}_{-0.018}$ from KiDS-1000, and $S_8 = 0.771^{+0.026}_{-0.027}$ from HSC-Y3. These parameter determinations are consistent with fits to the Planck Cosmic Microwave Background dataset, albeit with $1.5-2\sigma$ lower values in the $S_8-\Omega_{\rm m}$ plane. We perform our analysis with a unified pipeline tailored to the requirements of each cosmic shear survey, which self-consistently determines cosmological and astrophysical parameters. We generate an analytical covariance matrix for the correlation data including all cross-covariances between probes, and we design a new blinding procedure to safeguard our analysis against confirmation bias, whilst leaving goodness-of-fit statistics unchanged. Our study is part of a suite of papers that present joint cosmological analyses of DESI-DR1 and weak gravitational lensing datasets.

We present a joint $3\times2$-pt cosmological analysis of auto- and cross-correlations between the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1) Bright Galaxy Survey (BGS) and Luminous Red Galaxy (LRG) samples and overlapping shear measurements from the KiDS-1000, DES-Y3 and HSC-Y3 weak lensing surveys. We perform our analysis in configuration space and, in addition to the cosmic shear correlation functions for each weak lensing dataset, we fit the tangential shear of the weak lensing source galaxies around DESI lens galaxies. Finally, we make use of the anisotropic BGS and LRG clustering information by fitting the full shape of the two-point correlation function multipoles measured over the full DESI-DR1 footprint, presenting the first full-shape analysis of DESI measurements in configuration space. We find that the addition of weak lensing information serves to improve, with respect to the clustering-only case, the measurements of the power spectrum amplitude parameters $\ln(10^{10}A_{\rm{s}})$ and $\sigma_{12}$ by $15\%$ and $36\%$, respectively. It also improves measurements of the linear bias of the lens galaxies by $15-20\%$, depending on the tracer. Our results show excellent consistency, regardless of the weak lensing survey considered, and are furthermore consistent with a companion analysis that fits $3\times2$-pt correlations including DESI projected clustering measurements, as well as the results published by the weak lensing collaborations themselves. Our measured values for weak lensing amplitude are $S_{8}^{\mathrm{DESI\times HSC}}=0.787\pm0.020$, $S_{8}^{\mathrm{DESI\times DES}}=0.791\pm0.016$, $S_{8}^{\mathrm{DESI\times KiDS}}=0.771\pm0.017$, which are $1.9\sigma-2.9\sigma$ below the $S_8$ value preferred by Planck. Finally, our clustering-only results are in good agreement with the Fourier space full-shape analysis of all DESI tracers.

We present constraints on cosmic structure growth from the analysis of galaxy clustering and galaxy-galaxy lensing with galaxies from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1. We analyze four samples drawn from the Bright Galaxy Survey (BGS) and the Luminous Red Galaxy (LRG) target classes. Projected galaxy clustering measurements from DESI are supplemented with lensing measurements from the Dark Energy Survey (DES), the Kilo-Degree Survey (KiDS), and the Hyper Suprime-Cam (HSC) survey around the same targets. Our method relies on a simulation-based modeling framework using the AbacusSummit simulations and a complex halo occupation distribution model that incorporates assembly bias. We analyze scales down to $0.4 \, h^{-1} \, \mathrm{Mpc}$ for clustering and $2.5 \, h^{-1} \, \mathrm{Mpc}$ for lensing, leading to stringent constraints on $S_8 = \sigma_8 \sqrt{\Omega_\mathrm{m} / 0.3}$ and $\Omega_\mathrm{m}$ when fixing other cosmological parameters to those preferred by the CMB. We find $S_8 = 0.794 \pm 0.023$ and $\Omega_\mathrm{m} = 0.295 \pm 0.012$ when using lensing measurements from DES and KiDS. Similarly, for HSC, we find $S_8 = 0.793 \pm 0.017$ and $\Omega_\mathrm{m} = 0.303 \pm 0.010$ when assuming the best-fit photometric redshift offset suggested by the HSC collaboration. Overall, our results are in good agreement with other results in the literature while continuing to highlight the constraining power of non-linear scales.

We estimate the source redshift distribution of current weak lensing surveys by applying the clustering-based redshift calibration technique, using the galaxy redshift sample provided by the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1). We cross-correlate the Bright Galaxy Survey (BGS), Luminous Red Galaxies (LRGs) and Emission Line Galaxies (ELGs) from DESI, within the redshift range $0.1 < z < 1.6$, with overlapping tomographic source samples from the Dark Energy Survey (DES), Kilo-Degree Survey (KiDS), and Hyper Suprime-Cam (HSC) survey. Using realistic mock catalogues, we test the stability of the clustering-redshift signal to fitting scale, reference-sample choice, and the evolution of source galaxy bias, and we explicitly model and marginalise over magnification contributions, which become non-negligible at $z \gtrsim 1$ due to the depth of the DESI ELG sample. We then compare the resulting bias-weighted redshift distributions to those calibrated using self-organising map (SOM) techniques, finding agreement within uncertainties for all surveys and tomographic bins. Our results demonstrate that clustering redshifts enabled by DESI's unprecedented spectroscopic sample provides a robust, complementary, and independent constraint capable of reducing one of the dominant systematic uncertainties in weak lensing cosmology.

The origins of the Martian moons Phobos and Deimos are highly debated, and hypotheses include formation from an impact-generated circum-Martian disk or from capture of asteroids. With the impact scenario, Deimos (or its precursors) were formed or were pushed out beyond the synchronous orbit of Mars. Moons interior to the synchronous orbit, including Phobos (or its precursors), would tidally evolve and resonances between these moons could potentially excite Deimos' orbit. This contradicts Deimos' present-day orbit of low eccentricity ($0.00027$) and moderate inclination ($1.8^\circ$ to the Laplace plane). Tidal dissipation within Deimos is too inefficient for eccentricity damping, and without alternative mechanisms, Deimos' present-day orbit places strong constraints on the evolution of any inner moons. We propose that a runaway collisional cascade called the "sesquinary catastrophe'' acts as a natural barrier that prevents Deimos from having a more excited orbit. Using N-body simulations with collisional fragmentation, we show that if Deimos was more excited, it would undergo a sesquinary catastrophe and break apart into a Roche-exterior debris disk. Using a measure of sesquinary orbital excitation called $q$, our simulations and previous works suggest that breakup occurs for $q \gtrsim 8$ on timescales of $\sim 10^{3-4}$ years. If Deimos was destroyed in a sesquinary catastrophe and re-accreted from a (likely collisionally) damped debris disk, it should be a porous sand-pile moon, consistent with its smooth surface. The sesquinary catastrophe can be applied to other Deimos-like planetary moons at $q \gtrsim 8$.

The Large Magellanic Cloud (LMC) provides a key laboratory for exploring the diversity of star formation and interstellar chemistry under subsolar metallicity conditions. We present the results of a hot core survey toward 30 massive protostellar objects in the LMC using the Atacama Large Millimeter/submillimeter Array (ALMA) at 350 GHz. Continuum imaging reveals 36 compact sources in total, among which line analyses identify 9 hot cores and 1 hot-core candidate, including two newly identified sources. We detect CO, HCO+, H13CO+, HC15N, HC3N, SiO, SO, SO+, NS, SO2, 34SO2, 33SO2, CH3OH, 13CH3OH, HCOOH, HCOOCH3, CH3OCH3, C2H5OH, H2CCO (tentative), and hydrogen recombination lines from hot cores. CH3OCH3, a complex organic molecule larger than CH3OH, is detected for the first time in a hot core outside the LMC bar region. All hot cores show stronger emission in the high-excitation SO line compared to non-hot-core sources, suggesting that its strong detection will be useful for identifying hot-core candidates in the LMC. Chemical analysis reveals a spread of more than two orders of magnitude in CH3OH abundances, with some sources deficient in COMs. In contrast, SO2 is detected in all hot cores, and its abundance shows a good correlation with rotational temperature. The hot cores without CH3OH detections are all located outside the LMC bar region and are characterized by either high luminosity or active star formation in their surroundings. A combination of locally low metallicity, active star formation in the vicinity, and high protostellar luminosity may jointly trigger the COM-poor hot core chemistry observed in the LMC.

We present the discovery of a large-scale, limb-brightened outflow, extending at least 30 kpc above and below the star-forming disk of the edge-on galaxy ESO 130-G012 (D = 16.9 Mpc). Partially obscured by Galactic foreground stars and dust, this optically unremarkable, low-mass galaxy reveals one of the largest known hourglass-shaped outflows from the full extent of its bright stellar disk. The outflow was discovered in 944 MHz radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) obtained as part of the "Evolutionary Map of the Universe" (EMU) project. Its height is at least 3x that of the stellar disk diameter (~10 kpc), while its shape and size most resemble the large biconical, edge-brightened FUV and X-ray outflows in the nearby starburst galaxy NGC 3079. The large-scale, hourglass-shaped outflow of ESO 130-G012 appears to be hollow and originates from the star-forming disk, expanding into the halo with speeds close to the escape velocity before likely returning to the disk. Given ESO 130-G012's modest star formation rate, the height of the outflow is surprising and unusual, likely made possible by the galaxy's relatively low gravitational potential. Follow-up observations are expected to detect hot gas inside the bipolar outflow cones and magnetic fields along the X-shaped outflow wings. Neutral gas may also be lifted above the inner disk by the outflow.

The earliest stages of star formation are highlighted by complex interactions between accretion, outflow, and radiative processes, which shape the chemical and physical environment of the emerging protostar. James Webb Space Telescope observations of the low-mass, low-luminosity Class 0 protostar IRAS 16253-2429 reveal a central compact source. This object exhibits a rich mid-infrared emission spectrum of OH pure rotational lines and $\rm CO_2$ ro-vibrational lines. Unusually for a young stellar object, it has no mid-infrared line emission from $\rm H_2O$ to match the other molecules. We demonstrate that the emitting OH molecules arise from UV photodissociation of $\rm H_2O$ in its second absorption band at $\lambda = 114-145$ nm, and that the OH emission is a fluorescent cascade starting with highest-excitation rotational states. This situation offers the opportunity of using the infrared OH spectrum to measure the UV flux from the central protostar. Thereby we determine the disk-star accretion rate to be $3 \times 10^{-10} \ M_\sun \ {\rm year^{-1}}$, and demonstrate that the system luminosity arises mostly from the protostar's photosphere rather than from accretion luminosity. The result is in accord with the measured outflow rate of IRAS 16253-2429 and lies within the outflow/accretion-flow rate trend often inferred for protostars; and with episodic accretion as the dominant mechanism by which this protostar has grown.

During its first phase, from 2004 up to the end of 2012, the H.E.S.S. (High Energy Stereoscopic System) experiment observed the extragalactic skies for more than 2700 hours. These data have been re-analysed in a single consistent framework, leading to the derivation of a catalog of 23 sources. In total, about 5.7% of the sky was observed, allowing for several additional studies to be conducted: source variability, extragalactic gamma-ray background light, and comparison with the Fermi-LAT catalogues. In this contribution, we discuss these results and present the high-level data (catalogs, maps) released to the astrophysical community.

The HESS Galactic Plane Survey (HGPS), published in 2018, presented a decade of very-high-energy (VHE) gamma-ray observations along the Galactic plane. This study was accompanied by the release of several maps in FITS format, offering a detailed view of the region. The flux upper-limits from these HGPS maps can be compared to the high-energy (HE) spectra of sources catalogued by the Fermi-LAT in the same region. For some sources, extrapolating the Fermi-LAT flux into the VHE range predicts flux values exceeding the upper-limits set by HESS. In this work, we present the results of this comparison and highlight the sources that are of particular interest for future VHE observations.

Quasar photometric redshifts are essential for studying cosmology and large-scale structures. However, their complex spectral energy distributions cause significant redshift-color degeneracy, limiting the accuracy of traditional methods. To overcome this, we introduce LSTM-MDNz, a novel end-to-end deep learning model combining long short-term memory networks (LSTM) with mixture density networks (MDN). The model directly uses multi-band photometric fluxes and associated errors as wavelength-ordered sequential inputs, eliminating the need for manual feature engineering while enabling simultaneous point estimation and probability distribution function (PDF) prediction of quasar redshifts. We integrate data from four major sky surveys-SDSS, DESI-LS, WISE, and GALEX-to assemble a sample of over 550,000 spectroscopically confirmed quasars ($0 \leq z_{\mathrm{spec}} \leq 5$) across 14 ultraviolet to infrared bands for model training and testing. Experimental results show that using all 14 bands yields optimal performance, with a normalized median absolute deviation ($\sigma_{\mathrm{NMAD}}$) of 0.037 and an outlier rate ($f_{\mathrm{out}}$) of 3.5\% on the test set. These values represent reductions of 29\% and 56\%, respectively, compared to the commonly adopted SDSS+WISE band set. Probability integral transform ($\mathrm{PIT}$) and continuous ranked probability score ($\mathrm{CRPS}$) analyses confirm that the predicted PDFs align closely with the true redshift distribution. Band-ablation experiments further highlight the essential role of ultraviolet and infrared data in alleviating color degeneracy and reducing systematic bias. This study demonstrates the effectiveness of multi-band fusion in improving quasar photo-z accuracy and offers a ready-to-use estimation framework for future surveys like LSST, CSST, and Euclid.

Steady-state solutions to the Navier-Stokes equations are known to admit solutions that are singular at the sonic point. Consequently, inviscid solar wind models require special treatment of the solution near the sonic points, and this has proven to be a significant impediment to efficient modeling of the solar wind. In this paper we revisit the governing hydrodynamic equations for the expanding solar wind, with the inclusion of the classical (Newtonian) viscous stress , and we show how this inclusion eliminates the singularities that emerge from the inviscid equations. This result has been previously reported and used to generate solar wind profiles from initial conditions in the asymptotic limit; however, those studies did not include realistic treatments of the inner corona, and generally rejected the prospect of extrapolating solutions outward from the Sun into the heliosphere. Here, we expand this method to include external heating and optically thin radiative losses and show that solutions can be computed from initial conditions near the solar surface, thereby capturing the entire range of scales from below the transition region to the outer heliosphere in a single solution. Our approach is to cast the steady-state Navier-Stokes equations as a system of five coupled, ordinary differential equations (ODEs), which we solve using conventional methods, without any special treatment of the governing equations in the vicinity of the sonic point. The representative solutions that we present here demonstrate the utility and efficiency of this extrapolation method, which is considerably more realistic than commonly used analytical or empirical models. This method provides a direct approach to generating accurate solar wind profiles subject to observationally motivated initial conditions near the solar surface, at a fraction of the computational cost of comparable relaxation-based models.

Open clusters serve as laboratories to study and evaluate stellar evolution and Galactic chemical evolution models. Chemical peculiarities, such as lithium-rich giants, are rarely observed in these stellar systems. This work focuses on eight red giants (#005, #028, #034, #053, #087, #121, #126, and #190) previously reported as members of the Galactic cluster IC 2714. We conducted a detailed investigation using high-resolution spectroscopy, supplemented with data from the Gaia DR3 catalog. Besides deriving the cluster's fundamental parameters, we provide the most thorough chemical characterization of IC 2714 to date, reporting the abundance of 23 species, including light elements (Li, C, N, O), odd-Z elements (Na, Al), $\alpha$-elements (Si, Ca, Ti, Mg), iron-peak elements (Sc, Cr, Ni), $\textit{s}$-process-dominated elements (Y, Zr, Ba, La, Ce, Nd) and $\textit{r}$-process elements (Sm, Eu). We also present the carbon isotopic ratios $^{12}$C/$^{13}$C for the first time for seven stars. One particular star (#087) exhibits a high lithium abundance ($\log \varepsilon$(Li)$_{\rm NLTE}$ = $+$1.54 dex) and a slightly higher projected rotational velocity ($v \sin i$ = 6.7 km s$^{-1}$). Our results suggest that the analyzed stars are in the core-helium-burning phase of evolution, where the most lithium-rich giants are found. Combining astrometric probabilities and chemical abundances, we conclude that two giants (#028 and #034) might not be cluster members.

How a seemingly `dead' host galaxy provides fuel for its active galactic nuclei (AGN) remains an unresolved problem. Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we present a new high-sensitivity atomic-hydrogen (HI) observation toward the nearby elliptical galaxy NGC 4278 and its adjacent region. From the observation, we found that external gas accretion from a neighbouring galaxy fuels the low-luminosity AGN in NGC 4278 through tidal interactions. The accreted gas entering NGC 4278 exhibits a rotating gas disk. And the accreted galaxy has been gas-poor and has an HI to stellar mass ratio of about 0.02. Due to the process of gas accretion, it is likely that relativistic jets are generated in the AGN of NGC 4278. The emission of TeV gamma rays in NGC 4278 is likely to be associated with the newly accreted HI gas.

In recent years, graph neural networks (GNNs) have shown tremendous promise in solving problems in high energy physics, materials science, and fluid dynamics. In this work, we introduce a new application for GNNs in the physical sciences: instrumentation design. As a case study, we apply GNNs to simulate models of the Laser Interferometer Gravitational-Wave Observatory (LIGO) and show that they are capable of accurately capturing the complex optical physics at play, while achieving runtimes 815 times faster than state of the art simulation packages. We discuss the unique challenges this problem provides for machine learning models. In addition, we provide a dataset of high-fidelity optical physics simulations for three interferometer topologies, which can be used as a benchmarking suite for future work in this direction.

High-resolution observations of Class II protoplanetary disks frequently reveal annular structures that may indicate the presence of embedded planets. In this study, we model brightness profiles and geometries for 16 dusty disks in the DSHARP observations to identify asymmetric substructures, including possible planet-induced signatures. We find no compelling evidence for circumplanetary emission in these systems. We identify a possible one-armed spiral in Elias 27; while previous studies report an $m=2$ spiral, the morphology of the newly identified spiral agrees with a spiral for a possible protoplanet. Although non-detection of circumplanetary emission is consistent with low expected luminosities, the absence of additional dust spirals except for Elias 27 may constrain properties of potential embedded planets given their theoretical detectability. The analyses further suggest that spiral amplitudes and phases correlated with gap and ring locations in WaOph 6 and IM Lup, appearing as deflections of spirals in images. Five disks exhibit strong residual asymmetries attributable to the vertical extent of their dust layers. We find that the brightness temperature of inner disks ($1$-$5$ au) declines with stellar age on a $\sim 3$ Myr timescale, while the total flux shows no clear decreasing trend. This trend is consistent with the presence of pressure bumps that retain large dust grains at outer radii, while allowing the inner disks to disperse. The universality of disk dispersal timescale at millimeter wavelengths, observed in both the extended DSHARP disks and the compact disks from previous demographic surveys, may constrain timing of planet formation, including habitable planets.

The statistical analysis of fast radio burst (FRB) samples from repeaters may suffer from a band-limited selection effect, which can bias the observed distribution. We investigated the impact of this selection bias on the energy function through simulations and then applied our analysis to the particular case of FRB 20220912A. Our simulations show that, in the sample of bursts observed by the Five hundred meter Aperture Spherical Telescope (FAST), assuming a unimodal intrinsic energy distribution, the band selection effect alone is insufficient to produce a bimodal energy distribution; only the bimodal central frequency distribution can achieve this. The bursts' energy of FRB 20220912A that primarily fell within the observing band showed no significant correlation with the central frequency. In contrast, bursts with higher central frequency tend to exhibit narrower bandwidth and longer duration. The distribution of the intrinsic energy can be modeled as a log-normal distribution with a characteristic energy of $8.13 \times 10^{37}$ erg, and a power-law function with the index of $1.011 \pm 0.028$. In contrast to the initial energy function reported by \cite{2023ApJ...955..142Z}, the low-energy peak vanishes, and the high-energy decline becomes steeper, which implies the low-energy peak is an observational effect. The bimodality of the energy distribution seems to originate from the intrinsic radiation mechanism.

The large instantaneous sensitivity, a wide frequency coverage and flexible observation modes with large number of beams in the sky are the main features of the SKA observatory's two telescopes, the SKA-Low and the SKA-Mid, which are located on two different continents. Owing to these capabilities, the SKAO telescopes are going to be a game-changer for radio astronomy in general and pulsar astronomy in particular. The eleven articles in this special issue on pulsar science with the SKA Observatory describe its impact on different areas of pulsar science. In this lead article, a brief description of the two telescopes highlighting the relevant features for pulsar science is presented followed by an overview of each accompanying article, exploring the inter-relationship between different pulsar science use cases.

Most of the pulsar science case with the Square Kilometre Array (SKA) depends on long-term precision pulsar timing of a large number of pulsars, as well as astrometric measurements of these using very long baseline interferometry (VLBI). But before we can time them, or VLBI them, we must first find them. Here, we describe the considerations and strategies one needs to account for when planning an all-sky blind pulsar survey using the SKA. Based on our understanding of the pulsar population, the performance of the now-under-construction SKA elements, and practical constraints such as evading radio frequency interference, we project pulsar survey yields using two complementary methods for a number of illustrative survey designs, combining SKA1-Low and SKA1-Mid Bands 1 and 2 in a variety of ways. A composite survey using both Mid and Low is optimal, with Mid Band 2 focused in the plane. We find that, given its much higher effective area and survey speed, the best strategy is to use SKA1-Low to cover as much sky as possible, ideally also overlapping with the areas covered by Mid. In our most realistic scenario, we find that an all-sky blind survey with Phase 1 of the SKA with the AA* array assembly will detect $\sim10,000$ slow pulsars and $\sim 800$ millisecond pulsars (MSPs) if SKA1-Mid covers the region within $5°$ of the plane, while higher latitudes will be covered with SKA1-Low. The yield with AA4 is $\sim 20\%$ higher. One could increase these numbers by increasing the range covered by SKA1-Mid Bands 1 and 2, at the cost of a considerably longer survey. The pulsar census will enable us to set new constraints on the uncertain physical properties of the entire neutron star population. This will be crucial for addressing major SKA science questions including the dense-matter equation of state, strong-field gravity tests, and gravitational wave astronomy.

Because of their extreme stellar densities, globular clusters are highly efficient factories of X-ray binaries and radio pulsars: per unit of stellar mass, they contain about 1000 times more of these exotic objects. Thus far, 345 radio pulsars have been found in globular clusters. These can be used as precision probes of the structure, gas content, magnetic field, and dynamic history of their host clusters; some of them are also highly interesting in their own right because they probe exotic stellar evolution scenarios as well as the physics of dense matter, accretion, and gravity; one of them (PSR~J0514$-$4002E) might even be the first pulsar - black hole system known. Deep searches with SKA-MID and SKA-LOW will only require one to a few tied-array beams, and can be done during early commissioning of the telescope, before an all-sky pulsar survey using hundreds to thousands of tied-array beams is feasible. Even a conservative approach predicts new discoveries only with the core of SKA-MID AA*, and the full AA* and eventually AA4 is expected to increase the number of discoveries even more, leading to more than doubling the current known population. This offers a great opportunity for early SKAO pulsar science, even before all the collecting area is in place. On the other hand, a more optimistic prediction calls for a 4-5 times growth of the population, leading to a total of about 1700 pulsars to be detectable with SKA-MID AA4 configuration in all Galactic GCs visible by SKA telescopes. Thus, a dedicated search for pulsars in globular clusters will fully exploit the best possible natural laboratories to study many branches of physics and astrophysics, including properties of dense matter, stellar evolution, and the dynamical history of the Galactic globular cluster systems.

The detection of a pulsar closely orbiting our Galaxy's supermassive black hole - Sagittarius A* - is one of the ultimate prizes in pulsar astrophysics. The relativistic effects expected in such a system could far exceed those currently observable in compact binaries such as double neutron stars and pulsar white dwarfs. In addition, pulsars offer the opportunity to study the magneto-ionic properties of Earth's nearest galactic nucleus in unprecedented detail. For these reasons, and more, a multitude of pulsar searches of the Galactic Centre have been undertaken, with the outcome of just seven pulsar detections within a projected distance of 100 pc from Sagittarius A*. It is currently understood that a larger underlying population likely exists, but it is not until observations with the SKA have started that this population can be revealed. In this paper, we look at important updates since the publication of the last SKAO science book and offer a focused view of observing strategies and likely outcomes with the updated SKAO design.

The known population of non-accreting neutron stars is ever growing and currently consists of more than 3500 sources. Pulsar surveys with the SKAO telescopes will greatly increase the known population, adding radio pulsars to every subgroup in the radio-loud neutron star family. These discoveries will not only add to the current understanding of neutron star physics by increasing the sample of sources that can be studied, but will undoubtedly also uncover previously unknown types of sources that will challenge our theories of a wide range of physical phenomena. A broad variety of scientific studies will be made possible by a significantly increased known population of neutron stars, unravelling questions such as: How do isolated pulsars evolve with time; What is the connection between magnetars, high B-field pulsars, and the newly discovered long-period pulsars; How is a pulsar's spin-down related to its radio emission; What is the nuclear equation of state? Increasing the known numbers of pulsars in binary or triple systems may enable both larger numbers and higher precision tests of gravitational theories and general relativity, as well as probing the neutron star mass distribution. The excellent sensitivity of the SKAO telescopes combined with the wide field of view, large numbers of simultaneous tied-array beams that will be searched in real time, wide range of observing frequencies, and the ability to form multiple sub-arrays will make the SKAO an excellent facility to undertake a wide range of neutron star research. In this paper, we give an overview of different types of neutron stars and discuss how the SKAO telescopes will aid in our understanding of the neutron star population.

The SKA telescopes will bring unparalleled sensitivity across a broad radio band, a wide field of view across the Southern sky, and the capacity for sub-arraying, all of which make them the ideal instruments for studying the pulsar magnetosphere. This paper describes the advances that have been made in pulsar magnetosphere physics over the last decade, and details how these have been made possible through the advances of modern radio telescopes, particularly SKA precursors and pathfinders. It explains how the SKA telescopes would transform the field of pulsar magnetosphere physics through a combination of large-scale monitoring surveys and in-depth follow-up observations of unique sources and new discoveries. Finally, it describes how the specific observing opportunities available with the AA* and AA4 configurations will achieve the advances necessary to solve the problem of pulsar radio emission physics in the coming years.

The ionised media that permeate the Milky Way have been active topics of research since the discovery of pulsars in 1967. In fact, pulsars allow one to study several aspects of said plasma, such as their column density, turbulence, scattering measures, and discrete, intervening structures between the neutron star and the observer, as well as aspects of the magnetic field throughout. Such sources of information allow us to characterise the electron distribution in the terrestrial ionosphere, the Solar Wind, and our Galaxy, as well as the impact on other experiments involving pulsars, such as Pulsar Timing Arrays. In this article, we review the state-of-the-art in plasma research using pulsars, the aspects that should be taken into consideration for optimal plasma studies, and we provide future perspectives on improvements enabled by the SKA.

The magnetic field structure of the Milky Way can offer critical insights into the origin of galactic magnetic fields. Measurements of magnetic structures of the Milky Way are still sparse in far regions of the Galactic disk and halo. Pulsars are the best probes for the three-dimensional structure of the Galactic magnetic field, primarily owing to their highly polarized short-duration radio pulses, negligible intrinsic Faraday rotation compared to the contribution from the medium in front, and their widespread distribution throughout the Galaxy across the thin disk, spiral arms, and extended halo. In this article, we give an overview of Galactic magnetic field investigation using pulsars. The sensitive SKA1 design baseline (AA4) will increase the number of known pulsars by a factor of around three, and the initial staged delivery array (AA*) will probably double the total number of the current pulsar population. Polarization observations of pulsars with the AA* telescopes will give rotation measures along several thousand lines of sight, enabling detailed exploration of the magnetic structure of both the Galactic disk and the Galactic halo.

Produced by the interaction between the ``pulsar wind'' powered by the rotational energy of a neutron star and its surroundings, the study of pulsar wind nebulae (PWNe) provides vital insight into the physics of neutron star magnetospheres and ultra-relativistic outflows. Spatially-resolved studies of the continuum and polarized radio emission of these sources are vital for understanding the production of $e^\pm$ in the magnetospheres of neutron stars, the acceleration of these particles to $\gtrsim10^{15}~{\rm eV}$ energies, and the propagation of these particles within the PWN as well as the surrounding interstellar medium. The significant improvements in sensitivity, dynamic range, timing capabilities offered by the Square Kilometer Array have the potential to significantly improve our understanding of the origin of some of the highest energy particles produced in the Milky Way.

Binary (and trinary) radio pulsars are natural laboratories in space for understanding gravity in the strong field regime, with many unique and precise tests carried out so far, including the most precise tests of the strong equivalence principle and of the radiative properties of gravity. The Square Kilometre Array (SKA) telescope, with its high sensitivity in the Southern Hemisphere, will vastly improve the timing precision of recycled pulsars, allowing for a deeper search of potential deviations from general relativity (GR) in currently known systems. A Galactic census of pulsars will, in addition, will yield the discovery of dozens of relativistic pulsar systems, including potentially pulsar -- black hole binaries, which can be used to test the cosmic censorship hypothesis and the ``no-hair'' theorem. Aspects of gravitation to be explored include tests of strong equivalence principles, gravitational dipole radiation, extra field components of gravitation, gravitomagnetism, and spacetime symmetries. In this chapter, we describe the kinds of gravity tests possible with binary pulsar and outline the features and abilities that SKA must possess to best contribute to this science.

Matter inside neutron stars is compressed to densities several times greater than nuclear saturation density, while maintaining low temperatures and large asymmetries between neutrons and protons. Neutron stars, therefore, provide a unique laboratory for testing physics in environments that cannot be recreated on Earth. To uncover the highly uncertain nature of cold, ultra-dense matter, discovering and monitoring pulsars is essential, and the SKA will play a crucial role in this endeavour. In this paper, we will present the current state-of-the-art in dense matter physics and dense matter superfluidity, and discuss recent advances in measuring global neutron star properties (masses, moments of inertia, and maximum rotation frequencies) as well as non-global observables (pulsar glitches and free precession). We will specifically highlight how radio observations of isolated neutron stars and those in binaries -- such as those performed with the SKA in the near future -- inform our understanding of ultra-dense physics and address in detail how SKAO's telescopes unprecedented sensitivity, large-scale survey and sub-arraying capabilities will enable novel dense matter constraints. We will also address the potential impact of dark matter and modified gravity models on these constraints and emphasise the role of synergies between the SKA and other facilities, specifically X-ray telescopes and next-generation gravitational wave observatories.

Pulsar timing arrays (PTAs) are ensembles of millisecond pulsars observed for years to decades. The primary goal of PTAs is to study gravitational-wave astronomy at nanohertz frequencies, with secondary goals of undertaking other fundamental tests of physics and astronomy. Recently, compelling evidence has emerged in established PTA experiments for the presence of a gravitational-wave background. To accelerate a confident detection of such a signal and then study gravitational-wave emitting sources, it is necessary to observe a larger number of millisecond pulsars to greater timing precision. The SKAO telescopes, which will be a factor of three to four greater in sensitivity compared to any other southern hemisphere facility, are poised to make such an impact. In this chapter, we motivate an SKAO pulsar timing array (SKAO PTA) experiment. We discuss the classes of gravitational waves present in PTA observations and how an SKAO PTA can detect and study them. We then describe the sources that can produce these signals. We discuss the astrophysical noise sources that must be mitigated to undertake the most sensitive searches. We then describe a realistic PTA experiment implemented with the SKA and place it in context alongside other PTA experiments likely ongoing in the 2030s. We describe the techniques necessary to search for gravitational waves in the SKAO PTA and motivate how very long baseline interferometry can improve the sensitivity of an SKAO PTA. The SKAO PTA will provide a view of the Universe complementary to those of the other large facilities of the 2030s.

A Bayesian analysis of the astrophysical $S$ factor for the $^{12}\mathrm{C}+^{12}\mathrm{C}$ fusion reaction is presented, based on available experimental information at carbon--carbon relative energies $E \gtrsim 2~\mathrm{MeV}$, including direct measurements, indirect Coulomb-renormalized Trojan Horse Method (THM) results, and recent inverse-kinematics data. The Bayesian inference is performed on the quantity $\log_{10}S^{*}(E)$ rather than on $S^{*}(E)$ itself, which naturally accommodates the wide dynamic range of the data and leads to approximately Gaussian uncertainties. The logarithm of the astrophysical factor is parametrized by a quadratic polynomial in energy, and the posterior distribution of the fit coefficients is determined using a weighted Bayesian regression. From this posterior, a global median $S^{*}(E)$ curve is constructed, and the associated covariance matrix is used to define a low/medium/high (LO/MED/HI) band corresponding to a $68\%$ credible interval. Particular emphasis is placed on the extrapolation below $E_{\mathrm{cm}}=2~\mathrm{MeV}$, where the fusion reaction rate is most relevant for stellar carbon burning. At $E_{\mathrm{cm}}=1.5~\mathrm{MeV}$, the posterior distribution yields $S_{\mathrm{global}}^{*}(1.5~\mathrm{MeV})= \left(2.13^{+0.01}_{-0.01}\right)\times10^{16}\,\mathrm{keV\,b}, $ corresponding to a $68\%$ credible interval. The extracted result is consistent with recent inverse-kinematics measurements and with Coulomb-corrected Trojan Horse Method constraints, providing a tightly constrained estimate of the $^{12}\mathrm{C}+^{12}\mathrm{C}$ fusion $S$ factor in the energy region relevant for stellar carbon burning.

The minimum variability timescale (MVT) is a key observable used to probe the central engines of Gamma-Ray Bursts (GRBs) by constraining the emission region size and the outflow Lorentz factor. However, its interpretation is often ambiguous: statistical noise and analysis choices can bias measurements, making it difficult to distinguish genuine source variability from artifacts. Here we perform a comprehensive suite of simulations to establish a quantitative framework for validating Haar-based MVT measurements. We show that in multi--component light curves, the MVT returns the most statistically significant structure in the interval, which is not necessarily the fastest intrinsic timescale, and can therefore converge to intermediate values. Reliability is found to depend jointly on the MVT value and its signal-to-noise ratio ($\mathrm{SNR}_{\mathrm{MVT}}$), with shorter intrinsic timescales requiring proportionally higher $\mathrm{SNR}_{\mathrm{MVT}}$ to be resolved. We use this relation to define an empirical MVT Validation Curve, and provide a practical workflow to classify measurements as robust detections or upper limits. Applying this procedure to a sample of Fermi-GBM bursts shows that several published MVT values are better interpreted as upper limits. These results provide a path toward standardizing MVT analyses and highlight the caution required when inferring physical constraints from a single MVT measurement in complex events.

JWST has uncovered a substantial population of Massive ($M_{\star} \gtrsim 10^{10 }\mathrm{M_{\odot}}$), Quenched Galaxies (MQGs) in the early Universe ($z \gtrsim 2$), whose properties challenge current galaxy formation models. In this paper, we examine this population of MQGs within the new COLIBRE cosmological hydrodynamical simulations. We report number densities and stellar mass functions in broad agreement with the latest observations. The predicted quenching and formation timescales are qualitatively consistent with observational inferences. Leveraging the state-of-the-art physics in COLIBRE, the model predicts that MQGs have dust and $\rm H_{2}$ fractions more than $1$ dex lower than their massive star-forming counterparts; while their sizes and kinematics remain broadly similar. We further explore the processes driving galaxies to become massive and quenched in COLIBRE, identifying active galactic nucleus (AGN) feedback as the primary quenching mechanism. Compared to star-forming galaxies of similar mass, MQGs host more massive black holes (BHs) and exhibit higher star formation efficiencies. These differences arise from their environments, particularly at local ($\rm 0.3\,cMpc$) to intermediate scales ($\rm 1.0\,cMpc$) before quenching, where overdense regions are associated with enhanced gas inflows, higher BH accretion and, hence, feedback power. We find that about $55\%$ of MQGs survive as the main progenitors of $z=0$ galaxies when they are selected at $z=3$, although up to $55\%$ experience rejuvenation episodes. Our results provide robust predictions for MQGs, show that tensions with observations are reduced when an effective observational uncertainty is forward-modelled, and clarify the mechanisms behind their origin.

We exploit a new sample of around 400 bright dusty galaxies from the ALMA CHAMPS Large Program, together with the rich JWST multi-band data products in the COSMOS field, to explore and validate new selection methods for identifying dusty star-forming galaxies (DSFGs). Here, we present an effective empirical selection criterion based on a newly defined parameter: I_star = log(M_star) x log(SFR). Incorporating the F277W-F444W color as a second parameter further improves the purity of the selection. We then apply this method to the COSMOS2025 catalog to search for fainter dusty galaxy candidates below the ALMA CHAMPS detection limit and, through a stacking technique, identify a population of high-redshift (z=6-8) DSFGs with an average flux density of$S_1.2mm = 0.15uJy and a space density of ~6E-6 Mpc^-3. This faint population seems to have been missed by most of the previous submillimeter/millimeter surveys, and ground- and space-based UV-to-NIR surveys. Finally, we discuss the possibility of an evolutionary connection between the z > 10 UV-bright galaxies recently discovered by JWST, the faint dusty z=6-8 galaxies identified here, and the population of z=3-5 massive quiescent galaxies, potentially linked as progenitor-descendant populations based on their abundance, redshifts, and stellar masses.

The common envelope (CE) interaction between an expanding giant star and a compact companion typically leads to a rapid orbital decay, ending in either a merger or the formation of a close binary. However, the existence of post-red giant and post-asymptotic giant branch binaries with separations of 100 to 800 Rsun challenges this standard picture, as these systems appear to have experienced strong interactions without undergoing a classic CE inspiral. In this work, we investigate the effect of high mass ratio, q = M2/M1, on the CE inspiral using three-dimensional hydrodynamical simulations performed with the smoothed particle hydrodynamics code PHANTOM. The primary is a 0.88 Msun, 90 Rsun red giant branch star, while the companion masses span q = 0.68 to 1.5. Higher mass ratios lead to wider post-CE separations, with a maximum of approximately 40 Rsun. The pre-CE mass transfer phase is longer for larger companion masses, and for q greater than or equal to 1 the inspiral becomes significantly more stable, broadly consistent with analytical expectations. This phase is not fully converged with respect to numerical resolution, and higher resolution simulations are expected to further increase its duration and stability. Although higher q systems show enhanced mass loss through the L2 and L3 Lagrange points, we find that circumbinary discs are more likely to form from fallback of bound envelope material. Fallback times are short, of order a few hundred years, and fallback radii lie well outside the binary, between 0.5 and 5 au, where discs are expected to spread efficiently through viscous torques. While high mass ratio systems produce wider post-interaction separations, these remain smaller than those observed. In contrast, fallback-formed discs have properties consistent with observed circumbinary discs.

The coming decades will see gravitational-wave (GW) astronomy expand decisively into the mHz-Hz frequency range, opening access to a population of compact binaries that are currently invisible or only detectable moments before merger. The Lunar Gravitational Wave Antenna (LGWA) concept is designed to probe this gap, enabling continuous observation of compact binaries over months to years prior to coalescence, and detecting sources inaccessible to both space-based mHz detectors and current ground-based >10 Hz facilities. This new GW window fundamentally alters the landscape of time-domain multi-messenger astronomy. Rather than reacting to mergers after the fact, LGWA enables predictive, scheduled electromagnetic (EM) follow-up, transforming how compact-object mergers, their environments, and their astrophysical channels are studied. However, fully exploiting LGWA discoveries requires EM capabilities that do not exist today and are unlikely to be available by the 2030s, particularly for wide-area, rapid, spectroscopically rich follow-up at optical and near-infrared wavelengths. This White Paper identifies the key science cases enabled by LGWA that motivate new ground-based capabilities in the 2040s.

With the installation of next-generation phased array feed (PAF) receivers on radio telescopes, there is an urgent need to develop effective and computationally efficient radio frequency interference (RFI) mitigation methods for large-scale surveys. Here we present a new RFI mitigation package, called mRAID (multi-beam RAdio frequency Interference Detector), which uses the eigenvalue decomposition algorithm to identify RFI in cross-correlation matrix (CCM) of data recorded by multiple beams. When applied to high time-resolution pulsar search data from the Five-hundred-meter Aperture Spherical Radio Telescope (FAST), mRAID demonstrates excellent performance in identifying RFI over short timescales, thereby enhancing the efficiency of pulsar and fast radio burst (FRB) searches. Since the computation of the CCM and the eigenvalue decomposition for each time sub-integration and frequency channel are independent, the process is fully parallelisable. As a result, mRAID offers a significant computational advantage over commonly used RFI detection methods.

The gravitational-wave signal from the neutron star-black hole (NSBH) merger GW200105 is consistent with this binary having significant orbital eccentricity close to merger. This raises the question of how an eccentric NSBH might form. Compact object binaries that form via isolated binary star evolution should radiate away any orbital eccentricity long before their gravitational-wave signal enters the sensitive frequency range of the LIGO-Virgo-KAGRA detector network. Meanwhile, dynamical environments -- which can be conducive to mergers on eccentric orbits -- produce very few NSBHs. We estimate the minimum measurable eccentricity of NSBHs at 10 Hz orbit-averaged gravitational-wave frequency, $e_{\mathrm{min},10}$, finding that for GW200105, GW200115, and GW230529-like systems, $e_{\mathrm{min},10}$ is O(0.01). For a GW190814-like unequal-mass binary with significant higher-order mode content, $e_{\mathrm{min},10}=0.003$; this is an order of magnitude lower than when higher-order modes are not present. For dominant-mode signals from eccentric binaries with $m_2=1.5$ M$_\odot$ and a range of total masses from $3\,{\rm M_\odot} \leq M \leq50\,\rm M_\odot$, we find $0.008\leq e_{\mathrm{min},10}\leq0.022$. The relationship between $M$ and $e_{\mathrm{min},10}$ is linear when the binaries are non-spinning. When the binaries are maximally spin-precessing, $e_{\mathrm{min},10}$ decreases as mass ratio becomes more unequal. We estimate the sensitivity of a quasi-circular templated search to a population of NSBHs from field triples, finding that we recover only 27% of our simulated population. Finally, we show that if ~1/3 of present NSBH detections are measurably eccentric, then all of them are consistent with an isolated field triple origin.

Understanding the quenching of star formation in central galaxies remains a core challenge in galaxy evolution. Two decades ago, the concept of halo quenching was introduced as a dominant mechanism, positing that massive central galaxy quenching is governed by the thermodynamics of gas predominantly influenced by dark matter halos. However, a vastly increasing body of observational evidence consistently indicates that quenching correlates predominantly with central properties like velocity dispersion, bulge mass, and black hole mass. When these properties are controlled, halo mass appears to show weak influence, supporting AGN feedback as the primary mechanism. A recurring key issue, however, is that these studies rely on halo masses derived via abundance matching (AM). Direct observational measurements from weak lensing, satellite kinematics, and galactic dynamics reveal that AM systematically overestimates halo masses of star-forming centrals while underestimating those of passive ones. To accurately assess the true role of halo mass. we re-estimated halo masses for SDSS groups; the resulting halo mass function and stellar-to-halo mass relations (SHMRs) for both populations match theoretical predictions and weak lensing measurements. Using these improved masses, we find direct observational evidence that halo mass is the dominant factor in quenching central galaxies, with a clear threshold at $M_{h}\sim 10^{12.1}M_\odot; $. By applying a simple correction to AM data using weak lensing-derived SHMRs, we demonstrate that previous claims regarding the dominance of central properties stem primarily from systematic biases in AM halo masses. Our results suggest that the significance of AGN feedback is primarily manifested in halos above this mass threshold, in galaxies already primed for quenching. In other words, AGN feedback appears to become effective in halos above this mass threshold.

Solar energetic protons (SEPs) in different energy channels from 10 to above 100 MeV are analyzed and their relationship to solar and geomagnetic activity is investigated. We performed temporal association analysis between the SEPs, solar flares (SFs), coronal mass ejections (CMEs) and geomagnetic storms (GSs) that occurred during solar cycles 23 and 24. The energy dependencies between the SEPs and the strength of the space weather activity are evaluated and presented.

We present a spectropolarimetric study of the nearby M4.5V exoplanet host star YZ Cet, based on near-infrared observations obtained with the SpectroPolarimètre InfraRouge (SPIRou) at the Canada--France--Hawaii Telescope. We detect striking changes in the large-scale magnetic field strength and geometry over the course of just a few stellar rotations, a level of short-term global magnetic field evolution rarely reported in M dwarfs. We modeled the temporal variation of the longitudinal magnetic field using a Gaussian regression process, which allowed us to robustly determine the stellar rotation period and quantify the evolution timescale of the magnetic field. Independent Zeeman Doppler Imaging reconstructions of the two epochs confirm a significant reconfiguration of the star's global magnetic strength and topology. The detection of a weaker, complex, axisymmetric magnetic field (mean $|B| \sim 201$~G), which changes into a stronger, non-axisymmetric, dipole-dominated field (mean $|B| \sim 276$~G) over a few rotation cycles, is in contrast to results from similar fully convective M-dwarf stars. YZ Cet is known to exhibit polarized radio bursts potentially driven by auroral radio emission from star--planet interaction (SPI). By combining our magnetic maps with recent radio observations, we refine the constraints on the magnetic field strength of the innermost planet, YZ Cet b. These results underscore the importance of monitoring stellar magnetic variability to interpret multi-wavelength SPI signatures and to characterize the magnetospheres of potentially habitable exoplanets.

Thousands of X-ray sources have been detected in the Galactic center (GC), most believed to be cataclysmic variables (CVs). As a potential probe of the old stellar population, in particular CVs, the existence and detectability of novae in the GC remain elusive, due to the prohibitive extinction toward the GC and their relatively low occurrence rate. Nova remnants evolving in the characteristic hot ($T\sim{10^{6}~\rm K}$) and dense ($n_e\sim{10~\rm cm^{-3}}$) interstellar medium in the GC may shed light on recent novae and provide useful insight on the GC ecosystem. In this work, we perform hydrodynamical simulations of putative nova remnants in the GC environment and calculate their time-dependent multiwavelength emission to estimate the detectability. Among 79 models sampling the nova parameter space (primarily ejecta mass and velocity), 6, 44, and 51 modelled nova remnants are detectable at their X-ray, radio, and Paschen-$\alpha$ maximum, respectively, for existing {\it Chandra}, VLA, and HST observations of the GC. The predicted peak luminosities are $\sim10^{32}~\rm erg~s^{-1}$, $\sim10^{31}~\rm erg~s^{-1}$, and $\sim10^{36}~\rm erg~s^{-1}$ in these three bands and the detectable window ranges from weeks to notably hundred years. By specifying a CV population of the nuclear star cluster, we estimate the probability of detecting at least one remnant to be 20\%, 8\%, and 18\% in X-rays, radio, and Pa$\alpha$. The nova remnant would be best resolved in the X-ray band. Our study highlights the potential for detecting nova remnants through further observations, leveraging JWST and the potentially forthcoming AXIS and SKA.

This paper proposes four new methods to decontaminate spectra of stars and galaxies resulting from slitless spectroscopy used in many space missions such as Euclid. These methods are based on two distinct approaches and simultaneously take into account multiple dispersion directions of light. The first approach, called the local instantaneous approach, is based on an approximate linear instantaneous model. The second approach, called the local convolutive approach, is based on a more realistic convolutive model that allows simultaneous decontamination and deconvolution of spectra. For each approach, a mixing model was developed that links the observed data to the source spectra. This was done either in the spatial domain for the local instantaneous approach or in the Fourier domain for the local convolutive approach. Four methods were then developed to decontaminate these spectra from the mixtures, exploiting the direct images provided by photometers. Test results obtained using realistic, noisy, Euclid-like data confirmed the effectiveness of the proposed methods.

The ELT will provide groundbreaking science across a wide range of areas, including small habitable-zone exoplanets; however, true Earth analogs in the habitable zones of Sun-like stars are generally beyond the reach even of the ELT, due to the extreme contrast ratio and small angular separation between the planet and star. Here, we note that the combination of ELT and a space-based starshade would provide the contrast needed to observe potentially tens of Earth analogs, as well as other planets. This would yield the scientific basis needed for addressing central scientific questions regarding the frequency and distribution of habitability and life in the Universe. The huge aperture of ELT, combined with a contrast otherwise only reachable in space, opens up scientific avenues that are unmatched by any other existing or foreseen facility. ESO could conceivably collaborate with ESA (and others) to facilitate a starshade mission suitable for synergy with the ELT, as well as to prepare the ELT instrumentation in order to maximize its potential for synergy with a starshade.

The aim of the present report is to outline the observational methods used to determine the solar origin - in terms of flares and coronal mass ejections (CMEs) - of the in situ observed solar energetic protons. Several widely used guidelines are given and different sources of uncertainties are summarized and discussed. In the present study, a new quality factor is proposed as a certainty check on the so-identified flare-CME pairs. In addition, the correlations between the proton peak intensity and the properties of their solar origin are evaluated as a function of the quality factor.

The semi-forbidden CIII] $\lambda\lambda$1907,1909 doublet is a key tracer of high-ionization emission in the early universe. We present a study of CIII] emission in galaxies at z=5-7, using publicly available JWST/NIRSpec prism data from programs including CEERS, JADES, RUBIES and CAPERS. We built a sample of 61 CIII]-emitting galaxies, and we classified them as star-forming or active galactic nuclei (AGN) host galaxies using (1) rest-frame UV and optical emission-line diagnostic diagrams, and (2) the presence/absence of broad Balmer emission lines. The UV diagnostics are based on the combination of the rest-frame equivalent width (EW) of CIII] versus CIII]/HeII $\lambda$1640, and the EW of CIV versus CIV/HeII $\lambda$1640. For optical diagnostics, we employ the OHNO diagram, which combines [OIII] $\lambda$5007, H$\beta$, [NeIII] $\lambda$3869, and [OII] $\lambda\lambda$3727,3729- and we find it has a low efficiency on separating AGN from SFG. We find that half of the sources in our sample (29 out of 61 galaxies) exhibit at least one secure indication of AGN activity while 13 are potential AGNs based on the CIII] diagnostic. Physical properties, including stellar mass and star formation rate, are derived through spectral energy distribution modeling with Bagpipes. Our analysis reveals that JWST is uncovering a population of strong CIII] emitters at high redshift (5<z<7) with a median EW of 22.8 A. This EW is higher than that of a control sample of CIII] emitters at redshift 3<z<4 with a median EW of 4.7 A. We find that for the same range of Muv, the CIII] EW increases by $\sim$0.67 dex from 3<z<4 to 5<z<7, indicating strong redshift evolution in the line's strength. Finally, we identify five sources in our sample as Little Red Dots (LRDs); while four of these have already been identified as LRD in the literature, one is presented here for the first time.

Simulating one-dimensional stellar evolution models with MESA, we show that removing the outer inflated envelope of a mass-accreting evolved stripped-envelope star, like a Wolf-Rayet (WR) star, substantially moderates the stellar expansion during accretion at high-mass accretion rates. We study the accretion onto a star via an accretion disk, which launches jets that remove the high-entropy outer layers of the inflated envelope. This is the `jetted mass removal accretion scenario.' By manually removing the entire hydrogen-rich envelope from a red supergiant, we build a hydrogen-deficient WR stellar model with a mass of 6.03Mo and a radius of 0.67Ro. We then accrete mass onto it at a high rate. We mimic the real process of simultaneous mass addition near the equatorial plane and jet-induced mass removal from the outer envelope by dividing the accretion period into hundreds of pulses: in the first half of each pulse, we add mass, and in the second, we remove a fraction of this mass. The removal of tens of percent from the added mass decreases the stellar expansion by a factor of about 2-5. Our results show that WR stars can maintain a deep potential well and not expand much while accreting mass at high rates. This allows the formation of an accretion disk and the liberation of large amounts of gravitational energy. Our results strengthen models of intermediate-luminosity optical transients, such as luminous red novae, in which a non-degenerate star accretes at high rates and launches jets that power the transient event.

Habitable Worlds Observatory (HWO) will be NASA's flagship space telescope of the 2040s, designed to search for life on other planets and to transform broad areas of astrophysics. NASA are seeking international partners, and the UK is well-placed to lead the design and construction of its imaging camera - which is likely to produce the mission's most visible public impact. Early participation in the mission would return investment to UK industry, and bring generational leadership for the UK in space science, space technology, and astrophysics.

Characterizing temperate exoplanet atmospheres remains challenging due to their small size and low temperatures. Recent JWST observations provide valuable data, but their interpretation has led to diverging conclusions. Complementary approaches combining laboratory experiments and photochemical modeling are essential for constraining atmospheric chemistry and interpreting observations. We aim to identify chemical pathways governing the formation and evolution of neutral species and to assess their sensitivity to key parameters such as C/O ratio and metallicity. Our approach combines experimental and numerical simulations on H2-rich gas mixtures representative of sub-Neptune atmospheres, spanning a wide range of CH4, CO, and CO2 mixing ratios. A cold plasma reactor simulates out-of-equilibrium upper-atmospheric chemistry. A 0D photochemical model reproduces reactor conditions, guiding interpretation of key pathways and abundance trends. We observe the formation of both reduced and oxidized organic compounds. In CH4-rich mixtures, hydrocarbons form efficiently through methane chemistry, correlating with CH4 concentration and agreeing with models. In more oxidizing environments, particularly CO2-rich mixtures, hydrocarbon formation is inhibited by complex reaction networks and oxidative losses. Oxygen incorporation enhances chemical diversity and promotes formation of oxidized organic compounds of prebiotic interest (H2CO, CH3OH, CH3CHO), especially in atmospheres containing both CH4 and CO2. Atmospheres containing CH4 and CO, which balance carbon and oxygen supply without excessive oxidative destruction, favor efficient production of hydrocarbons and oxidized compounds. Out-of-equilibrium chemistry plays a key role in the diversification and organic complexification of temperate exoplanet atmospheres.

We investigate the connection between strong, blended Ly$\alpha$ absorption systems (SBLAs) and $\approx1000$ Ly$\alpha$ emitting galaxies (LAEs) at $z\gtrsim3$ in 28 quasar fields from the MUSE Analysis of Gas around Galaxies (MAGG) survey. Selecting SBLAs as spectral regions with transmitted flux $-0.05<F<0.25$ over $\approx138\text{ km s}^{-1}$ bins, we find a strong correlation with LAEs within a projected distance of $R\le300\rm\,kpc$ and line-of-sight velocity separation of $|\Delta v|\le300 \text{ km s}^{-1}$. The association rate increases significantly with decreasing flux, a trend that persists also at smaller separations ($R<100$ kpc). A two-dimensional cross-correlation analysis confirms significant clustering of LAEs around SBLAs, while no such clustering is seen for spectral regions with $F>0.25$. The correlation appears to also depend on the width of the spectral window used to identify SBLAs, with a larger window yielding a stronger signal. Our analysis confirms that SBLAs serve as probes of the CGM at the interface between the Ly$\alpha$ forest and the optically-thick Lyman limit systems. The significant dependence of the LAE-SBLA cross-correlation on the spectral binning used to select these absorbers motivates future tests of the current SBLA framework as a tracer of halos.

With the increasing sensitivity of modern radio interferometers, it has become important to image objects larger than the field of view while optimising sensitivity and image fidelity. We present a coherent visibility plane direction-dependent imaging, calibration and mosaicing framework. Our simulations and application to real MeerKAT data show that this joint deconvolution and primary beam correction approach, coupled with direction-dependent calibration, allows for deeper mosaics with greater fidelity and increased accuracy of recovered flux densities and spectral indices, especially beyond the half-power beam width. Our best-case mosaic produces precise flux values within a 6% uncertainty and spectral indices within 20\% throughout the imaged area, and is fully complete out to twice the radii and half the flux density than the image plane equivalent. The application to archival wideband MeerKAT 1283 MHz data produces the deepest high-resolution image of the Shapley Supercluster Core, with a sensitivity of 3.6 $\mu$Jy/beam within the primary beam at a 7$^{\prime\prime}$ resolution, constituting a $\sim$ 50% increase in dynamic range over the image plane counterpart, and a fluxscale that is consistent within 10% across the entire field of view. The compute time for the direction-dependent visibility plane mosaic was comparable to the sum of the times needed to perform direction-dependent calibration on the individual pointings. Our results suggest that visibility plane mosaicing with its capability for deeper deconvolution could improve the efficiency of deep and wide surveys, particularly for on-the-fly mapping and studies of low surface brightness sources, and could form the basis of future calibration pipelines for SKA-scale instruments.

The post-merger phase of binary neutron star (BNS) mergers encodes valuable information about the equation of state (EOS) of supranuclear matter. Extracting this information from the analysis of the post-merger waveforms remains challenging due to the high-frequency limitations of current detectors. Future third-generation observatories, such as the Einstein Telescope (ET) and NEMO, will have the sensitivity required to resolve post-merger signals with high fidelity. In this work, we apply CLAWDIA, our recently developed sparse dictionary learning (SDL) framework, to classify different EOS models using only the post-merger gravitational-wave emission of simulated BNS mergers available in the CoRe database. Our dataset comprises five EOS models representative of a broad range of neutron star properties. The SDL framework is optimised under realistic detection conditions by injecting signals into simulated noise matching the sensitivity curves of ET and NEMO. Our results show that classification is primarily driven by the dominant post-merger frequency, $f_2$, which encodes EOS-dependent information. At a modest signal-to-noise ratio of 5, our method achieves $F_1$ scores of $0.76$ for ET and $0.70$ for NEMO, with performance improving for higher signal-to-noise ratios. The reliability and generalisation capabilities of the model are assessed with additional tests, including the classification of an EOS not included in the training dataset and the analysis of detector-specific biases.

The Legacy Survey of Space and Time (LSST) that will be carried out by the NSF-DOE Vera C. Rubin Observatory promises to be the defining survey of the next decade, supplying unprecedented access to the night sky to static science- and time-domain science-focused researchers alike. Maximizing the output of the broad remit of Rubin Observatory science requires a non-trivial survey strategy. For time-domain science, the most promising strategy designed so far is a rolling survey strategy, whereby a subset of the full LSST survey area is observed at higher rate compared with the nominal rate dictated by weather conditions and the observatory's technical constraints. This strategy is now the baseline approach for the LSST as a whole. Focusing on static science (galaxy clustering and weak lensing), we study how these time-domain-optimized rolling strategies affect the depth uniformity at intermediate years of the survey. We characterize the amount of survey area at high risk of being lost in static-science analyses of a baseline rolling LSST dataset due to an insufficient combination of survey contiguity and uniformity. At intermediate data releases, nearly half of the survey could be lost for static science, decreasing the Dark Energy figure of merit by approximately 40\%. We describe additional metrics focused on key analysis tasks, such as photometric redshifts and galaxy clustering. We propose a new strategy that returns the survey to uniformity at key release years, enabling use of the full survey area and restoring our metrics to the values they would have in a non-rolling cadence without loss of time domain data relative to a rolling survey with the same number of rolling cycles. This work has informed the third round of optimization of the survey strategy, and the new uniform rolling strategies have been incorporated into the baseline strategy.

We explore the effects of massive neutrinos on the cosmic web using the FLAMINGO simulations. We classify the cosmic web into voids, sheets, filaments, and clusters, and find that massive neutrinos affect the environment by decreasing the volume occupied by clusters and voids. We find that increasing the neutrino mass shifts the volume-weighted density distribution towards higher densities and leads to a more narrow density distribution, which we interpret as neutrinos delaying structure formation. We construct the minimum spanning tree (MST) graph from the subhaloes, adopting a number density chosen to match that expected for DESI-like observations. We show that most MST edges lie in filaments, approximately 70% throughout different simulations, which we link to its sensitivity to neutrino mass. We also link the MST's edge length signal at different scales to different cosmic web environments, with clusters dominating the signal at small scales, voids at longer scales, and filaments at intermediate scales. The strong correlation between MST edges and cosmic web environments reinforces the MST's potential to be used as a classifier for large-scale structure in galaxy surveys. We compare the effects of baryonic physics and massive neutrinos and find that each produces distinct signatures in MST edge lengths. This analysis is performed in 3D space, using the true positions of subhaloes and not accounting for redshift space distortions. Nevertheless, these results emphasise the MST's capability to go beyond two-point statistics, motivating future applications to real observational data.

Winds are an important ingredient in the evolution of X-ray binary (XRB) systems, particularly those at high accretion rates such as ultra-luminous X-ray sources (ULXs), because they may regulate the accretion of matter onto the compact object. We aim at understanding the properties of ULX winds and their link with the source spectral and temporal behavior. We performed high-resolution X-ray spectroscopy of the variable source NGC 55 ULX-1 to resolve emission and absorption lines as observed with XMM-Newton at different epochs. Optically-thin plasma models are used to characterise the wind. We confirmed and thoroughly strengthened previous evidence of outflows in NGC 55 ULX-1. The presence of radiative recombination signatures and the ratios between the fluxes of the emission lines favours photoionisation balance and low-to-moderate densities, which confirm that the lines originate from classical XRB disc winds. An in-depth parameter space exploration shows line emission from a slowly moving, cool, and variable plasma perhaps associated with a thermal wind. Mildly-relativistic Doppler shifts (about -0.15c) associated with the absorption lines confirm, at higher confidence, the presence of powerful, radiatively-driven, winds. The comparison between results obtained at different epochs revealed that the wind responds to the variability of the underlying continuum and these variations may be used to understand the actual accretion regime and the nature of the source.

The observations of B1237+25 at a frequency of 111 MHz were analyzed. For the first time in the normal radiation mode a new component in the central region in the average profile was detected. This component is manifested in all modes of pulsar emission: quiet-normal (QN), flare-normal (FN) and in the abnormal mode (AB). The subpulse drift is observed in the QN mode only in the first and last components of the average profile. The normal mode is interrupted by nullings and transitions into the abnormal AB mode. In the AB mode, the structure at the edge of the outer cone is destroyed, the distance between the inner and outer cones is almost doubled, and the distance between the inner cone and the central region is this http URL of our data has shown that the components of the outer and inner cones of the average profile are formed by an ordinary mode of radio emission (O-mode) and form a single cone radiation of the pulsar. The central components of the average profile (wide and narrow) are formed by an extraordinary mode (X-mode). Estimates of the height of the radiation output from the central region (X-mode) and the cone radiation (O-mode) are obtained: 80~km and 370~km, respectively. A microstructure with a time scale of $\tau_\mu\le0.5$~$\mu s$ has been detected. This time scale corresponds well to the time of the development of a spark discharge in the polar cap. For this value $\tau_\mu$, the height of the vacuum gap should be $h_p\le750$~cm. Based on the steepness of the individual pulse's trailing edge at the longitude of the first component, a limit was obtained on the value of the $\gamma$ factor of the relativistic secondary plasma: $\gamma\ge$260. The dependence of the distance between the components of the outer and inner cone of radiation on the frequency is the same and corresponds to a power law with an exponent of -0.16.

Compact Groups (CGs) of galaxies are dense systems where projected separations are comparable to their optical diameters. A subset - non-isolated CGs - are embedded within major structures. Using multi-band S-PLUS data, we analyse galaxies in 122 non-isolated CGs within more massive systems such as larger groups and clusters. We compare them to galaxies in the host structures, hereafter surrounding group galaxies. Structural parameters were obtained with MorphoPLUS, a pipeline for multi-wavelength Sérsic profile fitting. Dividing galaxies into early (ETG), transition, or late types (LTG), we find: (1) Non-isolated CGs host higher quenched fractions and more ETGs, especially for stellar masses $\log(M/M_\odot) > 10.2$, than surrounding groups. (2) Sérsic indices increase with wavelength for all morphological types in both environments, whereas effective radii show a stronger morphology-dependent behaviour - ETGs become more compact towards redder bands, while LTGs exhibit flatter $Re(\lambda)$ trends. Environmental differences remain weak, with only a modest enhancement of the gradients for ETGs in non-isolated CGs. (3) Transition galaxies in CGs show a concentrated $R_e$-$n$ distribution and faint-end bimodality, consistent with ongoing morphological transformation absent in surrounding groups. (4) Phase-space analysis indicates that some CGs in clusters are projection artefacts, while others are genuine dense systems at various infall stages, from recent arrivals to ancient remnants. These results show that galaxies in non-isolated CGs follow distinct evolutionary paths compared to their surrounding groups galaxies, suggesting that the compact configuration plays a unique role beyond the influence of the larger-scale environment.

Based on the finding of unusual chemical abundance ratios of N-emitters, which resemble those of globular cluster (GC) stars, their compactness, high ISM densities and other properties, it has been suggested that N-emitters could indicate the formation sites of globulars. A recent statistical study of the N-emitter population has quantified the frequency $f_N$ of these rare objects and their redshift evolution (Morel et al. 2025). Using these results we here test if N-emitters trace the formation of GCs and use the observed cosmic star-formation rate density evolution to predict the cosmological evolution of the GC population with time, their age distribution, and the total present-day stellar mass density formed in globulars. The predicted age distribution of GCs strongly resembles the typical asymmetric observed distributions in the Galaxy and ellipiticals, with a peak at $\sim 11.5-12$ Gyr and a longer tail extending to younger ages. We derive a total stellar mass density formed in N-emitters down to redshift zero of $(2-7) \times 10^5$ M$_{\odot}$ Mpc$^{-3}$, which matches within a factor $\sim 2$ the observed fraction of stellar mass found in the GC population at $z=0$. These results provide additional indirect arguments supporting the hypothesis that N-emitters could represent sign-posts of a short phase of GC formation.

In this work, we investigate the dynamical survival of short-period inner planets during the high-eccentricity tidal migration of companion exterior giant planets. Using a combination of analytic arguments and N-body simulations including equilibrium tides and general relativistic precession, we find the boundary in parameter space where an inner companion can remain dynamically stable. We find that survival requires a periastron separation exceeding roughly 14 mutual Hill radii at closest approach. Below this threshold, secular eccentricity exchange, orbit crossing, and/or tidal evolution can lead to the destruction of the inner planet. We apply our methodology to the current exoplanet sample and find that none of the known systems containing a short-period giant and an inner companion could have assembled via high-eccentricity tidal migration. However, warm Jupiters with larger periastron distances ($q_{\mathrm{out}} \sim 0.05-0.08$ AU, corresponding to final observed semi-major axis values $a_{\mathrm{out}} \sim 0.10-0.16$ AU) can allow the survival of short-period inner planets while potentially also circularizing on $\lesssim 1$ Gyr timescales. Our results provide a framework for distinguishing disk migration from tidal migration in observed multi-planet systems containing close-in gas giants.

We investigate a perturbation-level modification of symmetric teleparallel gravity of the form $f(Q)=F(Q)+M\sqrt{Q}$ and assess its ability to ease the $\sigma_8$ tension. The square-root term leaves the background expansion unchanged while modifying the effective gravitational coupling, providing a pure decoupling between background cosmology and structure-growth evolution. Using the latest redshift-space distortion data, including DESI DR2 Full-Shape measurements, we constrain $M$ and $\sigma_8$ across three representative backgrounds: $\Lambda$CDM, an $H_0$-tension-reducing model, and a DESI-motivated dynamical dark energy scenario. In all cases, the square-root correction suppresses growth and can reconcile $\sigma_8$ with Planck at the $1\sigma$ level, with the strongest improvement occurring in the $H_0$-tension-oriented background. A residual degeneracy between $M$ and $\sigma_8$ remains, indicating that future multi-probe analyses combining lensing and full-shape clustering will be required to determine whether the $\sqrt{Q}$ term represents a genuine signal of modified gravity.

We built a simple toy model of a core-collapse supernova (CCSN) ejecta composed of two shells, an outer low-mass spherical shell and an inner elongated massive shell, and show that it can reproduce the evolution of the photospheric radius of SN 2024ggi, Rph(t). During the first week, the larger spherical shell, the S-shell, forms the photosphere. As the shell expands and becomes increasingly transparent, the photosphere moves inward along the mass coordinate, although it grows in size. When the photosphere reaches the long axis of the elongated inner shell, the E-shell begins to contribute to the photosphere, ultimately comprising the entire photosphere. The simple toy model explains the transition of Rph(t) from being concave (decreasing slope) to convex (increasing slope). A single-shell model predicts only concave behavior. The structure of a spherical shell with an inner elongated shell is motivated by the morphologies of several CCSN remnants whose structures have been attributed to multiple pairs of jets in the framework of the jittering jets explosion mechanism (JJEM). The deduced multiple-shell ejecta of SN 2024ggi in this study, and of SN 2023ixf in an earlier study, as well as studies of the polarization of SN 2024ggi, are better compatible with the JJEM than with the neutrino-driven mechanism. Our study supports the growing evidence that the JJEM is the primary explosion mechanism of CCSNe.

The search for sources of high-energy astrophysical neutrinos can be significantly advanced through a multi-messenger approach, which seeks to detect the gamma rays that accompany neutrinos as they are produced at their sources. Multi-messenger observations have so far provided the first evidence for a neutrino source, illustrated by the joint detection of the flaring blazar TXS 0506+056 in highenergy (HE, E > 1 GeV) and very-high-energy (VHE, E > 100 GeV) gamma rays in coincidence with the high-energy neutrino IceCube-170922A, identified by IceCube. Imaging atmospheric Cherenkov telescopes (IACTs), namely FACT, H.E.S.S., MAGIC, and VERITAS, continue to conduct extensive neutrino target-of-opportunity follow-up programs. These programs have two components: followup observations of single astrophysical neutrino candidate events (such as IceCube-170922A), and observation of known gamma-ray sources after the identification of a cluster of neutrino events by IceCube. Here we present a comprehensive analysis of follow-up observations of high-energy neutrino events observed by the four IACTs between September 2017 (after the IceCube-170922A event) and January 2021. Our study found no associations between gamma-ray sources and the observed neutrino events. We provide a detailed overview of each neutrino event and its potential counterparts. Furthermore, a joint analysis of all IACT data is included, yielding combined upper limits on the VHE gamma-ray flux.

Tidal disruption events (TDEs) occur when a star crosses the tidal radius of a black hole (BH) and is ripped apart, providing a novel and powerful way to probe dormant BHs over a wide mass range. In this study, we present our late-time observations and comprehensive multi-wavelength analyses of an extraordinary TDE at the center of a dwarf galaxy, which exhibited successive flares in the optical, X-ray, and radio bands. Notably, we discovered an unexpected high-state X-ray plateau phase following the peak until the present time. Along with its reported prolonged rise lasting at least 550 days, these unique characteristics are consistent with the scenario of a TDE caused by an intermediate-mass black hole (IMBH) with a mass of approximately $(1-6) \times 10^5$ solar masses. Furthermore, scaling relations derived from the host-galaxy properties indicated a similar BH mass in concert. This discovery highlights the invaluable role of TDEs in the search for elusive IMBHs.

Future direct-imaging missions, such as the Large Interferometer for Exoplanets (LIFE), aim to observe thermal emission from potentially habitable planets to characterize their surface environments and search for signs of life. Previous studies of directly imaged Earth-like planets have mainly examined the signatures of atmospheric composition, often using one-dimensional models, while the effect of horizontal temperature gradients has received limited attention. Because a pronounced horizontal temperature gradient may signal the absence of a global ocean, we investigate its detectability through thermal-emission direct imaging. Adopting Teegarden's Star b (zero-albedo equilibrium temperature $\sim 280$~K) as a benchmark, we compute three-dimensional atmospheric structures with and without a global ocean using the ROCKE-3D general circulation model and simulate geometry-dependent thermal emission spectra. We show that the temperature gradients that disfavor a global-ocean scenario manifest in both orbital phase variation and spectral shape of the snapshot spectra. The phase variation is more readily detectable: one-day integrations with LIFE at two orbital phases would reveal flux variations in no-ocean cases with 1-10~bar atmospheres, depending on background atmospheric composition. Shapshot spectra provide complementary diagnostics of global temperature contrast, the running brightness temperature of the continuum and detailed absorption band shapes, but require integration a few times longer. These three-dimensional effects, if neglected, can bias interpretations based on one-dimensional models. We also assess their detectability for other nearby exoplanets. Our results highlight the need to incorporate three-dimensional atmospheric structures when characterizing rocky exoplanets, both to constrain surface conditions and to avoid misinterpretation of spectral data.

We investigated the James Webb Space Telescope photometric color classification of mid-infrared (MIR) selected galaxies at high redshifts, toward cosmic noon. The aim of the present work is to obtain a z-dependent mid-infrared (MIR) photometric galaxy classification tool based on broad spectral emission and absorption lines using the JWST Mid-Infrared Instrument (MIRI) and its broadband filters. We used the largest Spitzer MIR spectral database to obtain synthetic photometry in the JWST/MIRI filters. We formed MIRI filter combinations to trace the strong polycyclic aromatic hydrocarbon (PAH) emission features and the 9.7 micron silicate feature in seven redshift windows from z = 0.25-2.10. Results. We present z-dependent MIRI color-color plots that separate active galactic nuclei (AGN), star-forming galaxies (SFGs), and silicate absorption-dominated galaxies up to z$\sim$2. We applied the photometric MIR colors to the largest ($\sim$34 arcmin square) MIRI survey called the Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES), to identify AGN, SFGs, and Si-absorption dominated galaxies out to substantial redshifts. Our JWST/MIRI SFGs sample includes galaxies with total IR luminosities of $10^{9.2} \sim 10^{11.9} L_{\odot}$ at 0.9 $\leq$z < 1.57. The majority of them are consistent with the z$\sim$1 main sequence. We also identified the first examples of z$\sim$1 galaxies with deep silicate absorption.

We develop an effective description of the local cosmic environment, namely, for redshift $z \lesssim 0.1$, to quantify the bias induced by local structure on cosmological observables. Our approach models the metric of the nearby Universe as a superposition of multi-structured $\Lambda$-Szekeres patches, calibrated against the HAMLET peculiar velocity and density field reconstructions of Cosmicflows-4. From this framework we compute the fully inhomogeneous and anisotropic quasilocal expansion field predicted by our model, and use it to assess the impact of local structure on estimates of $H_0$. For this purpose we analyse low-redshift Type Ia supernovae from the Pantheon+ catalogue. We find that accounting for the local structure increases the Hubble tension, yielding a shift in the best-fit value of the Hubble constant of order $\Delta H_0 \approx 0.5\ \mathrm{km\,s^{-1}Mpc^{-1}}$.

We present several closed-form expressions of useful mass distributions. These include the potentials and accelerations of circular rings and arcs, the potentials of uniform density rings and arcs at arbitrary eccentricities, and the potentials and accelerations of rings and arcs when the mass is time-averaged over a Kepler orbit. We show that these expressions can be expressed, often simply, in terms of elliptic functions of complex arguments. We show that in a few limiting cases, the expressions are entirely real. We expect that these expressions will allow for more rapid modeling in many areas of celestial mechanics.

Future telescopes such as the Large Interferometer For Exoplanets (LIFE) will enable mid-infrared characterisation of the atmospheres of nearby rocky exoplanets. Whilst 4D spatial and temporal variations of Earth as an exoplanet are below spectroscopic detection limits, such variability is planet-specific. We investigate LIFE's ability to detect 4D variability in the atmospheres of tidally locked exoplanets. We create daily synthetic LIFE observations of Proxima Centauri b in a 1:1 and an eccentric 3:2 spin-orbit resonance (SOR), using LIFEsim on spectra from daily 3D climate-chemistry model output of an aquaplanet with Earth-like composition. Hemispheric distributions of temperature, clouds, and chemical species determine spectral signatures and variability with orbital phase angle. Such variability dictates the extent to which parameters can be reliably inferred from snapshot spectra at arbitrary viewing geometries. In the 1:1 SOR, MIR spectra vary significantly with viewing geometry and indirectly probe atmospheric circulation. Nightside temperature inversions generate O3, CO2, and H2O emission features, though these lie below LIFE's detection threshold, and instead O3 features disappear at certain phase angles. In contrast, the 3:2 SOR yields a more homogeneous atmosphere with weaker phase variability but enhanced bolometric flux due to eccentric heating. Phase-resolved LIFE observations confidently distinguish between the SORs and capture seasonal O3 variability for golden targets like Proxima Centauri b. In case of abiotic O2/O3 build-up, the O3 variability presents a potential false positive scenario. Hence, LIFE can disentangle different spin-orbit states and resolve 4D atmospheric variability, enabling daily characterisation of the 4D physical and chemical state of nearby terrestrial worlds. Importantly, this characterisation requires phase-resolved rather than snapshot spectra.

We report on the very first radio detection associated with the peculiar hourglass-morphology X-rays surrounding 3C 294 at z=1.8. Using International Low Frequency Array (LOFAR) data at 144 MHz and Chandra data at 0.3-6 keV, we find that the co-spatial diffuse radio and X-ray emission is well described by synchrotron and inverse-Compton processes by the same electron population. Through modelling of this rare low-energy plasma, we find that the most defining property of the electrons up-scattering CMB photons at this redshift is very low electron Lorentz factors ($\gamma_{\text{max}}\ll 10^{4}$ and $\gamma_{\text{break}}\lesssim 10^{3}$) in the lobe: deep low frequency (<150 MHz) observations are critical to the detection of radio lobes at high redshift. The physical conditions imply a total energy in the diffuse emission significantly greater than that implied by the temperature of the protocluster gas: 3C 294 is one of the most powerful known radio-loud systems in a dense protocluster environment. Through resolved spectral analysis of archival radio data up to 15 GHz, we find evidence that the inner hotspots are due to restarted activity, while the outer hotspots remain energetic, suggesting a rapid duty cycle while the jet precesses. This allowed the low-energy aged plasma driving the X-rays to remain spatially distinct from the high-energy plasma. Together, our results promise a revelation of AGN-related radio emission at high redshift using future low-frequency arrays such as SKA-LOW.

We report the discovery of variable $\gamma$-rays up to petaelectronvolt from Cygnus X-3, an iconic X-ray this http URL $\gamma$-ray signal was detected with a statistical significance of approximately 10 $\sigma$ by the Large High Altitude Air Shower Observatory (LHAASO).Its intrinsic spectral energy distribution (SED), extending from 0.06 to 3.7 PeV, shows a pronounced rise toward 1 PeV after accounting for absorption by the cosmic microwave background this http URL detected month-scale variability,together with a 3.2$\sigma$ evidence for orbital modulation, suggests that the PeV $\gamma$-rays originate within, or in close proximity to, the binary system this http URL observed energy spectrum and temporal modulation can be naturally explained by $\gamma$-ray production through photomeson processes in the innermost region of the relativistic jet, where protons need to be accelerated to tens of PeV energies.

Mauve is a low-cost small satellite developed and operated by Blue Skies Space Ltd. The payload features a 13 cm telescope connected with a fibre that feeds into a UV-Vis spectrometer. The detector covers the 200-700 nm range in a single shot, obtaining low resolution spectra at R~20-65. Mauve has launched on 28th November 2025, reaching a 510 km Low-Earth Sun-synchronous orbit. The satellite will enable UV and visible observations of a variety of stellar objects in our Galaxy, filling the gaps in the ultraviolet space-based data. The researchers that have already joined the mission have defined the science themes, observational strategy and targets that Mauve will observe in the first year of operations. To date, 10 science themes have been developed by the Mauve science collaboration for year 1, with observational strategies that include both long duration monitoring and short cadence snapshots. Here, we describe these themes and the science that Mauve will undertake in its first year of operations.

Stellar flares are an intense stellar activity that can significantly impact the atmospheric composition of the surrounding planets and even the possible existence of life. During such events, the radiative energy of the star is primarily concentrated in the optical and X-ray bands, with the X-ray flux potentially increasing by tens or even hundreds of times. Einstein Probe (EP) detected a new X-ray transient EP J2322.1-0301 on 27 September 2024. Its spatial localization shows a high positional coincidence with the nearby high proper motion K-type star PM J23221-0301. Follow-up X-ray observations confirmed the flux enhancement of the source, while optical spectroscopic monitoring revealed time-variable features, particularly the disappearance of the H-alpha emission line. This X-ray flare is consistent with a characteristic fast-rise-exponential-decay (FRED) light curve, with a rise timescale of 1.4 ks, a decay timescale of 5.7 ks, and a total duration of about 7.1 ks. The peak luminosity in the 0.5-4.0 keV energy band reached about 1.3 x 10^31 erg s^-1, with a total energy release of about 9.1 x 10^34 erg, consistent with the empirical energy correlations observed in magnetic-reconnection-driven stellar flares, as inferred from the multitemperature plasma structure and H-alpha-X-ray energy correlation. This discovery underscores EP's capability in understanding stellar magnetic activity via observing stellar transients.

In the context of future Venusian missions, it is crucial to improve our understanding of Venus upper atmosphere through 3D modeling, notably for spacecraft orbit computation. This study compares three General Circulation Models (GCMs) of the Venusian atmosphere up to the exosphere: the Venus Planetary Climate Model (Venus PCM), the Venus Thermospheric Global Model (VTGCM) and the Tohoku University GCM (TUGCM), focusing on their nominal simulations (e.g. composition, thermal structure and heating/cooling rates). Similarities and discrepancies among them are discussed in this paper, together with data-models comparison. The nominal simulations analyzed in this study fail to accurately reproduce the daytime observations of Pioneer Venus, notably overestimating the exospheric temperature. This is linked to an underestimation of the atomic oxygen (O) abundance in the three GCMs, and suggests the need of additional O production in the thermosphere. The selection of solar spectrum is also the main reason for the discrepancies between the models in terms of temperature dependence on solar activity. A list of recommendations is proposed aiming at improving the modeling of Venus' upper atmosphere, among them: 1. Standardize the EUV-UV solar spectrum input. 2. Update the near-infrared heating scheme with Venus Express-Era data. 3. Reassess Radiative cooling schemes. 4. Investigate the underestimated atomic Oxygen abundance.

The mysterious Galactic Center Excess of gamma rays could be explained by a large population of millisecond pulsars hiding in the Galactic bulge, too faint to be detected as individual high-energy point sources by the Fermi Large Area Telescope, as well as too fast and too dispersed to be detected in shallow radio pulsation surveys. Motivated by an innovative candidate selection method, we aim at detecting millisecond pulsars associated with the Galactic Center Excess by carrying deep radio pulsation searches towards promising candidates detected in the inner Galaxy, in X rays by Chandra, and in radio or gamma rays by the Very Large Array or Fermi. We conducted deep radio observation and follow-up campaigns with MeerKAT, the Murriyang and the Green Bank telescopes towards 9 X-ray candidate sources. We here report the detection of two new millisecond pulsars, including a black widow candidate, towards the Galactic bulge: PSRs J1740-2805 and J1740-28. These discoveries double the number of MSPs discovered within the innermost 2 degree from the Galactic center.

The Serpens OB2 association (l ~ 18.5 deg, b ~ 1.9 deg, d = 1950 +/- 30 pc) is a large star-forming complex ~65 pc above the Galactic midplane, with a clumpy, elongated structure extending ~50 pc parallel to the plane. We analyse a sample of probable association members, including OB stars and low-to-intermediate-mass young stellar objects (YSOs) from the SPICY catalogue. While both populations are found throughout the association, the OB stars lie preferentially on the side nearest the Galactic plane, while the YSOs are generally younger and more strongly clustered around molecular-cloud clumps detected in 13CO MWISP data. Using Gaia DR3 proper motions to probe the association's internal kinematics, we find aligned stellar velocities on length scales <2 pc, two-point velocity statistics that show increasing velocity differences and predominantly divergent motions at larger separations, and distinct velocities for star clusters within the association. Finally, the association exhibits gradual, but statistically significant global expansion perpendicular to the Galactic plane, with a spatial gradient of 0.10 +/- 0.02 km/s/pc. The expansion of the H II region Sh 2-54, powered by the association's OB stars, may be accelerating the star-forming cloud away from the plane given the system's geometry, plausibly inducing the vertical stellar velocity gradient. The clumpy stellar distribution, correlated velocities on small scales, and increasingly divergent motions on larger scales are consistent with an initial velocity field inherited from a turbulent molecular cloud modified by stellar feedback. Ser OB2 demonstrates that the multi-scale expansion of an OB association can begin even while star formation is still ongoing throughout the complex.

Simulation-based inference (SBI) with neural posterior estimation (NPE) provides rapid X-ray spectral fitting in both Gaussian and Poisson regimes by learning approximate parameter posteriors from simulations. We investigate auto-encoders for compressing high-resolution X-ray spectra, motivated by newAthena X-ray Integral Field Unit (X-IFU), and use likelihood-based importance sampling to refine NPE outputs. Our auto-encoder maps spectra to a low-dimensional latent space and is trained with a custom loss equal to the Cash statistic (C-stat) between simulated and reconstructed spectra. A neural density estimator is then trained on the latent representations. Both models are trained in multiple rounds: at each round, new simulations are drawn from a truncated proposal concentrated around the observation, improving efficiency as the proposal contracts. After NPE convergence, we apply likelihood-based importance sampling to correct the learned posterior. To assess information retention, we train a diagnostic network that predicts the original spectral parameters from the latent space, and we also train a network to learn the likelihood directly to accelerate importance sampling. On X-IFU-like simulations, the auto-encoder and multi-round NPE outperforms PCA and hand-crafted spectral summaries in accuracy and robustness. After importance sampling, the resulting posteriors are statistically indistinguishable from those obtained with nested sampling. On a standard laptop, the full pipeline (simulation, compression, inference, correction) delivers 10x speedups. We further demonstrate the approach on XRISM/Resolve and on lower-resolution NICER and XMM-Newton EPIC-pn data, confirming applicability across instruments and resolutions. Overall, NPE on compressed spectra paired with likelihood-based importance sampling offers an exact yet efficient alternative for Bayesian X-ray spectral fitting.

The existence of a population of millisecond pulsars in the Galactic bulge is supported, along with other evidence, by the Fermi GeV excess, an anomalous {\gamma}-ray emission detected almost 15 years ago in the direction of the Galactic center. However, radio surveys searching for pulsations have not yet revealed bulge millisecond pulsars. Identifying promising bulge millisecond pulsar candidates is key to motivating pointed radio pulsation searches. Candidates are often selected among steep-spectrum or polarized radio sources, but multiwavelength information can also be exploited: The aim of this work is to pinpoint strong candidates among the yet unidentified X-ray sources. We investigated the multiwavelength counterparts of sources detected by the Chandra X-ray observatory that have spectral properties expected for millisecond pulsars in the Galactic bulge. We considered that ultraviolet, optical, and strong infrared counterparts indicate that an X-ray source is not a bulge pulsar, while a radio or a faint infrared counterpart makes it a promising candidate. We identify a large population of more than a thousand X-ray sources without optical, ultraviolet, or strong infrared counterparts. Among them, five are seen for the first time in unpublished radio imaging data from the Very Large Array. We provide the list of promising candidates, for most of which follow-up pulsation searches are ongoing.

We present new visible-wavelength spectroscopic observations of the Patroclus-Menoetius binary system in the Jupiter Trojan population. Motivated by previously published spectra from different instruments that showed evidence of significant longitudinal variability, we obtained two spectra spanning 440-680 nm at near-opposite rotational phases with the Gemini Multi-Object Spectrograph on the Gemini South telescope during the late 2024 apparition. The same solar analog was used for both observations to remove one source of inconsistency. We measured spectral slopes of 2.51% $\pm$ 0.05%/100 nm and 8.13% $\pm$ 0.05%/100 nm at the two different rotational phases. The first of these measurements was serendipitously obtained during an occultation of Menoetius by Patroclus. Although the statistical significance of the spectral slope discrepancy persists even after considering possible systematic errors stemming from differences in slit position angles and air masses between the asteroid and solar analog exposures, we consider this report of variability to be tentative. We briefly explore several scenarios that could explain the measured spectral slope variability. Additional follow-up observations are necessary to definitively confirm and characterize any inhomogeneities across the surface, which will have major implications for the 2033 flyby of Patroclus-Menoetius by the Lucy spacecraft.

Current observations of neutron stars and measurements of gravitational waves only provide constraints on the zero temperature ($T=0$) equation of state (EoS) of dense matter. The detection of the post-merger gravitational-wave signal from a binary neutron star merger would additionally provide access to finite-temperature properties of the EoS which contain more information about the composition and the interactions of dense matter than the cold EoS alone. In particular deconfined quark matter may be probed by its characteristic finite temperature effects. This is especially the case for color-superconducting phases, in which the quasiparticle contribution to the thermal pressure is exponentially suppressed at low temperatures. Here we develop a new finite $T$ framework to model the thermal EoS for dense quark matter based on the cold quark matter EoS which is useful for numerical relativity simulations. We test the validity of the framework against a three-flavor NJL mean-field calculation, both with and without diquark pairing. We find that even for the complicated phase diagram of the NJL model including multiple different phases the framework is accurate to the few percent level for temperatures up to $T\sim 50\,$MeV.

We report the unprecedented Ly$\alpha$ properties of AMORE6, an extremely metal-poor ($12+\log({\rm O/H}) < 6$), low-mass ($M_\star = 4.4\times10^{5}\,M_\odot$), and ultra-compact (effective radius $\lesssim30$ pc) dwarf galaxy at $z=5.7253$, gravitationally lensed by the cluster A2744. A prominent, narrow, and nearly-symmetric Ly$\alpha$ emission line is detected at the systemic redshift (the latter traced by H$\beta$, from JWST/NIRCam slitless spectroscopy), with rest-frame $EW=150\pm10$ Å, $\rm FWHM=58\pm1$ km s$^{-1}$, and a slight asymmetry, resulting in a $\rm \sim10\%$ flux excess in the red wing of the line. The negligible velocity offset from systemic ($dv = 4\pm67$ km s$^{-1}$, $3\sigma$ uncertainty), together with the sharpness and symmetry of the profile, indicates minimal radiative transfer effects implying a neutral hydrogen column density consistent with an optically thin medium, compatible with a non-zero ionizing photon escape fraction. If indirect spectral diagnostics calibrated at $z<4.5$ remain the only viable tools to identify LyC leakers during reionization, then based on its strongest indicator (Ly$\alpha$), AMORE6 stands out as one of the most compelling LyC-leaking candidates yet discovered in the reionization epoch.

Our study shows that the more a dark energy model allows the equation of state to become phantom-like toward the past, the better it fits current observations which also may lead to stronger statistical tension with $\Lambda$CDM. In particular, models with greater flexibility in their redshift evolution, such as power-law type forms, are mildly preferred over more restrictive parametrizations with similar parameter dimensionality. For the dataset including DESY5, the deviation from $\Lambda$CDM reaches a significance of $6.26\sigma$. These results indicate that present data are sensitive to enhanced phantom behaviour and provide clear guidance for constructing phenomenologically viable dark energy parametrizations.

The Wide-field Infrared Transient Explorer (WINTER) is a new near-infrared time-domain survey instrument installed on a dedicated 1-meter robotic telescope at Palomar Observatory. The project takes advantage of the recent technology advances in time-domain astronomy, robotic telescopes, large-format sensors, and rapid data reduction and alert software for timely follow up of events. Since June of 2023, WINTER robotically surveys the sky each night to a median depth of J_AB = 18.5 mag, balancing a variety of science programs including searching for kilonovae from gravitational-wave alerts, blind surveys to study galactic and extragalactic transients and variables, and building up reference images of the near-infrared sky. The project also serves as a technology demonstration for new large-format Indium Gallium Arsenide (InGaAs) sensors for near-infrared photometry without cryogenic cooling. WINTER's custom camera combines six InGaAs sensors with a novel tiled fly's-eye optical design to cover a >1 degree-squared field of view with 1 arcsecond pixels in the Y-, J-, and shortened-H-band filters (0.9 - 1.7 micron). This paper presents the design, performance, and early on-sky science of the WINTER observatory.

The discovery of exoplanets in binary star systems-now numbering about 850 of the nearly 4,600 known exoplanet systems-raises questions about whether observational bias or stellar companions inhibit planet formation. While most studies on terrestrial planet formation assume planar configurations, wide binaries likely feature random inclinations, potentially disrupting planet-forming disks. This study explores the evolution of embryo-planetesimal disks in S-type motion in misaligned binary systems, focusing on the stage after the gas phase when terrestrial planet formation begins and gravitational interactions dominate. Using our GPU-accelerated N-body code GANBISS, we simulate disks with 2,000 planetesimals and 25 planetary embryos, studying the influence of the planetesimals on the evolution of the embryos and tracking their growth through collisions. After the simulations, we analyse collision outcomes with an analytical model. Moreover, for certain inclined binary configurations, we compare dynamically excited (perturbed by the secondary star) with cold disks in inclined configurations, as the distribution after the gas phase in misaligned binaries remains unclear. Our simulations reveal two key outcomes: (i) embryos migrate slightly inward in misaligned systems, and (ii) The initial large oscillations in embryos' inclinations and nodes around the respective values of the secondary star dampen over time. Collision analysis shows distinct differences: planar systems favour accretive collisions, while inclined configurations exhibit more destructive events. These findings underscore the sensitivity of planet formation dynamics to binary star alignment and initial disk conditions.

Accurate modelling of the effective point spread function (ePSF) is essential for high-precision photometry and astrometry, particularly in undersampled imaging regimes. In this work, we build on a well-established ePSF modelling framework and its commonly used open-source Python implementation and demonstrate that several simple but effective modifications to existing ePSF modelling routines can significantly improve model accuracy. We use synthetic ePSFs to generate simulated datasets of stellar images, allowing us to evaluate the accuracy of ePSF models and determine the scale of the pixel-phase errors in resulting flux and position measurements. We systematically investigate how specific modelling choices affect ePSF accuracy, and evaluate the influence of oversampling, interpolation, gridpoint estimation, smoothing, star-sample distribution, and dithering on photometric precision. We apply our refined ePSF modelling routine to images from the Global Jet Watch observatories, demonstrating its improved ability to recover an accurate ePSF for real astronomical images. Our findings highlight the importance of tailoring the modelling approach to the specific characteristics of the instrument and detector, as well as to the nature of the available imaging data used to construct the ePSF model. These results provide practical guidance for optimising ePSF construction, thereby improving the reliability of photometric and astrometric measurements.

Recent cosmological measurements are hinting that dark energy may evolve, with its equation of state, $w_\mathrm{DE}$, even showing oscillatory patterns. In this work, we employ a model-independent approach to jointly reconstruct $w_\mathrm{DE}$ and the sum of neutrino masses, $\sum m_\nu$, adopting the PCHIP method with seven fixed nodes in which we allow the two parameters to vary. We employ CMB, Baryon Acoustic Oscillations and Supernovae Ia data to constrain the values of $w_\mathrm{DE}$ and $\sum m_\nu$ at each node. We conduct three different analyses in which we reconstruct $w_\mathrm{DE}$: one with fixed $\sum m_\nu=0.06~\mathrm{eV}$; one in which we allow $\sum m_\nu$ to vary, and one in which we also reconstruct $\sum m_\nu$ using the PCHIP method. We find the dark energy equation of state to be consistent with the cosmological constant scenario, except when including DESI data and allowing for phantom crossing, where we find a $95\%$ CL deviation from $w_\mathrm{DE}=-1$ around $z\sim1.2$. For neutrino masses, we obtain looser constraints when focusing on phantom dark energy, that show further early and late relaxation when reconstructing the mass via the PCHIP method.

An excess emission has been observed by Spitzer in the [3, 5] micron range of the SNR 1987A spectrum. It is generally argued that this excess could be due to the presence of warm amorphous carbon dust in the equatorial ring (ER) around the supernova, but the proposed models all have problems. This prompted us to present an alternative view on the interpretation of the Spectral Energy Distribution (SED) of SNR 1987A from the near-IR wavelengths to the radio frequencies (from 3 micron up to 1.4 GHz), between 6000 and 8000 days after outburst. We argue that the origin of that excess could be attributed instead to a free-free emission. We show that under very specific conditions (the free-free is self-absorbed at a cut-off frequency imposed by the mass of the emitting region), it could be produced by collisional heating of the gas. We then discuss the time evolution of the various components of the SED. We establish a linear relationship between the growth of the warm carbon dust mass and that of the silicates dust during the analyzed period. Finally, we build the Spitzer light curves and we show that our models reproduce the observations pretty well, although our study clearly favors the free-free case. In conclusion, we argue that the free-free model provides a formally very good description of the data, however the model does require some very specific parameter choices, and results in an unusually low temperature for the ionized gas.

We show that atomic clock measurements provides an exceptionally sensitive Solar System probe of scalar tensor dark energy. By connecting variations in Newton's constant and differential clock drifts to the dynamics of a single dark energy scalar, we derive a direct constraint on the present day equation of state and our results force any locally coupled scalar dark energy into a very slowly rolling regime, giving the strongest bounds on the equation of state parameter. This is independent of potential shape or kinetic structure and rules out broad classes of canonical and non canonical models, leaving only near Lambda CDM behavior or fully decoupled fields as viable late time scalar dark energy, thereby leaving cosmological constant and minimally coupled scalar field models as the most consistent dark energy regimes. We also use results from Lunar Laser Ranging and photon trajectories to further strengthen our the depth of our constraints.

Dual active galactic nuclei (DAGN) mark an observable stage of massive black hole (MBH) pairing in galaxy mergers and are precursors to the MBH binaries that generate low-frequency gravitational waves. Using the large-volume ASTRID cosmological simulation, we construct DAGN catalogs matched to current (COSMOS-Web, DESI) and forthcoming (AXIS, Roman) searches. With realistic selection functions applied, ASTRID reproduces observed dual fractions, separations, and host-galaxy properties across redshifts. We predict a substantial population of small-separation (<5 kpc) duals that current surveys fail to capture, indicating that the apparent paucity of sub-kpc systems in COSMOS-Web is driven primarily by selection effects rather than a physical deficit. By following each simulated dual forward in time, we show that dual AGN are robust tracers of MBH mergers: ~30-70% coalesce within $\lesssim 1$ Gyr, and 20-60% of these mergers produce gravitational-wave signals detectable by LISA. Duals accessible to AXIS and Roman are the progenitors of ~10% of low-redshift LISA events and ~30% of the PTA-band stochastic background. Massive green-valley galaxies with moderate-luminosity AGN, together with massive star-forming hosts containing bright quasars at $z>1$, emerge as the most likely environments for imminent MBH binaries. These results provide a unified cosmological framework linking dual AGN demographics, MBH binary formation, and gravitational-wave emission, and they identify concrete, high-priority targets for coordinated electromagnetic and GW searches in upcoming multi-messenger surveys.

We performed a search for fast X-ray transients (FXTs), with durations longer than one second and less than one day, through data of the Wide Field Camera (WFC) instrument onboard the BeppoSAX X-ray observatory collected between June 1996 and April 2002. (..) We focused our search on gamma-ray bursts (GRBs), X-ray flashes (XRFs), X-ray flares from high-mass X-ray binaries and stellar flares, while Type-I and II X-ray bursts from Galactic neutron stars were excluded. 149 such fast transient events were detected. 63 of these are new to the literature. 38 flares are identified with 22 nearby stars. Three stars have never been seen flaring before in X-rays or optical (NLTT 51688, GR Dra and UCAC4 255-003783). We find that the MeV transient GRO J1753+57 is most likely the same object as GR Dra rather than an AGN as previously thought. Eleven flares were detected from known high-mass X-ray binaries with irregular wind accretion (four of which are of the subclass of supergiant fast X-ray transients). 100 GRBs were identified of which 24 have not been published before. We classify 37% of the X-ray detected GRBs as XRFs with relatively large X-ray to gamma-ray flux ratio, gamma-rays being measured with the BeppoSAX Gamma Ray Burst Monitor. The duration/spectral hardness distribution of all FXTs is bimodal, separating the group roughly in transients shorter and longer than 1 ksec and with relatively hard and soft spectra, respectively. We identify the 'short' FXTs as GRBs and XRFs and the `long' FXTs as flares from nearby late-type stars and X-ray binaries. The BeppoSAX-WFC FXT sample is found to be consistent with the one observed by Einstein Probe, when the sensitivity of the two instruments is taken into account.

We present the first targeted centimeter-band radio observations of two recently-discovered exoplanet systems that are prime candidates for magnetic star-planet interaction (SPI): TOI-540 and SPECULOOS-3. The targets were selected due to the small orbital separation of their known planets, as well as for indications of stellar magnetic activity, given that for SPI radio emission may be strongest when a sufficiently magnetized star hosts a close-in planet. The deep, multi-epoch Very Large Array (SPECULOOS-3) and MeerKAT (TOI-540) observations yield non-detections, with $3\sigma$ limits of $\lesssim 7.5$ $\mu$Jy ($4-8$ GHz) and $\lesssim 30-80$ $\mu$Jy ($0.8-1.7$ GHz), respectively. For SPECULOOS-3 b we rule out observable SPI for most of its orbit, while for TOI-540 b we sample a narrower range, around planetary transit. We model possible planetary magnetic field strength constraints for both systems, and conclude that our observations are sensitive enough to sample SPI emission in these systems if present and directed at us, even for a planetary field of only $\sim 1$ G.

Gravitational waves (GWs) are unique messengers as they travel through the Universe without alteration except for gravitational lensing. Their long wavelengths make them susceptible to diffraction by cosmic structures, providing an unprecedented opportunity to map dark matter substructures. Identifying lensed events requires the analysis of thousands to millions of simulated events to reach high statistical significances. This is computationally prohibitive with standard GW parameter estimation methods. We build on top of state-of-the-art neural posterior algorithms to accelerate the lensed inference from CPU days to minutes with DINGO-lensing. We showcase its capabilities by reanalyzing GW231123, the most promising lensed candidate so far, and find that its statistical significance cannot exceed 4$\sigma$. We observe that 8% of GW231123-like nonlensed simulations favor lensing, which could be explained by the self-similarity of short-duration signals. Still, 58% of GW231123-like lensed simulations have larger support for lensing, showing that higher detection statistics are possible. Although GW231123 exposes the challenges of claiming the first GW lensing detection, our deep-learning methods have demonstrated to be powerful enough to enable the upcoming discovery of lensed GWs.

Matter-induced neutrino flavor mixing (the Mikheyev-Smirnov-Wolfenstein, or MSW, effect) is a central prediction of the neutrino mixing framework, but it has not been conclusively observed. Direct observation of the energy-dependent MSW transition in the solar electron-neutrino survival probability would solve this, but backgrounds have been prohibitive. We show that our new technique for suppressing muon-induced spallation backgrounds will allow JUNO to measure the MSW transition at $>$4$\sigma$ significance in 10 years. This would strongly support upcoming multi-\$1B next-generation long-baseline experiments and their goals in cementing the neutrino mixing framework.

A pulsar, i.e., a spinning neutron star, with a deformation could emit gravitational waves continuously. Such continuous waves, which have not been detected yet, will be very useful to study gravitational physics and to probe the extreme physics of neutron stars. While typically such waves from a pulsar are estimated considering an overall stellar ellipticity, there can be multiple irregularities or mountains in the stellar crust that the gravity of the star cannot smooth. In this paper, we consider this realistic situation and compute the strain, power, torque and the pulsar spin-down rate due to multiple mountains supported by the stellar crust. Here, we consider astronomically motivated mountain distributions and use the Brans-Dicke theory of gravity which has three polarization states: two tensors dominated by the time-varying quadrupole moment and one scalar dominated by the time-varying dipole moment. We also give the limiting results for general relativity.

Interpreting multi-messenger signals from neutron stars and black holes requires reliable general-relativistic magnetohydrodynamics (GRMHD) simulations across rapidly evolving high-performance-computing platforms, yet key algorithms are routinely rewritten within infrastructure-specific numerical-relativity codes, hindering verification and reuse. We present the General Relativistic Hydrodynamics Library (GRHayL), a modular, infrastructure-agnostic GR(M)HD library providing conservative-to-primitive recovery, reconstruction, flux/source and induction operators, equations of state, and neutrino leakage through an intuitive interface. GRHayL refactors and extends the mature IllinoisGRMHD code into reusable pointwise and stencil-wise kernels, enabling rapid development and cross-code validation in diverse frameworks, while easing adoption of new microphysics and future accelerators. We implement the same kernels in the Einstein Toolkit (Carpet and CarpetX) and BlackHoles@Home, demonstrating portability with minimal duplication. Validation combines continuous-integration unit tests with cross-infrastructure comparisons of analytic GRMHD Riemann problems, dynamical Tolman-Oppenheimer-Volkoff evolutions, and binary neutron-star mergers, showing comparable or improved behavior over legacy IllinoisGRMHD and established Einstein Toolkit codes.

Axion-like particles can be abundantly produced through scattering processes in the cores of neutron stars (NSs). If they are ultralight ($m_a \lesssim 10^{-4}$ eV), then they can efficiently convert to detectable photons in the external NS magnetospheres, and if they are heavy ($m_a \gtrsim 1$ eV), then they can decay into photons before reaching Earth. In this work, we search for the resulting X-ray signatures from both of these channels summing over the $\textit{cosmological}$ NS population. We compare the predicted axion-induced X-ray signal to the cosmic X-ray background today as measured by a number of instruments such as NuSTAR, HEAO, Swift, and INTEGRAL. We model the axion-induced signal using NS cooling simulations and magnetic field evolution models. We find no evidence for axions and derive strong constraints for both ultralight and heavy axion scenarios, covering new parameter space for the axion-photon and axion-nucleon couplings. Our results rule out the axion-explanation of the Magnificent Seven X-ray excess from nearby isolated NSs.

The inert doublet model is a two-Higgs-doublet extension of the standard model that provides a minimal and versatile framework for frozen-out dark matter. Assuming standard cosmology, if the dark matter mass ranges between approximately 120 GeV and 500 GeV then it turns out to be underabundant, as gauge interactions render its annihilation too efficient. In this work, we show that this mass window becomes allowed in cosmological scenarios where dark matter freeze-out occurs during a period with a stiff equation of state, $w > 1/3$, such as kination. This predictive setup satisfies all current experimental constraints while remaining within the reach of upcoming detection efforts.

It was recently shown that, in a binary coalescence, the greybody factor of the remnant black hole modulates the post-merger ringdown signal. In this work, we demonstrate that a simple four-parameter model based on the greybody factor accurately reproduces the frequency-domain amplitude of a large set of comparable-mass, aligned-spin numerical relativity waveforms from the SXS catalog, achieving mismatches of order ${\cal O}(10^{-5})$ and improving existing models by roughly two orders of magnitude. We also identify the optimal initial frequency for applying the model in the frequency domain and provide analytical fits of the model parameters in terms of the progenitor masses and aligned spins. Our results pave the way for new consistency tests of the ringdown phase, complementary to traditional black hole spectroscopy.

Thermal distribution functions can only be of the Fermi-Dirac or Bose-Einstein types, whereas distorted spectra encompass any possible deviations from these shapes. It is fruitful to devise parameterizations of these distortions with only a few parameters which depend on the physical system considered. A method proposed by Stebbins consists in describing a distorted spectrum as a sum of thermalized spectra with a distribution of temperatures, the moments of which are the parameters of interest. After revisiting and extending this approach by working at the level of the number density distribution instead of the standard spectrum, we build another method which consists in describing the distorted spectrum by a polynomial modulating a reference thermalized spectrum. The distortion parameters are then the coefficients of a decomposition on a suitable orthonormal polynomial basis. We advocate that the latter is computationally easier and allows to describe a wide range of distortions. With this formalism, we efficiently describe the standard distortions of the cosmological backgrounds of neutrinos and photons, and we obtain model-independent constraints on nonstandard distortions of these cosmological relics.

Symbolic regression (SR) has emerged as a powerful method for uncovering interpretable mathematical relationships from data, offering a novel route to both scientific discovery and efficient empirical modelling. This article introduces the Special Issue on Symbolic Regression for the Physical Sciences, motivated by the Royal Society discussion meeting held in April 2025. The contributions collected here span applications from automated equation discovery and emergent-phenomena modelling to the construction of compact emulators for computationally expensive simulations. The introductory review outlines the conceptual foundations of SR, contrasts it with conventional regression approaches, and surveys its main use cases in the physical sciences, including the derivation of effective theories, empirical functional forms and surrogate models. We summarise methodological considerations such as search-space design, operator selection, complexity control, feature selection, and integration with modern AI approaches. We also highlight ongoing challenges, including scalability, robustness to noise, overfitting and computational complexity. Finally we emphasise emerging directions, particularly the incorporation of symmetry constraints, asymptotic behaviour and other theoretical information. Taken together, the papers in this Special Issue illustrate the accelerating progress of SR and its growing relevance across the physical sciences.

The production of $^{60}$Fe is crucial for nucleosynthesis in massive stars and supernovae. In this work, by using the microscopic EP+IPM (exact pairing plus the independent-particle model) for the nuclear level density (NLD) and extended EP+PDM (exact pairing plus phonon damping model) for the $\gamma$-ray strength function (gSF), we re-evaluate the substantial enhancement of $^{60}$Fe production recently reported in {\it A. Spyrou et al., Nat. Comm. {\bf 15}, 9608 (2024)}, which was attributed to an unexpectedly large Maxwellian-averaged cross section (MACS). Our analysis demonstrates that this enhancement indeed originates from the choice of NLD, which, despite being constrained to reproduce the total NLD and gSF data, lacks a reliable spin dependence, a critical input for Hauser-Feshbach calculations of nuclear reaction rate. In contrast, our predictions yield a significantly lower MACS, calling the claimed enhancement into question. In particular, our approach highlights the microscopic nature of the low-energy enhancement of the gSF, the so-called upbend resonance, which arises from strong particle-particle ($pp$) and hole-hole ($hh$) excitations that emerge only at finite temperature, thereby further reinsisting on the invalidity of the Brink-Axel hypothesis in this low-energy region. Overall, our study reopens the question on the long-standing problem of $^{60}$Fe production in massive stars.

The dipole response of a nuclear system, characterized by its photon strength function (PSF), is a key ingredient of many applications of nuclear structure, ranging from nuclear reactor design and nuclear waste transmutation to astrophysical models of nucleosynthesis and stellar evolution. While the majority of those applications require the knowledge of PSF of mid-mass and heavy nuclei, there is now renewed interest in $E1$ strength distributions of light nuclei in the framework of the PANDORA project, which aims at an understanding of the mass distribution of ultrahigh-energy cosmic radiation (UHECR).UHECR is of extragalactic origin and its interaction along the travel path is dominated by photoabsorption of cosmic background radiation boosted to the Giant Dipole Resonance (GDR) energy region in the center-of-mass system. Thus, systematic knowledge of the photoabsorption cross sections in light nuclei and of their subsequent particle decay is required. The purpose of this work is to enhance the database of available theoretical evaluations of PSF of light nuclei that are necessary in the studies of UHECR propagation. We employ the Configuration Interaction Shell Model (CI-SM) approach to provide predictions of $E1$ dipole response for $p$ and $sd$-shell nuclei, with mass number $A$ between 7 and 40. Theoretical predictions are compared to available data and to existing predictions from phenomenological and microscopic models. Finally, the impact of using of CI-SM PSF on the predicted propagation of a $^{40}$Ca UHECR source is studied.

We examine the motion of an electron constrained to follow a magnetic field line near a primordial sub-stellar mass black hole. Earlier studies treated the problem in flat (Minkowski) spacetime, yielding qualitatively correct results and introducing a light cylinder (LC), a hypothetical surface where the linear velocity of rotation equals the speed of light. However, this picture changes significantly when gravity is included. By analyzing the electron's dynamics in the Schwarzschild metric, we obtain a modified light cylinder (MLC) whose geometry no longer resembles a cylinder. We then determine the maximum energies attainable by the electrons under the limiting effects of inverse Compton scattering, curvature radiation, and synchrotron radiation.

Gravitational waves can be gravitationally lensed by massive objects along their path. Depending on the lens mass and the lens-source geometry, this can lead to the observation of a single distorted signal or multiple repeated events with the same frequency evolution. We present the results for gravitational-wave lensing searches on the data from the first part of the fourth LIGO-Virgo-KAGRA observing run (O4a). We search for strongly lensed events in the newly acquired data by (1) searching for an overall phase shift present in an image formed at a saddle point of the lens potential, (2) looking for pairs of detected candidates with consistent frequency evolution, and (3) identifying sub-threshold counterpart candidates to the detected signals. Beyond strong lensing, we also look for lensing-induced distortions in all detected signals using an isolated point-mass model. We do not find evidence for strongly lensed gravitational-wave signals and use this result to constrain the rate of detectable strongly lensed events and the merger rate density of binary black holes at high redshift. In the search for single distorted lensed signals, we find one outlier: GW231123_135430, for which we report more detailed investigations. While this event is interesting, the associated waveform uncertainties make its interpretation complicated, and future observations of the populations of binary black holes and of gravitational lenses will help determine the probability that this event could be lensed.

This work investigates the nonlinearity of the power-law model of F(T) gravity, highlighting the inability of the Boltzmann solver CLASS to handle nonlinear models. As a workaround, a second-order Taylor expansion is applied to the nonlinear field equations, under the assumption that the extra degree of freedom n, which quantifies deviations from the currently favored cosmological model (\Lambda CDM), remains sufficiently small to preserve the key properties of the \Lambda CDM model. The validity of the Taylor expansion is supported by supernova data indicating n \leq 0.05, for which the power spectrum can be accurately computed within CLASS with a negligible truncation error.

We investigate the nuclear equation of state (EoS) for isospin-asymmetric matter using a new set of RMF interactions with the $\sigma$-$\delta$ and $\omega$-$\rho$ mixing, referred to as the OMEG family. These interactions are optimized so as to reproduce both terrestrial nuclear measurements and astrophysical constraints extracted from NICER and GW170817. The $\sigma$-$\delta$ mixing softens the nuclear symmetry energy and pressure around twice the saturation density, which enables relatively small neutron-star radii and tidal deformabilities while keeping the nuclear EoS sufficiently stiff at high densities to support $2M_{\odot}$ neutron stars. We find that the curvature parameter, $K_{\textrm{sym}}$, plays an important role in realizing the soft-to-hard behavior of the nuclear EoS, and the astrophysical data favor small or even negative values of $K_{\textrm{sym}}$.

We derive an extended expression for the relaxation time of a barotropic Israel-Stewart (IS) fluid using the non-linear causality constraint, and propose a new formulation for modeling causal viscous dissipation in barotropic fluids. With this generalized relaxation time, the non-linear IS equation simplifies to a first-order non-linear expression connecting bulk viscous pressure and energy density, which remains valid in any homogeneous and isotropic spacetime. In the case of spatially flat Friedmann universe, adopting this extended relation in the generalized non-linear IS theory, provides new class of analytical solutions in both, the linear, and the non-linear regimes. We also find that, the resulting effective equation of state in the linear regime naturally reproduces the generalized polytropic form which is often introduced phenomenologically in literature. Resulting dynamical implications are investigated and the constraints necessary for ensuring an acceptable evolutionary behavior for the fluid are determined. A detailed dynamical system analysis of the coupled Einstein-Israel-Stewart (EIS) system is also performed. Finally, we solve the coupled EIS equations numerically, and show that the model can support a transient Hubble slow-roll expansion phase with a smooth exit to a radiation-dominated universe, which is challenging to obtain in standard inflationary models.

In this study, we show that, in the background of the primordial magnetic field, the CME effect can significantly amplify the chiral chemical potential sourced by the CP violation near the bubble walls during the first-order electroweak phase transition. This effect can lift the generated baryon asymmetry by several orders, and make it possible to explain the baryon asymmetry of the Universe with a CPV in the fermion sector far beyond the limitation of the electron dipole moment.

For axions present during inflation, it has been shown that a non-minimal coupling $\xi_\sigma$ of the inflaton to gravity worsens isocurvature bounds, while a non-minimal coupling $\xi_\rho$ of the radial Peccei-Quinn field can alleviate them. We analyze the simultaneous presence of both couplings and determine when one effect dominates the other, in both the metric and Palatini formulations of gravity. The two tendencies interpolate smoothly, but introducing a non-minimal inflaton coupling reduces the viable interval of $\xi_\rho$ in which isocurvature bounds can be alleviated while avoiding backreaction on the inflationary dynamics. We illustrate our findings in Palatini Higgs inflation and Starobinsky inflation.

An interesting feature of a cosmological phase transition can be a stage of exponential expansion (supercooling). The modified expansion history and the entropy injection at reheating, can affect the final energy fraction of dark matter. In this paper, we revisit the calculation of the freeze-out and freeze-in dynamics, showing additional effects on top of the standard dilution factor if the dark matter production is completed during the supercooling stage. We show for the first time how these effects can be particularly interesting for direct detection, as the parameter space for WIMP-like candidates shifts from excluded to allowed regions, and freeze-in candidates get closer to experimental reach. A phenomenological motivation to consider supercooling is the associated gravitational wave background. The implications of a finite-duration reheating stage, when the equation of state is close to matter-domination, are a peculiar low-frequency spectrum, and its shift to lower frequencies. These effects are a complementary test of the dynamics that we study for dark matter production, and remarkably can link direct detection of dark matter and gravitational wave astronomy.

We propose a new mechanism for inflationary model building in the framework of metric-affine gravity. Such a mechanism involves an inflaton non-minimally coupled with the Holst invariant. If the non-minimal coupling function has a zero point and it is very steep at that same point, then the canonically normalized inflaton potential always features an exponential plateau, regardless of the shape of the original inflaton potential. The inflationary predictions in such a region are equivalent to the ones of Starobinsky inflation.

We investigate the velocity-dependent one-scale (VOS) model to the case of one cosmic F-string and two D-strings as color flux tubes in pure Spin($4N$) gauge theory. We analytically calculate the scaling string density as a function of the reconnection probabilities, and confirm our results with numerical calculations. We also determine the timescale at which the string density reaches the scaling regime, and find that for certain values of the reconnection probability, the scaling time can become extremely large, by many orders of magnitude. This leads to a characteristic suppression signature of the gravitational-wave signal at high frequencies, which may become observable in the frequency range of future interferometric gravitational-wave observations.

Axions and axion-like particles generically appear in extensions of the Standard Model. While many searches assume only a single axion species, there may exist a whole spectrum of multiple such fields. We develop general formulas for axion-photon oscillations in the presence of multiple axions and analyze the implications for experimental searches, including light-shining-through-a-wall experiments, helioscopes and haloscopes. We demonstrate that axion multiplicity can qualitatively alter observational signatures, particularly through coherence and interference effects. Multiple axions can not only enhance signals compared to single axion scenarios, but also suppress them. We show that variations of experimental parameters and searches allow identifying contributions of multiple axions and obtaining information about their properties.