Papers for Monday, Dec 15 2025

Interstellar ices are mainly composed of amorphous solid water (ASW) containing small amounts of hypervolatiles, such as O2, whose diffusion-limited reactions play a key role in space chemistry. Although O2 is an important precursor molecule present during the early stages of ice formation, its surface diffusion in ASW remains poorly constrained. In this study, we experimentally investigate the surface diffusion and the entrapment efficiency of O2 in porous ASW under astrophysically relevant conditions. Experiments were conducted in an ultra-high vacuum chamber and monitored using infrared (IR) spectroscopy and quadrupole mass spectrometry. Diffusion coefficients were extracted through a novel approach applicable to IR-inactive molecules, by fitting the mass spectrometer signal during the isothermal phase with a Fickian model. These coefficients were then used to derive the diffusion energy barrier of O2 in ASW. Entrapment efficiencies were measured by analyzing the subsequent temperature-programmed desorption phase. We measured the surface diffusion coefficients at different temperatures (35 K, 40 K, 45 K) and water ice coverages (40 ML, 60 ML, 80 ML), yielding values on the order of 1E-16 to 1E-15 cm^2 s^-1. From these values, we derived a diffusion energy barrier of ED = 10 +/- 3 meV (116 +/- 35 K), corresponding to a chi ratio of about 0.1. Entrapment measurements revealed that a residual amount of about 20% of O2 remains trapped in the ASW matrix at the highest temperatures investigated. This work demonstrates that the surface diffusion of IR-inactive molecules can be experimentally quantified using mass spectrometry. Our findings show that O2 exhibits a low diffusion barrier, indicating high mobility in interstellar water ices. Moreover, we suggest these water ices likely retain a residual fraction of hypervolatiles entrapped within their structure.

Many stars -- if they could be imaged with enough angular resolution -- would exhibit features expected from theory but not possible to extract from spectra. We may group these by increasing complexity as follows. First, smooth variations in brightness across the surface, resembling solar limb darkening but much more prominent and involving more processes in stars with fast spin or external tides. Next, there are periodic features: not only oscillations, but also convective cells and starspots, which appear to transit across a star as its spins, and exoplanets that really do transit across the star. Then, there are transients like flares. Current optical interferometers provide synthetic apertures of a few hundred metres and angular resolutions down to about nanoradian ($\simeq 0.2\,$milliarcsecond), enough to resolve some of the above features on the nearest upper main-sequence stars, giants and supergiants. Ongoing projects aims to km-scale synthetic apertures, enough to measure the radius of the nearest white dwarf. In this White Paper we briefly discuss what could be observed with synthetic apertures over $\sim20\,$km -- resolving detail on white dwarfs at the level currently possible on supergiants.

Inflationary models predicting abundant primordial black holes (PBHs) and large amplitude of scalar-induced gravitational waves (SIGWs) often rely on amplified fluctuations over limited scales, typically driven by phase transitions, particle production, or departures from slow-roll evolution. While the power spectrum of these models has been extensively studied, higher-order correlations are much less understood. Motivated by the complex physics involved and the fact that PBH and SIGW formation are both sensitive to non-linearities, we present a detailed study of the bispectrum as the leading non-linear effect in these scenarios. We refer to the scale- and shape-dependence of the bispectrum collectively as its geometry; and define a scale-dependent shape correlator to disentangle the two dependencies. Generally, we find that for the scales most affected by phase transitions and particle production (including the power spectrum peak), the bispectrum is strongest near the equilateral configuration, while non-attractor phases tend to produce correlations near the squeezed configuration. We further propose a simplified bispectrum estimator, resembling local-type non-Gaussianity but with scale-dependent amplitude, that captures the main features of the full bispectrum. As an implication of our results, we show that incorporating the bispectrum significantly broadens the range of scales with a substantial probability of large smoothed density contrasts compared to linear analysis. This suggests that non-linearities can alter not only PBH abundance and SIGW amplitude but also their mass and frequency spectra. In particular, and in contrast with the usual assumption, our results hint that the second-highest peak of the power spectrum may produce more PBHs than the highest peak.

The galaxy cluster XLSSC 122 is a rare system at $z = 1.98$, hosting surprisingly evolved member galaxies when the Universe was only one-third of its present age. Leveraging deep JWST/NIRCam imaging, we perform a weak-lensing analysis and reconstruct the cluster's mass distribution, finding a mass peak that coincides with both the X-ray peak and the position of the brightest cluster galaxy. We obtain mass and concentration estimates of $M_{200\rm c}=1.6 \pm 0.3\times 10^{14}~M_{\odot}$ and $c_{200 \rm c}=6.3 \pm 0.3$, respectively, in agreement with recent strong-lensing estimates. The high concentration in particular motivates tests against empirical and simulation-derived mass--concentration relations. Placing our weak-lensing mass map in the context of Chandra X-ray data, MeerKAT radio imaging, ALMA+ACA/ACT Sunyaev-Zel'dovich (SZ) mapping, and a new JWST intracluster light measurement, we identify consistent NE--SW elongation across datasets and a pronounced offset along the same axis between the SZ and mass/X-ray peaks, pointing to significant merger activity. XLSSC 122 thus serves as a JWST pilot study for high-$z$ lensing, demonstrating the telescope's unique ability to map cluster mass distributions at $z\sim2$ and motivating a uniform sample of analogous systems with joint lensing, X-ray, SZ, and radio data to probe cluster assembly at Cosmic Noon.

Hot, X-ray emitting atmospheres pervading galaxy clusters are rich in metals, which have been synthesised and released by asymptotic giant branch (AGB) stars, core-collapse supernovae (SNcc) and Type Ia supernovae (SNIa) over cosmic history. This makes the intracluster medium (ICM) an ideal astrophysical system to constrain its chemical composition, hence ultimately understand metal production and enrichment on megaparsec scales. In this work, we take advantage of the unprecedented ~5 eV resolution offered by XRISM/Resolve to measure the chemical composition of the core of the bright, nearby, and metal-rich Centaurus cluster (287 ks) with unprecedented accuracy. We use these measurements to provide constraints on the stellar populations having enriched the cluster core. We derived the Fe abundance and its relative Si/Fe, S/Fe, Ar/Fe, Ca/Fe, Cr/Fe, Mn/Fe, and Ni/Fe ratios. We completed this high-resolution view with N/Fe, O/Fe, Ne/Fe, and Mg/Fe ratios obtained with XMM-Newton/RGS archival data. Similarly to the core of Perseus, we find that nine out of our 11 measured abundance ratios are formally consistent with the chemical composition of our Solar System. However, the (super-solar) N/Fe and (half-solar) Mg/Fe ratios significantly differ from Perseus and/or other systems, thus provide tension with the picture of a fully solar composition ubiquitous to all systems. In addition, possible uncertainties in O/Fe and Ne/Fe with atomic codes highlight the need for studying more systems at high spectral resolution to assess (or rule out) the universality of the ICM composition in clusters' cool cores. Combinations of (AGB+)SNcc+SNIa yield models can reproduce our observed X/Fe ratios in all cases. However, whether two distinct populations of SNIa are needed depends on the weight of our RGS measurements. We also briefly discuss the possibility of a multi-metallicity gas phase in this respect.

Modeling binary neutron star merger (BNSM) ejecta evolution requires simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. We systematically investigate different treatments of nuclear heating, thermalization, and opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Ejecta from long-term numerical-relativity simulations are evolved to ~30 days using a 2D ray-by-ray approach. We compare simplified heating-rates, thermalization prescriptions, and gray opacities with in-situ nuclear networks (NN) that track energy deposition, and include a composition-dependent thermalization scheme and frequency-dependent, atomic-physics-based opacities. Coupling NN and hydrodynamics affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters the abundance evolution and produces a narrower second $r$-process peak and a third peak shifted to higher mass numbers. Nuclear heating back-reaction delays and reddens the early emission. A constant thermalization underestimates the early luminosity and overestimates the late emission. Analytical opacities yield dimmer and redder kilonovae at early times ($t\lesssim$ hour) and a prolonged emission at $t\gtrsim5$ days. Resolving the first hundreds of milliseconds of hydrodynamics is essential for robust nucleosynthesis calculations, and composition-dependent thermalization and frequency-dependent, atomic opacities are needed to accurately capture the ejecta temperature and kilonova brightness and color evolution. Analytic nuclear-power fits with simplified thermalization and opacities can reproduce the density and temperature evolution of the ejecta. [Abridged].

We quantify departures from hydrodynamical and centrifugal equilibrium in the gas discs of low-mass ($10^{10.75}<M_\mathrm{200c}/\mathrm{M}_\odot<10^{11}$) galaxies from the COLIBRE cosmological hydrodynamical simulations. We evaluate the full Eulerian acceleration balance in the midplane and show that disequilibrium is widespread: equilibrium-based circular velocity estimates typically have errors of $\geq 10$ per cent ($\approx 75$ per cent of midplane gas by mass). Disequilibrium is strongest and the largest associated errors occur in the inner few kiloparsecs that are crucial for constraining the dark matter density profile. Correcting the circular velocity to account for pressure and convective terms does not reliably improve its recovery in strongly perturbed systems where time-dependent forces dominate the residual acceleration budget. Stellar feedback, self-gravitating gas clumps and AGN energy injection account for most strong local perturbations, and large-scale gravitational asymmetries act as a scaffold for disequilibrium. We classify gas discs into coherent, perturbed, and slow/erratic rotators and show that this classification correlates with galaxy properties like mass, morphology and tracers of recent feedback. A majority of galaxies in our sample would be unsuitable for standard rotation curve analysis. Much of the observed diversity in the shapes of dwarf galaxy rotation curves may stem from non-equilibrium gas motions rather than diversity in mass profiles - resolving the discrepancy is then first and foremost a problem in gas dynamics.

We study the prevalence of thin galaxies as a function of stellar mass in the range $10^7<M_{\star}/\rm{M_\odot}<10^{11}$ using data from the GAMA, DESI, ALFALFA and Nearby Galaxy catalogs. We use the distribution of projected axis ratios, $q$, to infer the abundance of intrinsically flat galaxies needed to reproduce the observed abundance of highly elongated systems in projection. We find that as many as $40\%$ of galaxies in the mass range $10^9<M_{\star}/\rm{M_\odot}<10^{10}$ are intrinsically flatter than $1$:$5$ (i.e., $c/a<0.2$), a fraction that rises to $\sim 80\%$ for $c/a<0.3$. Although the incidence of thin galaxies decreases towards lower and higher $M_{\star}$, they are still quite common in dwarfs: $\sim 30\%$ and $\sim 65\%$ of $\sim 10^8 ~ \rm{M_\odot}$ galaxies are inferred to be intrinsically flatter than $c/a=0.2$ and $0.3$, respectively. A comparison of these results with several state-of-the-art cosmological hydrodynamical simulations (TNG50, FIREbox, Romulus25) reveals a distinctive lack of thin simulated dwarfs. In particular, there are no $M_{\star} < 10^9 ~ \rm{M_{\odot}}$ simulated galaxies flatter than $c/a=0.2$, in clear contrast with observational samples. This discrepancy likely reflects limitations in resolution and in the treatment of baryonic physics, suggesting that our understanding of the mechanisms regulating the formation of disk galaxies less massive than the Milky Way is still quite incomplete. Our results present a clear challenge to current numerical models of dwarf galaxy formation, which future models should attempt to meet.

