Main

Current solar flare predictions often lack precise quantification of their reliability, resulting in frequent false alarms, particularly when dealing with datasets skewed towards extreme events. To improve the trustworthiness of space weather forecasting, it is crucial to establish confidence intervals for model predictions. Conformal prediction, a machine learning framework, presents a promising avenue for this purpose by constructing prediction intervals that ensure valid coverage in finite samples without making assumptions about the underlying data distribution. In this study, we explore the application of conformal prediction to regression tasks in space weather forecasting. Specifically, we implement full-disk solar flare prediction using images created from magnetic field maps and adapt four pre-trained deep learning models to incorporate three distinct methods for constructing confidence intervals: conformal prediction, quantile regression, and conformalized quantile regression. Our experiments demonstrate that conformalized quantile regression achieves higher coverage rates and more favorable average interval lengths compared to alternative methods, underscoring its effectiveness in enhancing the reliability of solar weather forecasting models.

With rapid improvements in the assembly of large samples of galaxy clusters, we are approaching the ability to study clusters at $z\gtrsim2$. Evolutionary studies comparing these distant clusters to the clusters in our local universe depend heavily on the reliability of low-redshift cluster samples, most of which are subject to X-ray selection effects, biasing them to relaxed, cool core clusters. Here, we introduce the Cluster Evolutionary Reference Ensemble At Low-$z$ (CEREAL) sample, composed of Chandra X-ray observations of 169 galaxy clusters that have been selected from the Planck Sunyaev-Zel'dovich catalog. CEREAL has a simple and well-understood selection function, spans an order of magnitude in mass at $z\sim0.15$, and has uniform, high-resolution X-ray follow-up. We present the full sample and provide results based on X-ray surface brightness properties, finding significantly more non-cool core systems than in X-ray-selected samples. We use surface brightness concentration (c$_\mathrm{SB}$) as a proxy for cool core strength and centroid shift ($w$) to measure dynamical state. Over the full sample, we find a cool core (c$_\mathrm{SB} > 0.075$) fraction of $0.39_{-0.04}^{+0.04}$, a strong cool core (c$_\mathrm{SB} > 0.155$) fraction of $0.13_{-0.03}^{+0.03}$, and a dynamically relaxed ($w<0.01$) fraction of $0.42_{-0.04}^{+0.04}$. We find no mass dependence in the fraction of clusters that appear relaxed or have cool cores. We quantify the rarity of X-ray-bright central point sources (L$_\mathrm{nuc,~2-10~keV} > 10^{43}$ erg s$^{-1}$), finding them to be intrinsically rare ($0.7_{-0.5}^{+1.2}$\% of massive, low-z clusters) with a notable increase in occurrence rate at the centers of cool cores.

We study the neutrino distributions that arise in a simulation of a neutron star merger that uses a Monte Carlo (MC) neutrino transport scheme. In a snapshot taken 1 ms after merger, we calculate relevant observables to test when neutrinos behave like a thermalized gas, and when a free-streaming picture is more appropriate. We find that in hot, dense regions where neutrino-matter interactions are frequent, MC neutrino and antineutrino distributions are consistent with thermalized neutrinos. In moderately warm regions, where neither approximation is expected to hold, we find significant departures from the predictions of the thermalized-neutrino approximation, particularly for the (anti)neutrino average opacity and net rate of absorption per baryon, even when average energies appear approximately thermal. At lower temperatures, MC results approach the free-streaming limit. Our results demonstrate that energy-averaged agreement with thermalized-neutrino assumptions does not guarantee accurate weak interaction rates. Non-equilibrium aspects of the neutrino distribution are therefore crucial for neutrino-mediated microphysics such as composition evolution in the early post-merger phase.

Although restoration of chiral symmetry is predicted by quantum chromodynamics to take place at high baryon density, most modeling of neutron star interiors disregards a chiral phase transition. We model neutron star cores with a parity doublet model, which allows for dynamical chiral symmetry restoration and predicts the appearance of the parity partners of nucleons and hyperons at large densities, as well as deconfined quark matter. We study the thermal evolution of neutron stars, focusing for the first time on the impact of Urca processes involving the parity partners in neutron star cooling simulations. We find that Urca processes for the parity partners of the nucleons significantly affect the thermal evolution of massive stars and allow for improved agreement with observed surface temperature and ages.

We study structure formation in alternative cosmological models constrained by background observations, including $\Lambda$CDM, wCDM, the Chevallier-Polarski-Linder parametrisation and a flexible Chebyshev expansion of the dark energy equation of state. The models are constrained using baryon acoustic oscillations, cosmic microwave background, cosmic chronometers and strong lensing measurements. Using the best-fitting parameters, we generate cosmology-dependent initial conditions and perform N-body simulations to analyse the matter power spectrum, halo mass function and halo density profiles. Although all models remain broadly consistent with $\Lambda$CDM at the background level, differences in the physical matter density $\Omega_{0m}h^2$ and in the expansion history $H(z)$ lead to distinct growth histories that are amplified by non-linear evolution. We find a clear hierarchy in the power spectrum amplitude and in $\sigma_8$, with the Chebyshev and CPL models exhibiting enhanced small-scale power, earlier halo formation at $z\gtrsim2$ and a migration of excess toward higher masses at late times. The wCDM model displays milder and partially compensating effects driven by its different expansion history. When expressed in terms of the scaled radius $r/R_{200c}$, halo density profiles show a high degree of universality across cosmologies, indicating that internal halo structure is largely governed by the same gravitational dynamics. These results demonstrate that even modest background-level variations in $w(z)$ can translate into coherent non-linear signatures, highlighting the constraining power of large-scale structure observables in extended dark energy models.

We present the first constraints on multipole-dependent cosmic birefringence using CMB polarization data from the BK18 dataset, which combines observations from BICEP2, Keck Array, and BICEP3 at frequencies of 95, 150, and 220 GHz. Photon coupling to an axion-like field leads to the rotation of CMB polarization, inducing non-zero EB cross-correlations. We show that a multipole-dependent rotation beta(l) imprints a distinct signature in the polarization spectra that can be constrained. Specifically, we consider an Early Dark Energy (EDE) scenario in which a pseudoscalar field couples to photons through a Chern-Simons interaction, generating a polarization rotation with multipole dependence. We introduce a phenomenological beta(l) as a step function, obtaining constraints on the step function size consistent with zero, with uncertainties less than 0.15 degrees (68% CL). In addition, using multi-frequency EE, BB, and EB cross-spectra, along with robust BICEP/Keck foreground treatment and likelihood framework, we derive constraints on the axion-photon coupling amplitude g for several choices of EDE parameters. For the baseline best-fit value f_EDE = 0.087 from the Planck 2018 analysis, we obtain g = 0.11 +/- 0.37 (68% CL), consistent with previous limits.

Aims: We aim to characterise the multiscale statistical properties of the reconstructed velocity and density fields of the nearby universe, identify possible scaling regimes, quantify intermittency, and assess indications for the transition toward large-scale homogeneity within the range probed by current data. Methods: We analyse the CosmicFlows4 three-dimensional velocity and density-contrast cubes using absolute structure functions of arbitrary order, $q$. The analysis is performed within a volume extending to $z \lesssim 0.08$ ($\simeq 350~\mathrm{Mpc}$ $h^{-1}$). Structure function scaling exponents $\zeta(q)$ are estimated from configuration-space statistics. Intermittency is characterised using the Universal Multifractal formalism, and probability density functions of increments are examined. Results: Two regimes are detected. Small separations are dominated by reconstruction smoothing and show nearly linear $\zeta(q)$ behaviour. At larger separations, a scaling regime appears with $\zeta_\rho(1)\simeq0.3$ ($D_\rho\approx3.7$) and $\zeta_v(1)\simeq0.4$. The correlation function follows $\xi(r)\sim r^{-1.4}$ over $[45,250]~\mathrm{Mpc}\,h^{-1}$, implying $D_2\simeq1.6$. Non-linear $\zeta(q)$ and Lévy-stable increment PDFs indicate intermittency and strong non-Gaussianity. Velocity increments show a systematic negative skewness suggestive of a cascade-like asymmetry associated to amplification of negative, compressive gradients.

We present a spectroscopic analysis of XRISM and NuSTAR observations of the black hole X-ray binary GS 1354–64 during its 2026 outburst. A total number of 3.5 million photons are collected by the microcalorimeter Resolve on board XRISM, providing an unprecedented high-resolution view of the iron line profile. A clean broad iron line is found in the data, without significant narrow features. Modeling the broad iron line with relativistic reflection from the inner accretion disk suggests a rapidly spinning black hole (a>0.98) in the system. Measurements of the disk inclination angle from the reflection method are model-dependent. This work demonstrates the power of X-ray microcalorimeters in studying the inner accretion flow and constraining black hole parameters.

QSOs are essential for investigating the structure and evolution of the Universe. Historically, their identification has been concentrated in the northern hemisphere, primarily due to the sky coverage of major astronomical surveys. The QUBRICS survey, started in 2019 to address this asymmetry, has identified more than 1300 new bright (i<19.5) high-redshift (2.5<z<6) QSOs in the southern sky. We aim to quantify, using an independent QSO sample, the completeness and recall of the QUBRICS QSO selection methods, based on XGB (eXtreme Gradient Boosting) and PRF (Probabilistic Random Forest), since completeness is a fundamental metric for ensuring the statistical robustness of QSO-based cosmological investigations. A subset of Gaia DR3 sources with low-resolution spectra was analyzed, obtaining a sample of 3501 QSOs. To determine how many QSOs were correctly identified as candidates, we crossmatched this independent sample with the datasets used for selection: 894 QSOs with z>2.5 fell within the XGB dataset footprint, of which 152 were unclassified and thus eligible for completeness testing. Similarly, 675 QSOs with z>2.5 were within the PRF dataset footprint, including 69 unclassified objects. The XGB correctly identified as candidates 136 (89%) of the 152 QSOs with z>2.5 present in its dataset as unclassified objects. The PRF correctly identified as candidates 46 (66%) of the 69 QSOs with z>2.5 present in its dataset as unclassified objects. These findings confirm the high efficiency of the QUBRICS selection methods (recall=89%) and provide the completeness estimate for spectroscopically confirmed QSOs (82%), necessary for cosmological studies using QUBRICS data. This work also provides reliable redshifts for 1223 new QSOs (median redshift z=2.1 and magnitude G=17.8), that will help improve the performance of future selections.

3C 286, a well-known radio calibrator, exhibits the stability in both of total flux density (FD) and polarization parameters. However, its stable and luminous interstellar radio signal may encounter interplanetary scintillation (IPS) due to density irregularities in the solar wind within the heliosphere. In this work, we analyze high-time-resolution observations of 3C 286 obtained with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) from 2019 to 2023. Our analysis reveals that IPS affects the polarized flux densities of the Stokes I, Q, and U parameters, whereas Stokes V shows no detectable IPS-induced variations. The IPS variations detected in Stokes I are synchronous with those in Stokes U, while those in Stokes Q exhibit greater randomness. The crosscorrelation function (CCF) results indicate no time delay between Stokes I and U, but a delay of approximately 2.8 seconds between Stokes I and Q. This suggests that the different polarized radio emissions of 3C 286 originate from distinct emission regions, specifically the core and the southwestern jet. Furthermore, the projections of the radio core and jet component onto the scintillation screen at 1 AU yield a solar wind plasma speed of $\sim 637$ km/s.

Azimuthal arcs in millimeter continuum emission from protoplanetary disks are often attributed to dust-trapping vortices, but definitive observational confirmation of vortices remains lacking. We present sub-0.1" resolution ALMA continuum observations of the HD 34282 disk at 0.9, 1.3, 2.1, and 3.1 mm. These observations resolve a bright azimuthal arc superposed on a compact double-gap, triple-ring morphology, most clearly at shorter wavelengths, and enable us to probe the physical origin of the arc. It exhibits a lower spectral index than the surrounding rings, consistent with enhanced grain growth and/or higher dust surface density of a dust-trapping vortex. Its azimuthal width decreases with increasing wavelength, consistent with tighter confinement of larger grains, or lower optical depths at longer wavelengths. These observations probe dust with Stokes numbers St < 0.03. Vortex models predict negligible peak shifts in this regime, consistent with the 1.3 to 3.1 mm data. At 0.9 mm, however, the arc peak is offset by 15 +/- 4 degree in the direction of disk rotation relative to longer wavelengths, and the near-side ring emission is locally dimmer compared to the far-side, likely reflecting optical-depth or temperature effects. These observations are consistent with azimuthal dust trapping, potentially associated with a vortex-induced pressure maximum.

Fast radio bursts (FRBs) provide a sensitive probe of ionized baryons through their dispersion measure (DM). In addition to slowly evolving cosmological terms, at least two repeaters now show clear secular DM-decrease episodes: FRB 20190520B and FRB 20121102 , supporting a dense, dynamically evolving local environment. We adopt a \emph{forward-modeling} approach and use time-dependent 1D SNR simulations for a young magnetar embedded in SN ejecta, combining single-star and binary-stripped progenitors with HD+NEI calculations to follow shock structure, ionization, and electron density. The shocked region contributes only limited DM ($\lesssim10\,{\rm pc\,cm^{-3}}$), while the dominant time-varying component is the unshocked ejecta, whose early behavior follows ${\rm DM}\propto t^{-\alpha}$ with $\alpha\simeq1.8$–$1.9$. Although shocked-region DM is small, shock-amplified magnetic fields can still generate substantial RM; in our shock-only RM framework, only the $11\,M_\odot$ SS model reproduces the FRB 20121102 RM evolution. Binary-stripped progenitors generally yield smaller DM than single-star models at fixed $M_{\rm ZAMS}$, with composition-dependent mean molecular weights introducing non-monotonic mass trends. Matching the observed ${\rm dDM}/{\rm d}t$ of FRB 20190520B (and the late-stage slope of FRB 20121102), we infer local SNR DM contributions of tens to hundreds ${\rm pc\,cm^{-3}}$. We also find GHz escape is allowed in most models, with $\tau_{\rm ff}=1$ typically reached by $t_{\rm esc}\lesssim70$ yr; for weakly ionized ejecta, the source can be nearly transparent from very early times. These results support a young CCSN/SNR origin for a substantial fraction of ${\rm DM}_{\rm source}$ and highlight that physically consistent local-environment modeling is essential for robust FRB cosmological DM inferences.

