Papers for Wednesday, Dec 03 2025

We present a geometrical approach that provides a non-perturbative technique, allowing the standard FLRW observer to evaluate a measurable, scale-dependent distance functional between her idealized FLRW past light cone and the actual physical past light cone. From the point of view of the FLRW observer, gathering data from sources at cosmological redshift $\widehat{z}$, this functional generates a geometry--structure--growth contribution ${\Omega_\Lambda(\widehat{z})}$ to the FLRW cosmological constant ${\widehat\Omega_\Lambda}$. This redshift--dependent contribution erodes the interpretation of ${\widehat\Omega_\Lambda}$ as representing constant dark energy. In particular, ${\Omega_\Lambda(\widehat{z})}$ becomes significantly large at very low $\widehat{z}$, where structures dominate the cosmological landscape. At the pivotal galaxy cluster scale, where cosmological expansion decouples from the local gravitation dynamics, we get ${\Omega_\Lambda(\widehat{z})/\widehat\Omega_\Lambda}\,=\,\mathscr{O}(1)$, showing that late--epoch structures provide an effective field contribution to the FLRW cosmological constant that is of the same order of magnitude as its assumed value. We prove that ${\Omega_\Lambda(\widehat{z})}$ is generated by a scale-dependent effective field governed by structure formation and related to the comparison between the idealized FLRW past light cone and the actual physical past light cone. These results are naturally framed in the mainstream FLRW cosmology; they do not require exotic fields and provide a natural setting for analyzing the coincidence problem, leading to an interpretative shift in the current understanding of constant dark energy.

We report the discovery of two broad-line X-ray AGNs cid_414 and cid_947 at z~3 that exhibit prominent He I$\lambda$10830+ Pa$\gamma$ emission and absorption, identified from the JWST Cycle 3 large GO treasury program COSMOS-3D using NIRCam F444W grism spectroscopy. Additional UV/optical line measurements (e.g., Ly$\alpha$, Si IV, C IV) come from complementary COSMOS-field spectroscopy. Both sources are robustly detected in the mid-infrared, with detections in MIRI F1000W for both AGNs and an additional detection in MIRI F2100W for cid_414, indicating the presence of hot dust emission. The source cid_947 shows a higher He I$\lambda$10830 absorption column density and X-ray-inferred $N_{\rm H}$, and displays strong outflow signatures in He I, Si IV, and C IV with velocity offsets exceeding 5000 km/s. The source cid_414 shows a narrow Ly$\alpha$ emission line with luminosity $\log L_{\rm Ly\alpha}=42.49\pm0.01~\mathrm{erg~s^{-1}}$ and a higher intrinsic 2-10 keV X-ray luminosity. Host-galaxy decomposition and multi-component SED fitting indicate that cid_947 hosts a more massive black hole but lower star formation rate than cid_414. From simplified photoionization modeling, we infer that the dense absorbing gas has a characteristic size comparable to the nuclear broad-line region and is likely kinematically coupled to the obscuration associated with the dust torus. He I$\lambda$1083 absorption has also been identified in several compact little red dots at similar redshifts. Together with the two AGNs reported here, these findings suggest that dense circumnuclear gas are plausibly prevalent at high redshift and plays an important role in regulating AGN obscuration and black hole--host co-evolution.

The "Reduced Speed of Light" (RSL) approximation is commonly used to speed up radiative transfer calculations in cosmological simulations. However, it has been shown previously that the RSL approximation leads to photon non-conservation in some regimes. I show that these missing photons can be counted exactly for some numerical schemes. Adding them back into a simulation, however, is a much harder task. I show one example of such a scheme, which achieves sub-percent accuracy on simple tests. Unfortunately, the scheme performs much worse on semi-realistic simulations of cosmic reionization, leading to a faster overlap and significant errors in the point-to-point comparison of the RSL radiation field with the reference simulation that maintains the full speed of light for the radiative transfer.

James Webb Space Telescope (JWST) has uncovered a new population of compact objects that show a unique V-shaped spectral energy distribution (SED) in the UV and optical wavelength range. These so-called "little red dots" (LRDs) often exhibit broad Balmer emission lines, indicative of the presence of active galactic nuclei (AGNs). They generally lack detection of X-ray, radio, and mid-IR radiation, which is fundamentally different from typical AGNs. Various models, including super-Eddington-accreting black holes enshrouded in dense and dust-poor gas, have been proposed to explain these features. However, the nature of LRDs remains debated, and their evolutionary fate is unclear. Here we report two unusual LRDs at redshift z = 2.868 and 2.925 that are in a transitional phase towards typical AGNs. Their V-shaped SEDs, compact optical morphology, and broad emission lines satisfy the defining criteria of LRDs. On the other hand, they exhibit intense X-ray, radio, and mid-IR radiation that is much stronger than previously known LRDs. These hybrid properties suggest that the dense gas envelope around their central black holes is dispersing, allowing high-energy photons and radio emission to escape. Meanwhile, a dust torus is forming. This finding provides a direct insight into the nature of LRDs and indicates that at least some LRDs will evolve into AGNs or quasars at later times.

Identifying the sources of high-energy (TeV-PeV) astrophysical neutrinos is crucial for studies in both astrophysics and particle physics. Despite extensive searches for more than a decade, which revealed several individual potential sources and only one potential source class, the origins of these neutrinos remain largely unresolved; thus, more source classes should be investigated. In this work, we conduct the first search for high-energy neutrino emission from a new source class, supermassive black hole binaries (SMBHBs), which are also theoretically motivated. We perform an unbinned maximum-likelihood-ratio analysis on our constructed catalog of 693 SMBHB candidates and 10 years of IceCube public data. Our results show positive correlations, with higher significance in more physically motivated scenarios and the highest significance at 3.0$\sigma$. In addition, we also study potential connections between SMBHBs' high-energy neutrino and nano-Hz gravitational-wave emissions, the latter being the main target of pulsar timing arrays. Our results provide the first evidence of SMBHBs being high-energy neutrino emitters.

Interstellar objects (ISOs), particularly those with cometary activity, provide unique insight into the primordial physical and chemical conditions present during the formation of the planetary system in which they originated. Observations in the sub-mm regime allow for direct measurements of several parent molecules released from the comet nucleus into the coma. Here we present observations of the third ISO, 3I/ATLAS, with the `Ū`ū heterodyne receiver on the James Clerk Maxwell Telescope (JCMT), which targeted emission from HCN($J = 3 - 2$) and CO($J = 2 - 1$). Our observations, taken between 16 July 2025 and 21 July 2025 (UT), when 3I/ATLAS was at a heliocentric distance between 4.01 and 3.84 au, provide the earliest sub-mm constraints on its activity. We do not detect HCN or CO in these epochs, with 3$\sigma$ upper-limits on the production rates of $Q(HCN) < 1.7 \times 10^{24}$ s$^{-1}$ at $r_h = 4.01 - 3.97$ au and $Q(CO) < 1.1 \times 10^{27}$ s$^{-1}$ at $r_h = 3.94 - 3.84$ au, respectively. We combine this HCN limit with later JCMT observations of HCN to constrain its temporal evolution. Fitting the HCN detections with a $Q(HCN) \propto r_h^{-n}$ model and accounting for the upper-limits yields $n = 12.7^{+6.9}_{-2.5}$. This slope is steeper than those of typical Solar System comets, but consistent with the production rate slopes measured for other species in the coma of 3I/ATLAS.

The Milankovitch cycles of Earth result from gravitational interactions with other bodies in the Solar System. These interactions lead to slow changes in the orbit and angular momentum vector of Earth, and correspondingly influence Earth's climate evolution. Several studies have shown that Mars may play a significant role in these Milankovitch cycles, such as the 2.4 Myr eccentricity cycle related to perihelion precession dynamics. Here we provide the results of a detailed dynamical analysis that explores the Earth Milankovitch cycles as a function of the Martian mass to quantify the extent that Mars influences variations in Earth's orbital eccentricity, the longitude of perihelion, the longitude of the ascending node, and obliquity (axial tilt). Our results show that, although the 405 kyr long-eccentricity metronome driven by $g_2$ (Venus) and $g_5$ (Jupiter) persists at all Mars masses, the $\sim$100 kyr short-eccentricity bands driven by $g_4$ (Mars) lengthen and gain power as Mars becomes more massive, consistent with enhanced coupling among inner-planet $g$-modes. The 2.4 Myr grand cycle is absent when Mars approaches zero mass, reflecting the movement of $g_4$ with the Martian mass. Meanwhile, Earth's obliquity cycles driven by $s_3$ (Earth) and $s_4$ (Mars) lengthen from the canonical $\sim$41 kyr with increasing Mars mass, relocating to a dominant 45--55 kyr band when the mass of Mars is an order of magnitude larger than its present value. These results establish how Mars' mass controls the architecture of Earth's climate-forcing spectrum and that the Milankovitch spectrum of an Earth-like planet is a sensitive, interpretable probe of its planetary neighborhood.

We provide a new measurement of the expansion history of the Universe at $z=0.54$ with the cosmic chronometers (CC) method, exploiting the high-quality spectroscopic VLT/MUSE data for three galaxy clusters in close-by redshift bins: SDSS J2222+2745 ($z=0.49$), MACS J1149.5+2223 ($z=0.54$), and SDSS J1029+2623 ($z=0.59$). The central one, MACS J1149.5+2223, hosts the well-known supernova 'Refsdal', which allowed for $H_0$ measurements via time delay cosmography (TDC). This represents the first step for a self-consistent probe combination, where different methods are applied to the same data sample. After selecting the most passive and massive cluster members (38 CCs), we derive their age and physical parameters via full spectrum fitting. We use the code Bagpipes, specifically modified to remove the cosmological prior on ages. On average, the CC sample shows super-solar metallicities $Z/Z_{\odot} = 1.3 \pm 0.7$, low dust extinction $A_{\rm{V}} = 0.3 \pm 0.3$ mag and to have formed in short bursts $\tau = 0.6 \pm 0.2$ Gyr. We also observe both an ageing trend in redshift and a mass-downsizing pattern. From the age-redshift trend, implementing the CC method through a bootstrap approach, we derive a new $H(z)$ measurement: $H$($z$=0.542) = $66_{-29}^{+81}$ (stat) $\pm$13 (syst) km/s/Mpc. We also simulate the impact of increased statistics and extended redshift coverage, finding that $H$($z$) uncertainties can be reduced by up to a factor of 4 with $\sim$100 CCs and a slightly broader redshift range (d$z\sim$0.2).

We explore the impact of "slow-light" radiative transfer - i.e., general relativistic radiative transfer (GRRT) calculations in which the simulated fluid evolves while light rays are propagating through it - in general relativistic magnetohydrodynamic (GRMHD) models of the M87 jet. Because the plasma in the jet-launching region is accelerated to relativistic velocities, and because the jet in M87 is nearly aligned with the line of sight (offset by ~17 degrees), a slow-light treatment is important for accurately modeling the observable structure. While fast-light images exhibit prominent helical or loop-shaped features in the jet - which we associate with narrow bundles of magnetic field lines - these features become stretched and smoothed-out in slow-light images. Our slow-light images instead exhibit a double-edged, cone-like morphology that is more consistent with observations of M87 than corresponding fast-light images. We find that the radius at which the plasma transitions from sub-relativistic to relativistic velocities is imprinted on slow-light images via a transition from loop-dominated at small distances from the black hole to edge-dominated at a larger distance, with the loop-edge transition occurring at larger distances for lower black hole spins. The jet image dynamics also vary with black hole spin, with low-spin models producing jets that exhibit substantial "wobbling", while high-spin models produce jets that are straighter and more stable in time. The spin-dependent jet morphology and variability are revealed by slow-light imaging because slow-light effects become more enhanced as the plasma velocity becomes more relativistic, and because the plasma acceleration is itself a strong function of the spin.

We present ClumPyLen, a Python-based simulator designed to produce realistic mock observations of strongly lensed, high-redshift, clumpy star-forming galaxies. The tool models galaxy components such as disks, bulges, and spiral arms using Sérsic profiles, and it populates them with stellar clumps whose properties are sampled from physically motivated distributions. ClumPyLen includes the effects of gravitational lensing through user-provided deflection angle maps and simulates realistic observational conditions by accounting for instrumental effects, Point-Spread-Function convolution, sky background, and photon noise. The simulator can support a wide range of filters and instruments; here we focus on HST/ACS, HST/WFC3-IR, and JWST/NIRCam. We demonstrate the capabilities of the code through two examples, including a detailed simulation of the z = 6.145 source Cosmic Archipelago lensed by MACS J0416.1-2403. The simulated images closely match the morphology and limiting magnitudes of real observations. ClumPyLen is designed to explore the detectability of stellar clumps in terms of mass and size, especially in the low-mass regime, and it allows the study of clump blending effects. Thanks to its modular design, the code is highly adaptable to a wide range of scientific goals, including lensing studies, galaxy evolution, and the generation of synthetic datasets for machine learning or forward modeling applications.

