Main

Residual agreement in a kinetic or moment equation does not automatically identify every higher-order closure variable entering a nonequilibrium shock. We formulate this issue as an observability problem for the fourth-order closure content of monatomic normal shocks and follow it through a hierarchy of collision models and diagnostics. The kinematic part of the result is independent of the collision operator: the one-dimensional heat-flux budget observes the projected fourth-order channel $S=R^{\cl}_{xx}+\Delta/3$, not the tensorial R26-level moment $R^{\cl}_{xx}$ separately from the scalar fourth-order excess $\Delta$. The observation map therefore has a one-dimensional null space, so a heat-flux residual can be small while the split between tensorial anisotropy and isotropic tail intensity remains wrong. A DVM-consistent scalar-excess budget supplies the missing channel and gives the two-channel reconstruction $R^{\cl}_{xx}=S-\Delta/3$ without direct $R^{\cl}_{xx}$ data. Across BGK shocks at Mach 2–5, this reduces the active-zone $R^{\cl}_{xx}$ error from about $63$–$64\%$ to $2.4$–$4.1\%$. Sparse scalar-excess interpolation is used only as an information-reduction test: a representative 24-probe operating point gives $R^{\cl}_{xx}$ errors below $4.5\%$, and below $4.7\%$ with $1\%$ probe noise. Collision-model diagnostics then separate the invariant observation channel from the model-dependent source law. Shakhov changes the heat-flux relaxation to the correct Prandtl number but is neutral in the even $|\boldsymbol c|^4$ scalar-excess source; a direct discrete Shakhov channel check recovers $S$, $\Delta$ and $R^{\cl}_{xx}$ with errors $6.4\times10^{-4}$, $2.1\times10^{-7}$ and $1.0\times10^{-3}$, respectively.

The chaotic rotation of Saturn's moon Hyperion is a paradigmatic example of Hamiltonian chaos in a natural system. Although its tumbling motion is well established theoretically, identifying a robust observational signature of chaos from sparse and noisy astronomical time series remains a major challenge, making phase-space reconstruction techniques impractical under realistic conditions. In this work, we show that multifractal detrended fluctuation analysis (MFDFA) provides an effective alternative for detecting chaotic dynamics directly from photometric observations. Using historical ground-based light curves and synthetic datasets, we demonstrate that the intermittency associated with chaotic tumbling produces a broad multifractal singularity spectrum. While multifractality is a known feature of Hamiltonian chaos, we show that it can serve as a practical observational diagnostic when traditional chaos indicators fail because of sparse sampling. In particular, the multifractal spectrum remains detectable after realistic observational filtering and distinguishes chaotic tumbling from aliased regular rotation. By contrast, regular resonant rotation exhibits a significantly narrower spectrum, approaching the monofractal behavior expected for uncorrelated noise. For the observational data, we measure a broad spectral width consistent with the synthetic chaotic model, statistically distinct from surrogate datasets, and robust against finite time-series length. These results establish multifractal scaling as a viable observational signature of Hamiltonian chaos in sparse astronomical datasets, bridging nonlinear dynamics and planetary photometry.

Polarized images of a black hole encode the direction of electromagnetic energy flow near its event horizon. Measuring polarization from near-horizon emission can determine whether this energy flow is powered by the accreting plasma or the black hole spin. Here we consider the linear polarization of synchrotron radiation emitted from the base of horizon-threading field lines in a time-stationary, axisymmetric, and degenerate Kerr magnetosphere. We show that the observed polarization pattern displays universal behavior: it is completely determined by the black hole spin and observer inclination and is independent of the magnetic field geometry. We derive a simple analytic formula for this spin-dependent horizon polarization pattern. We find that this predicted pattern is also approached in time-averaged images from General Relativistic Magnetohydrodynamic simulations. Future observations with Very-Long-Baseline Interferometry at microarcsecond resolution could detect the trend of polarization toward the unique horizon value in M87*. Such observations may enable new measurements of black hole spin and provide evidence that magnetic field lines thread the horizon and extract spin energy via the Blandford–Znajek process.

The Simons Observatory will be extended by three Small Aperture Telescopes by 2027, increasing the total number of these instruments to six. We study the prospects for probing the reheating temperature and the inflaton coupling with this configuration, assuming a discovery of primordial gravitational waves in benchmark scenarios with a tensor-to-scalar ratio r=0.0036 or r=0.01. In popular plateau models of inflation, such an observation would fix the scale of inflation and enable determination of the order of magnitude of the reheating temperature and the inflaton interactions. For QCD-driven Warm Inflation the reheating temperature and inflaton coupling to gluons could, under optimistic assumptions, be measured with a precision of a few percent. Such a measurement would imply a clear prediction for complementary inflaton searches in axion experiments, paving the way toward probing the mechanism responsible for the initial conditions of the hot Big Bang in the laboratory.

Recent analyses of $f\sigma_8$ and weak-lensing data indicate that the linear growth rate at $z\lesssim 1$ may be lower than predicted by $\Lambda$CDM. This motivates models of dark matter in which large scale structure growth slows relative to $\Lambda$CDM at late times. We construct particle models in which dark matter experiences a drag with dark radiation that grows at late times, unlike conventional DM–DR interactions, which fade as the universe expands. A key ingredient is that the radiation interacting with the dark matter is produced at late times from dark matter decay. An explicit model, interacting Decaying Cold Dark Matter (iDCDM), adds two parameters beyond $\Lambda$CDM while leaving the background, BBN, and primary CMB intact. But it predicts a step-shaped suppression of the linear growth rate $f(k,z)$, a distinctive target for DESI, Euclid, and Rubin. Confronted with current data, iDCDM shows a modest preference over $\Lambda$CDM, driven by $f\sigma_8$, with $\Delta\chi^2$ between $-2.7$ and $-7.6$ depending on the assumed scaling of the drag with redshift and on neutrino masses. The decisive test will come from upcoming $k$- and $z$-resolved growth measurements.

Gas giants may shape the reservoir of solids and gas in the inner disk in which the small planets assemble. To test this possibility, we collect a sample of 43 exoplanetary systems containing 68 inner small planets (ISP) with both measured masses (1-20 M$_{\oplus}$) and radii (1-4 R$_{\oplus}$). After correcting for heterogeneous individual system sensitivities to distant gas giants, we calculate the gas giant occurrence rate in ISP systems P(GG$|$ISP) as a function of inner small planet density, envelope mass fraction (EMF), and core mass. While we find no significant difference between P(GG$|$ISP) given high/low small planet density, EMF, or core mass, we see hints of a trend when only looking at the metal-rich systems. Despite the substantial limitations due to small sample sizes, we find that gas giants in metal-rich systems are preferentially found with lower density planets with similar core masses. We find consistent hints of trends using larger samples of inner planets with measured radii divided across the radius valley or with measured masses divided across 10 $M_\oplus$. Our result is consistent with more metal-enriched disks catalyzing rapid core assembly and kickstarting the gas accretion early, while the muted difference in the outer giant occurrence rate with respect to core mass may indicate contamination by post-formation photoevaporation.

Classical novae are transient events powered by thermonuclear burning in a layer of hydrogen-rich material accreted by a white dwarf from its binary companion. Most classical novae reach optical maximum within  1 d, but a rare few rise far more slowly. We probe the envelope structure and ejection mechanism of the slowly-rising nova PGIR22akgylf with TESS photometry spanning 3 to 16 d after the nova discovery, supplemented by ground-based observations that cover its full  133 d ascent to maximum. We detect a 0.1802 +/-0.0012 d periodic brightness modulation with a peak-to-peak amplitude of  0.02 mag, identified with PGIR22akgylf via temporal and spatial coincidence. The period is stable over the two weeks of TESS coverage, suggesting an orbital origin. Whether this period corresponds to the full or half orbital period, it implies a dwarf donor companion. At the time of the TESS observations the nova was > 6 mag above quiescence (but still 4 mag below peak), so its light should be dominated by the expanding photosphere. We interpret the periodic signal as arising from the binary orbital motion distorting the nova envelope while its size remains comparable to the binary separation. This interpretation points to common-envelope interaction as a contributor to shell ejection in PGIR22akgylf and demonstrates that the slow-rise phenomenon is not exclusive to thermonuclear eruptions in symbiotic binaries, where the large orbital separation of the giant companion inhibits such interaction.

We develop a new grid-based multi-grid Poisson solver in numerical relativity. We report the performance of the multi-grid Poisson solver in the initial value problems for two-puncture black holes, a static spherical neutron star, a uniformly rotating neutron star in equilibrium, and a gravitationally collapsing massive star. As a demonstration, we conduct a numerical-relativity neutrino-radiation-transfer hydrodynamics simulation of the gravitational collapse of the $9M_\odot$ massive star in Ref. \cite{Aguilera-Dena:2020mfh} up to the core bounce. During the simulation, we employ the constraint-preserving regird prescription with the newly developed multi-grid Poisson solver to improve the resolution. It shows that the baryonic mass, the Arnowit-Deser-Misner (ADM) mass, and the ADM-like angular momentum are, respectively, preserved with $O(10^{-3})\%$ and $O(10^{-2})$–$O(10^{-1})$\% accuracy.