In this study, we investigate the dynamics of Intermediate-Mass Black Hole (IMBH) binaries within Nuclear Star Clusters (NSCs) that contain a population of stellar-mass black holes (BHs). We examine how these stellar and BH populations influence the dynamics of the IMBH binary and, in turn, how the evolving IMBH binary affects the surrounding stellar and BH populations. We conduct high-resolution $N$-body simulations of NSCs constructed based on observational parameters from two local dwarf galaxies: NGC205 and NGC404. For the first time, we achieve a star particle mass resolution of $1\rm\;M_{\odot}$ and a BH mass resolution of $10\rm\;M_{\odot}$. This level of resolution is crucial for accurately modeling the collisional dynamics of these dense systems. Including stellar-mass BHs within the stellar population significantly influences the IMBH binary dynamics, nearly doubling the sinking rate and halving the merger time. During the initial phase of the inspiral, the IMBH binary disrupts both the stellar and BH cusps. However, the BH cusp quickly regains its steep slope due to its shorter relaxation time and continues to dominate the evolution of the IMBH binary, despite being much less massive compared to the stellar component. We uncover an interesting mechanism in which BHs first efficiently extract energy from the IMBH binary and then transfer this energy to the surrounding stars, allowing the BHs to spiral back toward the center of the system and restart the process. Our results imply that, although stellar mass BHs are a minor component of a stellar population, they can significantly facilitate IMBH growth within NSCs via mergers. We also notice that these dense systems can potentially boost Intermediate Mass Ratio Inspirals (IMRIs) predominantly on radial orbits.

We present constraints on modified gravity from a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy) in combination with DESI measurements of baryon acoustic oscillations, eBOSS observations of redshift space distortions, and cosmic microwave background anisotropies from Planck. We study the Horndeski class of modified gravity models in an effective field theory framework employing a parameterisation that satisfies stability conditions by construction and, for the first time, present a cosmological analysis in this inherently stable parameter basis. Cosmic shear constrains the Horndeski parameter space significantly, matching or surpassing the CMB contribution. Adopting the de-mixed kinetic term of the scalar field perturbation, $D_{\rm kin}$, and the deviation of the Planck mass from its fiducial value, $\Delta M_*^2\equiv M_*^2-1$, as model parameters, we constrain their present values to be $\Delta \hat{M}_*^2=0.32^{+0.07}_{-0.21}$ and $\hat{D}_{\rm kin} = 3.74^{+0.69}_{-1.92}$, which deviate from general relativity at $1.5\sigma$ and $1.9\sigma$, respectively. We derive constraints on the structure growth parameter $S_8=0.813^{+0.008}_{-0.011}$, which is compatible with the $\Lambda$CDM constraint at $0.54\sigma$. We obtain the deviation of the effective Newtonian coupling from the GR value as $\Delta \mu_{\infty,{\rm eff}}=0.066\pm0.023$, corresponding to a $2.9\sigma$ significance. Although modified gravity provides a slightly better fit to the data, a model comparison shows only a weak preference for modified gravity at the $1.4\sigma$ level. When adopting a dynamical dark energy model of the background cosmology, the inferred modified gravity parameter constraints are stable with respect to a $\Lambda$CDM background, while a mild preference at $1.57\sigma$ for dynamical dark energy remains.

We constrain minimally extended cosmological models with the cosmic shear analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy) in combination with external probes. Due to the consistency of the KiDS-Legacy analysis with the Cosmic Microwave Background (CMB), we can combine these data sets reliably for the first time. Additionally, we use CMB lensing, galaxy redshift-space distortions, and baryon acoustic oscillations. We assess, in turn, the effects of spatial curvature, varying neutrino masses, and an evolving dark energy component on cosmological constraints from KiDS-Legacy alone and in combination with external probes. We find KiDS-Legacy to be consistent with the fiducial flat $\Lambda$CDM analysis with $c^2 \sum m_\nu\leq 1.5\,$eV, $w_0 = -1.0\pm 0.7$ and $w_a = -1.3^{+1.9}_{-2.0}$ while $\Omega_K = 0.08^{+0.16}_{-0.17}$ (1$\sigma$ bounds) with almost equal goodness-of-fit. $w_0w_a$CDM is not a significant improvement over $\Lambda$CDM when cosmic shear and CMB lensing are combined, yielding a Bayes factor $B = 0.07$. If all probes are combined, however, $B$ increases to 2.73, corresponding to a $2.6\sigma$ suspicousness tension. The constraint on $S_8 = \sigma_8\sqrt{\Omega_\mathrm{m}/0.3}$ is robust to opening up the parameter space for cosmic shear. Adding all external datasets to KiDS-Legacy, we find with $S_8 = 0.816 \pm 0.006$ in $\Lambda$CDM and $S_8 = 0.837 \pm 0.008$ in $w_0 w_a$CDM for all probes combined.

We present new JWST IFS observations of the active galaxy NGC 1068, combining Mid-IR and optical IFS data from MIRI and MUSE to characterize the multi-phase circumnuclear gas properties and its interaction with the AGN outflow and radio jet. MIRI data trace the multiphase gas emission up to 400 pc from the nucleus at 20--60 pc resolution, unveiling a clumpy ionized structure around the radio hot-spots and a rotating warm molecular disc. Innovative Mid-IR diagnostic diagrams highlight the role of the AGN as the main excitation source for the ionized gas in the entire MIRI field of view, consistent with optical diagnostics, and supporting the AGN-driven wind scenario. Density sensitive [NeV] and [ArV] Mid-IR transitions reveal high-density clumps (n_e > 10**4 cm**-3) along the edges of the jet and outflow, tracing gas compression by the expanding wind. We combined multi-cloud kinematic (MOKA) and photo-ionization (HOMERUN) modeling to characterize the ionized outflow properties and found that [OIV] traces an outflow 300 km/s faster than that inferred from [OIII], showing that the two lines originate from distinct gas components. This kinematic dichotomy is confirmed by the photoionization analysis, which requires a dust-poor component dominating the optical lines and a dust-rich component responsible for the Mid-IR emission. The Mid-IR-revealed dusty component carries a significantly larger ionized-gas mass than what can be inferred from optical lines alone, showing that most of the outflowing mass is hidden from classical optical diagnostics. Our modelling point to a two-stage acceleration scenario, with velocities up to ~2000 km/s, consistent with an energy-driven wind. Our findings indicates that the outflow entrains up to a few 10**6 solar masses of ionized gas and couples efficiently with the surrounding ISM, injecting turbulence and impacting the host-galaxy environment.

Data-driven stellar classification has a long and important history in astronomy, dating as far back as Annie Jump Cannon's "by eye" classifications of stars into spectral types still used today. In recent years, data-driven spectroscopy has proven to be an effective means of deriving stellar properties for large samples of stars, sidestepping issues with computational efficiency, incomplete line lists, and radiative transfer calculations associated with physical stellar models. A logical application of these algorithms is the detection of unresolved stellar binaries, which requires accurate spectroscopic models to resolve flux contributions from a fainter secondary star in the spectrum. Here we use The Cannon to train a data-driven model on spectra from the Keck High Resolution Echelle Spectrometer. We show that our model is competitive with existing data-driven models in its ability to predict stellar properties Teff, stellar radius, [Fe/H], vsin(i), and instrumental PSF, particularly when we apply a novel wavelet-based processing step to spectra before training. We find that even with accurate estimates of star properties, our model's ability to detect unresolved binaries is limited by its approx. 3% accuracy in per-pixel flux predictions, illuminating possible limitations of data-driven model applications.

The center of the Galaxy harbors a supermassive black hole, Sgr\,A*, which is surrounded by a massive star cluster known as the S-cluster. The most extensively studied star in this cluster is the B-type main-sequence S2 star (also known as S0-2). These types of stars are commonly found in binary systems in the Galactic field, but observations do not seem to detect a companion to S2. This absence may be attributed to observational biases or to a dynamically hostile environment caused by phenomena such as tidal disruption or mergers. Using a $N$-body code with first-order post-Newtonian corrections, we investigate whether S2 can host a stellar or planetary companion. We perform $10^{5}$ simulations adopting uniform distributions for the orbital elements of the companion. Our results show that companions may exist for orbital periods shorter than 100~days, eccentricities below 0.8, and across the full range of mutual inclinations. The number of surviving companions increases with shorter orbital periods, lower eccentricities, and nearly coplanar orbits. We also find that the disruption mechanisms include mergers driven by Lidov-Kozai cycles and breakups that occur when the companion surpasses the Hill radius of its orbit. Finally, we find that the presence of a companion would alter S2's astrometric signal by no more than $ 5\,\mu{\rm as}$. Current radial-velocity detection limits constrain viable stellar binary configurations to approximately 4.4\% of the simulated cases. Including astrometric limits reduces to 4.3\%. Imposing an additional constraint that any companion must have a mass $\lesssim 2\,M_{\odot}$ (otherwise it would be visible) narrows the fraction of undetectable stellar binaries to just 3.0\%.

Little red dots (LRDs), a population of active galactic nuclei (AGNs) recently discovered by JWST, show distinctive Balmer-transition features, including prominent Balmer absorption, pronounced Balmer breaks, and large equivalent widths of broad $\mathrm{H}\alpha$ emission, all of which indicate the presence of dense gas surrounding their central black holes. A further key property of LRDs is their large Balmer decrements with broad $\mathrm{H}\alpha/\mathrm{H}\beta$ line-flux ratios far exceeding the Case B recombination value. These ratios of $\mathrm{H}\alpha/\mathrm{H}\beta>3$ have often been interpreted as evidence for heavy dust extinction ($A_V\gtrsim 3$ mag), however such dust would inevitably produce strong near-to-mid infrared re-emission that is hardly seen in JWST/MIRI observations. To investigate the physical origin of these observed Balmer features, we perform radiation transfer calculations through dust-free, dense gas. We show that the observed large Balmer decrements ($\mathrm{H}\alpha/\mathrm{H}\beta$ and $\mathrm{H}\alpha/\mathrm{H}\gamma$) naturally arise from Balmer resonance scattering without invoking dust. At sufficiently high densities ($n_\mathrm{H} \gtrsim 10^{{8}-{10}}~\mathrm{cm^{-3}}$), the elevated multiple Balmer-line ratios converge to values that closely mimic dust reddening, explaining why LRD spectra resemble obscured AGNs. Furthermore, when the Balmer break and broad Balmer lines originate in the same dense gas, their strengths are physically linked, allowing us to constrain the density structure and infer a low broad-line region gas mass of $\sim O(10~M_\odot)$. Such a small gas reservoir would be enriched by even a single supernova, implying that LRDs with observed low-metallicity signatures likely experienced minimal star formation in their nuclei.