Selection effects, such as interstellar extinction and varying survey depth, complicate efforts to determine the gravitational potential - and thus the distribution of baryonic and dark matter - throughout the Milky Way galaxy using stellar kinematics. We present a new variant of the "Deep Potential" method of determining the gravitational potential from a snapshot of stellar positions and velocities that does not require any modeling of spatial selection functions. Instead of modeling the full six-dimensional phase-space distribution function $f\left(\vec{x},\vec{v}\right)$ of observed kinematic tracers, we model the conditional velocity distribution $p\left(\vec{v}\mid\vec{x}\right)$, which is unaffected by a purely spatial selection function. We simultaneously learn the gravitational potential $\Phi\left(\vec{x}\right)$ and the underlying spatial density of the entire tracer population $n\left(\vec{x}\right)$ - including unobserved stars - using the collisionless Boltzmann equation under the stationarity assumption. The advantage of this method is that unlike the spatial selection function, all of the quantities we model, $p\left(\vec{v}\mid\vec{x}\right)$, $\Phi\left(\vec{x}\right)$, and $n\left(\vec{x}\right)$, typically vary smoothly in both position and velocity. We demonstrate that this "conditional" Deep Potential method is able to accurately recover the gravitational potential in a mock dataset with a complex three-dimensional dust distribution that imprints fine angular structure on the selection function. Because we do not need to model the spatial selection function, our new method can effectively scale to large, complex datasets while using relatively few parameters, and is thus well suited to Gaia data.

Kappa Velorum (Markeb, HD 81188) is one of the brightest stars in the Southern sky and has long been known to be a single-lined spectroscopic binary. The binary mass function is large, $f(M)=1.15\ M_\odot$, suggesting that the bright (V=2.5) B2IV star may host a dark, compact object companion. We use VLTI GRAVITY observations to definitively test this possibility by directly resolving the binary. We detect a main sequence B star companion and rule out the compact object scenario. By combining the relative astrometric orbit and archival radial velocities, we report an updated precise characterization of the orbit (period $P=116.795\pm0.002$ d, eccentricity $e=0.1764\pm0.0004$, inclination $i=74.04\pm0.01^{\circ}$) and estimate the masses of the B stars. Using the original Hipparcos parallax measurement $\varpi = 6.05\pm0.48$ mas, we find $M_1 = 10^{+4}_{-2}\ M_\odot$ and $M_2 = 6.9\pm1.0\ M_\odot$. The uncertainties on the masses are primarily driven by the uncertain parallax, which we find is likely biased by the orbital motion. We use an archival UVES spectrum and MIST evolutionary tracks to refine our mass estimates. Finally, we discuss how interferometry and high-contrast imaging may be used to characterize other candidate star+compact object binaries, including those that will be discovered with Gaia DR4, as part of a larger effort to uncover the hidden population of black holes in the Milky Way.

We report on four Local Volume dwarf galaxies identified through our ongoing SEmi-Automated Machine LEarning Search for Semi-resolved galaxies (SEAMLESS): Hydrus A, LEDA 486718, Cetus B, and Sculptor 26, with the discovery of Hydrus A reported here for the first time. These four galaxies span a wide range of environments and evolutionary states. Hydrus A (MV = -9.39+/-0.20, D = 3.38-0.30+0.32 Mpc) and LEDA 486718 (MV = -11.62+/-0.08, D = 4.80+/-0.17 Mpc) are among the most isolated dwarfs known within 5 Mpc, while Cetus B (MV = -8.26+/-0.17, D = 3.32-0.23+0.25 Mpc) and Sculptor 26 (MV = -11.25+/-0.10, D =3.21+/-0.13 Mpc) lie < 2 Rvir of NGC 253. Hydrus A shows properties consistent with quenching driven by cosmic reionization, cosmic-web interactions, or internal feedback. LEDA 486718 is an isolated star forming dwarf. Cetus B appears quenched and morphologically disturbed, making it a low-mass satellite or backsplash candidate, while Sculptor 26 is red and seemingly gas-poor but displays signs of recent activity, consistent with a transitional evolutionary state. Together, these systems demonstrate the power of SEAMLESS for building a census of faint galaxies beyond the Local Group.

We present the abundances of 23 elements, including 11 heavy elements (Cu, Zn, Sr, Y, Zr, Ba, La, Ce, Nd, Sm, Eu) for up to 86 metal-poor (-2 < [Fe/H] < -1) subgiants. We use KORG, a state of the art spectral synthesis package, to derive 1D-LTE abundances from high-SNR and high-resolution spectra taken by the Large Binocular Telescope with the Potsdam Echelle Polarimetric and Spectroscopic Instrument. These precise spectra and abundance measurements minimize the impact of photon-noise (<0.06 dex), allowing us to robustly measure the intrinsic abundance scatter in [X/Fe]. After removing two stars with exceptional s-process enhancement, we find that the intrinsic scatter among the s- and r-process elements tends to be larger than for the lighter elements, with heavy element scatter ranging from 0.11 (Zn) to 0.27 (Eu) dex. Intrinsic abundance scatter could have multiple origins, including star-to-star variations in the ratios of nucleosynthetic sources as well as stochastic sampling of the progenitor supernovae properties, such as mass, rotation, and magnetic field strength. We explore the expected abundance scatter signature caused by stochastic sampling, finding that a fraction of both rapidly rotating CCSN and magnetorotationally driven SN are needed to reach the observed abundances and intrinsic scatter. This analysis is limited by the restrictive parameter spaces spanned by existing yield sets. A diverse, finely sampled grid of supernovae yields is needed to robustly model stochastic abundance scatter.

The most recently formed young stellar objects (YSOs) in active star forming regions are excellent tracers of their parent cloud motion. Their positions and dynamics provide insight into cluster formation and constrain kinematic decoupling timescales between stars and gas. However, because of their strong extinction and young age, embedded YSOs are mainly visible at infrared wavelengths and thus absent from astrometric surveys such as Gaia. We measured the proper motions of 6,769 sources toward the NGC 2024 cluster in the Flame Nebula using multi-epoch near-infrared observations from three ESO public surveys: VISIONS, VHS, and the VISTA/VIRCAM science verification program. Cross-validation of our results with Gaia using optically visible stars shows excellent agreement, with uncertainties on the same order of magnitude. For 362 YSO candidates identified from the literature, we derived proper motions on the order of <5 mas/yr with mean measurement uncertainties of ~0.22 mas/yr. This is the first homogeneous proper motion measurement of this quality for more than half of these stars. For Class I and flat-spectrum sources, our results provide a >13-fold increase in available proper motion measurements. We analyzed the positional and kinematic differences between YSO classes and confirmed a previously reported inside-out age segregation from younger to older stars, likely driven by an outward movement of older stars. No evidence of prolonged hierarchical assembly was found. Instead, the results support a rapid (<1 Myr) cluster collapse. This scenario also accounts for the observed slightly higher 1D velocity dispersion of Class I sources relative to Class flat objects. YSO radial velocities generally align with the gas velocities measured from 12CO(3-2), HNC(1-0), HCN(1-0), and show a weaker correlation with N_2H+(1-0). Some Class II and III objects appear to be already decoupling.

The rapid neutron-capture ($r$-process) element europium (Eu) is a valuable tracer of neutron star mergers and other rare nucleosynthetic events. The stellar spectroscopic survey GALAH's unique wavelength range and setup include the Eu absorption feature at $\sim6645$ Å for almost a million stars in the most recent Data Release 4 (DR4). However, DR4 also saw a decreased precision in reported Eu measurements compared to previous data releases. In this work, we use a convolutional neural network (CNN) to perform label transfer, wherein we use the GALAH DR4 spectra and stellar parameters to infer DR3 [Eu/H] abundances. This CNN is then applied to DR4 spectra without corresponding DR3 Eu abundances to develop a new, publicly available catalogue of [Eu/H] values for high signal-to-noise targets. We include [Eu/H] predictions for $118\,946$ stars, out of which $54\,068$ giants constitute our "golden sample" of high-confidence predictions, which pass stricter quality cuts and have a reported precision $\lesssim0.1$. To overcome the scarcity of training data in the low metallicity regime, we provide an additional catalogue of [Eu/H] abundances for metal poor ($\mathrm{[Fe/H]}<-1$) stars derived from synthesis of the Eu feature. Our "golden sample" can be combined with [Eu/H] values from GALAH DR3 to create a catalogue of over $100\,000$ vetted, high-quality abundances on a homogeneous scale. Moreover, we are able to reproduce known science results, including the elevated Eu abundances of accreted stars and previously observed Galactic chemical evolution trends. This catalogue represents one of the largest available samples of [Eu/H] abundances for high signal-to-noise targets.

We present the OASIS survey program discovery of a substellar companion orbiting the young A1V star HIP 71618, detected using precision astrometry from \textit{Gaia and Hipparcos} and high-contrast imaging with SCExAO/CHARIS and Keck/NIRC2. Atmospheric modeling favors a spectral type of M5--M8 and a temperature of $\sim$2700 $\pm$ 100 $K$. Dynamical modeling constrains HIP 71618 B's mass to be ${60}_{-21}^{+27}$ $M_{\rm Jup}$ or ${65}_{-29}^{+54}$ $M_{\rm Jup}$, depending on the adopted companion mass prior. It has a nearly edge-on, 11 au-orbit with a high eccentricity. HIP 71618 B will be located within Roman Coronagraph's dark hole region during the instrument's technological demonstration phase. A high signal-to-noise ratio detection of HIP 71618 B at 575 nm would demonstrate a 5-$\sigma$ contrast of 10$^{-7}$ or better. The system is also located within or very close to Roman's Continuous Viewing Zone -- near multiple candidate reference stars for dark-hole digging -- and its primary is bright ($V$ $\approx$ 5). The suitability of HIP 71618 as one potential Roman Coronagraph target for demonstrating the instrument's core requirement (TTR5) should motivate the timely, deep vetting of candidate reference stars.

We use the FLAMINGO galaxy formation model to quantify the impact of baryon-CDM isocurvature perturbations on galaxy formation in $\Lambda$CDM. In linear theory, these perturbations represent local, compensated variations in the ratio between the baryon and CDM densities; they freeze in amplitude at late times, with an RMS amplitude of $1.5\%$ on the Lagrangian scale of a $10^{11}\,\rm M_\odot$ halo ($0.85\, \rm{Mpc}$). Although such perturbations arise naturally within $\Lambda$CDM, most cosmological simulations and semi-analytic models to date omit them. These perturbations are strongly anti-correlated with the matter overdensity field such that halos form with baryon fractions below the cosmic mean, with earlier-collapsing halos exhibiting stronger baryonic suppression. To isolate the galaxy response, we analyse three hydrodynamical simulations with identical initial matter overdensity fields that: i) include isocurvature modes, ii) omit them, or iii) invert their amplitude. At $z=8$, isocurvature perturbations reduce the mean baryon fraction and star formation rates of resolved halos by $ 5\%$ and $ 12\%$, respectively, relative to the null-isocurvature case. These effects are almost independent of halo mass and diminish steadily with time, reaching $ 0.1\%$ and $ 1\%$ by $z=0$. We develop a model based on spherical collapse that accurately reproduces the mean baryon fraction suppression. As high-redshift observations become increasingly routine, incorporating isocurvature perturbations into simulations and semi-analytic models will be important for robust predictions of early galaxy and black hole formation in the JWST era.

We present global 3D GRMHD simulations of black hole (BH) accretion disks designed to investigate how MRI-driven dynamo action regulates jet formation and evolution. Unlike standard SANE/MAD setups that impose a coherent large-scale poloidal loop, our "sub-SANE" initial conditions use multiple same-polarity small-scale magnetic loops. Rapid reconnection erases magnetic memory and enables large-scale dynamo to emerge early from MRI turbulence. We perform two such sub-SANE simulations at different BH spins ($a = 0.5, 0.9375$) and compare them with conventional SANE runs. The sub-SANE disks show regular large-scale dynamo cycles with periods of about ten orbits. Decomposition of the induction equation shows that the turbulent dynamo term is stronger in 3D compared to 2.5D and balances advection in the saturated state, confirming sustained large-scale field generation. These dynamo-generated fields are advected inward with minimal time lag, producing correlated peaks in both poloidal and toroidal field strengths from $r_{\rm max}$ to the horizon. Early in the evolution, these peaks imprint directly onto the jet's electromagnetic energy flux, indicating that the jet mirrors the dynamo wave. Though jets form at early times, the sub-SANE runs eventually undergo jet shutdown. We show that this occurs when the magnetic field at the horizon loses coherence, as quantified by a decline in the signed-to-unsigned flux ratio $\mathcal{C}_{\rm BH}$ below $\approx 0.6$. In contrast, the SANE reference case with similar accretion rate and horizon magnetic flux maintains high magnetic coherence because its initial large-scale field persists, allowing its jet to survive. Our results show that both dynamo-driven field evolution and horizon magnetic-field coherence critically regulate jet longevity, establishing a direct dynamo-jet connection in GRMHD disks.