We report the first NuSTAR hard X-ray observations of the low-luminosity active galactic nucleus NGC 4278. The source is clearly detected beyond 10 keV with a hard X-ray spectrum consistent with a power law of photon index between 2.2 and 2.5 without evidence for a high-energy cutoff. The X-ray flux is low compared to the active state in 2021, but exhibits variability by a factor of  2 on a timescale of a month. We discuss the origin of the hard X-ray emission and explore its connection to gamma rays and high-energy neutrinos. We explain the X-ray data, including both quiescent and active states, using a radiatively inefficient accretion flow (RIAF) model with a variable accretion rate. We also show that TeV gamma rays cannot escape from the RIAF disk, and very high-energy gamma rays observed in LHAASO are likely to originate from outer regions such as jets and winds, which is consistent with our results favoring a magnetically arrested disk. We also discuss hidden neutrino emission from RIAFs together with possible connections to coronae of active galactic nuclei with standard, radiatively efficient disks.

Pulsar halos are produced by electrons and positrons diffusing in the interstellar medium around their parent pulsar wind nebulae. Recent observations by HAWC and LHAASO have revealed asymmetric morphologies in the halos surrounding Geminga and Monogem. The anisotropic diffusion model provides a natural explanation for such asymmetries, where the morphology is determined by the viewing angle of the mean magnetic field, the Alfvénic Mach number ($M_{\rm A}$), and the pulsar distance. In this work, we model the measured morphologies based on this framework and constrain the properties of interstellar magnetic turbulence. We find that the mean magnetic field orientations within the two halos are different, implying that they reside in different magnetic coherence regions, whereas the Alfvénic Mach numbers are relatively close ($M_{\rm A}\sim 0.2$). The results suggest a local magnetic field coherence length of approximately 100pc. Our study demonstrates that the morphology of pulsar halos serves as a powerful diagnostic tool for the properties of interstellar magnetic fields, highlighting the need for more accurate morphological measurements and sophisticated diffusion modeling in future studies.

We present that supernovae interacting with a dense Si and S-rich circumstellar matter like SN 2021yfj can originate from mergers of two white dwarfs. A C+O white dwarf accreting He from its non-degenerate He companion star can initiate a C burning frame at its surface propagating inward under certain conditions. Such a burning frame synthesizes intermediate mass elements such as Si and S, forming a hybrid WD with an outer Si+S-rich layer. After the He star companion becomes a white dwarf, the two white dwarfs can eventually merge. During the merger, the outer layers of the hybrid white dwarf can be tidally stripped, forming a dense Si and S-rich circumstellar matter. If a thermonuclear explosion is triggered in the merging white dwarfs, an explosion within a dense Si and S-rich circumstellar matter can be realized, resulting in SN 2021yfj-like events. We argue that the properties of SN 2021yfj can be reproduced by a dense Si and S-rich circumstellar matter having   0.3 Msun within which an explosion having kinetic energy of   4e50 erg and ejecta mass of   0.3 Msun occurred. These properties are consistent with the double white dwarf merger scenario. This scenario can naturally explain the existence of He observed in SN 2021yfj. Because white dwarf mergers can also lead to the formation of He and C+O dense circumstellar matter, some Type Ibn and Icn supernovae may also originate from a similar evolutionary path.

Fast radio bursts (FRBs) are millisecond-duration extragalactic radio transients, some of which are associated with compact persistent radio sources (PRSs), hinting at a physical connection. While several models have been proposed to explain PRSs and their connection to FRBs, direct observational tests remain limited. Here, we report for the first time a correlated trend between the long-term variation of the PRS flux density and the burst energetics of FRB 20190520B and FRB 20240114A, suggesting that both the PRS and FRB activity may be powered by a shared energy reservoir. We further examine additional repeaters with compact PRSs and find no clear correlation between PRS luminosity and burst activity, likely due to the limited observations. These results are consistent with scenarios in which both the PRS and FRB activity may be powered by a common energy reservoir, such as the magnetic or rotational energy of a magnetar.

The existence of dark matter is supported by multiple astrophysical observations, yet its particle nature remains unknown. The development of gravitational wave astronomy, especially with future space-based detectors such as LISA, provides new opportunities to study the interactions between dark matter and compact-object systems. This review summarizes the main dark matter candidates and their macroscopic distributions, and highlights three mechanisms through which dark matter can affect gravitational wave observations: (1) modifications to compact-object orbits and the dynamics of systems such as extreme mass-ratio inspirals, including dark matter spikes, dynamical friction, and potential perturbations; (2) gravitational lensing effects induced by the spatial distribution of dark matter, altering waveform amplitudes and phases; and (3) direct couplings between ultralight dark matter fields and detectors. As low-frequency gravitational wave detection techniques are proposed and continue to develop, these effects may offer a novel avenue for probing the properties of dark matter, and combining precise waveform modeling with multi-messenger observations could reveal insights into its microscopic structure.

A discrete tomography method has been developed that is able to reconstruct three-dimensional coronal mass ejection (CME) density structure. We test the method by producing synthetic coronagraph imagery for three CME events using the CORona–HELiosphere (CORHEL) model. We combine images from different numbers of observing spacecraft, ranging from three to seven, and we perform the method separately using polarimetric and non-polarimetric reconstructions, as a means to test their relative effectiveness. In all cases, we show that increasing the number of observing spacecraft reduces the mean relative absolute error (MRAE) between the simulated and reconstructed density. Furthermore, the MRAE is generally lower when using polarimetric reconstructions compared to non-polarimetric reconstructions. Methods applied to localise the CME front work well for all spacecraft configurations, and are improved when using polarimetric, over non-polarimetric, reconstructions. The presence of a CME front in the simulations can be identified with an accuracy of $(72\pm9)\%$, $(70\pm8)\%$ and $(52\pm12)\%$ for CME1, CME2 and CME3 via polarimetric reconstructions using only three spacecraft at L1, L4 and L5. The radial position of the CME front can be constrained to a high level of precision when using polarimetric reconstructions using the same three spacecraft; $0.003\pm0.004$ au, $0.004\pm0.005$ au and $0.005\pm0.004$ au for CME1, CME2 and CME3, respectively. We expect that at least four spacecraft are required in order to derive accurate information about 3D CME structure. We find no strong evidence of improvement when including out-of-ecliptic observers, but that their inclusion increases the volume of space within which the inversion can be performed.

We investigate the magnetic origin of the coronal mass ejection (CME) that occurred on November 28, 2023, at 19:50UT from active region (AR) 13500 located near the solar disk-center. The eruption was associated with an S-shaped sigmoidal structure formed by the inner AR polarities along a sheared polarity inversion line, while the outer polarities evolved through proper motions. During November 26-28, the AR exhibited a decrease in net magnetic flux while progressively injecting magnetic helicity and energy into the corona toward the eruption onset, highlighting the key role of helicity-injection in triggering eruptions. To simulate this magnetic evolution, we employed a data-driven magnetofrictional (MF) simulation starting 2.8 days prior to the eruption. The energy input for the model was constrained using the observed energy injection through an ad-hoc parameter. The initial potential-field configuration gradually evolved into a sheared-arcade and eventually developed into a twisted flux rope (FR) over the observed time-scale. Proxy emission maps based on electric currents show remarkable morphological agreement between the simulated and observed sigmoidal structure. The average FR-core twist increasingly builds-up leading the FR to initiate slow-rise motion of FR top from 50Mm until its eruption onset at 80Mm. Importantly, the ratio of current-carrying to total relative-helicity increased from 0.13 at FR formation to 0.30 at eruption, when the FR core entered the torus-unstable regime, suggesting an association between torus-instability and a threshold helicity ratio. These results demonstrate that data-driven MF simulations can successfully reproduce the evolving coronal magnetic configuration and may provide a robust tool for assessing the eruptive potential of ARs, particularly the helicity ratio.

We show that the flip-flop transitions in X-ray binaries (rapid cycling between different spectral states which are sometimes seen near the global state transition) show a series of analogies to the changing state phenomena (rapid changes in the emission line properties that seem to be driven by changes in the central engine) in active galactic nuclei (AGN). Specifically, (1) the timescales for the transitions scale approximately linearly with mass and (2) both phenomena occur at a few percent of the Eddington luminosity. Because most accretion physics is expected to be scale-free, it is likely that these represent two manifestations of the same phenomena. Demonstrating this would allow the use of a much wider range of observational techniques, on a much wider range of characteristic timescales, and provide a clearer pathway toward understanding these rapid transitions than is currently available. We discuss potential means to establish the connection more firmly, and to use the combination of the observational advantages of both classes of systems to develop a better understanding of the phenomenon.

Black hole-neutron star mergers, together with binary neutron star mergers, are key laboratories for neutron star physics. They enable us to probe merger dynamics imprinted in gravitational waves and potential electromagnetic counterparts. These systems link microphysics and macrophysics by placing constraints on the dense matter equations of state, potentially revealing the imprint of hadron-quark phase transitions, clarifying the role of neutrino irradiation in shaping the ejecta, its r-process nucleosynthesis, and kilonova emission, as well as assessing how magnetically driven instabilities affect mass ejection and possible electromagnetic signatures. Despite their importance, black hole-neutron star mergers remain relatively less studied and therefore not yet well understood, largely due to the lack of publicly available numerical relativity setups suitable for such investigations. In this work, we present a fully reproducible black hole-neutron star merger simulation performed exclusively using Einstein Toolkit thorns, targeting the detected event \texttt{GW230529}. The simulations are carried out at three resolutions with finest grid spacings of $162$, $222$ and $310$ meters to assess numerical robustness. The entire setup, from initial data to a parameter file with some of the analysis scripts, is publicly released as a new Einstein Toolkit gallery example and will be distributed as part of the Hypatia release, establishing a reference black hole-neutron star merger configuration within the Einstein Toolkit.

Catastrophic planetesimal disruptions offer a unique opportunity to study and characterize large planetesimal populations in exoplanetary systems that are not currently detectable by modern observatories. The unexpected discovery of a second collision event in the Fomalhaut system raises important questions about the planetesimal population and dynamical state inside the Fomalhaut main belt that led to two collisions in 20 years. We present a statistical model developed and applied to the archetypal Fomalhaut system to provide new constraints on the bulk properties of the planetesimals in Fomalhaut's main belt. Utilizing the constraints provided by the spatially resolved Fomalhaut cs1 and cs2 collision events, we retrieve the belt parameters that best reproduce the observed collision rate while remaining consistent with the system's age and dust mass. Our best-fit model suggests a total main belt mass of 200-360 $M_{\oplus}$, with the transition from a collisionally evolved to a primordial planetesimal population occurring at a radius of $115_{-10}^{+30}$ km and a maximum planetesimal radius of $380_{-202}^{+643}$ km. We estimate a catastrophic collision rate of $0.086_{-0.048}^{+0.067}$ collision events per year for planetesimals with radii $\ge$ 100 km in the region interior to the main belt. Our findings show that further observable collisions are likely, motivating continued monitoring of Fomalhaut and other nearby debris disks.

We present an analysis that investigates the $H_0$ and $S_8$ tensions by considering a dark energy model. The latter is a late-time model characterized by a future abrupt event known as the Little Rip (LR) model and characterised by one extra parameter, $\beta$, compared to the standard model, $\Lambda$CDM. To test this approach, we perform a statistical analysis by the MCMC method using the most recent observational data. We obtain a positive correlation in ($H_0$, $\beta$) plane. We also note that the Hubble tension is less than $3\sigma$ when using early measurements, i.e., Cosmic Microwave Background (CMB) data, and when combining it with Baryon Acoustic Oscillation (BAO) data, but it is no longer so when we combine early and late measurements (i.e. PantheonPlus (PP)). In addition, we test the model with DESI-DR2 combined with CMB and recent SNIa measurements. We notice that our model shifts toward the quintessence field. For a complete statistical analysis, we use the Akaike Information Criteria and Bayesian analysis of the evidence. According to Bayes factors, we find that the LR model provides an improved fit only to CMB data.

As the oldest building blocks of our Galaxy, globular clusters retain the archaeological footprint of the early stellar environments. Accurate absolute ages of globular clusters are required to interpret this ancient record. Existing dating techniques often produce precise but discordant ages, suggestive of systematic errors in excess of 1 Gyr. The James Webb Space Telescope (JWST) has unlocked a new dating method that leverages the cooling behavior of previously unobservable brown dwarf members. With a largely independent set of systematic errors, this new method provides a new consistency test for more established methodologies. I present a likelihood-based histogram-free method to derive globular cluster ages from multi-band JWST photometry of cluster members near and below the hydrogen-burning limit. By applying the method to a large set of simulated observations, I establish that formal age errors (i.e. errors based on measurement uncertainties alone) under 0.2 Gyr are attainable for nearby globular clusters. I also evaluate the significance of associated systematic effects, including the chemical heterogeneity of globular clusters (multiple populations), unresolved binary systems and uncertainties in brown dwarf cooling rates. As with other methods of age determination, systematic effects dominate the error budget (in selected cases, by over an order of magnitude), but may be reduced with more sophisticated analysis. Finally, I provide a lookup table for determining the number of observations, exposure times and temporal baselines required to estimate the age of a given cluster to a prescribed precision.

