高能暂现源 · 每日文献综述

The magnetic Penrose process (MPP) is one of the most efficient mechanisms for extracting rotational energy from a magnetized Kerr black hole (BH), enabling charged particles to reach very-high energies. We investigate this process in Sagittarius A* (Sgr A*), the supermassive BH at the Galactic center (GC), focusing on its ultra-high-efficiency regime, in which neutrons undergo beta decay inside the BH ergosphere. As a key novelty, we compute the neutron production spectrum in the accretion flow directly from the nuclear reaction kinematics and follow neutron trajectories in the Kerr spacetime to determine the population that reaches the ergosphere and undergoes the MPP. From this population, we derive the spectrum of accelerated protons and show that Sgr A* can accelerate them up to PeV energies, strengthening its interpretation as a candidate Galactic PeVatron. We further compute the gamma-ray and neutrino emission produced through hadronic interactions of the escaping protons in the Central Molecular Zone. The predicted gamma-ray fluxes exhibit distinctive spectral features that could provide an observational signature of the MPP and may constitute a non-negligible contribution to the very-high-energy emission detected by H.E.S.S. and HAWC. The associated neutrino fluxes remain below the diffuse Galactic component inferred by IceCube, but may still contribute to the high-energy emission from the GC. Remarkably, the predicted signals lie within the projected sensitivity of SWGO for all scenarios considered here and, for some models, only a factor of a few below the nominal CTAO sensitivity, while KM3NeT/ARCA and IceCube-Gen2 will provide complementary tests. Our results establish the MPP as an observable mechanism for extracting BH rotational energy, providing a direct connection between horizon-scale physics and multimessenger observations that can be extended to other magnetized BHs.

We present new observations of the metal-poor ($< 0.1$ $\rm Z_\odot$) Type II supernova 2023ufx and its host galaxy. The deep nebular spectrum, obtained roughly a year after explosion, has a triple-peaked [O I] emission profile suggesting an asymmetric explosion. Comparison to nebular spectra models suggests a zero-age main sequence mass of $ M_{\rm ZAMS} \sim 25-35$ $\rm M_\odot$, which is supported by the low [Ca II]/[O I] emission-line ratio of $\approx0.4$. The diminishing, broad, boxy H$α$ emission and flattening of the light curve in the late ($\sim2$ yr) photometry suggest a complex mass-loss history in the centuries to millennia before explosion. New optical and near-infrared imaging of the host galaxy confirms that it is a dwarf, with a stellar mass of $10^{6.6\pm0.1}$ $\rm M_\odot$ and a SFR of $10^{-2.5\pm0.1}$ $\rm M_\odot$/yr. Results from both SED fitting and galaxy stellar mass-metallicity scaling relations all lead to an environmental metallicity estimate of $0.02-0.07$ $\rm Z_\odot$. Taken together, these observations confirm that SN 2023ufx is the explosion of a very metal-poor, heavily stripped red supergiant.

The exceptional gravitational-wave event GW231123 exhibited unusually high total mass and near-extremal spins, while revealing strong waveform-dependent variations in its inferred source properties. Several follow-up studies favored a lensed interpretation, with reduced waveform systematics, raising the question of whether apparent lensing evidence could arise from waveform modeling systematics. Using GW231123-like numerical-relativity injections, we find that lensing degrees of freedom can reduce waveform discrepancies and yield non-negligible lensing support. However, this behavior is not systematic across injections, and cases where all waveform and lens models favor lensing while reducing waveform discrepancies are not found.

We present the Active Galactic Nuclei Database (AGN-DB), a comprehensive, multi-wavelength catalog compiled from more than 100 publicly available AGN catalogs and samples released by the end of 2025, spanning radio to $γ$-ray wavelengths. The database contains approximately 8.1 million unique sources, approximately 7.8 million of which remain after flagging stellar contaminants, and approximately 6.8 million of these are classified as AGN. Source cross-matching across catalogs is performed using Lyra, a Bayesian likelihood-ratio framework that jointly considers positional uncertainties, source densities, and photometric information to compute posterior match probabilities. The resulting catalog provides astrometric coordinates, redshifts, photometry, and classifications for each unique source. All multi-catalog provenance is preserved. For every property, we store the full array of values and originating catalog identifiers, enabling multi-epoch and multi-survey analyses. In this paper, we describe the AGN-DB pipeline, including the cross-matching methodology, and present the statistical properties of the v1.0 catalog. AGN-DB is designed to enable population studies, spectral energy distribution modeling, AGN classification, and variability analyses at an unprecedented scale. Its pipeline is designed to facilitate the integration of new catalogs, allowing AGN-DB to be updated regularly, with releases planned at least annually.

The non-extensive entropy models applied to the black-hole horizon are connected to the generalized Einstein--Hilbert action, $f(R)$, via the Wald entropy formalism. We demonstrate that quark star configurations characterized by the MIT bag model under this modified gravity approach conform to the observational lower and upper bounds established for compact objects linked to HESS J1731$-$347 and GW190814. Furthermore, we assess the effective energy conditions and corresponding speed of sound to evaluate the physical plausibility and stability of the stellar configurations. All non-extensive entropy models studied satisfied these conditions. Our findings show that the exotic compact objects offer a compelling and credible framework for exploring aspects of non-extensive entropy models, such as the Barrow, Tsallis-Cirto, and Rényi formalisms, and vice versa.

Long-duration gravitational-wave modelling must resolve fast orbital motion together with slow dissipative evolution while preventing small numerical errors from accumulating into secular phase drift. Here we ask whether the finite-time evolution map itself can be learned as an explicit, differentiable, structure-preserving object and then repeatedly composed through a complete inspiral. We construct three neural-flow architectures: a symplectic and slimplectic flow on Galley's doubled phase space, [SINFONIA-J0]; a Taylor-anchored flow, [SINFONIA-J1]; and a Magnusian flow that learns the finite-time dissipative correction in the interaction picture, [SINFONIA-J2]. Applied to a 2.5PN neutron-star inspiral, all three expose the same controlling mechanism: long-time accuracy is governed not by pointwise map error alone, but by its signed projection onto a single secular channel fixed by energy--angular-momentum balance. Encoding this structure allows the learned maps to remain accurate through $10^{2}$--$10^{5}$ window compositions to coalescence at timesteps of a full orbital period and beyond, reaching chained phase errors orders of magnitude below a benchmark slimplectic integrator at lower cost. The same secular structure can also be exploited for physics inference: when the channel is left unconstrained, the accumulated phase retains enough information to recover an un-modelled dynamical-friction-like force, both parametrically and as a learned function of separation. Network-off controls isolate the contribution of learning from the analytic structure already built into each map. These results establish a proof of concept for structure-preserving learned evolution maps as tools for fast long-duration integration and physics inference in gravitational-wave source modelling.

A new condition involving the complexity factor, a structure scalar arising from the orthogonal splitting of the Riemann tensor, is generated at the boundary of a spherically symmetric anisotropic fluid. This condition leads to the generation of a model, without invoking any further constraints on geometry or matter variables. The model generated is physically reasonable and suggests that there is a non-zero, positive minimum constraint on the complexity scalar value evaluated at the surface. This suggests that vanishing complexity is a globalized constraint which does not apply at the boundary, specifically for compact objects in the strong-gravity regime.

Massive-field perturbations on black-hole and wormhole backgrounds can form long-lived quasinormal ringing; this behaviour is not confined to brane-world spacetimes, but it is also not universal across compact-object geometries. Recent work on a brane-world black hole reported an unusual limiting behaviour in which the real oscillation frequency tends to zero as the scalar mass grows. We test how general this phenomenon is in the Bronnikov--Kim effective brane-world geometries, treating both the black-hole and wormhole sides of the same families. The comparison uses two complementary numerical treatments: Leaver's continued-fraction method for the black-hole side, after reducing the wave equation for the second Bronnikov--Kim family to rational form, and Chebyshev pseudospectral collocation for regular wormhole throats. The black hole of the second Bronnikov--Kim family shows the standard quasi-resonant pattern: the damping rates of the scalar branches considered here become very small, but the real frequencies remain finite. Representative wormholes from the first and second Bronnikov--Kim families also develop increasingly long-lived massive scalar modes, and extending the second-family wormhole scan to larger field mass reduces the damping by more than a factor of four relative to the massless case. Thus the robust brane-world effect found here is the appearance of long-lived massive-field ringing in both black holes and wormholes, whereas the vanishing of the real frequency is metric- or branch-dependent.