Von Zeipel-Lidov-Kozai (ZLK) oscillations, induced by bound, perturbative companions to white dwarfs, have been suggested as a dynamical mechanism that may contribute to white dwarf pollution. To trigger ZLK oscillations, however, a 3-body system must reach a sufficiently large mutual inclination between orbits. The occurrence of these high-mutual-inclination configurations can be curtailed by dissipative precession at the protoplanetary disk stage, which pushes exoplanet-hosting close binary systems toward preferential orbit-orbit alignment. In this work, we constrain the fraction of white dwarfs with binary companions that can undergo ZLK-driven pollution given the effects of dissipative precession. To accrete pollution via ZLK oscillations, a white dwarf binary system must be sufficiently inclined and the characteristic timescale of the oscillations must be sufficiently short to perturb material within the white dwarf's cooling age. Considering a sample of 4400 known white dwarf/main sequence binaries, we find that $50-70\%$ have favorable parameters for ZLK pollution, depending on the orbital separation of the polluting body. While the conditions for oscillations are favorable, the tendency for ZLK to result in massive but more infrequent polluters likely restricts the rates of ZLK-induced pollution among the observed population. In general, dissipative precession is a limiting factor in pollution rates for more closely separated binaries (initial separations $<500-800$~au), while ZLK timescale constraints are most limiting for wider binaries.

Tidal Disruption Events (TDEs) have long been hypothesized as valuable indicators of black holes, offering insight into their demographics and behaviour out to high redshift. TDEs have also enabled the discovery of a few Massive Black Holes (MBHs) with inferred masses of $10^4$--$10^6\,M_\odot$, often associated with the nuclei of dwarf galaxies or ultra-compact dwarf galaxies (UCDs). Here we present AT2022zod, an extreme and short-lived optical flare in an elliptical galaxy at z=0.11, residing within 3kpc of the galaxy's centre. Its luminosity and ~30-day duration make it unlikely to have originated from the host galaxy's central supermassive black hole (SMBH), which estimate to have a mass around $10^8\,M_\odot$. Assuming that the emission mechanism is consistent with known observed TDEs, we find that such a rapidly evolving transient could either be produced by a MBH in the intermediate-mass range or, alternatively, result from the tidal disruption of a star on a non-parabolic orbit around the central SMBH. We suggest that the most plausible origin for AT2022zod is the tidal disruption of a star by a MBH embedded in a UCD. As the Vera C. Rubin Observatory's Legacy Survey of Space and Time comes online, we propose that AT2022zod serves as an important event for refining search strategies and characterization techniques for intermediate-mass black holes.

EP240222a is the first intermediate-mass black hole (IMBH) tidal disruption event (TDE) captured in real-time with multi-wavelength observations and spectroscopic confirmation. However, its light curves deviate substantially from previous theoretical expectations. Motivated by these unique features, we have developed a novel model that successfully reproduces its peculiar evolution. Our model delineates five stages: (1) Initial Stage of inefficient circularization; (2) Slow-Rising Stage with a faint X-ray precursor disk fed by successive self-crossings; (3) Fast-Rising Stage, where delayed stream-disk interaction at momentum flux matching drives a sharp luminosity rise; (4) Plateau Stage with super-Eddington accretion, outflow, reprocessing, and a clear polar line-of-sight; and (5) Decline Stage of sub-Eddington accretion and ongoing reprocessing. Our fit indicates the disruption of a $M_* \approx 0.4~M_\odot$ main-sequence (MS) star with a penetration factor $\beta \approx 1.0$. Our model, which incorporates key TDE processes, establishes EP240222a-like light curves as typical IMBH-TDE signatures. The distinctive identifier is a slow rise in X-rays and a corresponding slow rise/quasi-plateau in the UV/optical, followed by a brighter, super-Eddington plateau in both bands, though other forms exist, such as the rapid rise from white dwarf (WD) disruptions over minutes to days.

We present the discovery of a superjovian planet around the young A5 star HIP 54515, detected using precision astrometry from the \textit{Hipparcos Gaia Catalogue of Accelerations} and high-contrast imaging with SCExAO/CHARIS from the recently-commenced \surveyname\ program. SCExAO/CHARIS detects HIP 54515 b in five epochs 0\farcs{}145--0\farcs{}192 from the star ($\sim$3--4 $\lambda$/D at 1.65 $\mu m$), exhibiting clockwise orbital motion. HIP 54515 b lies near the M/L transition with a luminosity of log(L/L$_{\rm \odot}$) $\sim$ -3.52 $\pm$ 0.03. Dynamical modeling constrains its mass and mass ratio to be ${17.7}_{-4.9}^{+7.6}$ $M_{\rm Jup}$ and ${0.0090}_{-0.0024}^{+0.0036}$ and favors a $\sim$25 au semimajor axis. HIP 54515 b adds to a growing list of superjovian planets with moderate eccentricities (e $\approx$ 0.4). Now the third planet discovered from surveys combining high-contrast extreme adaptive optics imaging with precision astrometry, HIP 54515 b should help improve empirical constraints on the luminosity evolution and eccentricity distribution of the most massive planets. It may also provide a key technical test of the Roman Space Telescope Coronagraph Instrument's performance in the low stellar flux, small angular separation limit and a demonstration of its ability to yield constrainable planet spectral properties.

To investigate the impact of nuclear reaction rate uncertainties on type I X-ray burst nucleosynthesis, comprehensive Monte Carlo simulations are performed with temperature-independent and -dependent variations in reaction rates using the REACLIB and STARLIB libraries, respectively. A total of 1,711 $(p, \gamma)$, $(p, \alpha)$, $(\alpha, p)$, and $(\alpha, \gamma)$ reactions are varied simultaneously, along with their inverse reactions, via detailed balance. For the first time, it is found that Monte Carlo sampling with larger perturbations to these reaction rates may lead to multi-peaked abundance distributions for some isotopes. These multi-peak structures arise not only from coupled reactions but also, in some cases, from single reactions. Our study also confirmed previously identified key reactions and provides more robust lists. These reactions deserve priority consideration in future study.

High-energy astrophysical sources such as active galactic nuclei, quasars, X-ray binaries, and gamma-ray bursts are powered by mechanisms that convert gravitational or rotational energy into radiation, jets, and relativistic outflows. Understanding the physics of these processes remains a major challenge. Black holes have traditionally served as the central engines behind such phenomena, with well established energy extraction mechanisms including the Penrose process, the Blandford-Znajek process, and the Banados-Silk-West mechanism. However, studies in general relativity indicate that, under certain conditions, gravitational collapse may lead to the formation of naked singularities or other horizonless compact objects, which could in principle allow more efficient energy extraction than classical black holes. This brief review summarizes recent progress on energy extraction mechanisms in naked singularity spacetimes. We examine the roles of rotation, electromagnetic fields, and particle interactions in shaping extraction efficiency and dynamics. Particular attention is given to negative energy orbits and ergoregion physics, which enable Penrose type and magnetic Penrose mechanisms without an event horizon. We also discuss collisional Penrose processes and particle acceleration near the singularity, emphasizing their potential astrophysical implications. By comparing extraction efficiencies and physical conditions in black holes and naked singularities, we highlight how the absence of a horizon fundamentally alters the dynamics of energy release. These results suggest that naked singularities may serve as natural laboratories for strong field gravity and as alternative engines for high-energy astrophysical phenomena in the era of multi-messenger observations.

The origin of the Sun's slow wind and its inherent variability remain unknown, but there is increasing evidence that interaction between closed and open magnetic flux in the corona plays a role. This paper studies the dynamic evolution of streamers and pseudostreamers with a particular focus on the release of plasma from the closed to the open field region. We employ a global magnetohydrodynamic model that extends from the solar chromosphere to 30 solar radii, and that extends previous interchange magnetic reconnection modelling by including self-consistent thermodynamics. We find that at both the helmet streamer and pseudostreamer there is a continual dynamic interaction between closed and open flux. At the helmet streamer, this takes the form of a ''breathing'' cycle in which the closed flux contracts and expands, and plasmoids are released along the heliospheric current sheet. The pseudostreamer exhibits a back-and-forth motion, driving interchange reconnection alternately on its opposite flanks. The resulting release of hot, dense plasma leads to density fluctuations in the open field that are significantly larger above the helmet streamer due to the persistence of the plasmoids there. Our model demonstrates that plasma is continually being released into the heliosphere from both streamers and pseudostreamers.

We present results of a spectroscopic campaign around the diffuse dwarf galaxy Crater 2 (Cra2) and its tidal tails as part of the Southern Stellar Stream Spectroscopic Survey ($S^5$). Cra2 is a Milky Way dwarf spheroidal satellite with extremely cold kinematics, but a huge size similar to the Small Magellanic Cloud, which may be difficult to explain within collisionless cold dark matter. We identify 143 Cra2 members, of which 114 belong to the galaxy's main body and 29 are deemed part of its stellar stream. We confirm that Cra2 is dynamically cold (central velocity dispersion $2.51^{+0.33}_{-0.30}\,{\rm km\,s^{-1}}$) and also discover a $\approx$7$\sigma$ velocity gradient consistent with its tidal debris track. We separately estimate the stream velocity dispersion to be $5.74^{+0.98}_{-0.83}\,{\rm km\,s^{-1}}$. We develop a suite of $N$-body simulations with both cuspy and cored density profiles on a realistic Cra2 orbit to compare with $S^5$ observations. We find that the velocity dispersion ratio between Cra2 stream and galaxy ($2.30^{+0.41}_{-0.35}$) is difficult to reconcile with a cuspy halo with fiducial concentration and an initial mass predicted by standard stellar mass$-$halo mass relationships. Instead, either a cored halo with relatively small core radius or a low-concentration cuspy model can reproduce this ratio. Despite tidal mass loss, Cra2 is metal-poor ($\langle \rm[Fe/H]\rangle=-2.16\pm0.04$) compared to the stellar mass$-$metallicity relation for its luminosity. Other diffuse dwarf galaxies similar to Cra2 in the Local Group (Antlia 2 and Andromeda 19) also challenge galaxy formation models. Finally, we discuss possible formation scenarios for Cra2, including ram-pressure stripping of a gas-rich progenitor combined with tides.

Stellar activity and planetary magnetospheres are powered by an underlying dynamo mechanism generated by magnetoconvection coupled with rotation. In astrophysical contexts, magnetoconvection typically occurs in parameter regimes that are currently inaccessible to direct numerical simulations (DNS). We investigate convective heat transfer in a magneto-convection and dynamo model under extreme parameter conditions, specifically high Rayleigh and Prandtl numbers, in a plasma flow with maximum kinetic helicity compatible with fast-rotating objects. Our approach to studying magneto-convection and dynamo mechanisms employs a simplified thermally driven shell model. Magnetic polarity reversals are obtained by including a pitchfork bifurcation term in the large-scale magnetic field equation, while nonlinear dynamics are described by a shell model formulation. The low computational cost of the model allows us to explore the asymptotic behavior of convective heat transfer in regimes beyond those reached by current DNS. Our results reveal that the Nusselt number $Nu$ -- a dimensionless measure of convective heat transport-- generally increases with turbulence, following a power-law scaling and showing a strong correlation with Ra and Pr. This relationship appears to be more pronounced than that observed in non-magnetized fluids, suggesting that magnetic fields may significantly enhance convective heat transfer. Despite the assumption to neglect spatial information such as density stratification -- an assumption that is necessary in the shell model approach -- our model captures the gross dynamical features of turbulent magnetoconvection in asymptotic regimes. It allows for a broad exploration of parameter space, indicating that magnetic fields may play a central role in modulating heat transport in stellar and planetary interiors.

