Priority

Radio continuum observations are widely used to study cosmic ray (CR) electron populations and transport processes in star-forming galaxies, but their interpretation relies on several simplifying assumptions. Here, we revisit three common assumptions: that some vertical radio profiles can be explained by CR advection alone, that radio spectra directly trace the galaxy-wide CR electron spectrum, and that bremsstrahlung and Coulomb losses are negligible for radio-emitting electrons. We model radio emission using time-dependent CR electron evolution in a magnetohydrodynamical simulation of an isolated Milky Way-mass galaxy. CR electron spectra are evolved self-consistently along Lagrangian tracer particles with the CREST framework, including injection at supernova remnants, advection with the gas, and spatially and temporally varying radiative losses. We compare these results to commonly adopted steady-state models. We find that advection-only transport in self-consistently driven galactic winds fails to reproduce the extended vertical radio intensity profiles observed in edge-on galaxies, despite reproducing the observed steepening of spectral indices with height. This is because slowly accelerating winds keep electrons in strong cooling environments for too long. Matching observed radio haloes with advection alone requires unrealistically high midplane wind velocities, implying that additional transport or re-acceleration processes are required. Although galaxy-integrated CR electron spectra at radio-emitting energies are similar across models, the resulting synchrotron spectra differ systematically because radio emission is biased toward young electrons in dense, strongly magnetised regions. Finally, we show that bremsstrahlung and Coulomb losses significantly shape radio spectra even when their loss rate is subdominant and therefore cannot be neglected.

We investigate the connection between supermassive black holes (SMBHs), their host galaxies, and large-scale dark-matter halos using broad-line X-ray AGN from the XMM–XXL and Stripe 82X surveys, together with galaxies from VIPERS and SDSS/Stripe 82. Building on the homogeneous host-galaxy catalogue presented in Paper I, we test whether AGN with a given black-hole mass, $M_{\rm BH}$, inhabit different large-scale environments from non-AGN galaxies with similar host properties. We first examine the empirical $M_{\rm BH}$–$M_{\star}$ relation of the AGN sample. We find a shallow trend with substantial scatter, likely driven by flux-limited selection effects and uncertainties in virial black-hole mass estimates. The ratio $M_{\rm BH}/M_{\star}$ decreases with increasing stellar mass, and AGN lying above and below the empirical relation show different median host properties, consistent with non-synchronous SMBH and stellar growth. We then divide the AGN into two black-hole mass bins, $8.0 \le \log(M_{\rm BH}/M_\odot) < 8.5$ and $8.5 \le \log(M_{\rm BH}/M_\odot) < 9.0$, and construct galaxy control samples matched in $M_{\star}$, SFR, and sSFR using a multivariate nearest-neighbour method. From AGN–galaxy cross-correlation functions, we infer the characteristic halo masses of AGN and matched galaxies. In the lower-$M_{\rm BH}$ bin, AGN occupy halos statistically indistinguishable from those of their controls. In the higher-$M_{\rm BH}$ bin, we find a mild indication that AGN may reside in somewhat more massive halos, with a difference of about 0.4 dex, although still consistent within the uncertainties. If confirmed with larger samples, this would suggest that halo-scale processes become important mainly at the highest $M_{\rm BH}$.

GD-1 is among the longest, coldest stellar streams in the Milky Way, making it an ideal target for probing dark matter substructure through dynamical heating. We present a catalog of 608 spectroscopically confirmed GD-1 members from the first three years of Dark Energy Spectroscopic Instrument (DESI) observations. This constitutes the largest homogeneous spectroscopic sample of GD-1, doubling the number of members previously available only through heterogeneous compilations combining multiple surveys with different systematics. Using these data, we derive updated stream tracks in sky position, proper motion, and radial velocity that extend over $100^\circ$ of the stream. We apply a Gaussian mixture model to decompose the stream into a dynamically cold thin component ($\sigma_V = 2.49\pm 0.28$ km s$^{-1}$, width $= 0.23\pm0.01^\circ$) and a kinematically hot cocoon ($\sigma_V = 6.13\pm0.75$ km s$^{-1}$, width $= 2.18\pm0.17^\circ$). The cocoon contains $\sim30\%$ of members and its velocity dispersion is consistent with $\sim11$ Gyr of heating by cold dark matter subhalos. We also detect a large proper motion dispersion ($41.36\pm4.98$ km s$^{-1}$) along the stream direction in the cocoon component. This feature indicates a significant line-of-sight distance spread in the cocoon, and its origin will be further explored in a forthcoming paper. These measurements demonstrate the power of DESI spectroscopy for characterizing the multi-component phase-space structure of stellar streams and constraining small-scale dark matter substructure.

The James Webb Space Telescope (JWST) has revealed a puzzling population of massive black holes in the first billion years, many of which are over-massive compared to their hosts (obese black holes), and reside in metal-poor hosts, posing a challenge for theoretical models at these early epochs. In this work, we compare the observational properties of astrophysically-seeded black holes using the DELPHI semi-analytic model and cosmologically-seeded primordial black holes (PBHs) using the PHANES analytic model. We explore the growth of light ($\sim 100 M_\odot$) and heavy ($\sim 10^{3-5}M_\odot$) seeds through mergers and accretion (both Eddington-limited and at super-Eddington rates) in the astrophysical scenario; PBHs (seeded between $10^{0.5-6}M_\odot$) only grow through accretion at sub-Eddington rates. Comparing to observables at $z \sim 5-10$, the only model that can be ruled out is the one where we allow Eddington-limited accretion onto light seeds. The observed high values of the black hole mass-stellar mass relation ($0.3-1$) can be reproduced by both PBHs and heavy seeds accreting at super-Eddington rates. However, only the PBH and Eddington-limited heavy seeding models can simultaneously reproduce the observed black hole masses (${\rm M_{bh}}$), stellar masses ($M_*$), and extremely low host metallicities ($Z \leq 0.01 Z_\odot$) inferred at $z \sim 7-10$. Crucially, we find PBHs show decrease in the black hole mass-stellar mass ratio with increasing halo mass at all redshifts, contrary to any astrophysical black hole model. Selecting systems at $z \sim 7$ with ${\rm M_{bh}}/M_* > 0.1$ and bolometric luminosities $\sim 10^{44-46} {\rm erg~s^{-1}}$ that show a negative black hole to stellar mass ratio and reside in $10^{9-11}M_\odot$ halos offer a promising clustering-based discriminant of PBH seeding models.