This paper reports recent discoveries from the Globular Clusters GMRT Pulsar Search (GCGPS) survey, which aims to uncover pulsars in the globular clusters (GCs) of the Milky Way using the upgraded Giant Metrewave Radio Telescope (uGMRT). Utilising the Band-4 (550$-$750 MHz) and Band-3 (300$-$500 MHz) receivers, the survey targets GCs accessible to uGMRT ($-53^\circ\,<\,\delta\,<\,-17^\circ$), excluding the declination range that can be covered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The survey focuses on GCs that have not previously been searched with comparable sensitivity in these radio frequencies. In this paper, we present the discovery of the five MSPs in two GCs, $-$ NGC~6637 (M69) and NGC~6681 (M70), each hosting MSPs identified here for the first time. Observations of M69 led to the discovery of two MSPs: J1831$-$3220A (M69A) and J1831$-$3220B (M69B), both of which we localize with arcsecond precision using interferometric imaging. Observations of M70 resulted in three new MSPs: J1843$-$3217A (M70A), J1843$-$3217B (M70B), and J1843$-$3217C (M70C). Although direct imaging did not yield precise localizations for these MSPs, we provide initial estimates based on uGMRT beam forming and imaging analysis. Additionally, we present preliminary imaging results for other observed GCs, and in cases of non-detections, we report upper limits on pulsed emission based on the rms noise levels in the image plane.

We present a novel self-consistent theoretical framework to characterize the formation, evolution, and merger sites of dynamically-formed black hole binaries, with a focus on explaining the most massive events observed by the LIGO-Virgo-KAGRA Collaboration. Our approach couples the galaxy formation model GAMESH with cluster population synthesis codes to trace the cosmic evolution of globular clusters simultaneously with mergers of massive black holes. Our reference model, which includes prescriptions for both cluster formation and disruption depending on properties of specific galaxies, accurately reproduces the observed age-mass distribution of the Milky Way globular clusters. We find that approximately 30% of the globular clusters observed in our galaxy's halo may have originated from satellite galaxies of the Milky Way. We confirm that hierarchical black hole mergers provide a significant contribution to the formation of black holes in and above the pair-instability mass gap. However, quantifying their contribution is challenging, as different population synthesis codes yield divergent results in terms of black hole mass function and merger rates. Furthermore, we characterize the host galaxies where massive black holes form in terms of their dark matter, stellar mass, and metallicity. Ultimately, we demonstrate that the merger and birth rate densities of binary black holes increase with redshift till z = 5. This cosmic evolution is a crucial signature with significant implications for future detectors like the LISA, the Einstein Telescope and Cosmic Explorer, which will be capable to probe the high-redshift Universe.