This work presents a study of the evolution of the stellar accretion rates of pre-main-sequence intermediate-mass stars. We compare the accretion rate of the younger intermediate-mass T Tauri stars (IMTTSs) with the older Herbig stars into which they evolve. We find that the median accretion rate of IMTTSs (1.2$\times$10$^{-8}$ M$_{\odot}$ yr$^{-1}$) is significantly lower than that of Herbig stars (1.9$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$). This increase stands in stark contrast with canonical models of disk evolution that predict that the stellar accretion rate declines with age. We put forward a physically plausible scenario that accounts for the systematic increase of stellar accretion based on the increase of the effective temperature of the stars as they evolve towards the zero-age main sequence. For example, the temperature of a 2M$_{\odot}$ star will increase from 4900~K in the IMTTS phase to 9100~K during the Herbig phase. Thus, the luminosity of the far ultraviolet (FUV) radiation will increase by orders of magnitude. We propose that this increase drives a higher stellar accretion rate. The scenario we propose to account for the increase in the stellar accretion rate solves the lifetime problem for Herbig disks because the increasing stellar accretion rates require lower initial disk masses to account for present-day disk masses. This work highlights the importance of the role FUV radiation has in driving the accretion rate, predicts a large population of pre-main-sequence non-accreting A stars, and has implications for interpreting disk morphologies that may serve as signposts of embedded gas giant planets in Herbig disks.

We use the cyclic figure of merit to determine the likelihood of achieving cyclic deconvolution for 312 pulsars with sub-40 ms spin periods across 14 different telescope-observing frequency combinations. We find that the optimal frequency range for achieving cyclic deconvolution for most pulsars is between $\sim$80$-$300 MHz, making low frequency observatories like uGMRT, LOFAR, and MWA the best-suited instruments for the technique. Moreover, we find that, as quantified by the total number sources with sufficient cyclic merits that are observed within the full deconvolution regime, uGMRT is likely the best current instrument for cyclic spectroscopy among the nine telescopes we considered, with LOFAR being the second best, and future telescopes like SKA-Low serving as the most likely eventual successors. The relatively high cyclic merit of the Crab Pulsar in the frequency ranges considered for GBT, MWA, LOFAR, and uGMRT suggests that some faster-spinning canonical pulsars may be able to achieve cyclic deconvolution, and we discuss potential follow-up analyses on other non-recycled pulsars. We conclude by advocating for near real-time cyclic spectroscopy backends to be considered for current and upcoming low frequency telescopes to increase the accessibility of this technique.

M-dwarfs are the most dominant stars in the Galaxy. Their interiors and atmospheres exhibit complex processes including dust condensation, convective feedback, and magnetic activity-driven heterogeneity. Standard stellar characterization methods often struggle to capture these coupled effects. Part I of this series introduced SPHINX I, a validated grid of self-consistent radiative-convective model atmospheres and spectra for M-dwarfs with up-to-date molecular opacities suitable for early-to-mid M-dwarfs. Here, we present SPHINX II, which extends the model grid to cover mid-to-late type M-dwarfs, including both gray and physically motivated condensate cloud treatments and shorter convective mixing lengths. We validate SPHINX II using 39 benchmark FGK+M binary systems observed with SpeX IRTF (Mann et al. 2014) and apply it to 32 mid-to-late-type M-dwarfs from the SpeX Prism Library. SPHINX II yields improved fits that are statistically consistent with empirical benchmarks, achieving precisions of 0.078 dex in metallicity and 0.13 dex in C/O. Across the model grid, condensate cloud mass peaks between 2100-2400 K, decreasing sharply toward both cooler and hotter temperatures. We find the onset of the cloud-free regime around 2900 K, and below 2100 K, we see formation of deep/buried clouds. As a case study, we also model Trappist-1 and show that even mass-limited silicate grains subtly modify its emergent spectrum, suppressing near-infrared flux and reddening the mid-infrared slope via shallow cloud formation near 1e-2 bar. In sum, SPHINX II provides an improved framework for constraining the fundamental properties of mid-to-late M-dwarfs.

We present a three-dimensional numerical model of a solitary spherically polarized Alfvén wave packet -- an Alfvénon, characterized by open field-line topology and magnetic field reversals, resembling the switchbacks observed by Parker Solar Probe to be a nearly ubiquitous feature of turbulence in the inner heliosphere. Direct magnetohydrodynamic simulations of the constructed Alfvénon demonstrates remarkable stability, confirming its nature as an exact, nonlinear solution of the ideal MHD equations.

Active region recurrent jets are manifestations of episodic magnetic energy release processes driven by complex interactions in the lower solar atmosphere. While magnetic flux emergence and cancellation are widely recognized as key contributors to jet formation, the mechanisms behind repeated magnetic reconnection remain poorly understood. In this letter, we report a sequence of nine recurrent jets originating from active region AR 12715 during its decay phase, where the jet activity was associated with a complex distribution of fragmented magnetic flux. Non-linear force-free field (NLFFF) extrapolations reveal the presence of low-lying, current-carrying loops beneath overarching open magnetic fields near the jet footpoints. These magnetic structures were perturbed by (i) emerging flux elements and (ii) interactions between oppositely polarized moving magnetic features (MMFs). To interpret these observations, we compare them with 3D radiative MHD simulation from the Bifrost model, which reproduce jet formation driven by interacting bipolar MMFs, leading to subsequent flux cancellation in the photosphere. Our results emphasize the critical role of MMF-driven flux interactions in initiating and sustaining recurrent jet activity in active regions.

The formation channels of binary neutron stars (BNSs) remain uncertain. The detection of GW190425 by LIGO/Virgo/KAGRA (LVK) suggests a subpopulation of massive BNSs, possibly formed through unstable "case BB" mass transfer with short merger delays. We investigate whether next-generation detectors such as the Einstein Telescope (ET) can identify and characterise such subpopulations. Using the latest LVK constraints on the BNS merger rate, we generate mock ET catalogues containing a mixture of light and heavy subpopulations. The redshift distribution of each subpopulation is modeled as the convolution of the cosmic star formation rate with a time-delay distribution: heavy BNSs have fixed short delays, while light BNSs follow power-law delays with indices -0.5,-1,-1.5. Hierarchical Bayesian analyses are then performed on catalogues of 100-5,000 events. With hundreds of detections from ET, we will be able to establish that the total mass distribution is bimodal. A few thousand events are sufficient to disentangle the redshift distributions of the two subpopulations for moderate time-delay indices (-0.5 or -1). For steeper indices (-1.5), the differences are more subtle, requiring larger catalogues, beyond what we could explore given our computational resources. Next-generation detectors should enable the detection of multiple BNS subpopulations and their redshift evolution, providing valuable insight into their formation pathways.

A complete overview of the stellar, gas and dust contents of galaxies is key to understanding their assembly at early times. However, an estimation of molecular and atomic gas reservoirs at high redshift relies on various indirect tracers, while robust dust mass measurements require multi-band far-infrared continuum observations. We take census of the full baryonic content of the main-sequence star-forming galaxy HZ10 at $z=5.65$, a unique case study where all necessary tracers are available. We present new ALMA Band 10 ($\lambda_\mathrm{rest}=50\mu$m) and Band 4 ($300\mu$m) observations towards HZ10, which combined with previously taken ALMA Band 6 through 9 data ($70-200\mu$m) constrains its dust properties. We complete the baryonic picture using archival high-resolution [CII] observations that provide both a dynamical mass and molecular and atomic gas mass estimates, a JVLA CO(2-1)-based molecular gas mass, and JWST metallicity and stellar mass measurements. We detect continuum emission from HZ10 in Bands 10 and 4 at the $3.4-4.0\sigma$ level, and measure a dust temperature of $T_\mathrm{dust} = 37_{-5}^{+6}$K and dust mass $\log(M_\mathrm{dust}/M_\odot) = 8.0 \pm 0.1$. Leveraging the dynamical constraints, we infer its total gas budget, and find that commonly used [CII]-to-H$_2$ and [CII]-to-HI conversions overpredict the gas mass relative to the dynamical mass. For this reason, we derive a [CII]-to-total ISM mass (atomic + molecular) conversion factor, which for HZ10 corresponds to $\alpha_\mathrm{[CII]}^\mathrm{ISM} = 39^{+50}_{-25}M_\odot L_\odot^{-1}$. We also find that HZ10 falls below the local scaling relation between dust-to-gas ratio and metallicity, suggesting inefficient ISM dust growth. These results demonstrate a powerful synergy between ALMA and JWST in disentangling the baryonic components of early galaxies, paving the way for future studies of larger samples.

Detecting Earth twins remains challenging because their shallow, long-period transits are difficult to distinguish from background noise. Motivated by the challenge, we developed Segmented-Polynomial-fitting Least Squares (SPLS), a new algorithm that simultaneously fits planetary transits and background trends using a segmented double polynomial model. Prior to signal detection, the optimal polynomial order for the trend component is selected using Bayes factor-based model comparison. During the periodogram search, the Signal Detection Efficiency metric is used to assess signal significance. The algorithm is accelerated by a three-step approximation and nonlinear parameter sampling tailored to SPLS. We compare the performance of SPLS with traditional detrending-detection approaches across different orbital periods, signal-to-noise ratios (SNR), planet radii and stellar noise levels on an injection-recovery test. When detecting signals with periods between 10 and 480 days and SNRs below 9, SPLS achieves at least a 22.6% higher true positive rate than other methods at the same 10% false positive rate. Using the threshold determined from the Receiver Operating Characteristic curve analysis, our method also recovers the most true signals while yielding the fewest false positives among all injected samples, and reaches a 97% recovery fraction in Kepler confirmed single-planet systems. The tests demonstrate that SPLS improves the detection of transiting planets, particularly for low-SNR, long-period signals. It offers the potential for finding Earth twins in future applications to data from Kepler, TESS, and upcoming PLATO and Earth 2.0 missions.

Fast Radio Bursts (FRBs) have become an indispensable tool for studying the Universe's ionisation properties, as well as its cosmological parameters. This is achieved by analysing their diffuse dispersion measure (${\rm DM}_{\rm diff}$) as a function of redshift. However, this requires an integration along the line-of-sight, which is time-consuming. In this work, we derive an analytical approximation formula for ${\rm DM}_{\rm diff}$ for flat, $\Lambda$CDM and $w$CDM universes. We show that our approximation works well for the ranges $0.01 \leq z \leq 2$, $0.2 \leq \Omega_m \leq 1.0$ and $-3.0 \leq w \leq -0.5$, with relative fractional error to a numerically evaluated ${\rm DM}_{\rm diff}$ always smaller than $3.5\ \%$, in the worst case scenario, and in most cases smaller than $0.5\ \%$. Moreover, the approximation is more than $15$ ($2$) times faster than the numerical solution of $\Lambda$CDM ($w$CDM). Therefore, we hope that it could be a useful tool for the FRB community.

M31 UCXB-1 is one of the brightest X-ray point sources in the bulge of M31, with a peak X-ray luminosity $ L_{\mathrm{0.5-10 \: keV}} = 2.9^{+0.2}_{-0.2} \times 10^{38} \: \mathrm{erg} \: \mathrm{s}^{-1} $. Both XMM-Newton and Chandra observations have detected an eclipsing signal with a period of about 465 seconds from this source, and we note that the periodic signal is detected exclusively during the source's high-luminosity states. This signal probably originates from its orbital motion, therefore it is an ultra-compact X-ray binary (UCXB) candidate with the highest X-ray luminosity. Our theoretical analyses show that M31 UCXB-1 is in good agreement with the luminosity-orbital period relation ($ L_{\mathrm{2-10 \: keV}}-P_{\mathrm{orb}} $) of the black hole/neutron star--white dwarf (BH/NS--WD) UCXB system. Moreover, our spectral analyses indicate that the primary in M31 UCXB-1 is more likely to be a BH rather than an NS. The results show that M31 UCXB-1 is a BH--WD system, with the shortest orbital period, the possibly strongest gravitational wave emission, and the most massive white dwarf among the known UCXBs.

