Papers for Tuesday, Mar 03 2026

\textit{PeVatrons} are the extreme galactic accelerators capable of producing PeV particles. Recent observation of Large High Altitude Air Shower Observatory have detected UHE photons ($\geq$ 100 TeV) from 43 galactic sources. Detection of UHE photons demands the presence of at least PeV particles in the acceleration site. Although the exact nature of most of the sources are still unknown, a large fraction of these sources have spatial association with pulsar wind nebula. In this work we investigate the acceleration mechanism in pulsar wind nebula by following a magnetohydrodynamics approach. Current study relates the MHD flow solution in immediate downstream with the particle spectrum and spectral energy distribution of photons. Our study shows that MHD description in the PWN environment reduces the parameter space and most of the parameters can be constrained in terms of a single parameter, \textit{the magnetization parameter} $\sigma$ only. Considering the effect of $\sigma$, we show that in low $\sigma$ environment pulsar wind nebula can produce PeV particles. We have also investigate the role of turbulence in the nebular region in acceleration of particle to PeV energy. Current study shows that both low $\sigma$ environment and turbulence environment is favorable for acceleration of particles up to PeV energy. We have also tested our model in two different LHAASO detected PeVatron 1LHAASO J1848-000u and 1LHAASO J1929+1846u.

A range of speculative space ventures envisage the use of ultra-large structures for the collection and reflection of light. Given the length-scale of such structures they cannot be considered as point masses for the calculation of gravitational and radiation pressure forces. Using a simplified model it will be demonstrated that ultra-large reflectors in static equilibrium levitating above a central star (so-called stellar engines) are always unstable if the reflector comprises a uniform disc. However, if the reflector has a non-uniform mass distribution, specifically a ring supporting a reflect or, a stellar engine can in principle be passively stable. Moreover, while it can be shown that static swarms of reflectors levitating above a central star (so-called Dyson bubbles) are unstable, in principle they can become passively self-stabilizing if arranged about the star as a dense cloud. While such ventures are clearly speculative, understanding the orbital dynamics of ultra-large structures, and in particular the conditions for passive stability, can provide insights into the properties of potential technosignatures in search for extraterrestrial intelligence studies.

The last evolutionary stages of massive black hole binaries prior to coalescence is dominated by the emission of gravitational waves, which will be probed by the future Laser Interferometer Space Antenna. If gas is present around the two black holes, however, the associated electromagnetic emission can provide additional information about the binary properties and location before the merger event. For this reason, a proper characterisation of the electromagnetic emission during these phases is of fundamental importance, and requires a detailed description of the gas dynamics close to the event horizon of the two black holes, only achievable via numerical simulations. Within this context, we present the implementation of the Superposed Kerr-Schild dynamic metric in the relativistic scheme in the meshless code GIZMO. Our code can now simulate black hole binaries approaching merger with high computational efficiency and accuracy, taking into account relativistic effects on the gas. To validate our implementation, we perform two tests. First, we explore the case of a relativistic Bondi flow around a binary, finding very good agreement with numerical relativity simulations. Then we explore the case of an inviscid relativistic circumbinary disc, comparing our results with a similar simulation run assuming Newtonian gravity. In this second case, we find moderate differences in the mass accretion rate and in the inflow dynamics, which suggest that the presence of a non-Keplerian potential and of apsidal precession in the orbiting gas trajectories may produce stronger shocks and boost angular momentum transport in the disc. Our work highlights the importance of accounting for relativistic corrections in accretion disc simulations around black hole binaries approaching merger, even at scales much larger than those currently probed by numerical relativity simulations.

We investigate the origins of quasi-periodic eruptions (QPEs) in galactic nuclei using global three-dimensional meshless finite-mass (MFM) simulations. By modeling stellar and black-hole impactors traversing accretion disks under various inclinations and surface densities, we evaluate their consistency with the observed properties of QPEs. Stellar impacts produce highly asymmetric bipolar ejecta with forward outbursts dominating by over an order of magnitude in energy and luminosity due to the star blocking downstream flow and creating a low-density wake. This shock-compression mechanism often renders backward events unobservable, implying one detectable burst per orbit, and challenging the standard assumption of two bursts. It also fails to explain alternating long-short recurrence patterns and places several sources near or within twice the tidal disruption radius for solar-mass stars, raising severe stability concerns. Whereas a stellar-mass black hole (sBH) gravitationally focuses and heats disk gas extending from its Bondi radius $R_{\rm B}$ to its Hill radius $R_{\rm H}$ during an impact, yielding nearly symmetric ejecta with mild contrasts. This gravitational-drag mechanism generates higher energy budgets at low inclinations due to enhanced mass accumulation. We suggest an ad hoc effective interaction radius $ R_{\rm eff} \simeq 0.5\, R_{\rm B}^{1/3} R_{\rm H}^{2/3} $ to quantify this trend. Our semi-analytical model confirms that sBH-disk collisions can power the full QPE energy range ($10^{44}$-$10^{48}$ erg), naturally accounting for periodicity, asymmetry, durations and diversity.

We report the results of a radial velocity monitoring program in the Hyades region, carried out at the Center for Astrophysics over a period of more than 45 yr. Nearly 12,000 spectra were gathered for 625 stars brighter than $V \approx 14.5$, of which 55% are members or possible members of the cluster. New or updated spectroscopic orbital solutions are presented for more than 100 members and non-members, including several triple systems. In a few cases we incorporate available astrometry. The frequency of binaries in the Hyades with periods up to $10^4$ days is determined to be $40 \pm 5$%, after corrections for incompleteness. This is marginally higher than in other open clusters. The orbital period and eccentricity distributions are found to be similar to those of solar-type binaries in the field. The mass ratio distribution is essentially flat, or slightly rising toward mass ratios of unity. We revisit the determination of the tidal circularization period, obtaining a longer $P_{\rm circ}$ value of $5.9 \pm 1.1$ days compared to the previous estimate of 3.2 days, still somewhat short of the value expected if most or all of the action of tides happens during the pre-main-sequence phase. We estimate a line-of-sight velocity dispersion of $0.21 \pm 0.05$ km s$^{-1}$ within 5.5 pc of the cluster center (approximately the half-mass radius) and a larger dispersion beyond that distance. Our velocity measurements are accurate enough to clearly reveal the signatures of gravitational redshift and convective blueshift among the dwarfs and giants in the Hyades.

We provide a complete census of the baryons in early galaxies to investigate the phases in which gas and stars reside, their corresponding budgets, depletion times, and stellar return fraction as a function of redshift and stellar age. We use the ColdSIM hydrodynamical time-dependent non-equilibrium chemistry simulations and perform a detailed analysis of the cold, warm, hot, and stellar phases for both bound structures (galaxies/CGM) and the diffuse IGM. We investigate in depth the cold HI and H2 components, explicitly computed in our simulations, and their relations with host mass, SFR, metallicity and depletion times. We also provide observational insights and discuss the implications for stellar mass functions, PopIII star formation and changes in the IMF. We find that cosmic gas prior to reionisation is mostly cold, while at later epochs the warm phase becomes dominant due to enhanced star formation activity and increasing UV reionising radiation. Stellar return fractions at these times are ~0.15-0.20, a factor of two lower than the values usually adopted. Cold, warm, and hot gas masses as well as HI and H2 components show increasing trends with mass and SFR, while depletion times decrease down to 0.01-0.1 Gyr with a weak metallicity dependence. The resulting star formation efficiency remains at the level of a few per cent and gas-to-star fractions decline with mass, influenced by local feedback and environment. Our findings are consistent with ALMA, VLA and IRAM surveys at later epochs, including ALFALFA, xCOLDGASS, GASS, xGASS, EDGE-CALIFA, PHIBBS, and ASPECS. Gas phases are quantitatively related to the underlying stellar populations and can be used to infer unknown quantities. In the appendix we provide fit functions describing the trends of the stellar return fraction, the main sequence, phase mass relations, gas-to-star fractions and depletion times.

We calibrate a physically motivated, super-Eddington eruptive mass-loss prescription for red supergiants (RSGs) using Local Group stellar populations. Building on MESA models that add eruptive mass loss with a free scaling parameter $\xi$, we generate stellar evolution tracks and isochrones, and synthesize mock populations at metallicities of $Z/Z_\odot=0.2,\ 0.4$, and $1.0$. We compare model luminosity functions to observations of RSGs in the SMC, LMC, and M31, restricting to $3.5<\log T_{\rm eff}/K<3.75$ and $\log(L/L_\odot)>4.5$. By-eye fits to the observations yield values of $\xi_\mathrm{SMC}=0.0-0.05$, $\xi_\mathrm{LMC}=0.1$, and $\xi_\mathrm{M31}=0.35$, implying a positive, linear trend between the strength of eruptive mass-loss and metallicity. This calibrated eruptive mass loss prevents stars with initial masses $\gtrsim 20~M_\odot$ from evolving to become red supergiants, with implications for the mass spectrum of core-collapse progenitors, compact remnants, early supernova interaction signatures, and the spectral energy distributions of unresolved galaxies.

Recent observations revealed that distant galaxies have bursty star formation histories, regulated by stellar or active galactic nuclei (AGN) feedback and gas inflows. According to theoretical models, feedback preferentially removes metal-rich gas, while subsequent starbursts are triggered by mergers and newly-accreted gas that is generally less enriched than the galaxy interstellar medium (ISM). Therefore, gas-phase metallicity holds key insights into the baryonic processes shaping early galaxies. We present the first NIRSpec/IFU study of spatially resolved ISM properties in the MACS0647-JD system ($z=10.17$). The system consists of two stellar components detected in NIRSpec/IFU and NIRCam photometry. The main component ($\log (M_{\ast}/M_{\odot}) =7.77 \pm 0.09$; $12+ \log(O/H)=7.89 \pm 0.11$) is more massive and significantly more metal-rich compared to its companion ($\log (M_{\ast}/M_{\odot}) =7.42 \pm 0.07$; $12+ \log(O/H)=7.47 \pm 0.14$), suggesting an older stellar population and a prolonged chemical enrichment history. We find that the H$\gamma$ line emission centroid is offset by $0.1^{\prime \prime}$ (150 pc in the source plane) from the stellar continuum centroid; the latter coincides with the location of the main stellar component. This offset provides possible evidence of a merger-driven starburst in this system. By comparing the spatial distribution of the metallicity, velocity dispersion, and the burstiness of star formation history, we infer the presence of turbulent, metal-poor gas outside the stellar components. This metal-poor, dynamically unstable gas is likely responsible for the increase in the recent star formation in the north-east region of the system.

We present the results of a multi-wavelength study of a population of X-ray bright ($\rm log(F_{0.2-12 \ keV})>-12.5$), non-$\gamma$-ray detected high and extreme high synchrotron peak (HSP, EHSP; $\rm log(\nu_{\rm peak,\ Hz})>16$) BL Lacs to $i$) put stronger constraints on the synchrotron peak location and shape and $ii$) model their expected behaviour in the very high-energy band. First, we performed an X-ray spectral analysis, using XMM-Newton, Chandra, Swift-XRT, and eROSITA data, and fitting the spectra using both a power law and a log parabola model. Out of 78 sources in the initial sample, 17 were best described by a log parabola model, a result that supports a scenario where the synchrotron peak falls in the X-ray band. Among these 17 sources, we further selected the 10 objects dominated by the jet emission, with no significant contamination of the host galaxy. We performed a $\gamma$-ray analysis of \lat\ data for these objects, obtaining upper limits providing information on their flux in the 100 MeV - 300 GeV energy range. We then modelled the broadband SED of these objects with JetSeT using two models: one assuming a log parabola for the electron distribution and the other one with a broken power law electron distribution, using parameters consistent with those describing the emission of the prototypical EHSP 1ES 0229+200. We found the models to be generally consistent with the available multi-wavelength detections and upper limits. Furthermore, they confirmed that a subsample of sources could display relevant emission in the TeV energy range, even potentially reaching the threshold for detectability by the Cherenkov Telescope Array Observatory.

We quantify the impact of finite catalog size, or "catalog variance," on current gravitational-wave population analyses. The distribution of merging binary black holes is commonly reconstructed via hierarchical Bayesian inference, with uncertainties reported as credible intervals. Such intervals are conditioned on the specific realization of the observed events and are therefore themselves subject to variability arising from the finite size of the catalog. We estimate this "uncertainty on the uncertainty" using statistical bootstrapping applied to data segments containing both detected events and sensitivity injections. Applying this framework to GWTC-4, we find that the inferred population distributions exhibit substantially broader uncertainties than those obtained in a standard single-catalog analysis. In particular, the $\sim 35\,M_\odot$ peak in the primary-mass distribution is largely absorbed by statistical fluctuations once catalog variance is taken into account. Unlike other studies that rely on simulating catalogs by assuming an underlying population, this work provides the first data-driven assessment of the uncertainty intrinsic to the observed gravitational-wave catalog. Accounting for catalog variance is important for drawing robust astrophysical conclusions from gravitational-wave data, avoiding inferences driven by a particular finite realization rather than genuine population features.

Hot subdwarf variables (sdBVs) that display large-amplitude ($>$1%), short-period variability, as a result of radial-mode pulsations, have recently become objects of interest as they show unique properties among the sdBV classes. Since the discovery of objects such as Balloon 090100001 and CS 1246, twelve more have been discovered in the Zwicky Transient Facility (ZTF) survey that display similar characteristics. However, due to lack of broad spectroscopic investigations, it remains unclear whether these objects constitute a distinct class of radial-mode dominant sdBVs that share common atmospheric and fundamental properties. Here we aim to spectroscopically define these peculiar sdBVs as a population. We collected low-resolution spectroscopy on a sample of sdBVs discovered in the ZTF survey, including time-series observations. We fitted the spectra to a grid of theoretical models to determine their mean effective temperature, surface gravity and helium abundance and any corresponding variability. We then use these properties to estimate the mass, radius and luminosity using a spectral energy distribution fitting method. We show that the resulting properties are similar to the radial-mode dominant sdBVs, Balloon 090100001 and CS 1246, and that they are distinguishable from other similar radial-mode pulsators, such as blue large-amplitude pulsators. We find that these stars, on average, have mean effective temperatures of 28,300 K and surface gravity measurements of $\log\,g=5.56$, with changes in these parameters on the order of 1000 K and 0.10 dex, respectively. The location of these stars on the $T_{\textrm{eff}}$ -- $\log\,g$ plane places them on the boundary region between the low-amplitude, multi-periodic V361 Hya and V1093 Her stars, where the hybrid DW Lyn pulsators lie. The masses and radii of the majority of the sdBVs in our sample align with canonical-mass sdB properties.