We propose a unified model for jet formation applicable to active galactic nuclei, young stellar objects, and X-ray binaries. In this model, the binding energy released from the accretion disk is primarily stored as turbulence rather than being radiated away, leading to the formation of advection-dominated accretion flows. Near the central object, a thick accretion disk with funnel-like structures develops. Within the turbulent flows, the smallest stable blobs can be accelerated beyond the escape velocity through a mechanism involving the combined effects of inward pressure forces and angular momentum conservation. These rapidly moving blobs may exit through the funnels, collectively forming two opposing jets. This model predicts that jets originate from the innermost region of the thick disk surrounding the central object. It can be extended to account for jet formation in active galactic nuclei, young stellar objects, X-ray binaries, and other analogous astronomical systems.

A growing number of white dwarfs (WDs) exhibit one or more signs of remnant planetary systems, including transits, infrared excesses, and atmospheric metal pollution. WD 1054-226 stands out for its unique, highly structured, and persistent photometric variability. We aim to investigate the long-term stability and nature of the periodic signals observed in WD 1054-226 to better understand the origin and evolution of its transiting material. We analyse all available TESS light curves from Sectors 9, 36, 63, and 90 using Lomb-Scargle (LS), Box-Least-Squares (BLS), and Gaussian process (GP) periodogram analyses. We complement these with multiband, high-cadence ground-based photometry from LCOGT, MuSCAT2, ALFOSC, and ProEM to test for colour dependence and confirm the periodicities. We confirm the persistence of the previously-reported 25.01 h and 23.1 min periodicities over a six-year baseline. The 25.01 h signal shows some temporal evolution, while the 23.1 min dips are highly coherent on long timescales. A transient 11.4 h feature, previously reported, is detected only in early TESS sectors and is absent in recent data. No significant colour dependence is found in the ground-based observations. The stability of both the 25.01 h and 23.1 min signals indicates a long-lived, dynamically sculpted debris structure around WD 1054-226. The lack of colour dependence implies high optical depth, consistent with an opaque, edge-on debris ring rather than an optically thin dust population. This makes WD 1054-226 a key laboratory for testing models of remnant planetary systems around white dwarfs.

Realistic modeling of stellar spectra requires accurate radiative opacity coefficients. Owing to the fragmentary nature of existing data from rigorous quantum-mechanical calculations, photoionization coefficients based on the rigid-wavefunction approximation remain the only practical option for studies of magnetic white dwarfs. Although variants of this approach have been widely used in spectral analyses for decades, a complete and explicit treatment of degeneracy-level breaking has not previously been presented. In this work, we provide a comprehensive description of this procedure, including explicit expressions for the photoionization probability of individual bound-free transitions as functions of magnetic field strength and radiation polarization. We also evaluate the occupation numbers of bound states in a magnetized gas under ionization equilibrium, enabling the calculation of absolute photoionization opacities. Because high-lying atomic states are strongly perturbed by the magnetic field and ultimately dissolved, substantial modifications of the monochromatic absorption are found even for fields below 10 MG–a regime where fully rigorous quantum calculations are numerically demanding and have not yet been applied. Over a wide range of magnetic field strengths, pronounced dichroic features appear in the hydrogen continuum absorption.

Theoretical formation models and exoplanet detection surveys indicate that systems with multiple giant planets are common. We investigate how multiple super-thermal mass planets embedded in a circumstellar disk shape the dust distribution and examine the consequences for interpreting disk substructures and inferring planetary properties. We perform two-dimensional hydrodynamical simulations with a modified PLUTO code, treating dust as Lagrangian particles in a wide range of sizes. We analyze systems with two planets of different masses and orbital separations, comparing them to the single-planet scenario. We generate synthetic ALMA continuum maps using RADMC-3D and compute the relative impact velocities of dust particles to assess potential limitations to grain growth. Dust morphologies in multi-planet systems cannot be described as a simple superposition of single-planet gaps. Secular planetary perturbations can generate multiple dust traps and asymmetric structures, while also exciting significant eccentricities in dust particle orbits. As a consequence, the locations and widths of dust rings and gaps depend on the size of the particles, the masses of the planet, and the orbital configurations. Synthetic continuum images may hide gaps carved by multiple planets, thereby complicating the interpretation of observed substructures. In addition, eccentricities induced in dust orbits lead to stronger gas drag, reducing the Stokes number for a given particle size, and the enhanced relative velocities associated with eccentric orbits can further suppress grain growth, promoting fragmentation and replenishment of small dust grains.

The chiral plasma instability (CPI) has been invoked as a possible mechanism for generating primordial magnetic fields in the universe and ultrastrong fields in neutron stars. We investigate chiral dynamos where the chirality imbalance is pumped by a source on a timescale $t_0$ and show that the CPI rate $\gamma$ is limited to $\gamma_0/(1+{\cal Q}^2)$, where ${\cal Q}= (\gamma_0 t_0)^{1/3}$ and $\gamma_0$ corresponds to models with instantaneously created chirality imbalance $(t_0=0)$. We then find that chiral flipping with rate $\Gamma_{\mathrm f}$ hinders the chiral dynamo if $\Gamma_{\mathrm f} >\gamma_0/(1+{\cal Q}^2)$ and completely suppresses it if $\Gamma_{\mathrm f} >\gamma_0/(1+{\cal Q}^{3/2})$. Realistic $t_0$ typically give ${\cal Q}\gg 1$, which makes the dynamo greatly vulnerable to the suppression by chiral flipping. The suppression is strong in protoneutron stars and may be (barely) avoided near the electroweak transition in the early universe.

While space-based telescopes offer unparalleled precision for asteroseismology, ground-based observations remain crucial for identifying compact pulsator candidates and enabling their pulsational study through multi-colour photometry. The BlackGEM telescope array significantly enhances the detection and characterisation possibilities for compact pulsators that tend to be much fainter than dwarf or giant pulsators. Using BlackGEM multi-colour photometry of the hot pre-white dwarf PG 1159-035, we demonstrate its capability to detect short, multi-periodic pulsations with amplitudes down to a few milli-magnitudes. The primary aim of this study is to establish the feasibility of pulsation mode identification in hot subdwarfs and white dwarfs via mode amplitude-ratio analysis derived from BlackGEM multi-colour observations. Pulsation frequencies were extracted from our target using iterative pre-whitening analysis. To validate our data-driven mode identification concept using multi-colour photometry, we used the well-studied hot pre-white dwarf PG 1159-035, with previously identified pulsation modes, as a prototypical object that served for validation. The pre-whitening analysis using BlackGEM's standard q-, u-, and i-band light curves of PG 1159-035 revealed pulsation frequencies of l=1 and l=2 modes, consistent with values obtained from the literature. Subsequent amplitude ratio calculations resulted in discernible distributions for the l=1 and l=2 modes. The future assembly of more BlackGEM amplitude ratios for well-known white dwarfs with already identified modes will lead to density estimators suitable for identifying newly detected modes in known or as-yet-undiscovered pulsators. Our proof-of-concept study paves the way for large-scale asteroseismic analyses of optically faint compact pulsating stars using ground-based facilities, such as BlackGEM.

We report the discovery of a new long-period radio transient, ASKAP J142431.2-612611, with a 36 minute period, identified in the Australian SKA Pathfinder Evolutionary Map of the Universe survey. We detected pulsed emission from ASKAP J142431.2-612611 over a period of eight days during follow-up observations with the Australia Telescope Compact Array, after which the source appears to have switched off. No optical or near-infrared counterpart is detected in archival surveys or in targeted Gemini South FLAMINGOS-2 observations. During its active state, the source exhibits a stable pulse profile with fractional polarisation consistent with 100%, evolving from elliptically to linearly polarised and tracing a well-defined great-circle trajectory on the Poincaré sphere. We show that this behaviour is consistent with fully linearly polarised intrinsic emission modified by propagation through a linearly polarised birefringent medium. This discovery expands the known population of long-period transients and highlights the intermittent nature of their activity. We discuss the implications for proposed models of long-period transients and outline future observations needed to constrain the origin of their intermittency and polarisation properties.

We present a particle physics motivated realization of interacting dark energy in which a radiatively stable dark energy sector couples to weakly interacting massive particle dark matter through pure momentum exchange. The dark energy field arises as a pseudo-Nambu-Goldstone Boson from a complex scalar singlet charged under a softly broken global $U(1)_S$, while dark matter is identified with an inert scalar doublet stabilized by a discrete $Z_4$ symmetry. This symmetry structure allows renormalizable dark matter-dark energy portal operators; however, requiring the dark energy field to emerge as a radiatively stable pseudo-Nambu-Goldstone Boson necessitates their absence, leaving derivative interactions as the leading coupling. As a result, energy transfer between the dark sectors is absent at the background level, while momentum exchange modifies the evolution of cosmological perturbations. We implement the resulting interacting dark energy model self-consistently in the Boltzmann code CLASS and study its impact on the growth of structure. We find that, despite sizeable momentum exchange, the suppression of the clustering amplitude $\sigma_8$ saturates above the level required to fully resolve current low-redshift tensions. Our results demonstrate that symmetry-protected, momentum-exchange-only dark sector interactions possess an intrinsic limit on structure suppression, providing a theoretically controlled benchmark for interacting dark energy scenarios.

Astrophysical neutrinos from Active Galactic Nuclei (AGN) offer a unique window into high-energy particle acceleration in obscured environments. The nearby Type II Seyfert galaxy NGC 1068 is a compelling example, exhibiting evidence of a high-energy neutrino excess without an associated TeV $\gamma$-ray counterpart. This suggests that hadronic processes may occur within an inner, magnetically dominated region, where the TeV emission is suppressed by $\gamma \gamma$ absorption and reprocessed via electromagnetic cascades in the dense, obscured environment. Building on our framework, which establishes turbulence-driven magnetic reconnection as the main driver for particle acceleration in this source, we present a refined lepto-hadronic model based on de Gouveia Dal Pino & Lazarian (2005) and Kadowaki et al. (2015). In these proceedings, we adopt a conservative inner disk radius compared to our previous results, moving the acceleration region further from the innermost stable circular orbit. We estimate the high-energy neutrino emission from hadronic and photo-hadronic processes, constrained by the acceleration timescale for first-order Fermi acceleration within the turbulent current sheet. The estimated model reproduces the IceCube neutrino flux excess, providing an essential technical complement and validation for our forthcoming comprehensive publication.

The black hole X-ray binary (BHXRB) Swift J151857.0-572147 was discovered during its first outburst in March 2024. We review the archive of NICER observations from this outburst, focusing on the soft states. We select spectra for which the disk to total flux ratio exceeds 0.8 and the coronal scattering fraction fsc is less than 25%, conditions under which the accretion disk is expected to extend to the innermost stable circular orbit (ISCO) and remains geometrically thin. Through a continuum fitting analysis, we explore the dependence of spin on mass, inclination, and distance. We constrain the spin within the parameter space typical of stellar-mass black holes (sBHs), assuming a mass around 10 solar masses, inclination angles between 20 and 60 degree, and distances between 4 and 16 kpc. For fiducial parameters: a mass of 10 solar masses, a distance of 10 kpc, and an inclination angle of 40 degree, a moderate spin of approximately 0.7 is obtained. However, precise determination of the spin will require accurate measurements of these parameters. Our analysis provides a framework to infer the spin and estimate its uncertainties once more precise measurements of mass, distance, and inclination become available. As we demonstrate, lower inclination angles, greater distances, or larger black hole masses result in higher spin values.

Bipolar Light Bridges (BLBs) are bright regions located between sunspot umbrae of opposite magnetic polarity. They are typically characterized by strong magnetic fields and intense flows, which are believed to be closely associated with major solar flares. Despite their importance, their fine structure, formation and evolution remain poorly understood. In this work, we analyze the observations of a well-defined BLB obtained by the Goode Solar Telescope at the Big Bear Solar Observatory and the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. The high-resolution GST observations reveal that the BLB is composed of fine, penumbral filament-like structures with widths of approximately 100-150 km. The corresponding Doppler velocity maps present a stable pattern of spatially adjacent red- and blueshifted patches within the BLB throughout the 5.5-hour GST observation. HMI observations show that the BLB arises from the converging and shearing motions of sunspots with opposite polarities. Penumbral regions originating from different polarities gradually evolve and interact, eventually forming the BLB. The observed Doppler velocity pattern, characterized by red- and blueshifted patches, can be interpreted as a projection effect of the Evershed flow within the penumbrae. Therefore, we argue that the BLB is formed through the compression and stretching of penumbral structures from oppositely polarized sunspots.

Recent obeservations have revealed significant soft X-ray polarizations from several quiescent magnetars, including the intriguing $90^°$ polarization angle (PA) swing as a function of photon energy for some sources. We present a general semi-analytical framework for calculating energy-dependent soft X-ray polarization signatures from magnetars, consistently incorporating both QED vacuum resonance in the atmosphere and resonant Compton scattering (RCS) in the magnetosphere. Starting from the polarized radiative transfer equation for RCS and treating vacuum-resonance-induced mode conversion as an input, we employ a first-order approximation in RCS optical depth to evaluate the effect of different magnetospheric plasma density (which depends on magnetic twist), drift velocity and temperature, and viewing geometry on the observed radiation. Our analysis reveals that magnetic twist and plasma drift velocity are the critical parameters controlling the impact of RCS on both the absolute polarization degree and its variation across the soft X-ray spectrum. We find that sufficiently strong RCS can wash out the PA swing caused by vacuum resonance. Furthermore, in addition to the QED vacuum resonance effect, significant relativistic signatures arising from plasma drift velocity ($\beta_0 \gtrsim 0.5$) may introduce an extra $90^\circ$ PA swing in the spectrum. Our calculation framework, based on single-scattering approximation, bypasses the need for complex, multi-dimensional Monte Carlo simulations, providing an analytical pathway for modeling full-surface emission and rotational-phase-resolved radiation from magnetic neutron stars, in support of current and future X-ray polarization missions.