The role of large-scale pre-existing interplanetary structures, including preceding coronal mass ejections (CMEs) and stream interaction regions (SIRs), in shaping the shock acceleration environment for energetic electrons remains not fully understood. In this study, we investigate nine interplanetary shocks observed by the Solar Terrestrial Relations Observatory (STEREO) that are associated with significant MeV electron enhancements, as such enhancements are rarely observed at interplanetary shocks. We combine remote-sensing observations, drag-based modeling, and in-situ measurements to analyze the shock propagation through pre-existing interplanetary structures. Eight of the nine events are associated with a preceding slow or intermediate-speed CME, while six shocks are in-situ observed propagating within preceding ICMEs, indicating that large-scale upstream trapping structures are a common feature of these events. Further analysis identifies three distinct scenarios associated with enhanced electron acceleration: shocks propagating through preceding ICMEs, shock-SIR interactions, and direct injection of flare-accelerated electrons into SIRs. As a representative shock-in-ICME event, the 2012 January 29 low-$β$, quasi-perpendicular shock ($θ_{Bn}\sim87^\circ$) is further investigated using observations together with one-dimensional Monte Carlo test-particle simulations of a shock propagating into a large-scale upstream magnetic loop. The simulation suggests that the upstream loop prolongs electron residence near the shock and substantially enhances acceleration efficiency. These results demonstrate that large-scale interplanetary structures can precondition the upstream magnetic environment, providing favorable conditions for prolonged electron residence and efficient shock acceleration.

Tidal perturbations from nearby compact objects can drive resonant "kicks" in extreme-mass-ratio inspirals (EMRIs), imprinting potentially detectable phase shifts in LISA-band gravitational waveforms. We provide a systematic survey of these tidal resonances across EMRI parameter space in the weak-tide, three-body hierarchy, treating the perturber as stationary on the resonance timescale. For generic Kerr orbits we compute (i) resonance contours in $(p,e,x)$ that locate where each $(n,k,m)$ resonance is encountered, (ii) the associated resonance duration, and (iii) the jump amplitudes in the angular momentum $L_z$ and the Carter constant $Q$ (with $ΔE$ vanishing by time-translation symmetry in the stationary model). Combining jump amplitudes with resonance durations, we construct a practical ranking of the resonances most relevant for waveform modeling, while emphasizing that even modest resonances may still be important through their influence on the phase at which subsequent resonances are entered. All contour and jump data are publicly available.

The properties of compact stars with a strong first-order phase transition to quark matter and with an additional fluid of self-interacting bosonic dark matter (DM) are studied. We find that the inclusion of DM changes considerably the stability of mass-radius configurations relative to the naive one-fluid criterion. For compact star configurations with similar masses and different radii, so-called twin stars, the presence of DM removes the unstable segment between the hadronic and the hybrid branch, so that the stable mass-radius sequence becomes continuous after the onset of the phase transition to quark matter. We furthermore find stable ultra-compact objects (UCOs), defined by a total compactness $C = M_\text{tot}/R_\text{grav} \ge 1/3$. We observe two distinct classes of UCOs: a DM-halo class with $f_\text{DM} \gtrsim 0.9$, and a DM-core class at $f_\text{DM} \lesssim 0.02$. The two classes can be separated by the surface redshift of the normal matter, which reaches $z=0.73$--$0.77$ for the DM-core class and stays below $0.45$ for the DM-halo class. Finally, we find hybrid star solutions of 'ultimate twins' with similar mass and visible radius, but different dark matter content, leading to different tidal deformabilities and surface redshifts. Future X-ray and gravitational measurements of ultra-compact neutron stars with radii and masses outside the allowed neutron star range can thereby probe the presence and the properties of DM in addition to a first-order phase transition to quark matter.

A gravitational-wave (GW) detector records a linear mixture of local couplings under the assumption that extra polarizations enter geodesic deviation. Vacuum general relativity (GR) admits two transverse-traceless (TT) amplitudes. It remains to determine the largest set of couplings that can sit in that mixture, and whether a tensor-only arm-length result selects GR uniquely. The little group \(E(2)\) of a null four-momentum classifies the strain amplitudes \(\mathbf{p}=(p_{+},p_{\times},p_{x},p_{y},p_{b},p_{\ell})\), commonly written \(h_{P}\), which determine the electric tidal tensor along a ray. Geodesic deviation, recorded as differential arm length, therefore contains only those \(p_{P}\) that enter that tensor. Lorentz mixing at helicity \(\pm1\) supplies a gravito-magnetic (GEM) field that does not enter the tidal tensor. if GEM field is static, that does not propagate as a wave. A time-varying helicity-\(\pm1\) current sources a GEM wave that enters the mixture as a coupling, to be isolated by its measured quantity. For a radiation-zone wave that depends only on retarded time, \(\betag_{\perp}=\hat{\mathbf{k}}\times\partial_{t}(p_{x},p_{y})\), and the eight couplings are \(\Pvec=(\mathbf{p},\betag_{\perp})\). Here we adopt unified field equations on \(\mathbf{p}\) to clarify the origin of each component of \(\Pvec\); the measured quantity of each coupling then isolates the polarizations in that mixture, which favors identification of distinct polarizations and model tests.

In the course of selecting targets for a spectroscopic survey of faint planetary-nebula nuclei, we encountered two anomalous cases of late-type variable stars surrounded by nebulosities discovered by amateur astronomers. We obtained follow-up deep narrow-band images of the two nebulae: StDr Objet 14 (diameter 1.5 arcmin) and StDr 4 (diameter 2.5 arcmin). StDr Objet 14 is found to be a low-excitation nebula, bright only in Ha+[N II], and showing a unique double-ring "bull's-eye" morphology. StDr 4 exhibits a bow-shock structure, bright in Halpha and [S II], embedded in an extended emission nebula. Both central stars, respectively GR Cygni and [D75] 141, are found to be pulsating semiregular variable stars. Our spectroscopic observations show GR Cyg to be a cool carbon star of type C9,2, with emission lines of the Balmer series, [O II] and [O III], [Ne III], [S II], and Fe II. A synthetic narrow-band image of GR Cyg shows that the source of its strong [O II] 3727 A emission line is a knot displaced from the star by about 1.6 arcsec, indicative of a recent collimated ejection event. [D75] 141 is a M6 III red giant, with weak Balmer-series emission. Both stars show a blue excess, and GR Cyg was detected in the near-UV by the GALEX spacecraft. These data establish both stars as new symbiotic binaries, containing a late-type red giant and a hot companion accreting from its stellar wind. We speculate that the double rings surrounding GR Cyg were ejected during helium-shell thermal pulses (or conceivably from recurring nova outbursts), while the nebula around [D75] 141 may be the result of a collision of the symbiotic star's ejecta with an interstellar cloud. Our accidental discoveries of these two binaries because they are surrounded by faint nebulae supports the view that there is a significant population of "missing" symbiotic stars yet to be discovered.

We propose using supernova microlensing to search for ultracompact minihalos and other extended dark structures seeded by enhanced primordial small-scale power. Unlike conventional microlensing searches, which are often optimized for primordial black holes, cosmological supernovae are sensitive to lenses with physical sizes comparable to their Einstein radii (up to $\sim {\rm pc}\left(M/10^4M_\odot\right)^{1/2}$ scale), opening a qualitatively different region of dark object parameter space. We compute the extended object detection efficiency relative to point lenses. We then show how projected limits on the abundance of ultracompact minihalos can be mapped, in benchmark formation scenarios, onto constraints on the primordial power spectrum. With the Vera C. Rubin Observatory expected to deliver an unprecedented LSST supernova sample, supernova microlensing is poised to become a timely and complementary probe of small-scale structure and the physics of the early Universe.

We update the X-ray information given in Maselli et al. (2024) for five bright radio sources in the SMS4 catalog, thanks to additional observations with the X-Ray Telescope (XRT) on board the Neil Gehrels Swift Observatory (hereafter Swift) carried out through February 2026. X-ray emission from MRC B1754-597, previously based only on data from the eROSITA-DE DR1 catalogs, is now detected also by Swift, which yields a more precise positional uncertainty. Thanks to ~7 ks of additional exposure with Swift we are now able to detect MRC B1817-391, a source for which no X-ray counterpart is found in eROSITA-DE DR1. Another source that is not detected in the eROSITA-DE DR1 catalogs, and that was not observed by Swift earlier than 2025, is MRC B1827-360: with a ~6 ks X-ray observation, the source is detected. Finally, based on ~8 ks of additional exposure, we detect X-ray emission for MRC B2032-350, a source that lies out of the DR1 footprint. In contrast, additional exposure (~2 ks) just gives an upper limit on the X-ray emission of MRC B1814-519. The analysis of the extent and the hardness ratio of the four detected X-ray sources suggests the presence of soft, diffuse X-ray emission as expected from galaxy group cores or hot galaxy coronae. Using the positional uncertainty of the X-ray detections to constrain the search for counterparts at lower energies, we provide a new infrared/optical counterpart for MRC B1817-391 and confirm the counterparts previously reported in the literature for the three remaining sources.