Active galactic nuclei (AGN) provide energetic feedback necessary to `turn off' star formation in high-mass galaxies (M$_{\rm halo} \geq $ 10$^{12.5}$ M$_{\odot}$, $10.4 \leq \log(\frac{M_*}{M_\odot}) \leq 11$) as observed. Cosmic rays (CRs) have been proposed as a promising channel of AGN feedback, but the nature of CR feedback from AGN remains uncertain. We analyze a set of high-resolution simulations of massive galaxies from the Feedback in Realistic Environments (FIRE-3) project including multi-channel AGN feedback, explicitly evolving kinetic/mechanical, radiative, and spectrally-resolved CRs from the central black hole. Specifically, we explore different CR feedback and transport assumptions, calibrated to Milky Way local ISM constraints, and compare them to observed galaxy scaling relations. We find that all parameterizations explored self-regulate within agreement with observed galaxy scaling relations, demonstrating that CR injection efficiencies varied by $\sim$1.5 dex and locally-variable transport produce quenched galaxies with reasonable bulk properties; however, they feature orders-of-magnitude variant circumgalactic medium (CGM) gas properties. Our results indicate that multi-wavelength synthetic observations probing these varied halo properties from larger simulated samples in conjunction with observational comparisons may place novel constraints on how AGN physically quench star formation in massive galaxies.

Cislunar space spans from geosynchronous altitudes to beyond the Moon and will underpin future exploration, science, and security operations. We describe and release an open dataset of one million numerically propagated cislunar trajectories generated with the open-source Space Situational Awareness Python package (SSAPy). The model includes high-degree Earth/Moon gravity, solar gravity, and Earth/Sun radiation pressure; other planetary gravities are omitted by design for computational efficiency. Initial conditions uniformly sample commonly used osculating-element ranges, and each trajectory is propagated for up to six years under a single, fixed start epoch. The dataset is intended as a reusable benchmark for method development (e.g., space domain awareness, navigation, and machine learning pipelines), a reference library for statistical studies of orbit families, and a starting point for community-driven extensions (e.g., alternative epochs). We report empirically observed stability trends (e.g., a band near 5 GEO and persistence of some co-orbital classes including L4/L5 librators) as dataset descriptors rather than new dynamical results. The chief contribution is the scale, fidelity, organization (CSV/HDF5 with full state time series and metadata), and open availability, which together lower the barrier to comparative and data-driven studies in the cislunar regime.

We present a systematic study of 189 dwarf galaxies and their globular cluster (GC) systems in the Perseus cluster, based on deep Subaru Hyper Suprime-Cam imaging and Keck spectroscopy, supplemented by literature data. This constitutes the largest sample of dwarfs in a single galaxy cluster to date with simultaneous deep imaging, spectroscopic coverage, and GC measurements, while uniquely spanning a broad and continuous range of galaxy properties. We find an anti-correlation between GC specific mass and galaxy stellar mass for dwarfs in Perseus similar to observations in other clusters. At fixed stellar mass, dwarfs with lower surface brightness or larger effective radius tend to be more GC-rich -- suggesting either high GC formation efficiency in an earlier compact-galaxy phase, or less efficient GC disruption. The correlation between GC richness and axis ratio in Perseus is weaker than in other environments. We find some connection between GC richness and infall time, but not with the clear correlations found in Virgo, Coma, and cosmological simulations. More complete observations are needed to test for cluster-to-cluster variations in galaxy and GC evolutionary histories. This work demonstrates the potential of new wide-field imaging and spectroscopy surveys for understanding GCs and dwarf galaxies, and highlights the need for further work in theoretical modeling.

A cosmological origin of the magnetic fields in large scale structures of the Universe would require a non-negligible magnetic field in cosmic voids, which, however, remains undetected. Gamma-ray emission from gamma-ray bursts (GRBs) offers the opportunity to indirectly probe such an intergalactic magnetic field (IGMF), as gamma rays interact with cosmic radiation fields, producing electron-positron pairs, and initiate an electromagnetic cascade. The deflection of the pairs in the IGMF results in a time-delayed signal at GeV energies. Using observations with the Fermi Large Area Telescope of the GRB 221009A, we are able to derive the most stringent constraints to date from the non-observation of the cascade and rule out magnetic fields B < 2.5 x 10^{-17} G at 95% confidence level for a coherence length larger than 1 Mpc. Our results are comparable to limits obtained from blazar observations but do not suffer from assumptions on the duty cycle of the gamma-ray source or whether inverse-Compton scattering losses dominate over the development of plasma instabilities.

In the exascale computing era, handling and analyzing massive datasets have become extremely challenging. In situ analysis, which processes data during simulation runtime and bypasses costly intermediate I/O steps, offers a promising solution. We present libyt (this https URL), an open-source C library that enables astrophysical simulations to analyze and visualize data in parallel computation with yt or other Python packages. libyt can invoke Python routines automatically or provide interactive entry points via a Python prompt or a Jupyter Notebook. It requires minimal intervention in researchers' workflow, allowing users to reuse job submission scripts and Python routines. We describe libyt's architecture for parallel computing in high-performance computing environments, including its bidirectional connection between simulation codes and Python, and its integration into the Jupyter ecosystem. We detail its methods for reading AMR simulations and handling in-memory data with minimal overhead, and procedures for yielding data when requested by Python. We describe how libyt maps simulation data to yt frontends, allowing post-processing scripts to be converted into in situ analysis with just two lines of change. We document libyt's API and demonstrate its integration into two astrophysical simulation codes, GAMER and Enzo, using examples including core-collapse supernovae, isolated dwarf galaxies, fuzzy dark matter, the Sod shock tube test, Kelvin-Helmholtz instability, and the AGORA galaxy simulation. Finally, we discuss libyt's performance, limitations related to data redistribution, extensibility, architecture, and comparisons with traditional post-processing approaches.

Occultations, the covering up of one celestial body by another celestial body, have been used in astronomy for millennia to learn about the sun and moon. Since 2018, VERITAS has implemented a program to detect predicted asteroid occultations, where an asteroid covers up a star. VERITAS has attempted to observe over 100 occultations to date and successfully observed 20 occultations. With these occultations, VERITAS can directly measure the smallest angular diameters of any instrument or technique in the optical for stars between magnitude 9 and 13. Each angular diameter is measured by fitting the diffraction pattern observed by the central VERITAS pixel at the start and end of an occultation. Once a planned FADC upgrade is complete, VERITAS will begin a program to search for serendipitous occultations within its full field of view (3 deg). Serendipitous occultations of sub-km trans-Neptunian objects (TNOs) have the potential to constrain models of solar system formation. This work details how VERITAS predicts and observes occultations as well as the overall status of the asteroid occultation program and future steps for observing occultations of both asteroids and TNOs.

Context. Sunspot boundaries are commonly outlined by contours of the continuum intensity. However, their magnetic nature has not yet been fully characterised. Aims. We investigate the properties of the outer boundary of a long-lived sunspot to identify the magnetic property that defines it. Methods. We analysed the magnetic properties of AR NOAA11591 spot during its two passages across the solar disc, using SDO/HMI continuum intensity and magnetic field, and determined their contours to outline the outer boundary. Results. During the 1st disc passage, in which the sunspot is in its stable phase, the intensity contours at 0.9 of the mean quiet Sun intensity and isocontours of the magnetic field strength of 625G provide an almost perfect match between the two contours. With these thresholds, the time-averaged area of mismatch is minimised, yielding an average distance between the contours of 0.58 pixel, corresponding to less than 0.26 arcsec. During the 2nd disc passage, the spot shows clear signs of decay, and we find that the 0.9 intensity and 625G magnetic isocontours detach from each other, coupled to the disappearance of penumbra. In this super-equipartition area, granulation still operates. Conclusions. Based on a comparison with simulation data from our previous work, and in agreement with findings of other authors, we conclude that the outer boundary of stable sunspots is defined by an invariant magnetic field: the equipartition field. From the discrepancy between intensity and magnetic contours during the decaying phase of the sunspot, we surmise that alongside the well-established (magneto-)convective regimes of the photosphere - granular, penumbral, and umbral - a super-equipartition granular regime can be identified. In this regime, bright, but smaller granules occur where the magnetic field exceeds equipartition but remains sub-critical for convection suppression.

Machine learning is a field that has been growing in importance since the early 2010s due to the increasing accuracy of classification models and hardware advances that have enabled faster training on large datasets. In the field of astronomy, tree-based models and simple neural networks have recently garnered attention as a means of classifying celestial objects based on photometric data. We apply common tree-based models to assess performance of these models for discriminating objects with similar photometric signals, pulsars and black holes. We also train a RNN on a downsampled and normalized version of the raw signal data to examine its potential as a model capable of object discrimination and classification in real-time.

Big Bang Nucleosynthesis (BBN), as one of the earliest processes in the universe accessible to direct observation, offers a powerful and independent probe of the cosmic expansion history. With recent advances in both theory and observation, including efficient and flexible BBN codes, percent-level measurements of primordial deuterium and helium-4 abundances, refined measurements of nuclear reaction rates, and precise determinations of the baryon density from the cosmic microwave background, particularly keen insights can be gained from BBN. In this work, we leverage these developments to place model-independent constraints on deviations from the Standard Model expansion history during BBN. Using the latest abundance data, we apply principal component analysis to identify the most constrained and physically meaningful modes of expansion history variation. This approach allows us to impose the most general constraints on early dark energy during the epoch of BBN. We further examine whether general modifications to the expansion rate could alleviate the long-standing lithium problem. Our results demonstrate that BBN, sharpened by modern data and statistical techniques, remains an indispensable probe of dark energy and new physics in the early universe.

Nuclear star clusters (NSCs) surrounding supermassive black holes (SMBHs) are among the densest stellar environments in the universe. In these environments, collisions can shape the stellar mass function and produce exotic stellar populations. In this work, we investigate how stellar collisions couple with stellar evolution in the inner parsec of an NSC. We simulate the evolution of a sample of $1000$ $1$ $M_\odot$ stars embedded in a uniform cluster of dynamically relaxed 0.5 $M_\odot$ stars. Using COSMIC to evolve stellar properties in time , we track the mass, radius, and evolutionary state of the stars as they collide in the cluster. Our results show that most stars within $0.1$~pc of the SMBH experience a collision while on the main-sequence. However, outside of this distance, stars collide during the red giant phase, when the stellar radius increases dramatically. We find that the most common type of collision -- main-sequence or red giant -- over the lifetime of the cluster depends on the steepness of the stellar cusp, which determines the spatial distribution of the stars in the cluster. These results show that stellar evolution plays a fundamental role in shaping the collisional history of stars in nuclear star clusters. Lastly, we consider whether the closest known stars to the Milky Way's SMBH have experienced a collision. We estimate that several of the S-stars have a high probability of experiencing a collision over their main-sequence lifetime, perhaps with implications for their observed youth and properties.

We review the formation of large-scale heterogeneity in the solid Earth from the magma ocean phase to the present day, focusing on lower-mantle structure and evolution, as well as continental formation and its impact on interior-exterior volatile exchange. The solidification of the magma ocean sets the "initial condition" for solid-state mantle structure and evolution, partitioning volatiles between the surface and interior, and differentiating major and trace elements between mantle reservoirs. A possibly long-lived basal magma ocean is a potential reservoir for iron and incompatible trace elements with important implications for the present-day structure of the lowermost mantle and the distribution of heat-producing elements. With emergent plate tectonics and mantle convection, the production of oceanic crust is the dominant differentiation mechanism, and may contribute to lower-mantle heterogeneity. In turn, convection acts to re-mix heterogeneity. Thus, the structure of the lower mantle records the early differentiation of our planet, as well as long-term evolution. The formation of continental crust and lithosphere has also an important influence on deep-mantle structure and composition, because the distribution of continents can control the scale of plate tectonics, and thus affect mantle convection patterns. While the formation of continental crust and lithosphere, scenarios can be tested by self-consistent long-term evolution models. The exchange of water between the surface and mantle also plays an important role in the dynamics of Earth\' s deep interior. By using relative sea-level change as a proxy for the interior-exterior exchange of water, it is possible to constrain the deep-mantle water cycle. Based on a comprehensive review of these topics, we conclude with several open questions that should be addressed by future studies.

Parker Solar Probe (PSP) observes abundant circularly polarized ion-scale waves throughout the inner heliosphere. These waves are a signature of the interplay between plasma microinstabilities and turbulent dissipation. We perform a mission-wide statistical survey of ion-scale waves observed by PSP, investigating if the waves correspond to specific free energy sources in the measured proton velocity distributions. We find that left-handed waves (LHWs) are frequently observed, with the fraction of time they are observed increasing closer to the Sun, reaching $\sim$30\%. Right-handed waves (RHWs) are less frequently observed, with the associated time fraction decreasing closer to the Sun. The observed LHWs are generally consistent with parallel propagating ion cyclotron wave (ICW) storms that occur continuously for extended periods of time. Turbulent energy spectra are consistently steeper when LHW storms are observed; these wave storms mediate the spatial transport of the free energy associated with temperature anisotropy. The observed RHWs are generally consistent with oblique and parallel fast magnetosonic waves (FMWs), and their observation is well correlated with enhanced proton parallel heat flux, which quantifies the presence of secondary proton populations. Using observations and the SAVIC machine learning instability identification algorithm, we identify a threshold on the proton heat flux beyond which FMWs are likely to be driven unstable by the proton beams. We are thus able to associate trends in the observed ion-scale waves with known sources of free energy for Encounters 3 through 24 of the PSP's prime science phase.

Telescope bibliographies record the pulse of astronomy research by capturing publication statistics and citation metrics for telescope facilities. Robust and scalable bibliographies ensure that we can measure the scientific impact of our facilities and archives. However, the growing rate of publications threatens to outpace our ability to manually label astronomical literature. We therefore present the Automated Mission Classifier (amc), a tool that uses large language models (LLMs) to identify and categorize telescope references by processing large quantities of paper text. A modified version of amc performs well on the TRACS Kaggle challenge, achieving a macro $F_1$ score of 0.84 on the held-out test set. amc is valuable for other telescopes beyond TRACS; we developed the initial software for identifying papers that featured scientific results by NASA missions. Additionally, we investigate how amc can also be used to interrogate historical datasets and surface potential label errors. Our work demonstrates that LLM-based applications offer powerful and scalable assistance for library sciences.

Current wide-field surveys discover ~15 Jupiter-family comets (JFCs) each year, typically identified via visual detection of a dust coma or tail. The same surveys also discover many asteroids that have distant JFC-like orbits, but with no reported activity. We observed asteroids on Jupiter-crossing orbits beyond the depth of typical survey imaging using the 2.5 m Isaac Newton Telescope. We used deep imaging to observe 16 asteroids in this region, plus 7 known comets for comparison. Three asteroids (2011 WM183, (669525) 2012 XO144, and 2020 RX133) showed surface brightness profiles consistent with low-level activity, equating to ~19% of our total sample. We note that 2020 RX133 is a Jupiter Trojan. When we considered the heliocentric distance range of the asteroids at the time they showed activity, this fraction increased to 33% of the targets in the 3.16 \leq Rh \leq 4.56 au region, and therefore it is possible to infer that at least ~30 asteroids with Tj \leq 3.05 and in the 4.05 < a < 5.05 au parameter space may potentially exhibit low-level activity. We also estimated nuclear radii for the three active targets of r_n = 1.8 \pm 0.2 km, r_n \leq 0.8 km, and r_n \leq 0.5 km for (669525) 2012 XO144, 2011 WM183, and 2020 RX133 respectively. The median color index for the observed asteroids is (g-r)_{PS1} = 0.52 \pm 0.13, aligning with those expected for D-type asteroids.

We investigated the morphology and strength of magnetic fields in pillar-shaped structures at the boundaries of HII regions by combining three-dimensional radiation-magnetohydrodynamic, R-MHD, simulations with synthetic polarimetric observations. Our analysis focuses on the first pillar formed self-consistently in the simulation and is used as a proof of concept to test the applicability of the Davis-Chandrasekhar-Fermi (DCF) method under conditions dominated by external agents. The pillar arises as the ionization front compresses a dense clump, producing a magnetically aligned, elongated structure whose morphology and field configuration closely resemble those observed in real systems such as the pillars of M16. Synthetic dust-polarization maps at 850 $\mu$m reproduce the large-scale magnetic morphology of the simulated pillar, confirming that polarimetry is a reliable tracer of magnetic field geometry. However, magnetic field strengths inferred using DCF-type methods systematically overestimate intrinsic values from the simulations by factors $\gtrsim 2$. We attribute this discrepancy to the fact that field alignment and amplification are primarily driven by external gas pressure from the expanding HII region rather than internal turbulence, thus violating key assumptions of the DCF method. Our results highlight the need for caution when applying classical DCF-based analysis to pillars or other structures shaped by external compressions.

We investigate how rotating convection responds to the imposition of a latitudinally-varying heat flux at the base of the convective layer. This study is motivated by the solar near-surface shear layer, whose flows are thought to transition from a buoyancy-dominated regime near the photosphere to a rotation-dominated regime at depth. Here, we conduct a suite of spherical 3-D, nonlinear simulations of rotating convection that operate in either the buoyancy-dominated (high-Rossby-number, high-Ro) or rotation-dominated (low-Rossby-number, low-Ro) regime. At the base of each model convection zone, we impose a heat flux whose latitudinal variation is opposite to the variation that the system would ordinarily develop. In both the low- and high-Ro regimes, a strong thermal wind balance is sustained in the absence of forcing. With a larger flux variation, this balance becomes stronger at high latitudes and weaker at low latitudes. The resulting differential rotation weakens in response and, at sufficiently high forcing, its latitudinal variation reverses for both low- and high-Ro systems. At fixed forcing, there exists a Rossby number above which the convective flows efficiently mix heat laterally, and the imposed flux variation does not imprint to the surface. At sufficiently high-Ro, thermal wind balance is no longer satisfied. We discuss these results within the context of the Sun's near-surface region, which possesses a weakened differential rotation when compared to the deep convection, along with little-to-no variation of photospheric emissivity in latitude.

This work presents possible quasi-periodic oscillations in the Transiting Exoplanet Survey Satellite observations of blazars that are in the 157-month hard X-ray survey done by Swift's Burst Alert Telescope. We report observations from four sources, J1104.4+3812, J1654.0+3946, J0353.4-6830, and J1941.3-6216, that show at least 3$\sigma$ local significance in generalized Lomb-Scargle periodogram and weighted wavelet Z-transform methods. These results are also checked using a continuous autoregressive moving average analysis that also predicts the absent data tentatively using stochastic differential equations. Each of these four sources exhibits a nominal QPO signal frequency in the range of 0.5--1.1 d$^{-1}$, resulting in at least 5 putative cycles. However, when the number of frequencies examined and the number of sources examined are both taken into account, the global significances are reduced to $\approx2.2\sigma$ for J1654.0+3964 (Mrk 501), and to $\approx1.9\sigma$ for other sources. These QPOs are thought to arise from the processes within the relativistic jet. Plausible explanations include the kink instability, which arises due to current-driven instabilities in the plasma or the precession of substructures, or mini-jets, within the jets.

Cool Planet Imaging Coronagraph (CPI-C) on Chinese Space Station Survey Telescope (CSST) is proposed to direct image the cool planets around nearby solar-type stars (within 40 pc). The core scientific objective of CPI-C is to conduct high-contrast directly imaging surveys of exoplanets ranging in size from Neptune-like to Jupiter-like, located at separations of 0.5 to 5 AU from their host stars, and to perform systematic spectroscopic analysis of the detected planets through high-precision multi-band photometry. CPI-C employs a step-transmission apodization technique to suppress the diffraction noises from the telescope pupil and a precise phase correction technique to eliminate the speckle noises due to imperfections of the optical surfaces. The contrast requirement is better than $10^{-8}$ at an inner working angle (IWA) of $3-4\lambda/D$, in the visible wavelength from 600 nm to 900 nm. CPI-C will be the first space-based instrument capable of directly imaging the reflection light from the cool exoplanets in the visible wavelength enabling the measurement of key physical parameters such as the effective temperature, surface gravity, radius, mass, and other key parameters. The potential observation results will significantly contribute to further understand the formation and evolution mechanisms of planets, which will also lay a solid foundation for future confirmation of the Earth-twins in the next generation space flagship missions.

We investigate the correlation between star formation rate (SFR) surface density and gas surface density (known as the Kennicutt-Schmidt, KS, relation) at kiloparsec (kpc) scales across cosmic time ($0\le z \le 8$) using the COLIBRE state-of-the-art cosmological hydrodynamical simulations. These simulations feature on-the-fly non-equilibrium chemistry coupled to dust grain evolution and detailed radiative cooling down to $\approx 10$~K, enabling direct predictions for the atomic (HI) and molecular (H$_2$) KS relations. At $z\approx 0$, COLIBRE reproduces the observed (spatially-resolved) KS relations for HI and H$_2$, including the associated scatter, which we predict to be significantly correlated with stellar surface density, local specific SFR (sSFR), and gas metallicity. We show that the HI KS relation steepens for lower-mass galaxies, while the H$_2$ KS relation shifts to higher normalisation in galaxies with higher sSFRs. The H$_2$ depletion time decreases by a factor of $\approx 20$ from $z = 0$ to $z = 8$, primarily due to the decreasing gas-phase metallicity. This results in less H$_2$ and more HI being associated with a given SFR at higher redshift. We also find that galaxies with higher sSFRs have a larger molecular gas content and higher star formation efficiency per unit gas mass on kpc scales. The predicted evolution of the H$_2$ depletion time and its correlation with a galaxy's sSFR agree remarkably well with observations in a wide redshift range, $0\le z\le 5$.

The apparent anisotropies of galaxy clustering and 21-cm mapping in redshift space offer a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales through the Alcock-Paczynski (AP) effect and redshift-space distortions (RSD). Although improved theoretical models exist for anisotropic clustering, their applicability is limited by the non-perturbative smearing effect caused by the randomness of relative velocities. Here, we consider an alternative approach using the statistical power of cross-correlation between galaxy clustering and 21-cm line intensity mapping. Based on Fisher matrix analysis, fully incorporating nonlinear RSD, we estimate the benefit of combining both observables. We find that, for spectroscopy surveys like DESI combined with 21-cm line-intensity mapping surveys, the constraint on the growth of structure is improved by a factor of two relative to the galaxy auto-correlation, while the constraint on the cosmic expansion is only slightly improved. Crucially, such an observation can strongly constrain the neutral hydrogen (HI) content Omega_HI to a sub-percent level. This level of precision unlocks the potential of this method to probe post-reionization astrophysics with enhanced precision. It would far surpass existing constraints from stacked 21-cm emission (resulting in ~ 50%-100% uncertainties) and break the degeneracy between Omega_HI and the HI bias b_HI inherent in linear-regime power spectrum analysis. This cross-correlation approach effectively compensates for the loss of constraining power when using galaxy clustering alone.

Core-envelope coupling provides a reasonable explanation of the spin-down stalling of stars in open clusters, which was not predicted by classical gyrochronology. However, it remains an open question whether the coupling efficiency is constant or variable. M dwarfs, possessing thicker convective envelopes and thus longer coupling timescales than other late-type stars, are ideal objects for this investigation. In this letter, we present a new analysis using LAMOST and DESI spectra to construct the rotation--activity ($R_{\rm{HK}}^{'}$--Ro) relation for M dwarfs. The new relation consists of three distinct regimes of fast, intermediate, and slow rotation, closely matching the three sequences of gyrochronology, namely the ``Convective'' sequence, ``Gap'', and ``Interface'' sequence. Our study reveals, for the first time, a variable activity decay rate in the intermediate-rotation regime (i.e., the ``Gap'' region). This implies a varying core-envelope coupling efficiency, peaking towards the end of this region. It also coincides with the well-known stage of stalled stellar spin-down.

Tidal torque theory (TTT) predicts that galaxy angular momenta are imprinted by the early tidal field acting on their proto-structures, which are preserved through cosmic evolution and provide the potentially most precise measurement of the early universe. We test this prediction using the gas component of central massive elliptical galaxies, whose angular momenta respond sensitively to external torques. By comparing the observed gas angular momentum vectors with those predicted from the primordial density field reconstructed by ELUCID for the nearby universe, we conclusively detect a strong direction correlation with a significance of about $7\sigma$. These results provide the first robust observational confirmation of TTT and open a new window for cosmological measurements of neutrino mass and other cosmological parameters.

Large-scale intergalactic magnetic fields may contain a mixture of galactic and cosmogenic contributions, that can be probed via observations of $\gamma$-ray "echo" from sources at cosmological distances. While these contributions may be disentangled based on the difference in their redshift evolution, thus far indications of non-negligible magnetic field have been found only at low redshifts. This work aims at extending the intergalactic magnetic field constraints to redshifts $z\gtrsim 1$ using 17-year long all-sky observations of high-redshift active galactic nuclei with Fermi/LAT $\gamma$-ray telescope. Combing the Fermi/LAT measurements in the 0.1 GeV - 1 TeV energy range with the Monte Carlo simulations of the $\gamma$-ray "echo", it is shown that the zero field strength hypothesis at $z = [0.5; 3]$ is disfavoured at the $\approx 9\sigma$ significance level, yielding the lower limit of $B \gtrsim 1\times10^{-18}$ cG for the magnetic field correlation length above 1 Mpc. It is further shown that the same data put a lower limit on the volume-filling fraction of this field of $f\gtrsim 90\%$ in the redshift range considered. It is also demonstrated that the derived limits are not substantially affected by either source flux variability or the assumed $\gamma$-ray emission attenuation model. Implications of these limits for intergalactic magnetic field origin are discussed.

Fallback in core-collapse supernovae plays a central role in setting compact-remnant masses and may produce late-time emission. In hydrogen rich progenitors, the reverse shock arising at the hydrogen-helium interface has the potential to dramatically enhance fallback, yet its overall impact across a broad explosion-energy range has not been systematically quantified. Using one-dimensional hydrodynamic simulations for metal-poor progenitors with $M_{\rm ZAMS}=18$-$28\,M_\odot$ and models with and without hydrogen envelopes, we explore fallback over explosion energies of $10^{48}$-$10^{52}\,{\rm erg}$. We find a robust and universal mass-transition behaviour: when the explosion energy reaches only $2$-$3$ times the binding energy of the hydrogen envelope, the reverse shock returns to the centre and sharply increases the remnant mass by $\gtrsim 2\,M_\odot$. Above this threshold, the reverse shock escapes and hydrogen-rich and stripped-envelope progenitors yield nearly identical remnant masses. By normalizing the results with the envelope binding energy, we show that all progenitor models converge to a common fallback relation. We further provide a simple analytic prescription that connects explosion energy, hydrogen-envelope binding energy, and final compact-remnant mass. This relation provides an important link between progenitor properties and compact-remnant masses, and is useful for population-synthesis and galactic chemical-evolution studies.

We explore the roles of the three competitors, namely, gravity, turbulence, and magnetic fields, in controlling star formation (SF) within dense, massive clumps identified in the ATLASGAL survey. By examining scaling relations, virial parameters, and turbulent energy spectra, we evaluate the dynamical state of these clumps. We observe a weak velocity dispersion-size relation (sigma proportional to L^0.11), which is much shallower than the classical Larson-like relations, suggesting that turbulence does not mainly drive internal dynamics. The turbulent energy spectrum, E(k) proportional to k^-1.21, is also less steep than what is expected for both incompressible and compressible turbulence. We equally observe a decreasing trend in the virial parameter with increasing mass (alpha_vir proportional to M^-0.37), indicating that more massive clumps are increasingly gravitationally bound. These trends indicate an increasing relative dominance of gravity over turbulence at smaller scales, aligning with multiscale collapse scenarios; however, the absolute energy balance remains unquantifiable with the current data. Although magnetic fields are not directly measured, their potential influence is considered in the interpretation of pressure balance and dynamical support. Our findings imply that gravitational processes appear to primarily regulate the structure and evolution of massive clumps.

Nickel, one of the most stable elements alongside iron, is the most abundant heavy element beyond iron in cosmic rays. With DAMPE's excellent charge resolution and broad energy range, a high-precision energy spectrum provides valuable insights into the acceleration sources of heavy nuclei and their propagation through the interstellar medium. In this analysis, we report the direct measurement of cosmic-ray nickel spectrum from 10 GeV/n to 2 TeV/n with nine years of flight data. The nickel spectrum is consistent with a single power law with spectral index -2.60 +/- 0.03 from 40 GeV/n to 1 TeV/n. This work provides an accurate measurement of differential flux of nickel with kinetic energy extending to TeV/n for the first time.

Stripped-envelope supernovae (SESNe), including Type IIb, Ib, and Ic supernovae (SNe), originate from the explosions of massive stars whose outer envelopes have been largely removed during their lifetimes. The main stripping mechanism for the hydrogen (H) envelope in the progenitors of SESNe is often considered to be interaction with a binary companion, while that for the helium (He) layer is unclear. We conducted photometric and spectroscopic observations of the Type Ib SN 2021efd, which shows signs of interaction with H-free circumstellar material (CSM). Around 30 days after the r-band light curve (LC) peak, until at least ~ 770 days, its LCs display excessive luminosity compared to regular SESNe and at least three distinct peaks. The light curve evolution is similar to that of SN 2019tsf, whose previously unpublished spectrum at 400 days is also presented here. The nebular spectrum of SN 2021efd shows narrow emission lines (~ 1000 km/s) in various species. Our observations suggest that SN 2021efd is a Type Ib SN interacting with CSM with the following parameters: The estimated CSM mass, composition, and distribution are at least a few times 0.1 M_sun, H-free, and clumpy, respectively. Based on the estimated ejecta properties, we conclude that this event is a transitional SN whose progenitor was experiencing He-layer stripping at the epoch of the explosion, and was on the way to becoming a carbon-oxygen star from a He star. The estimated CSM properties suggest that the progenitor had some episodic mass ejections with a rate of ~ 0.001-0.01 M_sun/yr for the last decade and slightly smaller before this final phase at least from ~ 200 years before the explosion, for the assumed CSM velocity of 100 km/s. For the case of ~ 1000 km/s, the necessary mass-loss rate would be increased by a factor of ten, and the timescales decreased by a factor of ten.

The Sun periodically emits intense bursts of radio emission known as solar radio bursts (SRBs). These bursts can disrupt radio communications and be indicative of large solar events that can disrupt technological infrastructure on Earth and in space. The risks posed by these events highlight the need for automated SRB classification, providing the potential to improve event detection and real-time monitoring. This would advance the techniques used to study space weather and related phenomena. A dataset containing images of radio spectra was created using data recorded by the Compound Astronomical Low frequency Low cost Instrument for Spectroscopy and Transportable Observatory (e-Callisto) network. This dataset comprises three categories: empty spectrograms; spectrograms containing Type II SRBs; and spectrograms containing Type III SRBs. These images were used to fine-tune several popular pre-trained deep learning models for classifying Type II and Type III SRBs. The evaluated models included VGGnet-19, MobileNet, ResNet-152, DenseNet-201, and YOLOv8. Testing the models on the test set produced F1 scores ranging from 87\% to 92\%. YOLOv8 emerged as the best-performing model among them, demonstrating that using pre-trained models for event classification can provide an automated solution for SRB classification. This approach provides a practical solution to the limited number of data samples available for Type II SRBs.

Polarization is one of the fundamental natures of electromagnetic radiation. The detection of polarization or polarized photons from distant X-ray radiating systems (such as X-ray binaries (XBs), active galactic nuclei (AGN), pulsars, and stars) complements the timing, spectral, and imagining analysis to better understand the physical mechanisms taking place in these sources. Polarization has enhanced the understanding of the internal geometry of these systems and their vicinity. Polarized X-rays can be generated either directly through non-thermal physical processes in the presence of a magnetic field(B) or through the scattering of unpolarized thermal radiation within plasma structures such as an accretion disk. X-ray polarization can measure the two important independent parameters, the polarization degree (PD) and polarization angle (PA) of the X-ray photons. These parameters are crucial as they reveal the characteristics of particles in such a strong magnetic and gravitational field. In this chapter, we have discussed (i) the basic idea of polarization, (ii) some distant sources radiating polarized X-ray photons, (iii) missions dedicated to observing polarized X-ray photons, and (iv) recent breakthroughs and upcoming missions.

Understanding the presence and distribution of prebiotic precursors in the interstellar medium (ISM) is key to tracing the chemical origins of life. Among them, 4-oxobutanenitrile (\ch{HCOCH2CH2CN}) has been identified in laboratory simulations as a plausible intermediate in the formation of glutamic acid, a proteinogenic amino acid. Here, we report its gas-phase rotational spectrum, measured using two complementary techniques: chirped-pulse Fourier transform microwave (CP-FTMW) spectroscopy (2$-$18 GHz) and free-jet millimeter-wave (FJ$-$AMMW) absorption spectroscopy (59.6$-$80 GHz). Quantum chemical calculations revealed nine low-energy conformers, of which the TC conformer was assigned based on the measured spectra. The resulting spectroscopic parameters were used to search for the molecule in the ultradeep spectral survey of the G+0.693-0.027 molecular cloud, located in the Galactic Center. No signal attributable to 4$-$oxobutanenitrile was detected. A stringent upper limit to its column density was derived ($N<$ 4 $\times$10$^{12}$ cm$^{-2}$), corresponding to a molecular abundance of $<$ 2.9 $\times$10$^{-11}$ relative to H$_{2}$. This upper limit lies well below the observed abundances of simpler structurally related species containing $-$HCO and $-$CN groups, underscoring the challenge of detecting increasingly complex prebiotic molecules in the ISM and the need for future, more sensitive astronomical facilities.

Context. Comets, asteroids, planetesimals, free-floating planets and brown dwarfs, are continuously injected into the intra-cluster environment after expulsion from their host planetary systems or binary system. The dynamics of large populations of such free-floating comets (ffcs) in a star cluster environment is not yet fully understood. Aims. We investigate the dynamical evolution of comet populations in star clusters, and characterize the kinematics and ejection rates of ffc in a star cluster. Moreover, we determine whether a different initial energy distribution affects the mass segregation of the less massive population. Methods. We carry out simulations using the N-body code NBODY6++GPU-MASSLESS (Flammini Dotti et al. 2025), which allows fast integration of star clusters that contain large numbers of massless particles, to characterize the dynamics of populations of low-mass particles with sub-virial and super-virial distributions. Results. Comets do not participate in the mass segregation process, similarly to planet-size objects, regardless of their initial energy distribution. The latter is slightly changing the whole dynamical evolution at the start of the simulation. We only observe an initial relaxation or collapse of the objects for super-virial and sub-virial ratios, respectively. The external regions of the ffcs population tend to be pulled back in the cluster core at the end of the simulation, suggesting the gravitational pull of the stars is pulling them back in the core. This phenomenon occurs at later times if the system in virial equilibrium. Compared to less massive bodies, brown dwarfs experiences more mass segregation the inner regions tend to be more mixed with the stellar population

The dark ages 21-cm signal is a promising probe of the currently unobserved infant universe between the formation of the Cosmic Microwave Background around $z \approx 1100$ and the first galaxies around $z\approx 30$. A detection of the signal will help researchers understanding the nature of dark matter and dark energy, the expansion of the universe and any extensions to the concordance $\Lambda$CDM model that could explain the reported Cosmic Dawn 21-cm signal from EDGES and the Hubble tension. In this letter we take existing constraints on the $\Lambda$CDM cosmological model from two early time probes, Planck and WMAP, and two late time probes, DES galaxy lensing and clustering and Baryon Acoustic Oscillations, and propagate these through to constraints on the magnitude of the dark ages 21-cm signal. We constrain the magnitude and central frequency of the signal while methodically accounting for uncertainties in the cosmological parameters. We find that within the context of our modelling assumptions and the $\Lambda$CDM paradigm, the depth of the dark ages 21-cm signal is known to better than 1 mK and the central frequency to within 0.05 MHz.

We present James Webb Space Telescope (JWST) and Hubble Space Telescope (HST) observations of the counterpart of EP240801a, at z=1.67, the first fast X-ray transient (FXT) identified as an X-ray flash (XRF) by the Einstein Probe (EP) and Fermi-GBM. Our observations reveal strong photometric and spectroscopic evidence for an associated broad-lined Type Ic supernova (SN) SN 2024aihh, the most distant spectroscopically identified gamma ray-burst (GRB)-SN to date. The SN exhibits similar luminosity and light curve evolution to the proto-type GRB-SN 1998bw with an absolute magnitude of the SN at $\sim$23 days rest-frame of $M_{F140W} \approx$-19 mag. The SN is located in a host galaxy with complex morphology at a large ($\sim$6 kpc) offset in a region of relatively low surface brightness. The region around the SN has a modest star formation rate and is dominated by an intermediate mass-weighted age (1.4$\pm$0.3 Gyr) population, despite the apparent presence of a young, massive broad-lined Type Ic SN progenitor. These observations demonstrate that observations with HST and JWST can greatly extend the redshift range over which the GRB/FXT-SN connection can be studied, including in relatively low luminosity, X-ray rich events. They demonstrate little apparent evolution in the SN properties from local examples despite EP240801a originating from an epoch 10 Gyr ago.

Magnetic reconnection powers explosive releases of magnetic energy, heating and particle acceleration throughout the plasma universe. Knowledge of this universal process is vital to understanding the Heliosphere, as it plays a key role in solar flares, coronal mass ejections, coronal heating, solar wind acceleration, geomagnetic storms, and interactions between the solar wind and planetary magnetospheres. As such, reconnection underpins multiple science objectives of multiple future space missions. The UK plays a leading role in this international field, through a combination of in situ measurements from Earth's magnetosphere and the solar wind, observations of the solar corona and chromosphere, and world-class numerical simulations and theory. This white paper identifies: Nine priority science objectives for reconnection research in the next decade; Recommendations to guide investment in theory, simulations and infrastructure; Mission priorities and required measurements to ensure the UK maintains and improves its world-class credentials in reconnection science.

Models of interacting dark energy (DE) and dark matter (DM) involving pure momentum exchange are a promising avenue for resolving cosmological tensions. However, the behaviour of these interactions in the theoretically challenging limit where the DE equation of state, $w$, approaches $-1$ is not fully understood. We demonstrate that a generic feature of these models is a $w$-dependent velocity-locking mechanism, which systematically shifts the onset of matter power spectrum suppression to smaller scales as $w \rightarrow -1$. The suppression magnitude depends on the difference in fluid velocities. In this limit, however, the interaction's drag dominates over the DE pressure support and causes the DE velocity to track that of the DM fluid at larger scales. This mechanism provides a physical explanation for the weaker constraints found in the literature when $w\approx-1$ in models where the interaction strength does not explicitly depend on $w$. We also demonstrate that the common approximation of neglecting DE perturbations ($\delta_{\mathrm{DE}}=\theta_{\mathrm{DE}}=0$) fails in this limit. By artificially increasing the velocity difference between the fluids, this simplification incorrectly removes the $w$-dependent velocity-locking mechanism and erases the shift in power spectrum suppression to smaller scales. This leads to an overestimation of the constraining power of cosmological data on the interaction strength.

We present an original data analysis and a physical model that provide new insights into the origin of, and relationship between, two observables of the dusty, polarized Galaxy at intermediate and high latitudes: (i) the misalignment between HI filamentary structures and magnetic fields and (ii) the positive $TB$ correlation measured in Planck data suggesting parity violation in the interstellar medium (ISM). We confirm an observational link between the two effects and find that both are predominantly produced at large angular scales ($\geq 10^\circ$, multipoles $\ell \leq 20$) with a significantly stronger signal in the northern hemisphere. We propose a model in which filaments and magnetic fields appear misaligned in projection because they are sourced by cold and warm gas phases distributed in different proportions in the Solar neighborhood, from the wall of the Local Bubble to larger distances. These projection effects at large angular scales can produce coherent signatures that propagate to smaller scales in power spectra without invoking local, small-scale filament misalignment. Within this frame, HI filaments remain statistically aligned with the magnetic field in 3D, although with a projected scatter of tens of degrees that requires further investigation. The multiphase, geometrical model presented in this work is supported by Planck polarization data at 30 GHz, where synchrotron radiation dominates, and at 217 and 353 GHz, where dust emission dominates. Our analysis also incorporates starlight polarization measurements. The model introduced here suggests a new interpretation of two unexplained observables and emphasizes the role of the large-scale magnetized ISM in shaping polarized Galactic emission, which has important implications for both Galactic astrophysics and cosmological foreground characterization.

We present a cosmological zoom-in simulation targeting the high redshift compact progenitor phase of massive galaxies, with the most massive galaxy reaching a stellar mass of $M_{\star}=8.5\times 10^{10} \ M_{\odot}$ at $z=5$. The dynamics of supermassive black holes (SMBHs) is modelled from seeding down to their coalescence at sub-parsec scales due to gravitational wave (GW) emission by utilising a new version of the KETJU code, which combines regularised integration of sufficiently massive SMBHs with a dynamical friction subgrid model for lower-mass SMBHs. All nine massive galaxies included in this study go through a gas-dominated phase of early compaction in the redshift range of $z\sim 7-9$, starting at stellar masses of $M_\star\gtrsim 10^8\ \mathrm{M}_\odot$ and ending at a few times $M_{\star}\sim 10^9\ \mathrm{M}_\odot$. The sizes, masses and broad band fluxes of these compact systems are in general agreement with the population of systems observed with JWST known as `Little Red Dots'. In the compact phase, the stellar and SMBH masses grow rapidly, leading to a sharp decline in the central gas fractions. The outer regions, however, remain relatively gas-rich, leading to subsequent off-centre star formation and size growth. Due to the very high central stellar densities ($\rho_{\star}\gtrsim 10^{13}\,\mathrm{M_\odot/kpc^3}$), the SMBHs merge rapidly, typically just $\sim 4-35\ \mathrm{Myr}$ after the SMBH binaries have become bound. Combining KETJU with the phenomenological PhenomD model resolves the complete evolution of the GW emission from SMBH binaries through the Pulsar Timing Array frequency waveband up to the final few orbits that produce GWs observable with the future LISA mission.