Dwarf galaxies are ideal laboratories to probe the interplay between galaxy formation and the growth of black holes (BHs) in the early Universe. Mounting observational evidence reveals the presence of BHs in low-mass galaxies across cosmic time, with $\textit{JWST}$ uncovering a likely population of $\textit{overmassive}$ BHs at $2 \lesssim z \lesssim 11$. Simulations struggle to reproduce this high-redshift regime, motivating revisions to models of BH accretion and feedback from active galactic nuclei (AGN). To address this, we present high-resolution cosmological zoom-in simulations of a dwarf galaxy based on FABLE physics, introducing novel sink-based BH accretion models and relaxing the fiducial assumption of strong supernova feedback. BHs accrete more efficiently in the sink-based runs compared to the `traditional' Bondi-based counterparts, with AGN feedback leading to early, rapid quenching maintained by fast, hot and metal-enriched outflows. These outflows pollute the outer circumgalactic medium, yielding flat metallicity gradients down to $z=0$. We further assess the performance of two widely used virial estimators and find significant departures from the true dynamical mass, especially during the high-redshift dwarf assembly. Since our galaxy is dark-matter-dominated at all times and radii, BH growth, tied to the baryon cycle, shows no clear correlation with global dynamical properties. Efficient AGN feedback, produced by overmassive BHs relative to extrapolated local $M_\bullet - M_\star$ relations, indicates that dormant BHs residing in local, quenched dwarfs might be the relics of some of the high-redshift $\textit{JWST}$ BHs.

We measure the star formation history (SFH) of Local Group dwarf galaxy WLM using wide-area ($\sim4$ half-light radii) ground-based NIR imaging of bright ($M_{J}<-4.9$ mag) AGB stars. From our NIR CMD of 825 stars, we find that our recovered SFH is in excellent agreement with literature SFHs of WLM measured from much deeper CMDs ($M_{F090W}\sim+4.3$ mag) based on JWST imaging. We find good agreement in the qualitative shape of the SFHs as well as quantitative metrics such as the timescales for which 50% and 90% of the stellar mass formed with $\tau_{50, {\rm AGB}}=5.16_{-0.50}^{+2.07}$ Gyr ago and $\tau_{90, {\rm AGB}}=1.33_{-0.09}^{+0.11}$ Gyr ago versus $\tau_{50, {\rm JWST}}=5.29_{-0.28}^{+0.34}$ Gyr ago and $\tau_{90, {\rm JWST}}=1.42_{-0.01}^{+0.16}$ Gyr ago. The coarser precision of the AGB star-based values is driven by the low number of AGB stars. We also recover an age gradient that is in good agreement with the age gradient measured from JWST data, where we find the outer regions of WLM are on average older than the inner regions. We derive an age-metallicity relation (AMR) from the AGB star CMD fitting that is similar to the JWST-based AMR and is consistent with other reported metallicities in WLM (i.e., spectroscopy, RR Lyrae). From our wide-area AGB star map and SFH, we identify a stellar over-density ($M_*\sim2.0\times10^6~M_{\odot}$, $r_h\sim340$ pc) in WLM's northwestern outer disk. The over-density's SFH shows a burst of star formation $\sim8$ Gyr ago and its spatial location is near a known warp in WLM's H I. Despite WLM having long been considered an isolated galaxy, the mass, size, and age of this over-density are highly suggestive of an accreted dwarf galaxy. Overall, our findings illustrate the power of NIR observations of AGB stars for efficiently and accurately measuring SFHs and for identifying and characterizing substructures in nearby galaxies.

On the scale of the cosmic horizon, signatures that are unique to general relativity are concealed within the statistics of the large scale distribution of galaxies. These were thought to be beyond the reach of all but the most ambitious galaxy surveys, as they are substantially suppressed relative to standard redshift-space distortions. We show that the detectability of these higher-order relativistic effects can be dramatically enhanced by a sampling strategy that splits a galaxy catalogue into faint and bright subsamples and then combines their auto-bispectra. For current surveys such as DESI, this implies that this new signal will be detectable for the first time using our new strategy.

Thermal dust continuum polarimetry is a powerful indirect probe of magnetic field geometry in dense molecular clouds while at the same time providing information on the alignment of dust grains with the magnetic field. The leading theory of grain alignment, Radiative Torque Alignment (RAT), has been successful in explaining a variety of observations, including the loss of polarization fraction toward high column densities. One prediction of RAT is that an increase in grain alignment efficiency should be observed in the environments surrounding protostars, due to radiation from the embedded source. However, observational confirmation of this prediction remains scarce. In this study, we sought to test the theoretical prediction of enhanced grain alignment near protostars in the high-mass star forming region DR21 using 214 $\mu m$ SOFIA/HAWC+ observations. We investigated the correlation of the polarization fraction of dust emission, $p$, and the polarization angular dispersion, $S$, with respect to total intensity. We also probed intrinsic dust polarization properties using the product $S\times p$ as a proxy. We detected significant polarization fractions even at the highest intensities, where strong depolarization is typically expected. The polarization fraction-intensity trend flattens at $I > 1.6^{+0.3}_{-0.3} \times 10^4$ MJy/sr ($N_{H_2}$ $\sim 2\times 10^{23}$ ${cm}^{-2}$). We compared the observed trends with predictions from an analytical model of a centrally heated envelope surrounding an embedded luminous protostar. The predictions from the simple model agree well with the observed trends. Our results provide strong support for enhancement of grain alignment by local radiation from embedded sources.

We present new Hubble Space Telescope (HST) imaging of ionised filaments in the brightest group galaxy NGC 5044. These filaments extend several kiloparsecs and have widths of $\sim$50--120 pc, with some as narrow as those in cluster cores and others broader, reflecting the lower confining pressure in groups. Filament width ($W$) scales with ambient pressure ($P$) as $W \propto P^{-0.4}$. Combining HST, ALMA, and MUSE data, we measure column densities and magnetic field strengths. Equipartition fields decline from $\sim$40 $\mu$G at the centre to $\sim$20 $\mu$G at 5 kpc, about 2--3 times weaker than in clusters. Dynamical stability requires stronger radial fields ($\sim$10$^2$ $\mu$G), consistent with simulations and magnetic draping, though such high values exceed Faraday Rotation Measure limits. Turbulence and cosmic rays also contribute support. Group and cluster filaments are stable against gravitational collapse, and ultraviolet imaging reveals no star formation in NGC 5044 ($<$10$^{-3}$ M$_\odot$ yr$^{-1}$). NGC 5044 hosts an ionised gas core within its Bondi radius with $n_e \propto r^{-1}$ and filling factor $f \gtrsim 3 \times 10^{-3}$, that is connected to the extended filaments, suggesting a channel for gas inflow toward the black hole. Group and cluster filaments likely share a common origin, with magnetic fields and AGN feedback preserving their structure. Ambient pressure and dust survival regulate molecular gas formation. Lower-pressure groups favour broader, more diffuse filaments with sporadic molecular clumps and weaker dust shielding, whereas higher-pressure clusters host narrower strands with stronger molecular-ionised gas alignment. We predict that (i) filament width scales with ambient pressure, (ii) filament-coincident Faraday rotation structures emerge at $\leq 0.1$ kpc resolution, and (iii) molecular/ionised gas co-spatiality is weaker in groups than in clusters.

The Herbig Ae star AB Aurigae hosts a vast, low-inclination protoplanetary disk that exhibits a plethora of substructures, including the protoplanet candidate AB Aur b. We present M-band spectroscopic data taken with NASA IRTF from Feb 2024 covering multiple position angles that captured emission from an off-centered, low temperature, and compact source. Analysis of the ${}^{12}$CO $\nu=$1-0 low-J ro-vibrational emission line profiles and spectroastrometric signals localizes the source at around an orbital radius of 65 au and a position angle of 143$^\circ$. These coordinates are distinctly different from those of AB Aur b, which was not detected. Although there is no obvious explanation for the detected source, if we assume it was a circumplanetary disk, then its maximum temperature would be about 550 K and its maximum radius would be about 5 au. Our results alludes to a previously unknown companion that may be residing in the AB Aurigae system.

We use controlled $N$-body experiments to study the dynamical evolution of dwarf spheroidal galaxies (dSphs) embedded in dark-matter (DM) haloes containing a large population of dark subhaloes. We show that stellar orbits subject to stochastic force fluctuations irreversibly gain energy and expand toward a dynamical attractor characterized by a stellar half-light radius $r_{\rm half} \approx r_{\rm max}$ and a velocity dispersion $\sigma \approx 0.5\,v_{\rm max}$, where $v_{\rm max}$ is the peak circular velocity of the host halo at radius $r_{\rm max}$. This state is reached both in isolation and under tidal stripping, although tidal mass loss significantly accelerates the evolution. Assuming that the Milky Way (MW) dSphs have reached this state, we find that the inferred halo masses collapse onto narrow sequences as a function of $r_{\rm half}$. Under this assumption, MW satellites with $r_{\rm half} \lesssim 1\,\mathrm{kpc}$ follow the tidal tracks of cuspy haloes, while larger systems deviate in a manner consistent with cored DM profiles. Moreover, the mass--luminosity relation follows the slope expected from abundance matching, but with halo masses systematically lowered from their peak values at fixed luminosity. These results suggest that the structural diversity of dSphs is largely an evolutionary outcome driven by internal heating and tides, rather than by the conditions of star formation. This framework predicts that isolated, early-quenched dSphs should have systematically larger sizes than satellites, a prediction testable with upcoming surveys.

As part of its commissioning, the Vera C. Rubin Observatory observed several fields repeatedly for a month with ComCam, an instrument that uses the same hardware as the LSST camera but covers a smaller field of view. We photometrically classify AT 2024ahzi, a transient discovered by ComCam, as a Type IIP supernova (SN IIP) using both ComCam and DECam photometry. We find that the duration, luminosity, and color of AT 2024ahzi's photometric plateau are all consistent with those from a large sample of SNe II. By comparing its multi-band light curves to SN II models and analytic relations, we place constraints on the SN progenitor, explosion dynamics, and circumstellar environment. We argue that the progenitor has an extended density profile indistinguishable from a slowly accelerating CSM. We discuss how a similar workflow can identify and characterize future Rubin SNe II.

The cold gas content of post-starburst galaxies (PSBs) provides important insight into the mechanisms that drive rapid quenching, but a multiphase assessment of both the atomic and molecular gas in PSBs does not yet exist. We introduce the Ensemble of Multiphase Baryons Evolving in Rapidly-quenching Systems, or EMBERS, a homogeneously selected, nearly mass- and redshift-complete survey of the global atomic (HI) and molecular gas (H2) in PSBs, observed with the Five Hundred-metre Aperture Spherical Telescope (FAST) and the Institut de radioastronomie millimetrique (IRAM) 30m telescope. We present new CO(1-0) observations for 52 PSBs with the IRAM 30m, which, combined with 9 archival observations, gives a total H2 sample of 61, of which 58/61 have ancillary HI measurements. We detect CO(1-0) in 34/61 galaxies, corresponding to molecular gas fractions (fH2 = MH2/M*) ranging from two to 250 per cent. By comparing with a stellar-mass matched star-forming (SF) control sample from xCOLD GASS, we find that PSBs on average are 0.3-0.6 dex depleted in H2. However, considering both HI and H2, individual PSBs host diverse gas reservoirs ranging from gas-rich in both phases, elevated in one phase, or gas-poor, the latter of which is common at lower stellar mass. The existence of gas-normal and gas-depleted PSBs in both phases suggests that some PSBs may rejuvenate their star formation, but the rapid shutdown of star formation in others is likely terminal. Despite this diversity, the majority of EMBERS PSBs are gas-poor compared to SF controls, with the typical PSB hosting gas reservoirs intermediate to those found in star-forming and quenched galaxies.

We have conducted a search for radio pulsars in six Galactic stellar binary systems having unseen primary stars. All six systems have estimated primary masses in the range that could be consistent with neutron stars. We used the Green Bank Telescope at a center frequency of 350 MHz to search for dispersed periodicities and single pulses across a range of possible dispersion measures (DMs) and binary accelerations. No astrophysical signals were detected in our search. The estimated 400-MHz luminosity upper limits from the search are comparable to or smaller than the lowest values observed for almost all the known Galactic binary pulsars having cataloged 400 MHz radio luminosities. This implies that the systems we observed either do not harbor radio-emitting pulsars, contain pulsars that do not beam in our direction, or contain pulsars with luminosities that are significantly lower than this subset of the known Galactic binary pulsar population.

[...] We observed XTE J1810-197 between 29 September 2022 and 14 July 2023 with the radio telescopes at the Argentine Institute of Radioastronomy (IAR). We searched for single pulses in time series at a DM range of 100-400 pc cm-3 , with a threshold in signal-to-noise ratio (S/N) of 8. [...] We found 249 giant pulses at a DM mean value of 178.8$\pm$0.1 pc cm-3 . We measured peak flux densities up to 119 Jy, and fluences up to 58 Jy ms. We fitted a power law distribution to the flux density, obtaining an index of -4.0$\pm$0.3. We observed a maximum rate of approximately 15 pulses per hour on 20 February 2023, followed by an abrupt disappearance of transient radio emission, indicating a transition to a less active state. The brightest single pulses are limited to a $\sim$2$\%$ of the rotational phase and have similar fluence values to the reported intermediate FRB-like bursts of SGR 1935+2154. No significant X-ray activity in the MAXI data was detected during the radio observing period. This is the first study of single radio pulses of a magnetar using IAR data, showing the potential of the upgraded telescopes for investigating the transient radio sky. The properties of the single pulses detected here show the magnetar transient nature and capability to emit high-luminosity pulses. We compared the detected emission to FRB-like bursts and single pulses emitted by SGR 1935+2154. Even though the mechanism producing all the events should be coherent, the luminosity of the events, features on the dynamic spectra, and the difference between being phase confined or not, indicate that XTE J1810-197 presents GP emission, while SGR 1935+2154 only shows normal single pulses or FRB bursts. This could indicate that the conditions for producing each type of event differ.

Interacting coronal mass ejections (CMEs) result in complex heliospheric structures that can dramatically enhance their geoeffectiveness compared to isolated events. A striking example of such complex structures is that of the Mothers Day event, which occurred during 10-14 May 2024, leading to the strongest geomagnetic storm in decades. It was driven by at least ten interacting CMEs accompanied by flare and filament eruptions from the solar western hemisphere. We aim to understand its solar and heliospheric origins, variation in plasma and magnetic characteristics, and predict its geoeffectiveness. Furthermore, we focus on bringing to the community the challenges and limitations faced during themodelling of such a complex, yet important environment. We use multi-point remote observations of CMEs and employ different methods to estimate their speeds and geometry. These parameters drive the 3D magnetohydrodynamic simulations of CME evolution and heliospheric propagation in the framework of the EUHFORIA model. The best-performing simulation reproduced the arrival time of the storm and peak intensity with a lead of 2 hours, and estimated the storm strength with about 70% accuracy. This work underscores the requirements and challenges involved in accurately modelling extreme events. The major bottleneck was to constrain the CME input parameters that dictate the accuracy of the data-driven simulations. Predicting the impacts of the complex heliospheric structures during this event required the modelling of the CME-CME interactions accurately. As a community, we need to invest more in observational infrastructure while improving the speed and accuracy of our MHD forecasting models.

The small-scale properties of circumgalactic gas in ordinary galaxies drive its bulk properties: the mass loading of cold neutral gas in galactic outflows affects their bulk momentum; gas cooling processes on small scales affect the spatial distribution of gas in the cool (T~$10^4$K) circumgalactic medium (CGM). However, hydrodynamical simulations have yet to resolve the CGM on such small scales. Spectroscopy remains our primary probe of the small-scale CGM, with which sub-parsec scales are challenging to resolve. Fast radio bursts (FRBs)--microsecond to millisecond duration radio pulses--are temporally broadened ("scattered") by gradients in the electron density transverse to the line of sight, often generated by fluctuations on the smallest spatial scales. This makes FRB scattering a powerful, complementary, and scalable probe of the small-scale CGM. We show that the distribution of scattering timescales introduced by density fluctuations within a single, foreground halo--the tau distribution function, or TDF--is sensitive to the small-scale spatial morphology of the gas. The TDF is readily measurable and is analogous to areal covering factors reported in quasar absorption statistics. We compute the TDF in two regimes: scattering from a turbulent, volume-filling medium ("volumetric scattering") distributed along the line of sight; and scattering from discrete structures localized along the line of sight ("intermittent scattering"). Within these regimes, the TDF is sensitive to whether the cool gas comprises primarily spherical, filamentary (1D), or sheet-like (2D) structures. This work sets the stage for upcoming observations which will use hundreds of sight-lines through nearby halos to probe the small-scale CGM, and points out a novel science case for FRB detectors like MeerKAT, Parkes, FAST, and the DSA-2000, which are exquisitely sensitive over a narrow field of view.

The Askaryan Radio Array (ARA) is a neutrino experiment at the South Pole, designed to detect radio-frequency emissions produced by interactions of ultra-high energy (UHE) neutrinos with the Antarctic ice. The array consists of five autonomous stations, each equipped with deep in-ice antennas sensitive to both vertically and horizontally polarized radio signals. With nearly 30 station-years of livetime accumulated, ARA is now conducting its first comprehensive array-wide search for diffuse UHE neutrinos. This analysis is expected to deliver the most stringent constraints from any in-ice radio-based detector up to 1~ZeV and is capable of probing flux levels suggested by KM3NeT around 220~PeV. The results from this analysis marks a critical step toward establishing scalable techniques for next-generation detectors.

Identifying the host galaxies of fast radio bursts (FRBs), and comparing their redshifts and dispersion measures, has unlocked a new probe of the cosmological distribution of ionised gas. However the necessary optical observations to identify FRB hosts, and measure their redshifts, are becoming increasingly onerous as the detection rate of precisely localised FRBs increases. Here we analyse the ability of the Legacy Survey of Space and Time (LSST), being conducted by the Vera C. Rubin Observatory, to identify FRB host galaxies, and the utility of LSST photometric redshifts for FRB cosmology. By combining a model of FRB host galaxy r-band magnitudes, $m_r$, with predictions for the FRB z-DM distribution, we create a method to predict the $m_r(z)$ distribution for the host galaxies of FRBs detected by radio surveys. We then predict these distributions for the coherent modes of the Australian Square Kilometre Array Pathfinder (ASKAP) and MeerKAT. We find that even a single visit with Rubin will be able to identify 65% of FRB host galaxies detected by ASKAP's coherent upgrade, `CRACO'; while the final 10 year co-added images will identify 81% of those from MeerKAT's tied array beams. We also simulate the impact of using photometric redshifts for a simplified analysis to determine $H_0$, finding that estimated photo-z errors result in a decreased precision of only 7% on $H_0$ for ASKAP's CRACO system. The impact of missing dim FRB hosts, which are likely at higher redshifts, is more significant, and might degrade sensitivity to $H_0$ by 47%, or 62% when combined with photo-z errors. All told, Rubin's LSST will be an incredibly powerful survey for facilitating FRB cosmology, although supplemental observations may be useful for particularly dim and distant host galaxies.

We present the detection and validation of a small, temperate transiting exoplanet orbiting TOI-1080 every 3.9652482$^{+0.0000014}_{-0.0000015}$ days. The host is a quiet M4V star at 25.6 pc. The planet signal was first detected by TESS and validated using TESS and ground-based observations. By fitting the available light curves, the planet radius is measured to be 1.200+- 0.058 Rearth and its equilibrium temperature of 368$^{+12}_{-10}$ K. With NIRPS radial velocities, we are able to place a 3-sigma upper limit on the mass of TOI-1080 b of 10.7 Mearth. Our injection-recovery tests enable us to discard additional transiting planets in the TOI-1080 system with radii down to 0.9 Rearth and periods between 0.5 and 7.7 days, and planets with radii larger than 1.4 Rearth for periods up to 19 days. We demonstrate that it is highly amenable to characterisation of its mass and putative atmosphere. In particular, we find that TOI-1080 b is an exceptional target for the ongoing JWST+HST Rocky Worlds DDT programme, having a priority score that is higher than four out of nine targets currently being investigated by the programme. TOI-1080 b can be added to the sample of nearby benchmark planets accessible for detailed study with JWST.

This study presents findings from narrowband imaging of comet C/2006 P1 (McNaught) using the 3.6-metre New Technology Telescope (NTT) at La Silla, Chile. Observations commenced on January 27, 2007, 15 days after perihelion, and continued until February 4, with additional sessions from February 25 to 28. Imaging was conducted using the ESO Multi-Mode Instrument (EMMI) in both broadband (B, V, R) and six comet-specific narrowband filters (CN, C3, C2, NH2, blue and red continuum). Various image processing techniques were employed to enhance structural features, including azimuthal mean/median division and subtraction, azimuthal renormalisation, and division by inverse profile, as well as the Larson-Sekanina technique. These enhancements revealed dynamic coma structures, with jets transitioning from spiral patterns to linear or fan-like shapes over time. The consistency of morphological patterns across different processing methods validated their authenticity. The periodic recurrence and temporal evolution of CN coma features in narrowband images indicate a nucleus rotation period of 11.8 h, consistent with stable active regions and rotationally modulated outgassing near perihelion.