The spectra of an accretion disk around black holes are the basic diagnostic tool to enlighten the underlying flows and then black holes. Accretion flows around black holes, however, are controlled by parameters like the magnetic field, spin of the black hole, accretion rate and temperature of the flow. These quantities affect the (magneto)hydrodynamics of the flow thus consequently lead to variations in the spectrum. We first consider numerical steady state magnetohydrodynamic (MHD) solutions of magnetized accretion flows around black holes to study the dependence of the spectra on these disk properties. The spectrum exhibits strong dependence on the spin of the black hole, accretion rate, magnetic field and the electron temperature of the flow. Variations in these quantities influence the emission peaks and overall luminosity, which can be a tell-tale sign to extract physics of observed spectra. We further validate our results with general relativistic MHD (GRMHD) simulations using the standard and normal evolution (SANE) and magnetically arrested disk (MAD) vector potentials. We consider two black hole spins ($a=0.5$ and $a=0.9375$) to model the magnetic field configurations and study the resulting spectra by comparing MAD and SANE results. We find a large difference in the bolometric luminosities and the location of the emission peaks between SANE and MAD flows. Certain properties of the spectra, like, the ratio of synchrotron radiation to synchrotron self-Comptonization peaks in SANE and MAD, show drastically distinct features. The overall luminosity combined with such metrics can distinguish the magnetic field characteristics in astrophysical systems.

Accurate modelling of the covariance of clustering observables is essential to fully exploit current and future survey data, which is expected to constrain large-scale clustering signals with unprecedented precision. Computational costs of simulation-based estimates motivate analytical approaches, especially in light of the growing interest towards multi-tracer analyses and parity-odd signatures in two-point statistics, which respectively mitigate cosmic variance and probe relativistic projection effects on cosmological scales. In this work, we generalise previous theoretical results for the Gaussian covariance of multi-tracer power spectrum measurements, providing a general expression applicable to both real (even-parity) and complex (both even- and odd-parity) power spectra. We focus on a generic weighted estimator at first, and then showcase how our general formalism applies to Legendre power spectrum multipoles and two-dimensional power spectrum, recovering known limits in appropriate cases. We validate our predictions against Gaussian Monte Carlo simulations and investigate the structure of the covariance matrix, including the Hermitian properties of its imaginary part.

We describe spectroscopic tomographic weak lensing measurements (spectrotomography) for two rich clusters of galaxies, A1767 and A2065, based on extensive spectroscopy and Subaru/Hyper Suprime-Cam (HSC) imaging. These detections represent the first use of spectrotomography based on archival Subaru/HSC imaging. The measurements depend only on galaxies with spectroscopic redshifts reported here. The approach cleanly separates cluster members from the background and suppresses systematics that may be introduced by the use of photometric background redshifts. We detect the tomographic shear signals at $3.1\sigma$ (A1767) and $3.5\sigma$ (A2065). The shear signal amplitudes are consistent with the cluster dynamical (caustic) masses and they scale appropriately with source redshift. However, comparison with the first spectrophotometric detection of A2029 based on DECam imaging reveals some subtle potential systematic issues in deriving the shear signal for the relatively bright background galaxies used in the analysis. These issues may be important for understanding future more extensive applications of spectrotomography based on further Subaru imaging, as well as Euclid and LSST data. The total of three spectrophotometric detections (A1767, A2029, and A2065) sets the stage for broader application of the technique for unbiased cluster weak lensing mass determinations and potentially for a geometric cosmological test that is independent of other methods.

Solar wind stream interaction regions (SIRs) drive recurrent geomagnetic storms, yet most existing catalogs rely on expert inspection and simple thresholds that are subjective and can miss events with complex morphologies. We present SIREN (SIR Encoder Network), a lightweight Transformer based model for per timestep SIR detection from in situ solar wind observations. The model ingests sequences of 11 solar wind parameters, spanning magnetic field, velocity, and thermodynamic properties. With approximately 100,000 trainable parameters in a two layer encoder architecture, SIREN is trained using weighted binary cross entropy loss and a cosine annealing learning rate. Platt scaling is applied to produce well-calibrated detection probabilities. On a held-out test set of 102 events, the calibrated model achieves a ROC-AUC of 0.93, F1 score of 0.78, and true skill statistic of 0.67. Analysis of the self-attention weights confirms that the model concentrates on the SIR, grounding its decisions in the physically relevant portion of each sequence. Integrated Gradients attribution reveals a quantifiable feature hierarchy: proton density (24.3%) and magnetic field magnitude (21.6%) dominate, followed by temperature (13.9%) and bulk speed (12.1%). Notably, the transverse velocity component Vy and east-west flow angle together contribute 13-17%, identifying flow deflection as a consistent but previously under-quantified SIR signature. By producing continuous probabilities rather than binary labels, SIREN enables flexible threshold tuning for operational use and provides a template for compact, interpretable deep-learning systems in space weather.

High-mass X-ray binaries (HMXBs) provide a natural laboratory to study radiatively driven stellar winds under strong external X-ray irradiation. As the compact object moves along its orbit, the wind density and ionization structure vary with orbital phase, leaving characteristic signatures in X-ray emission and absorption features. The amplitude and morphology of this variability depend on the system geometry, including the orbital inclination (via line-of-sight and occultation effects) and the orbital eccentricity (via phase-dependent changes in the orbital separation). We present a computational framework that connects phase-resolved spectroscopic variability to the three-dimensional structure of irradiated winds, and that enables inference of physically meaningful wind–irradiation parameters within a Bayesian setting. We combine photoionization calculations with a parametric wind description and orbital geometry to construct three-dimensional maps of density and ionization. From these maps we compute orbital-phase-dependent diagnostics, accounting for geometric occultation and wind inhomogeneity through a clumping prescription. We then use Bayesian inference to compare model predictions with phase-resolved observables and to quantify parameter constraints and degeneracies. The framework reproduces the main orbital-phase-dependent trends expected for irradiated winds and yields robust constraints on the system ionization balance. While combinations of wind and luminosity parameters are well constrained, individual parameters can remain partially degenerate depending on the orbital configuration and data quality. This modular and computationally efficient approach provides a route to interpret HMXB phase variability in physical terms, and offers a foundation for future extensions toward forward spectral modeling and population-level applications.

We present an analysis of XMM-Newton and Chandra observations of 85 nearby main-sequence FGK stars with age estimates ranging from 0.2-12 Gyr. We measure quiescent 0.3-10 keV luminosities, variability metrics, and multi-temperature thermal plasma spectral parameters. Quiescent spectra are typically described by three characteristic plasma components ($kT\approx0.1$, 0.4, 0.8 keV); the fraction of flux from $T\ge7$ MK rises with X-ray surface flux, reaching $\sim$50% for $F_X\gtrsim10^6$ erg cm$^{-2}$ s$^{-1}$. We derive relations between emission measure-weighted coronal temperature and both $L_X$ and $F_X$, enabling temperature-informed count-rate conversions for faint sources. We quantify how bandpass conversions (ROSAT 0.1-2.4 keV vs. XMM-Newton 0.3-10 keV) depend on temperature, and show that inferred ROSAT-band $L_X$ broadly follows the canonical $t^{-1.5}$ decay, while the harder band exhibits increased scatter at $>$4 Gyr. Several stars show excess activity suggestive of age errors, inclination effects, or unresolved companions. Some of these "outlier" stars are potential direct imaging targets for the Habitable Worlds Observatory, and detailed characterization of these stars is needed to inform their likely influence on the atmospheric evolution of orbiting planets.

Strong lensing of massive black hole binaries can produce multiple gravitational-wave images with different magnifications and arrival times. LISA signals remain in band for months to years, allowing multiple lensed images to overlap during the inspiral stage and generate beat patterns. A singular isothermal sphere lens model is adopted to describe the lensing configuration, and two-image beat waveforms are constructed from massive black hole binary signals. To isolate the beat pattern itself, waveform mismatch is evaluated only during the overlapping inspiral stage before the coalescence of the first image, excluding contributions from the delayed merger peak of the second image. Using the HS-nod-SN (B+20) strong-lensing population, the occurrence rate of identifiable beat events is estimated, and Bayesian parameter estimation is performed with a beat template. Beat patterns are most readily identified when the lensing time delay is short and the delayed image has a relatively large magnification. Among 196 detectable two-image lensed events, 92 satisfy the temporal-overlap condition and 14 satisfy the beat-identification criterion, corresponding to an identifiable beat fraction of about 7\%. Posterior inference shows that the beat template can recover the lensing time delay and magnification parameters for a representative beat event. These results indicate that lensed beat patterns constitute a distinguishable subset of strongly lensed LISA events and provide a unique observational signature of strong lensing in the LISA band.

We identify a geometric obstruction to resolving the Hubble tension by combining early-time sound-horizon reduction with late-time smooth dark energy. Within $\Lambda$CDM, the BAO–SN matter-density gap $\Delta\Omega_m = 0.037$ is exactly invariant under the sound-horizon rescaling $\alpha \equiv r_s^{\rm mod}/r_s^{\Lambda{\rm CDM}}$, and late-time $w(z)$ deformations cannot eliminate this gap either: reconciling the two datasets requires \emph{opposite} deformations – phantom ($w < -1$) for BAO, quintessence ($w > -1$) for SN at $z < 0.5$ – an anti-alignment quantified by $\cos\theta = -0.97$ in $w(z)$ space. A full MCMC analysis of DESI DR2 BAO, Planck plik\_lite, and Pantheon+ bears this out: the optimal $\alpha^* = 0.992$ ($0.8\%$ $r_s$ reduction) brings the joint fit to $H_0 = 70.3 \pm 0.3\;\mathrm{km\,s^{-1}\,Mpc^{-1}}$, still $3.2\sigma$ below SH0ES, with the inter-dataset tension reduced but not removed. The obstruction reflects not a shortage of model freedom but an irreducible disagreement between probes. The deformation space $\{\alpha, \beta_{\rm damp}, w(z)\}$ already spans $93\%$ of the $\Omega_m$ response direction; nonetheless BAO and SN constrain $\Omega_m$ through independent channels and disagree, while the residual $H_0$ deficit, anchored by the local distance ladder, resides in the absolute distance scale that $w(z)$ reshapes but cannot rescale.