Narrow-line Seyfert 1 (NLSy1) galaxies are key laboratories for studying rapid supermassive black hole (SMBH) growth and active galactic nucleus (AGN) feedback at high accretion rates. We investigate the nearby NLSy1 Ton S180 with VLT-MUSE optical integral field spectroscopy to connect its nuclear accretion properties with the spatially resolved ionized gas and host-galaxy kinematics. We modeled the unresolved nuclear spectrum and applied a custom point spread function subtraction to recover the host-galaxy emission on kiloparsec scales. The nuclear spectrum requires a complex permitted-line decomposition and a blueshifted [O III] outflow component. Single-epoch estimators and the stellar velocity dispersion imply black hole masses in the range $\log(M_{\rm BH}/M_\odot) = 6.5 - 7.7$. Combined with the observed luminosity, this implies a dimensionless mass accretion rate of $\dot{M} / \dot{M}_{\rm Edd} = 4.1 - 980$, confirming the extreme accretion regime. The host galaxy shows a circumnuclear ring, an inner elongated structure consistent with a bar, and rotation-dominated gas and stellar kinematics. Simple inflow models do not significantly better reproduce the observed velocity field. We detect a resolved ionized outflow extending about 2 kpc west of the nucleus, with mildly blueshifted velocities (with a maximum of $v_{\rm max} \sim 340$ km s$^{-1}$). Its mass outflow rate is only $\sim 4.2 \times 10^{-4}\,M_\odot\,\mathrm{yr}^{-1}$, whereas the unresolved nuclear outflow reaches $v_{\rm max} \sim 1140$ km s$^{-1}$ and $\dot{M}_{\rm out} > 1.5\,M_\odot\,\mathrm{yr}^{-1}$. This contrast may reflect either weak ionized coupling from nuclear to galactic scales or different episodes of AGN activity over time. These results show that Ton S180 is undergoing rapid SMBH growth, while the observed ionized outflow remains confined to the inner few kiloparsecs and is weak on host-galaxy scales.

We present a phenomenological frequency-domain model for the gravitational-wave (GW) memory signal from nonspinning binary-black-hole mergers on quasicircular orbits. We develop separate amplitude and phase models for the dominant $(l,m)=(2,0)$ spherical-harmonic mode of the GW memory signal. The amplitude and phase models are built from superpositions of elementary and transcendental functions, which can be evaluated efficiently. Both portions of the model are calibrated against a numerical-relativity surrogate model over mass ratios from one to eight. Their accuracy is assessed by computing their mismatch with the numerical-relativity models using the Advanced LIGO sensitivity curve from the fourth observing run. The mismatches are of the order $10^{-4}\unicode{x2013}10^{-3}$ over the parameter space of total mass covered by LIGO. The resulting frequency-domain model is a more computationally efficient waveform model for the GW memory signal than a related time-domain model earlier produced by the authors. An open-access implementation of both the time-domain and frequency-domain models is provided in the Python package GWMemoryModel, which can be used for relevant analyses of nonspinning binary black holes.

Quadratic $f(R)$ gravity carries a massive scalar degree of freedom, the scalaron, whose finite range $λ$ screens its influence on the geometry outside a compact object. We show that this screening severs the exterior of a black hole from that of a neutron star of the same mass. The correction to the Schwarzschild metric is Yukawa-suppressed, falling as $e^{-r/λ}$ instead of polynomially in the coupling, and is thus invisible to any expansion in powers of that coupling; also, the exterior is intrinsically isotropic, so that inverting the temporal potential alone, as in single-potential solutions, fails to solve ]field equations. Applying the Gauss-Bonnet theorem to this geometry, we find the leading deflection angle to be exactly the general-relativistic $4GM/b$, the scalaron contributions cancelling identically in the combination that bends light. The first correction carries the unfamiliar signature $λ^{-1/2}b^{-3/2}e^{-b/λ}$, screened beyond the scalaron range and confirmed against direct quadrature to around 20% at $b=4λ$, improving to 7% at $b=12λ$. What survives is not a difference of degree but of kind. The three ingredients that deliver it, non-analyticity in the coupling, the intrinsically isotropic gauge and the pressure-weighted scalar charge, are here obtained within a single, self-consistent derivation for the first time. A static black hole carries no scalar hair, and lenses precisely as GR requires; a neutron star acquires a scalar charge weighted by the pressure supporting it against collapse, and does not. Weak lensing hence closes as a discriminant of the theory, whilst the horizon, as against a material surface, remains one in principle. The observational advantage lies not in bending angles at large $b$ but in the strong-field imaging of the photon sphere, and in the stellar interior, where the scalar charge is fixed by EoS.

We investigate the tidal forces experienced by massive particles in radial free fall in the quantum Oppenheimer--Snyder black hole, a loop quantum gravity correction to the Schwarzschild geometry that describes the exterior of a bouncing dust ball. Using an orthonormal tetrad adapted to a freely falling observer, we show that the radial and angular components of the tidal tensor reproduce their Schwarzschild counterparts at large distances but reverse sign in the interior. The angular zero remains confined between the Cauchy and event horizons and coincides with the latter only in the extremal limit. A particle released from rest does not reach the center, it stops at a turning point inside the Cauchy horizon, where the ratio of the tidal components is independent of the parameters of the solution and is fixed by the transverse equation of state of the effective source. Solving the geodesic deviation equations for two sets of initial conditions, we find that the radial deviation vector attains its maximum at the minimum of the metric function, exactly for one set of initial conditions and asymptotically for the other, and remains finite throughout, in contrast with the Schwarzschild case, in which it diverges at the singularity. This regularity originates in the bounce rather than in a regular core, and therefore protects timelike radial infall while leaving radial null geodesics unaffected. Examining the full range of the quantum parameter, we find that its sign determines whether the tidal sector possesses any of this structure, while its magnitude determines only whether that structure is hidden: beyond the extremal value the geometry becomes horizonless, and the entire deformation history, bounce included, is exposed to distant observers.

The heliospheric current sheet (HCS) is the boundary between open magnetic fields of opposite polarity in the solar wind that connect back to their respective open flux regions in the corona. The HCS is often described as a disc-like sheet warped by the combined effects of the inclination of the magnetic equator on the Sun, solar rotation, and solar wind expansion. The HCS is an extension of the helmet streamer belt-the closed-flux region resulting from the global dipole component of the observed photospheric magnetic field. In the absence of a strong dipole, a dominant quadrupole can result in two separate HCSs. Here we report on a set of idealized, axisymmetric solar wind simulations designed to study the transition of the single-HCS corona and inner heliosphere into a multiple-HCS configuration, run with the WindPredict-AW PLUTO MHD code. We systematically vary the relative strengths of the dipole/quadrupole contributions to the global field and show this results in the formation of a large-scale, closed-flux pseudostreamer that interacts with the ambient solar wind. We find that when the dipole and quadrupole magnitudes are approximately equal, our model corona can oscillate between a one- and two-HCS configuration on the timescale of days-to-weeks. We characterize the plasma dynamics of the transition of a global-scale pseudostreamer into a localized, secondary helmet-streamer belt with its own HCS surrounding a newly-opened, opposite-polarity coronal hole. We discuss the implications of our results for solar cycle evolution and the dynamic coupling of the highly structured corona and inner heliosphere.

We use Gaia DR3 astrometry and TESS time-series photometry to study NGC 1901, a sparse open cluster in the TESS Southern continuous viewing zone, that has been neglected due to its proximity on the sky to the Large Magellanic Cloud. We refine its membership list and measure rotation periods for 32 high-confidence rotating stars. By applying the gyro-interp empirical model, we calculate a gyrochronal age of 830 $\pm$ 50 Myr. Because NGC 1901 possesses a sub-solar metallicity of [Fe/H]$\approx -0.11$, its excellent agreement with empirical models suggests that temperature-period relations mitigate metallicity effects for Sun-like stars. Our results place NGC 1901 in the gap between the benchmark open clusters Praesepe (670 Myr) and NGC 6811 (1 Gyr). We find that the rotation sequence exhibits significant spin-down stalling in late K dwarfs, supporting models of core--envelope recoupling. Furthermore, our findings provide constraints for the cluster's structural evolution, including its recently identified tidal tail and corona. This study establishes NGC 1901 as a new benchmark for calibrating the next generation of gyrochronal models.