Gravitational wave (GW) detector LISA will observe near-coalescence extreme mass ratio inspirals (EMRIs), which typically form in galactic central accretion disks. Gas torques on EMRI will alter its GW-driven inspiral trajectory from the vacuum expectation, leading to potentially LISA-observable GW dephasing ($\Delta\psi_{\rm gas}$). Most studies compute $\Delta\psi_{\rm gas}$ for a thin, laminar disk, with negligible flow turbulence, where the disk exerts a fairly well-understood linear torque ($T_{\rm lin}$). However, these disks must be turbulent due to magneto-rotational instability in the inner regions. Hence, we present a proof-of-concept general, agnostic prescription for the turbulent torque ($T_{\rm turb}$) acting on an EMRI by modeling it as a Gaussian distribution around $T_{\rm lin}$, based on recent advances from a global hydrodynamical (HD) study. We compute $\Delta\psi_{\rm gas}$ for the “golden” circular EMRI with total source mass $M=10^6~{\rm M}_\odot$ and mass ratio $q=5\times10^{-5}$ in its final four-year evolution at redshift $z=0.276$ and signal-to-noise ratio (SNR) $=50$ by varying Eddington ratio ${\rm f}_{\rm Edd}$, turbulence normalization $C$ ($=~360$ in the aforementioned HD study), disk aspect ratio $h_0$, and turbo-viscous coefficient $\alpha$ in a reasonable parameters space. We find that for ${\rm f}_{\rm Edd}\gtrsim0.3$, $C\gtrsim300$, $h_0\gtrsim0.03$, and $\alpha\gtrsim0.1$, gas-induced dephasings are unobservable if only considering $T_{\rm lin}$ but could become detectable ($\Delta\psi_{\rm gas}>8/$SNR) if EMRIs exhibit chaotic migration due to turbulent gas flow. Hence, this work motivates running MHD simulations of accretion disks with embedded LISA EMRIs in the early in-spiral phase over long enough timescales to understand the evolution of their orbital elements and the imprint of the turbulent environment on their gravitational waveforms.

We report on the spectroscopic confirmation of overdense regions of massive quiescent galaxies (QGs) in the early Universe with JWST/NIRSpec. Based on data from the DeepDive NIRSpec program and archival data from the Dawn JWST Archive, we confirm three QGs in the vicinity of Jekyll & Hyde, a pair of massive QG and a dusty star-forming galaxy, at $z=3.71$ and two QGs around SXDS-27434 at $z=4.01$. According to the analysis of galaxy number density with photometric redshifts, Jekyll & Hyde (SXDS-27434) are in an overdense region, where the number density of galaxies is three times higher than the average in the COSMOS (SXDS) field. SED fitting suggests that most of the QGs follow similar star formation histories and have consistent formation and quenching epochs. The same trend is observed in other proto-clusters hosting QGs that were already identified by ground-based telescopes, indicating that the large-scale environment plays an important role in the formation of QGs. In addition, JWST spectra reveal a broad H$\alpha$ emission line from SXDS-27434 and faint emission lines from other three QGs, which are identified as AGN-driven based on their emission line ratios. The overdensity is also reproduced by the Illustris TNG300 simulation at $z=3.71$, in which the member QGs also have similar quenching epochs. These results suggest that large-scale structure may enhance merger activity and/or gas accretion and trigger AGN feedback, which simultaneously drives galaxy quenching in the overdensity.

The internal structure of dark matter halos on sub-galactic scales remains a key open question, particularly in the context of the core-cusp problem. Ultra-faint dwarf galaxies (UFDs), owing to their extreme dark matter dominance, provide a promising laboratory to probe these density profiles through stellar tracers. In this work, we assess the capability of the Chinese Space Station Telescope (CSST) to detect and characterize wide binary stars in the nearby UFD Segue 1, using mock observations. We generate mock binary populations based on our existing $N$-body simulations and incorporate realistic CSST observational conditions, including the expected deep-field limiting magnitude ($g \sim 27.5$ mag) and a photometric completeness of approximately $90\%$. The two-point correlation function (2PCF) of stellar pairs is used as a statistical tool to recover the binary fraction under these assumptions. We find that CSST can robustly detect wide binaries at the $3\sigma$ level for binary fractions as low as $f_b \gtrsim 0.01$, provided a stellar sample size of $N_{\mathrm{star}} \gtrsim 2300$. However, distinguishing between cusped and cored dark matter profiles is significantly more demanding, requiring $N_{\mathrm{star}} \gtrsim 6000$ and $f_b \gtrsim 0.1$ within $\sim 40\mathrm{kpc}$.

Various processes can induce long-lived overdense rings and arcs in protoplanetary and AGN accretion discs, such as the accumulation of gas at the outer edge of the dead zone, or the infall of material. Using the local approximation of dynamical friction, we investigate the orbital evolution of a low-mass highly-eccentric point-mass accretor (perturber) embedded in an isothermal disc hosting a density ring. We specifically consider the regime in which the eccentricity exceeds four times the disc aspect ratio. For prograde perturbers, orbits that cross the ring progressively circularize while their semi-major axes converge toward the ring radius. As a result, perturbers accumulate, forming a population ring superimposed on the gaseous ring. The ring therefore acts as a migration trap for these eccentric orbits. We also find that prograde orbits tangent to the ring, either at apocentre or pericentre, remain tangential throughout their evolution; perturbers confined to these trajectories experience the highest accretion rates. In contrast, retrograde perturbers always migrate inward. Once the semi-major axis becomes smaller than the ring radius, the eccentricity grows, but not enough for the orbit to intersect the ring again. We also discuss how feedback effects, such as jet launching and thermal torques, could modify the effective forces acting on the perturbers.

Models predict that chemical enrichment and gas redistribution should proceed rapidly once star formation begins, yet direct observational constraints at the earliest cosmic epochs have been scarce. Here we present evidence that metal-enriched gas in multiple ionic phases was already present around galaxies before the midpoint of cosmic reionization. Using JWST/NIRSpec rest-frame ultraviolet spectroscopy of three galaxies at redshifts $z=7.2-9.3$, we detect blueshifted metal absorption in all three systems; across the sample, the detected transitions span neutral, low-ionization, and high-ionization species, including O I, Si II, C II, Si IV, and C IV. These absorption features show velocity offsets of order $|\Delta v| \sim 50$–$250\,\mathrm{km\,s^{-1}}$, predominantly blueshifted relative to the systemic redshifts of the host galaxies derived from nebular emission lines. This ionic coexistence within a broadly shared velocity structure, together with the observed equivalent-width ratios, is consistent with outflowing or otherwise kinematically disturbed galaxy-associated gas, similar to that seen at lower redshift. The observations therefore indicate that metal-enriched gas associated with galaxies was already kinematically disturbed at very early times, requiring rapid metal production in the early generations of stars. These results show that key conditions for baryon cycling were established in at least a subset of luminous galaxies within the first several hundred million years of cosmic time, well before the completion of reionization.

We investigate whether the observed radio-infrared-$\gamma$-ray correlation in star-forming galaxies is a geometric effect rather than a signature of local cosmic-ray (CR) calorimetry. Using the GALPROP framework, we generate synthetic observations for external viewers from a grid of 3D Milky Way models with varied CR source, gas, interstellar radiation, and magnetic field distributions, all normalised to reproduce local CR data. We find that a tight, quasi-linear correlation arises naturally from line-of-sight integration through the extended, radially-structured disc, even when local calorimetry is absent. The correlation's properties depend strongly on viewing geometry, preserving its form under moderate inclination but breaking down in edge-on views where galactic components are stratified. We conclude that the correlation is primarily an emergent property of geometric projection, not local physics. This implies that its scatter is likely not random noise but a diagnostic of underlying galactic structure and viewing angle.