We present a data-driven investigation of Galactic diffuse emission. Using multi-frequency Planck radio/microwave maps (30-857 GHz) and Fermi-LAT gamma-ray data (50 MeV-814 GeV), we construct a nonlinear mapping between radio emission and gamma-ray intensity through supervised machine learning. Our models achieve high predictive accuracy (R^2 > 0.90 in the 0.1-10 GeV range), demonstrating that multi-frequency radio observations encode sufficient information to reconstruct both spatial morphology and spectral properties of diffuse gamma-ray emission. By analyzing model performance across different frequency bands and spatial regions, we identify high-frequency radio bands as the dominant predictor, providing direct empirical support for the hadronic origin of Galactic 0.1-10 GeV gamma rays, while low-frequency radio bands for the leptonic origin above 10 GeV. Residual maps reveal coherent large-scale structures, including Loop I and III, highlighting regions where standard interstellar emission models are incomplete or biased. Compared with the GALPROP model, our machine learning approach yields a higher R^2=0.95 and lower mean absolute relative error (14.7%) in the inner Galactic disk and the Galactic center region. Our results illustrate that machine learning serves as a physically interpretable tool for multi-messenger astrophysics, providing a data-driven baseline for separating non-standard emission components and deriving new constraints on cosmic-ray propagation and interstellar medium structure.

Synchrotron emission in high-energy transients is produced by relativistic electrons accelerated by shocks. As high-energy transients are often unresolved even on angular scales probed by very long baseline interferometry, it is difficult to obtain a full picture of the ejecta and circumstellar medium (CSM) properties that are probed by the radio synchrotron emission. Radio spectra of high-energy transients frequently show optically thick slopes shallower than the standard $F_\nu \propto \nu^{5/2}$ expected from synchrotron self-absorption (SSA) models, or broader spectra near the self-absorption frequency. Such deviations are often interpreted phenomenologically, without providing clear insights into the structure of the emitting region. We show how information on the homogeneity and symmetry of the emitting region can be directly inferred from SSA spectra, even when the source is unresolved. We discuss the circumstances under which inhomogeneities in the emitting region can change the spectrum below the self-absorption frequency, causing it to follow a different slope. We examine which parameters can be constrained from observations and which remain degenerate. We apply this method to the stripped-envelope supernova (SN) 2016coi and to the fast blue optical transient (FBOT) AT2018cow, showing that SSA spectra constrain the degree of inhomogeneity in these systems, providing strong evidence for inhomogeneities in the emitting region in the SN 2016coi, and asymmetry in the case of AT2018cow, and we infer the characteristics of the emitting region. When well sampled spectra are available, our method can be applied as a general, model-independent, inference method. This approach can be used to constrain inhomogeneities in a variety of unresolved high-energy astrophysical transients, including SNe, FBOTs, tidal disruption events and gamma-ray bursts.

Many details of the gas accretion phase during giant planet formation remain untested. We present new 0.2$\unicode{x2013}$0.7 $\mu$m UV-through-red optical imaging of the young, wide-orbit planetary-mass companion SR 12 c from the Wide Field Camera 3 (WFC3) instrument on board the Hubble Space Telescope. SR 12 c exhibits strong accretion-related continuum excess blueward of $\sim$5000 $\unicode{x212B}$ and clear signs of the Balmer jump at 3646 $\unicode{x212B}$. We derive a total accretion luminosity of 1.65 $\pm$ $0.19 \times 10^{-5} L_{\odot}$ and a mass accretion rate of 8 $\pm$ $2\times 10^{-12}$ M$_{\odot}$ yr$^{-1}$. Based on its mass and age, SR 12 c will not grow by an appreciable amount at its current accretion rate; it is at the end stages of assembly. No accretion variability is evident between the two epochs of the WFC3 observations spanning a month-long baseline, but the H$\alpha$ emission line strength decreases by 90% compared to the reported flux from five years earlier. Combined with previous observations of SR 12 c, we assemble one of the most complete spectral energy distributions of a young giant planet to date, spanning the UV through sub-mm wavelengths (0.2$\unicode{x2013}$880 $\mu$m). This adds SR 12 c to the small yet growing sample of planets with detailed accretion and disk constraints, which together are beginning to establish the diversity of timescales and physical processes governing the formation of giant planets.

The masses of galaxy clusters carry cosmological and astrophysical information. We develop a simulation-based inference pipeline to infer cluster masses from full projected phase-space information of member and interloper galaxies. Our method combines a permutation-invariant Deep Sets architecture with neural posterior estimation using normalizing flows, enabling the recovery of expressive posterior distributions. We train the model to predict residual corrections to the classical $M$–$\sigma$ relation, thus explicitly isolating information beyond velocity dispersion. Using the Uchuu-UniverseMachine simulation, we evaluate the method under both idealized (interloper-free) and realistic (cylindrical) observational setups. In the idealized case, our model reduces the scatter in mass estimates to as low as $\sim 0.1$ dex, representing a twofold improvement over the traditional $M$–$\sigma$ relation. In the cylindrical setup, we achieve comparable performance at the high-mass end ($> 10^{14.5}\,M_\odot/h$), demonstrating robustness against interloper contamination. We demonstrate that set-based simulation-driven inference provides a powerful and flexible framework for galaxy cluster mass estimation, enabling improved accuracy and reliable uncertainty characterization for upcoming large-scale surveys. Our model saturates the kinematic information content and thus suggests a baseline for future studies.

The gravitational wave background (GWB) produced by extreme-mass-ratio inspirals (EMRIs) serves as a powerful tool for probing the astrophysical and dynamical processes in galactic centers. EMRI systems are a primary target for the space-based detector LISA due to their long-lived signals and high signal-to-noise ratios. This study explores the statistical properties of the GWB from EMRI, focusing on the calculation methods for the GWB, the astrophysical distribution of EMRI sources, and the influence of key parameters, including the spin of supermassive black holes (SMBHs) and the masses of compact objects (COs). By analyzing these factors, we determine the distribution range of the characteristic strain of the GWB from EMRIs. We find that the final eccentricity distributions appear to have negligible effect on the intensity of the GWB due to rapid circularization before they become detectable and the spin of the SMBH enhances the GW characteristic strain by approximately 1$\%$ compared to cases without spin effects. The masses of COs can also significantly affect the characteristic strain of the GWB from EMRIs, with Black Hole (BH) as CO producing a GW signal intensity that is approximately one order of magnitude higher compared to cases where Neutron Star (NS) or White Dwarf (WD) are the COs.

Fast radio bursts (FRBs) are millisecond-duration radio flashes of unknown origin, observationally classified into repeating and apparently nonrepeating (one-off) populations. In this work, we use a statistical population approach to investigate the redshift evolution of one-off FRBs. We compare a pure star formation history (SFH) tracing model, phenomenological delayed models, physically motivated delayed models that correspond to binary neutron star related and neutron star age-window channels, and mixture models which is obtained when two physically motivated models are normalized separately and then combined as a weighted mixture. The samples in CHIME/FRB Catalog do not support an intrinsic event rate density that directly follows the SFH. The preferred model is mixture model corresponds to an effective mean delay time of $\bar{\tau}=1.426^{+0.032}_{-0.035}~\mathrm{Gyr}$. These results suggest that the current data may naturally explained by delayed, possibly multi-channel progenitor evolution for the one-off FRBs.

The $\Lambda$CDM cosmological model has been successful in explaining many astronomical observations. However, recent observations increasingly point to deviations from the standard $\Lambda$CDM framework. Among these, one of the most significant discrepancies is the \textit{Hubble tension}, which refers to the difference in values obtained for the Hubble constant $H_0$ from high-redshift measurement and local observation. To address this issue, numerous cosmological models and methodological approaches have been proposed. This review offers a concise overview of recent progress in resolving the Hubble tension. The combination of Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillations (BAO) and uncalibrated Type Ia supernovae data yields a value for $H_0$ that is significantly higher than the $\Lambda$CDM predication based on early-universe probes, even without incorporating local distance ladder constraints. This result indicates that the origin of the Hubble tension lies in new physics at low redshifts. Our findings suggest that although many unresolved systematics persist in current observations, they are insufficient to account for the magnitude of the current Hubble tension. This implies the likely existence of new physical mechanisms that have yet to be discovered.

Magnetic reconnection plays a central role in solar flares and coronal mass ejections. Identifying where reconnection is physically active within coronal magnetic field models is a key part of magnetic field analysis. We investigate the field line slippage rate as a physics-weighted proxy for three-dimensional reconnection in nonlinear force-free field (NLFFF) extrapolations. The slippage rate measures the instantaneous deviation of magnetic field lines from ideal evolution, due to non-ideal terms in Ohm's law, providing a direct link between magnetic geometry and reconnection physics. For NLFFFs, we show that the resistivity-induced slippage rate is governed by cross-field gradients of the field-aligned twist, thus establishing a clear connection between current structure and reconnection signatures. We further examine its relationship to the squashing factor $Q$, used to identify quasi-separatrix layers (QSLs). By deriving a scaling estimate, we demonstrate that strong magnetic squashing amplifies slippage only insofar as it produces small transverse length scales; large values of $Q$ alone do not guarantee significant reconnection. We apply this framework to a sequence of NLFFF extrapolations of NOAA active region 11158 spanning the X2.2 flare of 15 February 2011. The slippage rate reveals enhanced reconnection signatures associated with distinct phases of the active region's evolution. In comparison with the squashing factor, we show that the field line slippage rate provides a physics-weighted complement to QSL analysis, distinguishing between regions that are geometrically favourable for reconnection and those where reconnection is physically significant.