Quantifying the natal kick distribution of pulsars is essential for understanding supernova physics and binary evolution, yet measurements are historically limited by the lack of radial velocity data. Most previous studies rely on transverse velocities under the assumption of spatial isotropy. In this work, we reconstruct the intrinsic three-dimensional (3D) velocity distribution for a curated sample of 18 pulsars by explicitly incorporating the observational constraint of spin-velocity alignment. Using a hierarchical Bayesian framework that accounts for measurement uncertainties, we compare nine candidate velocity distribution models. We find that a Gamma distribution provides an adequate description of the inferred 3D velocities; however, the modest Bayes factor (1.65 relative to a single Maxwellian) indicates that the current data lack sufficient resolving power to discriminate decisively among the models considered. The Gamma model is characterized by a peak velocity of $237^{+67}_{-84}$ km s$^{-1}$. The reconstructed 3D velocities under alignment are systematically lower than those inferred under isotropy, indicating that projection effects can bias individual kick estimates high, while leaving the overall population scale largely unchanged within uncertainties. A complementary analysis of 465 pulsars with transverse velocity estimates favors a Log-Normal distribution for the full sample, while isolated young pulsars remain consistent with a Gamma-like profile. Our results underscore the importance of geometric assumptions in population inference and highlight the need for larger samples with improved distance and spin-axis measurements to place tighter constraints on natal kick physics.

Terraforming Mars can be evaluated with a small set of system-level feasibility constraints linking (i) target pressures and compositions to required atmospheric inventories, (ii) target surface temperatures to required radiative control authority, (iii) inventories and radiative agents to sustained industrial throughput and power over a build time, and (iv) persistence against collapse, escape, and geochemical sinks. We use transparent order-of-magnitude scalings to map proposed levers (endogenous CO$_2$ release, synthetic super-greenhouse gases, CO$_2$-H$_2$ CIA, engineered aerosols/nanoparticles, orbital mirrors/albedo modification, and regional solid-state greenhouse ``paraterraforming'') onto common metrics $\{M,\ \tau_{\rm IR}/\Delta F_{\rm TOA} \dot M,\ P\}$. We find: (1) human-relevant pressures imply exaton-class inventories, $M_{\rm atm}\simeq 4\pi R_{\rm Mars}^2 P_s/g_{\rm Mars}\sim 10^{17}$-$10^{18}$ kg; (2) accessible CO$_2$ plausibly provides $\lesssim 20$ mbar, yielding $\lesssim 10$ K warming under present insolation; (3) achieving $T_s$ ~ 250-273 K at current insolation requires an effective IR opacity target $\tau_{\rm IR,eff}\sim 2$--4 (uncertain at the ~30-50% level but not altering mass-scale conclusions); (4) breathable endpoints are dominated by O$_2$ and buffer-gas mass and by a minimum oxygenation work $\gtrsim 10^{25}$ J, implying $\dot M\sim 10^{7}$-$10^{8}$ kg\,s$^{-1}$ and multi-$10^2$ TW to PW-class average power for century-to-millennial build times. We conclude that regional habitability gains via paraterraforming are plausible on near-term industrial scales, whereas global transformation of Mars requires multi-century planetary industry and becomes credible only under conditions of (a) massive exogenous volatile supply or much larger discovered inventories, and (b) sustained high-authority climate control and retention against sinks and loss.

We report the design and functionality of the Murchison Widefield Array Particle Detector Array (MWA PDA), an array of eight particle scintillation detectors deployed to Inyarrimanha Ilgari Bundara, the Murchison Radio-astronomy Observatory (MRO). The purpose of the instrument is to identify cosmic ray extensive air showers (EAS) occurring over the core of the MWA radio telescope, and generate a trigger to allow radio data on the event to be captured and analysed. The system also acts as a pathfinder for a much larger instrument to be deployed in the core of the low-frequency component of the Square Kilometre Array, SKA-Low, by the SKA's ultra-high-energy particles science working group. Here, we describe the instrument and associated infrastructure, which has been verified to comply with the strict radio-frequency emissions requirements of the MRO, and was deployed in November 2024. We present calibration data, which demonstrates the ability of each detector to identify individual atmospheric muons at the expected rate, and we characterise the temperature dependence of the system. We describe a sample of 35,500 EAS identified using multi-detector coincidence over a 13-day period, and show how the detector data can be used to reconstruct the arrival directions and approximate energies of these events. We conclude that the PDA can reliably trigger on and reconstruct EAS contained within the $\sim 103 \times 90$ m$^2$ core region, arriving within 20$^{\circ}$ of zenith, at primary cosmic ray energies above $\sim 4$ PeV. We have also verified that the detector array can generate triggers, allowing the capture of radio data from the MWA correlator for offline analysis.

We investigate the role of the accretion history in shaping the depletion radius of dark matter halos using a large cosmological N-body simulation. We show that the inner depletion radius, rescaled by the virial radius, depends strongly on the recent mass accretion rate (MAR) measured over a dynamical timescale, while exhibiting only weak dependence on halo mass. While this dependence mirrors that of the splashback radius and the two radii are tightly correlated, the depletion radius exhibits a more nuanced response to the detailed accretion mode. Specifically, we find that the dependence on MAR steepens at lower redshifts, aligning with self-similar spherical collapse models yet contrasting with the behavior of the splashback radius. This redshift dependence is largely driven by dynamic events, as it diminishes significantly when halos undergoing recent major mergers are excluded. Furthermore, we identify a dichotomy in the drivers of the depletion radius. For slowly accreting halos, the MAR is the primary dependence, whereas for rapidly accreting halos, other properties (shape, spin, concentration, and formation time of the central subhalo) related to the anisotropic or perturbed accretion mode also play a significant role. These results establish the depletion radius as a sensitive physical probe of the detailed accretion history of dark matter halos, complementary to the splashback radius.

Using the TNG50 simulation, we determine observationally motivated metrics that can distinguish quenching galaxies from star forming galaxies for $M_{*} \geqslant 10^{9.5}~M_{\odot}$, based on the spatial distribution of their stellar populations. Quenching galaxies are not fully quenched but have low levels of ongoing star formation that decreases over time. The morphological metrics consider the concentration of star formation, size of the star forming disk, and characteristic radii that trace sharp truncations of star formation. These metrics can separate simulated quenching galaxies based on morphology into populations where star formation is suppressed inside-out and outside-in. Inside-out quenched galaxies are more likely to be the most massive galaxy within their halo in the field, while outside-in quenched galaxies are satellites residing in dense environments and begin quenching ${\sim} 1~\text{Gyr}$ after being accreted. Outside-in quenched galaxies typically take ${\sim} 1.5~\text{Gyr}$ to quench, and inside-out quenched galaxies can take up to ${\sim} 3.5~\text{Gyr}$, where the duration of quenching is a function of stellar mass. We find that each population of quenched galaxy experiences evolution of their morphological metrics, where the different quenched populations reside in unique locations in parameter space. Galaxies in the later stages of quenching are more easily distinguished than those in the early stages, when compared to star forming galaxies. In addition, inside-out quenched galaxies can be distinguished compared to outside-in quenched galaxies, and the progress through the quenching episode can be estimated for both populations. These results have broad implications for distinguishing quenching galaxies in large galaxy surveys.

2024 YR4 is a 40-100 meter-diameter asteroid and former Torino Scale 3 object which currently has a roughly 4% chance of impacting the Moon on 2032 December 22, an event which recent studies suggest could pose a hazard on Earth due to impact ejecta. We present a search for, and identification of, potential precovery observations of the virtual lunar impactor in Intermediate Palomar Transient Facility (IPTF) survey data, as well as other publicly accessible surveys, dating from 2016. These candidate detections, not accounting for any currently-undetected Yarkovsky forces, predict a perilune of 22001 +/- 49 km and a perigee of 277534 +/- 46 km (relative to the center of each respective body) representing an improvement of > 300 times in the approach distance uncertainty above the existing orbit solution and, if confirmed, decisively ruling out a lunar impact in 2032. Using a matched filter tuned to 2024 YR4's predicted appearance in each image, we find the detection to be significant at Pnull = 5x10-9. The resultant possible orbit solution should be easy to confirm during 2024 YR4's 2028 approach to Earth, potentially greatly reducing the effort required by the planetary defense community at large to characterize 2024 YR4 before its potential lunar impact.

Upcoming next-generation sky surveys will detect large number of faint objects with magnitudes larger than 25. When objects are crowded within a limited a field of view, blending becomes unavoidable. Blending leads to the omission of many sources during photometry in these fields, which cause an underestimates of tens of percent in crowded fields, and remains a major challenge for existing source-extraction techniques. Although artificial neural networks had shown promising results in the detection and classification in wide-field surveys, they often fail with severely blended stars. We developed a robust deep learning model, Astro-RetinaNet, based on the Retinanet algorithm to detect and classify blended sources in single-band astronomical images. After training and evaluating the performance of our network on simulated images, we find precision of 0.96, 0.89,0.70, 0.50,0.75 for single star, 2-star, 3-star, 4-star and 5-or-more star blending cases, respectively, with star number density $\sim$22000 stars per $\rm arcmin^2$. We compare our method's detection capability and completeness both on CSST simulated NGC 2298 images and HST observed M31 images. In crowded and non-crowded stellar fields of simulated NGC 2298, our results show that the model can recover $82\%$ and $95\%$ sources respectively at magnitude ($i$ band) of 25, while for SExtractor and Photutils the completeness reduces to $20\%, 59\%$ and $60\%, 88\%$ respectively. In the M31 case, as faint as 27 magnitude ($F814W$) in a crowded field, Astro-RetinaNet detects 2,224 sources, significantly outperforming Photutils and SExtractor by factors of 3.4 and 7.1, respectively.

In view of the current and increasing evidence of a running Hubble constant, we investigate its redshift dependence within the flat $\Lambda$CDM framework using a 20-bin analysis of the Master SNe~Ia Sample \citep{2025JHEAp..4800405D}, considering cases with and without very low-redshift data. For each case, we obtain best-fitting values of $H_0$ and $\Omega_{m0}$, and employ both logarithmic \citep{2025arXiv250902636L} and power-law \citep{2021ApJ...912..150D,2022Galax..10...24D,2025JHEAp..4800405D} parameterizations. The two parameterizations are consistent over the redshift range considered and coincide for low redshifts. To assess their behavior at earlier epochs, we extrapolate both forms to the Cosmic Microwave Background radiation (CMB) era ($z\simeq1100$), Big Bang Nucleosynthesis (BBN, $z\sim10^{9}$), and inflationary scales ($z\sim10^{20}$). The reconstructed Hubble constant remains nearly indistinguishable up to the CMB scale, diverges at the few-to-ten percent level around BBN, and differs more substantially when extrapolated to inflationary redshifts. A qualitative distinction emerges at very-high redshift: the logarithmic form predicts a vanishing of $\mathcal{H}_0^{\mathrm{Log}}(z)$ at finite $z$, while the power-law form, $\mathcal{H}_0^{\mathrm{PL}}(z)$, approaches zero asymptotically as $z \rightarrow \infty$. In future studies, independent high-redshift observations and extensions beyond $\Lambda$CDM, such as $f(R)$ modified gravity, could allow a comparative study of the two parameterizations beyond the SNe~Ia regime and their high-$z$ physical implications.