In this review we outline our knowledge on slow neutron captures, concentrating on its main part occurring during the final stages of stellar evolution for low or intermediate-mass stars when they evolve during the Asymptotic Giant Branch, or AGB, stars. We focus our attention on how, in this field, studies passed from a first era of inquiries based on nuclear systematics, to numerical nucleosynthesis computations performed in stellar codes. We then discuss how these last were forced, by observational constraints, to almost abandon, for the synthesis of nuclei between Sr and Pb, the rather naturally activated Ne22 neutron source (operating efficiently at T > 30 keV, and producing a neutron density N_n > 5 10^8 cm^-3). This implied considering the alternative reaction 13C(alpha,n)16O, that can be activated locally after each of the TDU. The mentioned crucial reaction occurs at T< 8 keV, in the time intervals separating two subsequent thermal pulses (TP). The layers where 13C(alpha,n)16O operates are characterized by a radiative equilibrium and their low temperature also yields low values for the neutron density (N_n < 10^7 cm^-3).

The potential of prebiotic chemistry to propagate on an exoplanet fundamentally depends on whether the atmospheric conditions can facilitate the production of prebiotic feedstock molecules. Photochemical simulations of exoplanet atmospheres can be used to explore this potential atmospheric synthesis, but require a comprehensive chemical network. We present the implementation of the CRAHCN-O network, constructed to simulate the formation of feedstock molecules such as HCN, H$_2$CO, and simple hydrocarbons, into the VULCAN photochemical kinetics code. We investigate the production of feedstock molecules driven by M-star radiation and compare these to predictions by the N-C-H-O network in VULCAN, for N$_2$-dominated atmospheres with C/O ratios between 0.5-1.5. Predicted abundances are similar for C/O${=}$0.5. Once CH$_4$ is included (i.e., for C/O${>}$0.5), the abundance profiles diverge in the photochemical regions. By analysing the attenuation of UV radiation, we find that hydrocarbon photochemical shielding causes the diverging profiles. CRAHCN-O accumulates C$_2$H$_6$, while N-C-H-O accumulates C$_4$H$_3$ and C$_3$H$_4$. Importantly, C$_2$H$_6$ is photochemically active whereas C$_4$H$_3$ and C$_3$H$_4$ are assumed inactive. With mixing ratios up to a few percent in CRAHCN-O, C$_2$H$_6$ shields CH$_4$ and CO$_2$ from photodissociation and weakens the destruction of HCN and H$_2$CO. Maximum HCN mixing ratios reach 1000 ppm with CRAHCN-O compared to only 3 ppm with N-C-H-O. Other feedstock molecules like HC$_3$N and C$_2$H$_2$ form more efficiently in N-C-H-O. The shielding mechanism and its impact on feedstock molecules persist for radiation from distinct M-star types. These results demonstrate the crucial role of chemical kinetics in understanding prebiotic processes in exoplanet atmospheres, including important considerations for the construction and applicability of chemical networks.

We report optical time-resolved photometric observations of a newly-discovered outbursting system, ASASSN-25dc. Its 8-mag amplitude, 40-day duration, 1-mag dip in the outburst plateau, and positive superhumps are characteristic of a dwarf nova superoutburst in a non-magnetic cataclysmic variable. We establish its stage-A and stage-B superhump periods as 0.059387(5) d and 0.058864(3) d, respectively. The negative superhump period derivative ($-$1.4(2)$\times10^{-5}$ cycle$^{-1}$) during the stage-B superhumps and the empirical relation indicate the mass ratio is 0.054(7), below the period bounce range. The long outburst decline timescale (35.2(1) d mag$^{-1}$) and small superhump amplitude ($\simeq$0.08 mag) observed in ASASSN-25dc are also seen in some period bouncer systems, but not seen in systems well before the period bounce. Despite its short superhump period and indicated small mass ratio, we find no evidence of the excitement of the 2:1 tidal resonance. Moreover, its outburst rise timescale (1.62(9) d mag$^{-1}$) is significantly longer than those measured at less than 0.4 d mag$^{-1}$ in other dwarf nova outbursts around the period minimum. Overall, an inside-out dwarf nova outburst from a massive disc in a system with a mass ratio around or even below the period minimum, but lacking the 2:1 tidal resonance, may explain all these observations. However, this challenges the existing models of dwarf nova superoutbursts, which do not predict these outburst properties in low-mass-ratio systems.

We propose that the recently discovered luminous hypersoft X-ray sources can be explained by accretion onto a binary black hole in a hierarchical triple system. For black hole masses $\sim 15 M_\odot$, the orbital separation of the internal binary might be $\sim 0.01 $ AU. If the donor provides $\gtrsim 10^{-8} M_\odot$ yr$^{-1}$, then the circumbinary accretion disc can explain the observed properties of the most luminous supersoft sources.

Event horizon-scale observations of low-luminosity black hole accretion flows favor magnetically arrested disks, characterized by dynamically important magnetic fields ($\beta\lesssim1$, where $\beta$ is the ratio of plasma thermal pressure to magnetic pressure) and a two-temperature transrelativistic plasma. Motivated by plasma conditions in the synchrotron-emitting regions of these models, we perform 2D particle-in-cell simulations of electron-ion plasmas with a realistic mass ratio, subject to continuous compression perpendicular to the mean magnetic field $\boldsymbol{B}_0$. Conservation of particle magnetic moments drives pressure anisotropy $P_{\perp}>P_{\parallel}$, triggering anisotropy-driven instabilities. For ion plasma beta $\beta_{i0}=0.5$ and ion temperature $k_{\text{B}}T_{i0}/m_i c^2=0.05$, the ion pressure anisotropy is regulated by the ion cyclotron instability, while the mirror mode influences the late-time electron anisotropy. Both species develop nonthermal components at high energies, consistent with stochastic acceleration by cyclotron-scale fluctuations. We characterize how the onset and time evolution of the plasma instabilities, as well as the resulting ion and electron anisotropies and energy spectra, vary with $\beta_{i0}$, $k_{\text{B}}T_{i0}/m_i c^2$, electron-to-ion temperature ratio $T_{e0}/T_{i0}$, and the compression rate. Increasing the thermal energy toward relativistic values raises the anisotropy thresholds for all instabilities observed in our simulations, allowing larger anisotropies to develop. For $T_{e0}/T_{i0}<1$, as expected in collisionless two-temperature accretion flows, the growth of mirror and whistler instabilities is delayed or suppressed, leading to increasingly adiabatic evolution of the electrons. Our findings can be used to inform global fluid models of black hole accretion.

As one of the most widely recognized high-quality standard stars, the Stetson standards have been extensively used as a photometric reference for calibrating other surveys. In this work, we present an independent validation and re-calibration of the Stetson standard star photometry using the BEST database. Based on typically 30,000-70,000 calibration stars per band, we find that the original Stetson photometry achieves field-to-field zero-point precisions of approximately 10–40 mmag in the $UBVRI$-band. In addition, significant spatially dependent magnitude offsets are detected within individual Stetson fields for all bands, with magnitudes exceeding 1\%, probably caused by the calibration errors in the Stetson photometry. After correcting those systematic errors, the agreement between the Stetson and BEST photometry is improved to $\sim$5 mmag for individual fields for $BVRI$-band. The re-calibrated photometry is further validated using the SCR standards, yielding agreement better than 10 mmag for individual stars in the $BVRI$ bands and confirming zero-point precisions of 2–4 mmag in the $BVI$ band. The precisions is further confirmed by checks using Gaia DR3 broad band colors. These results highlight the power of the BEST database for improving photometric calibration and suggest that, if feasible, it be incorporated into the calibration process of future releases of the Stetson standard catalog.

There is an ongoing discussion in the literature on the nature of long-period transients (LPTs), radio-emitting sources with periods ranging from hundreds to tens of thousands of seconds. Although some of these objects have been identified as white dwarf (WD) + M-dwarf binaries, this description currently does not fit the entire class. An example is GLEAM-X J162759.5-523504.3 (hereafter GLEAM-X J1627-5235), with a period of 1091 s, for which the lack of an optical counterpart disfavors the presence of such a binary system. In this case, GLEAM-X J1627-5235 could be interpreted as an isolated, massive, fast-rotating, and highly magnetized (  1e+9 G) WD pulsar. Its properties are consistent with a carbon-oxygen WD of mass  1.3 Msun and radius  2500 km, possibly supported by small-scale multipolar magnetosphere structures that keep it above the death line for WD-pulsars. We assess a double WD merger origin, modeling the post-merger rotational evolution under accretion, propeller, and magnetic braking torques. We find rotational age of  572 Myr for GLEAM-X J1627-5235, i.e., the post-merger time required to reach its observed period. This result is consistent with current optical upper limits for GLEAM-X J1627-5235 and support the WD pulsar interpretation for this source. We also discuss how the same model can apply to other LPTs.

Eclipsing close double white dwarf (WD) systems provide a unique opportunity to directly constrain hydrogen-envelope retention and test common-envelope (CE) evolution in low-mass stars, since they allow precise determinations of stellar masses and radii. We analyze J2102-4145, an eclipsing binary composed of two low-mass helium-core white dwarfs in a 2.4-hour orbit. By comparing the observed radii and effective temperatures with updated evolutionary models for CE evolution and stable Roche-lobe overflow (SRLOF), we confirm that both stars are helium-core white dwarfs. The primary, with a mass of 0.375 solar masses, is consistent with SRLOF models that retain thick hydrogen envelopes and sustain residual nuclear burning, whereas the secondary, with a mass of 0.314 solar masses, can only be reproduced by CE models in which the hydrogen envelope is almost completely removed. The inferred cooling ages (approximately 220 Myr for the secondary and between about 260 and 510 Myr for the primary, depending on the contribution of residual nuclear burning) support a formation sequence in which the primary formed first through SRLOF, followed by a CE phase that produced the compact secondary. Reconstruction of the CE energy budget yields progenitor and orbital parameters consistent with this evolutionary picture. The unusually small radius of the secondary requires an extremely thin hydrogen envelope, with a mass below about 10e-7 solar masses, well below the values predicted by standard bifurcation criteria. J2102-4145 therefore provides one of the strongest observational constraints on the hydrogen-envelope mass of post-CE low-mass white dwarfs and represents a benchmark challenge for current prescriptions of envelope ejection.

On the Sun, the inertial mode with the largest observed amplitude (rms velocity exceeding $10$ m/s) is the high-latitude mode with longitudinal wavenumber $m=1$. In two dimensions, on the sphere, linear theory predicts that this mode is unstable due to a shear instability associated with latitudinal differential rotation (fast equator, slower polar regions). We investigate the evolution of this instability numerically and theoretically. The nonlinear vorticity equation is solved using direct numerical simulations in the time domain. The only control parameter is the Ekman number $E$. For $10^{-3}\lesssim E< E_c \approx 1.5\times10^{-3}$, only the high-latitude $m=1$ mode is unstable. We extract its saturation amplitude as a function of $E$ and compare the results with predictions from two perturbative approaches in nonlinear stability theory. The simulations reveal a supercritical Hopf bifurcation. Near onset, the mode amplitude is well described by the Landau equation $d|A|/dt=\sigma_I |A|+\beta_I |A|^3$, with a positive linear growth rate $\sigma_I$ and a negative nonlinear coefficient $\beta_I$. The coefficient $\beta_I$ depends weakly on $E$, implying that the saturated amplitude scales approximately as $|A|\propto\sigma_I^{1/2}$. The equilibrium mode contains the $m=1$ fundamental and harmonics $m=2$ and $m=3$, whose amplitudes scale as $\sigma_I^{m/2}$. Saturation results from Reynolds stresses that smooth the latitudinal differential rotation. For $E=4\times10^{-4}$, consistent with solar-like turbulent viscosity, the saturated velocity reaches $28$ m/s, comparable to solar observations. These results should be interpreted cautiously, since in three dimensions the instability is baroclinic and involves different physics.

Most long gamma-ray bursts (GRBs) exhibit afterglows broadly consistent with external forward-shock emission, typically described by smooth broken power-law decays in the multiband light curve. However, a minority of well-sampled GRBs deviate from this behavior, including GRB 250129A. This object shows multiple late-time rebrightenings at X-ray and optical wavelengths. Rebrightenings are often attributed to energy injection from prolonged central engine activity, refreshed shocks from delayed shell collisions, density jumps in the ambient medium, or angular jet structure and viewing-angle effects. After analysing the prompt emission observed in gamma-rays and the near-infrared, we investigate the origin of X-ray and optical flaring episodes in GRB 250129A. Physical processes in the afterglow light curves were investigated using methods ranging from empirical fitting to Bayesian inference. The well-sampled flares and the connection between the prompt and afterglow emission allow us to test the consistency of the fireball model and alternative scenarios. Conducting the prompt and time-resolved analyses, we obtained an isotropic-equivalent energy of E_iso,gamma = (1.35 +/- 0.12) x 10^53 erg. By modeling the afterglow using an agnostic Bayesian framework (NMMA), we rule out both a single external-shock evolution and a one-time energy-injection scenario. Numerical calculations show that the rebrightening episodes are consistent with refreshed shocks from delayed collisions between relativistic shells. Based on the consistency between our analyses of the prompt and afterglow GRB 250129A data, we find that two statistically significant rebrightening episodes occur within 1.1 days post trigger and can be explained by a sequence of refreshed shocks. Temporally and spectrally rich GRB datasets such as the one presented in this work, provide a powerful means to test current modeling frameworks.