SN 2018bsz is the closest known stripped superluminous supernova (SLSN-I) to date, making it an ideal laboratory for investigating the physical mechanisms powering this class of extreme explosions. We present a multi-epoch X-ray spectroscopic study of SN 2018bsz based on four Chandra observations followed by one XMM observation, spanning 87 to 1253 days after explosion. The source is detected at all Chandra epochs and is also tentatively detected in the late XMM observation, although more uncertain due to nearby contaminating sources. Regardless of the XMM detection, this makes SN 2018bsz the second X-ray detected SLSN-I and the third X-ray detected SLSN overall. We explore potential power sources for the observed X-ray emission and find that a millisecond (ms) magnetar central engine underpredict most of the observed X-ray luminosities and fail to reproduce the relatively flat light curve. Accounting for ejecta absorption further increases the discrepancy. While asymmetries and magnetar-driven ionization could reduce the effective absorption, ionization breakout is expected years after our observational window. Instead, the observations are more readily explained by early-time interaction between the ejecta and the circumstellar medium, while the magnetar emission is absorbed by the ejecta. This scenario is supported by the flat temporal evolution, previous optical results, and inferred mass-loss rates which resemble those of stripped SNe that later evolve into interacting systems. Our results thus favor the scenario where SN 2018bsz is part of a distinct group of SLSN-I, where interaction is crucial for the strong emission.

The accretion of matter on a massive white dwarf (WD) can lead to repeated nuclear explosions on its surface over a timescale of years to decades. The seventh explosion of the recurrent symbiotic nova RS Ophiuchi (RS Oph) was recorded on August 8, 2021. In this paper, we examine its early evolution, from 9 hours before its optical maximum until day 42. We achieved our goal by modeling the spectral energy distribution (SED) using optical spectroscopy and simultaneous $BVR_{\rm C}I_{\rm C}$ photometry, supplemented by $JHKL$ photometry and ultraviolet spectroscopy from previous explosions in 2006 and 1985. Our SED models revealed an early stage of development of the ejecta bipolar structure, consisting of a flared, density-enhanced equatorial disk and low-density regions in bipolar directions. The comparability of the internal shocks' luminosity in the equatorial outflow, inferred from our model parameters, with the luminosity of the warm WD pseudophotosphere during its presence in the spectrum (until $\sim$day 42) confirmed that a significant part of its radiation originates from reprocessed shock emission. We explain the formation and evolution of the bipolar ejecta structure during RS Oph explosions by the rotation of the accreting WD. Such an ejecta structure provides a natural framework for the generation of strong internal shocks and thus $\gamma$-ray emission inside the ejecta.

The collision of comet SL9 with Jupiter in 1994 changed the chemical composition of the Jovian stratosphere for decades. New molecules were detected minutes after the impacts (HCN, CO, CS, etc.) and some are still present today. They were deposited in the stratosphere at pressures lower than 0.1 mbar and were most probably formed by shock-induced chemistry recombining Jovian and cometary material. However, the question of the origin of these molecules is still not completely understood. One way to address this open question is to determine the isotopic composition of the new molecules. Isotopic ratios have long been measured in the Solar System. They present a variety of values depending on the object or the molecule and therefore trace different reservoirs of material. Derivations of carbon and nitrogen isotopic ratios in HCN four years after SL9 showed atypical depletions in the heavier isotopes that had never been observed before in the Solar System. These results suggested an unusual cometary composition or an unknown fractionation mechanism in the hot and shocked air parcels. We aim to measure carbon and nitrogen isotopic ratios in HCN to shed light on the puzzling results of 1998. With Atacama Large Millimeter/submillimeter Array data from 2017 and radiative-transfer calculations, we derived the abundance of two HCN isotopologues, H13CN and HC15N, at pressures probed from 0.03 to 1.8 mbar. We find 12C/13C = $73\pm5$ and 14N/15N = $245\pm29$, respectively (0.76-0.87) and (0.80-1.00) times the terrestrial references, and (0.69-0.87) and (0.42-0.70) times the solar-Jovian bulk values. In contrast to the strong depletions reported in 1998, our values are instead compatible with an enrichment in the heavier isotopes relative to the Jovian bulk. We interpret these enrichments as the direct signature of the cometary contribution in HCN and/or as 23 years of chemical evolution.

Flare ribbons often display fine structures along their fronts that are commonly interpreted as signatures of intermittent reconnection dynamics including plasmoid formation in the flare current sheet. We introduce an automated method that detects and tracks the spiral- and wave-like imprints of these structures and as a proof of concept apply it to maps of magnetic field-line length from a high-resolution 3D eruptive-flare simulation. The workflow applies the correlation-dimension method, density-based clustering, and a minimum-area ellipse fit to summarise each feature. We show that across the simulated flare, the detected spirals remain locked to the ribbon's outward motion while drifting coherently along the ribbon. The two ribbons show opposite along-ribbon drift and motion away from their hooks in accordance with theoretical expectations, with instantaneous speeds of 10-800 km s^1, all well below the local surface Alfven speed. Occurrence, lifetimes, and mean magnetic flux of the features peak during the impulsive phase. The distribution of per-spiral mean unsigned flux shows a scale-free tail above roughly 6x10^18 Mx with a power-law exponent near 3.4. Together, these results show that bursty, plasmoid-mediated flare reconnection leaves a clear, measurable signature on the flare ribbons. The method provides a practical surface diagnostic of ribbon fine structure that can potentially be used to inform our understanding of three-dimensional magnetic reconnection in the flare current sheet.

Titan possesses a thick N$_2$-CH$_4$ atmosphere that makes it difficult to study its surface spectroscopically. The chemical composition of the solid surface of Titan thus remains very uncertain. By leveraging JWST's high sensitivity and large spectral coverage, we searched for any signature from Titan's surface in the broad and less explored 5-$\mu$m atmospheric window. We also investigated the JWST spectrum of Pluto which has a thin Titan-like atmosphere. We made selections of JWST NIRSpec and MIRI spectra around Titan's disk center and compared the NIRSpec average spectrum with a radiative transfer model including gas and haze opacity. We detected an unidentified absorption in both NIRSpec and MIRI spectra of Titan centered at 5.113 $\mu$m (1956 cm$^{-1}$) and 6-7% deep. The width of the feature is 0.024$\pm$0.0008 $\mu$m (9.2$\pm$0.3 cm$^{-1}$) in the NIRSpec spectrum recorded on the trailing side and is possibly 25% narrower in the MIRI spectrum of the leading side. This absorption most likely originates from the surface. We could not identify this signature among published laboratory spectra of ices relevant to Titan's atmospheric compounds but present a few plausible candidates. A 4-5% deep absorption is also present in the MIRI spectrum of Pluto but is about 3 times broader than on Titan's trailing side.

The accuracy of stellar population properties inferred through spectral energy distribution fitting hinges on the reliability of the underlying spectral models. Binary interactions are fundamental for massive star evolution, and ignoring their spectral contribution can lead to incorrect results. We use the POSYDON binary population synthesis code to generate spectral models of stellar populations that include binaries at solar metallicity. Our framework incorporates a collection of spectral libraries that is designed to address key outcomes of binary stellar evolution like Wolf-Rayet stars, stripped helium stars, and a treatment for stellar mergers. Our models confirm previous results showing that the inclusion of binary interactions has a significant effect on the UV and ionizing regime of the integrated spectrum. In particular we find that Wolf-Rayet stars and other massive stars dominate the production of ionizing radiation at earlier times, but after $\simeq$16 Myr stripped stars produced through mass transfer begin to dominate. Furthermore, we show that the production of ionizing He II photons is especially sensitive to the underlying population of stripped stars. While our results currently focus on high-mass stars ($\ge4~M_{\odot}$) at Solar metallicity, they provide the framework for binary spectral synthesis across a range of metallicities and masses and lay the foundation for calculations of the emergent emission-line spectra in the UV, optical, and IR regimes. We make the spectral models from this work publicly available for use in a format that can be integrated into fitting codes.

Following our earlier work on studying 3-dimensional nuclear recoil response to Galactic Weakly Interacting Massive Particles (WIMPs) in directional direct Dark Matter detectors, in this paper, we simulate 3-D coherent elastic neutrino-nucleus scattering (CEvNS) events induced by Solar B-8 neutrinos. Our numerical results show that, in contrast to the approximately fixed patterns of the WIMP-induced signals, the characteristic ring-like angular distributions of the nuclear recoil flux/energy of CEvNS events show clearly annual variations along the trajectories of the moving direction of incident Solar neutrinos in different celestial coordinate systems without experimentally distinguishable target dependence.

One effective way of learning about the production and properties of dark matter in the early universe is by extracting information about the primordial dark-matter phase-space distribution from the matter power spectrum. Several years ago a simple empirical formula was introduced which successfully reproduces most of the salient features of the primordial dark-matter phase-space distribution from the matter power spectrum – even in situations in which this distribution is non-thermal, multi-modal, or exhibits other complicated features. Continuing this line of research, we investigate the extent to which machine-learning techniques can improve upon this analytic approach. Interestingly, we find that a one-dimensional convolutional neural network not only succeeds in reconstructing the dark-matter phase-space distribution with greater accuracy, but can also be applied to a broader range of matter power spectra.

Cosmic voids and superclusters are among the largest structures in the Universe and provide complementary probes of the growth of large-scale structure and the underlying gravitational field. We present a comprehensive catalogue of the main local Universe cosmic web voids and knots, using the updated CosmicFlows-4++ Zone of Avoidance (CF4++ZOA) catalogue out to redshift z=0.1. To do this we use the V-web algorithm which provides a dynamically motivated map of the the local cosmic web, tracing the major expanding and converging regions of the local Universe. The robustness of the catalogue is assessed using the ensemble of Hamiltonian Monte Carlo realizations of the CF4++ZOA reconstruction. We additionally remove any structure that exceeds the survey boundaries and present catalogues of 37 voids and 42 knot regions within the reconstructed survey volume. The identified emptying regions (voids) have effective radii ranging from 13 to 38 h-1Mpc. The high density converging regions (knots) have volumes ranging from 10^4 to 3.3x10^5 h-3 Mpc^3.