We have applied our approach of using ''SNe Ia twins''in the Hubble flow to obtain distances to SNe Ia at z $>$ 0.015 and derive H$_{0}$. Our results, taking a single step between the low z domain and the Hubble flow, validate the three rung classical method. We find, however, that the full compilation of distances, both in Pantheon+ and in the Carnegie-Chicago Hubble Program (CCHP), contain some inaccurate values in the colors due to an underestimate of reddening by dust. This produces odd individual values for H$_{0}$ from single SNe Ia. On the average, those erroneous estimates do not affect the mean value of H$_{0}$, which is characterized by the bulk of well--modeled SNe Ia. Our sample of carefully addressed SNe Ia in the Hubble flow contains a dozen supernovae, for which the distances are determined with high accuracy. Three of these SNe Ia are of the Broad Line subtype and can be compared with SN 1989B in M66, a host galaxy with a unique convergence of the Cepheid distance determination and the Tip of the Red Giant Branch stars (TRGB) determination by the CCHP group. They point to a weighted average of H$_{0}$ $=$ 73.556 $\pm$ 2.084 (stat) km s$^{-1}$ Mpc $^{-1}$. There is as well a very good agreement on the distances to NGC 7250 and NGC 5643 between those derived with Cepheids by SH0ES and those derived with the use of J-Asymptotic Giant Branch stars (JAGB stars) by the CCHP, which makes them very good anchors. The sample of 12 SNe Ia gives a value of H$_{0}$ $=$ 72.833 $\pm$ 1.306(stat) $\pm$ 1.151 (sys) km s$^{-1}$ Mpc$^{-1}$, when anchored in Cepheids, and of H$_{0}$ $=$ 72.388 $\pm$ 1.272 (stat) $\pm$ 1.015 (sys) km s$^{-1}$ Mpc$^{-1}$, when anchored in JAGBs by the CCHP. We take a mean of the two values of H$_{0}$ and obtain H$_{0}$ $=$ 72.610 $\pm$ 1.289(stat) $\pm$ 1.085 (sys) km s$^{-1}$ Mpc$^{-1}$.

Lunar laser ranging (LLR) currently delivers mm-class tests of relativistic gravity and the lunar interior, but further gains are limited by photon-starved pulsed systems, array-induced pulse broadening, and atmospheric variability. Pushing into the sub-0.1mm regime is required to probe ultra-slow tidal dissipation, long-period librations, uHz-band gravitational waves. We develop a complete amplitude-modulated continuous-wave (AM-CW) framework for high-power LLR, including measurement models, estimators, a unified covariance description, and an observatory CONOPS. CW LLR is treated as RF phase metrology on a bright CW optical carrier: 2-way range and 1-way line-of-sight range-rate obtained from the phase and phase slope of a GHz-class modulation tone. The joint estimators are characterized by a covariance matrix decomposed into photon, atmospheric, instrumental, oscillator, nonlinearity contributions, mapped into subsystem requirements, including a closed-form range-rate bias from modulation-frequency slew across the ~2.56s round trip and synthetic-wavelength constraints on multi-tone nonlinearity. Using realistic photon budgets for a 1 kW, 1064 nm CW transmitter on a 1-m telescope ranging to 10cm corner-cube retroreflectors, we show that the photon-statistical floor for 2-way range at a GHz-class modulation frequency is at the level of a few-10 um in ~100s. With allocations to residual atmosphere and instrument, the resulting single-station 2-way precision is ~80 um, with range-rate sensitivities of ~0.1-1 um/s. Differential LLR between reflectors separated on the Moon by ~1000km can deliver 20-50 um and few (0.1-1) um/s differential range and range-rate. These precisions are realized across three representative regimes that span detected photon rates of (5-7) x 10^3-10^5 s^{-1} and are translated directly into hardware and CONOPS requirements for next-generation high-power AM-CW LLR.

It has been suggested that there are $\sim 10^5$ black hole-neutron star (BH-NS) binaries in our Galaxy. However, despite the effort of intensive radio search for decades, none of these binaries has been found to date. These binaries are regarded as a holy grail of astronomy because they can greatly improve our understanding about relativistic systems of compact objects and fundamental physics. In this article, we propose the existence of exotic BH-NS binaries which can open a new way in searching the missing BH-NS binaries in our Galaxy. By considering the possible dark matter density spikes formed around the primordial black holes in the BH-NS binaries, we show that extremely high temperature ($\sim 10^6$ K) could be maintained on the surface of the neutron stars due to effective dark matter capture. This interesting feature can also help reveal the nature of dark matter and possibly further improve the upper limit of dark matter scattering cross section well below $10^{-47}$ cm$^2$.

We present the redshift evolution of radial color gradients (in rest-frame ${\textit{U}} - {\textit{V}}$ and ${\textit{V}} - {\textit{J}}$) for galaxies in the range $0.5<$ z $<2.5$ and investigate their origin and dependence on stellar mass. We select $\sim 10,200$ galaxies with stellar masses $M_\star>10^{9.5}~{\text{M}}_\odot$ from publicly available JWST/NIRCam-selected catalogs. Using 2D Sérsic profile fits to account for PSF broadening, we perform spatially resolved SED fitting on HST and JWST/NIRCam photometry retrieving accurate rest-frame ${\textit{U}} - {\textit{V}}$ and ${\textit{V}} - {\textit{J}}$ color gradients within 2$R_\text{e, F444W}$. Star-forming galaxies generally exhibit negative ${\textit{V}} - {\textit{J}}$ color gradients that are strongly mass and redshift dependent. For massive star-forming galaxies ($M_\star>10^{10.5}~{\text{M}}_\odot$) at $z>1.5$ ${\textit{V}} - {\textit{J}}$ colors are $\approx 0.5$ mag redder within the effective radius than outside, on average. We find that, at all redshifts and across the entire stellar mass range, ${\textit{V}} - {\textit{J}}$ gradients strongly correlate with global attenuation ($A_V$), suggesting that they predominantly trace dust attenuation gradients. Edge-on galaxies are redder and have stronger gradients at all $z$, although the correlation weakens at higher $z$. The ${\textit{U}} - {\textit{V}}$ and ${\textit{V}} - {\textit{J}}$ color gradients in the quiescent galaxy population, in contrast, are weakly negative (from $\approx -0.1$ to $\approx- 0.2$ mag), though significant, and show little or no dependence on stellar mass, redshift or axis ratio. The implication is that quiescent galaxies must be largely transparent, with low $A_V$, and color gradients mostly attributable to stellar population gradients.

Recent observations show a persistent tension in the Hubble constant $H_0$, suggesting an incomplete understanding of cosmic expansion and local dynamics. Using the Cosmicflows-4 catalogue, we mapped the angular and radial variations of $H_0$ in radial shells with a distance modulus $\mu \in [29,36]$ (approximately corresponding to $[20,100]$ $h^{-1}$ Mpc) and equal-area sky patches, applied adaptive weighing, and fitted spherical harmonics up to the octupole. Our results reveal a clear, distance-decreasing dipole that remains coherent across shells, with subdominant higher-order multipoles, and the octupole fit capturing the main anisotropic features except in sparsely sampled or SDSS-dominated shells. The direction and amplitude of the dipole depend on whether the observed radial recessional velocities are corrected for peculiar velocities or not. If the correction is not applied, the dipole aligns with the major gravitational structures in the local universe. If it is applied, a global dipole still seems to be present, but the signal is much weaker and with much lower statistical significance. This decrease in the amplitude of the dipole supports the idea of a differential expansion rate in our universe, but does not clarify whether the origin is astrophysical or cosmological. Finally, we verify that, while this anisotropy could influence local measurements of the Hubble constant, its effect on the large-scale Hubble tension appears to be limited, as the distribution of galaxies hosting SNeIa, both used as calibrators to constrain $H_0$ and in the Hubble-flow, does not show a strong correlation with the dipole signal.

The heliospheric current sheet (HCS) is an important large-scale structure of the heliosphere, and, for the first time, the Parker Solar Probe (PSP) mission enables us to study its properties statistically close to the Sun. We visually identify the 39 HCS crossings measured by PSP below 50~\Rs~during encounters 6 to 21, and investigate the occurrence and properties of magnetic reconnection, the behavior of the spectral properties of the turbulent energy cascade, and the occurrence of kinetic instabilities at the HCS. We find that 82\% of HCS crossings present signatures of reconnection jets, showing that the HCS is continuously reconnecting close to the Sun. The proportion of inward/outward jets depends on heliocentric distance, and the main HCS reconnection X-line has a higher probability of being located close to the Alfvén surface. We also observe a radial asymmetry in jet acceleration, where inward jets do not reach the local Alfvén speed, contrary to outward jets. We find that turbulence levels are enhanced in the ion kinetic range, consistent with the triggering of an inverse cascade by magnetic reconnection. Finally, we highlight the ubiquity of magnetic hole trains in the high $\beta$ environment of the HCS, showing that the mirror mode instability plays a key role in regulating the ion temperature anisotropy in HCS reconnection. Our findings shed new light on the properties of magnetic reconnection in the high $\beta$ plasma environment of the HCS, its interplay with the turbulent cascade and the role of the mirror mode instability.

The Large High Altitude Air Shower Observatory (LHAASO) has recently reported five Galactic microquasars as Ultra-High-Energy (UHE) $\gamma$-ray emitters (> 100 TeV). Among these sources, the microquasar V4641 Sgr exhibits $\gamma$-ray emission up to $\sim$0.8 PeV, requiring the acceleration of particles to multi-PeV energies, as well as the hardest UHE spectrum. The mechanisms behind particle acceleration to such energies are not well understood. Furthermore, the limited multi-wavelength (MWL) information on this source appears contradictory, further complicating interpretation and suggesting that V4641 Sgr may represent a particularly unusual case. In this work, we present a detailed physical model of V4641 Sgr that combines first-principles simulations of stochastic (turbulent) particle acceleration with MWL emission modeling. We adopt a leptonic scenario in which electrons are accelerated via the second-order Fermi process driven by relativistic strong turbulence ($\delta B/B \sim 1$). The particle energization is simulated using a dedicated Monte Carlo framework STRIPE that incorporates the effects of intermittent energy gains and radiative losses. The resulting accelerated electrons produce UHE $\gamma$-rays through inverse Compton scattering on both the cosmic microwave background (CMB) and the interstellar radiation fields (ISRF). Our model is capable of reproducing key observational characteristics of the system, including particle acceleration to energies of tens of PeV, as well as the TeV-PeV $\gamma$-ray spectrum and the hard spectral index measured by LHAASO. Nonetheless, several aspects remain unresolved, highlighting the need for deeper observational coverage and further theoretical refinement.

Chemical models are routinely used to predict the atmospheric composition of exoplanets and compare it with the composition retrieved from observations, but little is known about the reliability of the calculated composition. We carried out a sensitivity analysis to quantify the uncertainties in the abundances calculated by a state-of-the-art chemical atmosphere model of the widely observed planets WASP-33b, HD209458b, HD189733b, WASP-39b, GJ436b, and GJ1214b. We found that the abundance uncertainties in the observable atmosphere are relatively small, below one order of magnitude and in many cases below a factor of two, where vertical mixing is a comparable or even larger source of uncertainty than (photo)chemical kinetics. In general, planets with a composition close to chemical equilibrium have smaller abundance uncertainties than planets whose composition is dominated by photochemistry. Some molecules, such as H2O, CO, CO2, and SiO, show low abundance uncertainties, while others such as HCN, SO2, PH3, and TiO have more uncertain abundances. We identified several critical albeit poorly constrained processes involving S-, P-, Si-, and Ti-bearing species whose better characterization should lead to a global improvement in the accuracy of models. Some of these key processes are the three-body association reactions S + H2, Si + O, NH + N, and N2H2 + H; the chemical reactions S + OH --> SO + H, NS + NH2 --> H2S + N2, P + PH --> P2 + H, and N + NH3 --> N2H + H2; and the photodissociation of molecules such as P2, PH2, SiS, CH, and TiO.

Fast Radio Bursts (FRBs) are energetic radio bursts that typically last for milliseconds. They are mostly of extragalactic origin, but the progenitors, trigger mechanisms and radiation processes are still largely unknown. Here we present a comprehensive analysis on 415 non-repeating FRBs detected by CHIME, applying manual filtering to ensure sample completeness. It is found that the distribution of fluence can be approximated by a three-segment power-law function, with the power-law indices being $-3.76 \pm 1.61$, $0.20 \pm 0.68$ and $2.06 \pm 0.90$ in the low, middle, and high fluence segments, respectively. Both the total dispersion measure (\text{DM}) and the extragalactic \text{DM} follow a smoothly broken power-law distribution, with characteristic break DM values of $\sim 703$ pc $\rm cm^{-3}$ and $\sim 639$ pc $\rm cm^{-3}$, respectively. The redshifts are estimated from the extragalactic \text{DM} by using the Macquart relation, which are found to peak at $ z \sim 0.6$. The isotropic energy release ($E_{\text{iso}}$) is also derived for each burst. Two-Gaussian components are revealed in the distribution of $E_{\text{iso}}$, with the major population narrowly clustered at $\sim 2.3 \times 10^{40} {\rm erg}$. The minor population have a characteristic energy of $\sim 1.6 \times 10^{39}$ erg and span approximately one order of magnitude. The distribution hints a near-uniform energy release mechanism for the dominant population as expected from some catastrophic channels, whereas the lower-energy component (potentially including repeat-capable sources) may reflect a broader diversity in FRB origins, emission mechanisms and evolutionary stages.