Calcium-aluminum-rich inclusions (CAIs) in carbonaceous chondritic meteorites are the oldest relics in the solar system. Notably, their radiogenic age feature a brief (100 kyr) condensation episode. In contrast, the reservoirs of the short-lived isotopes in CAIs, presumably supernovae or asymptotic giant stars, pollutes star-forming regions in giant molecular cloud complexes (GMC) over much longer (Myr) duration. Through a series of numerical simulations, we show here the possibility that, within an extended region (2$\sim$3 AU), nearly all “pre-solar” CAI-loaded grains in the infall clouds were sublimated and re-condensed during the early ($ \lesssim 10^5$ yr) infall and formation of class-0 disks. We adopt a set of initial conditions from a previous hydrodynamic simulation of the collapse of GMC and the formation of young stellar clusters. We analyze the evolution of the disk's thermal distribution and dynamical structure resulting from the interaction between circumstellar disks and infalling gas. Our follow-up simulations, with much higher resolution, show significant and rapid changes in the disk orientation and morphology due to the dynamic infall of external streamers. Warps and global spiral density waves commonly appear. They lead to intense dissipation which heats the gas to sufficiently high temperature to sublimate prior-generation CAIs. This solid-to-gas phase transition is followed by subsequent cooling and re-condensation. The CAI contained in the meteorites today could be the relics of the last episode of major infall onto class 0 disks.

The estimation of stellar atmospheric parameters for large-scale samples, particularly metal-poor stars, is a cornerstone of Galactic archaeology. In this work, we optimized a photometric filter design tailored to measuring stellar metallicities for very metal-poor stars with [Fe/H]$< -1$.The optimal configurations consist of a central wavelength $\lambda_{\rm c}$ = 3960 Angstrom with a bandwidth $\Delta\lambda$ = 80 Angstrom for giant stars, and $\lambda_{\rm c} $= 3920 Angstrom with $\Delta\lambda$ = 80 Angstrom for dwarf stars. By applying these optimized filters to synthetic photometry derived from Gaia XP spectra, we inferred metallicities for both populations. Both internal and external validations demonstrate high precision across a wide metallicity range: 0.18-0.19 dex for $-2 \le \rm [Fe/H] \le -1$, 0.23-0.33 dex for $-3 \le \rm [Fe/H] \le -2$, and approximately 0.39 dex for the most metal-poor regime, successfully extending down to $\rm [Fe/H] \approx -4$ for giant stars, $\rm [Fe/H] \approx -3.3$ for dwarf stars. Finally, we present a catalog of approximately 14.5 million metal-poor stars with robust $\rm [Fe/H]$ measurements, along with more than ten thousand red giant ultra metal-poor candidates with $\rm [Fe/H] < -4.0$, providing a valuable resource for exploring the early formation and chemical evolution of the Milky Way.

Most Active Galactic Nuclei (AGN) are Radio Quiet, with radio emission that may arise from star-formation activity, AGN-driven winds, weak jets, and coronal activity. Disentangling these mechanisms is challenging and requires detailed multi-wavelength investigation, but it is crucial for quantifying AGN feedback in galaxy evolution. We present a detailed radio investigation of 21 X-ray selected AGN in the Supermassive Black Hole Winds in X-Rays (SUBWAYS) sample (log Lbol = 44.9-46.3 erg/s, z=0.1-0.5), selected to systematically search for Ultra-Fast Outflows (UFOs). UFOs are detected in 30% of the targets, making the sample particularly well-suited for investigating the role and signatures of multi-scale outflows at different frequencies. We build the radio SED of the sources complementing our proprietary data, collected with the JVLA at 1.5 and 6 GHz, with images from LoTSS and other publicly available radio surveys between 150 and 1400 MHz. We investigate the role and occurrence of the aforementioned mechanisms, with particular interest in outflows and their possible relation with UFOs. We combined information on spectral indices, luminosities, and morphologies of the radio emission with properties derived in other wavebands, such as Star Formation Rate, X-ray luminosity, Eddington ratio or the UFO kinetic luminosity. All the sources are detected and are mostly consistent with RQ AGN. For 80% of the sources the data suggest the presence of an outflow (wind or weak jet). Interestingly, our results indicate that AGN with UFOs tend to have larger radio extension and a steep radio spectrum consistent with outflows. Moreover, the radio emission of the 6 UFO hosts is consistent with predictions from wind-driven shock models, possibly indicating a direct connection between the two phases. Alternatively, this may reflect physical conditions favouring the rise of both phenomena.

We present a new implementation of the SFUMATO code, called SFUMATO#, for solving self-gravitational radiation hydrodynamics problems using adaptive mesh refinement (AMR) with the CUDA/HIP programming frameworks. The code incorporates a multigrid solver for self-gravity, radiation transfer with M1 closure and reduced speed of light approximation, non-equilibrium chemistry, thermal evolution, and sink particle schemes. We develop new non-equilibrium chemistry and thermal solvers based on a linearized implicit method, whose accuracy is validated through a series of test problems by comparison with solutions obtained using the Newton-Raphson method. By incorporating the heat capacity of dust grains, the dust temperature can be evolved without iterative energy-balance calculations. From the perspective of computational cost, we demonstrate that adopting an increased pseudo dust heat capacity accelerates the chemistry solver while preserving accuracy, even when the value is increased by up to three orders of magnitude relative to the realistic value. In addition, we perform a suite of test problems to confirm the validity of the other components of our implementation. The code supports multi-GPU execution via MPI-based parallelization. We measure the strong-scaling performance of the hydrodynamics and self-gravity solvers on both uniform and AMR grids, as well as the overall code performance using a giant molecular cloud simulation. We find that the computational cost of the self-gravity solver increases with the number of MPI processes, indicating that efficient parallel performance is achieved only when the number of devices is chosen such that the cost of the self-gravity solver remains comparable to that of the other components.

We present ultra-deep ($\simeq 20-30$ hours), rest-frame UV spectroscopy with NIRSpec/JWST of two UV-bright galaxies at $z\sim 8.7$, CEERS-1019 and CEERS-1025 ($Z_{\rm neb} \simeq 0.1 Z_{\odot}$). The spectra reveal exceptionally strong P-Cygni profiles in wind lines (NV $\lambda$1240 and CIV $\lambda$1550) and significant broad and strong HeII $\lambda$1640 emission ($\rm EW\simeq 2-4$ A). We compare the observations with synthetic stellar population models at $Z_{\star} \simeq 0.1 Z_{\odot}$, both including and excluding very massive stars (VMS). Models including VMS provide a markedly improved fit to the data relative to non-VMS models ($\Delta$AIC and $\Delta$BIC $> 70$), which fail to reproduce the observed strengths of the wind features. A comparison with empirical spectra of VMS-dominated systems in the local Universe further supports this interpretation. The best-fit VMS models imply extremely young ages of the stellar populations ($\simeq 1.5-2.0$Myr) and high ionizing photon production efficiencies ($\log \xi_{\rm ion} [\rm Hz erg^{-1}] \gtrsim 25.8$), exceeding those inferred from non-VMS models by $\sim 0.1-0.2$ dex. These results provide evidence for an overabundance of VMS at high-$z$ with an IMF extending well beyond $100 M_{\odot}$, and highlight their potential role in shaping the rest-frame UV spectra, chemical enrichment, and ionizing output of galaxies in the early Universe.