We study the full evolution of a 1.313 Msun white dwarf star that descended from a 9 Msun main-sequence progenitor with an initial metallicity of Z=0.02. Using MESA r24.08.01, we calculate its entire evolution from pre-ZAMS to the WD cooling curve, including both the evolution through 139 thermal pulses and the post-AGB phase. The resulting remnant is an ultramassive H-deficient WD, for which the composition, in mass fraction, is 47.7% O16, 39.7% Ne20, 4.2% Mg24, 3.3% Na23 and 0.386% C12 – corresponding to a total mass of 5 x 10^-3 Msun of C –, surrounded by a 1.5 x 10^-5 Msun He layer. We also investigate the effects of fully suppressing the TP-SAGB stage by adopting a high mass-loss rate only after the second dredge-up, and find only minor differences in the final mass and composition. In addition, we calculate models with and without phase separation during the WD stage, estimating a cooling delay of only 16 Myr. This is the first ultramassive white dwarf sequence for which both the TP-SAGB and post-AGB stages are calculated and, to our knowledge, the most massive WD model from complete evolution for which cooling times and detailed abundance profiles are published

The observation of high-energy neutrinos from the direction of the nearby active galaxy, NGC 1068, was a major step in identifying the origin of high-energy cosmic neutrinos. The multimessenger data imply that high-energy neutrinos originate from the hearts of active galaxies which are opaque to GeV-TeV $\gamma$-rays. This realization is reinforced by the excess of neutrinos in the direction of NGC 4151 and Circinus Galaxy, other nearby active galactic nuclei (AGNs). Understanding the vicinity of supermassive black holes with electromagnetic radiation is often challenging due to uncertainties associated with the absorption of emission in these dense environments, and neutrinos can be used as a powerful probe of the inner parts of the active galaxies. Considering the five brightest neutrino-active galaxies, NGC 1068, NGC 4151, CGCG 420-15, Circinus Galaxy, and NGC 7469, we employ the measured neutrino spectra together with the sub-GeV $\gamma$-ray emission measured by the {\em Fermi} satellite to break the degeneracy and narrow in on the parameter space of neutrino emission from turbulent coronae of AGNs. We also study contributions of jet-quiet AGNs, whose properties are similar to NGC 1068 and NGC 7469, to the isotropic neutrino background flux, through exploring possibilities that the neutrino luminosity function may deviate from the X-ray luminosity function. Our results will help estimate the prospects for identifying additional neutrino-active galaxies and guide future targeted analyses.

Context. Ultra fast outflows (UFOs) from active galactic nuclei (AGN) are expected to lead to the formation of sub-relativistic strong shocks expanding in a dense circumnuclear medium, and thus have the potential for being efficient particle accelerators, and to be proficient sources of gamma rays and neutrinos. Aims. We investigate the detectability of a sample of nearby identified UFOs in gamma rays and neutrinos with current and next- generation instruments. Methods. We model the acceleration of particles at the strong shocks of UFOs, and estimate the associated gamma-ray and neutrino signal. We adopt our model to investigate the prospects for detection with current and next-generation observatories. Results. We find that several UFOs could be detectable in the very-high-energy (VHE) domain - for example, by the Cherenkov Telescope Array Observatory (CTAO)- even if they remain undetected by Fermi-LAT in the high-energy range. Detectability is favored for hard proton spectra (spectral index {\alpha} \lessim 3.9), high acceleration efficiencies, and amplified magnetic fields. Our results suggest that next-generation VHE observatories could detect the first gamma-ray signatures of AGN UFOs, providing a new probe of particle acceleration in sub-relativistic shocks

Spectral-sirens inference enables the extraction of cosmological parameters from gravitational-wave data alone, without electromagnetic counterparts or galaxy catalogs. We introduce new parametric mass functions for the binary black hole population that capture significant structure across the mass spectrum and are moderately favoured by Bayesian evidence over simpler models. Analysing the latest gravitational-wave transient catalog, GWTC-4.0, we show that powerlaws-only population models constrain the Hubble constant to $23\%$ precision, $H_0 = 53.3^{+14.0}_{-10.8} ~\rm km \,s^{-1} \,Mpc^{-1}$ at $68\%$ confidence level. This represents a $\sim 50\%$ improvement over the corresponding binary black hole-only analysis by the LIGO-Virgo-KAGRA collaboration, achieving precision comparable to their joint analyses including neutron stars and galaxy catalogs. We further test alternative cosmological models, establishing competitive constraints on modified gravitational-wave propagation, while bounds on the dark energy equation-of-state parameters remain uninformative. Projecting to future O5 observing run, we forecast substantial improvements in $H_0$ and modified propagation parameters with larger datasets at higher redshifts. Our results highlight the strong interplay between the black hole mass distribution and inferred cosmology, underscoring the need for suitable population models to fully exploit gravitational-wave data.

This study presents the trajectory design for a mission touring Saturn's Inner Large Moons (Rhea, Dione, Tethys, Enceladus, and Mimas) engineered to meet observational requirements, including full surface coverage, while ensuring low fuel consumption and compatibility with current power and propulsion technologies (radioisotope thermoelectric generators and Hall effect thrusters). The tour begins at Rhea and ends at Mimas, using a trajectory concept that alternates between extended observation phases around each moon and Saturn centered low-thrust spiral arcs to transition efficiently to the next target. The J2-perturbed Circular Restricted Three-Body Problem is adopted to design exploration paths, with halo orbits serving as staging points for heteroclinic and homoclinic loops that enable prolonged, repeated, and comprehensive surface reconnaissance (including critical regions such as Enceladus poles, where geological activity produces intense plumes). Stable and unstable hyperbolic invariant manifolds of the halo orbits act as departure and arrival gateways for propelled inter-moon transfers, modeled in an ephemeris-based framework including gravitational perturbations from the moons, the Sun, and Saturn's oblateness. The dynamical model setup is guided by a rigorous perturbation analysis to maximize computational efficiency while maintaining high fidelity trajectory design. A locally optimal guidance law minimizes propellant consumption. The proposed tour offers an alternative to traditional flyby missions, providing comparable total duration but greater observing time and reduced fuel requirements, and advances previous work by achieving both complete lunar surface coverage and high-fidelity modeling.

We adopt a general relativistic ray-tracing approach to study the shadows and polarization images of spherically symmetric Gauss-Bonnet (GB) black holes enveloped by geometrically thick accretion flows. Specifically, we adopt a phenomenological RIAF-like model and an analytical Hou disk model. In the RIAF-like model, increasing the GB coupling parameter $\lambda$ reduces both the size and brightness of the higher-order image, while increasing $\theta$ alters the shape of the higher-order image and obscures the horizon's outline. The main difference between isotropic and anisotropic emission is that the latter produce distortion of the high-order image in the vertical direction, leading to an elliptical morphology. For the Hou disk model, due to specific regions being geometrically thinner with the conical approximation, the high-order images are narrower with the increase in $\lambda$ than the RIAF model. While increasing $\theta$ enhances the brightness of the direct images outside the higher-order images, but hardly changes the size of the higher-order images, which is in sharp contrast to the RIAF model. Meanwhile, the Hou disk produces polarization patterns that trace the brightness configuration and are affected by $\lambda$ and $\theta$, reflecting the intrinsic structure of spacetime. These results illustrate that intensity and polarization in thick-disk models provide probes of GB black holes and near-horizon accretion dynamics.

Previous work on the gravitational-wave background generated in a two-scalar-field cosmological model, in which warm inflation and the dark sector are unified within a single framework, has shown that the resulting spectrum could be potentially detectable by planned next-generation gravitational-wave observatories. In this work, we extend this analysis to the weak dissipation regime of warm inflation, highlighting how the features of the inflationary scenario play a crucial role in the production of gravitational waves. The full gravitational-wave energy spectrum is calculated using the formalism of continuous Bogoliubov coefficients. By comparing our results with those obtained in the strong dissipation regime and with the sensitivity curves of future detectors, we find that the weak dissipation regime improves the prospects for observational detection.

Energy dissipation and particle acceleration in the collisionless magnetotail plasma remain incompletely understood. While Kinetic Alfvén Waves (KAWs) are widely hypothesized to mediate these processes, observational characterization of their spectral properties and dissipation signatures in magnetotail boundary layers remains limited. We report observations of KAW turbulence and parallel electric fields ($E_{\parallel}$) in the outer Plasma Sheet Boundary Layer (PSBL) using high-resolution burst-mode data from the Magnetospheric Multiscale (MMS) mission. For a crossing event on May 31, 2017, we identify broadband KAW turbulence characterized by a normalized electric-to-magnetic field ratio $\mathcal{R} = |\delta E_{\perp}|/(v_A|\delta B_{\perp}|) = 2.5 \pm 1.2$ exceeding the MHD limit, a spectral break near ion scales, a steep kinetic-range spectral slope ($\alpha = -3.48 \pm 0.13$), and low magnetic compressibility ($C_{\parallel} \approx 0.03$). We observe impulsive parallel electric field structures (up to 15 mV/m) and large-amplitude density fluctuations (up to 68\%) during intervals of enhanced wave activity. The steep spectral slope, steeper than theoretical predictions for undamped KAW cascades ($-7/3$ to $-8/3$), is consistent with substantial energy removal from the cascade at kinetic scales. The near-zero correlation between the $E_{\parallel}$ waveform and density fluctuations ($r \approx -0.03$) suggests that the observed $E_{\parallel}$ structures are not straightforwardly organized by compressive density variations, consistent with dissipation through direct wave–particle interaction. Attribution to a specific damping channel (e.g., Landau damping) is not uniquely constrained by the present diagnostics. These observations support collisionless damping of KAW turbulence at kinetic scales in the intermediate-beta, outer PSBL of the terrestrial magnetotail.

We show that the hidden symmetry and separability characteristic of Kerr geometry arise as an inevitable \emph{local} consequence of the Einstein equations for rotating spacetimes. Without assuming separability, algebraic speciality, Killing-Yano symmetry, or global boundary conditions, we analyze general stationary and axisymmetric geometries in a locally non-rotating orthonormal frame. Imposing a minimal physical requirement, the local equilibrium condition, we find that the mixed Einstein equations enforce a rigid projective alignment between the radial and angular sectors. This result does not rely on vacuum: the mixed equations are insensitive to the detailed form of a comoving stress-energy tensor. Consistency then requires equality of Schwarzian derivatives, leading to a universal classification of local solutions into Möbius, exponential, and trigonometric branches. Global regularity generically excludes the trigonometric branch, leaving precisely the Kerr-type sector and the emergence of Petrov type D structure. This provides a structural precursor to Kerr uniqueness, demonstrating that the kinematical core of Kerr geometry is encoded locally within the Einstein equations.

Black holes are the most efficient scramblers in nature. By mapping the instantaneous mass and angular momentum of two spinless black holes in a quasi-circular binary onto those of an effective Kerr black hole, we demonstrate that the final state of the merger remnant corresponds with remarkable accuracy to the configuration that renders null geodesics unstable at the highest possible rate. This suggests a deep connection between the properties of black holes resulting from binary mergers and their unstable null orbits.

We derive new constraints on light CP-even scalars using old gamma-ray observations in the direction of SN1987A by the Solar Maximum Mission (SMM) satellite. Light scalars can be abundantly produced in the supernova core via the nucleon bremsstrahlung process, can stream out of the supernova-environment and decay into photons – either primary photons or secondary photons from lepton-antilepton pairs – thus leading to a gamma-ray signal. From the non-observation of excess photon flux by SMM after the detection of the neutrino burst from SN1987A, we set new constraints on the mixing angle of the CP-even scalar with the Standard Model Higgs boson.

The slopes of the linear relations between sunspot and white light (WL) facular areas at the onset of sunspot Cycles 12-21 correlate well with the amplitudes of those cycles between 1878-1980 (Brown and Evans, 1980). We use continuum images from the SOHO Michelson Doppler Imager and SDO Heliospheric Magnetic Imager to show that the relation holds also for Cycles 24 and 25. The amplitudes of Cycles 12-21 and 24 calculated using this relation agree with the observed amplitudes to within +/- 4% rms. It also enabled us in 2022 to correctly predict a larger Cycle 25 than estimated by the International Prediction Panel, 3 years before maximum. The technique offers an objective, physically based predictor of cycle amplitudes 3-4 years ahead of their maxima, given a stable source of continuum full disk photospheric images.