We analyse imaging (EUV, X-ray) and spectral (radio, X-ray) data obtained by ground based and space instruments on board space missions both on Earth (Fermi, Hinode, Solar Dynamics Observatory) and solar orbit (Solar Orbiter, STEREO-A), which provide a multi-directional view on the same event. The combination of EUV and X-ray images and X-ray spectra allowed us to identify hot loops in the vicinity of the filament before its eruption. We interpreted their interaction with the rising filament as a signature of an arcade-to-rope reconnection geometry. The subsequent EUV brightening within the filament revealed helical structure of the erupting rope. We explained co-temporal radio slowly positively drifting bursts as a result of beam acceleration within the magnetic rope and propagation along the helical structure. Corresponding X-ray spectra were consistent with a thermal origin. The filament rising was accompanied by co-temporal normal and reverse drift type III radio bursts. We interpreted them as a signature of a reconnection event and estimated electron density at the reconnection site. Further untwisting of the helical structure led to formation of a quasi-circular EUV structure seen from Earth and STEREO-A. Its occurrence was co-temporal with a unique tangle of radio U- and inverse U-bursts. We proposed that several accelerated beams propagate within that complex structure and generate the burst tangle. During the start of the flare hard X-ray emission was concentrated near the filament leg only suggesting predominant propagation of the beams towards its rooting. We collected multi-wavelength observations indicating interaction of the erupting magnetic flux rope with the overlying arcade and internal magnetic reconnection inside the rising flux rope.

The well-known correlation between stellar metallicity and planet occurrence is strongest for giant planets, but weaker for smaller planets, suggesting that detailed elemental patterns beyond [Fe/H] may be relevant. Using abundances from the fourth data release of the GALAH spectroscopic survey, we analyzed 104 host stars with 141 confirmed transiting planets. We divide planets at 2.6 Earth radii, the theoretical threshold radius above which planets are unlikely to be pure-water worlds. We find that large-planet hosts are enriched by approximately 0.2 dex in iron and show a possible excess of highly volatile elements (C, N, O), though these measurements are affected by observational limitations, whereas small-planet hosts exhibit an enhanced contribution of the classical rock-forming elements (Mg, Si, Ca, Ti) relative to iron, corresponding to a modest [Rock/Fe] offset of 0.06 dex, which is statistically significant, with a p value of 10^{-4}. These offsets remain significant for alternative radius cuts. A matched control sample of non-planet-host stars shows only weak and mostly statistically insignificant similar trends, confirming that the stronger chemical signatures are linked to the planetary characteristics. As our study relies on transiting planets, it mainly probes short-period systems (periods shorter than 100 days). These results refine the planet-metallicity relation, highlighting the role of the relative balance between iron, volatiles, and rock-forming elements in planet formation.

We present a systematic survey of strange mode pulsations in Cepheids using MESA and for the linear stability analysis, MESA RSP. Our model grid spans $2-15\,M_\odot$ in mass and [Fe/H] $= -0.95$-$0.17$ ($Z = 0.0015$-$0.0200$) in metallicity, with four convective overshoot prescriptions. Strange modes were identified in a relatively small fraction ($5-12.5$ %) of models, occurring at $n_\mathrm{pg} = 5-9$, with $n_\mathrm{pg} = 6$-$7$ as the most frequent radial modes. No unstable solutions were identified beyond $n_\mathrm{pg} = 9$, in contrast to earlier studies reporting strange modes at $n_\mathrm{pg} = 10-12$. We quantified the duration of the instability crossing phase ($\tau_\mathrm{IS}$), the strange mode phase ($\tau_\mathrm{s}$), and their ratio $\mathcal{P}_\mathrm{s} = \tau_\mathrm{s} / \tau_\mathrm{IS}$. Toward higher masses, both $\tau_\mathrm{IS}$ and $\tau_\mathrm{s}$ decrease, yet their ratio shows no systematic trend with mass in models that include convective core overshoot. The absolute timescales for strange modes remain short, typically $\tau_\mathrm{s} \sim 10^{4.5}$-$10^{6}$ years, while $\tau_\mathrm{IS}$ is often an order of magnitude shorter, implying that these stars may spend a larger fraction of their life in the strange mode phase than in the instability strip itself. The extended duration of the strange mode phase may enhance the detectability of strange mode pulsators, provided that observational precision is sufficient to capture their low-amplitude variability. The predicted periods ($0.6-6.3$ days) are well covered by a single $27-$day TESS sector, making strange mode pulsators potentially detectable with current space-based photometry, although blending with nearby sources may pose challenges.

Type-2 quasars (QSO2s) are AGN seen through a significant amount of dust and gas that obscures the central supermassive black hole and the broad line region. Despite this, recent mid-infrared spectra of the central 0.5-1.1 kpc of five QSO2s at z~0.1, obtained with the MRS module of JWST/MIRI, revealed 9.7, 18, and 23 micron silicate features in emission in two of them. This indicates that the high angular resolution of JWST/MIRI now allows us to peer into their nuclear region, exposing some of the directly illuminated dusty clouds that produce silicate emission. To test this, we fitted the nuclear mid-infrared spectrum of the QSO2 with the strongest silicate features, J1010, with two different sets of torus models implemented in an updated version of the Bayesian tool {\tt BayesClumpy}. These are the CLUMPY and the CAT3D-WIND models. The CAT3D-WIND model is preferred by the observations based on the marginal likelihood and fit residuals, although the two torus models successfully reproduce the spectrum by means of intermediate covering factors ($\rm C_T=0.45\pm^{0.26}_{0.18}$ and $\rm C_T=0.66\pm^{0.16}_{0.17}$ for the CLUMPY and CAT3D-WIND models) and low inclinations ($\rm i=50^\circ\pm^{8^\circ}_{9^\circ}$ and $\rm i=13^\circ\pm^{7^\circ}_{6^\circ}$). Indeed, four of the five QSO2s with JWST/MIRI observations, including J1010, are in the blowout or ''forbidden'' region of the Eddington ratio-column density diagram, indicating that they are actively clearing gas and dust from their nuclear regions, leading to reduced covering factors. This is in contrast with Seyfert 2 galaxies observed with JWST, which are in the ''permitted'' regions of the diagram and show 9.7 micron silicate features in absorption. This supports a scenario where the more luminous the AGN and the higher their Eddington ratio, the lower the torus covering factor, driven by radiation pressure on dusty gas.

The cosmological distribution and formation rate of compact astrophysical objects such as fast radio bursts (FRBs) are typically assumed to be proportional to a linear combination of cosmological star formation rate and stellar mass. In the literature, a template for star formation rate, which is just a function of redshift, is typically used. In this work, we point out the importance of galaxy stellar mass function which captures the host galaxy information of observed FRBs as well as the redshift evolution of galaxy stellar mass. Using this information, we find that FRB formation efficiency per stellar mass has to be more efficient (by a factor of $\approx 3$) than previously calculated, in order to reproduce the observed volumetric rate of FRBs at $z=0$. We show that cosmological population studies of FRBs have to include host galaxy information along with its redshift evolution in order to obtain unbiased results. This consideration is also applicable to other transients, e.g. gamma-ray bursts and merging binary black hole events. We show that our approach may open up the possibility to infer the host galaxy stellar mass of merging binary black holes with a detection of few thousand gravitational wave events.

Blazars are promising targets for neutrino astronomy, as highlighted by IceCube's identification of TXS 0506+056 as a cosmic neutrino source candidate. High-frequency-peaked BL Lacs (HBLs) stand out due to their characteristic electromagnetic high-energy emission properties, which makes them as well promising candidates for the production of high-energy neutrinos. Such neutrinos could be detected by KM3NeT/ARCA, a next-generation deep-sea Cherenkov detector under construction in the Mediterranean Sea. This contribution presents the results of a binned likelihood stacking analysis investigating high-energy neutrino emissions from a subset of HBLs. A total of 232 sources have been selected from the 3HSP catalogue and targeted in the analysis of KM3NeT/ARCA data taken with several detector configurations, specifically in the period in which from 6 to 21 detection units were deployed. The analysis is based on neutrino emission models of blazars developed using the LeHa-Paris numerical code, which simulates the spectral energy distribution of jetted active galactic nuclei.