JWST has revealed a population of super-luminous early galaxies with a volume density in excess of most expectations. The spectra reveal diverse properties: while some reveal strong emission lines characteristic of galaxies in the midst of strong bursts, others show weak emission lines that could reflect old stellar populations, large escape fractions, or post-burst star formation histories. Through the JWST Cycle 4 large program SPURS, we have obtained ultra-deep (29 hr) rest-frame UV spectroscopy of a z=9.3 super-luminous ($M_{\rm UV}=-21.66$) galaxy with large assembled stellar mass (1.6$\times$10$^9$ $M_\odot$) and extremely weak emission lines (H$\beta$ EW $\approx25$ Å). The strong stellar wind features and rest-optical line ratios suggest the galaxy is already significantly enriched, with a metallicity of 0.4–0.7 Z$_\odot$. The interstellar absorption lines reveal outflows ($v\simeq -161$ km s$^{-1}$) with a large neutral gas covering fraction, suggesting that the weak emission lines are not due to large escape fractions. The combination of the Balmer break, weak emission lines, and stellar wind features constrains the star formation history, indicating a recent burst of star formation lasting 10–20 Myr followed by a downturn over the last 10 Myr. The observations suggest that $z\gtrsim 9$ weak emission line galaxies such as this source can be explained by stochastic star formation, provided that the downturns in star formation are recent (i.e., <10 Myr prior to observation). The ultra-deep grating spectrum enables the IGM damping wing to be characterized, decoupling the effects of local absorption. The smooth Ly$\alpha$ break indicates that this source, one of the most massive galaxies known at z>9, is likely situated in a small ionized bubble ($0.29_{-0.09}^{+0.11}$ pMpc), as is common at large neutral hydrogen fractions ($\bar{x}_{\rm HI}=0.81_{-0.21}^{+0.14}$).

We report the discovery of DESI-HVS1, a hypervelocity star (HVS) candidate identified from DESI DR1 spectroscopy and Gaia DR3 astrometry. DESI-HVS1 is an old, low-mass, metal-poor F-type star with a mass of $0.8\,M_\odot$, an age of $\sim14.1$ Gyr, and $\mathrm{[Fe/H]}=-1.6$. It is located at a heliocentric distance of $3.77^{+0.39}_{-0.36}$ kpc and has a Galactocentric total velocity of $523^{+46}_{-47}\,\mathrm{km\,s^{-1}}$, marginally exceeding the local escape speed, corresponding to an unbound probability of $P_{\rm ub} \sim 50\%$. Backward orbit integrations show that DESI-HVS1 had a closest approach to the Galactic Centre (GC) of $0.40^{+0.23}_{-0.11}\,\mathrm{kpc}$, with a velocity of $682^{+22}_{-35}\,\mathrm{km\,s^{-1}}$ and a flight time of $12.89^{+0.92}_{-0.74}\,\mathrm{Myr}$. The reconstructed orbit exhibits a clear perigalactic turning point and only a single crossing of the Galactic midplane ($P_{\rm cross} > 0.95$). These properties suggest that DESI-HVS1 is most naturally explained by the Hills mechanism, although alternative scenarios cannot be entirely ruled out. Its discovery provides the first strong evidence for an old, low-mass HVS candidate consistent with a GC origin, indicating that the apparent dominance of young, massive GC-origin HVSs is likely a consequence of observational selection effects.

We investigated the high-mass star formation activity in a subregion of the Sagittarius E star-forming complex, centered at (l,b) = (358.69 deg, 0.03 deg), where infrared and radio sources trace a prominent U-shaped structure that has not been identified in previous studies. We used radio continuum data from the Global View on Star Formation (GLOSTAR) survey, which is a wide-band radio (4-8 GHz) survey of the Milky Way that combines data from the Karl G. Jansky Very Large Array and the Effelsberg 100 m telescope. Using BLOBCAT source extraction software, we identified 49 compact radio sources. Based on multiwavelength associations and spectral index estimates, we identified GLOSTAR counterparts to 27 previously confirmed HII regions, detected radio emission from 3 WISE "radio-quiet" candidates, and report 5 new HII region candidates. The derived physical properties indicate that most are relatively evolved HII regions. We find around 50 cold dust clumps, predominantly toward the south and southeast. Mid-infrared flux-ratio maps ([4.5]/[3.6]) show localized shock enhancements along the arc and adjacent clumps, and 15 clumps exhibit SiO emission with broad components indicative of shocks. Together with CO data, the SiO velocity components delineate a continuous (>100 km/s) velocity bridge that links the far dust-lane inflow to the central molecular zone (CMZ) stream. The largest concentration of clumps and compact HII regions lies at this interface. These combined diagnostics favor a scenario in which bar-driven cloud-cloud collision at the far dust-lane-CMZ interface compressed the gas and triggered the observed high-mass star formation.

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for cosmic ray transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a $20000 M_{\odot}$ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inwards and compress the gas into higher density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the SFE in the CRT simulation including stellar wind CRs is 43 \% higher than the non-CRT simulation. The IMF is also top heavy in comparison, with a slope above 1 $M_{\odot}$ that is shallower by $\sim 20$ \%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16 \% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by $\sim 10$ \%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the galactic center.

(Abridged) We aim to investigate the inner regions of large and massive disks orbiting T Tauri stars, thought to be progenitors of systems with wide-orbit planets and possible cases of halted pebble drift. We analyze the MIRI spectra of three disks from the MINDS program: V1094 Sco, DL Tau, and IM Lup. The spectra reveal a striking diversity. V1094 Sco and DL Tau exhibit the highest C$_2$H$_2$/H$_2$O flux ratio in the MINDS sample of T Tauri disks. In V1094 Sco, even cold C$_4$H$_2$ is seen. In contrast, the IM Lup spectrum is dominated by O-bearing species. No one-to-one correspondence is found between the gas in the outer disk, as traced by the C$_2$H/C$^{18}$O flux ratio, and that of the inner disk as traced by the C$_2$H$_2$/H$_2$O flux ratio. To explain these results, we propose a scenario based on a toy model of halted pebble drift. We show that a volatile C/O ratio close to unity and low C and O abundances in inner disks arise only if: (1)  95$\%$ of the icy grains are blocked in the outer disk, (2) the outer disk is chemically evolved, and (3) the gas in the outer disk has had time to reach the inner disk. DL Tau and perhaps V1094 Sco would be the rare examples for which all these conditions are met. Therefore, a high C$_2$H$_2$/H$_2$O flux ratio in pebble-rich disks would have a different origin than proposed for very-low mass stars, for which fast drift of O-rich pebbles would eventually leave a C-rich inner disk. We also show for the first time that the disks with high C$_2$H$_2$/H$_2$O flux ratio exhibit a prominent silica dust component, a result found in four disks published so far (V1094 Sco, DL Tau, CY Tau, DoAr 33). We propose that the reformation of dust at the sublimation front of silicates in a gas with super-solar (but below unity) C/O ratio leads to a silica stoichiometry (SiO$_2$). In turn, silica is a promising diagnostic of the C/O ratio in the inner disks.