The relatively recent revelation of the high occurrence rate of binary interactions, especially in intermediate- and high-mass systems, has prompted multiple investigations into their asteroseismic imprints. The near-core region just outside the convective cores of mass-accreting early-type main-sequence stars in binaries has been theorised to be sensitive to assumptions about mixing (notably semiconvection) and accretion physics. In turn, the predicted asteroseismic properties depend strongly on the physical properties of this near-core region. We explore how robust the previously identified asteroseismic imprints of mass accretion are to changes in semiconvective mixing. Using one-dimensional stellar structure and evolution models, this parameter study shows the dominant effect of convective boundary mixing on rejuvenation and the post-accretion asteroseismic properties. The recovered seismic imprint, largely robust to variations in semiconvective mixing efficiency, changes drastically when convective boundary mixing is not included in the models. We find that the post-accretion thermal relaxation is key in determining the final near-core structure and the asteroseismic imprint of accretion. We reaffirm the potential of Fourier transforms of period spacing patterns to quantify the effects of different near-core mixing and accretion-rate assumptions on asteroseismic signals. Overall, this work highlights the sensitivity of the asteroseismic imprint of accretion not only on stellar structure and evolution modelling assumptions, but also on the accretion physics. The logical next step is to arrive at a more general picture of the asteroseismic imprint of mass transfer by exploring its properties in a multi-dimensional parameter study including single- and binary-star assumptions.

We present a simple universal computational algorithm for computing compositional phase diagrams of rocks and their melts at given temperature and pressure. It makes use of the mathematical concept of the convex hull of a set of points in the space spanned by the composition and the Gibbs free energy. All the complexities of determining the stability or separation of phases, the localization and orientation of tie lines, as well as the determination of characteristic points, curves and surfaces such as the solidus, liquidus, solvus, and the eutectic/peritectic points etc, are taken care of by the algorithm that computes the convex hull, supplemented with an algorithm to physically classify the resulting simplices. For the convex hull computation, the publicly available Qhull package can be used, which is available in SciPy. This makes this method accessible and intuitive for a broad set of scientific and educational applications. Although the method is not practical for systems of a large number of components, it is remarkably stable and efficient for systems of up to four. We present our implementation of the method as a publicly available Python package.

We analyze Maxwell theory coupled to an axion domain wall as a spectral boundary value problem. We find that a finite-width axion domain wall generically supports a localized normalizable chiral electromagnetic mode with linear, gapless dispersion. This mode arises from helicity-dependent coupling sourced by the axion gradient: one polarization experiences an effective attractive potential and forms a bound state, while the opposite polarization is repelled. The existence of this chiral surface photon is robust over a wide regime of wall structures and axion masses. Our result shows that axion domain walls generically support a localized chiral photon that has been missed in previous analyses.

We study the ability of the upcoming Laser Interferometer Space Antenna (LISA) to constrain gas torques acting on extreme-mass-ratio inspirals (EMRIs) when these are embedded in accretion disks, using recently developed relativistic models for the binary-disk interaction. Using a fully Bayesian setup, we find that, contrary to previous forecasts based on Newtonian results, these observations can provide simultaneous estimates of the disk surface density and the accretion rate (or, equivalently, its total luminosity) without the need for an electromagnetic counterpart. Our analysis also indicates that simpler measurement constraints based on the linear-signal (Fisher matrix) approximation are not valid for these systems. For typical EMRI observations, the torque amplitude can be constrained to within  10%, strengthening the prospect of probing accretion physics at (sub)microparsec scales, deep in the strong-field gravity regime and complementing electromagnetic observations. This also strengthens LISA's ability to help answering questions such as how massive black holes grow and coevolve with their host galaxies and, by helping to identify the EMRI's host galaxy through cross-correlation with AGN catalogues, to improve the use of these sources as (dark) sirens for cosmology.

We propose a generic, largely inflaton-potential-independent mechanism in which a light axion spectator, initialized near the hilltop of its potential, reshapes inflationary observables through purely gravitational multi-field dynamics. During inflation the axion is frozen and the background follows an effectively single-field trajectory. After inflation ends, the axion rolls, inducing a turn in field space and transient tachyonic phases of the isocurvature mode. The resulting “tachyonic encore” occurs entirely on super-horizon scales. These phases generate a nearly scale-invariant enhancement of the curvature power spectrum, suppressing the tensor-to-scalar ratio and shifting the scalar tilt to a weighted combination of adiabatic and entropic tilts at horizon crossing. We show that these effects can reconcile otherwise disfavored inflaton potentials with current CMB constraints. The same dynamics predict local non-Gaussianity, $f_{\rm NL}^{\rm loc.}\sim \mathcal{O}(1)$, within reach of upcoming surveys.

While vertically-sheared stratified flows have been studied extensively, their horizontally-sheared counterparts have received considerably less attention. Yet, horizontal shear instabilities remain active even when the mean Richardson number is large or even formally infinite, and can drive turbulence in strongly stratified (low Froude number) flows at sufficiently high Reynolds number. In this work, we combine linear theory with direct numerical simulations to investigate two pathways to turbulence in low Froude / high Reynolds number horizontally-sheared flow with no initial vertical shear. In the first pathway, vertical shear emerges directly from vertically-modulated eigenmodes of the primary horizontal shear instability, and becomes unstable to secondary small-scale Kelvin-Helmholtz (KH) instabilities on the buoyancy scale at sufficiently large buoyancy Reynolds number $Re_b$. In the second pathway, a vertically-invariant eigenmode of the primary horizontal shear instability initially dominates, causing the background flow to evolve nonlinearly into a long-lived time-dependent two-dimensional (columnar) vortical flow. The vortices are subsequently unstable to secondary three-dimensional hyperbolic instabilities from which vertical shear emerges, which is finally unstable to tertiary small-scale KH instabilities on the buoyancy scale at sufficiently large $Re_b$. This shows that the emergence of vertical shear driving small-scale KH instabilities is an inevitable by-product of horizontal shear instabilities in strongly stratified flows at sufficiently large $Re_b$. However, we also find that the two pathways excite different ranges of vertical scales, which results in different peak mixing efficiencies.

Grand unified theories (GUTs) predict the overproduction of magnetic monopoles, leading to the so-called monopole problem, which is often addressed by cosmological inflation that dilutes their abundance. However, if inflation occurs before the GUT symmetry breaking, monopoles are produced afterwards and the problem persists. This motivates the exploration of alternative mechanisms. We propose a simple solution based on the Langacker–Pi mechanism within an $SU(5)$ GUT framework with symmetry breaking into its special subgroup. In particular, after the gauge symmetry is broken to the Standard Model (SM) gauge group $SU(3)_C\times SU(2)_L\times U(1)_Y$ by the vacuum expectation value of an adjoint scalar, it is further broken to $SO(3)_C\times SO(2)_L$. This structure is naturally realized by introducing a symmetric tensor scalar, a singlet scalar, and multiple singlet fermions. During this intermediate phase, the monopoles become connected to antimonopoles by cosmic strings, which enhances their pair annihilation and reduces their abundance. Subsequently, the symmetry is restored to the SM gauge group. The restoration transition can be either first-order or second-order, depending on the model parameters. In the case of a first-order phase transition, a stochastic gravitational-wave (GW) signal is generated. For a certain region of the parameter space, the resulting signal can lie within the sensitivity of future GW experiments.

Geoneutrinos, electron antineutrinos produced by radioactive decays of heat-producing elements (HPEs) within the Earth, provide unique insights into Earth's interior and heat budget since their first detection in 2005 by KamLAND. Conventional geoneutrino detectors currently provide integrated global information and lack the capability to spatially resolve structures deep within the Earth. Here, we evaluate the ability of angular-sensitive geoneutrino detectors to distinguish between homogeneous and heterogeneous mantle models, focusing on Large Low Shear Velocity Provinces (LLSVPs). Our results show that LLSVPs enriched in Th and U yield a distinct flux of geoneutrinos with distinctive angular patterns. An oceanic site above the Pacific LLSVP is considered a particularly favorable detector location. The Ocean Bottom Detector (OBD) project aims to leverage this spatial resolving advantage by deploying a kiloton-scale liquid scintillator detector directly on the ocean floor, enabling unprecedented sensitivity for mantle geoneutrino detection. These findings demonstrate the critical role of combining geophysical and geochemical data to guide detector site selection, ultimately improving constraints on Earth's internal heat and the HPE distribution.

We constrain the cosmological constant $\Lambda$ using astrometric and spectroscopic observations of the S2, S1, and S14 stars orbiting Sgr A$^*$. The stellar motion is modelled by numerically integrating timelike geodesics in Schwarzschild-de Sitter spacetime, including relativistic redshift and time-delay corrections. Orbital and spacetime parameters are inferred using a Bayesian MCMC analysis. The resulting posterior distributions place upper bounds on the magnitude of $\Lambda$ at the Galactic Center (GC). Combining the independent constraints from the S2, S1, and S14 orbits yields upper bounds of $\Lambda \lesssim 6.9\times10^{-48} \mathrm{m}^{-2}$ at 68\% credibility and $\Lambda \lesssim 1.0\times10^{-38} \mathrm{m}^{-2}$ at 95\% credibility.