Among more than 1000 known fast radio bursts (FRBs), only five sources - FRBs 20121102A, 20190520B, 20201124A, 20240114A and 20190417A - have confirmed associations with persistent radio sources (PRS). The observed quasi-steady emission is consistent with synchrotron radiation from a composite of magnetar wind nebula (MWN) and supernova (SN) ejecta. Using a phenomenological model that incorporates simplified treatments of the nebular dynamics and particle acceleration, we compute the synchrotron flux by solving kinetic equations for energized electrons, accounting for electromagnetic cascades of electron-positron pairs interacting with nebular photons. Within the framework of our model, the rotation-powered scenario requires a young neutron star (NS) with age $t_{\rm age}\approx 20\,{\rm yr}$, dipolar magnetic field $B_{\rm dip}\approx (3-5)\times10^{12}\,{\rm G}$ and initial spin period $P_i\approx 1.5-3\,{\rm ms}$ in an ultra-stripped SN progenitor to account for emissions from FRBs 20121102A and 20190520B. In contrast, FRB 20201124A requires $t_{\rm age}\approx 10\,{\rm yr}$, $B_{\rm dip}\approx 5.5\times10^{13}\,{\rm G}$ and $P_i\approx 10\,{\rm ms}$ in a conventional core-collapse SN progenitor. For the magnetar-flare-powered model, NS aged $t_{\rm age} \approx 25\,/40\,{\rm yr}$ in a USSN progenitor and $t_{\rm age} \approx 12.5\,{\rm yr}$ in a CCSN progenitor explains the observed flux for FRB 20121102A/20190520B and FRB 20201124A, respectively. Finally, we estimate a minimum NS age $t_{\rm age,min} \sim 1-3\,{\rm yr}$ based on the near-source plasma contribution to observed DM, and $t_{\rm age,min} \sim 6.5-10\,{\rm yr}$ from the absence of radio signal attenuation.

A Chern-Simons interaction between a pseudo-scalar field and a U(1) gauge field results in the generation of a chiral gravitational wave background. The detection of this signal is contrasted by the fact that this coupling also generates primordial scalar perturbations, on which strong limits exist, particularly at CMB scales. In this study, we propose a new extension of this mechanism characterized by a non-canonical kinetic term for the pseudo-scalar. We find that a decrease of the sound speed of the pseudo-scalar field highly suppresses the sourced scalar with respect to the sourced tensor modes, thus effectively allowing for the production of a greater tensor signal. Contrary to the case of a canonical axion inflaton, it is in this case possible for the sourced tensor modes to dominate over the vacuum ones without violating the non-Gaussianity constraints from the scalar sector, which results in a nearly totally polarized tensor signal at CMB scales. We also study the extension of this mechanisms to the multiple field case, in which the axion is not the inflaton.

The cosmological Mass Varying Neutrino (MaVaN) model is considered, where the interaction between a fermionic field and a scalar field with a Ratra-Peebles potential via a Yukawa coupling is investigated. Observational constraints on the flat and non-flat MaVaN models, as well as on the standard $\Lambda$CDM model, are derived from 32 $H(z)$ measurements using MCMC analysis. Comparison with the $\Lambda$CDM model using the criteria $\Delta\chi^2_{\rm min}$, $AICc$, and $BIC$ shows no statistically significant improvement for MaVaN models. Deviations in the expansion history remain well below $1\sigma$, indicating that the $H(z)$ data alone do not provide sufficient constraining power to distinguish MaVaN models from the $\Lambda$CDM model. The non-flat MaVaN model reduces the tension between the $H_0$ value inferred from $H(z)$ data and the Planck CMB measurement from $\sim 2\sigma$ in the $\Lambda$CDM framework to $\sim 1.1\sigma$. The discrepancy with the SH0ES measurement of $H_0$ is also reduced to below $1\sigma$, primarily due to the large uncertainties of the $H(z)$ data.