The Dyson Minds 2025 Workshop, held at the Center for Brains, Minds & Machines at MIT and organized by Penn State, MIT, and The Ultraintelligence Foundation, brought together researchers in astrophysics, engineering, artificial intelligence, computer science, and philosophy to examine "Dyson Minds" – large-scale post-biological intelligences powered by energy harvested from supermassive black holes (SMBHs). Building on the ideas of F. J. Dyson (1960, 1966) and I. J. Good (1966), participants explored the physical, engineering, behavioral, and observational consequences of civilizations embodied as machinery operating near the universe's most powerful energy sources. The workshop aimed to develop new observational strategies capable of detecting signatures of such systems. Despite the highly cross-disciplinary scope, discussions centered on how a Dyson Mind might be constructed, how it might behave, and how those factors would shape strategies for the search for extraterrestrial intelligence. Key themes included the thermodynamic, mechanical, and stability limits of Dyson swarms; the trade-offs between power availability and communication latency in distributed minds; and how observability changes depending on whether Dyson Minds act as coherent entities or as loosely coordinated collectives. Across these topics, the consensus was that details of architecture and behavior strongly influence observational signatures. A major recommendation was to apply anomaly-detection methods to archival datasets, including those from WISE, JWST, and the Event Horizon Telescope, to identify unusual sources potentially overlooked by standard reduction pipelines. By integrating insights from multiple disciplines, the meeting advanced concrete, observation-focused strategies for future technosignature searches around SMBHs.

Hydrogen atom tunneling likely plays a substantial role in the gas-phase chemistry of astrochemical environments. To determine the potential effect that it has on the chemical modeling of aromatic molecules, we screened the kida.uva.2024 network, and our own expanded network to find reactions which could be significantly accelerated by hydrogen atom tunneling in the ISM. In total, 64 reactions were identified. The hydrogen abstraction reactions from H$_{2}$ to four key interstellar radicals (C$_{2}$H, OH, CN, and NH$_{2}$) were studied further using newly calculated potential energy surfaces and RRKM analyses to determine rate coefficients for a temperature of 10 K and a density of 2 $\times$ 10$^{4}$ cm$^{-3}$. Despite having low rate coefficients of 1.66 $\times$ 10$^{-15}$, 8.17 $\times$ 10$^{-16}$ and 3.15 $\times$ 10$^{-16}$ $\mathrm{cm^{3}\,s^{-1}}$ the C$_{2}$H, OH, and CN reactions are competitive in the ISM, due to large overall rates caused by the high abundance of molecular hydrogen. The calculated value for the NH$_{2}$ reaction, however, was much smaller and found to be inefficient at ISM conditions. The possible effects of all other considered reactions were studied with simulations using calculated collision limit rate coefficients. Upper and lower bounds were then placed on modeled aromatic abundances using the most significant reactions. Due to the dependence of calculated aromatic abundances on reactions involving c-C$_{6}$H$_{5}^{+}$ and the recent questions surrounding its reactivity, we also explored the abundance variations caused by reactions leading to or involving c-C$_{6}$H$_{5}^{+}$.

Accurate modeling of ion-molecule reaction networks is essential for understanding the chemical evolution of planetary ionospheres, particularly for giant planets where proton-transfer chains drive atmospheric composition. However, predicting reaction rates in these ultracold environments remains a challenge due to the non-trivial interplay between vibrational dynamics and quantum tunneling. In this work we present a chaos-diagnostic framework that integrates multireference electronic structure theory, Adiabatic Gauge Potentials (AGP), and Random Matrix Theory (RMT) to characterize the microscopic dynamics of proton transport. Using the formation of H+3 and the proton-bound cluster H+5 as representative model systems relevant to Jovian atmospheres, we demonstrate that the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity. This diagnostic approach offers a generalizable, data-driven metric for identifying vibrationally gated pathways in complex astrochemical networks, providing a theoretical basis for refining kinetic models of planetary and interstellar plasmas

Dark matter search strategies have started advancing towards the neutrino fog. In this regard, compact objects such as neutron stars have already demonstrated their ability in probing such low DM-nucleon cross-sections from dark matter induced effects. In the optically thin limit, effect of dark matter self-interaction becomes relevant and may assist the capture and thermalization of dark matter inside stars, imparting observable changes on neutron star temperatures. The resulting radiation although weak can be potentially detected by the James Webb Space Telescope and upcoming Thirty Meter Telescope and the European Extremely Large Telescope. Observation of cold neutron stars accompanied by advancements in direct detection probes would provide stringent constraints or a smoking-gun signature for dark matter self-interactions. The potential detection of a neutron star with surface temperatures $\sim (1000 - 1200)$ K in the optically thin limit can push the bounds on asymmetric dark matter self-interaction cross-section to approximately two orders of magnitude more stringent than the bullet cluster.

3C 454.3 is a flat spectrum radio quasar (FSRQ) known for its high variability across the electromagnetic spectrum, showing structural and flux variability in its pc-scale jet, and correlated variability among frequency bands. This study aims to identify the structure, dynamics, and radiative processes common to the innermost regions of the blazar 3C 454.3. We investigate whether any jet component can be associated with $\gamma-$ray emission and variability. We analyze the relationship between the variable $\gamma-$ray emission and pc-scale jet properties in 3C 454.3 by combining $\gamma-$ray data spanning twelve years with contemporaneous VLBA multi-epoch images at 15 and 43 GHz. Spearman rank correlation tests are conducted to determine if the flux variability of any jet component is associated with $\gamma-$ray variability. Core emission at 43 and 15 GHz strongly correlates with $\gamma-$ray emission. The 43 GHz core (Q0) contributes around 37$\%$ of the observed $\gamma-$ray variability, while the 15 GHz core (K0) accounts for 30$\%$. A quasi-stationary component at 43 GHz, at a projected distance of 4.6 pc, correlates with the $\gamma-$ray flux, accounting for 20$\%$ of its emission between 2016 and 2021. We found a mobile component (Q3 between 2010.18 and 2011.16) at 43 GHz with a projected distance between 0.8 and 2.3 pc and apparent velocity of $\beta_{app} = 9.9 \pm 1.1$ c, accounting for approximately 28% of the $\gamma-$ray emission. The observed simultaneous variability in emission regions beyond the central parsec strongly suggests synchrotron self-Compton (SSC) as the primary mechanism for $\gamma-$ray production in these regions. Our findings demonstrate the existence of multiple $\gamma-$ray emission regions within the blazar jet but also suggest that some of these regions are non-stationary over time.