We study supercooled first-order phase transitions in a supersymmetric hidden sector with a spontaneously broken $U(1)_X$, focusing on the frequency range of the Einstein Telescope and Cosmic Explorer. Along the D-flat direction the tree-level quartic vanishes, so the barrier is generated radiatively by soft SUSY-breaking splittings. In the $\overline{\rm DR}$ scheme the gaugino mass $M_{\tilde\lambda}$ sets the barrier depth, while the soft scalar mass $m_0$ stabilizes the broken vacuum. For $M_{\tilde\lambda}/v_X\simeq0.05$–$0.23$, the predicted signal reaches $\Omega_{\rm GW}h^2\sim3\times10^{-10}$ near the percolation boundary. The observable amplitude depends sensitively on the portal coupling $\delta$ through the hidden-to-visible temperature ratio at percolation: for a cold initial hidden sector the signal rises from the ET floor at $\delta=10^{-6}$ to $\Omega_{\rm GW}h^2\simeq7\times10^{-11}$ as the sectors approach thermal contact at $\delta=10^{-4}$, while a hotter initial hidden sector gives a large signal already for weak portal coupling. We follow this evolution with an 11-variable Boltzmann system that separates the cold nucleating exterior from the reheated true-vacuum interior; reheating mainly enters through the energy budget and redshift factors. The same hidden sector can reproduce $\Omega_{\rm CDM}h^2=0.12$ through relativistic dark-quark freeze-out followed by entropy dilution from hidden-Higgs decay, with $m_q\simeq30$–$800\;$keV and $N_{\rm eff}\lesssim{\rm few}\times10^{-5}$.

We study the impact of imposing an infinite-future condition on Gaussian-process (GP) reconstructions of $H(z)$ from 37 cosmic-chronometer measurements. Implementing the asymptotic limit $H(-1)=0$ expected in constant-$w$CDM with $w>-1$ as a pseudo-point at $z=-1$ with tunable uncertainty $\sigma_{-1}$ lowers the GP-inferred Hubble constant from $H_0=68.71\pm6.08$ to $H_0\simeq(64.67$–$65.86)\pm(4.45$–$4.85)~{\rm km\,s^{-1}\,Mpc^{-1}}$. The resulting $H_0$ remains within $\sim0.32$–$0.61\,\sigma$ of the \textit{Planck} $\Lambda$CDM value, while the separation from representative local distance-ladder measurements increases to $\sim1.45$–$1.83\,\sigma$. A scan over $\sigma_{-1}$ shows that the shift is governed by the effective pseudo-point weight, interpolating between the hard-condition and data-dominated limits. Finally, constant-$w$CDM Markov Chain Monte Carlo (MCMC) fits to the same $H(z)$ data with $H_0$ fixed to the GP-inferred values show that $\Omega_m$ is only weakly affected, whereas lower $H_0$ shifts $w$ toward less negative values; allowing curvature broadens constraints and remains consistent with $\Omega_k=0$.

Tomographic weak lensing surveys contain intrinsic symmetries that depend solely on the geometric structure of the Universe. These symmetries can be revealed through null tests and verifying their validity provides constraints on cosmological parameters that govern the background evolution-particularly the redshift dependence of the angular diameter distance. This forms the foundation of the tomographic cosmic shear nulling test introduced in this work. We describe how this test can be implemented, what aspects of cosmology it can constrain, and its specific efficiency in doing so. We also assess its sensitivity to astrophysical effects-such as magnification bias and reduced shear corrections-as well as to observational systematics, including errors in the mean redshift of source bins. Our results show that, in a survey configuration comparable to that of Euclid, this null test can yield complementary constraints on key cosmological parameters such as ${\Omega_{\rm m}, {\rm w}_0}$. However, due to its subdominant constraining power compared to standard 3x2pt analysis and through the identification of a required precision of order $10^{-3}$ on the mean redshift of the bins, we conclude that nulling would better be used as a photometric redshift calibration probe or consistency check. In combination with standard weak lensing and galaxy clustering analysis, it would then offer a promising route to better control systematics and improve the precision of future cosmological measurements.

We present updated observational constraints on the spatially flat $\phi$CDM model, where dark energy is described by a minimally coupled scalar field $\phi$ with an inverse power-law potential $V=V_0 \phi^{-\alpha}$. Using Planck 2018 CMB temperature, polarization (P18), and lensing power spectra (lensing), along with a compilation of non-CMB data including baryon acoustic oscillation, type Ia supernova, Hubble parameter, and growth rate measurements, we constrain $\phi$CDM and $\phi$CDM+$A_L$ models where $A_L$ is the CMB lensing consistency parameter. The scalar field parameter $\alpha$, which governs dark energy dynamics, is more tightly constrained by non-CMB data than by CMB data alone. For the full dataset, we obtain $\alpha = 0.055 \pm 0.041$ in the $\phi$CDM model and $\alpha = 0.095 \pm 0.056$ in the $\phi$CDM+$A_L$ model, mildly favoring evolving dark energy over a cosmological constant by $1.3\sigma$ and $1.7\sigma$. The Hubble constant is $H_0=67.55_{-0.46}^{+0.53}$ km s$^{-1}$ Mpc$^{-1}$ in the $\phi$CDM model, consistent with median statistics and some local determinations, but in tension with other local determinations. The constraints for matter density and clustering amplitude ($\Omega_m = 0.3096 \pm 0.0055$, $\sigma_8 = 0.8013_{-0.0067}^{+0.0077}$) of the flat $\phi$CDM model statistically agree with $\Lambda$CDM model values. Allowing $A_L$ to vary reduces tensions between CMB and non-CMB data, although we find $A_L = 1.105 \pm 0.037$, $2.8\sigma$ higher than unity, consistent with the excess smoothing seen in Planck data. Model comparison using AIC and DIC indicates that the $\phi$CDM model provides a fit comparable to $\Lambda$CDM, with the $\phi$CDM+$A_L$ slightly preferred. Overall, while the $\Lambda$CDM model remains an excellent fit, current data leave open the possibility of mildly evolving quintessence-like dynamical dark energy.

We present optical + near-infrared (NIR) + mid-infrared (MIR) observations of the normal Type Ia supernovae (SN Ia) 2022aaiq and 2024gy in the nebular phase, continuously spanning 0.35-28 microns. Medium-resolution JWST spectroscopy reveals novel narrow ($v_{\mathrm{FWHM}}<1500$ km s$^{-1}$) [Ni II] 1.94 and 6.64 micron cores in both events. The MIR [Ni II] 6.64 micron line exhibits a distinct narrow core atop a broader base, indicating a central enhancement of stable Ni. This structure points to high central densities consistent with a near-Chandrasekhar-mass ($M_{\text{Ch}}$) progenitor or a high-metallicity sub-$M_{\text{Ch}}$ progenitor. From detailed line-profile inversions of SN 2024gy, we derive emissivity profiles for stable iron-group elements (IGEs), radioactive material, and intermediate-mass elements (IMEs), revealing spatially distinct ejecta zones. The [Ni III] 7.35 micron line shows a shallow-to-steep slope transition - a "broken-slope" morphology - that matches predictions for delayed detonation explosions with separated deflagration and detonation ashes. We also reanalyze and compare to archival JWST spectra of SN 2021aefx and the subluminous SN 2022xkq. From the stable Ni luminosities, we infer that SN 2024gy produced  5-10 times more stable Ni mass than SN 2022xkq, favoring a near-$M_{\text{Ch}}$ scenario for SN 2024gy and sub-$M_{\text{Ch}}$ scenario for SN 2022xkq. These results demonstrate that resolved line profiles, now accessible with JWST, provide powerful diagnostics of explosion geometry, central density, and progenitor mass in SN Ia.

The evolved cosmological matter density field is fully determined by the initial matter density field at fixed cosmological parameters. However, the two-dimensional cosmological projected matter density field, relevant for weak-lensing and photometric galaxy studies, is fully determined by the initial projected matter density field only at the linear order. At non-linear order, the entire volume of initial matter contributes. We study a model for the evolved projected density field that is deterministic in the initial projected density fields and probabilistic in the effects of the remaining modes in the initial conditions. We write down predictions for the mean evolved projected field model using Lagrangian perturbation theory. We run a suite of small $N$-body simulations with fixed projected initial conditions and measure the statistical properties of the ensemble of evolved projected fields. Measurements and theory are in good agreement and show that the information on the initial projected fields is exponentially suppressed on non-linear scales. We implement this approach in a likelihood code and use Hamiltonian Monte-Carlo sampling to show that initial fields can be reconstructed even in the presence of non-trival mask features.

In collisionless astrophysical plasmas, turbulence mediates the partitioning of free energy among cascade channels and its dissipation into ion and electron heat. The resulting ion heating is often anisotropic, with ions observed to be preferentially heated perpendicular to the local magnetic field; understanding the mechanisms responsible for this heating is a key step in understanding the evolution of such plasmas. In this paper, we use the framework of quasi-linear theory to compute analytically the heating rates of ions interacting with turbulent, large-scale Alfvénic fluctuations. We show how the imbalance of the turbulence (the difference in energies between Alfvénic fluctuations travelling parallel and antiparallel to the magnetic field) modifies the spatiotemporal spectrum of these fluctuations, allowing the heating mechanism to smoothly transition between stochastic heating in balanced turbulence and cyclotron-resonant heating in imbalanced turbulence. The resultant heating rate is found to have a general form regardless of the level of imbalance, exhibiting a suppression related to the conservation of the ions' magnetic moment at small turbulent amplitudes and recovering previous empirical results in a formal calculation. The results of this work help to consolidate our qualitative understanding of ion heating within astrophysical plasmas, as well as yielding specific quantitative predictions to analyse simulations and observations.