The puzzling bimodality displayed by blazars on the broad-band spectral indices ($α_{\rm ro}-α_{\rm ox}$) plane is still an open question. To investigate its physical origin, we compiled a sample comprising 136 flat-spectrum radio quasars (FSRQs), 64 low-synchrotron-peaked BL Lac objects (LBLs), 39 intermediate-synchrotron-peaked BL Lac objects (IBLs), and 105 high-synchrotron-peaked BL Lac objects (HBLs). Our results show that (FSRQs+LBLs)$\rightarrow$IBLs$\rightarrow$HBLs follow a \reflectbox{$\angle$}-shaped evolutionary sequence on the $α_{\rm ro}-α_{\rm ox}$ plane, with the average 5\,GHz radio Doppler factor, jet power, and accretion Eddington ratio decreasing along this \reflectbox{$\angle$}-shaped track. Additionally, the distributions of $α_{\rm ox}$ exhibit only slight differences among blazar subclasses, whereas those of $α_{\rm ro}$ and $α_{\rm rx}$ differ significantly. Our results can be naturally explained if the distribution of blazars on the $α_{\rm ro}-α_{\rm ox}$ plane is regulated by the accretion rates and the Doppler beaming effect across different bands ($δ_i^{q_i}/δ_j^{q_j}$). However, $δ_i/δ_j$ plays a pivotal role in the distribution, while jet shape ($q=2(3)+α_ν$) and spectral indices ($α_ν$) play a minor role. Our results suggest that FSRQs+LBLs are strong-jet sources coupled with radiatively efficient accretion, and their jets are accelerated to a larger distance from the core and still maintain relativistic speeds down to the radio region; HBLs are weak-jet sources coupled with radiatively inefficient accretion, and their jets may begin to decelerate in the optical region and exhibit sub-relativistic speeds as they propagate to the radio region. IBLs are regarded as transitional objects between FSRQs+LBLs and HBLs.

We investigate the dynamics of horizon soft degrees of freedom associated with near-horizon BMS supertranslations in an evaporating Vaidya-Schwarzschild black hole spacetime. Considering the dynamic nature of the supertranslation parameter, we derive its effective action directly from the Einstein-Hilbert action. The resulting Goldstone sector is intrinsically coupled to the evolving black hole dynamics, with the mass function entering the evolution of the Goldstone modes while the Goldstone configuration contributes to the dynamical evolution of the black hole mass. Azimuthal periodicity and regularity at the poles select the physically admissible angular sector and yields a mass-dependent $selection\ rule$ on the azimuthal modes. Consequently, the spectrum of horizon-supported Goldstone modes evolves with the black hole mass; as the horizon shrinks during evaporation, the horizon is left with a progressively reduced set of lower-order modes. The evaporation thus induces a dynamical filtering of the horizon soft sector, arising intrinsically from the black hole mass dynamics. These results provide an effective gravitational framework linking near-horizon BMS symmetry, dynamical Goldstone modes, and black hole evaporation, suggesting a direct connection between macroscopic horizon evolution with the microscopic organization of the gravitational soft degrees of freedom residing on the horizon.

Black holes, with their strong gravitational fields, provide an important testing ground for theories of gravity beyond General Relativity. Among the many proposed alternatives, considerable recent work has focused on scalar-tensor theories in which a scalar field couples to higher-curvature terms. The black hole solutions that arise in such theories can differ significantly from the Schwarzschild and Kerr solutions of General Relativity. Characteristic properties of these black holes include instabilities, shadows, and gravitational wave spectra, which can be used to constrain the couplings of the underlying theories.

We characterise the accretion rate properties of radio-detected AGN by combining deep radio and X-ray observations of the Boötes field. We used deep international LOFAR telescope observations to identify 2840 radio-AGN across $0.3 < z \leq 2$, divided into samples that are complete in radio luminosity. We further split this sample into four different classes: radio-quiet AGN (RQ-AGN), high-excitation radio galaxies (HERGs), and low-excitation radio galaxies (LERGs) hosted by star-forming (SF-LERGs) and quiescent galaxies (Q-LERGs). Through performing X-ray stacking, we determined the average X-ray luminosities, $L_{\rm{X, 2-10\,keV}}$ and average specific X-ray luminosities, $λ_{sL_{X}}$ (X-ray luminosity scaled by the stellar mass; a proxy for the accretion rate). We studied how these X-ray properties depend on radio luminosity, stellar mass, and redshift for each of the four AGN classes. We found that the LERGs, regardless of their star-formation activity, show significantly lower $L_{\rm{X, 2-10\,keV}}$ and $λ_{sL_{X}}$ than both HERGs and RQ-AGN across all redshifts. The average X-ray luminosities for the HERGs, RQ-AGN, and SF-LERGs typically increase with redshift, which may be associated with the increased cold gas fractions at earlier times, resulting in more enhanced black hole accretion. The average X-ray luminosities show weak-to-no correlation with radio luminosity at a given redshift, suggesting that the physical processes producing the X-ray and radio emission may not be coupled on spatial and/or temporal scales.

The moving lens (ML) effect is a secondary anisotropy of the cosmic microwave background (CMB) generated by the transverse motion of gravitational potentials, providing a direct probe of the large-scale cosmic velocity field. Its detection relies on cross-correlating a CMB map with the transverse galaxy momentum field, constructed from the galaxy overdensity and a velocity field reconstructed from it. Restricting the reconstruction to large-scale modes strongly suppresses contamination from small-scale astrophysical foregrounds while preserving the ML signal. In this work, we develop a theoretical framework for the ML estimator and its foreground contamination. We show that the signal and foregrounds have a distinct dependence on the direction of the long-wavelength mode represented by the reconstructed galaxy velocity. In the squeezed limit enforced by the velocity reconstruction filter, the bispectrum sourcing the foreground correlation becomes independent of this direction, causing its leading contribution to vanish after angular averaging. In turn, the ML signal survives by matching this directional dependence, with its amplitude reduced only by the filtered velocity variance. Using the halo model, we derive the leading corrections beyond the squeezed limit and show that the residual contamination remains parametrically suppressed. Finally, we model the cross-correlation exactly in the curved sky, and show that it is sourced solely by the longitudinal component of the galaxy momentum field, in the form of a spin-1 $E$-mode, with all other contributions either strongly suppressed on small scales or exactly zero.

The LZ experiment has reported a single 248 keV nuclear recoil event in their $2.84$ ton-yr exposure on a xenon target. The absence of lower nuclear recoil events is naturally explained by inelastically scattering dark matter (DM), for which the nuclear recoil energy spectra are shifted upwards from zero. Excitingly, the observed rate and spectra are consistent with a thermal $1.1$ TeV higgsino that explains all of DM with a mass splitting $δ\sim 370 - 490$ keV between the two lowest-lying Majorana states, with the range depending on the underlying DM velocity distribution. In this work, however, we show that higgsino DM with these properties would generically produce more higher-energy recoil events in a high-energy sideband along the LZ $S1c$ energy axis, while the LZ experiment reports no events in that bin. This suggests a possible tension which could be further explored by the LZ Collaboration, as the acceptance in that bin is not public. We further show that a non-thermal $\sim$500 GeV higgsino may evade the sideband constraint and that future experiments using heavy targets can potentially test the higgsino interpretation of the event either way.

We study the shape of three-dimensional and projected dark-matter halo profiles extracted from cosmological $N$-body simulations in $Λ$CDM and non-standard cosmologies, using the \texttt{DUSTGRAIN-PF} and \texttt{DEMNUni} suites. The models considered include massive neutrinos, $f(\mathcal{R})$ gravity, and dynamical dark energy. By comparing density, mass, velocity-dispersion, and excess-surface-density profiles up to $5\,r_{500{\rm c}}$, we quantify the differential imprint of non-standard physics on halo structure in view of \textit{Euclid} cluster WL studies. Our main analysis is performed at $z=1.1$, a high-redshift regime where the weak-lensing signal-to-noise starts to degrade, providing a conservative stress test for detectability; for \texttt{DUSTGRAIN-PF} we additionally analyse $z=0.5$ and $z=0.3$ snapshots. In low-mass haloes ($M_{\rm 200c}<7\times10^{13}\,M_\odot$), $f(\mathcal{R})$ gravity produces deviations of order $10\,\%$ in projected and three-dimensional profiles, especially in the outskirts where screening is less efficient. Massive neutrinos partially reduce this signal, reflecting the competition between free streaming and fifth-force-enhanced growth. Dynamical dark energy and massive-neutrino cosmologies generally induce smaller, few-percent deviations, with the largest effects again found in low-mass haloes and at large radii. Under simplified assumptions for \Euclid WL, detecting such profile differences at $z=1.1$ requires stacks of $\sim10^5$ haloes, while a few thousands haloes may be sufficient at $z\lesssim0.5$. This further calls for the need of integrating such precise modelling of non-standard effects -- along with other observational effects -- in any likelihood involving \textit{Euclid} WL masses to avoid non-negligible systematic biases. Concentration--mass relations show weaker cosmology dependence, typically at the $\sim5\,\%$ level. [...]