Galaxy mergers are known to trigger bursts of central star formation, which should therefore lead to stellar mass growth in their inner regions. However, observational measurements of this `burst mass fraction' are scant. Here, we assemble a large sample of  14,000 post-coalescence galaxies that have recently completed their merger-induced star formation, and compare various measurements of central stellar mass with a matched control sample. Specifically, we quantify (at fixed redshift, star formation rate and total stellar mass) the stellar mass enhancement within a fixed angular aperture (Delta M_{star,fibre}) and in the galactic bulge (Delta M_{star,bulge}), finding burst mass fractions of 10 – 20%. 61 galaxies in our sample are at z<0.05 and have integral field unit data from the Mapping Galaxies at Apache Point (MaNGA) survey, allowing further kpc-scale assessment of excess stellar mass and radial gradients. When assessed within apertures defined in units of kpc we again find a   15 – 20% excess of stellar mass in the central regions of the post-mergers compared with matched controls. However, within apertures defined in units of effective radius this stellar mass enhancement increases to 40%, suggesting that the relative structure/size of the galaxy is important for regulating the location of the merger induced star formation. Moreover, we find that these stellar mass enhancements are spatially extended, out to  7 kpc or around 1 R/R_e, although the small sample size of the MaNGA overlap limits our radial sampling. Our work represents the first direct measurement of merger-induced stellar mass growth that is independent of stellar population modelling, or fitting light profiles, demonstrating significant and extended mass build-up in late stage post-mergers.

Ultra-high-energy photons play an important role in probing astrophysical models and beyond-Standard-Model scenarios. We report updated limits on the diffuse photon flux using Telescope Array's Surface Detector data collected over 14 years of operation. Our method employs a neural network classifier to effectively distinguish between proton-induced and photon-induced events. The input data include both reconstructed composition-sensitive parameters and raw time-resolved signals registered by the Surface Detector stations. To mitigate biases from Monte Carlo simulations, we fine-tune the network with a subset of experimental data. The number of observed photon candidates is found to be consistent with the expected hadronic background, yielding upper limits on photon flux $\Phi_\gamma(E_\gamma > 10^{19} \text{eV}) < 2.3 \cdot 10^{-3} $, and $\Phi_\gamma(E_\gamma > 10^{20} \text{eV}) < 3.0 \cdot 10^{-4} $ $ (\text{km}^2 \cdot \text{sr} \cdot \text{yr})^{-1} $.

The third data release of the Gaia mission (Gaia DR3) has enabled large-scale searches for dormant black hole and neutron star binaries with stellar companions at AU-scale separations. A recent study has proposed thousands of dormant black hole and neutron star binary candidates using summary statistics from Gaia DR3 by simulating and fitting Gaia observables. In this work, we perform broadband spectral energy distribution (SED) fitting from the optical to the infrared for 1,328 candidates, incorporating GALEX ultraviolet photometry to assess the presence of hidden hot companions. We quantify ultraviolet excess by comparing observed near-ultraviolet fluxes with single-star SED predictions and further test whether excesses can be explained by non-degenerate stellar companions for sources exhibiting moderate excess. We additionally examine the Galactic kinematics of the sample to identify systems potentially affected by natal kicks during compact-object formation. By combining the ultraviolet and kinematic diagnostics, we identify 182 sources as the highest-priority candidates for follow-up observations, in which 19 are black hole candidates with fit_companion_mass $\geq$ 3 $M_\odot$.

High-$z$ quasars are believed to reside in massive dark matter haloes (DMHs), suggesting that they reside in galaxy overdense regions. However, previous observations have shown a range of environments around them. These fields have been limited to luminous quasars ($M_{1450}\lesssim-25$), for which photoevaporation may hinder galaxy formation in their vicinity. Here, we present Subaru/Hyper-Suprime Cam observations of the environments of four low-luminosity quasars ($-24<M_{1450}<-22$) at $z\sim6.18$, which are expected to have a smaller photoevaporation effect. We detect Lyman $\alpha$ emitters (LAEs) with narrowband NB872 imaging, and measure the local LAE overdensity. One quasar (J0844$-$0132) resides in an overdense region ($\delta_\mathrm{LAE}=1.97\pm0.40$), whereas the other three fields are consistent with no overdensity. These results hold over the proximity zone of each quasar, suggesting that the diverse environment around quasars is independent of photoevaporation. We find no significant correlation between the LAE overdensities and the characteristics of host galaxies and supermassive black holes. Our quasars have host stellar mass measurements from JWST, allowing us to compare them with the LAE overdensity around galaxies without quasar activity with comparable stellar masses. We find that the LAE overdensity in the J0844$-$0132 field is stronger than that of galaxies with similar stellar mass at $z\sim6$.

The dynamical nature and formation mechanism(s) of galactic spiral arms remain long-standing problems in astrophysics. Most theoretical work is based on analytic calculations or idealised simulations, which has yielded several theories of spiral structure. The radial profile of the spiral arm rotation speed - the pattern speed - is a key observable prediction of these theories. However, observations that infer spiral pattern speeds reveal a mixed picture with no clear consensus. Here, we expand on theoretical efforts by examining the pattern speed profiles in the Auriga Superstars set of high-resolution cosmological magnetohydrodynamic simulatons of Milky Way-mass spiral disc galaxies. These simulations combine galaxy formation in a cosmological environment with the high dynamical fidelity afforded by an $\sim 800$ $\rm M_{\odot}$ star particle resolution, giving $\sim 100$ million star particles in the disc. We show that several different spiral arm theories are realised among our simulations, including large-scale kinematic density waves, manifold spirals, dynamic (co-rotating) spirals, and overlapping modes. In particular, we demonstrate that a strong tidal interaction leads to clear kinematic density waves, and that manifold spirals are present in a strongly-barred galaxy. Interestingly, we find that the same galaxy may show qualitative evolution of their spiral pattern speed profiles, indicating that the nature of spiral arms can evolve on potentially sub-Gigayear timescales. Our results demonstrate that in the absence of a strong external encounter or a strong bar, galactic spiral structure is highly transitional and complex with no clear long-lived underlying wave.

Stellar microlensing surveys are a powerful tool for probing dark matter in the form of planetary and stellar mass compact objects (COs), in particular primordial black holes (PBHs). Under standard assumptions, current observations exclude COs in the mass range $10^{-11} \lesssim M/M_{\odot} \lesssim 10^{4}$ making up all of the dark matter. We provide an overview, aimed at theorists working on PBHs, of the history, theory, observational status, and future prospects of the field.

We investigate how misalignment between the core angular momentum and the large-scale magnetic field affects protostellar outflows, and whether a single protostellar system can drive multiple outflow components. We perform three-dimensional nonideal magnetohydrodynamic simulations of magnetized rotating cores, focusing on the formation of a protostar, a circumstellar disk, and magnetically driven outflows. The initial angle between the core angular-momentum vector and the magnetic field is systematically varied from $0^\circ$ to $90^\circ$. All models launch a classical magnetocentrifugal disk wind (DW) roughly along the local disk normal. For large misalignment, the system also develops a spiralflow (SF) component that propagates parallel to the disk plane. In a representative model with a $60^\circ$ misalignment, the outflow transitions from a DW-dominated to an SF-dominated phase, with the SF becoming more massive and more extended than the DW, and the two components intermittently coexisting. Across the model suite, the maximum mass and size ratios of SF to DW, as well as the relative lifetimes of the two components, increase for misalignment angles $\gtrsim60^\circ$. We propose that secondary, misaligned outflows (or their fossil remnants) observed in some protostellar systems can be interpreted as the SF component, while the main bipolar outflow traces the DW from the same misaligned system.

The structures and dynamics of the interstellar medium are governed by a combination of self-gravity, external gravity, and various sources of ordered and random motions on different spatial scales. This paper uses ALMA CO (3-2) observations at 0.1" $\approx$ 5 pc resolution to examine the scale dependence of molecular gas structure and dynamics in the central molecular zone (CMZ) of a nearby galaxy, NGC 3351. We use the dendrogram technique to characterize hierarchical molecular gas structures spanning two decades in spatial scales and measure their size, gas mass, and velocity dispersion. Their size-linewidth relation shows a power-law slope of 0.58, comparable to measurements for CMZs in other galaxies and suggestive of significant contribution from ordered motion on large scales. We further decompose the observed velocity dispersion in each gas structure into ordered versus random motions. The former appears stronger in gas structures at $\gtrsim$ 30 pc while the latter becomes more dominant at $\lesssim$ 30 pc. Modulo uncertainties with the CO-to-H$_2$ conversion factor, the estimated gravitational free-fall time is comparable to the crossing time of ordered motions for structures on all spatial scales, and both becomes longer than the crossing time of random motions at small, $\lesssim$ 10 pc scales. Our results highlight the varying sources and drivers of gas motions on different spatial scales in the CMZ of a Milky Way-like galaxy.

The recent detections of thioformyl cyanide (HCSCN) and propynethial (HCSCCH) in TMC-1 provide critical insights into the interstellar sulfur inventory, yet their sequestration and survivability on dust grain mantles remain poorly constrained. Here, we present a computational study of the site-specific adsorption of HCSCN and HCSCCH on amorphous solid water (ASW), modelled via water clusters (H2O)n, n = 6-16, at the wB97X-D/def2-TZVP level of theory, corroborated by QTAIM topological analyses and TD-DFT vertical excitations. Our results reveal a highly heterogeneous binding environment, with desorption energies spanning 1500 to 4900 K. Strongly bound cavity sites induce significant Stark shifts in the C=S stretching modes. Crucially, while the ice matrix exerts a negligible solvatochromic shift on UV transition wavelengths, deeply bound CN-cavity configurations exhibit a pronounced hyperchromic enhancement of the oscillator strength. Implementing these site-specific parameters into the UCLCHEM gas-grain code demonstrates that these species undergo a gradual thermal desorption profile rather than a singular sublimation event. Furthermore, the hyperchromic effect establishes a Survival Paradox: while deeply trapped populations are thermodynamically shielded against thermal desorption, they simultaneously possess enhanced UV absorption cross-sections, rendering them vulnerable to photodissociation by the interstellar radiation field prior to sublimation.

When a radio jet is partially optically thick in the launching region, its apparent compact core may display frequency-dependent positional shifts. High-precision astrometric measurements of core shifts enable astronomers to pinpoint the jet's origin and place tight constraints on the magnetic field. BL Lacertae, the archetypal BL Lac object, hosts a highly variable and well-collimated jet. To independently constrain its innermost core shifts, we conducted very long baseline interferometric (VLBI) observations at 8.4, 12.4, 15.2, 23.6, and 43.2 GHz. By exploiting a nearby (13.3 arcmin) steep-spectrum calibrator (NVSS J220340+420839) through inverse phase-referencing VLBI astrometry, we detect nearly unbiased two-dimensional core shift measurements with state-of-the-art precisions of 5-30 $\mu$as, which are significant at $>3\sigma$ confidence. The core shift between 8.4 and 43.2 GHz reaches 250 $\mu$as. The apparent core shifts scale with frequency as $\nu^{-1/k_r}$, implying the existence of an optically thick region in the upstream of jet. The derived core-shift index, $k_r\!=\!1.18^{+0.59}_{-0.34}$, is consistent, within uncertainties, with the canonical $k_r\!=\!1$ expected under energy equipartition between the jet particle and magnetic field energy densities, while allowing for modest deviations given that BL Lacertae was captured in a flaring state.

The surface and subsurface of worlds beyond Mars remain largely unexplored. Yet these worlds hold keys to fundamental questions in planetary science - from potentially habitable subsurface oceans on icy moons to ancient records preserved in Kuiper Belt objects. NASA's success in Mars exploration was achieved through incrementalism: 22 progressively sophisticated missions over decades. This paradigm, which we call Planetary Exploration 2.0 (PE 2.0), is untenable for the outer Solar System, where cruise times of a decade or more make iterative missions infeasible. We propose Planetary Exploration 3.0 (PE 3.0): a paradigm in which unvisited worlds are explored by a single or a few missions with radically adaptive space systems. A PE 3.0 mission conducts both initial exploratory science and follow-on hypothesis-driven science based on its own in situ data returns, evolving spacecraft capabilities to work resiliently in previously unseen environments. The key enabler of PE 3.0 is software-defined space systems (SDSSs) - systems that can adapt their functions at all levels through software updates. This paper presents findings from a Keck Institute for Space Studies (KISS) workshop on PE 3.0, covering: (1) PE 3.0 systems engineering including science definition, architecture, design methods, and verification & validation; (2) software-defined space system technologies including reconfigurable hardware, multi-functionality, and modularity; (3) onboard intelligence including autonomous science, navigation, controls, and embodied AI; and (4) three PE 3.0 mission concepts: a Neptune/Triton smart flyby, an ocean world explorer, and an Oort cloud reconnaissance mission.

The Lunar Surface Electromagnetics Experiment (LuSEE-Night) is a joint NASA-DOE-ESA low-frequency radio telescope that will reach the lunar far side in 2027. The unknown dielectric properties of the subsurface at the LuSEE-Night landing site impose the most significant limitation for precision instrument calibration, as reflections from the lunar subsurface can change the primary beam at the 10-20% level. Simulations of these effects have provided insight and concern, showing that the lunar subsurface modeled as a lossy dielectric can absorb a large amount of the power of the sky signal. While this absorption may not strongly impact the signal-to-noise ratio in a sky-noise-dominated regime, it could complicate the beam pattern and make the signal more difficult to model and interpret. We have simulated the far-field properties of the LuSEE-Night beam for varying dielectric profiles of the lunar subsurface. We find that varying the properties of the lunar subsurface has the most significant impact around the antenna resonance, impacting its amplitude, position and width. Conversely, changing the properties of the foreground impacts the data across the band. We use a Bayesian inference pipeline to jointly estimate parameters of a galactic foreground model and dielectric properties of the lunar subsurface around the LuSEE-Night landing site and find that parameters of both the galaxy and subsurface properties can be estimated jointly. While the modeling is somewhat idealized, we believe that the results are largely robust owing to the fact that spectral variations for plausible subsurface and galaxy models have very different spectral signatures.

We report the closed-form expression for Hill's surfaces in the circular restricted three-body problem. The solution $\phi(r,\theta)$, derived in the primary-centric spherical coordinate system, is deduced from a cubic equation delivering at most two roots on each side of a separatrix. The famous patterns (tadpole, horseshoe and peanut shapes, Roche lobes and Hill's quasi-spheres) are exactly produced.

Context: ASPIICS is a giant-baseline visible light solar coronagraph, which relies on the millimetric positioning performance of the precision formation flying Proba-3 mission of the European Space Agency. Proba-3 was launched on 5 Dec 2024, and since then ASPIICS observes the solar corona with the field of view (1.1-3) R_sun. Aims: Diffraction, in particular diffraction of solar disk light on the external occulter, is known to provide a major source of straylight in coronagraphs. We aim to analyze diffracted light visible in ASPIICS images, compare it with the analytical-numerical diffraction model reported earlier, and fine-tune the model. Methods: We compare diffraction effects visible in ASPIICS data with simulated diffraction images; in particular, we compare the geometrical properties and the radiometric signal. The properties of the diffraction described in previous works suggest how to fine-tune the model in order to achieve a better correspondence with the observations. Results: Early ASPIICS observations, where diffraction is pronounced, fully confirm all the qualitative properties of diffracted light suggested by the model. After fine-tuning of the model we see quantitative correspondence of the level of 30\% – 50\%, depending on the configuration. Conclusions: The performed analysis allows (a) to validate our analytical-numerical model and justify the assumptions, and (b) to estimate the contribution of the diffracted light in the ASPIICS images. In the majority of the field of view the diffracted light is two orders of magnitude below the coronal signal.

We have succeeded in operating a transition-edge sensor (TES) spectrometer and evaluating its performance at the SPring-8 synchrotron X-ray light source. The TES spectrometer consists of a 240 pixel National Institute of Standards and Technology (NIST) TES system, and 220 pixels are operated simultaneously with an energy resolution of $4$ eV at 6 keV at a rate of about 1 c/s/pixel. The tolerance for high count rates is evaluated in terms of energy resolution and live time fraction, leading to an empirical compromise of about 2 x 10^3 c/s/all pixels with an energy resolution of 5 eV at 6 keV. By utilizing the TES's wide-band spectroscopic capability, simultaneous multi-element analysis is demonstrated for a standard sample. We conducted X-ray absorption near-edge structure (XANES) analysis in fluorescence mode using the TES spectrometer. The excellent energy resolution of the TES enabled us to detect weak fluorescence lines from dilute samples and trace elements that have previously been difficult to resolve due to the nearly overlapping emission lines of other dominant elements. The neighboring lines of As K alpha and Pb L alpha2 of the standard sample were clearly resolved and the XANES of Pb L alpha2 was obtained. Moreover, the X-ray spectrum from the small amount of Fe in aerosols was distinguished from the spectrum of a blank target, which helps us to understand the targets and the environment. These results are the first important step for the application of high resolution TES-based spectroscopy at hard X-ray synchrotron facilities.

Systematic errors in the parameter estimation (PE) of gravitational wave (GW) mergers can arise from various sources, including waveform systematics, noise mischaracterization, data analysis artifacts, and other unknown factors. In this study, we analyze selected events from the first three observing runs of the LIGO-Virgo-KAGRA (LVK) collaboration. We choose events that have been flagged in various studies as potentially affected by systematic errors. Here, we reanalyze these events using a couple of parametric models developed in previous work that incorporate uncertainties in both the phase and amplitude of the GW waveform. In this data-driven approach, we apply sufficiently broad priors on the uncertainty parameters to account for potential systematic errors. Our findings show that the proposed method effectively reduces systematic errors, even those arising from data artifacts, such as glitches occurring near a signal and the deglitching process in GW frame files. Similarly, inconsistent results from different waveform models become much more consistent in our framework. One noteworthy event we examine is GW191109\_010717, which is particularly interesting due to its anti-aligned spin properties. We report that, within our framework, the event still exhibits anti-aligned spin characteristics, but the inference results become consistent across raw and deglitched frame files, as well as across the waveform models used for this event (IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4). A similar trend is observed for the event GW200129\_065458, which previously yielded a high, but inconsistent precession parameter among different waveform models. In contrast, we observe a non-zero and consistent value of $\chi_{p}=0.60^{+0.31}_{-0.33}, 0.58^{+0.30}_{-0.29}$ and $0.56^{+0.31}_{-0.28}$ for the IMRPhenomXPHM, IMRPhenomXO4a, and NRSur7dq4 waveform models, respectively.

Ringdown signals from perturbed black holes (BHs) offer a clean window into BH spacetime, strong-field gravity, and fundamental physics. Presently the quasi-normal modes of stellar-mass BH ringdowns have been successfully extracted in the ground-based gravitational wave (GW) observations. Looking ahead, the future space-borne observatories will listen to the ringdowns from massive BH binary coalescences more loudly and resolve multiple modes to unprecedented precision, which calls for efficient approaches to mitigate the sharply increasing computational burden. We develop a practical ringdown analysis pipeline for space-borne detectors by implementing FIREFLY, a novel acceleration algorithm validated in ground-based detectors, and for the first time demonstrate its compatibility and effectiveness with the time-delay interferometry (TDI) observables. With high fidelity, we achieve a $\sim 200$-fold speedup for a simulated ringdown signal including six modes, providing a viable and scalable route for multi-mode ringdown analysis in the space context. This new approach has sound statistical interpretation and is extensible to other GW sources in band.

A major challenge in extracting information from current and upcoming surveys of cosmological Large-Scale Structure (LSS) is the limited availability of computationally expensive high-fidelity simulations. We introduce calibrated Neural Quantile Estimation (NQE), a new Simulation-Based Inference (SBI) method that leverages a large number of approximate simulations for training and a small number of high-fidelity simulations for calibration. This approach guarantees an unbiased posterior regardless of approximate simulation accuracy, while achieving near-optimal constraining power when the approximate simulations are reasonably accurate. As a proof of concept, we demonstrate that cosmological parameters can be inferred at field level from projected 2-dim dark matter density maps up to $k_{\rm max}\sim1.5\,h$/Mpc at $z=0$ by training on $\sim10^4$ Particle-Mesh (PM) simulations with transfer function correction and calibrating with $\sim10^2$ Particle-Particle (PP) simulations. The calibrated posteriors closely match those obtained by directly training on $\sim10^4$ expensive PP simulations, but at a fraction of the computational cost. Our method offers a practical and scalable framework for SBI of cosmological LSS, enabling precise inference across vast volumes and down to small scales.

The approaches to searching for axion-like signals based on pulsars include observations with pulsar timing arrays (PTAs) and pulsar polarization arrays (PPAs). However, these methods are limited by observational uncertainties arising from multiple unknown and periodic physical effects, which substantially complicate subsequent data analysis. To mitigate these issues and improve data fidelity, we propose the Artificial Pulsar Polarization Arrays (APPA): a satellite network comprising multiple pulsed signal transmitters and a dedicated receiver satellite. To constrain the axion-photon coupling parameter $g_{a\gamma}$, we generate simulated observations using Monte Carlo methods and investigate the sensitivity of APPA using two complementary approaches: Likelihood analysis and frequentist analysis. Simulations indicate that for the axion mass range of $10^{-22}-10^{-18}$ eV, APPA yields a tighter upper limit on $g_{a\gamma}$ (at the 95\% C.L.) than conventional ground-based observations, while also achieving superior detection sensitivity. Moreover, a larger spatial distribution scale of the satellite network corresponds to a greater advantage in detecting axions with lighter masses.

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

Due to the limited generalization and interpretability of deep learning classifiers, The final vetting of rare celestial object candidates still relies on expert visual inspection–a manually intensive process. In this process, astronomers leverage specialized tools to analyze spectra and construct reliable catalogs. However, this practice has become the primary bottleneck, as it is fundamentally incapable of scaling with the data deluge from modern spectroscopic surveys. To bridge this gap, we propose Spec-o3, a tool-augmented vision-language agent that performs astronomer-aligned spectral inspection via interleaved multimodal chain-of-thought reasoning. Spec-o3 is trained with a two-stage post-training recipe: cold-start supervised fine-tuning on expert inspection trajectories followed by outcome-based reinforcement learning on rare-type verification tasks. Evaluated on five rare-object identification tasks from LAMOST, Spec-o3 establishes a new State-of-the-Art, boosting the macro-F1 score from 28.3 to 76.5 with a 7B parameter base model and outperforming both proprietary VLMs and specialized deep models. Crucially, the agent demonstrates strong generalization to unseen inspection tasks across survey shifts (from LAMOST to SDSS/DESI). Expert evaluations confirm that its reasoning traces are coherent and physically consistent, supporting transparent and trustworthy decision-making. Code, data, and models are available at this https URL .