Interacting Dark Energy (IDE) models offer a promising avenue to explore possible exchanges of energy and momentum between dark matter and dark energy, providing a dynamical extension of the standard $\Lambda$CDM paradigm. Such interactions modify the growth of cosmic structures, imprinting distinctive signatures on the matter power spectrum that can be tested through large-scale structure (LSS) observations. In this work, we compute the one-loop corrections to the matter power spectrum in IDE models. We then reinterpret these results within the standard framework of the Effective Field Theory of Large-Scale Structure (EFTofLSS), which provides a consistent description of mildly non-linear scales and allows for reliable comparisons with observational data. We investigate two commonly studied forms of the coupling function, $Q$, namely $Q = \xi \mathcal{H} \rho_{\rm m}$ and $Q = \xi \mathcal{H} \rho_{\rm DE}$, and introduce a novel interaction term, $Q = \Gamma \, \rho_{\rm m} \, \rho_{\rm DE} \, \theta_{\rm m}$, characterized by the non-linear coupling constant $\Gamma$, which links the interaction strength to the velocity divergence of dark matter. This coupling function is proposed to isolate the effects solely of the IDE model on mildly non-linear scales. Using Full-Shape (FS) measurements of the galaxy power spectrum from BOSS DR12, we constrain the interaction rate $\Gamma$, the cosmological parameters, and the bias parameters. We find $\Gamma = 0.0039 \pm 0.0082$, which is highly consistent with the $\Lambda$CDM model. This work opens the possibility of testing IDE models at mildly non-linear scales, potentially providing new insights for this class of models beyond the standard $\Lambda$CDM framework.