Time-delay cosmography provides an independent probe of the Hubble constant ($H_0$), but exploiting the high-resolution imaging of the James Webb Space Telescope (JWST) requires strong-lens models that capture complex source-galaxy morphologies and the point spread function (PSF) with high fidelity. Flexible source and joint PSF reconstruction are computationally expensive, motivating more rigid parameterizations that can bias the inferred time-delay distance $D_{Δt}$. We present ALPACA, a semi-automatic, GPU-accelerated, differentiable strong-lens modeling pipeline built on Herculens, Starred, and JAX. ALPACA jointly reconstructs a pixelated correlated-field source, the lens mass, and the PSF within a single likelihood, reducing per-system runtimes from weeks to less than 15 hours. We validate the pipeline on realistic mock JWST/NIRCam observations built from morphologically complex COSMOS galaxy sources. Holding an imperfect PSF (RMSE = 12% from truth) fixed biases $D_{Δt}$ by up to 14%, whereas joint PSF reconstruction reduces this bias to the percent level. We further show that the reconstructed PSF core can be marginalized over jointly with the lens parameters during posterior sampling. Propagating the per-lens posteriors of a ten-lens mock sample through hierarchical inference, we recover the input cosmology without significant bias. ALPACA shows that rigorous control of source and PSF systematics is computationally affordable, supporting reliable percent-level precision in the JWST era.

The Dark Energy Spectroscopic Instrument is conducting a redshift survey, mapping the universe in three dimensions to measure the history of cosmic expansion. As well as extragalactic objects, it is targeting millions of Milky Way stars, including white dwarfs. Using Data Release 1 we assemble the largest catalogue of spectroscopically-confirmed white dwarfs to date, whose straightforward selection function enables statistically-robust population studies. We visually inspect and fit models to spectra of 63968 objects, finding 44409 white dwarfs. We present their spectral classifications, atmospheric parameters and radial velocities. We assess survey completeness and uniformity, identify potential spectral contamination caused by flux from nearby sources entering the fibre, and assign confidence scores to our classifications to facilitate selection of statistical and observational samples. We present spectra representing most white dwarf classes, common and exotic. We conclude with recommendations and warnings regarding the use of the catalogue.

Constraints on the effective viscosity of the intracluster medium (ICM) come almost exclusively from imaging: the power spectrum of X-ray surface brightness fluctuations, and the morphology of Kelvin-Helmholtz rolls at cold fronts. Both indicate a viscosity below the isotropic Spitzer value. XRISM/Resolve opens a complementary, purely kinematic route, in which the viscosity is inferred from the thickness of the viscous shear layer that develops at a sliding contact discontinuity. We apply it to the prototypical cold front of A3667. The interface is the classical Rayleigh problem, in which momentum diffuses to a depth d in a time T, giving a kinematic viscosity ν~ d^2/(4T). The XRISM velocity map shows a line-of-sight velocity jump of 535 km/s across the front, completed within the 100 kpc-wide region immediately inside it, the only region with an anomalously broad line. Hence d ~< 100 kpc. Adopting interaction times set by the sizes of the interacting regions divided by this velocity, we obtain dynamic viscosities of 5.2 x 10^3 and 3.3 x 10^3 g/cm/s for the cool and the hot side of the front, or 6.4 and 0.72 times the Spitzer value. Because d is limited by the angular resolution of XRISM, and the viscosity scales as its square, these are upper limits. The kinematic constraint is consistent with, and independent of, the imaging results.

White dwarf asteroseismology aims to constrain the interior structures of pulsating white dwarf stars by fitting observed pulsation periods to those computed for stellar models. However, there remain several obstacles for achieving reliable results with the model grid-fitting approach. We simulate data for a parameter recovery experiment in one and two dimensions (mass and effective temperature) to demonstrate improved methodologies that address challenges for characterizing degenerate asteroseismic solutions. We show that interpolating model periods onto a finer grid can adequately resolve asteroseismic solutions that would be missed with sparse model grids. Incorporating absolute magnitude from Gaia astrometry into the statistical fitting is shown to reduce solution degeneracy. We define a seismic solution as a solution to the mode identification problem, and we show that fitting to consistent model pulsation modes can isolate each candidate solution in degenerate solution space. A criterion based on $χ^2$ for what should be considered a reasonable fit is adopted, and we identify all combinations of model periods that meet this criterion, not only those nearest to the measured periods. Finally, we demonstrate how fitting Gaussians to probability distributions allows for the robust characterization of each candidate solution, including uncertainties. Our accurate characterization of the degenerate solution landscape is supported by the comparison of our solutions to a direct marginalization of the likelihood function. These fitting approaches can be generalized to higher dimensions, where there are more than two free parameters.

Cool-gas inflows, required to sustain star formation, have been fundamental in simulations yet remained observationally elusive. Using DESI spectroscopy of ~30,000 galaxies, we identify coherent inflowing gas (~100 km/s) in 20-50% of the sample, yielding a population-level census of gas flows. We uncover a striking inversion: inflows are detected in quiescent galaxies, whereas star-forming systems are dominated by gravitationally bound outflows. At fixed age, galaxies with inflows, outflows, or no/weak flows share similar masses, environments, and structures, indicating that these properties do not differentiate flow states. Instead, gas-flow state is linked to stellar population age and recent evolutionary history, consistent with age-dependent gas flows in two regimes. In some star-forming galaxies, elevated star formation surface densities drive outflows that recycle on ~0.5 Gyr timescales, consistent with a galactic fountain. In quiescent systems, low-level ``drizzling'' inflows persist, consistent with slowly cooling enriched halo gas and weak radio-mode nuclear activity. Broad gas-phase metallicity distributions---and absence of a pristine dilution signature---indicate that detected inflows are predominantly recycled or enriched. Detectability is modulated by dust, ionization, and geometry: in star-forming disks, inflowing gas lies near the disk plane and is obscured or ionized, while outflow hosts exhibit higher dust and metal content. As star formation declines, cold-outflow signatures weaken, and recycled or slowly cooling gas is more readily detected as inflow. Post-starburst galaxies provide snapshots of this transition. Our results resolve the scarcity of observed inflows, provide evidence for widespread gas accretion and recycling in present day galaxies, and establish an observational framework linking gas flows to star formation, chemical evolution, and galaxy structure.

H I gas traces the large-scale structure and provides the primary fuel for star formation. High-$z$ protoclusters are ideal laboratories to study how H I gas is accreted and consumed during the formation of the most massive structures in the Universe. However, much remains unknown about the distribution and physical state of their H I gas. We examine a rare configuration in which a protocluster candidate is located in front of a quasar at $z=3.09$. Our spectroscopic campaign confirms a protocluster at $z=3.079$: however, no corresponding strong H I absorption is found in the background quasar spectrum. Instead, we serendipitously discover a prominent H I absorption feature at $z\sim 3.01$, offset by $\sim 60$ cMpc from the centre of the protocluster. Spanning an exceptionally broad velocity range of $\sim 2000$ km s$^{-1}$ ($\sim 40$ cMpc), this absorption is decoupled from the confirmed member galaxies. Detailed kinematic modelling reveals this absorption comprises five distinct components rather than a single cloud. Moreover, one of these components exhibits a super-solar metallicity ($[\mathrm{O/H}] = +1.19^{+0.91}_{-0.78}$). We propose two physical scenarios for this unique system: (1) an additional, hidden massive protocluster along the line of sight, and/or (2) metal-rich outflows and metal-poor inflows driven by a single massive galaxy. The discovery highlights that while protoclusters are not universally associated with strong H I absorption, targeting the strong H I absorbers may serve as a beacon for uncovering massive, metal-rich protoclusters or complex gas kinematics in the early Universe.

The standard linearly perturbed Friedmann-Lemaître-Robertson-Walker (FLRW) framework provides an accurate description of the Universe on large scales but cannot fully capture the non-linear inhomogeneities that characterize the local Universe. In this work, we investigate whether an exact inhomogeneous solution of the Einstein field equations can provide a more realistic description of the nearby cosmic environment. To this end, we employ the quasi-spherical Szekeres class I solution with $Λ$ and focus on its axisymmetric subclass, motivated by recent observations indicating an axisymmetric expansion rate of the local Universe. We determine light propagation by numerically integrating the Sachs optical equations, thereby obtaining the luminosity distance and its evolution with redshift. We further compute the multipoles of the covariant cosmographic parameters and apply the formalism to an attractor-repeller configuration representative of the local matter distribution. We find that the resulting covariant cosmographic parameters are consistent with current observational constraints. Finally, we derive the gravitational potential and peculiar velocity field generated by the Szekeres density distribution, establishing the connection between the exact relativistic solution and its Newtonian counterparts.