We present global, three-dimensional general relativistic magnetohydrodynamic simulations of accreting black holes that incorporate pair physics. Pairs are modeled as a passive scalar that maintains a constant temperature. For high accretion rate models, we observe a maximum pair fraction of $\sim \mathcal{O}(0.01)$, consistent with those inferred from some X-ray binaries, and identify a `pair void' extending to a few gravitational radii from the black hole. Pair fractions peak in the midplane just outside the plunging region and within a thin strip at the base of the corona. For moderate to high accretion rate models, pairs are near equilibrium close to the disk midplane, where the scattering optical depth is high and pair equilibrium timescales are short, and could be comparable to the Coulomb collision timescale. This suggests the possibility of a pair-regulated coronal temperature. In contrast, the upper corona and jets, where the scattering optical depth is relatively low and pair equilibrium timescales are long, are populated with pairs that may exceed their equilibrium value by orders of magnitude. These pairs are transported by advection from the disk, which dominates over local pair processes. This result highlights advection as a significant source of pair injection, which may be relevant for certain X-ray binaries exhibiting $\gamma$-ray signatures. The pair density along the magnetically dominated poles exceeds the Goldreich-Julian density in some models.

The redshifted 21-cm signal is a unique probe of the early universe, particularly the Epoch of Reionization (EoR). While the 21-cm power spectrum has been the primary statistic for parameter inference, it fails to capture the non-Gaussian information in the signal, motivating the use of higher-order statistics such as the bispectrum. We perform a rigorous cross-simulation validation to infer the mean neutral hydrogen fraction ($\bar{x}_{HI}$) by training a Bayesian neural network on 21cmFAST simulations and applying it to mock observations generated by ReionYuga code. Our analysis spans six redshifts and includes realistic SKA system noise and cosmic variance, calculated from 50 statistically independent realizations. We find that the bispectrum adds useful information, but the improvement is moderate, with constraints tightened by $\sim 1.4\times$ the power-spectrum only case. The cross-simulation analysis also identifies a persistent systematic discrepancy between inferred and true values that often exceeds the statistical uncertainties, implying that modeling uncertainty remains the dominant limitation. Our results, therefore, indicate that the highly stringent constraints obtained in same-code validation studies may be overly optimistic, and mitigating cross-model systematics is crucial for robust parameter inference in the SKA era.

Compact binaries with orbital periods shorter than about 7 days show an absence of transiting planets, a feature known as the ``circumbinary planet desert". The physical mechanism behind this desert remains unclear. We investigate its origin by simulating the long-term dynamics of multi-planet circumbinary systems with evolving inner binaries. Our simulations are based on the single-averaged secular equations that average only over the binary orbital period and fully incorporate planet-planet interactions. When an eccentric binary decays via tides, an outer planet can be captured into resonance advection in eccentricity, a state in which its apsidal precession locks with that of the binary, driving extreme eccentricity growth. While such growth can occur in a binary-single planet system, the parameter space is limited and may not necessarily induce instability. In a multi-planet system, however, the excited orbit inevitably crosses those of its neighbors, which triggers violent planet-planet scatterings and produces collisions or ejections. Crucially, these mutual gravitational interactions amplify the ``localized" instability of a single planet into a system-wide chain reaction, drastically reshaping the orbital architecture and potentially clearing out the inner regions of planetary systems. Our results suggest that the resonance-induced instability provides a natural explanation for the observed circumbinary planet desert.

Galaxy clusters imprint a distinctive signature on the cosmic microwave background through the thermal Sunyaev-Zel'dovich (SZ) effect, which enables to study the intracluster plasma distribution and makes them powerful cosmological probes. We present the first Sardinia Radio Telescope (SRT) detection of the SZ effect in the galaxy cluster MACS J1752+4440 at 18.6 GHz, with a resolution of 0.9'. We detected a decrement in brightness toward the cluster centre, which we attributed to the thermal SZ effect. We modelled the signal using a spherically symmetric $\beta$ model for the electron density distribution and we employed a Bayesian retrieval to estimate the core radius, central electron density, and $\beta$ parameter of the cluster. We found values consistent with expectations for a galaxy cluster of the mass of MACS J1752+4440: a core radius of ($160 \pm 30$) kpc, a central electron density of ($2.5^{+0.7}_{-0.5} \cdot 10^{-3}$) cm$^{-3}$, and $\beta$=$0.6\pm0.1$. The mean Compton-$y$ parameter within a radius of 3.5' is $(2.6 \pm 0.3) \cdot 10^{-5}$, higher than the value reported by Planck, which is coherent considering the different resolution of the instruments and the modelling adopted. This work demonstrates the potential of the SRT to detect the onset of the SZ decrement at low frequencies, providing higher angular resolution than current all-sky surveys and enabling an improved reconstruction of the SZ decrement profile and the plasma distribution in the intracluster medium.

Context. Various simplified models have been investigated to understand the complex dynamical environment near irregular asteroids. We propose a generalized dipole-segment model (GDSM) to describe the gravitational fields of elongated bodies. The proposed model extends the dipole-segment model (DSM) by including variable pole masses and a connecting rod while also accounting for the spheroidal shape of the poles instead of assuming point masses. Methods. A nonlinear optimization method was employed to determine the model parameters, which minimizes the errors between the equilibrium points predicted by the GDSM and those obtained using a more realistic approach, such as the polyhedron model, which is assumed to provide the accurate values of the system. The model was applied to three real irregular bodies: the Kuiper belt objects Arrokoth, Kleopatra, and comet 103P/Hartley. Results. The results show that the GDSM represents the gravitational field more accurately than the DSM and significantly reduces computational time and effort when compared with the polyhedron model. This reduction in computational complexity does not come at the cost of efficiency. This makes the GDSM a valuable tool for practical applications. The model was further employed to compute heteroclinic orbits that connect the unstable triangular equilibrium points of the system. These trajectories, obtained from the intersections of the stable and unstable manifolds, represent natural pathways that enable transfers between equilibrium regions without continuous propulsion. The results for Arrokoth, Kleopatra, and 103P/Hartley are consistent and validate the GDSM as an accurate and computationally efficient framework for studying the dynamical environment and transfer mechanisms around irregular small bodies.

With the rapid advancements in observational technologies and the widespread implementation of large-scale sky surveys, diverse electromagnetic wave data (e.g., optical and infrared) and non-electromagnetic wave data (e.g., gravitational waves) have become increasingly accessible. Astronomy has thus entered an unprecedented era of data abundance and complexity. Astronomers have long relied on unimodal data analysis to perceive the universe, but these efforts often provide only limited insights when confronted with the current massive and heterogeneous astronomical data. In this context, multimodal data fusion (MDF), as an emerging method, provides new opportunities to enhance the value of astronomical data and deepening the understanding of the universe by integrating information from different modalities. Recent progress in artificial intelligence (AI), particularly in deep learning (DL), has greatly accelerated the development of multimodal research in astronomy. Therefore, a timely review of this field is essential. This paper begins by discussing the motivation and necessity of astronomical MDF, followed by an overview of astronomical data sources and major data modalities. It then introduces representative DL models commonly used in astronomical multimodal studies, the general fusion process as well as various fusion strategies, emphasizing their characteristics, applicability, advantages, and limitations. Subsequently, the paper surveys existing astronomical multimodal studies and datasets. Finally, the discussion section synthesizes key findings, identifies potential challenges, and suggests promising directions for future research. By offering a structured overview and critical analysis, this review aims to inspire and guide researchers engaged in DL-based MDF in astronomy.

Plumes have been proposed to channel MHD waves and the solar wind into the heliosphere. High-speed propagating disturbances (PDs), though well detected in plumes, cannot yet be clearly assigned to MHD waves or to mass flows. Additionally, plume bases as observed in the extreme ultraviolet are riddled with small-scale transients that could be related to the PDs. We study three plumes within an equatorial coronal hole observed by the EUV High Resolution Imager of the Extreme Ultraviolet Imager on board Solar Orbiter. The properties of the small-scale brightenings at the plume bases are investigated to interpret their nature and possible relation with PDs. We process images with the Difference of Gaussians method to highlight the target brightenings, which are further identified with two different approaches. In the 30-min observation, 50 brightenings are visually selected, which also help set thresholds for automatic detection, where we find 451 brightenings. Their properties, including velocities on the plane of sky (PoS), are analyzed statistically. Potential field extrapolation based on the magnetic field data from the Polarimetric and Helioseismic Imager on board Solar Orbiter is used for correcting the PoS velocity to the real velocity along the magnetic field. We observe that the majority of the base brightenings are small-scale, short-lived, and slightly elongated at the plume bases. They display intricate movements, with most exhibiting velocities in the PoS of less than 10 km/s. Their 3-dimensional velocities are found to be substantially lower than (and difficult to reconcile with) the speeds of PDs. A direct link between base brightenings and PDs remains inconclusive. We propose two possibilities for base brightenings: they may be related to wave-driven Type I spicules or originate from interchange reconnections. Further investigation is required to validate these hypotheses.

We present a detailed analysis of the trajectory of the interstellar comet 3I/ATLAS, focusing on its non-gravitational acceleration (NGA) parameters and their uncertainties. Orbital solutions are computed with models that implement symmetric, time-offset, and asymmetric radial dependence of the outgassing law relative to perihelion. We assess solution robustness through multiple data-selection strategies and comparison with independent determinations. The radial and normal NGA components (A_1 and A_3) are broadly consistent across all configurations, whereas the transverse component (A_2) is more sensitive to data selection, parameter correlations, and orbital phase coverage. Models with asymmetric radial decay slopes marginally improve the fit, but they also introduce additional degeneracies, contributing to systematic uncertainties. The magnitude of the NGA scaled to 1 au constrains the nucleus size of 3I. While our total acceleration estimates agree well with that of JPL's Small-Body Database solution, inclusion of systematic modeling effects implies a significantly larger uncertainty in the inferred radius of 3I.

The dynamical effects on the stellar motion produced by the Galactic central bar and the spiral arms perturbations are investigated separately and compared. The stars from the Gaia DR3 catalog are selected in the region of observable completeness, which we estimate as $\sim$1 kpc from the Sun. We apply the 2D model of the Galactic potential consisting of three axisymmetric components, the disk, the bulge, and the dark matter halo, and two non-axisymmetric components, the central bar and the spiral arms. The stellar dynamics is studied using analytical and numerical techniques, such as Hamiltonian topology analysis, the construction of dynamical maps on the representative planes, dynamic spectra, and Poincaré sections. We identify the main dynamical features in the solar neighborhood (SNd), the corotation (CR) and Lindblad resonances (LRs). By assuming that the main moving groups (MGs) in the SNd originate from the resonances, we compare their locations, structures, and intensities with the theoretical predictions and provide a description of the process involved in the formation of the MGs. In addition, we explore parametric planes by adjusting the values of the pattern rotation speed $\Omega_{p}$ with the positions of the MGs, for both the spiral arms and bar models, and conclude that the spiral arms model shows better results when compared to those of the bar, under the hypothesis of the dynamical origin of MGs.

Cosmic explosions play a critical role in a broad range of astrophysical fields. Although considerable progress has been made to understand the explosive engines and their progenitors, many of the details are not well understood. One of the most powerful electromagnetic probes of the explosive mechanism and the stellar progenitor is the first burst of photons emitted from this blastwave as it exits the stellar photosphere, known as shock breakout (SBO). Our understanding of SBO has evolved considerably in the past decade. Shock heating as the blastwave propagates through the star and circumstellar material can drastically alter this emission producing a much broader range of potential SBO signals than that predicted by standard analytical approaches. Here we present a semi-analytic approach to model this diverse SBO emission, focused on thermal Bremsstrahlung radiation, which more accurately captures the complexities in Nature over previous treatments. We calculate a range of signals for a range of supernova and gamma-ray burst types. Our models demonstrate how we can use these signals to place constraints on the nature of the explosive engines and better understand the role SBO can play in prompt gamma-ray bursts. We study the implications of these results to historic observations, Einstein Probe transients, and in the context of proposed missions. We find that stripped envelope events can be detected serendipitously with survey telescopes, but type Ia and II SBO detections require fast-pointing X-ray observations in response to early warning alerts from gravitational wave or neutrino detectors.

Non-local models of stellar convection can account for mixing effects in regions adjacent to convectively unstable layers and for changes to the mean temperature structure caused by free, buoyancy driven convection. The physical completeness of such models, however, depends on how third order correlations, which characterize the non-local transport processes, are expressed in terms of second order correlations and the stellar mean structure. Physical arguments and 3D hydrodynamical simulations were used to develop and test new closure relations for the skewness of the vertical velocity and temperature fields and third order cross-correlations to improve the predictive capabilities of non-local models of convection used in stellar astrophysics and in other disciplines such as meteorology. The structural form of the closure correlations was developed by a series of physical arguments and their accuracy was evaluated through self-consistency tests based on 3D hydrodynamical simulations for the Sun and a DA type white dwarf. The new closure relations derived for the skewness of vertical velocity and temperature fields provided improvements of up to an order of magnitude compared to previous models. This allows releasing the full potential of closure relations for the vertical velocity and temperature cross-correlations previously proposed in meteorology as well as the construction of new, more reliable models for the third order moments of vertical velocity and temperature in non-local models of turbulent convection. The new models for the skewness and third order cross-correlations of vertical velocity and temperature permit the construction of non-local models of turbulent convection which remove, among others, several major short-comings of three equation non-local convection models that are based on the downgradient approximation.

The accretion-ejection activities of black holes play a vital role in shaping the Universe. Bright and recurrent black hole X-ray binaries are ideal objects for studying accretion physics across a wide range of accretion rates, providing insights into the understanding of their supermassive counterparts. This short review summarizes X-ray techniques capable of measuring accretion geometry, our current understanding, and open questions. In particular, X-ray spectroscopic studies indicate that the accretion disk can extend close to the innermost stable circular orbit in the bright hard state. Some hints of disk-corona-jet connections are also discussed.

The axion, a well-motivated hypothetical particle arising in extensions of the Standard Model, can be produced copiously within the hot, compact cores of white dwarf stars. The shape of the white dwarf luminosity function (WDLF) is a powerful tool for constraining theoretical particles that would imply an additional cooling channel in white dwarfs. In this work, and for the first time, we use the 100-parsec Gaia DR3 white dwarf sample and compare it with theoretical predictions. We have simulated synthetic populations of white dwarfs using a population synthesis code based on Monte Carlo techniques, incorporating realistic observational errors, and based on state-of-the-art white dwarf models that incorporate the anomalous cooling caused by the presence of axions. Axion bremsstrahlung emission rates were implemented using the latest theoretical calculations. We find that, for the brightest white dwarfs in the sample ($M_{\mathrm{Bol}} < 10$), the $\chi^2$ statistic is largely insensitive to the assumed stellar formation rate (SFR), which is typically the dominant uncertainty in modeling the Galactic-disk WDLF. The resulting $\chi^2$ analysis disfavors a sizable additional cooling contribution. This conclusion contrasts with earlier studies in which axion-electron couplings in the range $0.7 \times 10^{-13} < g_{ae} < 2.1 \times 10^{-13}$ provided mildly improved fits to the Galactic-disk WDLF. We attribute the discrepancy to simplifying assumptions in previous modeling and to the substantially improved observational quality of the 100-pc Gaia DR3 sample. We obtain the upper limit $g_{ae} < 1.68 \times 10^{-13}$ ($95\%$ C.L.), which is among the strongest available.

We investigate the influence of nonlinear gravity corrections, arising from the Eddington-inspired Born-Infeld (EiBI) theory on wave dynamics, stability, and energy transport processes in polytropic, viscous, and turbulent solar plasmas. Analytical and numerical analyses of the Jeans-normalized quadratic dispersion relation demonstrate that both the EiBI gravity parameter $(\chi)$ and the relative polytropic sound speed $(\beta)$ independently regulate oscillation frequencies, growth rates, phase velocities, perturbation energy partitioning, and outward acoustic energy flux. Positive $\chi$ systematically elevates oscillation frequencies, phase velocities, and outward energy flux level by $\sim$10% relative to the Newtonian predictions, while larger $\beta$ enhances them by up to 55%, thereby promoting wave propagation and efficient acoustic transport. Conversely, negative $\chi$ strengthens gravitational binding and increases damping rates by $\sim$40%, particularly for the \textit{g}-modes. Energy partitioning analyses reveal that the EiBI corrections fundamentally restructure the kinetic-electrostatic-gravitational energy balance. While the Newtonian gravity contributes negligibly ($<$4%), nonzero $\chi$ channels up to one-third of oscillation energy into gravitational modes. The modal surface flux calculations further confirm that only the \textit{p}-modes drive outward energy transport (amplification for $\chi>0$, suppression for $\chi<0$). A direct comparative analysis with four years of SDO/HMI Doppler velocity observations demonstrate a robust theoretical agreement for $\chi=3\times10^7$ m$^5$kg$^{-1}$s$^{-2}$, providing the first empirical constraint on the solar EiBI gravity through helioseismology. The findings offer a rigorous framework for advancing our understanding about solar plasma stability, helioseismic signatures, and ambient atmospheric energy transport processes.

The discovery of many Earth-like planets has renewed interest in whether life and technological civilizations exist elsewhere. The Search for Extraterrestrial Intelligence (SETI) seeks evidence for technological civilizations via technosignatures across the electromagnetic spectrum. Here, focusing on artificial radio emissions with extremely narrowband signals, we model Earth as a distant, unresolved source and simulate its narrowband transmissions as observed with current and forthcoming radio facilities. Planetary rotation induces small but coherent Doppler drifts (maximum fractional shift of order $10^{-6}$) that imprint a characteristic, time-varying pattern on the spectrum. We develop a forward-inverse framework that exploits this modulation: adopting a population-weighted model for terrestrial transmitters, we compute time-resolved spectra and then apply a new inversion method that reconstructs the underlying transmitter distribution from the temporal pattern of fractional frequency offsets. In noise-added tests, the method recovers the low-order spherical-harmonic structure of the map and retrieves major population centers despite the north-south degeneracy of unresolved observations. The recovered distribution is expected to correlate with continents, climate zones, and population density. This approach moves SETI beyond mere detection, enabling quantitative cartography of a civilization's activity and inference of host-planet properties through sustained, time-resolved spectroscopy.

We present chemical abundance measurements for 7,302 red giant branch stars within the half-light radius (~5') of $\omega$ Centauri ($\omega$ Cen), derived from MUSE spectra using the neural network model DD-Payne. DD-Payne effectively identifies spectral features of C, N, and O for [Fe/H]>-1.0 dex; Mg for [Fe/H]>-1.5 dex; and Na, Ca, and Ba for all metallicities. By combining these measurements with previous high-resolution studies, we create the most comprehensive picture of $\omega$ Cen's rich chemical evolutionary history. For the first time, we map elemental variations across the entire chromosome diagram, which is widely used to identify multiple populations. We analyze the median chemical abundance trends as functions of age and metallicity for different subpopulations. The DD-Payne measurements of [C/Fe], [N/Fe], and [O/Fe] extend literature trends to higher metallicities and show continuous abundance-metallicity relations, with [(C+N+O)/Fe] increasing steadily with [Fe/H]. [Ca/Fe] and the s-process element [Ba/Fe] also increase with metallicity across all populations. For [Ba/Fe], the chemically enhanced (P2) populations are more enriched than primordial (P1) and the intermediate (Im) populations. Furthermore, [N/Fe] correlates strongly with stellar age while [Ca/Fe] and [Ba/Fe] exhibits a weaker age dependence. Using these abundance-metallicity-age relations, we evaluate different formation scenarios of $\omega$ Cen proposed in the literature. Our study demonstrates that combining MUSE with machine learning enables large-sample stellar abundance measurements in crowded cluster cores, overcoming the limitations of fiber-fed spectroscopy for studying multiple stellar populations and their evolutionary histories.

We investigate the velocity field derived from HI measurements of the irregular galaxy ESO 358-60 using the Velocity Ring Model (VRM) method. This technique, which assumes a coplanar disk, allows us to reconstruct coarse-grained maps of both radial and tangential velocity components from the observed line-of-sight velocity field. Such maps reveal that tangential motions dominate the inner regions, while radial motions become increasingly significant toward the outskirts. This kinematic behavior contrasts with that inferred from the Tilted Ring Model (TRM), which suggests that radial motions are more prominent in the intermediate disk and negligible in the outskirts and detects a pronounced warp of approximately $20^\circ$, with the inner disk nearly edge-on and the outer regions inclined by approximately $60^\circ$. In contrast, the VRM analysis finds that the disk exhibits a bar-like structure in its central regions. This interpretation is further supported by the intensity and velocity dispersion maps. To test the origin of the TRM-derived warp, we construct a toy model based on the TRM results and analyze it with the VRM technique, finding evidence that the warp is likely an artifact arising from the TRM's assumptions. Finally, we estimate the galaxy's mass using both the standard dark matter halo model and a dark matter disk (DMD) model, where all mass lies in the disk plane. The DMD yields a total mass approximately three times lower and provides a slightly better fit to the rotation curve.

Wind accretion in binary systems is commonly described using the Bondi-Hoyle-Lyttleton (BHL) formalism. However, its standard implementation fails in the slow-wind regime, where the wind velocity of the donor star ($v_\mathrm{w}$) is comparable to or smaller than the orbital velocity of the accretor ($v_\mathrm{o}$). Tejeda & Toalá recently proposed a geometrical correction to the BHL formalism that accounts for the wind aberration caused by the binary's orbital motion, which tilts the accretion cylinder and reduces its effective cross-section. Here we present a suite of smoothed particle hydrodynamic simulations performed with PHANTOM to test wind accretion in binary systems operating in this slow-wind regime. We explore circular configurations and directly measure mass accretion efficiencies from the simulations. Our results confirm that the standard BHL prescription systematically overestimates accretion rates for $v_\mathrm{w}/v_\mathrm{o} < 1$, while the geometrically corrected model reproduces the simulated efficiencies with remarkable accuracy. A key finding is that the velocity relevant for accretion estimates is not the value derived from the unperturbed stellar wind, but the local gas velocity measured upstream of the accretor. The gravitational potential of the accretor perturbs the flow, altering the effective relative velocity and modifying the accretion efficiency, particularly for compact orbits. These results provide strong numerical support for the geometrically corrected framework and establish a physically motivated basis for modeling wind-fed accretion in interacting binaries, including symbiotic systems.

We report on a giant stellar flare from the RS CVn-type binary UX Arietis, detected with the Monitor of All-sky X-ray Image (MAXI) and followed by a 12-day optical spectroscopic campaign using the 3.8~m Seimei Telescope. The flare released $5 \times 10^{37}$~erg in X-rays (0.1--100~keV) and $(2$--$6) \times 10^{36}$~erg in the H$\alpha$ line, placing it among the most energetic events of its kind. The H$\alpha$ light curve showed sinusoidal modulation atop an exponential decay, consistent with reappearance of the flaring region due to binary rotation. At orbital phase 0, when the primary star is farthest from the observer, 40\% of the H$\alpha$ flux was obscured, while at phase 0.5 the full emission was visible. This suggests the H$\alpha$ emitting region is located at a relatively low latitude and is comparable in size to the stellar disk. Radial velocity modulation implies that the region lies at $\sim19\,R_{\odot}$ from the system's rotation axis, farther out than the stellar limb at $14.4\,R_{\odot}$. Photometric monitoring with the Chuo-university Astronomical Telescope revealed a large low-latitude starspot covering $\sim25\%$ of the surface. These findings are consistent with a scenario in which the flare occurred above the starspot, and the H$\alpha$-emitting plasma was magnetically confined in a loop extending at least $5\,R_{\odot}$ above the stellar surface. From the MAXI data and assuming a radiatively cooling plasma, the electron density and volume are estimated to be $10^{10}$~cm$^{-3}$ and $1 \times 10^{35}$~cm$^3$, respectively. If cubic in shape, this corresponds to $7\,R_{\odot}$, consistent with the H$\alpha$ region height. These results provide direct constraints on the geometry of the plasma and its spatial relationship with the starspot in one of the most energetic stellar flares ever observed.

The 21-cm signal from neutral hydrogen serves as a critical tool for unraveling the astrophysical processes that shaped cosmic dawn and the epoch of reionization. We explore the usually overlooked impact of re-scattering of 21-cm photons during and after reionization, similarly to cosmic microwave background photons. This scattering affects the observed brightness temperature by mixing the original signal with light scattered into the line of sight from other regions, effectively at the mean 21-cm brightness temperature. This gives a small but significant effect. We show that it attenuates the fluctuations in a 21-cm image by $4-7\%$, while reducing the 21-cm power spectrum by a scale-independent $7-13\%$ during cosmic dawn and reionization. Incorporating this correction is vital for precisely comparing theoretical predictions with observations from experiments such as NenuFAR, LOFAR, and HERA, and the upcoming Square Kilometre Array.

We study outflow from an unmagnetized, shocked accretion disk around a non-rotating super-massive black hole using multidimensional hydrodynamics simulation with radiative cooling. We aim to investigate whether such shocked accretion flow can launch sustained collimated bipolar outflow reaching out to thousands of gravitational radii even in the absence of magnetic field and if yes, what terminal velocity can they achieve. We present the results of a few simulations of geometrically thick accretion flow with increasing specific angular momentum on a vertically elongated cylindrical domain. We show that bipolar outflow from a region very close to the black hole is originating and propagating vertically out to our simulation domain boundary at around 2651 Schwarzschild radius. The outflow attains a terminal velocity with a maximum value found to be 0.14c and the outflow rate depends on the angular momentum value of the accreting material. We also compute the self-Comptonized bremsstrahlung spectra for all the disk-jet runs.

Context. This paper presents a detailed description of the Decametre Pulsar and Transient Survey of the Northern Sky that was carried out in 2012-2017 using the world's largest radio telescope at decametre wavelengths - UTR-2 in Ukraine. This extensive survey covers the northern sky from declination -10 to +80 deg , with a temporal resolution of 8 ms, and explores dispersion measures up to 30 pc/cc. Aims. The major advantage of the decametre wavelength range is a comparatively wide band, in which the dispersive delay due to the interstellar plasma reaches hundreds of seconds, giving us the opportunity to determine the dispersion measure with a very high accuracy. This enables us to discover new transients, while avoiding data contamination from numerous weak signals of a different nature. Methods. The drift-scan survey in 5-beam mode of UTR-2 was carried out at night time. To cover the entire sky along the right ascension, the duration of the sessions was more than 12 hours at a time close to the autumn and spring equinoxes (to obtain the same conditions for the interference situation). 90 degrees along the declination were covered by five beams in 40 days (each equinox). Results. We discovered 380 individual transient signals with dispersion measures significantly differ from those of known sources. We determined the parameters of each single transient signal. We show that they cannot be explained by ionospheric scintillations. Repeated observations have shown that some detected transient signals are repetitive and are thus likely to originate from pulsars or rotating radio transients. Key words. Stars: neutron - pulsars: general - Methods: data analysis - Methods: observational - Astronomical databases: Surveys

Young stellar moving groups offer unique opportunities to investigate the early evolution of stellar and planetary systems. In continuation of an ongoing effort to age-date compelling planetary systems, we provide an in-depth age analysis of KELT-20, a young A-type star that hosts a well-aligned ultra-hot Jupiter. This system poses a useful case study to investigate migration mechanisms at early stages of evolution. Using Gaia DR3 data, we identify 77 stars with proper motions consistent with KELT-20, including 19 with measured radial velocities that enable full 3D kinematic confirmation. Using isochronal analyses, gyrochronology, photometric variability, and stellar activity indicators, we converge on an age of $58 \pm 5\,$Myr. This constraint provides critical insights into the dynamical processes shaping hot Jupiter formation.

The disagreement between early and late Universe electromagnetic measurements of the Hubble constant, $H_0$, known as the Hubble tension, highlights the need for independent and complementary probes. Gravitational-wave events have recently emerged as such a probe for constraining cosmological parameters. $H_{0}$ inference using these events relies on sky localization and luminosity distance estimates, both of which can be significantly improved for strongly lensed events with appropriate lens modeling. In this context, we propose utilizing strong lensing of dark sirens, gravitational-wave events without identified electromagnetic counterparts, in combination with strong lensing of galaxies as a novel method for measuring $H_0$. The constant is inferred from the luminosity distances of these lensed dark sirens and the redshifts of their host galaxies, combining information from individual events to obtain statistically stronger constraints when multiple events are available. We adopt a simulated galaxy catalog, \texttt{MICECATv2}, as the basis for simulating strong lensing of galaxies and to provide the redshift information of host galaxy candidates required to infer $H_0$. We also examine the impact of galaxy catalog incompleteness on the resulting $H_0$ inference. Our results demonstrate that using only 8 strongly lensed dark sirens, analyzed with a dedicated galaxy-galaxy lensing catalog, can improve the precision of $H_{0}$ by roughly 50\% compared to 250 unlensed events.

Understanding the formation pathway for close-companion protostars is central to unraveling the processes that govern stellar multiplicity and very early star formation. We analyze a large sample of 51 Class 0/I close-companion protostellar systems, of which 38 show detectable outflows, yielding 42 measured outflows used in our analysis. We use ALMA observations of 11 systems in Perseus and 40 systems in Orion. These companions formed either directly at these small scales ($\lesssim 500$ au separations) via disk fragmentation or at larger scales ($> 1000$ au separations) via turbulent fragmentation followed by inward migration. Because of differences in formation mechanism, the former is expected to have preferentially aligned disks and outflows, whereas the latter is expected to show no preferred alignment. The relative prevalence of these formation pathways remains uncertain, yet it is critical to forming a comprehensive picture of star formation. We examine the distribution of position angles of companion protostars relative to the position angles of their molecular outflows. The outflow, as traced by $^{12}$CO ($J=2\rightarrow1$), is a useful proxy for the angular momentum of the system, expected to be orthogonal to the binary orbital plane. We use a simple model to account for random sampling of inclination and orbital phase in each system, finding that the observations are consistent with a distribution in which the outflows are preferentially orthogonal to the companions. Based on this analysis, we suggest disk fragmentation is the dominant formation pathway for close-companion protostellar systems.

We present ultra-high-resolution (0.03"~230 pc) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the hyperluminous dust-obscured galaxy W2305-0039 at z=3.111, targeting the CO J=7-6 and J=11-10 lines. The CO(11-10) emission is extremely compact and exhibits anomalously high excitation relative to CO(7-6) within the central <500 pc. X-ray-dominated region models successfully reproduce this excitation, providing strong evidence for intense X-ray irradiation by a deeply obscured active galactic nucleus (AGN), while photodissociation-region models fail to match the observed ratio. Forward modeling of the nuclear CO(11-10) position-velocity diagram yields a dynamical black-hole mass of log(M$_{\rm BH}$/M$_{\odot}$) = 8.3$^{+0.7}_{-0.6}$ and an intrinsic gas velocity dispersion of $277~^{+16}_{-14}$ km s$^{-1}$. Combined with the AGN luminosity from infrared spectral energy distribution decomposition, these measurements imply a highly super-Eddington accretion state with $\lambda~_{\rm Edd}~\gtrsim 4$. Our results provide dynamical evidence that the most rapid phases of black-hole growth can occur within a compact, heavily obscured nuclear region. Extending ALMA beyond its current 16 km maximum baselines will be essential for pushing such dynamical measurements to tens-of-parsec scales and resolving the black-hole sphere of influence in massive galaxies at $z \gtrsim 6$.

Time-domain photometric monitoring is essential for characterizing cometary evolution, particularly for rare interstellar objects with limited observing opportunities. We aimed to characterize the pre-perihelion photometric behavior and dust activity of the interstellar comet 3I/ATLAS, and to test the capability of the Black Hole Target and Observation Manager (BHTOM) platform and telescope network for coordinated high-cadence non-sidereal observations. We obtained 70 days of time-series photometry of 3I/ATLAS from 2025 July 4 - September 11 using 16 telescopes and 1554 images. The data were processed and calibrated with the BHTOM pipeline. High-cadence, multi-band imaging was used to measure the rotation period and color evolution, while the dust activity was quantified via Afp measurements. We present a pre-perihelion light curve of 3I/ATLAS from Rh = 3.18 - 2.19 au, which exhibited a steady increase of ~3 magnitudes with no evidence of anomalous behavior. We measured a rotation period of P_rot = 15.98 +/- 0.08 h. The relative dust production increased from A(0)fp ~600 - 1100 cm, and the upper limit on the dust mass-loss rate increased from \leq 217 kg/s to \leq 328 kg/s. We measured an activity index of n = -1.24 +/- 0.02, consistent with a well-developed dust coma. The colors were statistically non-changing, with only a weak, non-significant tendency for 3I/ATLAS to become bluer at 3.5 > Rh > 2.2 au.

Compton-thick active galactic nuclei (CT-AGNs), defined by column density $\mathrm{N_H} \geqslant 1.5 \times 10^{24} \ \mathrm{cm}^{-2}$, are so heavily absorbed that their X-ray emission is often feeble, even undetectable by X-ray instruments. Nevertheless, their radiation is expected to be a substantial contributor to the cosmic X-ray background (CXB), predicting that CT-AGNs should comprise at least $\sim$30\% of the total AGN population. In the Cosmological Evolution Survey (COSMOS), the identified CT-AGN fraction falls far below theoretical expectations, indicating that a substantial population of CT-AGNs is hidden due to their low photon counts or their flux below the current flux limits of X-ray instruments. This work focuses on identifying CT-AGNs hidden in mid-infrared (MIR)-selected AGNs. First, we selected a sample of 1,104 MIR-selected AGNs that were covered but individually undetected by X-ray. Next, we reduced the X-ray data in the COSMOS and analyzed multiwavelength data in our sample to derive the key physical parameters required for CT-AGN identification. Using MIR diagnostics, we first find out 7 to 23 CT-AGN candidates. Their subsequent X-ray stacking analysis reveals a clear detection at $>3\sigma$ significance in the soft band and only $>1\sigma$ significance in the hard band. We fit the stacked soft- and hard-band fluxes with a physical model and confirm that these sources are absorbed by Compton-thick material. However, CT-AGNs constitute only 2.1\% (23/1104) of our sample, significantly below the fraction predicted by CXB synthesis models, indicating that a considerable population of CT-AGNs remains missed by our selection. A comparison of host-galaxy properties between CT-AGNs and non-CT-AGNs reveals no significant differences.

In scientific and engineering domains, modeling high-dimensional complex systems governed by partial differential equations (PDEs) remains challenging in terms of physical consistency and numerical stability. However, existing approaches, such as physics-informed neural networks (PINNs), typically rely on known physical fields or coefficients and enforce physical constraints via external loss functions, which can lead to training instability and make it difficult to handle high-dimensional or unobservable scenarios. To this end, we propose PhysFormer, a generative modeling framework that is self-consistent at both the data and physical levels. PhysFormer leverages a low-dimensional, physically interpretable latent space to learn key physical quantities directly from data without requiring known high-dimensional physical field parameters, and embeds the physical process of radiative flux generation within the network to ensure the physical consistency of the generated spectra. In high-dimensional, degenerate inversion tasks, PhysFormer constrains generation within physical limits and enhances spectral fidelity and inversion stability under varying signal-to-noise ratios (SNRs). More broadly, this approach shifts the physical processes from external loss functions into the generative mechanism itself, providing a physically consistent generative modeling paradigm for complex systems involving unknown or unobservable physical quantities.

We investigate the infrared (IR) properties of a sample of local star-forming galaxies using the WISE data. Focusing on the 20 confirmed strong Lyman continuum (LyC) leakers ($f_\mathrm{esc}>5\%$) included in this sample, we find that the IR detection separates these strong LyC leakers into two populations. The IR-undetected strong leakers in our sample exhibit high [OIII]5007/[OII]3725,3727 (O32) ratios, blue UV slopes, and low stellar masses, consistent with the classical density-bound scenario where the entire ISM is highly ionized. However, IR-bright strong leakers display unexpectedly low O32 ratios while maintaining a substantial escape fraction of $f_{\text{esc}} \sim 12\%$. They also have higher stellar masses and redder UV slopes, similar to weak or non-leakers, along with low nebular extinction, suggesting that the LyC photons are likely to escape through low-column-density channels in a high-density, ionization-bound clumpy medium. Morphologically, these two populations echo the compact starburst and merger-driven LyC leakers observed at $z \sim 3$, indicating that both escape pathways coexist across cosmic time. Our results demonstrate that efficient LyC escape is not limited to the lowest-mass, dust-poor dwarfs but can also occur in more massive, dusty environments, highlighting the potential contribution from dusty systems and posing new questions about the actual dominance of different galaxy populations during the Epoch of Reionization.

Little red dots (LRDs) draw extensive attention because of their unique observational characteristics and apparent overabundance in the early Universe, raising new insights into early black hole formation and growth. Early studies show that LRDs exhibit weak variability in broad-band photometry and emission-line fluxes, suggesting a preference for super-Eddington accretion or disfavouring an AGN origin. However, the cadence of the current data, and therefore, the resulting light curves for LRDs, is limited, preventing us from placing strong constraints on their variability. Based on Zwicky Transient Facility (ZTF) light curves with a baseline of $\sim6$ years, we here study the optical variability of seven previously reported local analogs of LRDs at $z \sim 0.3$, offering an insight into LRDs from a low-redshift sample. Three out of seven local analogs show excess variances on all three bands of their light curves, and two of them can be fitted with the damping random walk model, supporting their AGN origins for the variability. The remaining sources show weak variance in at least one band, but no detectable variability at the current sensitivity level, exhibiting $\rm SF_\infty$ upper limits consistent with estimates from high-redshift (high-$z$) LRDs. Their non-detection of variability is likely due to the large photometric uncertainty. As an implication, by simulating long baseline light curves with the variability amplitude of local analogs and adopting JWST observation cadence, we investigate the limitation of the variability amplitude estimate for LRDs. Our mock observations imply that the current constraints on LRDs' variability are probably underestimated. This underestimation might be induced by the short temporal baseline of observations, as well as the intrinsic scatter of the empirical $M_{\rm BH}-\tau$ relation.

High-redshift Lyman continuum emitters (LCEs) are crucial for understanding how galaxies ionize the neutral hydrogen in the epoch of reionization. However, detected LCEs at $z > 4$ are quite rare. Here we report an LCE candidate at $z = 4.444$, dubbed LCEz4-M1, which is the highest-redshift LCE to date with LyC detections confirmed in two independent data sets. The redshift is determined from the $Ly\alpha$ emission line detected in the VLT/MUSE spectrum. The LyC signal is detected independently in the Hubble Space Telescope (HST) F435W image and the VLT/MUSE spectrum at significances of $\simeq 3.7~\sigma$ and $\simeq 2.8 - 3.0~\sigma$, respectively. The centroid of the LyC emission is closely aligned with the rest-frame optical continuum traced by James Webb Space Telescope (JWST) imaging, with an offset of $\simeq 0''06$ (0.4 kpc in physical scale). Based on HST/ACS F435W photometry and MUSE spectroscopy, we infer LyC escape fractions of $f_{esc}(F435W) = 0.38_{-0.15}^{+0.25}$ and $f_{esc}(MUSE) = 0.33_{-0.13}^{+0.22}$. Using the combined JWST and MUSE data set, we characterize the physical properties and morphology of LCEz4-M1. The galaxy is compact and lies in the starburst regime, with a high star formation rate surface density of $\Sigma_{\rm SFR} \simeq 7~M_{\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, consistent with conditions that can drive strong feedback and outflows. The feedback may generate low-column-density pathways in the interstellar medium that facilitate LyC escape. While we find no clear evidence for an ongoing major merger, the presence of a faint companion ($\sim 0''05$) detected in the F277W band suggests a potential minor interaction. This is also consistent with LCEz4-M1 residing in an overdense environment, where elevated interaction rates and dynamical perturbations are expected.

We study the role of domain walls and their relics in the very early Universe within the framework of the mode-matching technique. Domain walls formed during the spontaneous breaking of discrete symmetries are modelled as a short lived contribution to the background energy density, leading to a controlled deviation from standard slow-roll inflation. Solving the modified Friedmann equation, we obtain a smooth, time-dependent Hubble parameter that asymptotically approaches the standard inflationary value. We analyze the evolution of scalar perturbations using the gauge-invariant Mukhanov variable and perform mode matching across the transition between the domain wall affected phase and standard inflation. We find that only modes exiting the horizon near the transition experience a modified evolution, while modes far from the transition remain unaffected. As a result, the primordial scalar power spectrum exhibits localized, scale-dependent deviations from scale invariance, while preserving the overall success of inflation. These results demonstrate that unstable domain-wall relics can leave subtle imprints on primordial fluctuations, potentially providing a probe of early-Universe phase transitions.

Compared to the commonly used planet occurrence rates, the multiplicity distribution of planets can be more useful in constraining the formation and evolution pathways of planetary systems. This work follows an earlier work of Zhu (2022) and derive the intrinsic multiplicity distribution of giant planets (with masses above Saturn mass) from two independent radial velocity (RV) surveys. In particular, we find that $(7.8\pm1.4\%, 2.3\pm1.2\%, 0.5^{+0.8}_{-0.3}\%)$ of Sun-like stars in the HARPS sample have $(1, 2, 3)$ giant planets within 10 au, whereas $(7.3\pm2.8\%, 7.2\pm2.3\%, <1.3\%, 1.0^{+1.0}_{-0.6}\%)$ of Sun-like stars in the California Legacy Survey (CLS) have $(1, 2, 3, 4)$ giant planets within 10 au. Here we have further cleaned the CLS sample and removed planet detections that were not discovered in the survey mode. The total fraction of Sun-like stars with giant planets within 10 au from the two samples are $10.6\pm1.2\%$ and $15.8\pm2.1\%$, respectively, and the difference may be accounted for by their different metallicity distributions. We briefly discuss the theoretical implications of our results. In particular, the inferred giant planet multiplicity distribution is inconsistent with most of the proposed theoretical models involving planet-planet scatterings, which predict either too many giant planets or too many systems with multiple giants.

Recent papers by Ansari et al. (2024, Science Advances 10, eadn4650) and Richardson et al. (2025, arXiv eprint 2504.01455) have suggested that global warming of the Martian surface ('terraforming') by 35 K to sustain local habitats above the melting point of water could be achieved through the injection of engineered aerosols into the Martian atmosphere. Using the MarsWRF 3D Global Climate Model, we investigate how artificial warming of Mars through engineered aerosol release would affect the planetary water cycle and the distribution of the major surface ice reservoirs. Within our model framework, every 20 K of global warming induces a tenfold increase in atmospheric water vapour content due to sublimation of H2O ice from the North Polar Cap. This increases the potency of cloud radiative feedbacks which induces nighttime warming (~5-10 K) at low latitudes, but daytime cooling (up to 40 K) in the winter midlatitudes. Water is transferred from the edge of the North Polar Cap to the South Polar Cap and there is minor destabilisation of shallow northern midlatitude subsurface ice. As a result, seasonal sublimation of H2O ice from the South Pole has an increased impact on the global water cycle. These changes persist on Mars at least decades after loading of the atmosphere with engineered aerosols ceases. Our model is limited by the gaps in our knowledge of present-day Martian weather and climate, and of the microphysics and radiative properties of candidate warming agents. Much more data is therefore needed before warming Mars could become feasible.

Braking of pulsars, or the law of their spin deceleration, is a manifestation of the combination of various processes occurring in the magnetospheres of neutron stars and their internal structural dynamics. The interaction of superfluid neutron vortex lines with superconducting magnetic flux tubes in the neutron star core plays a significant role in rotational evolution and magnetic field evolution through spin glitches and magnetic flux expulsion. In this study, the effect of this interaction on the temporal variation of the braking index of pulsars is investigated. The variation of the deviation from the $n=3$ value predicted by the generally accepted magnetic dipole formula with respect to physical parameters such as the magnetic field and the age of the neutron star has been elaborated, and applications have been made to pulsars with reliably measured braking indices.

SPACE-Timers are a lightweight hierarchical profiling framework for C++ designed for modern high-performance computing (HPC) applications. It uses a stack-based timing model to capture deeply nested execution patterns with minimal overhead, representing runtime behaviour as a tree of timing nodes with precise attribution. The framework provides structured reports with recursive aggregation, detection of unaccounted time, and compact visual summaries of runtime distribution, supporting both quick inspection and detailed analysis. It also includes checkpointing and error detection mechanisms. SPACE-Timers supports multiple profiling backends, including NVTX, ITT, ROCtx, and Omnitrace, and integrates with the MERIC runtime system to enable energy-aware optimisation. Its successful use in OpenGadget3 demonstrates its effectiveness for large-scale scientific applications.

We investigate the spectro-polarimetric properties of the long-duration GRB~220107A, which exhibited two distinct emission episodes separated by a 40 s quiescent gap, to test whether such multi-episode bursts show evidence for evolution in their underlying radiation mechanisms. We analyzed prompt emission data from AstroSat/CZTI, Fermi/GBM, and Konus-Wind, performing spectro-polarimetric analysis for each emission episode. The time-integrated polarization analysis shows no significant detection (PF$ < 38 \%$, $2\sigma$). Time-resolved analysis reveals clear spectral evolution between the two episodes, with episode 1 exhibiting a hard low-energy photon index and episode 2 showing substantial spectral softening ($\alpha \sim -0.72$). Regarding polarization: Episode 1 shows a low polarization upper limit (< 52\%), consistent with expectations for photospheric emission dominated by quasi-thermal Comptonization in a baryon-rich outflow. Episode 2 also shows overall low polarization (PF$ < 55 \%$, $2\sigma$), though sliding-window analysis yields a marginally elevated signal (PF$= 70 \pm 30\%$, BF = 2.8) between T0+76 to T0+88 s. The robust spectral softening between episodes could arise from sub-photospheric dissipation, optically thin synchrotron radiation in small-scale magnetic fields, or if the tentative polarization enhancement proves intrinsic, it would favor synchrotron emission in large-scale ordered magnetic fields. The spectral evolution of GRB 220107A, combined with our polarimetric constraints, demonstrates the diagnostic potential of time-resolved spectro-polarimetry for constraining GRB prompt emission physics. We present GRB 220107A as a test case illustrating how future higher sensitivity observations could discriminate between competing emission models for multi-episode bursts. Our results emphasize both the promise and current limitations of prompt phase polarimetry.

Recent JWST surveys have revealed two puzzling high-redshift phenomena: (1) an unexpectedly large abundance of flattened, disk-like galaxies at z > 3, and (2) a rare population of compact, extremely red sources at z ~ 4-9 ("Little Red Dots," LRDs) that often show V-shaped SEDs and very broad Balmer lines. These findings lack a consensus interpretation and have motivated models ranging from dusty starbursts to obscured AGN and more exotic scenarios. We propose that both phenomena are linked by a simple geometric consequence of a third clue: mounting evidence from structural modeling and axis-ratio statistics indicates that many low-mass galaxies at z > 3 are intrinsically prolate (cigar-like), not oblate rotation-supported disks. In this picture, a substantial fraction of the flattened, disk-like morphologies reported at z > 3 arise from side-on and intermediate-angle projections of prolate systems, while the rare near end-on views appear extremely compact and high-surface-brightness, and are preferentially reddened by the maximal line-of-sight column, naturally matching key elements of LRD selection. The expected fraction of near end-on systems, $P_{\text{end}} = 1 - \cos(\theta_{\text{max}})$, is ~1-3% for $\theta_{\text{max}}$ ~ 10°-15°, consistent with LRD demographics in wide JWST fields. This orientation-based framework does not exclude AGN or starburst activity; rather, it explains LRD rarity as an orientation effect and provides a natural route to the large columns of gas/dust and scattering depths inferred in recent dense-gas and electron-scattering interpretations of LRD spectra, without fine-tuned new physics. The model makes falsifiable predictions to validate or rule out this geometric interpretation.

The evolution of planet-forming disks and the processes of planet formation influence each other, and both of them are possibly impacted by the local environment. Extensive high-resolution imagery of disks across space and time is the best tool for determining their evolution. We compiled a comprehensive list of disk-bearing young stars with near-IR high-contrast images available. The sample sums up to 268 sources, including 51 targets with no prior publications, which makes this study the largest of its kind and the most extensive release of IR disk images to date. Our census reveals very diverse disk and ambient morphologies. Disks in Lupus are bright, in Chamaeleon are faint, in Corona Australis and Taurus are frequently surrounded by ambient emission. Disks experience an abrupt increase in IR brightness between 2 Myr and 5 Myr. The earliest IR disk cavities around single stars arise after 2-3 Myr explaining why are young disks faint in the near-IR, and determining which disks can live longer. Well-known, high-longevity disks (>8 Myr) are always bright. Ambient material is detected in more than 20% of young sources but the fraction drops with time. We find a clear correspondence for the presence of ambient material with the stellar variability, near-IR excess, and mass accretion rate as well as, in turn, with spirals and shadows in disks. Half of the disks with ambient material show spirals while none of them show rings. We therefore propose that the spirals and the disk warps responsible for shadows are generally induced by late infall from the medium, and that this also affects the stellar accretion. The emerging picture proves the fundamental role of the environment for the disk evolution and planet formation.

Stellar mass, size, and star formation rate (SFR) are fundamental properties that encode the structural and evolutionary states of galaxies. Observations reveal a mass-SFR-size relation whereby galaxies become more compact both above and below the ridge of the star-forming main sequence (SFMS), linking galaxy structure to star formation activity. We investigate this relation by comparing galaxies from two cosmological hydrodynamical simulations, EAGLE and TNG100, with observational samples from SDSS and CANDELS over three redshift intervals (0 < z < 0.2, 0.5 < z < 1.5, and 1.5 < z < 2.5). Both simulations reproduce the observed trend that galaxy sizes decrease with increasing offset away from the SFMS. This trend, however, weakens and is not detected in the observational sample at 1.5 < z < 2.5, likely due to increased measurement uncertainties. In contrast, the trend persists in both simulations up to z = 2.5. Across all redshifts, EAGLE predicts a stronger size dependence on SFMS offset than observed, whereas TNG100 exhibits a weaker dependence. We discuss how this mass-SFR-size relation can be understood in terms of different time variability in star formation rate across the SFMS.

Analysis of binary-lens microlensing events typically requires intensive computation because of the multimodal and complex posterior distributions. With the recent development of the JAX-based differentiable binary-lensing modeling package microlux, we present an analysis of two microlensing events with planet/brown-dwarf candidates, KMT-2025-BLG-1314 and KMT-2025-BLG-1392. Both events exhibit the "Close/Wide" degeneracy, and KMT-2025-BLG-1314 suffers from the "Planet/Binary" degeneracy and a recently recognized "Point/Finite" degeneracy among the planetary solutions. For KMT-2025-BLG-1314, the binary mass ratio is $\log q \sim -3.5$ for the planetary solutions and $\log q > -1.5$ for the binary solutions, while for KMT-2025-BLG-1392, we find $\log q \sim -1.3$. We show that for the analysis of KMT-2025-BLG-1314, Hamiltonian Monte Carlo (HMC), enabled by microlux, provides robust parameter inference and outperforms traditional Markov chain Monte Carlo (MCMC) methods in the presence of bimodal posteriors.

We benchmark the semi-analytical isothermal Jeans model against a high-resolution isolated N-body simulation that follows a self-interacting dark matter (SIDM) halo into deep core collapse. The model accurately reproduces the density evolution through much of the collapse phase, although it does not capture the sharp rise in central velocity dispersion during collapse. When applied to strong gravitational lensing observables, such as the projected mass and logarithmic density slope of SIDM halos, the Jeans model tracks the simulated evolution more closely than the parametric approach in the deep collapse regime. Our results demonstrate that the isothermal Jeans model provides a reliable and computationally efficient description of SIDM halo evolution.

Quantifying tidal effects on giant planets has recently made significant advances, thanks in particular to the Cassini space probe. During its thirteen-year orbit around Saturn, numerous measurements from different instruments made it possible to characterize fundamental parameters such as Saturnś Love number $k_2$ and quality factor $Q$ at different frequencies. In this article, we summarize the various measurements and methods that have allowed to arrive at such a result, as well as the extrapolations that can be deduced for other systems. More generally, the state of the art concerning the four giant planets of the Solar System is presented, as well as the case of exoplanets.

Giant planets and brown dwarfs play a crucial role in star and planet formation, as they are situated at the boundary between planets and stars with uncertain formation mechanisms. Previous observational searches for the formation boundary were hampered by the lack of large unified samples of wide-orbit giant planets and substellar companions. A combined analysis of radial velocity and astrometry mitigates this problem and has significantly enlarged the sample. Here we present a rigorous statistical analysis of the sample of 55 giant planets, brown dwarfs and low-mass stellar companions orbiting FGK stars. We quantitatively analyze the occurrence rates of brown dwarfs and identify a distinct brown dwarf desert at approximately $30\,M_\mathrm{J}$, with no evidence of disappearance up to 20 au. Unlike previous studies that predicted a declining planet occurrence rate beyond the water-ice line, we identify a new population of giant planets and low-mass brown dwarfs in this region. The metallicity and eccentricity trends in our sample suggest that these are the consequences of two different formation scenarios. Our combined population synthesis model successfully accounts for the observed brown dwarf desert, supporting the dual formation hypothesis.

We present new long-slit stellar kinematics for a sample of 21 Brightest Cluster Galaxies (BCGs) and triaxial Schwarzschild models for 16 of these objects using our orbit modelling code SMART. The new kinematics obtained with the Large Binocular Telescope (LBT) is complemented with high-resolution photometry from HST or new AO-assisted ground-based observations also obtained at LBT and combined with wide-field imaging from the Wendelstein Observatory. These data enable robust modeling from the innermost regions - where the Supermassive Black Hole dominates the potential - to larger radii, where stars and dark matter (DM) are the primary mass contributors. As already discussed in a companion paper, we discovered 8 Ultramassive Black Holes (UMBHs, with mass $> 10^{10}$ M$_\odot$) in this BCG sample, more than doubling the number of galaxies with dynamically detected UMBHs. We show that the DM halos display a wide variety of geometries. Purely kinematical results include low central velocity dispersion with increasing profiles towards the outskirts, and the discovery of one Kinematically Decoupled Core.

Scattering polarization signals in the Sr I 4607 Å spectral line are among the strongest originating from the solar photosphere, offering a powerful diagnostic of tangled magnetic fields in the 3--300 G range via the Hanle effect. However, measuring them with sub-arcsec resolution remains a significant challenge. We analyze spatially resolved quiet-Sun observations of these signals performed with the Visible Spectropolarimeter (ViSP) at the Daniel K. Inouye Solar Telescope (DKIST) and identify its current observational limits. We present high-resolution, high-precision spectropolarimetric observations in a spectral window including the Sr I 4607 Å line at various limb distances. We apply consistent instrumental corrections across all spectral lines, enabling the adjacent lines to serve as reliable references. At a limb distance of $\mu = 0.74$, the signal-to-noise ratio is low but sufficient in the total linear polarization map to directly reveal sub-arcsec structures in the Sr I line for the first time, which can be attributed to scattering polarization. Disk-center measurements are still dominated by noise related to the current limitations of the observational setup. By combining high spatio-temporal and spectral resolution with exceptional polarimetric precision, DKIST enables measurements of solar photospheric scattering polarization at fine scales. However, current signal-to-noise limitations still hinder direct detection of disk-center scattering polarization and must be addressed before further progress can be made.

Gamma-ray bursts (GRBs) are among the most energetic explosive phenomena in the universe, and their peak energy ($E_{\rm p}$) is a key physical quantity for understanding the prompt emission mechanism. However, due to the limited energy coverage of the Swift satellite, a large fraction of Swift GRBs lack reliable measurements of the peak energy. Therefore, developing an accurate and efficient method to predict $E_{\rm p}$ is of great importance. In this work, we propose a method based on the SuperLearner framework that integrates multiple supervised machine learning algorithms to predict $E_{\rm p}$ of Swift/BAT GRBs. We use the Swift/BAT observational data from December 2004 to September 2022 as training features, and adopt the peak energies of 516 GRBs jointly detected by Swift and either Fermi/GBM or Konus-Wind as training labels. After training and testing multiple supervised models, the final SuperLearner ensemble yields a more robust and reliable predictive model. In 100 iterations of 5-fold cross validation, the predicted $E'_{\rm p}$ values show a tight correlation with the observed $E_{\rm p}$, with an average Pearson correlation coefficient of $r = 0.72$. Compared with previous Bayesian estimates, our model provides predictions that are likely closer to the true values. Based on the trained model, we further predict the peak energies of 650 Swift GRBs, significantly increasing the number of GRBs with known peak energies and providing new statistical support for constraining GRB emission mechanisms and energy origins.

We report spatially resolved measurements of cool gas traced by Mg II and Fe II absorption in the circumgalactic medium (CGM) of a star-forming galaxy at $z\sim1$ (G1). The fortuitous alignment of a background gravitational arc at z$\sim$2.4 provides seven closely spaced ($\sim$6 kpc) transverse sightlines along the minor axis of G1, probing its CGM out to $\sim$50 kpc. This geometry allows us to detect a galactic-scale outflow simultaneously in down-the-barrel and transverse directions, where blue-shifted Mg II absorption is detected along both types of sightlines, revealing a large-scale, collimated wind. We measure blue-shifted line-of-sight velocities of $v_{\mathrm{los}}$ $\sim$ 62 - 239 km s$^{-1}$ and line-of-sight velocity dispersions $\sigma_{\mathrm{los}}$ $\sim$ 53 - 133 km s$^{-1}$, suggesting a structure dominated by bulk motion. De-projection of $v_{\mathrm{los}}$ along the minor axis indicates that the outflow material barely approaches the escape velocity and is likely to be gravitationally bound to G1. We constrain an outflow opening angle $\theta_c\sim$ 18$^\circ$ - 25$^\circ$, and a mass outflow rate of $ \dot{M}_{\mathrm{out}}$ $\gtrsim$ 0.06 $M_\odot$ yr$^{-1}$, corresponding to a mass loading factor $\eta$ $\gtrsim$ 0.004, estimated within $\sim$10 - 50 kpc ($\sim$ 0.05 - 0.3 $R_\text{vir}$) of the galaxy centre. Our measurements, combined with previous arc tomography data along the major axis, indicate that normalizing impact parameters by galaxy B-band luminosity substantially reduces scatter in the established anti-correlation between Mg II equivalent width and impact parameter, while also diminishing possible excess of Mg II equivalent width towards the minor axis.

The protosolar helium mass-fraction is a key ingredient of solar, planetary models and enrichment laws. However, the assumed values often rely on simplified descriptions of the transport of chemicals in solar models. They are also based on the inferred helium mass fraction in the solar convective envelope, which is itself sensitive to uncertainties in the solar equation of state. We update the reference protosolar helium abundance by including the effects of macroscopic mixing at the base of the convective zone and more recent determinations of the helium mass fraction in the convective envelope. We combine results from our inversions to spectroscopic abundances, as well as literature values to provide a robust interval of the current helium mass fraction in the convective zone. We combine this measurement to models including light element depletion to provide an udpated protosolar helium abundance. We show that macroscopic mixing at the base of the envelope is key to infer protosolar helium. We find a revised interval of primordial helium mass fraction of 0.27575 +- 0.00315 slightly lower than previous estimates when combining our latest estimate of surface helium mass fraction and spectroscopic abundances. We find that the effects of macroscopic mixing are partially compensated by an increase in the inferred solar helium mass fraction in recent studies. We also derive more precise estimates based on various reference works in the litterature. Using the usual surface helium mass fraction, the primordial helium mass fraction drops to 0.2669 +- 0.00415 due to the inclusion of macroscopic mixing. The dominant source of uncertainty is found the surface helium abundance inferred from helioseismic constraints and more specifically, the impact on the equation of state of the solar material on this inference result.

Dust polarization induced by aligned non-spherical grains acts as an important tool to trace the magnetic field (B-field) morphologies and strengths in molecular clouds and constrain grain properties and their alignment mechanisms. The widely accepted grain alignment theory is the alignment induced by RAdiative Torques (RATs). In this work, we investigate grain alignment mechanisms in a massive, quiescent and filamentary Infrared Dark Cloud G16.96+0.27 using thermal dust polarization observation with JCMT/POL-2 at 850 $\mu$m. We observe the so-called phenomenon of polarization hole attributed to the decrease in polarization fraction in denser regions of higher total intensity and gas density. Our study finds that B-field tangling effect is minimal to cause the polarization hole, and the dominant factor is the reduction in grain alignment efficiency in denser regions, consistent with RAT mechanism. To test RAT theory, we calculate various quantities describing grain alignment, including minimum size of aligned grains, magnetic and magnetic relaxation parameter, and show that RAT mechanism can explain observational data. Our study also reveals evidence for magnetically-enhanced RAT (M-RAT) mechanism required to explain the observed high polarization fractions of above 10 % in the outer regions of the filament. Finally, we perform detailed modeling of thermal dust polarization using $\mathrm{DustPOL\_py}$ based on M-RAT theory and find that the modeling could successfully reproduce the observational data when maximum grain size is around 0.45 $\mu$m accompanied by an increase in grain axial ratio, along with the consideration of variations in the magnetic field's inclination angle with the line of sight.

Counter-Rotating Disk (CRD) galaxies have two co-spatial stellar disks rotating in opposite directions, and provide a rare opportunity to study external gas accretion and dynamical assembly processes. In the 16th data release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, only 64 CRDs were visually identified. Using this as a training sample, we developed an automated pre-selection method that reduces the number of galaxies requiring visual inspection by removing systems unlikely to host counter-rotation. Applying this method to MaNGA Data Release 17, we identified 126 confirmed CRDs and an additional 143 candidate galaxies, more than doubling the MaNGA CRD sample. With this extended sample, we analyzed their Baldwin, Phillips, and Terlevich (BPT) emission-line diagrams and compared them with a matched control sample of early-type galaxies (ETG). We found no statistically significant difference in photoionization sources between CRDs and the ETG control sample, indicating emission-line diagnostics cannot solely be used to identify counter-rotating disks, nor do they correspond to a distinct present-day photoionization signature. Our method facilitates efficient discovery of CRDs in large spectroscopic surveys, enabling stronger statistical studies of their formation and evolution.

In this work, we investigate the General Relativistic Entropic Acceleration (GREA) framework, in which late-time acceleration emerges from entropy production associated with the cosmological horizon, and compare its performance with the standard $\Lambda$CDM description of the Universe. We first confront GREA with current background observations, including baryon acoustic oscillations, type Ia supernovae, compressed CMB information, and cosmic chronometers, with particular emphasis on the geometric horizon parameter $\sqrt{-k}\eta_0$. We then introduce a phenomenological extension of the theory by allowing for an additional dark energy component, $\Omega_{de}$, enabling the recovery of a $\Lambda$CDM-like expansion history as a limiting case. We perform a Bayesian parameter inference and model comparison analysis using both current data and mock datasets representative of future surveys, including SKAO, LSST, and ET. While current data statistically prefer $\Lambda$CDM when compressed CMB information is included, GREA remains competitive for low-redshift combinations. Forecasts indicate that gravitational wave standard sirens are expected to enhance the ability to discriminate between entropic-driven and dark-energy-driven expansion scenarios, and to identify the underlying cosmological model favored by the data.

Atomic carbon ([CI]) is a key species in the carbon chemistry of the interstellar medium (ISM). Using the Submillimeter Wave Astronomy Satellite (SWAS), we conducted a [CI]($^3$P$_1$--$^3$P$_0$) 492 GHz survey covering approximately 4 deg$^2$ of the L1688 and L1689 regions in the $\rho$ Oph molecular cloud, achieving a spatial resolution of 4.25$\hbox{$^{\prime}$}$. The derived [CI] column densities, N([CI]), range from 4.85 $\times$ 10$^{14}$ cm$^{-2}$ to 6.29 $\times$ 10$^{17}$ cm$^{-2}$, corresponding to an abundance ratio N([CI])/N($H_2$) of 2.24$\times$ 10$^{-7}$ to 2.39$\times$ 10$^{-4}$, with a median value of 1.8$\times$ 10$^{-5}$. Combining observations with photodissociation region (PDR) modeling, we find that [CI] abundance varies less than CO in regions with UV intensity G$_0$ $> 16$ and N(H$_2$) $<$ 4.6 $\times$ 10$^{21}$ cm$^{-2}$, suggesting [CI] is a more reliable tracer of molecular hydrogen in low-density, high-radiation environments where the [CI]-to-CO transition occurs. Utilizing [CI] as direct H$_2$ tracer, the CO-dark gas fraction is estimated to be 0.43 , meaning that 43% of the total cloud mass will be missed by conventional calculation based on CO observations but can be calibrated by [CI] emission. The [CI] line widths are systematically broader than those of $^{13}$CO, possibly due to contributions from atomic carbon. These findings provide key insights into Galactic [CI] emission and the carbon cycle evolution in the interstellar medium. Future high-sensitivity [CI] ($^3$P$_1$--$^3$P$_0$) surveys with the Chinese Survey Space Telescope (CSST) will significantly advance our understanding of the carbon cycle evolution.

We prove the absence of scale-invariant one-loop corrections to the superhorizon curvature perturbations from small-scale (potentially enhanced) scalar perturbations in a general inflationary setup, including the transient ultra-slow-roll scenario. We demonstrate this by analyzing the symmetry structure of an in-in effective field theory for the soft curvature perturbations, and by explicitly performing one-loop calculations, integrating out hard modes in the soft limit of external momenta. The dilatation symmetry, respected by a counter term necessary for the tadpole cancellation, guarantees the cancellation of scale-invariant corrections.

Upcoming imaging surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will enable high signal-to-noise measurements of galaxy clustering. The halo occupation distribution (HOD) is a widely used framework to describe the connection between galaxies and dark matter haloes, playing a key role in evaluating models of galaxy formation and constraining cosmological parameters. Consequently, developing robust methods for estimating this statistic is crucial to fully exploit data from current and future galaxy surveys. The main goal of this project is to extend a background subtraction method to estimate the HOD with more photometry-based information in preparation for the clustering analysis of the upcoming LSST data and to enable the study of the HOD with significantly improved statistical power. We evaluate the performance of the method using a mock galaxy redshift survey constructed from the cosmoDC2 catalogue. We implement an extension of the background subtraction technique to utilize information from photometric galaxy surveys. To identify the centres of galaxy groups, we implement an iterative centroiding approach (Central Galaxy Finder). We evaluate the impact of each step in our pipeline, including group size estimation from luminosity and purity, and completeness on group identification, along with the influence of observational systematics such as the use of photometric redshifts and halo mass uncertainties. We demonstrate the validity of the proposed method using a mock galaxy catalogue, recovering the HOD from cosmoDC2 over the absolute magnitude range $M_r = -20.0$ to $-17.0$ and halo masses up to $10^{15}\, \mathrm{M_\odot}$. We present key performance metrics to quantify the precision and reliability of the group finder and the resulting HOD measurements.

After the recent detection of solar equatorial Rossby waves, a renewed interest has been brought to the study of gravito-inertial waves propagating in the convective envelope of solar-type stars. In particular, the ability that some of these envelope gravito-inertial modes have to couple with the ones trapped in the radiative interior might open new windows to probe the deep-layer dynamics of solar-type stars. The possibility for such a coupling to occur is particularly favoured in pre-main sequence (PMS) solar-type stars. Indeed, due to the contraction of the protostellar object, they are able to reach large rotation frequencies before nuclear reactions are ignited and magnetic braking becomes the driving mechanism for their rotational evolution. In this work, we therefore study the coupling between the envelope inertial waves and the radiative interior g modes in PMS stars, focusing on the case of prograde dipolar modes. We consider the case of 0.5 Msun and 1 Msun PMS models, each with three different scenarios of rotational evolution. We show that, for stars that have formed with a sufficient amount of angular momentum, this coupling can occur in frequency ranges that are accessible to space-borne photometry, creating inertial dips in the period spacing pattern. With an asymptotic analysis we characterise the shape of these inertial dips to show that they depend on rotation and on the stiffness of the convective-radiative interface.

We present an analytical instrument model of the TMS radiometer, a pseudo-correlation system designed for absolute sky-temperature measurements through a continuous comparison between sky and reference-load signals. The goal of this work is to quantify and understand the impact of instrumental non-idealities that are intrinsic to absolute radiometric measurements. We used a combination of the Jones-matrix formalism, to describe non-ideal signal mixing effects in the different instrument components, and the Friis formalism, to account for noise-temperature contributions which introduce specific offset terms. The model includes all components from the cryostat entrance window onwards and allows us to propagate realistic losses, return losses, and noise figures through the system. We find that non-ideal mixing effects, such as intensity-to-polarization and sky-to-load leakage, are expected to appear at the percent level, but they could be in principle calibrated out. The total TMS intensity response exhibits a frequency-dependent offset with a band-averaged level of 6.9 K, dominated by the cryostat window, and with smaller contributions from the infrared filter and orthomode transducers. Under the required TMS thermal stability of 1 mK over one hour, the resulting variation in the output signal can reach up to 91.3 mK, which sets the fundamental limit on the absolute sky-temperature accuracy. In contrast, relative spectral measurements across the 10-20 GHz band are stable at the few-microkelvin level, consistent with the instrument target sensitivity of ~ 10 Jy/sr. This work provides a detailed and quantitative assessment of the systematic effects affecting absolute radiometric measurements and establishes a robust framework for the calibration and performance optimization of the TMS instrument.

We present a multi-band study of three symbiotic binaries using combined ground- and space-based monitoring that spans up to 14 years. These datasets enable a systematic investigation of variability on intermediate timescales (tens of days) and the detection of shorter-period signals. All systems display coherent photometric modulations that are distinct from the orbital cycles. In AX Per, a dominant 75-day signal and its 37-day harmonic are identified, which we interpret as pulsations of the cool giant. CI Cyg exhibits a stable modulation between 70 and 74 days, which likely arises from a combination of pulsation and circumstellar or disk-related variability. For Z And, we confirm a persistent modulation between 55 and 60 days, consistent with semiregular pulsations of the cool component. Additionally, space-based data reveal further short-period variability, including coherent signals at 26.7 and 66.6 minutes in Z And and CI Cyg, respectively, and a quasi-periodic modulation near 0.95 days in AX Per. These detections suggest the presence of rapid activity driven by accretion or rotation, superposed on the intermediate timescale behavior. Our results show that the observed variability in these symbiotic binaries reflects the combined effects of cool giant pulsation, circumstellar or disk activity, and possible rotation of the hot component. The multi-timescale behavior revealed here offers new constraints on mass transfer and activity cycles in interacting binaries.

How substructures and disk properties affect dust evolution and the delivery of solids and volatiles into planet-forming regions remains an open question. We present results from tailored dust evolution modeling of the AGE-PRO ALMA large program, a sample of 30 protoplanetary disks spanning different evolutionary stages. Visibility fitting of the AGE-PRO ALMA data (at 1.3\,mm) reveals that approximately half of the disks exhibit radial substructures. Combined with stellar properties, disk inclinations, and gas mass estimates from CO isotopologues and N$_2$H$^+$, this well-characterized set of disks provides an ideal testbed to constrain dust evolution models across different ages and disk morphologies. Using the dust evolution code \texttt{DustPy}, we simulate dust evolution in each disk under four model configurations, varying two key free parameters: the turbulent viscosity ($\alpha = 10^{-4}, 10^{-3}$) and fragmentation velocity ($v_{\rm{frag}} = 1 \mathrm{m\,s^{-1}}, 10 \mathrm{m\,s^{-1}}$). Pressure traps are incorporated by perturbing the gas surface density based on the continuum intensity profiles, and synthetic observations generated with \texttt{RADMC-3D} are compared to these profiles. While no single model fits all disks, nearly half are best reproduced by the configuration with low turbulence and low fragmentation velocity ($\alpha = 10^{-4}, v_{\rm{frag}} = 1\,\mathrm{m\,s^{-1}}$). Models of smooth disks underpredict dust mass, possibly indicating unresolved substructures. Pebble fluxes into inner disk regions correlate more strongly with disk age than with the presence of substructures, highlighting time-dependent dust transport as a key factor in shaping inner disk composition. Our results also provide a comparative baseline for interpreting multiwavelength and JWST water vapor observations.

The cross spectrum encodes the correlated variability between two time signals. In recent years, the cross spectrum has been used to study astronomical sources, particularly in the field of X-ray timing. In the literature, it has been common to either simultaneously fit the real and imaginary components of the cross spectrum, or fit the phase and magnitude. Until now, a full discussion of the statistical distribution of the cross spectrum has been missing from the astronomical literature. In this paper, we present a derivation of the full statistical distribution of a cross spectrum between two time series, showing that it follows an asymmetric Laplace distribution. We further provide the probability distribution function for a cross spectrum random variable, along with the marginal distributions for many quantities. We also relate the cross spectrum to the power spectra of the constituent time series. This work will enable the cross spectrum to be used more accurately as a probe of physical processes such as accretion onto black holes and neutron stars.

Recent PTA analyses reporting evidence for a nanohertz common-spectrum process motivate targeted tests of whether any anisotropic component of the stochastic gravitational-wave background (SGWB) is correlated with the nearby large-scale structure (LSS), as anticipated for an astrophysical background dominated by supermassive black hole binaries. We present the first Bayesian PTA likelihood analysis that embeds an externally observed, full-sky galaxy-survey LSS template directly as an overlap-reduction-function (ORF) component. Using the 2MASS Photometric Redshift (2MPZ) galaxy catalog, we construct low-multipole LSS--correlated ORF templates in two redshift slices ($0<z\le0.1$ and $0.1<z\le0.2$) and model PTA cross-correlations as $\Gamma_{ab}=\Gamma^{\rm HD}_{ab}+\sum_i \epsilon_i\,\Gamma^{\rm LSS(i)}_{ab}$, where $\epsilon_i$ quantifies the amplitude of an SGWB component whose angular correlations project onto the fixed 2MPZ LSS templates. Applying this framework to the NANOGrav 15-year dataset, we find no statistically significant evidence for an LSS-correlated component: $\epsilon_i$ is consistent with zero in both single-bin and two-bin analyses (e.g., $\epsilon_1=0.20^{+1.68}_{-1.66}$ and $\epsilon_2=-0.11^{+2.04}_{-1.83}$; 68\% credible intervals), and Bayes factors favor the isotropic Hellings--Downs hypothesis ($\mathcal{B}_{{\rm HD+LSS}_1,{\rm HD}}=0.40$, $\mathcal{B}_{{\rm HD+LSS}_2,{\rm HD}}=0.43$, $\mathcal{B}_{{\rm HD+LSS}_{1+2},{\rm HD}}=0.11$). We therefore place upper limits on any 2MPZ-traced, LSS-correlated contribution to the SGWB at $z<0.2$. More broadly, our framework provides a reproducible pathway for incorporating observed LSS information into PTA anisotropy searches and naturally motivates extensions to finer redshift tomography and next-generation PTA datasets.

We construct constant rest-mass sequences of equilibrium models of differentially rotating neutron stars which resemble binary neutron star post-merger remnants. For a more realistic description of the post-merger remnant, we impose that each model carries approximately $95\%$ of the angular momentum that a binary system with the same total rest-mass has at the moment of merging, based on an empirical relation informed from neutron star merger simulations. We account for equation of state effects by employing two distinct microphysical descriptions for high density matter. We dynamically evolve the equilibrium models with a three-dimensional general relativistic hydrodynamics code that employs the conformal flatness approximation. We investigate the connection between the occurrence of the instability and the existence of corotation radii within the stellar configurations and determine the instability window for both equation of state sequences. The occurrence of low-$T/|W|$ instabilities leads to pronounced gravitational wave emission in the range $0.13 \lessapprox\beta \lessapprox 0.2$, while models outside this range exhibit less pronounced features in the gravitational wave spectrum. The prominence of gravitational wave emission is primarily determined by $\beta$, while the equation of state seems to have a more minor effect. We present correlations between the strength of the gravitational wave emission associated with the instability and properties of the equilibrium models. Stellar configurations modelled by different equations of state display differences in the timescales over which the various dynamical features develop, as well as whether they exhibit a pronounced $m=1$ deformation. Potential relations between the instability growth timescales and properties of the stellar models are studied.

Map-space $E$/$B$ decompositions of linear polarization are attractive for foreground and CMB analyses because they isolate the $B$-family patterns that contaminate primordial tensor searches from $E$-family patterns that trace coherent Galactic structures. However, the $E$/$B$ transform is non-fully-local and induces apparent spectral complexity in projected fields even when the underlying sky is spectrally simple in $\underline{P}=Q+iU$. We quantify this effect for synchrotron emission. We introduce a complex-parameter description of the frequency dependence of $\underline{P}$, its spin-preserving projections $\underline{P}_E$ and $\underline{P}_B$, and the scalar $\underline{S}=E+iB$, using complex log--Taylor and moment expansions (with simple transformation rules under $E$/$B$ projection) and linking their coefficients to spectral-index variations, line-of-sight mixing, synchrotron ageing, and Faraday effects. Using a toy model and a PySM template, we find that scalar combinations, especially $|E|$ and $|B|$, acquire the largest induced complexity, while $\underline{S}$ is less affected but lacks a directly interpretable amplitude and angle. By contrast, $\underline{P}_E$ and $\underline{P}_B$ retain a clear geometric meaning and exhibit only moderate spectral distortions, while satisfying the closure relation $\underline{P}=\underline{P}_E+\underline{P}_B$ (which extends to all spectral orders in the moment formalism). Finally, with three frequency channels, we compare low-order spectral truncations and propose diagnostics to test whether the data favour a single power law in $P$ or independent power laws in $(P_E,P_B)$. This work is intended to be of practical relevance for both Galactic science and CMB $B$-mode analyses and lays the conceptual foundation for a series of papers applying the framework to observational data.

The coming data releases of Gaia are expected to result in an upheaval of exoplanet science, in particular for long period giant planets ($0.2 {\rm M_{J}} \leq M\leq25 {\rm M_{J}}$). One class of exoplanets which Gaia will help investigate is circumbinary planets. Using the current knowledge of the circumbinary exoplanet population as well as expectations for the Gaia sensitivity, we investigate the impact Gaia will have on our understanding of circumbinary planets. We compare our results to a pre-launch estimate, the main differences arising from a better understanding of the circumbinary planet population, which result in a lower expected yield than previously predicted, though still significant compared to the known population. We make a rough yield estimate, with conservative detection criteria and parameter-space cuts, predicting in the 10s - 100s of detections in Gaia DR4. More importantly, we show how the yield estimate varies strongly with different assumptions on the injected circumbinary population, showing Gaia's sensitivity to the mass and orbital period distribution of circumbinary planets. We find that Gaia circumbinary exoplanet detections will be biased towards planets closer to the instability zone surrounding the binary, due to the larger number of binaries on wider orbits and the limited timespan of Gaia. We also assess the impact Gaia will have on known circumbinary systems, one being that it may resolve the question of reliability of the claimed planets orbiting post-common-envelope binaries, with Gaia DR5 being sensitive to between 3 and 11 out of 32 such planet candidates.

Atmospheric characterization of exoplanets using ground-based high-resolution transmission spectroscopy has traditionally focussed on large and close-in planets, such as hot Jupiters. In this work, we aim to extend this technique to smaller and more temperate planets by studying the atmospheric composition of the temperate super-Earth planet L 98-59 d ($\sim$$1.5\,\mathrm{R_{\oplus}}$; $\sim$$1.9\,\mathrm{M_{\oplus}}$). Using high-resolution transmission spectra obtained using IGRINS on the Gemini-South telescope, we demonstrate the feasibility for atmospheric characterization of super-Earths using ground-based facilities, and confirm the previous tentative JWST inference of hydrogen sulfide (\ce{H2S}) in the atmosphere of L 98-59 d at $\lesssim$3.9$\,\sigma$ ($B\sim390$). This is the first ground-based inference of a molecular species in the atmosphere of a super-Earth planet, and reveals the sensitivity of spectrographs on 8\,m-class telescopes to the atmospheric characterization of such planets. By exploring a grid of atmospheric models, we find that the data favors a cloud-free atmosphere with an abundance of H2S corresponding to $\sim$1-10\,$\times$ solar metallicity. We additionally place constraints on the atmospheric abundances of other molecular species. Assuming cloud-free models, super-solar abundances for CH4 and NH3 are ruled out at 3.6$\sigma$ and $4.6\sigma$, respectively. Our results are consistent with previous suggestions that L 98-59 d is a super-Earth with possible disequilibrium production of \ce{H2S} driven by volcanic outgassing from the surface. Future studies combining multiple observations with different facilities may be able to further constrain the atmospheric composition of this planet. This work underscores the promise of atmospheric characterization of super-Earth exoplanets using high-resolution spectroscopy with ground-based facilities.

The construction of an accurate 3D model of the Milky Way center is necessary to understand inflow processes that drive its overall evolution, and to compare our Galactic nucleus to other galaxies' nuclei. A main point of contention is the line-of-sight location of sources observed toward the central 10 pc of the Galaxy, including recent star formation (the Sgr A East supernova remnant and Sgr A HII regions) and copious gas (the 50 and 20 km/s molecular clouds, the Circumnuclear Disk, and the Sgr A West ionized "minispiral" that encircles the central supermassive black hole, Sgr A*). Some models place all of these structures within a radius of 5 pc from Sgr A*, while others place the 20 and 50 km/s clouds at a distance of at least 30 - 50 pc away from Sgr A* along the line of sight. We present new radio and millimeter observations of the molecular gas toward the central ~10 pc, from which we have constructed an alternative 3D model that is consistent with both prior radio observations and orbital gas kinematics. Our model places the 20 km/s cloud, 50 km/s cloud, and Sgr A East more than 10 pc in front of Sgr A*. While this model does not conclusively rule out a connection between the 50 and 20 km/s clouds and the circumnuclear disk, we argue that prior evidence for these connections is tenuous, especially given the complex spatial and kinematic overlap of structures along the line of sight.

Implementing the axion concept in the context of 3-3-1 extensions of the Standard Model (SM) leads to richer properties than in the simpler axion setups, and related to the Dirac neutrino seesaw mechanism. In this way the smallness of neutrino masses, the strong CP problem, the nature of dark matter and the number of families all have a common origin. Besides having an enhanced coupling to photons, our revamped axion can also be distinguished from DFSZ and KSVZ axions through its couplings to fermions. The latter lead to interesting phenomenological consequences, including flavor-changing axion-emitting two-body K, B and D meson decays.

We review theoretical frameworks in which small neutrino masses arise radiatively through interactions with a dark sector that also accounts for cosmological dark matter (DM). A prototype is provided by scotogenic schemes, that extend the inert Higgs doublet model to include dark fermions. We outline their key features and limitations, discussing the advantages of the revamped scotogenic extension. The phenomenological signatures of fermionic and bosonic scotogenic dark matter are discussed, along with scoto-seesaw models that merge scotogenic and seesaw mechanisms. We also consider scenarios where the dark sector seeds a low-scale seesaw. These frameworks can accommodate dark matter as Weakly or Feebly Interacting Massive Particles (WIMPs or FIMPs). While hidden dark sector models are inherently difficult to exclude, visible dark sector schemes should be confirmed--or ruled out--by forthcoming dark matter, collider, and lepton flavor violation studies.

Bayesian evidence evaluation becomes computationally prohibitive in high dimensions due to the curse of dimensionality and the sequential nature of sampling-based methods. We introduce SunBURST, a deterministic GPU-native algorithm for Bayesian evidence calculation that replaces global volume exploration with mode-centric geometric integration. The pipeline combines radial mode discovery, batched L-BFGS refinement, and Laplace-based analytic integration, treating modes independently and converting large batches of likelihood evaluations into massively parallel GPU workloads. For Gaussian and near-Gaussian posteriors, where the Laplace approximation is exact or highly accurate, SunBURST achieves numerical agreement at double-precision tolerance in dimensions up to 1024 in our benchmarks, with sub-linear wall-clock scaling across the tested range. In multimodal Gaussian mixtures, conservative configurations yield sub-percent accuracy while maintaining favorable scaling. SunBURST is not intended as a universal replacement for sampling-based inference. Its design targets regimes common in physical parameter estimation and inverse problems, where posterior mass is locally well approximated by Gaussian structure around a finite number of modes. In strongly non-Gaussian settings, the method can serve as a fast geometry-aware evidence estimator or as a preprocessing stage for hybrid workflows. These results show that high-precision Bayesian evidence evaluation can be made computationally tractable in very high dimensions through deterministic integration combined with massive parallelism.

First-order cosmological phase transitions (PT) can take place in a dark sector at relatively late times between the big-bang nucleosynthesis and recombination epochs. Because bubble nucleation is stochastic, the PT completes at different times in different regions of the Universe. This fluctuation sources a curvature perturbation whose (dimensionless) power spectrum ${\cal P}_\zeta(k)$ features a universal infrared tail, independent of the microscopic details of the PT. Even in the absence of any non-gravitational interaction between the dark sector and the Standard Model, these additional curvature perturbations at small scales impact a variety of observables. We derive new constraints on dark sector phase transitions using {\it Planck}, baryon acoustic oscillation (BAO), Lyman-$\alpha$ observations, spectral distortion limits from FIRAS, constraints on early reionization, and the existence of ultra-faint dwarf galaxies.

Cosmological time dilation (CTD) serves as a fundamental probe of cosmic expansion, historically verified through the characteristic (1+z) broadening of Type Ia supernova (SNe Ia) light curves. However, significant tensions arise when extending this test to other astrophysical regimes. While discrete, event-based transients such as Gamma-Ray Bursts (GRBs) exhibit large scatter in interred time-dilation signatures, analyses of stochastic variability in persistent sources, specifically Quasars (QSOs), frequently yield null results. I demonstrate that these discrepancies stem from a previously overlooked distinction between discrete geometric clocks and continuous thermal emission, presenting a resolution within the framework of Generalized Cosmological Time (GCT). The central premise relies on strictly distinguishing global coordinate time, characterized by a generalized lapse function, from the local proper time measured within gravitationally bound systems. We propose that the progenitors of transients, specifically SNe Ia and GRB central engines, are effectively shielded from background time evolution due to strong gravitational binding and environmental decoupling. Consequently, they act as standard clocks tracing pure geometric path dilation, obeying \tau_{\rm obs} \propto (1+z)^{1+b/4}. Conversely, the lack of dilation in QSOs is derived as a consequence of observing persistent thermal accretion disks at fixed wavelengths, introducing an intrinsic selection effect (\tau_{\rm intr} \propto (1+z)^{-2}) that masks the cosmological signal. This framework reconciles the diverse behaviors of transient and persistent sources without modifying local physical laws.

We investigate Extended Geometric Trinity of Gravity at both classical and quantum cosmological levels using the minisuperspace approach. Adopting Noether symmetries to select viable models, we examine metric-affine theories of gravity, in particular the extensions of General Relativity, Teleparallel Equivalent General Relativity and Symmetric Teleparallel Equivalent General Relativity, and show that the equivalence among these different formulations can be restored by including in the Lagrangian the divergence terms that relate their respective geometric invariants to the Ricci scalar. Exact cosmological solutions are derived and compared in the different models.

Laser-ranged satellites have demonstrated exceptional effectiveness in high-precision verification of General Relativity but also in the accurate inversion of geophysical parameters such as Earth tidal parameters. Due to the extremely weak frame-dragging signal, after utilizing the orbital symmetry of LARES 2 and LAGEOS satellites to cancel most of the nodal precession caused by Earth's oblateness, precise modeling of incompletely symmetric Earth tidal perturbation patterns becomes the core prerequisite for effectively extracting this signal. This study, based on Kaula's perturbation theory and Lagrange equations, investigates the perturbations of LARES 2 and LAGEOS satellite orbital nodes and inclinations caused by Earth tides. From the tidal perturbations computed from 402 earth tide constituents, 83 tidal perturbations with significant amplitudes were selected by threshold based on the RMS of overlap orbit differences of the two satellites, the asymmetric characteristics of tidal perturbations between the two satellites were quantitatively analyzed. The traditional amplitude threshold method for individual tidal perturbations is limited, as the coherent superposition of minor tidal constituents and frequency-orbit resonance lead to the total effect exceeding the threshold. These results provide important support for high-precision tests of General Relativity and the refinement of satellite orbital dynamics modeling.

We study single-field inflation in which the inflaton is identified with the lightest scalar (dilaton) excitation of a confining gauge theory. The inflaton potential is not postulated: it follows from the pure effective Gluodynamics Lagrangian tightly constrained by the trace anomaly and the associated infinite tower of Ward identities, yielding a Coleman--Weinberg form with a logarithmic term fixed by nonperturbative condensates. After coupling to gravity via a non-minimal interaction $\xi\,\varphi^2 R$, the Einstein-frame potential develops a plateau consistent with current CMB observables. In the large-$\xi$ limit the model approaches the standard plateau attractor, while the Migdal--Shifman(MS) logarithmic structure induces a controlled, testable deformation governed by $A/\lambda$ across the CMB window. We quantify the resulting shifts in $(n_s,r)$ and the running analytically and confirm them with numerical scans over $(\xi,\lambda,A,\mu)$, making the departure from the attractor both microphysically motivated and observationally predictive.

The master student at Lomonosov University in Moscow E.P. Velikhov formulated 1959 the theory of magnetorotational instability, which dominates current astrophysics. A meteoric career made him later the science, nuclear and disarmament advisor to Gorbachev and Yeltsin. This article describes the interactions between Velikhov and the PROMISE team from Potsdam and Dresden-Rossendorf when it came to experimentally testing the theory in the laboratory in the 2000s. At an MHD conference in Catania, Sicily, he offered the public the use of small, transportable Russian nuclear power plants anywhere in the world until the fusion machines currently under development had finally solved the energy problem.

The heating of the solar wind is a key to understand its dynamics and acceleration process. The observed radial decrease of proton temperature in the solar wind is slow compared to the adiabatic prediction and it is thought to be caused by turbulent dissipation. To generate the observed 1/R decrease, the dissipation rate has to reach a specific level which varies in turn with temperature, wind speed, and heliocentric distance. We want to prove that MHD turbulent simulations can lead to the 1/R profile. We consider here the slow solar wind, characterized by a quasi-2D spectral anisotropy. We use the EBM (expanding box model) equations, which incorporate into 3D MHD equations the expansion due to the mean radial wind, allowing to follow the plasma evolution between 0.2 and 1 AU. We vary the initial parameters which are: Mach number, expansion parameter, plasma beta, and properties of the energy spectrum as the spectral range and slope. Assuming turbulence starts at 0.2 AU with a Mach number equal to unity, with a 3D spectrum mainly perpendicular to the mean field, we find radial temperature profiles close to 1/R in average. This is done at the price of limiting the initial spectral extent, corresponding to the small number of modes in the inertial range available, due to the modest Reynolds number reachable with high Mach numbers.

Aims. We investigated plasma turbulence in the context of solar wind. We concentrated on properties of ideal second-order magneto-hydrodynamic (MHD) and Hall MHD invariants. Methods. We studied the results of a two-dimensional hybrid simulation of decaying plasma turbulence with an initial large cross helicity and a negligible magnetic helicity. We investigated the evolution of the combined energy and the cross, kinetic, mixed, and magnetic helicities. For the combined energy and the cross, kinetic, and mixed helicities, we analysed the corresponding Kármán-Howarth-Monin (KHM) equation in the hybrid (kinetic proton and fluid electron) approximation. Results. The KHM analysis shows that the combined energy decays at large scales. At intermediate scales, this energy cascades (from large to small scales) via the MHD non-linearity and this cascade partly continues via Hall coupling to sub-ion scales. The cascading combined energy is transferred (dissipated) to the internal energy at small scales via the resistive dissipation and the pressure-strain effect. The Hall term couples the cross helicity with the kinetic one, suggesting that the coupled invariant, referred to here as the mixed helicity, is a relevant turbulence quantity. However, when analysed using the KHM equations, the kinetic and mixed helicities exhibit very dissimilar behaviours to that of the combined energy. On the other hand, the cross helicity, in analogy to the energy, decays at large scales, cascades from large to small scales via the MHD+Hall non-linearity, and is dissipated at small scales via the resistive dissipation and the cross-helicity equivalent of the pressure-strain effect. In contrast to the combined energy, the Hall term is important for the cross helicity over a wide range of scales. The magnetic helicity is scantily generated through the resistive term and does not exhibit any cascade.

Gauss's method of orbit determination (OD) and its variants are among the most popular initial state estimation techniques for astronomers and engineers alike. However, owing to its assumptions regarding the two-body problem, Gauss's method is inapplicable in the cislunar domain, where three body effects dominate. We introduce a hybrid Particle Gaussian Mixture filtering method, a purely recursive probabilistic orbit determination framework based on a combination of the Markov Chain Monte Carlo based Particle Gaussian Mixture-II (PGM-II) and Particle Gaussian Mixture-I (PGM-I) filters. This method enables us to fuse probabilistic information with angles-only observations from terrestrial telescopes for short and long-term cislunar target tracking. We demonstrate this technique on an important cislunar orbit regime.

Next-generation gravitational-wave detectors will make black-hole ringdown an increasingly sensitive probe of small departures from General Relativity in the strong-field regime. This motivates obtaining high-precision predictions of gravitational effective field theory, as spectral shifts can be quite small. Here we perform a focused stress test of the modified-Teukolsky framework by designing two null diagnostics. First, we consider an action with redundant operators that must produce zero first-order vacuum QNM shifts. Second, we exploit a Ricci-flat identity relating two physical cubic Riemann to test such a relation is satisfied by the ringdown spectra obtained. We compute the shifts using two independent numerical approaches: the eigenvalue-perturbation and generalized continued-fraction (Leaver-type) methods. Both null tests are passed across multiple multipoles and overtones, and the control-operator results agree in magnitude with the benchmark values reported in Ref. [1]. These validations support using the framework for obtaining accurate precitions for robust strong-field tests, with straightforward extensions to rotating backgrounds and coupling with matter fields.

We present $\mathcal{H}\mathtt{-EFTCAMB}$, the official successor to $\mathtt{EFTCAMB}$. The original $\mathtt{EFTCAMB}$ is designed as a consistent and numerically stable implementation of the effective field theory (EFT) of dark energy in the Einstein-Boltzmann code $\mathtt{CAMB}$. On top of this, $\mathcal{H}\mathtt{-EFTCAMB}$ introduces a new Horndeski module that supports computing cosmology for an arbitrary input covariant Horndeski Lagragian. $\mathcal{H}\mathtt{-EFTCAMB}$ supports both mapping the Horndeski theory to an EFT lagrangian to solve in the EFT framework as well as directly solving for the scalar field equations of motion derived from the covariant Lagrangian. The latter approach also works for the cases when the Horndeski field experiences turn-overs, e.g. oscillation, where the EFT approach breaks down. The Horndeski module has been validated by comparing internally with existing models in the original $\mathtt{EFTCAMB}$ and externally with $\mathtt{hi\_class}$. $\mathcal{H}\mathtt{-EFTCAMB}$ features a flexible Python wrapper that is seamlessly integrated into the widely utilized cosmological sampler $\mathtt{Cobaya}$. \heft~is publicly available and serves as a comprehensive tool for testing gravity against the precision data from current and next-generation surveys.

We present ImpCresst, a Geant4-based Monte Carlo tool to simulate backgrounds from natural and cosmogenic radionuclides, and calibration signals in solid-state detectors and their response to it. It is tuned for a fast-evolving and heterogeneous detector environment with a focus on physics at the keV range. This tool was originally developed and validated by the CRESST collaboration; however, its flexibility and configurability make it suitable for other experiments with similar requirements. Key features of ImpCresst include the dynamic geometry implementation directly from CAD files, ROOT-based data persistency of the whole event topology and automatic metadata annotation for data provenance, and interfaces to various particle generators, particularly for radiogenic and cosmogenic radionuclides. It includes also a newly developed particle generator for radioactive bulk and surface contaminations which is completely independent of any user defined confinement volumes. The auxiliary tool CresstDS applies detector-specific energy and time resolution based on a user-provided data set of empirical parameterization. We discuss also how to manage an ImpCresst based workflow in an HPC environment based on Apptainer and nextflow.