The search for free-floating planets (FFPs) is a key science driver for upcoming microlensing surveys like the Nancy Grace Roman Galactic Exoplanet Survey. These rogue worlds are typically detected via short-duration microlensing events, the characterization of which often requires analyzing noisy, irregularly-sampled observations. We present a pipeline for this task using simulation-based inference. We use a Transformer encoder to learn a compressed summary representation of the raw time-series data, which in turn conditions a neural posterior estimator. We demonstrate that our method produces accurate and well-calibrated posteriors over three orders of magnitude faster than traditional methods. We also demonstrate its performance on KMT-BLG-2019-2073, a short-duration FFP candidate event.

We study the spatially resolved star formation, gas ionisation, and outflow properties of 1813 active galactic nuclei (AGNs) from the MaNGA survey, which we classify into infrared (IR), broad-line (BL), narrow-line (NL), and radio (RD) AGNs based on their mid-infrared colours, optical spectra, and/or radio photometry. We also provide estimations of AGN power at different wavelengths. AGN incidence is found to increase with stellar mass following a power-law, with the high-mass end dominated by RDAGNs and the low-mass end dominated by NLAGNs. Compared to their mass-matched non-AGN counterparts, we find that IRAGNs, BLAGNs, and NLAGNs on average show enhanced specific star formation rates, younger stellar populations, and harder ionisation towards the centre. RDAGNs, in contrast, show radial profiles similar to quiescent galaxies. [OIII] outflows are more common and stronger in BL/IRAGNs, while RDAGNs on average show no outflow features. The outflow incidence increases with [OIII] luminosity, and the features in BL/IRAGNs on average extend to ~2 kpc from the nuclei. We further discuss a possible evolutionary sequence of AGNs and their host galaxies, where AGNs with strong emission lines or dust tori are present in star-forming galaxies. Later, young compact radio jets emerge, the host galaxies gradually quench, and the AGN hosts eventually evolve into globally quiescent systems with larger radio jets that prevent further gas cooling.