Blazars are a class of active galactic nuclei, supermassive black holes located at the centres of distant galaxies characterised by strong emission across the entire electromagnetic spectrum, from radio waves to gamma rays. Their relativistic jets, closely aligned to the line of sight from Earth, are a rich and complex environment, characterised by the presence of strong magnetic fields over parsec-scale lengths. Owing to their cosmological distance from Earth, these sources serve as ideal targets to probe non-standard gamma-ray propagation. In particular, axion-like particles (ALPs) could be detected through their coupling to photons, which enables ALP-photon conversions in external magnetic fields, leading to distinct signatures in the blazars' gamma-ray spectra. In this work, we estimate the potential of the Cherenkov Telescope Array Observatory (CTAO) to constrain the ALP parameter space by simulating observations of two bright blazars, Mrk 501 and PKS 2155$-$304. We obtain projected $2\sigma$ exclusion regions, demonstrating that CTAO will be able to consistently improve present limits thanks to its greater energy resolution and point-source sensitivity with respect to present ground-based gamma-ray telescopes. In addition to the standard statistical technique based on the likelihood ratio test, we further demonstrate the application of a new method based on machine learning classifiers, which may help in reducing the effect of systematic model-dependent uncertainties in future ALP searches.

We used \texttt{MIST} isochrone fitting and a dedicated grid of stellar evolution models computed with MESA to constrain the ages of the components of the $\gamma$ Persei binary system. While individual stars can be matched to the models at specific metallicities, no joint isochrone solution reproduces both the observed masses and evolutionary states. The stellar evolutionary tracks calculated by MESA reveal a clear evolutionary mismatch. The primary component of the system is in a post-main-sequence phase consistent with the red giant branch or red clump. In contrast, the lighter secondary component lies near the turn-off point of the main sequence or is in the early phase of the subgiant branch. This discrepancy can be overcome by assuming that the $\gamma$ Persei system was born as a triple and the primary component is a rejuvenated star formed through a merger of a close-by pair of main-sequence stars. We show that the merger must have occurred no later than a few hundred megaryears after system formation, and the progenitor masses of the merging stars are restricted by a combination of stars that fall within a narrow band in the $(M_{1,a},M_{1,b})$ plane, corresponding to $M_{1,a}\simeq0.9$–$2.1\,M_\odot$ and $M_{1,b}\simeq2.3$–$2.5\,M_\odot$.

Polarimetric interferometry is a powerful tool for probing both black hole accretion physics and the background spacetime. Current models aimed at explaining the observed multiwavelength flares in Sgr A* often assume hotspots moving on geodesic, Keplerian orbits. In many scenarios, though, a hotspot may instead follow an inspiraling trajectory, potentially transitioning into a plunge toward the black hole. In this work, we present a general framework to simulate the polarized emission from generic equatorial inspiraling hotspots in Kerr spacetime using a parametric four-velocity profile. This parametrization defines a continuous family of flows, ranging from Cunningham's disk model (fixed radius orbits outside the innermost stable circular orbit and plunging motion within the innermost stable circular orbit) to purely radial motion, thereby extending the standard assumptions. Within this framework, we show that inspiral motion produces a distinctive observational signature: a precessing, unwinding evolution of the polarimetric Stokes Q-U looping pattern, in sharp contrast with the closed Q-U loops associated with stable orbits at a fixed radius. We then explore how the morphology of these signatures depends on black hole spin, observer inclination, and magnetic-field configuration. The presented model can be applied to current and near-future interferometric observations of linear polarization, offering a new avenue to probe the physics of matter spiraling inward and the relativistic velocities of plunging plasma.

Extremely high-peaked BL Lac objects - also named extreme blazars - are among the most energetic and persistent extragalactic accelerators in the Universe, defined by a synchrotron emission peaking above $10^{17}$ Hz in X-rays. Such emission is then reprocessed and produces radiation extending deeply into very-high-energy (VHE, energy E>100 GeV) gamma rays. Observations in this energy band - optimally investigated by the Imaging Air-Shower Cherenkov telescopes - are crucial for probing the physical processes that drive their extreme behavior. This study extends our investigation of extreme blazars in the VHE gamma-ray range, providing a second new mini-catalog of sources observed by the MAGIC telescopes. We report on the monitoring of seven targets between 2017 and 2025, including four newly observed sources and three that have been part of long-term observation campaigns, for a total of approximately 338 hours of observations. The analysis of MAGIC data reveals two new VHE detections of extreme blazars, along with three additional sources showing hints of VHE emission. Joint observations of MAGIC and the first Large-Sized Telescope (LST-1) also confirmed a new VHE extreme blazar. Our results are complemented by simultaneous multiwavelength observations in other energy bands, including optical-UV, X-rays, and high-energy gamma rays (100 MeV<E<100 GeV). We confirm typical behavior of extreme blazars, such as a modest variability and a “harder-when-brighter” trend in X-rays across the sample. This new set increases the population of extreme blazars and their broadband analysis confirms the physical properties of these extreme sources.

We study the optical depth to reionization in a cosmological setting that includes backreaction from matter inhomogeneities, using the Buchert averaging formalism. We construct a spacetime model consisting of multiple inhomogeneous domains, hereafter referred to as the backreaction model, characterized by a set of parameters. We first examine how these parameters influence the computation of the optical depth to reionization, $\tau_{reion}$. Next, we carry out a Markov Chain Monte Carlo (MCMC) analysis based on the PantheonPlus+SH0ES Type Ia supernova sample to infer the best-fit values of the model parameters, and then use these to evaluate $\tau_{reion}$. We obtain $\tau_{reion} = 0.0581^{+0.0105}_{-0.0096}$ (68$\%$ confidence limits). This result indicates that, when PantheonPlus+SH0ES data are used to constrain the model parameters, our backreaction model yields a value of $\tau_{reion}$ that aligns more closely with observational estimates than the value predicted by the standard cosmological model. We further demonstrate that the backreaction model leads to a modest reduction of the Hubble tension, while avoiding the need for exotic or non-standard physics.

Measuring the chemical and isotopic compositions of gas giants and brown dwarfs provides insights into their formation pathways and birth environments. 2MASS J0249-0557 c is an L2-type planetary mass companion ($\sim 12 M_{\mathrm{Jup}}$) orbiting a pair of brown dwarfs in the $\beta$ Pic young moving group and vicinity. Its mass places it at the intersection of planets and brown dwarfs, making it an interesting target for constraining formation pathways at the planet-brown-dwarf boundary. Using high-resolution spectroscopic data of the planet acquired with CRIRES+ mounted on VLT, we conduct atmospheric retrieval with the radiative transfer code \texttt{petitRADTRANS} and the nested sampling tool PyMultiNest. We retrieve a C/O ratio of $0.57\pm0.01$, a metallicity of [M/H] = $0.18\pm0.05$, and a $^{12}$CO/$^{13}$CO ratio of $95^{+23}_{-17}$. We also retrieve atmospheric compositions for two benchmark brown dwarfs in the $\beta$ Pic YMG, 2MASSI J0443+0002 and SIPS J2000-7523, using CRIRES+ data and find consistent compositions. Together with 2MASS J0249-0557 c's wide separation from its host, its compositional consistency with benchmark brown dwarfs supports gravitational collapse in a star-like manner as its most likely formation mechanism. These results deliver a homogeneous comparison of three substellar members in the $\beta$ Pic YMG and vicinity. Their solar-like abundances provide a baseline for exoplanet members in the same moving group, such as $\beta$ Pic b, 51 Eri b, and AF Lep b, whose host stellar compositions are difficult to measure. Future comparisons of atmospheric compositions among this moving group offer the potential to distinguish between formation mechanisms for its planetary members.

We investigate the spectrum and polarization of radiation emerging from a subcritical X-ray pulsar using self-consistent radiation-hydrodynamic simulations of an accretion channel in a strong magnetic field. The polarized radiative transfer in the channel above the hot spot is simulated for the two normal modes, taking into account resonant Compton scattering in a strongly magnetized plasma and the effects of vacuum polarization. We show that the deceleration of the accreting matter in the subcritical regime is mainly governed by resonant scattering. Our simulations provide the velocity profiles of the plasma flow and demonstrate that vacuum polarization dominates over plasma birefringence, enhancing both the cyclotron spectral feature and the radiative deceleration of the plasma. We also find that the energy of the cyclotron feature increases with accretion luminosity, indicating a positive correlation consistent with previous observational results and theoretical interpretation.

We present a compressed Gaussian likelihood for the Planck CMB low-$\ell$ E-mode polarization data, constructed from the SRoll2 likelihood which provides the tightest constraint on the reionization optical depth $\tau$ to date. The non-Gaussian form of CMB low-$\ell$ TT and EE likelihoods makes them incompatible with Fisher matrix analyses that require an analytic Gaussian $\chi^2$, such as the Fisher-bias formalism and Fisher forecasts. We show that the $\chi^2$ of an offset log-normal likelihood takes a Gaussian form in the log-transformed power spectrum amplitudes, and can therefore serve as a proxy for the true Gaussian likelihood of this variable in Fisher matrix analyses, without any explicit change of variables. Building on this, we compress the SRoll2 likelihood into a small number of piecewise offset log-normal functions and validate it against the full SRoll2 likelihood via MCMC combined with Planck and ACT DR6 data, finding excellent agreement across all $\Lambda$CDM parameters and in extended cosmological models. We further demonstrate that Fisher matrix uncertainty estimates from our compressed likelihood agree well with the full MCMC posteriors. We release our compressed likelihood planck-gaussian-lowl, a lightweight Python package incorporating the compressed low-$\ell$ TT likelihood from previous work, allowing a straightforward incorporation of the Planck CMB low-$\ell$ data into any Gaussian-likelihood-based analysis. The package is publicly available at \href{ this https URL }{ this http URL }.