We study the stability of collective neutrino oscillations using a distributional approach motivated by the statistical mechanics of Kuramoto synchronization. Treating the ensemble of neutrino flavor polarization vectors in the thermodynamic limit $N\to\infty$, we derive an exact nonlinear Fokker--Planck (continuity) equation for the one-body distribution $F(\vec{\mathbf{S}},ω,t)$ on the flavor sphere. This equation admits a two-parameter family of azimuthally symmetric stationary solutions, whose stability we analyze by linearizing around them. The resulting eigenvalue condition determines the growth or decay rate of small perturbations from \emph{any} initial distribution -- not merely from a state close to full flavor coherence -- thereby going significantly beyond the conventional linear stability analysis of collective modes. In special limits the condition reproduces known synchronization thresholds in the two-beam model, providing a non-trivial check of the framework. We present analytical results for the eigenvalue equation and explore stability phase diagrams for physically relevant frequency distributions.

Supermassive binary black holes (SMBBHs) are expected to be surrounded by circumbinary disks (CBDs) which affect the binary through gravitational forces and accretion. It has been reported that the binary can experience ``preferential accretion'' where one black hole (BH) out-accretes the other for hundreds of orbits, but this asymmetry has yet to be fully described or understood. In this work, we utilize a suite of 80 SMBBH hydrodynamical simulations with varying mass ratios ($q_b$) and eccentricity ($e_b$) in order to robustly delineate the behavior of preferential accretion, determine its relationship to the structure of the CBD, and study its observational consequences. We characterize the accretion-rate ratio $λ(t) \equiv \dot{M}_2(t)/\dot{M}_1(t)$ and the mass-ratio rate of change $\dot{q}_b \equiv d/dt(M_2/M_1)$ across the suite. We confirm that the secondary tends to out-accrete the primary ($λ\geq 1$), and find this preference to be strongest for low-$e_b$, low-$q_b$ binaries and increasingly time-variable toward high $e_b$. We also find that (i) the time-variability of $λ$ tracks the precession of the CBD, (ii) there can be sub- and super-Eddington accretion in a single binary, and (iii) the gas-driven approach toward equal mass becomes particularly slow for highly eccentric, high $q_b$ binaries, suggesting that some binaries may not reach $q_b=1$ within the $30\,\mathrm{Myr}$ lifetime of a quasar and therefore allowing LISA to constrain the accretion history of SMBBHs. Our findings also suggest that periodically flickering jets are a potential observable signature of many binaries.

We generalize the apparent-horizon thermodynamic construction of [1] to a FLRW universe with non-zero spatial curvature, non-interacting cold dark matter and holographic-type dark energy. For nonzero curvature, the scale factor is an additional geometric variable in the horizon equation of state. The usual criticality conditions are thus not well defined until a closure prescription for the curvature sector is provided. We introduce a dimensionless curvature variable and restrict the thermodynamic variations to slices of constant value of this variable. For each such slice, a positive holographic coupling and nonlinear powers larger than one guaranty the existence of a unique positive critical point. The critical specific volume, temperature and pressure are shifted by the spatial curvature, while the critical ratio and the mean-field critical exponents remain unchanged. It thus rescales the critical quantities, without generating a new local universality class. We construct a Helmholtz potential out of the physical thermodynamic volume and the entropy conjugate to the rescaled horizon temperature. In the quadratic holographic model, the parametric coexistence curve which is equivalent to the Maxwell construction is derived from the equality of the Gibbs free energies of the competing branches. The associated latent heat disappears at the critical endpoint, resulting in a global first-order coexistence in the fixed-curvature ensemble. Whether a physical FLRW trajectory intersects the critical or coexistence locus is a separate dynamical question because the dimensionless curvature variable generally evolves during cosmological expansion.

We derive the consistency conditions governing the propagating content of linear shear-free perturbations of Friedmann-Lemaître-Robertson-Walker cosmologies in metric $f(R)$ gravity. The divergence of the shear-free constraint is shown to be the total momentum-conservation equation and therefore does not imply geodesic flow. Its projected time derivative instead supplies a nontrivial integrability condition. Scalar, vector, and tensor sectors must then be tested separately. Shear-freeness removes the independent transverse-traceless electric-magnetic Weyl pair, so no ordinary $+$ or $\times$ tensor gravitational wave survives. The scalaron is not removed algebraically, but its four-variable harmonic system must remain in the time-dependent consistent subspace generated by the shear-free constraint and all of its time derivatives. For a geodesic shear-free congruence in the spatially flat expanding de Sitter patch, this excludes every nonzero fixed comoving scalar harmonic, including the healthy $R+αR^2-2Λ$ model with $α,Λ>0$. Vector modes obey a curvature-dependent global eigenvalue condition whose right-hand side is negative for positive-density matter with $f'>0$ and $1+w>0$; hence no nonzero vector harmonic survives on that branch. The familiar coasting $R^3$ example evades this sign obstruction only because it lies on an $f'<0$ branch. Thus tensor radiation is absent in the imposed shear-free sector, whereas scalar radiation is not excluded kinematically but remains a model- and background-dependent constraint-closure question.

The James Webb Space Telescope (JWST) has unveiled the existence of numerous z>10 bright ultraviolet (UV) galaxies, potentially challenging galaxy formation models in a Lambda Cold Dark Matter (LCDM) universe. Modifications to star formation and stellar feedback models have been suggested to alleviate the tension. However, the fundamental challenge is to design a galaxy formation model that simultaneously reproduces observations from z=0 to those at the highest redshifts. Here, we present predictions from the Shark semi-analytic model of galaxy formation, which is tuned to reproduce the z=0 universe. We show that the same model is capable of reproducing the current UV luminosity function constraints even up to z=17 without the need to invoke variations in the baryon physics model. This model is also capable of reproducing reasonably well the stellar mass function evolution from z=0 to z=10 and the cosmic star formation rate (SFR) density at 0<z<15, again with the same baryon physics and parameters. The key to the success of Shark is the triggering of starbursts by violent disk instabilities included in the model, which leads to an increase in the variance of the UV luminosity with increasing halo mass at z>10. We demonstrate that without this mechanism, galaxies are not bursty enough in the model to reproduce the observations at z>10.

The Trans-Neptunian objects (TNOs) formed together with the planets from a flat disc of gas and dust but now orbit the Sun mostly on inclined, eccentric orbits. One explanation for the TNOs' orbits is a close flyby of another star. One with a mass of $M_p =$0.8 Msun at $q_p =$ 110 au and $i_p =$70° fits the observations. However, such close encounters were more frequent when the Sun was young and still part of its birth cluster. Assuming that the flyby happened then, we use numerical $N$-body simulations to model how the TNOs' orbits changed due to interactions with the planets over the 4.56 Gyr since the Sun formed. We find that the cold Kuiper belt region lost about 80% and the hot Kuiper belt population about 40% of its initial population. About 7% --8% of the TNOs were injected into the planet region, but almost all (99%) were ejected afterwards. By contrast, the orbital parameters of distant Sedna-like objects ($q >$60 au) remained basically unchanged. Surprisingly, in sum, these changes improve the fit to the observed TNO population even more, strengthening the argument for a close flyby to the Solar System. The remaining differences concern the location of the cold population and the inclination distribution of the Sedna population. We discuss possible reasons and steps to resolve these issues.

The cosmic web contains most of the matter in the Universe, and its filamentary structure is closely related to the properties of dark matter. Compared with the cold dark matter (CDM) model, warm dark matter (WDM) suppresses the formation of low-mass halos and produces smoother cosmic filaments. In this work, we explore whether filamentary 21cm emission can be used to distinguish between these two dark matter models at $z=4$ and $2.5$. We generate intrinsic 21cm brightness temperature maps from cosmological hydrodynamical simulations and produce mock observations by convolving them with Gaussian beams of different angular resolutions. We then compare the predicted filament signals with the expected sensitivities of SKA1-Low and SKA2-Low. We find that the intrinsic 21cm emission clearly reflects the different filament morphologies in the CDM and WDM models, particularly at $z=4$. However, these differences rapidly disappear as the beam size increases, especially for beam sizes larger than 6 arcsec. Although larger beams improve the observational sensitivity, they also erase the small-scale structures that encode the dark matter information. As a result, neither SKA1-Low nor the planned SKA2-Low can directly resolve individual filaments with sufficient sensitivity to distinguish between the two models. We therefore conclude that direct imaging of cosmic filaments and constraining dark matter through their morphology are unlikely to be feasible with current and planned SKA-Low facilities.