The multi-band HSC-CLAUDS survey comprises several sky regions with varying observing conditions, only one of which, the COSMOS Ultra Deep Field (UDF), offers extensive redshift coverage. We aim to exploit a complete sample of labeled galaxies from the COSMOS UDF at i<25 (z<~5) to train a convolutional neural network (CNN) and infer more accurate photometric redshifts in the other regions than those currently available from SED-fitting methods. To address the severe domain mismatch problem we observed when applying the trained CNN to regions other than the COSMOS UDF, we developed an unsupervised adversarial domain adaptation network that we grafted onto the CNN. The method is validated by three tests: the predicted redshifts are compared to the spectroscopic redshifts that are available for limited samples of mostly bright galaxies; the predicted redshift distributions of the entire galaxy population of a given field in several intervals of magnitude are compared to those of the COSMOS UDF, assumed to be representative; the redshifts predicted for a sample of galaxies selected by narrow-band filter observations sensitive to [OII] emitters at z~1.47 are compared to those of confirmed [OII] emission line galaxies. The results show successful domain adaptation: the network is able to transfer its redshift classification capability learnt from the COSMOS UDF to other regions of HSC-CLAUDS. Accuracy varies depending on magnitude and redshift, following that of the labels we used, but far exceeds that of currently available photometric redshifts. The catalogs of CNN redshifts we inferred for the XMM, DEEP2 and ELAIS fields and for the remaining COSMOS region (~4 million sources in total at i<25) are made public.

Context: This paper presents new results from the ESO Large Programme Looking into the faintEst WIth MUSE (LEWIS). The sample consists of low-surface brightness (LSB) and ultra-diffuse galaxies (UDGs) located inside 0.4 R$_{vir}$ of the Hydra I cluster. Integral field spectroscopy is acquired for 24 galaxies with the MUSE spectrograph mounted on the Very Large Telescope (VLT). Aims: Our main objective is to analyse possible correlations between the environment and the integrated stellar population properties. Methods: For each galaxy, we extract the 1D stacked spectrum in an aperture of one effective radius and adopt previously published stellar kinematics to derive age, metallicity and [Mg/Fe] through a full spectral fitting technique. Results: We find that the analysed LEWIS sample has a mean metallicity of [M/H] = -0.9 dex, a mean age of 10 Gyr, comparable to results of UDGs in other clusters. According to the projected phase-space, galaxies can be classified into two groups: early infallers, with slightly higher metallicities ([M/H]$_{early} = -0.8 \pm 0.1$ dex), and late infallers, with slightly lower values ([M/H]$_{late} = -1.0 \pm 0.1$ dex). Two types of galaxies are found in the early-infall region: half of them have metallicities consistent with the dwarf mass-metallicity relation, while the other half show higher values. Moreover, they display different timescales for stellar mass assembly. Conclusions: We identified different classes of UDGs in the Hydra I cluster, which suggest different formation mechanisms taking place. Almost all of the UDGs and LSBs in the cluster are consistent with the puffed-up dwarf formation scenario, having dwarf-like metallicities and being consistent with the dwarf mass-metallicity relation. In the innermost regions of the cluster, where metal-richer UDGs lie, tidal or environmental effects might have influenced their formation and evolution.