To mitigate the severe information loss arising from widely adopted linear scale cuts in constraints on modified gravity parameterisations with Weak Lensing (WL) and Large-Scale Structure (LSS) data, we introduce a novel alternative method for data reduction. This Principal Component Analysis (PCA)-based framework extracts key features in the matter power spectrum arising from nonlinear effects in a set of representative gravity theories. By performing the analysis in the space of principal components, we can replace sweeping `linear-only' scale cuts with targeted cuts on the transformed data vector, ultimately reducing parameter bias and significantly tightening constraints. We forecast constraints on a minimal parameterised extension to $\Lambda$CDM which includes modifications to the growth of structure and lensing of light ($\Lambda$CDM$+\mu_0+\Sigma_0$) using mock Stage-IV data for two simulated cosmologies: the $\Lambda$CDM model and Extended Shift Symmetric (ESS) gravity. Under the assumption of a Universe defined by $\Lambda$CDM and General Relativity, our method offers constraints on $\mu_0$ a factor of 1.65 tighter than traditional linear-only scale cuts. Crucially, our approach also provides the necessary constraining power to break key degeneracies in modified gravity without relying on $f\sigma_8$ measurements, introducing a promising new tool for the analysis of present and future WL and LSS photometric surveys.

Diffusive shock acceleration and diffusion propagation are essential components of the standard cosmic ray model. These theories are based on extensive observations of high-energy solar processes, providing substantial direct evidence in the MeV energy range. Although the model is widely and successfully used to explain high-energy cosmic phenomena, direct validation has been elusive. The multi-wavelength spectra and angular profile measurements of the Geminga pulsar wind nebula and its pulsar halo, particularly the precise spectral observations by HAWC and LHAASO-KM2A in recent years, offer a rare opportunity to test these theories with cosmic rays energies between 1 TeV and several hundred TeV. These observations are expected to elevate the direct testing of theoretical models from multi-MeV to sub-PeV energies. In this work, a method is developed to test the diffusive shock acceleration and diffusion propagation model between one and several hundred TeV energies through the latest spectral and morphological data of the Geminga region from HAWC and Fermi-LAT. Our results show that the theories of diffusive shock acceleration and diffusion propagation are consistent with experimental observations. However, the published morphological data adopted rather wide energy bins and currently do not allow a high precision test of the inferred energy dependent diffusion coefficient by observed energy spectra with DSA theory. It is anticipated that future HAWC and LHAASO-KM2A observations will yield higher-precision results, and the confirmation of a rapidly increasing diffusion coefficient above 100 TeV would serve as important evidence supporting the diffusive shock acceleration and diffusion propagation theory. Similar tests would be both important and valuable for other models.

Gravitational waves (GWs) distort galaxy shapes through the tidal effect, offering a novel avenue to probe the nature of gravity. In this paper, we investigate how extra GW polarizations beyond those predicted by general relativity imprint observable signatures on galaxy shapes. Since galaxy shapes are measured as two-dimensional images projected onto the celestial sphere, we present three-dimensional statistical quantities of the projected tidal tensor sourced by the tensor perturbation. We show that the presence of extra polarization modes modifies both the amplitude and angular dependence of the correlation functions. Furthermore, we identify a distinct observational channel for probing parity violation in helicity-two and helicity-one modes. In particular, we show that if they propagate at different speeds, galaxy surveys can disentangle the source of parity violation. Our findings establish a theoretical framework for using upcoming large-scale galaxy surveys to test modified gravity theories through the polarization content of GWs.

Nanohertz gravitational waves (GWs) from supermassive binary black holes (SMBBHs), detectable via pulsar timing arrays (PTAs), offer a novel avenue to constrain dark energy. Based on cosmological simulations and semi-analytic galaxy formation models, this study explores the detectability of individual nanohertz SMBBH sources using next-generation PTAs and their potential for constraining dark energy under an optimistic scenario considering only the presence of white noise. By constructing light-cone SMBBH populations across hardening timescales ($\tau_H = 0.1/5/10$Gyr) and computing signal-to-noise ratios (SNR), we find advanced PTAs can resolve $10^2$–$10^3$ sources with SNR $> 8$ (primarily at $z < 1$ with chirp masses of $10^8$–$10^{10}M_{\odot}$). If electromagnetic counterparts can be identified, optimal configurations ($\sigma_t = 50$ns, $N_p = 1000$, $T_{\text{obs}} = 30$yr with$ \tau_H \leq 5$Gyr) could constrain the dark energy equation-of-state (EoS) parameter $w$ to $\Delta w \sim 0.023$–$0.048$, where the constraints only exhibit weak dependence on $\tau_H$ within $0.1$–$5$Gyr. If only $10\%$ of GW sources have detectable electromagnetic counterparts, constraints weaken to $\Delta w = 0.075$ ($\tau_H = 0.1$Gyr) and $\Delta w = 0.162$ ($\tau_H = 5$Gyr) under the most optimal parameter configuration. What's more, conservative PTAs ($N_p = 500$, $\sigma_t = 100$–$200$ns) with additional $30$-year data accumulation could double resolvable source counts and improve $\Delta w$ precision by $\sim 40\%$.

The excess of muons observed in ultra-high-energy cosmic-ray air showers relative to simulation predictions, known as the muon puzzle, provides indirect evidence of our incomplete understanding of high-energy hadronic interactions. An unambiguous resolution requires that each proposed solution be directly tested through cosmic-ray and collider experiments probing hadronic interactions. In this work, we develop a framework to assess the strangeness enhancement scenario, wherein an increased yield of kaons relative to pions boosts muon production, which can connect a model prediction and cosmic-ray and collider measurements. Using the \textsc{MCEq} air-shower simulation package, we first identify the key phase-space regions of hadronic interactions that drive muon yields in this scenario. Subsequent analysis demonstrates that a strangeness enhancement starting at $10^6-10^7~\mathrm{GeV}$ can consistently explain the latest cosmic-ray experiments and requires substantial enhancement at the Large Hadron Collider (LHC) energy. Furthermore, evaluating the required precision for LHC measurements, assuming Pierre Auger Observatory muon measurements and forthcoming kaon-to-pion ratio data from LHC Run 3, reveals that these experiments can robustly constrain the majority of the scenario's parameters. In particular, achieving 10.8\% precision on the kaon-to-pion ratio at LHCb and 8.4\% at FASER is sufficient to test the strangeness enhancement scenario over its viable parameter space. These upcoming experimental results will provide the first direct constraints on strangeness enhancement as a potential resolution of the muon puzzle.

The hot gas in the outskirts of galaxy cluster-sized halos, extending around and beyond the virial radius into nearby accretion regions, remains among one of the least explored baryon components of the large-scale cosmic structure. We present a stacking analysis of 680 galaxy clusters located in the western Galactic hemisphere, using data from the first two years of the SRG/eROSITA All-Sky Survey. The stacked X-ray surface brightness (SB) profile reveals a statistically significant signal extending out to 2 r_200m ( 4.5 Mpc). The best-fit SB profile is well described by a combination of terms describing orbiting and infalling gas, with a transition occurring around r_200m. At this radius, the gas number density corresponds to a baryon overdensity of about 30. By integrating the density profile out to r_200m, we inferred a gas fraction exceeding the universal baryon fraction, assuming a typical halo concentration. However, correcting for possible clumping effects reduces the baryon fraction by more than 20%. Additionally, we examined the distribution of hot gas in massive clusters in the IllustrisTNG simulations, from the halo center to the accretion region. This analysis reveals differences in radial gas profiles depending on whether the direction points toward voids or toward nearby cosmic filaments. Beyond r_200m, the density profile along the filament direction exceeds that along the void direction. This pattern aligns with the observed transition radius between the one-halo and two-halo terms, suggesting that r_200m is the approximate radius marking the location at which cosmic filaments connect to galaxy clusters. Meanwhile, comparisons of the gas density and gas fraction profiles between the observation and the IllustrisTNG simulation suggest that the feedback processes in the stacking sample are more efficient at distributing gas to large radii than the IllustrisTNG model.

Next-generation large-scale structure spectroscopic surveys will probe cosmology at high redshifts $(2.3 < z < 3.5)$, relying on abundant galaxy tracers such as Ly$\alpha$ emitters (LAEs) and Lyman break galaxies (LBGs). Medium-band photometry has emerged as a potential technique for efficiently selecting these high-redshift galaxies. In this work, we present clustering analysis of medium-band selected galaxies at high redshift, utilizing photometric data from the Intermediate Band Imaging Survey (IBIS) and spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We interpret the clustering of such samples using both Halo Occupation Distribution (HOD) modeling and a perturbation theory description of large-scale structure. Our modeling indicates that the current target sample is composed from an overlapping mixture of LAEs and LBGs with emission lines. Despite differences in target selection, we find that the clustering properties are consistent with previous studies, with correlation lengths $r_0\simeq 3-4\,h^{-1}$Mpc and a linear bias of $b\sim1.8-2.5$. Finally, we discuss the simulation requirements implied by these measurements and demonstrate that the properties of the samples would make them excellent targets to enhance our understanding of the high-$z$ universe.

Without making any assumption on the underlying geometry and metric of the local Universe, we provide a measurement of the expansion rate fluctuation field using the Cosmicflows-4 and Pantheon+ samples in the redshift range $0.01 < z < 0.1$ ($30 \,h^{-1}\,\mathrm{Mpc} < R < 300\,h^{-1}\,\mathrm{Mpc}$). The amplitude of the anisotropic fluctuations is found to be of order a few percent relative to the monopole of the expansion rate. We further decompose the expansion rate fluctuation field into spherical harmonic components and analyze their evolution with redshift across the studied redshift range. At low redshift, the dipole is clearly dominant, with an amplitude of $\sim (2.2 \pm 0.15)\times 10^{-2}$, significantly larger than the higher–order modes. As redshift increases, the dipole amplitude steadily decreases, reaching roughly half its value in the highest redshift bin investigated. The quadrupole is also significant, at about half the dipole amplitude, and persists across all redshift bins, with no clear decreasing trend, although uncertainties grow at higher redshift. A nonzero octupole is also detected at low redshift. The dipole, quadrupole, and octupole components are found to be aligned, exhibiting axial symmetry around a common axis ($l = 295^\circ,\, b = 5^\circ$). We interpret the observed fluctuations in the expansion rate within the framework of covariant cosmography. Our results indicate that the multipoles of the expansion rate fluctuation field are primarily driven by a strong quadrupole in the covariant Hubble parameter, together with dipole and octupole contributions from the covariant deceleration parameter. These few parameters suffice to reconstruct the luminosity distance with high precision out to $z \sim 0.1$, in a manner that is model–independent, non–perturbative, and free from assumptions about peculiar velocities.

Recent observational evidence of axially symmetric anisotropies in the local cosmic expansion rate motivates an investigation of whether they can be accounted for within the Lemaître-Tolman-Bondi (LTB) framework with an off-center observer. Within this setting, we compute the exact relativistic luminosity distance via the Sachs equation and compare it with the approximate expression obtained from the covariant cosmographic approach (including Hubble, deceleration, jerk and curvature parameters). This comparison allows us to identify the regimes in which the covariant cosmographic method remains reliable. In addition, we compare the LTB relativistic distance for small inhomogeneities with the corresponding result derived from linear perturbation theory (LPT) in the standard cosmological model. This analysis establishes a precise correspondence between the LTB and LPT approaches, offering a consistent dictionary for the interpretation of the observed anisotropies of the large-scale gravitational field. This analysis will be instrumental in interpreting expansion-rate anisotropies, facilitating investigations of the local Universe beyond the FLRW framework with a fully non-perturbative metric approach.

Hydrodynamic escape can strip the envelopes of close-in exoplanets, but most observations of atmospheric mass loss to date have been confined to planets orbiting K and M dwarfs. A growing body of detections of atmospheric escape from planets orbiting early-type stars indicates that they may have significantly stronger and more extended outflows than planets orbiting cooler stars. However, it is unclear whether this limited sample of planets is representative of all gas giants orbiting early-type stars. Motivated by this question, we initiated the first dedicated survey of atmospheric escape from gas giants orbiting F stars in order to understand how their distinct radiation environments shape planetary outflows. We observed ten transits of six planets in an ultra-narrowband filter centered on the metastable helium line using Palomar/WIRC. We report strong ($>3\sigma$) detections of atmospheric escape for WASP-12 b and WASP-180 A b, tentative ($>2\sigma$) detections for WASP-93 b and HAT-P-8 b, and non-detections for WASP-103 b and KELT-7 b. We fit these measurements with a 1D Parker wind model to derive corresponding mass-loss rates, and combine our results with literature measurements to obtain an updated picture of mass loss from planets orbiting early-type stars. Our results indicate that the observed variation in mass-loss rates can be explained by a combination of Roche filling factor and XUV luminosity, and disfavors NUV-driven escape models.

The Polarimetric and Helioseismic Imager (SO/PHI) onboard Solar Orbiter is a spectropolarimeter scanning the Fe I line at 617.3 nm, providing data of the solar photosphere. The same line is sampled by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) and many other on-ground instruments. In this paper, we aim at assessing the consistency between line-of-sight (LoS) velocity measurements from the two instruments. Reliable measurements of up and down flows from SO/PHI are crucial and unique when Solar Orbiter is facing the far side of the Sun. Also, a combination of measurements from two vantage points to study horizontal flows must rely on consistent observations. For this purpose, we compare the LoS velocity measured by SO/PHI's High Resolution Telescope (SO/PHI-HRT) and SDO/HMI on 29 March 2023, when Solar Orbiter was crossing the Sun-Earth line at a distance of 0.39 au from the Sun. Because such co-alignments are rare, this configuration offered an almost unique opportunity to directly compare data products from both telescopes. The data are aligned and remapped to allow a pixel-by-pixel comparison of the whole time series of 4 hours length. Temporal and spatial variations are considered for a direct combination of the measurements. The LoS velocity distributions are evaluated and a clear linear relation is found between the two instruments with a slope of 0.96 and a correlation of 92%. We find that the signals form at similar heights, with a separation of 7$\pm$14 km, which is larger than previous estimates. A close-up look at the penumbra of a sunspot and its Evershed flow is presented. We conclude that the signals inferred by SO/PHI-HRT and SDO/HMI show very good agreement and high correlation when instrumental effects and large-scale velocities on the Sun are properly accounted for.

We propose a predictive $SO(10)$ Grand Unified Theory (GUT) framework for cosmic inflation in the Palatini $\mathcal{R}^2$ formulation of gravity. In this model, a GUT Higgs field both drives inflation and induces intermediate-scale symmetry breaking, thereby linking primordial cosmology, gauge unification, and topological defect formation. A partial inflationary phase of $N_I \sim 10$–$17$ $e$-folds following monopole formation can dilute magnetic monopoles to abundances $Y_M \sim 10^{-35}$–$10^{-27}$. The model yields Cosmic Microwave Background (CMB) predictions of $0.955 \lesssim n_s \lesssim 0.974$, accommodating the tensions between Planck-BICEP ($n_s \approx 0.965$) and Planck+ACT ($n_s \approx 0.971$) via $\phi < M$ and $\phi > M$ branches repectively. The predicted tensor-to-scalar ratio $r \lesssim 8\times10^{-4}$ lies within current observational constraints and is accessible to forthcoming experiments, including the Simons Observatory and LiteBIRD. The resulting correlations between the unification scale $M_U$, the inflationary observables $(n_s, r)$, and proton-decay lifetimes highlight a complementarity between CMB measurements and proton-decay searches, with regions of parameter space testable in forthcoming experiments such as Hyper-Kamiokande and DUNE.

We investigate how the cosmological Euler equation can be tested in the presence of viscous dark matter, violations of the equivalence principle (EP), and modifications of gravity, while relying on minimal theoretical assumptions. Extending the previous analysis, we generalize the observable $E_P$, which quantifies EP violation, to $\tilde{E}_P$, discuss the degeneracy between bulk and shear viscosities and EP-violating effects, and explicitly show that the EP can still be tested in the small-viscosity limit. In addition, we identify a model-independent observable, $C_{\rm vis,0}$, which characterizes the present-day dark matter viscosity and can be measured from relativistic galaxy number counts by cross-correlating two galaxy populations. We perform forecasts for three forthcoming Stage-IV surveys: DESI, Euclid, and SKA Phase 2 (SKA2), and find that $C_{\rm vis,0}$ can be tightly constrained, at the level of $\mathcal{O}(10^{-6})$ or better in all cases. Among these surveys, SKA2 provides the tightest constraint, with a $1\sigma$ uncertainty of $1.08 \times 10^{-7}$ on $C_{\rm vis,0}$.

It is usually taken for granted that mutual synchronisation of a tidal two-body system is attained through tidal recession, assuming the reduced Hill sphere is not reached. However, synchronisation can be achieved also via tidal approach, provided the Roche limit is not crossed. For each of the two scenarios, we derive the condition under which the evolving synchronicity radius catches up with the tidally evolving orbit. We consider the two scenarios for the Pluto-Charon system and examine the impact-origin hypothesis of Charon's formation against capture. Based on geophysical evidence, we propose that capture appears more likely. Motivated by this conclusion, we investigate the capture scenario, wherein the orbital evolution of Charon starts at a higher altitude than present and undergoes tidal descent, both analytically and numerically. We also consider the possibility that Pluto's initial prograde spin underwent a reversal by a tidally approaching retrograde Charon. Depending on the initial conditions, we observe temporary locking of Charon into higher spin-orbit resonances (3:2 to 7:2) during the first 0.5 Myr of the system's evolution. Owing to a greater initial separation between the partners, the power dissipated in each of them turns out to be much lower than in the case of tidal recession of bodies of the same internal structure. The greater initial separation also results in lower tidal stress, which may explain the absence of tidally generated fracture patterns.

Gravitational waves (GWs) are unique messengers as they travel through the Universe without alteration except for gravitational lensing. Their long wavelengths make them susceptible to diffraction by cosmic structures, providing an unprecedented opportunity to map dark matter substructures. Identifying lensed events requires the analysis of thousands to millions of simulated events to reach high statistical significances. This is computationally prohibitive with standard GW parameter estimation methods. We exploit DINGO-lensing, a deep-learning algorithm that accelerates the inference from CPU days to minutes to thoroughly reanalyze GW231123, the most promising lensing candidate to date. By performing more than 200,000 simulations with 3 different waveform models, we find that its statistical significance is below 4$\sigma$ and the event cannot be claimed as lensed. We observe that 8% of GW231123-like nonlensed simulations favor lensing, which could be explained by the self-similarity of short-duration signals. Still, 58% of GW231123-like lensed simulations have larger support for lensing, showing that higher detection statistics are possible. We show that analyzing simulations with different waveform models only lowers the significance, highlighting the relevance of waveform systematics. Although GW231123 exposes the challenges of claiming the first GW lensing detection, our deep-learning methods have demonstrated to be powerful enough to enable the upcoming discovery of lensed GWs.

We develop an exact framework for neutrino decoherence in power-law correlated turbulent matter, as encountered in core-collapse supernovae. Employing the Nakajima–Zwanzig projection technique, we derive an exact non-Markovian master equation for the neutrino density matrix. For kernels \( K(t) \propto t^{-\nu} \), the red-noise sector in our convention corresponds to \( \nu < 0 \), while \( \nu=1 \) is the white-noise boundary. To treat ultraviolet singularities for \( \nu \geq 1 \) without spoiling the fractional structure, we use a renormalization prescription based on Hadamard finite parts and analytic continuation. The exact Laplace-space solution for the survival probability is obtained. In the high-density matter basis relevant to supernovae, the solution is expressed through Mittag-Leffler functions, establishing a direct link to anomalous diffusion phenomena. For red spectra (\( \nu < 0 \)), the memory integral corresponds to a higher-order fractional operator. Our work clarifies how spectral index, renormalization scale, and decoherence efficiency interrelate, providing a complete analytical description and practical tools for supernova neutrino simulations. The fractional calculus formulation reveals fundamental mathematical connections between neutrino flavor evolution and other systems governed by long-range temporal correlations.

Accurate modelling of redshift-space distortions (RSD) is essential for maximizing the cosmological information extracted from large galaxy redshift surveys. In preparation for the forthcoming analysis of the Euclid spectroscopic data, we investigate three approaches to modelling RSD effects on the power spectrum multipoles of mock H$\alpha$ emission line galaxies. We focus on two one-loop perturbation theory models – the effective field theory (EFT) and velocity difference generator (${\rm VDG_ \infty}$) – which differ in their treatment of the real-to-redshift space mapping on small scales, and a third approach, the BACCO emulator, which adopts a hybrid strategy combining perturbation theory with high-resolution N-body simulations. We assess the ability of these models to recover key cosmological parameters, including the expansion rate $h$, the cold dark matter density parameter $\omega_{\rm c}$, and the scalar amplitude $A_{\rm s}$, across four redshift bins spanning $0.9 \leq z \leq 1.8$. In each bin, we find that ${\rm VDG_ \infty}$ and BACCO outperform the EFT model across all scales up to $k_{max} \lesssim 0.35 h\,Mpc^{-1} $. While BACCO saturates in constraining power at intermediate scales and higher redshift, the ${\rm VDG_ \infty}$ model continues to improve parameter constraints beyond $k_{max} \gtrsim 0.30 h\,Mpc^{-1}$. The EFT model, although robust on large scales, exhibits significant parameter biases for $k_{max} \gtrsim 0.25 h\,Mpc^{-1}$, limiting its applicability to Euclid-like H$\alpha$ samples. Among the full perturbation theory-based models, the enhanced treatment of small-scale RSD effects in ${\rm VDG_ \infty}$ improves cosmological parameter constraints by up to a factor of two.

We present the statistical validation and full photometric characterization of the exoplanet candidate TOI 7475.01 (TIC 376866659), detected by the TESS mission. Using a custom pipeline combining natural flux preservation with robust BLS detection, we identified a transit signal with a period of $3.2538$ days and a depth of $\sim 4600$ ppm. To rule out false positives, we performed centroid analysis, spatial contamination checks using Gaia DR3, and a statistical validation using TRICERATOPS. Our results show a Signal-to-Noise Ratio (SNR) of 294.13 and a False Positive Probability (FPP) of $\approx 0$. Based on the clean spatial environment, stable centroids, and high statistical probability, we validate TOI 7475.01 as a planetary companion. We subsequently performed a Bayesian transit fit using juliet with Dynesty nested sampling, and derived physical parameters via Monte Carlo error propagation. We find a planetary radius of $R_p = 1.18^{+0.39}_{-0.40}\ R_\mathrm{Jup}$ and an equilibrium temperature of $T_\mathrm{eq} = 1455^{+77}_{-56}$ K, consistent with a hot Jupiter classification. The planet mass is estimated at $M_p \approx 2.2\ M_\mathrm{Jup}$ (MAP) via the Chen&Kipping mass-radius relation; radial velocity follow-up is required for a precise measurement. The impact parameter remains unconstrained ($b = 0.46 \pm 0.34$), a limitation of single-band TESS photometry that future CHEOPS observations could resolve.

The Euclid galaxy survey will provide unprecedented constraints on cosmology, but achieving unbiased results will require an optimal characterisation and mitigation of systematic effects. Among these, the intrinsic alignments (IA) of galaxies are one of the dominant contaminants of the weak lensing (WL) and galaxy-galaxy lensing (GGL) probes. In this work, we assess IA modelling choices for Euclid DR1 3x2pt analyses by comparing the performance of the two most commonly used IA models, nonlinear alignment (NLA) and tidal alignment tidal torquing (TATT), along with several variations. Our analyses combine three perspectives: i) the constraining power on the IA and cosmological parameters for each IA model, ii) the bias that results when the IA analysis model differs from the model used to generate the synthetic data vector, and iii) the degeneracies between IA and photometric redshift (photo-z) nuisance parameters. Among the IA models analysed, the redshift-dependent TATT model (zTATT) provides the most flexible description of IA, with a similar constraining power compared to simpler IA models, making it a well-motivated choice for Euclid DR1 3x2pt analyses.

We present a multitask surrogate for neutron-star equations of state (EoSs) that delivers \emph{distribution-free}, certified uncertainty via split conformal prediction (CP) and its Mondrian variant. The surrogate ingests a six-parameter piecewise-polytropic representation $(\log_{10}p_1,\Gamma_1,\Gamma_2,\Gamma_3,\rho_1,\rho_2)$ – with fixed transition densities $\rho_1$ and $\rho_2$ – and jointly performs (i) validity classification under physical/observational constraints and (ii) regression of $M_{\max}$, $R(M_{\max})$, $R_{1.4}$, and $\Lambda_{1.4}$. Trained on a balanced set of $40{,}000$ EoSs, the model attains near-perfect discrimination (AUC $\approx 0.997$) and sub-percent relative errors for masses and radii, with few-percent error for tidal deformability. Across $\alpha\in[0.05,0.25]$, empirical coverages closely track $1-\alpha$ for both Standard and Mondrian CP; in conservative regimes, Mondrian yields narrower average physical widths at comparable coverage. To our knowledge, this is the first application of class-conditioned (Mondrian) conformal calibration to neutron-star EoS surrogates, enabling efficient, reproducible, and uncertainty-aware inference; the framework is readily extensible to functional targets (e.g., full $R(M)$ curves).

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

We propose a new, general form for a pseudo-Newtonian gravitational potential (PNP), expressed as a series of Paczyński-Wiita-like functions with the addition of increasing negative powers of $r$ with arbitrary coefficients. We present a procedure for determining these coefficients to construct a custom PNP that replicates key features of Schwarzschild geodesics for a test particle near a black hole. As an example, we construct potentials set to reproduce (I) the presence of an innermost stable circular orbit at the $r=6$ (geometric units), with the correct infall velocity for small deviations (on the geodesic universal infall), (II) the periapsis advance at large distances, and (III) the presence of a marginally bound circular orbit with specific angular momentum $L=4$, and the periapsis advance of parabolic orbits close to it. We compare the performance of our examples against the Paczyński-Wiita potential and other existing potentials. Finally, we discuss the limitations and advantages of our formulation.

Spatial curvature is one of the fundamental cosmological parameters that is routinely constrained from observations. The forward modelling of observations, in particular of large-scale structure, often relies on large cosmological simulations. While the so-called separate universe approach allows one to account for the effect of curvature on the expansion rate in small sub-volumes, the non-Euclidean geometry is harder to accommodate. It becomes important when observables are computed over large distances, e.g. when photons travel to us from high redshift. Here we present a fully relativistic framework to run cosmological simulations for curved spatial geometry. The issue of consistent boundary conditions is solved by embedding a spherical cap of the curved spacetime into a hole within a flat exterior, where it can undergo free expansion. The geometric nature of gravity is made explicit in our framework, allowing for a consistent forward modelling of observables inside the curved patch. Our methodology would also work with any Newtonian code to a good approximation, requiring changes only to the initial conditions and post-processing.

We investigate the propagating modes of New General Relativity (NGR) in second-order linear perturbations in the Lagrangian density (first-order in field equations). The Dirac-Bergmann analysis has revealed a violation of local Lorentz invariance in NGR. We review the recent status of NGR, considering the results of its Dirac-Bergmann analysis. We then reconsider the vierbein perturbation framework and identify the origin of each perturbation field in the vierbein field components. This identification is mandatory for adequately fixing gauges while guaranteeing consistency with the invariance ensured by the Dirac-Bergmann analysis. We find that the spatially flat gauge is adequate for analyzing a theory with the violation of local Lorentz invariance. Based on the established vierbein perturbative framework, introducing a real scalar field as matter, we perform a second-order perturbative analysis of NGR with respect to tensor, scalar, pseudo-scalar, and vector and pseudo-vector modes. We reveal the possible propagating modes of each type of NGR. In particular, we find that Type 3 has stable five propagating modes, \textit{i.e.}, tensor, scalar, and vector modes, compared to five non-linear degrees of freedom, which results in its Dirac-Bergmann analysis; the linear perturbation theory of Type 3 is preferable for applications to cosmology. Finally, we discuss our results in comparison to previous related work and conclude this study.

We report on a novel set of very-high-frequency quasi-periodic oscillations (VHFQPO's) in the context of compact, non-singular horizonless objects. Focussing on the static, spherically symmetric case we utilize metrics of non-singular black holes that are accompanied by a regulator length scale $L > 0$. The choice $L \gtrsim GM$ generically removes the horizon from these metrics leading to compact, horizonless but non-singular objects. This generically guarantees the existence of a stable orbit at small radii $r \ll r_\text{ISCO}$, independent of the angular momentum of the massive particle. Crucially, the absence of a horizon allows the resulting VHFQPO's to escape to infinity, spanning the range from 1kHz ($M = 10M_\odot$) to 25 kHz ($M = 2M_\odot$). Within the paradigm of non-singular spacetime geometries, the absence of such VHFQPO's from X-ray binary spectra implies the presence of a horizon around the central, compact object.

In this work we consider the scattering between non-relativistic particles with different finite sizes. We first calculate their interaction potential and apply the partial wave method to obtain their scattering cross section. Our findings show that the particle size can significantly affect the scattering between non-relativistic particles. Then we apply such a study to direct detection of puffy dark matter. We find that the finite size of the target nucleus may introduce non-perturbative effects that differ from the scenario of point-like dark matter. For large-size dark matter particles, this non-perturbative regime in the dark matter nucleus scattering cross section effectively disappears; while for small values of the size-to-range ratio in the scattering process, a significant non-perturbative regime can maintain. Finally, for the direct detection of nugget-type puffy dark matter with a small number of constituent particles, we find that the stability conditions for the formation of bound-state dark matter can provide constraints on the dark matter nucleus scattering cross section.

A gravitational-wave background (GWB) arising from the superposition of numerous unresolved gravitational-wave signals has yet to be detected. Potential contributing sources to such a background include compact binary coalescences (CBCs) and core-collapse supernovae (CCSNe). In this work, we place upper limits on the gravitational-wave energy emitted by CCSNe using cross-correlation measurements made with Advanced LIGO and Advanced Virgo detectors during their third observing run. Specifically, we obtain a $95\%$ credibility upper limit of $0.01~ {M_\odot c^2}$ while accounting for the contribution from CBC sources to a GWB. This result improves on previous constraint obtained from initial LIGO data by approximately two orders of magnitude. We also explore the detection prospects of third-generation ground-based detectors such as the Einstein Telescope and Cosmic Explorer for both individual CCSNe events and the GWB. Our results show that single events are likely to be detected prior to the GWB.

The recent observations of the Dark Energy Spectroscopic Instrument (DESI) indicated the possibility that the dark energy equation of state parameter $w$ might change from $w<-1$ to $w>-1$ when the redshift $z\sim 0.5$, which is called the inverse phantom crossing. In this paper, we investigate the possibility of the phantom crossing, and we construct realistic models realizing the crossing in the framework of the scalar–Einstein–Gauss-Bonnet gravity and ghost-free $f(\mathcal{G})$ gravity. We also investigate the scenario of the apparent phantom crossing, where dark matter energy density decreases more slowly than usually expected, which might explain the DESI observations. In the scenarios developed, the energy conditions are not violated by any component of the cosmic fluid. In the framework of the apparent phantom crossing, we also propose a new scenario, where the particle corresponding to the scalar field in the scalar–Einstein–Gauss-Bonnet gravity is dark matter. The mass of the particle might increase due to the coupling with the Gauss-Bonnet invariant, which makes the decrease of the dark matter energy density slower. This last scenario may suggest that the inverse phantom crossing might be related to the transition from the decelerating expansion of the Universe to the accelerating expansion.

We consider the gravitational Noether-Ward identities for the evolution of general metric perturbations on quantum matter backgrounds. In this work we consider Einstein's gravity covariantly coupled to a massive, non-minimally coupled, quantum scalar field in general curved backgrounds. We find that each term in the equation of motion for gravitational perturbations satisfies its own Noether-Ward identity. Even though each term is non-transverse, the whole equation of motion maintains transversality. In particular, each counterterm needed to renormalize the graviton self-energy satisfies its own Noether identity, and we derive the explicit form for each. Finally, in order to understand how the Noether-Ward identities are affected by the definition of the metric perturbation, we consider two inequivalent definitions of metric perturbations and derive the Noether-Ward identities for both definitions. This implies that there are Noether-Ward identities for every definition of the metric perturbation.

The theory of $f(R)$ gravity with constant curvature (i.e. constant scalar curvature) admits rotating and charged black hole solutions obtained from the Kerr-Newman-(A)dS metrics of general relativity through appropriate rescalings of the metric parameters. In this paper, we focus on the Kerr-Newman-de Sitter case and present a unified analytic treatment of the horizon structure and its physical properties, allowing for a transparent comparison between general relativity and $ f(R)$ gravity with constant curvature. We solve the quartic equation determining the horizon locations and derive closed analytic expressions for the horizon radii. Focusing on extremal configurations, we obtain analytic formulas for the squared rotation parameter $ a^2 $ and the inverse square of the curvature radius $ l^{-2} $ as functions of the horizon location and the electric charge. For generic values of these parameters, the extremality conditions are non-universal, reducing to the familiar Kerr-Newman bound only in the limit of vanishing background curvature. We identify an ultra-extremal configuration in which $ a^2 $ attains its maximal value at zero charge and decreases monotonically to zero as the charge approaches its limiting value, while $ l^{-2 }$ increases correspondingly. As an illustrative example, we show that black holes with charge $ q=M/2 $ necessarily possess a minimum rotation, which emerges naturally as an intersection point in our analytic description of $ a^2 $ and $ l^{-2 }$, when embedded in a universe characterized by a critical value of $ l^{-2} $ (equivalently, the scalar curvature or the cosmological constant). Finally, we demonstrate that when the mass satisfies $ M^2= (a^2+q^2)(1-a^2/l^2)$, the quartic horizon equation factorizes, leading in the extremal regime to a chiral-like horizon structure that allows only the outer-cosmological horizon merger.

We investigate the impact of gravitational wave (GW) dephasing due to gas accretion on the subtraction of massive black hole (MBH) binary signals over 4 yr of LISA data in the context of the global-fit. Based on state of the art predictions for the population of merging MBHs, we show that imperfect subtraction with vacuum waveform templates leaves a GW residual with an SNR of $3.2^{+5.4}_{-1.9}\times \sqrt{f_{\rm Edd} \langle \dot n \rangle/(20\, {\rm yr}^{-1})}$, where $f_{\rm Edd}$ is the typical Eddington ratio and $\langle \dot n \rangle$ the mean merger rate of LISA MBH binaries. We characterize the dependence of the residual on key population hyper-parameters and provide a simple fitting function. {Finally, we consider the distinguishability of the residual as a stochastic signal by comparing with Bayesian power-law sensitivity curves, while discussing several additional detection and mitigation strategies. Overall, our analysis indicates that the residual is unlikely to be confidently distinguishable from instrumental noise while nevertheless being likely to bias the inference of other signals.

The hyperon puzzle refers to the challenge of reconciling the existence of hyperons in neutron star cores and the observed high masses of neutron stars. The recent discovery of PSR J0952-0607 ($2.35\pm0.17 M_{\odot}$) has intensified this challenge. Existing solutions fail to achieve such a high mass, and often predict unrealistically fast cooling that is at odds with observations. Here, we propose a novel solution to the hyperon puzzle. Using the Dyson-Schwinger equation approach, we incorporate the quantum many-body effects caused by strong baryon-meson interactions into the equation of state for cold baryonic matter and find it stiff enough to support a maximum hyperon-star mass of $M_{\mathrm{max}} \approx 2.59 M_{\odot}$, which can explain all the observed high neutron-star masses. The resulting proton and hyperon fractions are remarkably low, thus the nucleonic and hyperonic direct Urca processes are significantly suppressed. As a result, fast cooling typically does not occur in ordinary neutron stars.

Through 2D3V PIC simulations of freely decaying sub-ion turbulence, intermittent localized regions with $\mathbf{E} \cdot \mathbf{B} \neq 0$ are found to be statistically associated with reductions in the magnitude of magnetic helicity while evolving in the early electron-scale interaction phase. Motivated by this behavior, we propose a source-compensated, history-dependent helicity density that satisfies an exact local balance identity by construction, enabling Saffman-type two-point correlation integrals which, under standard flux-decorrelation assumptions, can exhibit intermediate-scale plateaus that are roughly time-independent. In our simulations we demonstrate such plateaus to remain approximately invariant even as the usual Saffman helicity integral plateau value $I_H$ evolves during the early kinetic stage. Under approximate single-scale self-similarity, the plateau behavior of the magnetic integral is consistent with the 2D decay constraint $BL \sim \text{const}$. For initially net-helical configurations, we observe rapid development of mixed-signed magnetic helicity patches and a decrease of the global fractional helicity, such that the decay over the kinetic interval is again most consistent with the cancellation-dominated scaling constraint.

Black holes are some of the most interesting objects in the universe. While they first arise in the complicated behavior of general relativity, the physical laws ruling their behavior are surprisingly simple. For example, one of the core facts about black holes is that their area never decreases, much like the entropy in thermodynamics. In this note directed at introductory physics students and their instructors, we use this similarity to understand properties of black hole physics using standard techniques from an undergraduate course in thermal physics. We explore the never-decreasing nature of black hole area to obtain bounds on the energy emitted in a black hole merger (a calculation originally done by Hawking). We show how this allows us to think of black holes in manners very similar to heat engines, and how these ideas have been used in modern gravitational wave observatories to test general relativity. This allows a research-level topic to be discussed in introductory physics lectures.

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