The external environment where protoplanetary disks are embedded regulates disk evolution. External irradiation, which heats and evaporates gas, as well as frequent stellar encounters, which truncate disks, can reduce disk sizes and masses. Disk surveys in nearby star-forming regions with various environments can help us understand how external environments impact disk evolution. The Serpens star-forming region, which is thought to be dense but not highly irradiated by massive stars, is an ideal laboratory to study the dynamical effects on disk properties. We aim to study how dynamical interactions modify disk dust masses by comparing Serpens with other nearby star-forming regions through their disk masses and 3D stellar densities. We survey 321 young stellar objects from Class I to III in Serpens using ALMA at 0.25 arcsec. We measure disk dust masses from millimeter observations under the optically thin assumption and recompute the 3D stellar density using Gaia data. We apply the same method to other nearby star-forming regions for direct comparison with Serpens. The cumulative disk dust mass distribution of Serpens is similar to those of similarly aged nearby star-forming regions, such as Lupus and Taurus. Along with re-assessment of the 3D stellar density, it suggests that Serpens is not likely to have experienced strong tidal truncation capable of producing lower disk dust masses. We also find that the low disk dust mass tension in Ophiuchus and Corona Australis can be resolved when only disks that have been explicitly identified as members of young (1-2 Myr) sub-clusters are considered. Disks without Gaia identifications, especially in Ophiuchus, are spatially clustered in the 2D projected sky plane and less massive than nearby disks with Gaia identifications, possibly tracing the effects of tidal truncation in denser environments at earlier evolutionary stages.

High-resolution cross-correlation spectroscopy (HRCCS) is a powerful technique for detecting molecules whose individual spectral lines are too weak to be identified directly, but its sensitivity is limited by the availability and quality of high-resolution opacity data. Many molecules of atmospheric and astrobiological interest lack complete line lists, restricting traditional template-based searches. In this work, we use Titan as a controlled testbed to develop and validate a new cross-section-based methodology for HRCCS template construction. We analyse K-band CRIRES+ observations of Titan (1.99 - 2.48 μm), and compute cross-correlation functions using both line-by-line and cross-section-based templates. Our analysis recovers known hydrocarbons such as methane (CH4) and acetylene (C2H2), and yields the first HRCCS detection of ethane (C2H6) with a significance of SNRpeak = 5.17 {\pm} 0.07. The ethane detection was made possible exclusively through cross-section-based templates, as no high-resolution line list currently exists for this molecule. These results demonstrate that cross-section-based template construction is a practical and powerful strategy for extending HRCCS to molecules that currently lack reliable line lists, and establish Titan as a benchmark for calibrating molecular detection techniques that can be applied to both solar system and exoplanet atmospheres. Future applications of this approach to other ground-based high-resolution spectrographs, as well as to JWST's highest-resolution modes and next-generation facilities such as the ELT, could significantly expand the inventory of molecules detectable in planetary atmospheres.

We propose a manifestly conformally invariant off-shell definition of the Weyl tensor in Galilean geometry. We also propose definitions for the electric and magnetic parts thereof. For the latter to vanish, observers with specific kinematical properties need to exist. While this is guaranteed by the Newton--Cartan equation, it might not be true for other Galilean invariant theories. Therefore, we propose that imposing the existence of observers for which the off-shell magnetic part is zero should be a necessary condition for a Galilean invariant theory to be called `Newtonian', in the spirit of the `Newtonian' condition, introduced by Trautman, in standard Newton--Cartan gravity. As a side result, we also show that there exists a unique Galilean boost-invariant connection that can be built from a Galilean structure and a choice of Coriolis field, even when the clock form is not closed. No extra structure is needed, contrasting with the standard approach used to construct boost-invariant connections which introduces a mass gauge field.

It has long been assumed that a particle---having crossed a black hole event horizon unentangled with another particle in the exterior universe---can no longer dynamically entangle with that particle. Here, we demonstrate within a gravitational retarded-potential model that entanglement can in fact be created from scratch between freely falling spatial superpositions across an event horizon. Yet, extracting this state radially requires non-inertial deceleration, triggering soft-graviton bremsstrahlung, and establishing a strict dephasing bound in terms of entangling phase, $Γ\ge\frac{729}{160π}Φ$. This decoheres the entangled state. By contrast, equivalent macroscopic optical masses allow local tangential harvesting of entanglement via quantum erasure, thus revealing a remarkable geometric duality: Spacetime irreversibly degrades entanglement the moment the localized mass is dragged away from the horizon, while allowing the transverse teleportation of its entangled state to infinity.

In this work, we constructed a new metric describing a rotating cosmological black hole with electric and magnetic charges in the presence of quintessence dark energy, where the cosmological constant and quintessence are associated with external matter. Starting from a Schwarzschild-type metric and considering all energy contributions in the function $f(r)$, we introduced rotation through the Newman--Janis algorithm using a unique and general complexification rule, avoiding possible ambiguities. Assuming that the metric satisfies the Einstein field equations with matter, we calculate the event horizon condition, the total stress--energy tensor $T_{μν}$, and the Kretschmann and Ricci scalars. We also analyze the influence of the cosmological constant $Λ$ and quintessence parameter $α$ on the event horizon and ergosphere. In the limit $α\rightarrow 0$, the scalar curvature differs from the vacuum result, $R \neq -4Λ$, while the singularity region remains unchanged. Finally, we investigate the angular and rotational velocities of a test particle, as well as the energy conditions, through the energy density and pressures required for the existence of this black hole solution.

Minimal electroweak dark matter (DM) extended by a Majorana (M) and a Dirac (D) multiplet coupled through Higgs interactions provides a predictive framework for high-energy inelastic nuclear recoils. We show that this setup can account for the type of event recently reported by the LUX-ZEPLIN (LZ) Collaboration. At the custodial point, electroweak symmetry breaking induces a neutral-state splitting $δ$ and an off-diagonal $Z$ interaction. For the coupled multiplets $3_M2_D$, $5_M4_D$, $7_M6_D$, $9_M8_D$, $11_M10_D$, and $13_M12_D$, both the splitting and the leading inelastic interaction are universal at fixed $(m_χ,y)$, independent of the electroweak representation. The observed high recoil energy points to splittings of a few hundred keV, while the fixed $Z$-mediated rate allows us to determine a mass-dependent $δ_{\rm LZ}(m_χ)$ by requiring one expected inelastic event in the LZ exposure, with a two-sided $90\%$ Poisson band. This defines a universal region in the $(m_χ,y)$ plane. Intersecting it with the representation-dependent thermal relic trajectories selects a benchmark for each multiplet. The corresponding electroweak representation then predicts a correlated loop-induced elastic spin-independent signal at lower recoil energies.

Primordial non-Gaussianity (PNG) offers a powerful test of inflationary physics beyond the simplest single-field scenarios. In large-scale structure, biased tracers respond to local-type PNG with a distinctive scale-dependent signature that is strongest on large scales, but single-tracer measurements are fundamentally limited by cosmic variance. Multitracer methods can mitigate this limitation; however, many existing proposals require splitting tracers using a wide range of halo masses which are often not available, or using secondary halo properties that are not directly observable. Here we develop and validate an alternative strategy based on splitting a tracer sample by its large-scale dark matter environment. We derive a theoretical framework for the $f_{\rm NL}$ response of environmentally selected tracers and test it in $N$-body simulations, confirming consistency between separate-universe predictions and direct power-spectrum inference. To enable observational implementation, we show that a simple linear-theory reconstruction of the dark matter environment from halos yields unbiased $f_{\rm NL}$ constraints in $\mathtt{Quijote}$ simulations. Forecasts for the DESI LRG sample indicate that environmental multitracing can improve constraints on local PNG by a factor of 2-3 relative to conventional single-tracer analyses for high density tracers, comparable to gains achieved when halo formation time is assumed to be perfectly known. We develop an analytic model explaining why the gains are smaller for a halo environment split and why separate-universe predictions fail in this case. Unlike approaches that require modeling assembly bias or inferring secondary halo properties, our method relies only on linear theory. These results establish environmental multitracer techniques as a promising route toward $σ(f_{\rm NL})\sim 1$ with upcoming galaxy surveys.