While convection has been known to play a key role in stars for many decades, its implementation in one-dimensional stellar evolution codes still represents a major uncertainty today. The purpose of this work is to investigate the impact of initial mass dependent convective boundary mixing (CBM), often referred to as overshooting, on the frequency and type of nuclear burning shell interactions that occur in low metallicity massive stars and the subsequent effect on their fates. Two grids of models were calculated using the Modules for Experiments in Stellar Astrophysics (MESA) code and a 22-isotope nuclear network, each with a different strength of CBM applied. One grid uses the typical CBM value for diffusive overshooting used in literature whereas the other grid uses CBM values guided by the results of 3D convection simulations. Interactions between the carbon, neon and oxygen shells (C-Ne-O) are common throughout both grids. The higher CBM grid also exhibits more frequent H-He and He-C interactions at lower initial masses than in the lower CBM grid. Several models also undergo multiple interaction events during evolution. While future work will be needed to fully assess the impact of the new CBM and the interactions it leads to, one expects interesting effects like unusual nucleosynthesis including more common or enhanced i- and gamma-process nucleosynthesis. Furthermore, SN precursors and a significant change to the pre-SN structure are also expected, with many models not having the commonly expected onion-ring like structure and having a different explosion probability.

The Wesenheit function is widely used to reduce the effects of interstellar reddening in distance measurements. Its construction, however, relies on the assumption of a universal extinction curve and on fixed values of the total-to-selective extinction ratio, Rv. Recent studies have shown that Rv varies significantly across the Milky Way and between different galaxies, raising concerns about systematic biases in Wesenheit magnitudes and period-Wesenheit relations. In this work, we discuss the impact of non-universal extinction on Wesenheit indices by combining the Rv-dependent extinction curve with a grid of stellar atmosphere models. We compute the integrated extinction in optical and near-infrared passbands, derive Rv-dependent R coefficients for multiple Wesenheit indices, and examine how changes in Rv propagate into Wesenheit magnitudes and Cepheid distances in our Galaxy. We find that the R coefficients in the Wesenheit functions vary strongly with Rv. For classical Cepheids in the Milky Way disk, variations of Rv within the typical observed range (2.6-3.6) can lead to substantial differences in the Wesenheit function, reaching +-0.7 mag from the mean for the Gaia-based Wesenheit index W_G and resulting in distance errors of almost 40%. Near-infrared Wesenheit indices are much less sensitive to Rv changes. Our results clearly show that accounting for variable Rv is essential when applying period-Wesenheit relations, particularly in the optical regime, or that near or mid infrared based distances should be used. While we present this effect for classical Cepheids, it applies to all pulsating stars for which period-Wesenheit relations are used to infer distances.

We present deep JWST/NIRSpec G395M spectroscopy of GLIMPSE-16043, a promising $z\sim6$ Pop III candidate originally identified through NIRCam photometry as having weak [OIII]$\lambda\lambda4959,5007$ emission. Our follow-up reveals clear [OIII] emission, ruling out a genuine zero-metallicity nature. However, the combination of the measured line fluxes and photometry indicates that its spectral energy distribution requires an extraordinarily strong Balmer jump ($-1.66 \pm 0.47$ mag) and H$\alpha$ equivalent width ($3750\pm1800$ Å), features that cannot be reproduced by current stellar+nebular or pure nebular photoionization models. The only models approaching the observations to almost within $1\sigma$ involve a hot ($T_{\rm eff}\!\simeq\!10^{4.7}$ K) single blackbody embedded in a low-$T_{\rm e}$ nebular environment, suggestive of scenarios such as a tidal-disruption event or a microquasar with strong disk winds. This cautions that photometric Pop~III selections are vulnerable to contamination when the rest-frame optical continuum is undetected. Motivated by this, we refine the photometric Pop III selection criteria to exclude the locus of extreme Balmer-jump objects. The revised criteria also recover the recently reported spectroscopic candidate AMORE6, demonstrating that the updated selection preserves sensitivity to genuine Pop III-like sources while removing key contaminants. Applying the refined criteria across legacy survey fields and five newly released CANUCS lensing cluster fields, we revisit the Pop III UV luminosity function and estimate the Pop III cosmic star-formation rate density to be $\approx[10^{-6}$--$10^{-4}]$~$M_{\odot}$~yr$^{-1}$~cMpc$^{-3}$ at $z\simeq6$--7, falling in the range of current theoretical predictions.

Cold dark matter (CDM) evolves as a collisionless fluid under the Vlasov-Poisson equations, but N-body simulations approximate this evolution by discretising the distribution function into particles, introducing discreteness effects at small scales. We present a physics-informed neural network approach that evolves CDM fields without any use of N-body data or methods, using a Kolmogorov-Arnold network (KAN) to model the continuous displacement field for one-dimensional halo collapse. Physical constraints derived from the Vlasov-Poisson equations are embedded directly into the loss function, enabling accurate evolution beyond the first shell crossing. The trained model achieves sub-percent errors on the residuals even after seven shell crossings and matches N-body results while providing a resolution-free representation of the displacement field. In addition, displacement errors do not grow over time, a very interesting feature with respect to N-body methods. It can also reconstruct initial conditions through backward evolution when sufficient final-state information is available. Although current runtimes exceed those of N-body methods, this framework offers a new route to high-fidelity CDM evolution without particle discretisation, with prospects for extension to higher dimensions.

It has recently been proved that, for constant density stars, there is a critical value $\Lambda^{*}=1$ for the dimensionless density parameter $\Lambda\equiv 4\pi R^2\rho_{\text{max}}$ of the star above which the asymptotically measured travel time $T_{\text{s}}$ along a semi-circular trajectory that connects two antipodal points on the surface of the star is {\it shorter} than the travel time $T_{\text{c}}$ along the (shorter) straight-line trajectory that connects the two antipodal points through the center of the compact star [here $\{R,\rho_{\text{max}}\}$ are respectively the radius and the maximum density of the compact astrophysical object]. This intriguing observation provides a nice illustration of the general relativistic time dilation (redshift) effect in highly curved spacetimes. One expects that generic compact astrophysical objects whose dimensionless density parameters are smaller than some critical value $\Lambda^*$ would be characterized by the `normal' relation $T_{\text{c}}\leq T_{\text{s}}$ for the travel times between the two antipodal points. Motivated by this expectation, in the present paper we prove, using analytical techniques, that spherically symmetric compact astrophysical objects whose dimensionless density parameters are bounded from above by the model-independent relation $\Lambda\leq\Lambda^*={3\over2}[1-({{2}\over{\pi}})^{2/5}]$ are always (regardless of their inner density profiles) characterized by the normal dimensionless ratio $T_{\text{c}}/T_{\text{s}}\leq1$.

For axions that couple to nucleons, the presence of dense nuclear matter can displace the axion from its vacuum minimum, sourcing large field gradients around neutron stars (and, more generally, compact objects). These gradients, which we refer to as axion hair, couple to the local background magnetic field, inducing a large voltage drop near the surface of the star; here, we demonstrate that the presence of axion hair decouples local near-field particle acceleration in the open magnetic field line bundle from the rotational frequency of the pulsar itself. This is significant as the non-observation of old slowly-rotating pulsars is attributed to the fact the rotationally-induced electric fields are not strong enough to sustain $e^\pm$ pair production. In this work, we review the evidence for the existence for `pulsar death', i.e. the threshold at which $e^\pm$ pair production (and thus, by association, coherent radio emission) ceases, and demonstrate using both semi-analytics and particle-in-cell simulations that the existence of axion hair can dramatically extend pulsar lifetimes. We show that the non-observation of extremely old, slowly rotating, pulsars allows for a new probe of light QCD and CP-violating axions. We also demonstrate how the observation of emission from both poles of pulsars with nearly orthogonal rotational and magnetic axes, as seen e.g. in PSR J1906+0746, can be used to set competitive limits on CP-violating axion-nucleon interactions.

The Laser Interferometer Lunar Antenna (LILA) presents a novel concept for observing gravitational waves from astrophysical sources at sub-Hertz frequencies. Compared to the Earth, the seismic environment of the moon, while uncertain, is known to be orders of magnitude lower, opening the possibility for achieving this sub-Hz band. This band fills the gap between space-based detectors (mHz) and Earth-based detectors (10 Hz to a few kHz). The initial version of LILA, known as LILA Pioneer, calls for non-suspended optics, relying on the moon's resonant modes to respond to gravitational waves. However, the follow-on design, LILA Horizon, requires suspensions to realize in-band free floating test masses and to filter the residual seismic background. This paper will establish baseline designs for these suspensions for different assumptions of the seismic background.

We present a new algorithm for numerical magnetohydrodynamics on staggered meshes preserving $\nabla \cdot B = 0$. Our algorithm is based on the constrained transport method and supports both cell-based adaptive mesh refinement and temporal substepping. We handle resolution changes directly on the logically Cartesian grid without needing interpolation or projection between nested or neighboring grids, nor coupling the solution between refinement levels.

Recently, a combinatorial approach to discrete, finite, and irreversibly aggregating systems has been progressively developed. In this work, we review its achievements up to the present moment, focusing on the practical aspects and discussing its limitations. First, we present the assumptions and combinatorial foundations of the approach, which are based on direct counting of the system states, in contrast to the previous approaches of Smoluchowski and Marcus--Lushnikov. A method to obtain combinatorial expressions for the average number of clusters of a given size and the corresponding standard deviation is described by solving the simplest example of a constant kernel. Then, we extend consideration to a number of kernels (e.g., additive, product, linear--chain, condensation), which were recently solved by explicitly finding the number of internal states of the cluster of a given size. Next, we show that theoretical predictions for any given kernel may be obtained with no need to find an explicit solution but using a recursive expression. We exploit this opportunity to present the use of combinatorial expressions to solve kernels related to the real processes of aerosol growth and planetesimal formation. At this point, a comparison to numerical results appears. Other potential application fields are indicated, including dust agglomeration and polymer growth. Finally, issues related to the varying precision of the theoretical predictions are summarized. In the last section, we propose open problems.

Three theoretically plausible techniques to developing a fractional scalar field cosmological model are pointed in this paper; the time-dependent kernel weighted action being then selected. Upon this choice, we proceed to establish a fractional cosmological model in $n$ dimensions considering the FLRW metric and a generalized version of the Sáez-Ballester (SB) theory. Our study focuses on the following purposes. Firstly, to investigate the fundamental gravitational structural features of the model, we analyze the dynamical behavior of the field equations, the fulfillment of the Bianchi identities, the associated conservation laws, and the application of the second Noether theorem at the background and first-order perturbation levels. Moreover, the model's distinguishing characteristics and theoretical differences from the corresponding standard scenarios are also investigated. Secondly, we aim to obtain exact analytical solutions and analyze the time evolution of key cosmological quantities, considering the integration constants influence, and the number of spacetime dimensions, and the fractional parameter effects. Furthermore, the model's predictions are compared with those of the corresponding standard models and observational data. Lastly, we propose new ideas to further generalize our model, with a focus on constructing an effective potential and investigating the conditions under which bounce solutions may emerge.