Although measurements of the Cosmic Microwave Background (CMB) are consistent with a nearly scale-invariant primordial spectrum of adiabatic perturbations, in which the energy densities of different components (radiation, baryons, and dark matter) fluctuate proportionally, there could also exist isocurvature perturbations, in which density fluctuations of the individual components differ from the adiabatic mode. Cold dark matter isocurvature (CDI) perturbations with a variety of spectral tilts generated in pre-inflationary axion models provide one such example. In this article, we present the most updated constraints on these axion CDI perturbations using the latest CMB anisotropy measurements from Planck, the Atacama Cosmology Telescope (ACT), and the South Pole Telescope (SPT). We study both fixed spectral indices with values ranging from red- to blue-tilted spectra as well as the case with a free index. We find that the constraint on the spectral index gets moderately improved with the combined datasets compared to Planck alone, while the bounds on the isocurvature amplitudes for the fixed spectral indices we consider do not get tighter. We also discuss the theoretical implications of our constraints, in particular for models giving rise to blue-tilted spectra.

Cassini states correspond to the equilibria of the spin axis of a body when its orbit is perturbed. They were initially described for planetary satellites, but the spin axes of black hole binaries also present this kind of equilibria. In previous works, Cassini states were reported as spin-orbit resonances, but actually the spin of black-hole binaries is in circulation and there is no resonant motion. Here we provide a general description of the spin dynamics of black hole binary systems based on a Hamiltonian formalism. In absence of dissipation the problem is integrable and it is easy to identify all possible trajectories for the spin for a given value of the total angular momentum. As the system collapses due to radiation reaction, the Cassini states are shifted to different positions, which modifies the dynamics around them. This is why the final spin distribution may differ from the initial one. Our method provides a simple way of predicting the distribution of the spin of black hole binaries at the end of the inspiral phase.

In this paper, we present 52 new numerical-relativity (NR) simulations of black-hole-neutron-star merger (BHNS) mergers and employ the data to inform TEOBResumS-Dalí: a multipolar effective-one-body model also including precession and eccentricity. Our simulations target quasicircular mergers and the parameter space region characterized by significant tidal disruption of the star. Convergent gravitational waveforms are produced with a detailed error budget after extensive numerical tests. We study in detail the multipolar amplitude hierarchy and identify a characteristic tidal signature in the $(\ell,m)=(2,0)$, and $(3,0)$ modes. We also develop new NR-informed models for the remnant black hole and for the recoil velocity. The numerical data is then used to inform next-to-quasicircular corrections and the ringdown of TEOBResumS-Dalí for BHNS. We show an overall order of magnitude improvement in the waveform's amplitude at merger and more consistent multipoles over our older TEOBResumS-GIOTTO for BHNS. TEOBResumS-Dalí is further validated with a new 12 orbit precessing simulation, showing phase and relative amplitude differences below $\sim 0.5$ (rad) throughout the inspiral. The computed mismatches including all the modes lie at the one percent level for low inclinations. Finally, we demonstrate for the first time that TEOBResumS-Dalí can produce robust waveforms with both eccentricity and precession, and use the model to identify the most urgent BHNS to simulate for waveform development. Our new numerical data are publicly released as part of the CoRe database.

We propose a boomerang mechanism for the explanation of the excess radio background detected by ARCADE 2. In an early stage of the Universe, at a temperature $T$ in the range $\sim 0.1\,{\rm keV}$–$1\,{\rm MeV}$, a fraction of relic neutrinos is resonantly converted into dark neutrinos by mixing induced by a pre-existing lepton asymmetry. Dark neutrinos decay much later into a dark-standard photon state and a dark fermion, with a lifetime longer than the age of the Universe, as required by a solution to the excess radio background. This scenario circumvents the upper bound on the neutrino magnetic moment but still implies a testable lower bound.

In this paper, we compute the inflationary trispectrum of primordial gauge fields generated through the scalar and tensor exchanges in models with spectator $U(1)$ gauge fields which are kinetically coupled to the inflaton. Focusing on the connected four-point autocorrelation function of gauge fields, we derive exact analytical expressions for the full trispectrum of both electric and magnetic fields using the in-in formalism and cosmological diagrammatic rules, and explore their respective contributions in specific momentum configurations. For the scalar exchange, we find that the trispectrum signal in the equisided configuration grows with the exchange momentum and reaches its maximum in the flattened limit. However, in the counter collinear limit, we show that the non-linearity parameter associated with the trispectrum scales quadratically with the corresponding parameter of the cross-correlation bispectrum of magnetic fields and curvature perturbations, thereby establishing a hierarchical relation between the higher- and lower-order correlation functions. For the tensor exchange, the trispectrum displays a richer angular dependence, reflecting the sensitivity to the orientation of the momentum quadrilateral with respect to the tensor polarisation, producing characteristic angular modulations in the trispectrum. Detecting such angular signatures in future high-precision cosmological observations would provide a novel window into tensor-mediated interactions in the early universe.

Scattering black holes spin up and gain mass through the re-absorption of orbital angular momentum and energy radiated in gravitational waves during their encounter. In this work, we perform a series of numerical relativity simulations to investigate the spin-up and mass-gain for equal-mass black holes with a wide range of equal initial spins, $\chi_{\rm i}\in[-0.7,0.7]$, aligned (or anti-aligned) to the orbital angular momentum. We also consider a variety of initial momenta. Furthermore, we explore a range of incident angles and identify the threshold between scattering and merging configurations. The spin-up and mass-gain are typically largest in systems with incident angles close to the threshold value, large momenta, and negative (i.e. anti-aligned) initial spins. When evaluated at the threshold angle, we find that the spin-up decreases linearly with initial spin. Intriguingly, systems with initial spin $\chi_{\rm i}=0.7$ sometimes experience a spin-down, in spite of an increase in the black-hole angular momentum, due to a corresponding gain in the black-hole mass. Across the simulation suite, we find a maximum spin-up of $0.3$ and a maximum increase in the black-hole mass of $15\%$.

A dense neutrino plasma can exhibit collective flavor evolution caused by neutrino–neutrino refraction. Recently, a new class of exact nonlinear inhomogeneous solutions was discovered: single-wave (SW) solutions of the fast flavor system. The key property is that the flavor occupation numbers remain homogeneous, whereas the field of flavor coherence varies spatially with a single wave vector. The equations of motion for this structure resemble those of a collection of classical spins, in analogy with the homogeneous slow and fast flavor cases. In contrast, the SW system is not integrable (it does not possess Gaudin invariants) so that, while two-beam pendulum solutions are inevitable, they do not extend to a multi-angle system. We develop a taxonomy of all known nonlinear collective flavor solutions, explaining the overlap between categories and their differences.

Flavor instabilities in dense neutrino media trigger exponential growth of flavor waves, yet their nonlinear saturation remains poorly understood. We examine a simple proxy for this effect in the form of a single-wave solution of an axially symmetric fast flavor system. When the angular crossing is shallow and the growth rate of the instability correspondingly small, the flavor wave primarily affects resonant neutrinos that move in phase with it. The evolution of these resonant neutrinos becomes periodic, undergoing cycles of full flavor reversal. They feed power into the unstable wave, and subsequently return to their initial state, draining power back out. This new flavor pendulum captures the dynamics of weak, nearly monochromatic fast flavor instabilities. Since weakly unstable distributions always exhibit a narrow range of unstable wavenumbers, our model likely describes the earliest development of a flavor instability when it first appears. When the instability is not weak, the linear phase of a single-wave excitation does not connect to a regular nonlinear solution, unless the angle distribution consists of only two beams.

We extend Maxwell electrodynamics with a Chern–Simons coupling to a dark sector top form sourced by domain walls. Cosmic birefringence can arise from a distinct mechanism in which photon polarization is rotated when crossing vacuum interfaces, rather than through adiabatic propagation in a background field. Ultrathin walls induce a finite, frequency-independent polarization rotation generated by an electromagnetic Chern-Simons term localized at the interface. The effect persists even in the absence of ultralight axions or other propagating scalar degrees of freedom. For phenomenologically viable parameters, such walls can generate cosmic microwave background polarization rotation at the level $\Delta\vartheta \sim 10^{-3}$ rad, providing a signature of the topological structure of the dark-sector vacuum.

One of the fundamental challenges of quantum gravity is to understand how the microscopic degrees of freedom of the cosmological horizon shape the evolution of the Universe. One possible approach to this problem is based on the Barrow–Tsallis entropy. This entropy accounts for both quantum gravitational effects and the nonextensive effects inherent in any long-range interaction. By employing an inverse cosmographic reconstruction of the model parameters, we derive a relation between the Barrow parameter, which encodes the microscopic deformation of the horizon geometry, and the Tsallis parameter, which characterizes macroscopic nonextensivity. Within the IR–UV correspondence, this relation determines the scaling of the microscopic length uncertainty in terms of the current cosmographic parameters and demonstrates how long-range nonextensive effects alter the standard Karolyhazy-type scaling. We also used our method for finding the parameters of cosmological models to evaluate the feasibility of using fractional derivatives to describe the late evolution of the Universe. The resulting relationships are exact. Therefore, the uncertainty in the relationship between the model parameters depends only on the current uncertainty in the values of the cosmographic parameters.

The spontaneous breaking of $SO(10)$ via flipped $SU(5)$ to the Standard Model yields a novel scenario in which the superheavy topologically stable GUT monopole carrying a single unit ($2\pi/e$) of Dirac magnetic charge emerges from the merger of a confined but topologically distinct monopole-antimonopole pair that are pulled together by a string. The $SO(10)$ breaking via the subgroup $SU(4)_c\times SU(2)_L\times SU(2)_R$, following a similar reasoning, produces a topologically stable monopole that carries two units ($4\pi/e$) of Dirac charge. We explore the cosmological consequences of this scenario by assuming that the monopoles and strings experience a limited number of inflationary $e$-foldings, before re-entering the horizon and ultimately forming a network of quasi-stable strings bounded by monopole-antimonopole pairs. We identify regions of the parameter space that yield an observable number density of the GUT monopole from the collapse of the appropriate string segments. The gravitational waves emitted by these quasi-stable cosmic strings lie in the Hz to kHz range, which can be tested in a number of proposed and ongoing experiments.