A class of supernovae (SNe) termed `SN Ic-CSM' are characterized by late-time emergence of narrow emission lines of elements formed in the oxygen core of a massive star. A popular scenario is the interaction of the SN ejecta and O-rich circumstellar medium (CSM), i.e., Circumstellar Interaction (CSI). Uncovering the progenitor system of SNe Ic-CSM plays a critical role in understanding the final evolution of a massive star to a bare C+O star. In this Letter, we present observations of SN 2022esa which we show is an SN Ic-CSM. Surprisingly, a stable periodicity of ~32 days is found in its light-curve evolution with a hint of a slowly increasing period over ~200 days. We argue that the main power source is likely the interaction of the SN ejecta and O-rich CSM, while the energy input by the post-SN eccentric binary interaction within the SN ejecta is another possibility. In either case, we propose a massive Wolf-Rayet (WR) star as the progenitor, in a WR-WR or WR-BH (black hole) binary that will eventually evolve to a BH-BH binary. Specifically, in the CSI scenario, the progenitor system is an eccentric binary system with an orbital period of about a year, leading to the observed periodicity through the modulation in the CSM density structure. We also show that some other objects, superluminous SN I 2018ibb (a pair-instability SN candidate) and peculiar SN Ic 2022jli (the first example showing stable periodic modulation), show observational similarities to SNe Ic-CSM and may be categorized as SN Ic-CSM variants. Complemented with a large diversity in their light-curve evolution, we propose that SNe Ic-CSM (potentially linked to SNe Ibn/Icn) are a mixture of multiple channels that cover a range of properties in the progenitor star, the binary companion, and the binary orbit.

Fast Radio Bursts (FRBs) are millisecond-duration radio pulses originating from cosmological distances, as indicated by their large dispersion measures. While numerous FRBs have now been localised to their host galaxies, a distinct class of compact electromagnetic counterpart, a Persistent Radio Source (PRS), has also been identified in some cases. Currently, only three, and possible a fourth repeating FRBs (FRB20121102A, FRB20190417A, FRB20190520B, and FRB20240114A) have confirmed associations with a PRS. Insight into progenitors, local environments, and the evolution of FRBs can be clarified by characterising these PRSs. In this work, we present 2 detected candidate PRSs using MeerKAT radio telescope data and one non-detection (as part of a larger study involving 25 FRB positions). Both FRB20221106 and FRB20181112 were found to have a host galaxy, and whether the detected radio continuum emission comes from the host galaxy or PRS is still an open question. High-resolution observations from a telescope such as e-MERLIN are required to resolve this question. If a compact PRS is detected, this telescope will provide the size, and investigate the flux variability and spectral shape of this compact PRS. Lastly, in the case of FRB220190102, which was observed over two epochs, no radio continuum was detected. However, a flux upper limit is provided for both epochs.

We update the status of the spectropolarimetric campaign dedicated to characterise the magnetic field properties of a sample of known exoplanet-hosting stars included in the current target list of the $Ariel$ mission. We analysed spectropolarimetric data collected for 15 G-M type stars with Neo-Narval, HARPSpol, and SPIRou to assess the detectability of the large-scale magnetic field. For three stars we reconstructed the magnetic field topology and its temporal evolution via Zeeman-Doppler imaging (ZDI). Such reconstructions were then used to perform 3D MHD simulations of the stellar wind and environment impinging on the hosted exoplanets. We detected the magnetic field of six stars. Of these, we performed ZDI reconstructions for the first time of TOI-1860 and DS Tuc A, and for the second time of HD 63433, providing temporal information of its large-scale magnetic field. Consistently with previous results on young ($\rm\sim 50-100~Myr$) solar-like stars, the large-scale magnetic field is moderately strong (30-60 G on average) and complex, with a significant fraction of magnetic energy in the toroidal component and high-order poloidal components. We found the orbit of TOI-1860 b to be almost completely sub-Alfvénic, the orbits of DS Tuc A b and HD 63433 d to be trans-Alfvénic, and the orbits of HD 63433 b and c to be super-Alfvénic. We obtained marginal detections of the magnetic field for TOI-836 and TOI-2076, and detections for TOI-1136, but the number of observations is not sufficient for magnetic mapping. A magnetic star-planet connection can occur for most of TOI-1860 b's orbit. This can happen more sporadically for DS Tuc A b and HD 63433 c given the lower fraction of their orbit in the sub-Alfvénic regime. The orbit of HD 63433 c is nevertheless more sub-Alfvénic than previously simulated owing to the temporal evolution of the stellar magnetic field.

The chemical abundance of the interstellar medium sets the initial conditions for star formation and provides a probe of chemical galaxy evolution models. However, unresolved inhomogeneities in the electron temperature can lead to a systematic underestimation of the abundances. We aim to directly test this effect. We use the SDSS-V Local Volume Mapper to spatially map the physical conditions of the Trifid Nebula (M 20), a Galactic H II region ionized by a single mid-type O star, at 0.24 pc resolution. We exploit various emission lines (e.g., Hydrogen recombination lines and collisionally excited lines, including also faint auroral lines) and compute spatially resolved maps of [O II] and [S II] electron densities; [N II], [O II], [S II], [S III] electron temperatures; and the ionic oxygen abundances. We find internal variations of electron density that result from the ionization front, along with a negative radial gradient. However, we do not find strong gradients or structures in the electron temperature and the total oxygen abundance, making the Trifid Nebula a relatively homogeneous H II region at the observed spatial scale. We compare these spatially resolved properties with equivalent integrated measurements of the Trifid Nebula and find no significant variations between integrated and spatially resolved conditions. This isolated H II region, ionized by a single O-star, represents a test case of an ideal Strömgren sphere. The physical conditions in the Trifid Nebula behave as expected, with no significant differences between integrated and resolved measurements.

We examine the connection between galaxy mergers and the triggering of active galactic nuclei (AGNs) using a sample of 614 type 1 AGNs at $z<0.07$, along with a control sample of inactive galaxies matched to the AGNs for comparison. We used tidal features, detected in deep images from the DESI Legacy Imaging Survey, as direct evidence of recent mergers. We find that the fraction of type 1 AGN hosts with tidal features ($f_T$) is higher for AGNs with higher luminosities and (to a lesser extent) more massive black holes. Specifically, $f_T$ rapidly increases from $0.05\pm0.03$ to $0.75\pm0.13$ as the luminosity of the $[\mathrm{O\,{\scriptstyle III}}]$ $\lambda$5007 emission line ($L_{[\mathrm{O\,{\scriptstyle III}}]}$), an indicator for bolometric AGN luminosity, increases in the range $10^{39.5}\lesssim L_{[\mathrm{O\,{\scriptstyle III}}]}/(\mathrm{erg\,s}^{-1}) \lesssim10^{42.5}$. In addition, $f_T$ increases from $0.13\pm0.03$ to $0.43\pm0.09$ as black hole mass ($M_\mathrm{BH}$) increases in the range $10^{6.0}\lesssim M_\mathrm{BH}/M_{\odot}\lesssim10^{8.5}$. The fraction $f_T$ also increases with the Eddington ratio, although the trend is less significant compared to that with $L_{[\mathrm{O\,{\scriptstyle III}}]}$ and $M_\mathrm{BH}$. The excess of $f_T$, defined as the ratio of $f_T$ for AGNs to that of their matched inactive counterparts, exhibits similar trends, primarily increasing with $L_{[\mathrm{O\,{\scriptstyle III}}]}$ and weakly with $M_\mathrm{BH}$. Our results indicate that, in the local Universe, galaxy mergers are the predominant triggering mechanism for high-luminosity AGNs, whereas they play a lesser role in triggering lower-luminosity AGNs. Additionally, strong events, such as galaxy mergers, may be more necessary to activate massive black holes in more massive galaxies due to their lower gas fractions.

Very high-energy (VHE) detected flat-spectrum radio quasars (FSRQs) are relatively few despite being the most persistent bright MeV-GeV sources. Focusing on VHE emission, we investigated the spectral and temporal properties of VHE-detected FSRQs using 14-year Fermi-LAT data. All are highly variable (flux-amplitude$>$100) with VHE detection associated with brighter flux states and relatively harder spectra. Above a flux limit, flux anti-correlates with spectral index, exhibiting a bluer-when-brighter trend. The low-flux state spectral energy distributions (SEDs) for all resembles a power-law, while high-flux and VHE-associated states resemble a log-parabola, accompanied by an almost nil (PKS0736+017, PKS1510-089) to marginal (4C+21.35, 3C279) to significant (B21420+32, TON0599, PKS1441+25, S30218+35, PKS0346-27, OP313) MeV-GeV peak-upshift -- more prominent in high-redshift sources. For no/marginal peak-upshift, the VHE emission is consistent with external Comptonization of infrared photons (EC-IR) driven primarily by a power-law continuation of the particle spectrum to higher energies. For those with a significant MeV-GeV peak-upshift, PKS0346-27 and OP313 shows peak-upshift in the synchrotron spectrum, and thus VHE is EC-IR origin, while for others without synchrotron-peak upshift, we attribute the VHE to an HBL-like component with a Compton-Dominance (CD) like FSRQs, with VHE driven primarily by particle spectrum continuation. In some, even high-state SEDs seem to require an HBL-like component. Thus, VHE activities in FSRQs mainly result from particle spectrum continuation, aided by spectral transition or a new HBL-like component with FSRQ-like CD. Such spectral changes naturally brighten the GeV-VHE flux, overcoming extragalactic background light absorption without requiring extraordinary brightening under the traditional EC-IR scenario than normally exhibited.

A very recent and meticulous timing effort suggests that an eccentric millisecond pulsar (eMSP) binary, namely PSR J1618-3921, is likely to be a part of a hierarchical triple (HT) system with a $0.6M_{\odot}$ companion in a $\sim 300$yr orbit. We investigate observational implications of the proposed HT scenario for PSR J1618-3921 and our ability to constrain the scenario. We model the MSP-Helium White Dwarf binary to be a part of bound point-mass HT, while incorporating the effects due to the quadrupolar interactions between the inner and outer binaries, along with dominant order general relativistic contributions to the periastron precession of the inner binary. If the proposed HT system is indeed undergoing Kozai oscillations at the present epoch, the orbital eccentricity ($e$) would be expected to decrease, while the rate of periastron advance ($\dot{\omega}$) would correspondingly increase, for plausible ranges of the HT parameters. Furthermore, the fractional variations in $e$ are anticipated to be at the level of a few parts in $10^{5}$-a magnitude that is substantially larger than the current measurement precision of $e$. We find that, for this eccentric MSP binary, the HT configurations that minimize the temporal evolution of orbital eccentricity and argument of periastron are mutually incompatible. This indicates that the continued high-precision timing of PSR J1618-3921-when analyzed within the framework introduced here-should place stringent limits on the presence and properties of a potential third body in the system.

The Extragalactic Background Light (EBL) is the accumulated light emitted throughout the history of the universe, spanning the UV, optical, and IR spectral ranges. Stars and dust in galaxies are expected to be the main source of the EBL. However, recent direct measurements performed beyond Pluto's orbit (less affected by foregrounds than those performed from the Earth) hint at an EBL level in the optical band larger than the one expected from the integrated contribution of known galaxy populations. One approach that could solve this controversy uses Very High Energy (VHE) photons coming from sources at cosmological distances. These photons can interact with the EBL producing electron-positron pairs, a process that leaves an imprint on the observed VHE spectrum. This technique, however, requires assumptions on the intrinsic spectrum of the source, which can compromise the robustness of EBL constraints. In this contribution, we used Monte Carlo simulations and archival data of the MAGIC telescopes to study the impact that the assumptions adopted in the literature have on the estimate of the EBL density, and how using more generic ones would modify the results. Our results show how the EBL density constraints obtained highly depend on the intrinsic spectral shape assumed for the source. We have studied two different methods to reduce the assumptions on the intrinsic spectral shape to get more robust results. This will be especially important for upcoming with new VHE facilities, where systematic uncertainties are expected to play a more significant role compared to statistical ones.

In dense environments, mechanisms like ram pressure stripping (RPS) and gravitational interactions can induce similar morphological features in galaxies, distinguishable only through detailed study of their stellar properties. While RPS affects recently formed stars by displacing the gas disk from which they form, gravitational interactions perturb stars of all ages rather similarly. We present the first observational test of the Size-Shape Difference (SSD) measure, a novel approach validated for simulated galaxies, that quantifies morphological differences between young and intermediate-age stellar populations to distinguish RPS from gravitationally interacting galaxies. We analyze 67 galaxies from the GASP survey using spatially-resolved star formation histories derived using SINOPSIS. In our fiducial model, we compare stellar populations in two age bins (t < 20 Myr and 20 Myr <= t < 570 Myr) to calculate SSD values. The sample includes confirmed RPS cases with different stripping intensities, as well as undisturbed and gravitationally interacting galaxies. We find that extreme cases of RPS show SSD values ~3.5x higher than undisturbed and gravitationally interacting galaxies (56(+24/-15) as compared to 16(+6/-2) and 16(+6/-3), respectively), confirming simulation predictions. This enhancement reflects RPS-induced asymmetries: youngest stars are compressed along the leading edge and/or displaced into the extended tails of cold gas, while older populations remain undisturbed. In contrast, gravitational interactions perturb all stars uniformly, producing lower SSD values. SSD robustly distinguishes strong RPS cases, even adopting different age bins. This holds even without correcting for disk inclination, or when single-band imaging are used to trace stellar distributions. This makes SSD a promising tool to select RPS candidates for spectroscopic follow-up in upcoming surveys.

We present details regarding the construction of a composite spectrum of quasar (QSO) absorption line systems. In this composite spectrum we identify more than 70 absorption lines, and observe oxygen and hydrogen emission features at a higher signal-to-noise ratio than in any previous study. As the light from a distant quasar travels towards an observer, it may interact with the circumgalactic medium environment of an intervening galaxy, forming absorption lines. In order to maximize the signal of these absorption lines, we have selected a sample of 238,838 quasar spectra from the second data release of the Dark Energy Spectroscopic Instrument (DESI), each identified to have absorption lines resulting from such an interaction. By stacking these spectra in the restframe of the absorption, and calculating a median composite spectrum, we are able to isolate and enhance these absorption lines. We provide a full atlas of all detected absorption and emission lines as well as their fit centroids and equivalent width values. This atlas should aid in future studies investigating the compositions and physical conditions of these absorbers.

We present new MIRI F560W, F770W and F1000W imaging of GN-z11, extending the previous rest-frame coverage from 0.38 to 0.86$\mu$m. We report significant detections (14$\sigma$) in the F560W and F770W images, and a marginal detection (3.2$\sigma$) in F1000W. Here, we analyse its SED combining new MIRI imaging data with archival NIRSpec/Prism and MRS spectroscopy, and NIRCam imaging. The continuum emission shows a flat energy distribution, in f$_\nu$, up to 0.5$\mu$m, compatible with the presence of a mixed stellar population of young (4$\pm$1 Myr) and mature (63$\pm$23 Myr) stars that also account for the [O III], H$\beta$ and H$\alpha$ emission lines. The continuum at rest-frame 0.66$\mu$m shows a 36$\pm$3% flux excess above the predicted flux for a mixed stellar population, pointing to the presence of an additional source contributing at these wavelengths. This excess increases to 91$\pm$28% at rest-frame 0.86$\mu$m, although with a large uncertainty due to the marginal detection in the F1000W filter. We consider that hot dust emission in the dusty torus around a type 2 AGN could be responsible for the observed excess. Alternatively, this excess could be due to hot dust emission or to a photoluminiscence dust process (Extended Red Emission, ERE) under the extreme UV radiation field, as observed in some local metal-poor galaxies and in young compact starbursts. The presence of a type 1 AGN is not supported by the observed SED since high-z QSOs contribute at wavelengths above rest-frame 1$\mu$m, and an additional ad-hoc red source would be required to explain the observed flux excess at 0.66 and 0.86$\mu$m. Additional deep MIRI imaging covering the rest-frame near-IR are needed to confirm the flux detection at 10$\mu$m with higher significance, and to discriminate between the different hot dust emission in the extreme starburst and AGN scenarios with MIRI imaging at longer wavelengths.

We investigate the LiF spectrum up to $7800$ cm$^{-1}$ above the first dissociation limit. The ab initio calculations of the adiabatic potentials and other molecular functions are performed in a wide range of interatomic separations, $r=1$-17 bohr. We consider the model of two interacting electronic states including both the bound states and the resonances of two kinds, the tunneling resonances and the predissociative ones. The Born-Oppenheimer potentials are modeled with use of two functions imitating the diabatic terms whereas the coupling between them was set constant equal to the half of the minimum separation between the adiabatic terms, and then we define the final diabatic terms and the final diabatic coupling via the adiabatic potentials and the angle of the adiabatic-to-diabatic basis rotation obtained by integration of the nonadiabatic coupling matrix element. The energies of the bound states, as well as the positions and widths of the resonances are calculated. The observed transition frequencies are reproduced with the standard deviation of 0.0009 cm$^{-1}$ for $^7$LiF, 0.0006 cm$^{-1}$ for $^6$LiF, and within the experimental uncertainties for the most of the lines. The line lists for the bound-bound $X$-$X$ rovibrational transitions are calculated for quantum numbers $v\le50,\Delta v\le15,J\le170$ ($J\le200$ for the 0-0 and 1-0 bands).

Physics beyond the Standard Model may give rise to first-order phase transitions proceeding via the nucleation of vacuum bubbles, whose subsequent collisions generate gravitational waves (GWs). Their detection would open the possibility of investigating the universe in its first instants. If the transition is slow enough, such that it completes with the nucleation and collision of only two bubbles, the resulting GW signal is linearly polarized. We show that in this case triangular interferometers such as LISA and the Einstein Telescope could be able to not only measure the magnitude of the GW but also establish its linear polarization. This would give a strong hint about the origin of the signal.

We study a holographic D3/D7 system, whose dilaton profile has been phenomenologically adjusted in the infrared. The model is used to describe a deconfined yet massive quark phase of QCD at finite density, concluding that the equation of state of such a phase can be stiff enough to support exotic dense stars as massive as 2 solar masses. Nucleons are modeled phenomenologically using the this http URL EFT baryon phases. For the stiff phenomenological baryon phases the transition to the quark phase is weakly first order allowing for stable quark cores. We also find that holographic baryons, modeled as wrapped D5-branes, provide unrealistic pressures (in the homogeneous approximation) and have to be discarded. We compute the mass vs. radius relation and tidal deformability for these hybrid stars. Contrary to a large number of other holographic models, this holographic model indicates that quark matter could be present at the core of heavy compact stars and may be used to explore the phenomenology of such objects.

The hot quark matter created in heavy-ion collision experiments can exhibit strong temperature and chemical-potential gradients, which in turn can generate electric fields through thermoelectric effects. In this work, we investigate two relevant thermoelectric coefficients - the thermopower (Seebeck coefficient) and the Thomson coefficient - of two-flavor quark matter using the Kubo formalism and the Nambu--Jona-Lasinio model as an effective description of dense, finite-temperature QCD. The required two-point equilibrium correlation functions are evaluated using the Matsubara formalism of thermal field theory, applying a 1/Nc expansion to the relevant multi-loop Feynman diagrams. We employ previously derived quark spectral functions obtained from one--meson-exchange diagrams above the Mott transition temperature. Our numerical results show that both thermoelectric coefficients increase approximately linearly with temperature and decrease with increasing chemical potential. We also estimate the magnitude of the electric fields that can be generated in heavy-ion collisions by thermal gradients via the Seebeck effect.

We explore extensions of type I seesaw framework with a scalar mediator ($\Phi$) connecting to a complex scalar dark field ($S$), and right handed neutrinos ($N_i$), with an aim to correlate neutrino mass generation, leptogenesis, and dark matter. $\mathcal{Z}_4\times CP$ turns out to be a phenomenologically viable choice of the extended symmetry, which can accommodate a dimension five effective interaction $\bar{l}_L^\alpha \tilde{H}\Phi N_i$, involving the SM lepton isodoublet ${l}_L$, and Higgs $H$; prohibiting the canonical Yukawa term $\bar{l}_L^\alpha \tilde{H} N_i$. The $\mathcal{Z}_{4}$ symmetry is spontaneously broken via the vacuum expectation value (VEV) of the $\Phi$ filed, which directly affects neutrino mass generation and leptogenesis; while the $CP$ symmetry stabilises one component of $S$, making it a viable dark matter candidate. The discrete symmetry breaking creates domain wall, which needs to be annihilated before the over-closure of the Universe. This paves the way for gravitational wave signal associated with the model set up, which probes the symmetry breaking scale, and indirectly connects to the other phenomena.

We investigate the axial electromagnetic quasinormal modes of a static, asymptotically Anti--de Sitter (AdS) black hole sourced by a nonlinear electrodynamics model of Plebański type. Starting from the master equation governing axial perturbations, we impose ingoing boundary conditions at the event horizon and normalizable (Dirichlet) behavior at the AdS boundary. Following the approach of Jansen, we recast the radial equation into a linear generalized eigenvalue problem by using an ingoing Eddington--Finkelstein formulation, compactifying the radial domain, and regularizing the asymptotic coefficients. The resulting problem is solved using a Chebyshev--Lobatto pseudospectral discretization. We compute the fundamental quasinormal mode frequencies for both the purely electric ($Q_m=0$) and purely magnetic ($Q_e=0$) sectors, emphasizing the role of the nonlinearity parameter $\beta$ and the effective charge magnitude $Q$. Our results show that increasing either $\beta$ or $Q$ raises both the oscillation frequency $\omega_R$ and the damping rate $-\omega_I$, leading to faster but more rapidly decaying ringdown profiles. Nonlinear electrodynamics breaks the isospectrality between electric and magnetic configurations: magnetic modes are systematically less oscillatory and more weakly damped than their electric counterparts. For sufficiently large $\beta$ and small $Q_m$, the fundamental mode becomes purely imaginary ($\omega_R \approx 0$), in agreement with the absence of a trapping potential barrier in this regime. These findings reveal qualitative signatures of nonlinear electromagnetic effects on black hole perturbations and may have implications for high-field or high-charge astrophysical environments.

In the framework of $F(\mathcal{R},\tilde{\mathcal{R}})$ Einstein-Cartan gravity with an action depending both of the Ricci scalar and the so-called Holst-invariant curvature we consider models that include cubic terms of the latter in the action and study their inflationary behavior. These terms can have a considerable effect either positive or negative in relation to the agreement with present observational data, depending on parameters. In parameter regions where the quadratic models fail to produce results consistent with observational data, the presence of these additional cubic terms can lead to compatible predictions.

Modified Newtonian Dynamics (MOND) is a framework of theories that adjust Newton's laws of gravity to explain effects such as galactic rotation anomalies, offering an alternative to dark matter. This essay examines the justification of MOND by assessing its inter-theoretical relationship to established theories across relevant scales, in particular its connection to Newtonian gravitation. We argue that MOND fails a key condition for a theory's justification--what we call 'reduction-wise justification'--since it does not adequately reduce to Newtonian gravity in a fully non-arbitrary way. More precisely, despite satisfying the standard formal criteria for successful limiting reduction, MOND does not properly reduce to Newtonian gravitation because of (i) the absence of a fundamental theoretical framework to justify the interpolating function introduced in MOND and (ii) the lack of a unified mathematical structure working across all scales, independent of Newtonian theory. Hence, the case study of MOND provides crucial results for the general debate on inter-theoretic reduction in science: MOND's failure as a case of reduction highlights important limitations in standard accounts of limiting reduction. We respond by proposing a more refined framework for limiting reduction that introduces two additional criteria to better distinguish successful from pathological reductions. More broadly, this case illustrates how analysing reduction-wise justification can serve as a powerful tool for evaluating the validity of novel theories that are not yet empirically established.

Gravitational-wave data analysis relies on accurate and efficient methods to extract physical information from noisy detector signals, yet the increasing rate and complexity of observations represent a growing challenge. Deep learning provides a powerful alternative to traditional inference, but existing neural models typically lack the flexibility to handle variations in data analysis settings. Such variations accommodate imperfect observations or are required for specialized tests, and could include changes in detector configurations, overall frequency ranges, or localized cuts. We introduce a flexible transformer-based architecture paired with a training strategy that enables adaptation to diverse analysis settings at inference time. Applied to parameter estimation, we demonstrate that a single flexible model -- called Dingo-T1 -- can (i) analyze 48 gravitational-wave events from the third LIGO-Virgo-KAGRA Observing Run under a wide range of analysis configurations, (ii) enable systematic studies of how detector and frequency configurations impact inferred posteriors, and (iii) perform inspiral-merger-ringdown consistency tests probing general relativity. Dingo-T1 also improves median sample efficiency on real events from a baseline of 1.4% to 4.2%. Our approach thus demonstrates flexible and scalable inference with a principled framework for handling missing or incomplete data -- key capabilities for current and next-generation observatories.

$\alpha$-attractor models naturally appear in supergravity with hyperbolic geometry. The simplest versions of $\alpha$-attractors, T- and E-models, originate from theories with non-singular potentials. In canonical variables, these potentials have a plateau that is approached exponentially fast at large values of the inflaton field $\varphi$. In a closely related class of polynomial $\alpha$-attractors, or P-models, the potential is not singular, but its derivative is singular at the boundary. The resulting inflaton potential also has a plateau, but it is approached polynomially. In this paper, we will consider a more general class of potentials, which can be singular at the boundary of the moduli space, S-models. These potentials may have a short plateau, after which the potential may grow polynomially or exponentially at large values of the inflaton field. We will show that this class of models may provide a simple solution to the initial conditions problem for $\alpha$-attractors and may account for a very broad range of possible values of $n_{s}$ matching the recent ACT, SPT, and DESI data.

In this work we study the possibility of modeling the dark matter content in galaxies as a core-halo model consisting of self-gravitating, self-interacting fermions. For the core of the halo, the dark matter fermions are degenerate, while for the halo we have considered two possibilities: the fermions have thermalized as a perfect fluidor they will follow a standard cold dark matter Navarro-Frenk-White profile. The core density profile is obtained by solving the Tolman-Oppenheimer-Volkoff equations, and their properties are determined by the fermion mass, the central density and the interaction strength. The mass of the fermion and the strength of the fermion self-interaction is fixed by doing a $\chi^2$ analysis to fit that fit the rotational curves of Low Surface Brightness galaxies. It was found that the fermion mass should be in the range $38.73~\rm{eV}< m_{f} < 42.11~\rm{eV}$ and the interparticle strength in the range $269.69 < y <348.48$ at $68$ C.L. in order to reproduce the rotational curves adequately, in the case when the halo is modeled as a thermalized ideal gas. Similar values are obtained if the halo is modeled following a Navarro-Frenk-White case, namely $41.54 ~\rm{eV} < m_{f} <49.87 ~\rm{eV}$ and $5606.06< y < 17484.84$. Once fixed the values of the mass of the fermion $m_f$ and the interaction strength $y$, we tested the core-halo model with data from the Milky Way and the SPARC database. We have found good agreement between the data and the core-halo models, varying only one free parameter: the central density. Thus a single fermion can fit hundreds of galaxies. Nevertheless, the dark matter halo surface density relation or the halo total mass and radius depend strongly on the model for the halo.

We study the structural and dynamical consequences of introducing a distributional density profile inside a neutron star, representing a massless, topological shell located at an arbitrary radius. We incorporate this effect into the structure of neutron star and construct equilibrium sequence for several realistic equations of state. Radial stability is examined through the Sturm-Liouville formulation of the $\ell=0$ perturbation equation, supplemented with a jump condition and imprinting distinct features on the fundamental $f$-mode spectrum. We find strong, non-monotonic variations in the mode frequency relative to standard no-shell models. Using first-principles scaling relations, we estimate various gravitational wave observables such as the damping time, quality factor, luminosity and characteristic strain. These observables are then compared with the sensitivity of Advanced LIGO, and third-generation detectors such as the Einstein Telescope and Cosmic Explorer. Our results demonstrate that internal topological shells can leave potentially observable signatures in the oscillation and gravitational wave properties of neutron stars.

We investigate the challenges of constructing neutron star (NS) solutions in $f(\mathbb{Q})$ gravity, highlighting the importance of treating the affine connection as an active, dynamical component of the theory. We begin by clarifying under what conditions standard simplifications -such as the coincident gauge or General Relativity (GR)-like connections- inadvertently lead to GR behavior, even in non-trivial $f(\mathbb{Q})$ models. Building on previous work in black hole (BH) spacetimes, we adapt the formalism to NS and extend it to non-vacuum configurations. Focusing on two representative models, $f(\mathbb{Q}) = \mathbb{Q} + \alpha \mathbb{Q}^2$ and $f(\mathbb{Q}) = \mathbb{Q}^\beta$, our analysis suggests that, under standard regularity assumptions, solutions with Maclaurin/Laurent-type series recover GR dynamics, pointing to more intricate structures as the likely seat of beyond-GR effects, and reflecting the constraints imposed by the connection's dynamics on the asymptotic behavior of genuinely beyond-GR solutions. We then formulate the problem as a boundary value problem (BVP) and highlight the numerical pathologies that may arise, together with possible strategies to prevent them. This work aims to provide a concrete framework for future numerical studies and outlines the theoretical consistency conditions required to construct physically meaningful beyond-GR NS solutions in $f(\mathbb{Q})$ gravity.