We develop an analytical treatment of the emergent cosmological dynamics induced by local polynomial interactions in a deparametrised group field theory model, moving beyond the usual non-interacting approximation. For a single field mode with quartic and sextic couplings we obtain a closed-form generalised Friedmann equation within a controlled Gaussian regime. The dynamics preserves the quantum bounce while generating effective dark matter and dark energy contributions as collective quantum-geometric phenomena, providing the first derivation of both components from a single model. Matching the resulting dark-energy-to-dark-matter density ratio to observations places concrete phenomenological constraints on the fundamental theory, linking the cosmic coincidence problem directly to the underlying quantum gravity dynamics. Our results extend to arbitrary even polynomial interactions, establishing a systematic dictionary between microscopic interaction orders and effective equations of state.

The recent analyses of DESI DR2 BAO, Planck CMB and various Supernovae datasets have shown preferences for evolving dark energy at various levels of significance, under frequentist model comparison tests. However, the same analysis done with a pure Bayesian model comparison test can lead to the opposite conclusion due to the impact of $w_0w_a$ priors. In contrast to these in-sample model comparison route, we approach the problem along a third axis -- which model better predicts unseen data. We use the Expected Log Predictive Density (ELPD) metric to score the predictiveness of $Λ$CDM and $w_0w_a$CDM models, using the leave-one-redshift-block-out cross-validation estimator. The comparison metric $Δ$ELPD is out-of-sample, unlike $Δχ^2_{\rm MAP}$ and Bayes Factor, and insensitive to prior width, unlike the Bayes factor. Aggregated, we find modest preferences for $w_0 w_a$CDM, primarily driven by the BAO. Decomposing the $Δ$ELPD score by redshift we find that a single BAO block, LRG2 ($z = 0.706$), supplies the entire BAO preference. Interestingly, the biggest outlier point, LRG1, contributes little to the model predictive scoring because both $Λ$CDM and $w_0w_a$CDM fail to predict the observation by equal amount. On the other hand, all SNIa datasets' $Δ$ELPD scores are consistent with $0$, indicating that at an out-of-distribution predictive level, the $w_0w_a$CDM tension is entirely driven by one BAO point, not SNIa.

We develop a branch-complete classical free-spectrum analysis of local parity-even symmetric teleparallel gravity containing all pure-gravity nonmetricity operators through four derivatives. After canonicalisation and integration by parts, the nonlinear action contains 116 independent operators distributed as $5+13+29+69$ among $Q^2$, $(\nabla Q)^2$, $Q^2\nabla Q$ and $Q^4$. Around Minkowski spacetime in the coincident gauge only the first two sectors contribute to the quadratic action, which reduces to four two-derivative and five four-derivative bilinears. We map the original coefficients to this nine-dimensional basis, construct the Hessian and momentum kernel, and classify local linear gauge symmetries. Unrestricted linearised diffeomorphism invariance imposes six independent relations and leaves a three-parameter family $(A,B,C)$. We also recover TDiff, Weyl, WTDiff and additional accidental symmetry branches. Barnes--Rivers decomposition, gauge fixing, exact inversion and conserved-source saturation reduce the physical exchange to spin-two and scalar factors $A+2Cq$ and $2A+(3B-2C)q$. The generic regular spectrum contains a massless graviton, a massive spin-two state and a massive scalar, with eight degrees of freedom. The massive spin-two residue is necessarily opposite to that of the healthy massless graviton. The regular classification leaves only GR/STEGR, $A>0$ with $B=C=0$, and the scalar extension, $A>0$, $B>0$, $C=0$, as fully healthy Minkowski loci. The latter has three degrees of freedom and mixed ultraviolet behaviour, with tensor exchange scaling as $k^{-2}$ and scalar exchange as $k^{-4}$. Thus, within this finite local metric theory, full four-derivative tensor suppression requires the finite massive spin-two pole with opposite residue.

In this work, we investigate exactly the dynamics of a relativistic spinless particle influenced by a static and uncharged black hole surrounded by a quintessence-like anisotropic fluid or Kiselev black hole. Considering two quintessence-like models as examples, we calculate the radial wave function and determine, in both cases, the quasispectrum of energy, as well as the Hawking radiation and temperature. We find that the stronger the influence of quintessence-like anisotropic fluid, the smaller the radiation observed outside the event horizon. Furthermore, the temperatures obtained in both cases are directly influenced by the quintessence-like anisotropic fluid, and we recover the values obtained in other contexts, such as those derived using the surface gravity. Finally, in the absence of quintessence, we recover in both cases the Hawking temperature $T_{H} = 1/(8πk_{B} M)$ corresponding to the Schwarzschild black hole.

In this work, we provide the first derivation of the line-intensity mapping (LIM) power spectrum super-sample covariance (SSC) from first principles, and also derive as a by-product the non-Gaussian in-box contributions to the covariance for the first time. Previous studies have typically modelled the LIM power spectrum covariance using either the Gaussian approximation or estimates obtained from mocks or the data itself, neglecting uncertainties related to whether the limited volume surveyed sits in a cosmological overdensity. This contribution, known as the SSC or, depending on the context, the field-to-field variance, cannot be estimated from the data, but it is crucial for a correct inference of global quantities, i.e., for ensemble-averaged parameters rather than the actual values just within the patch of the Universe observed. For our derivation, we employ a combination of the halo model and standard perturbation theory that allows us to capture the nonlinearity and non-Gaussianity of the covariance. After a successful validation of our predictions against painted N-body simulations, we explore different scenarios related to current and future LIM experiments, quantifying the relative importance of the non-Gaussian in-box and SSC. We find that the newly derived contributions to the LIM power spectrum covariance are crucial at intermediate and small scales, especially for cases in which the covariance is not dominated by instrumental noise. We find that the relative relevance of the SSC with respect to the other covariance contributions is roughly independent of the survey volume, but does depend on the specific response of the power spectrum to large-scale modes for each line and redshift. Therefore, the impact of the SSC will be increasingly significant for parameter inference from the current and the next generation high signal-to-noise LIM surveys.

Recently, a new Trojan Horse Method (THM) measurement of carbon-carbon fusion was reported by Li \textit{et al.} [Phys. Lett. B (2026) 140675]. The purpose of the present work is to demonstrate the breakdown of the plane-wave approximation used in the analysis of these data and the critical role of Coulomb distortions in the initial and final states. The reaction mechanism underlying the THM analysis of the $^{12}\mathrm{C}+{}^{12}\mathrm{C}$ fusion reaction using the $^{16}\mathrm{O}+{}^{12}\mathrm{C}\to α_s+α+{}^{20}\mathrm{Ne}$ reaction is investigated. Particular attention is paid to the spectator momentum distribution and to the dependence of the THM reaction amplitude on the relative carbon-carbon energy $E$. It is demonstrated that agreement with the measured spectator momentum distribution does not by itself validate the plane-wave approximation. Although the experimental momentum distribution can be reproduced, inclusion of Coulomb distortions in both the initial and final channels leads to an energy dependence of the THM amplitude that is completely different from the plane-wave result. Consequently, the energy dependence of the $^{12}\mathrm{C}+{}^{12}\mathrm{C}$ fusion cross section extracted from the THM data can be strongly distorted by the plane-wave treatment. It is concluded that the astrophysical factor extracted in the plane-wave analysis cannot be regarded as reliable and may lead to misleading conclusions concerning the low-energy $^{12}\mathrm{C}+{}^{12}\mathrm{C}$ fusion reaction.

Parametric component separation methods used in cosmic microwave background (CMB) analyses commonly assume instrumental noise to be uncorrelated between frequency channels. This approximation may become insufficient for experiments whose design naturally introduces cross-channel correlations, such as interferometric or multi-band bolometric systems. This work aims to provide a consistent framework to include frequency-correlated noise in Commander2 and to quantify its impact on component separation in controlled and instrument-motivated simulations. We generalized the amplitude sampling of the Commander2 framework to consistently handle dense frequency-frequency noise covariance matrices, modifying the likelihood evaluation and parallel architecture. The implementation is validated on low-resolution simulations of a simplified CMB-S4-like configuration. We applied the generalized framework to instrument-motivated simulations of the QUBIC bolometric interferometer, using representative noise models derived from end-to-end simulations. Validation tests show excellent convergence and agreement between sampled parameters and analytical maximum-likelihood estimates. Accounting for frequency correlations reduces the width of the maximum-likelihood residual distributions by up to 8% on average, and up to 11% for individual parameters, relative to the diagonal-noise approximation. In QUBIC-based simulations, accounting for an anti-correlated noise produces posterior widths reduced by about 60%, while a positive-correlation test increases them by about 40%. This shows that neglecting cross-channel correlations can lead to incorrect posterior estimates in either direction. The generalized implementation presented here enables the Commander2 software to consistently propagate such correlations into parameter estimation, providing a framework for assessing their impact in future CMB analyses.