Papers

We present a computational audit of the identifiability limits for structural corrections to the Radial Acceleration Relation (RAR) using a canonical subset of the SPARC database (N = 126). Rather than proposing a new dynamical law, we establish the observational conditions under which such a law would be mathematically recoverable. We resolve three methodological controls that heavily influence RAR interpretations. First, we isolate a -0.39 dex residual offset between the eight galaxies in the lowest-quality observational tier (Q = 3, mean -0.433 dex) and the remaining 118 (mean -0.040 dex). Of five predictions that beam smearing must satisfy, four fail: the radial profile plateaus at -0.331 +/- 0.028 dex rather than decaying to zero, and the offset does not scale with the number of resolution elements across a curve. Beam smearing is disfavoured as the primary driver, though the offset remains a data-quality signature rather than a physical one – Q = 3 in SPARC flags major asymmetries and strong non-circular motions, conditions under which a rotation curve does not trace the equilibrium potential. Second, we decouple the architectural limits of the dataset into three independently measured quantities – the no-model point scatter (sigma_M0 = 0.1860 dex), the residual floor after free per-galaxy intercepts (sigma_M3 = 0.1058 dex), and the propagated analytic error floor – together with the absorbable budget sqrt(sigma_M0^2 - sigma_M3^2) = 0.1530 dex derived from the first two. Finally, we provide an 8-cell protocol grid to reconcile Leave-One-Out (LOO) Mean Squared Prediction Error (MSPE) ratios. We demonstrate that reported structural-dynamical couplings must be evaluated with strict adherence to residual definitions (median vs. mean, signed vs. absolute) and baseline denominators to avoid adopting labeling artifacts as new physics.

We present a simple model with a n + 4-dimensional metric with compactified extra dimensions with a space dependent radius R. We show that changes in the profile of R are sourced by (visible) matter density and Newtonian potential. We can reproduce rotation curves of spiral galaxies and provide an expression for R. In order to test our hypothesis, we propose six experiments.

The Einstein Telescope (ET) is a proposed underground infrastructure in Europe to host future generations of gravitational-wave (GW) detectors. Two configurations are under investigation. The first is the two-site solution with an L-shape detector at each site similar to the LIGO configuration. The second configuration is the single-site triangular detector. While numerous analyses over the last decade have already proven that the distributed network of GW detectors is the scientifically preferred solution, it was also argued that the advantage might not be decisive enough to dismiss ET's triangle configuration. The modeling of aspects such as noise correlations in these configurations has only been treated superficially in the past. In this article, we introduce a time-domain simulation of environmental noises in the ET triangle using available data from ET's candidate sites. We simulate and analyze correlated environmental noise as well as coincident glitches produced by transients in the environmental fields. We find that these features pose an important limitation to the observational capabilities of the ET triangle and limit the utility of its null stream.

Recent detections of Galactic diffuse neutrinos by IceCube and $\gamma$-rays by LHAASO offer direct probes into the origin of Galactic cosmic rays. Conventional diffuse templates typically assume a single cosmic-ray injection spectrum across a wide energy range, without accounting for independent contributions from distinct accelerator populations. Here, we present a numerical framework that models Galactic diffuse neutrino and $\gamma$-ray emission that incorporates contributions from both microquasar and supernova remnant populations. By anchoring CR injection to local observations and utilizing high-resolution 3D target gas distributions, our model suggests multi-population contributions to the diffuse sky: escaped cosmic rays from supernova remnants dominate below $\sim10\text{ TeV}$, while those from microquasars become the primary driver at higher energies. Our predicted neutrino flux agrees well with the recent 12-year IceCube measurements, establishing a physically motivated baseline for the diffuse hadronic background while leaving room for unresolved point-like sources. This flexible framework provides testable predictions for current and future multi-messenger observatories, accommodating diverse accelerator populations and updated observational constraints.

High-resolution spacecraft images provide important astrometric constraints for orbit refinement, but measurements of resolved bodies are often limited by labor-intensive control-point selection and the difficulty of achieving consistent reductions over large image archives. We present an automated shape-model-based astrometric pipeline for Phobos and apply it to Mars Express Super Resolution Channel (SRC) images. For each exposure, a synthetic image is rendered from a high-resolution 3D shape model under the nominal spacecraft-target-Sun geometry. Feature correspondences between the observed and synthetic images are established using SuperPoint and SuperGlue, followed by RANSAC filtering. The matched synthetic-image keypoints are then associated with surface points through ray-shape intersection. The geometric adjustment fixes the adopted body orientation, spacecraft state, and corrected camera pointing and estimates only two effective plane-of-sky position offsets using the exact perspective-projection model. These offsets are used to derive the center-of-figure position of Phobos. We first test the method on an image set previously analysed with a control-point approach and obtain comparable astrometric performance. We then extend the analysis to a larger SRC dataset spanning 2007-2025 and obtain 1113 successful measurements. Relative to the JPL MAR099 ephemeris, the resulting observed-minus-computed residuals have mean values of 0.186 km in $\alpha \times cos(\delta)$ and 0.053 km in $\delta$, with corresponding standard deviations of 0.609 km and 0.583 km. These results demonstrate that the proposed pipeline provides a practical approach to large-scale, homogeneous astrometric reduction of archival spacecraft images of Phobos, with potential application to other resolved bodies.

The objective of this study was to characterize the human-made object that impacted the Moon on 5 August 2026, to estimate the dimensions of the crater it produced, and to compare the estimates with the crater subsequently measured by the Lunar Reconnaissance Orbiter (LRO). The impactor was the discarded upper stage of a SpaceX Falcon 9 rocket, which struck the surface near Einstein crater at a velocity of 2.43 km/s. The kinetic energy and linear momentum of the impactor were calculated, and the crater diameter was estimated using pi-group scaling relationships and measured rocket-body impacts as empirical analogs. The kinetic energy of the impactor was found to be 1.18 x 10^10 J, equivalent to approximately 2.8 tonnes of TNT. Compact-body pi-group scaling yielded a final rim-to-rim diameter of 46-57 m, and the empirical analogs yielded 26-27 m, leading to a pre-observation preferred estimate of 25 m within a range of 20-30 m. Between 11 and 12 August 2026, LRO imaged a crater approximately 18 m in diameter and less than 3 m deep. The results revealed that the compact-body scaling overpredicted the observed diameter by a factor of 2.6-3.2, whereas the empirical-analog method overpredicted it by 39-50%. When the nearest-mass analog was applied as a single crater rather than as a double crater, the estimated diameter was approximately 18 m, in close agreement with the observation. It was found that the hollow and elongated geometry of the stage, rather than a single equivalent bulk density, is the principal factor governing the crater size. This study provides a direct comparison of impact-scaling methods against a measured artificial lunar crater.

Superconducting cosmic strings (SCSs) are characterized by three parameters: the dimensionless string tension $G\mu$, the current amplitude $Y$, and the vector coupling, $\tilde{e}$ as recently explored in \cite{Rybak:2024our}. Using the Charge-Velocity-dependent One-Scale (CVOS) model, we compute the stochastic gravitational wave background (SGWB) emitted by current-carrying loops, including the suppression effect of vector radiation. While previous studies have shown that large $\tilde{e}$ suppresses the SGWB, allowing larger $G\mu$ to evade pulsar timing array (PTA) bounds \cite{Rybak:2024our}, we demonstrate that this parameter space is definitively constrained by an independent, astrophysically clean probe: the Dark Ages global $21$-cm signal. We recalculated the recently derived $21$-cm bounds \cite{Si:2025vsj} into the generic $(\tilde{e},Y,G\mu)$ parameter space and find that at the threshold where vector radiation emission starts dominating, any coupling $\tilde{e}\gtrsim 10^{-7}$ and $Y\simeq 0.67$ and $G\mu\gtrsim 10^{-13}$ injects sufficient ionizing radiation to completely erase the $21$-cm absorption signal at redshift $z\simeq89$. Greater current amplitudes impose more stringent constraints. This excludes a large parameter space considered viable. Our model-independent results highlight the complementarity between future lunar-based $21$-cm experiments and gravitational wave (GW) observatories.

Extreme mass-ratio inspirals (EMRIs) can magnify small departures from Kerr dynamics into appreciable gravitational-wave phase shifts accumulated over many orbital cycles. We exploit this sensitivity to investigate the imprint of a rotating black hole's acceleration on an EMRI waveform. The spinning C metric poses two obstacles to the standard Kerr flux framework: the spacetime is not asymptotically flat, and the acceleration breaks the reflection symmetry that supports exactly equatorial circular timelike orbits. For sufficiently small acceleration $AM$, we therefore formulate the calculation in an intermediate Kerr-like wave zone satisfying $M/r\ll1$ and $Ar\ll1$, and construct a near-equatorial circular orbit by examining its coupled radial–polar stability. We derive the separated point-particle source for the spin$-2$ radial Teukolsky equation, construct a regular normalized angular solution, solve the radial equation using the Sasaki–Nakamura transformation and the Green function method, and couple the resulting horizon and far-zone fluxes to the adiabatic evolution of stable near-equatorial circular orbits. The framework recovers the Kerr limit and reproduces the dominant $l=2$ Kerr fluxes with relative errors of order $10^{-7}$. Acceleration modifies both radiation reaction and the orbital frequency, producing a characteristic nonmonotonic accumulated dephasing. For $M=10^6M_\odot$, $m_s/M=10^{-5}$, $a/M=0.7$, and $AM=3\times10^{-7}$, the dominant-mode dephasing slightly exceeds $1$ rad over one year. Thus even weak acceleration can generate an order-radian secular phase imprint on long-duration EMRIs within the controlled regime of the present approximation.

Estimating spectral density matrices for future space-based gravitational-wave detectors such as LISA is challenging due to the long duration of the data and the correlated instrumental noise across multiple time-delay interferometry channels. In this work, we investigate a specialized mean-field stochastic gradient variational Bayes (SGVB) procedure for fast posterior approximation in long-duration multivariate spectral density estimation. Based on the existing eigenbasis representation of the blocked Whittle likelihood approach, the posterior model applies a Cholesky factorization to represent the inverse spectral density matrix and models the resulting frequency-dependent entries with cosine basis functions, with a discounted regularized horseshoe prior assigned to the basis coefficients. We test mean-field SGVB as a stand-alone posterior approximation by comparing it with Hamiltonian Monte Carlo (HMC) targeting the same posterior model. In simulation studies based on autoregressive and moving average processes, SGVB produces posterior median spectral estimates close to those obtained from HMC at substantially lower computational cost. We then apply the method to two one-year, fixed-delay, stationary, noise-only LISA simulations and show that the SGVB spectral density estimates are consistent with HMC and Welch estimates across the LISA analysis band, while requiring substantially less computation. These results demonstrate that SGVB provides a scalable Bayesian approach to spectral density estimation for long-duration multivariate LISA noise analysis.

We study the thermal history of the primordial Universe in the presence of non-zero lepton and baryon asymmetries. Considering different scenarios, we determine the equation of state (EoS) of the Universe from T = 10 GeV down to T = 1 keV, spanning the QCD transition, hadron gas phase and neutrino decoupling epochs. Using a combination of numerical codes, we track the cosmic trajectories of chemical potentials associated with the baryonic, leptonic and electric charges, and follow the evolution of lepton asymmetries including through the era where neutrino oscillations take place. Combining peak theory with numerical-relativity simulations of the collapse threshold, we show the EoS-induced modifications to the primordial black hole (PBH) mass spectrum. We determine the associated Gravitational Wave (GW) signal, showing how lepton asymmetries and a particular spectral index of curvature perturbations can be hinted at by current ground interferometer-based GW observations. Finally, we discuss constraints and positive evidence for PBHs.

We have used archival GALEX data to separate the cosmic ultraviolet background at the Galactic Poles into two components: the dust scattered light and an offset. We have modeled the dust-scattered light using a single- scattering model finding 1 sigma limits of 0.54 – 0.71 for the albedo (a) and 0.74 – 0.83 for the phase function asymmetry factor (g) at 1530 Å and 0.66 – 0.73 for a and 0.71 – 0.77 for g at 2360 Å, that is, the grains are moderately reflective and highly forward-scattering. The offsets are 277 – 284 photon units at 1530 Å and 513 – 520 photon units at 2360 Å. We have estimated other Galactic and extragalactic contributors to the offset finding that 161 +- 18 photon units is unaccounted for at 1530 Å and 335 +- 38 at 2360 Å. The offsets are constant over these regions with variances of 20 – 30 photon units.

The electromagnetic properties of a dilute unpolarized nonrelativistic neutron gas are described in the case when neutron-by-neutron scattering is negligible. The Gaussian (Maxwell-Boltzmann) and Fermi-Dirac one-particle density matrices for the neutron gas are considered, the typical space scale of variations of these density matrices being assumed to be much larger than the the typical space scale of variations of the electromagnetic field. The Green functions for the static effective Maxwell equations are obtained. The parameters of a nonequilibrium neutron gas are found where this gas possesses a ferromagnetic instability. The analog of Friedel oscillations is described. The properties of plasmon-polaritons in the neutron gas are revealed. It turns out that the dispersion law of longitudinal and transverse plasmon-polaritons has an infinite number of branches at a given momentum. The region of parameters where the plasmon-polaritons can be regarded as quasiparticles is found. The plasmon-polaritons on a Gaussian wave packet of a single electron are described. Their dispersion law proves to have an infinite number of branches at a given momentum. It is shown that there exist stable longitudinal plasmon-polaritons even on a single electron.

We investigate a proof-of-concept code designed to model the functionality of current generation astronomical fiber positioner instruments. This code accommodates an arbitrary number of fiber probes and astronomical targets. To establish the prerequisites for future fiber positioning technologies, we conduct a systematic analysis of both the observational and instrumentation parameter space (such as target density, fiber density, and exclusion and patrol radius). The tool assesses the efficiency and completeness of such a system by modeling the allocation of fibers to targets over multiple visits or observations, enabling the quantification of performance based on the specific parameters. This tool holds significance in delineating the technical specifications for future multi-object spectroscopy instrumentation. This will inform the development of new technologies for the next generation of telescopes using novel instrumentation techniques or the extension of existing technologies.

The gravitational wave background detected by pulsar timing arrays is sourced by a finite population of supermassive black hole binaries, and is therefore anisotropic. We ask what measuring that anisotropy can teach us about the population, using models that span a wide range of effective source counts, all normalized to the measured background amplitude. We find four things. First, the expected anisotropy is produced by the single brightest binary: a dipole at the level of the published NANOGrav 95% upper limit would require one source to supply about 60% of the power in the band. Second, no model that also reproduces the measured strain spectrum contains a source that bright. In every case the loudest binary stays below the NANOGrav continuous-wave upper limit at every frequency it covers, consistent with the joint search of the 15-year data, which finds no resolved source. Third, because the anisotropy is produced by one source, compressing the sky to an angular power spectrum discards the phase information that locates it. Such a search is never more sensitive than looking for the source directly, and is strictly worse once more than a dipole is kept. Fourth, the published upper limits on the angular power spectrum therefore reflect the analysis prior rather than the data: they coincide with the 95th percentile of the prior induced by the square-root spherical harmonic basis adopted in the analysis. What does constrain the population today is the shape of the strain spectrum. Rare, bright sources depress the median spectrum below its mean, and fitting the measured spectrum already disfavors mass functions dominated by binaries above about 10^10 solar masses.

The core-cusp problem is a discrepancy between simulations and observations of the slopes of inner dark matter density profiles. Splitting dwarf spheroidals by metallicity is a standard method to break the mass-anisotropy degeneracy. Four dynamical methods are applied to the 1339-star Tolstoy et al. (2023) VLT/FLAMES Sculptor catalog: spherical Jeans, GravSphere, the Walker and Peñarrubia (2011) two-population estimator, and a continuous $f(\boldsymbol{J},[\text{Fe}/\text{H}])$ action-based distribution function. This work finds $\Gamma=2.49^{+0.23}_{-0.20}$, with 99.6% of the posterior lying above $\Gamma=2$ (84.4% after $\lambda$ calibration), robust to membership cuts and metal-poor contamination; GravSphere gives $\gamma=0.46^{+0.40}_{-0.28}$. Controlled mocks that hold tracer geometry fixed while breaking the metallicity-kinematics link attribute $\sim77$–$83\%$ of the recovered bias to geometry rather than decomposition itself. The continuous distribution function encodes a metallicity–action coupling $k_J$ inaccessible to any split-based method; a direct action-space computation on 575 stars with Gaia proper motions finds $[\text{Fe}/\text{H}]$ declining by $\sim0.3\,\text{dex}$ across three decades in radial action. On the central-density diagnostic of Read et al. (2019), the same posteriors give $\rho_\text{DM}(150\,\mathrm{pc})=1.77\times10^8\,M_\odot\,\mathrm{kpc^{-3}}$, agreeing with their independent Sculptor measurement at $0.74\sigma$ and placing Sculptor in their cusp-like class. The shallow $\gamma$ of this work is an asymptotic quantity; the local slope at $150\,\mathrm{pc}$ is consistent with the Navarro–Frenk–White profile under the spherical models used here, though a recent axisymmetric analysis using resolved 3D velocities finds a shallow slope over that range and identifies an inclination degeneracy to which spherical models are blind.

In 1980, the Universe was made of stars, cosmology was the province of less than 100 astronomers, the Hubble constant was only known to within a factor of two, and the hot big bang was the “the standard model of cosmology.” Particle physics was a thriving enterprise concerned with the inner space of quarks, leptons and the strong, weak and electromagnetic interactions. The Standard Model was newly established and the aspirations were for a “grand unification” of the forces and particles. By 2000, the agendas of inner space and outer space were inextricably linked, and $\Lambda$CDM with its inflation, particle dark matter and dark energy – and a Hubble constant reliably known to better than 10% – was the new standard model. The big questions connecting inner space and outer space included the origin of ordinary matter, the identity of dark matter, the nature of dark energy, and an understanding of the origin of space, time and the Universe. Particle dark matter, the idea that the bulk of the mysterious dark matter is a new, long-lived or stable elementary particle, played a central in bringing together inner space and outer space, and today is still a hot topic in both fields.

[CII] emission is a powerful tool for studying the gas content in galaxies, especially at early stages of evolution. Given its low excitation potential, the [CII] 158 \mu m line is the main coolant of neutral gas in photodissociation regions and giant molecular clouds, hence a tracer of all phases of the interstellar medium. However, there is no consensus on what is the nature of the gas traced by this emission. We aim to provide insights into the physics of the [CII] emission and predictions for future observations. We used four pre-prepared simulations performed with a version of p-gadget-3, which includes the krome chemistry package. We implemented a postprocessing semi-analytical model to simulate the [CII] emission. We analysed an isolated Milky Way-mass-size galaxy, and three merger configurations: a co-rotating, counter-rotating, and perpendicular major mergers of two Milky Way-mass-size galaxies. We reproduce fundamental relations such as the [CII] luminosity and the star formation rate, L[CII]-SFR relation, in which mergers show an extra component for the [CII] emission not traced by the star formation alone. During quiescent or starburst stages, galaxies deviate from the L[CII]-SFR relation. We find that [CII] is a robust tracer of molecular gas, as well as atomic gas in both cold and warm phases on a global scale. Our results indicate that around 50 per cent of the total [CII] emission is traced by molecular gas, while 30 to 40 per cent is traced by atomic and ionised hydrogen regardless of the configuration. Our findings suggest that the nature of the [CII] emission depends strongly on the evolutionary history of galaxies, in which mergers and starbursts act as more efficient drivers of [CII] emission through collision with molecular hydrogen.

Binding energies on amorphous solid water regulate molecular desorption in molecular clouds and protoplanetary disks, but computing them from first principles across the heterogeneous binding sites of disordered ice surfaces is computationally prohibitive. We aim to develop a rapid, interpretable method for predicting binding energies on amorphous solid water that captures surface heterogeneity without requiring explicit electronic-structure calculations. We trained a machine-learning model using 27 geometric descriptors of the local ice-adsorbate environment. The model was trained on binding energies for 13 adsorbates on amorphous solid water clusters and evaluated through pooled leave-one-cluster-out validation across 15 ice surfaces. Its transferability was assessed on radical species excluded entirely from training. The model achieves a pooled leave-one-cluster-out coefficient of determination $R^2 = 0.90$ (the fraction of variance in the binding energies captured by the model, where 1 is a perfect fit), and a mean absolute error of 378 K across 15 ice surfaces. It also transfers to radical species withheld from training, with a pooled $R^2 = 0.79$, although its accuracy decreases when the unpaired electron contributes directly to the surface interaction. Our results demonstrate a route towards rapid, distribution-aware parameterisation of binding energies and desorption/diffusion rates in gas-grain astrochemical models, while preserving physical interpretability on heterogeneous ice surfaces.

Whether the slow solar wind originates from one coronal source or two distinct channels remains a central open question in heliophysics. Resolving this requires unsupervised separation of two populations that arrive at nearly the same bulk speed and differ mainly in heavy-ion composition. We present Solar-CDC, a self-supervised contrastive deep clustering (CDC) framework that maps plasma observables to a latent space via a Transformer encoder, optimizes a triplet margin loss, and updates pseudo-labels via $k$-means. Theoretically, we prove that neighborhood-preserving embeddings such as t-SNE and UMAP are fundamentally constrained. Preserving the neighbor graph leaves the cross-cluster cut fraction unchanged, and preserving all but a fraction $\varepsilon$ of its links moves that fraction by at most $\varepsilon$. Neither bound depends on the target dimension. A margin objective rewrites the graph and drives the cut fraction to zero. Empirically, on 30,602 Solar Orbiter observations, thirty combinations of dimensionality reduction and clustering peak at a silhouette of $0.454$, whereas Solar-CDC reaches $0.869$. Escaping the geometric bound alone does not guarantee physical validity: TriMap also optimizes triplets and reaches $0.824$, yet its clusters score below chance against the published composition taxonomy. Solar-CDC instead recovers clusters with mean charge-state ratios of $0.080$, $0.160$, and $0.400$, placing the intermediate population inside the window associated with coronal-hole boundaries. Even when the defining charge-state ratio is withheld from the inputs entirely, the model still recovers the taxonomy defined on it. Solar-CDC thus connects self-supervised representation learning to coronal source diagnostics. Importantly, a learning loss recovers physical populations only when driven by dynamically updated physically-aware clusters rather than distances.

To leading order, the gravitational-wave (GW) signal from binary neutron stars depends on the individual spins, $\chi_1$ and $\chi_2$, only through the effective spin parameter $\chi_{\rm eff}$. We present the first systematic investigation of individual-spin effects on ejecta, $r$-process nucleosynthesis, and kilonova emission, including comparisons at fixed total mass, mass ratio and $\chi_{\rm eff}$. We use numerical relativity ejecta from three total mass regimes with the finite-temperature, composition-dependent SFHo equation of state and neutrino emission and absorption. For $M_{\rm tot}=2.55\,M_\odot$ at fixed $\chi_{\rm eff}=0$, individual spins change the dynamical ejecta mass by a factor of ${\sim}45$, while the absolute $A\geq140$ yield spans more than two orders of magnitude and the lanthanide to light $r$-process mass ratio increases from ${\sim}2$ to ${\sim}70$. Prompt-collapse $4.10\,M_\odot$ models show heavy-element yield differences exceeding four orders of magnitude at $\chi_{\rm eff}=0$. The kilonova retains the individual-spin imprint, with peak brightness differences reaching ${\sim}0.9$ mag. At $40\,{\rm Mpc}$, all three fixed $\chi_{\rm eff}=0$ pairs remain above adopted depths at common epochs for all viewing angles, with same-epoch colour differences reaching ${\sim}1.5$ mag. The colour imprint persists when the simulation-derived secular ejecta are replaced by the same parametric disc outflow, indicating that disc mass differences are not the primary driver and that dynamical ejecta make an important contribution. Neutrino absorption systematically brightens the kilonova and shifts peak-associated colours blueward. These EM signatures can break the degeneracy between individual spins in the GW signal.

Stellar mass black holes accreting at a super-Eddington rate produce strong outflows and non-thermal radiation up to PeV energy. Recently, very high energy (VHE) emission was discovered from the environment of the galactic microquasar V4641 Sgr, which is known to experience super-Eddington outbursts. We present MHD simulations supplemented with a nonlinear kinetic Monte Carlo model that reproduces the observed radio, X-ray, and gamma-ray emission from the region over tens of parsecs surrounding V4641 Sgr. The reported multi-wavelength non-thermal emission can be explained as produced by PeV electrons accelerated by a system of strong shock waves created by the interaction between a narrow trans-relativistic polar jet and a quasi-spherical accretion-disk wind. The shock positions are associated with bright spots of the gamma-ray and X-ray radiation, whose spatial profiles can be tested with dedicated observations. The model predicts that a sizable fraction of the jet's kinetic energy is channeled into TeV-PeV protons, suggesting that such systems might be powerful sources of PeV-range cosmic rays during the accelerator activity cycle. In this model, the bow-tie shape structure discovered by MeerKAT in V4641 Sgr is associated with the radio emission from the termination surface of the quasi-spherical accretion-disk wind. The newly presented (from the SRG/eROSITA all-sky survey) and archival (INTEGRAL and Swift BAT) X-ray data are quantitatively consistent with the model, and simultaneously reproduce VHE measurements with H.E.S.S., HAWC, and LHAASO. Deeper X-ray observations will provide a test of the intermittent operation of the accelerator, in contrast with more persistent sources such as SS433/W50.

Investigating the correlation of flux variability in blazars across multiple spectral bands provides a key tool for constraining emission mechanisms and jet physical parameters. We study the gamma-ray and optical variability correlation of BL Lacertae on timescales of about 1 day during September 2022, when the source was within the TESS field of view, and over a longer period (April 2022 - September 2023) encompassing its gamma-ray flare, which finishes the series of flares during the high-activity phase. Our analysis confirms the Synchrotron Self-Compton mechanism as the origin of the gamma-ray emission, characterized by a zero time lag. For the longer interval, the Discrete Correlation Function (DCF) reveals a second statistically significant peak at the time lag of about 12 days. We demonstrate that this peak arises due to the presence of an additional flare in one of the bands, which accompanies a gamma-ray and optical flares correlated at zero time lag. We interpret the orphan flare origin in terms of the propagation of an injected blob of high-energy electrons along a jet with a decaying magnetic field, accounting for electron radiative cooling.

Coherent radio emission is observed from compact astrophysical objects with relatively weak magnetic fields, including millisecond pulsars and white dwarfs, the latter being proposed as possible sources of long-period radio transients. In such environments, the standard pair discharge mechanism—driven by curvature radiation and one-photon pair production—can fail because the low magnetic field strength suppresses photon conversion. We analyze a discharge mechanism that operates efficiently in weak-field magnetospheres, including two-vertex quantum electrodynamic processes: inverse Compton up-scattering of background photons followed by either two-photon or one-photon pair creation depending on the seed photon temperature. Using first-principles radiative particle-in-cell simulations incorporating exact QED cross sections, we demonstrate that these mechanisms robustly generate pair cascades in millisecond-pulsar magnetospheres and in hot white-dwarf environments. The resulting system exhibits limit-cycle behavior and generates electromagnetic field fluctuations capable of producing coherent radio emission, providing a natural explanation for radio activity in low-field compact objects.

We extend the finite-solid-angle general relativistic radiation transport method of C. J. White et al. (2023) to time-dependent spacetimes represented in ADM form. This Valencia-type solver retains the angular transport and local implicit matter coupling of the original HARM-type solver, but replaces its time-independent Kerr–Schild tetrad and conserved-variable normalization with an Eulerian formulation. A Cholesky-gauge spatial tetrad makes the frame and its derivatives algebraic functions of the ADM variables, providing smooth, metric-compatible angular transport. The resulting transport system couples directly to an evolving spacetime, agnostic to the gauge evolution equations. Our test suite, including flat and curved beams, radiation–fluid coupling, and time-dependent geometries, establishes the accuracy and robustness of this approach. We further apply the Valencia-type solver to a radiative circumbinary disk, illustrating its potential for multi-messenger studies of dynamical strongly gravitating systems.

To investigate the initial conditions of high-mass star-forming regions, we use SiO (2-1) emission to trace early shock-related kinematics toward sixteen 70 $\mu$m-dark and massive clumps with luminosity-to-mass ratios ($L/M$) $< 1\,L_{\odot}/M_{\odot}$, as part of the Multiwavelength Line-Imaging Survey of the 70 $\mu$m-dark and bright clouds (MIAO) project. Using ALMA observations at a spatial resolution of $\sim$0.06 pc and a velocity resolution of 0.21 km s$^{-1}$, we identify a total of thirty-seven outflows with a variety of morphologies. Outflow parameters were derived by integrating the HCO$^+$ (1-0) line wings, excluding the quiescent dense core component traced by H$^{13}$CO$^+$ (1-0). We find that outflow masses and velocities show moderate positive correlations with the masses of their driving cores. Owing to the high sensitivity of our observations, which yield longer projected outflow lengths compared to previous studies, the derived outflow dynamical ages span $\sim10^{3}$-$10^{5}$ yr. We detect six narrow-linewidth (0.6-1.4 km s$^{-1}$) and three broad ($>$ 2 km s$^{-1}$) SiO (2-1) features not associated with outflows driven by clearly identified protostars. Lacking coincident 3 mm dust continuum cores, their origins may be young outflows from undetected low-mass protostars, dissipating shocks, cloud-cloud collisions, or projection effects when the outflows lie close to the plane of the sky. The detection of these shocks and outflows in such extremely young environments demonstrates that protostellar activity has already begun.

In this work we explore the imprints of self-interacting dark matter (SIDM) on the structure of simulated galaxy groups and clusters in view of future optical surveys. We find that both baryonic and projection effects lead to a significant alleviation of the SIDM core-forming behavior on the mass and density profiles, resulting in very small ($\lesssim 5\%$) deviations from standard collisionless dark matter. The opposite trend emerges for our most massive halos, where SIDM tends to significantly enhance the central density instead; we compare those halos with the observed dark matter profiles of two galaxy clusters, MACS J1206 and Abell S1063, finding that SIDM may leave potentially detectable imprints on cluster cores when considering the small observational uncertainty of current kinematic mass measurements. Finally, we find that weak lensing observations provide a promising tool for testing SIDM on cluster scales: in an idealized scenario where systematics are under control, weak lensing shear profiles exhibit up to $\sim 20\%$ level deviations for $M_{200} \gtrsim 10^{13}\ \text{M}_\odot$ halos. Conversely, the projected concentration-mass relation is only slightly modified by our SIDM models. We also discuss the limitiations of our analysis and the impact of the simulation resolution, as well as possible extension of this work. Overall, our results provide a baseline for future SIDM studies: although line-of-sight projection and baryonic effects significantly reduce SIDM imprints on halo mass profiles, weak lensing and kinematic analyses in galaxy clusters may potentially allow for robust detection and constraining of dark matter self-interactions in massive structures from megaparsec down to kiloparsec scales.

We have discovered PSR J2108+5055, an extremely luminous but sporadic mode-changing rotating radio transient (RRAT), with a period of 0.495 s using the Canadian Hydrogen Intensity Mapping Experiment/Fast Radio Burst instrument (CHIME/FRB). We show that the brightest pulse from PSR J2108+5055 has a flux density of 3.1+/-1.0 kJy, corresponding to a spectral luminosity of (7.5+/-2.3)x10^23 erg s^-1 Hz^-1, assuming a dispersion-measure-derived distance of 1.6+/-0.3 kpc. We place an upper limit on the duty cycle of PSR J2108+5055 at 0.002% and measure a surface magnetic field of 5.8x10^11 G, a characteristic age of 1.2x10^7 yr, and a spindown luminosity of 2.2x10^32 erg s^-1. Using a simple model based on the Green Bank North Celestial Cap (GBNCC) survey, we place an upper limit on the number of PSR J2108+5055-like sources with similar peak flux densities in the GBNCC survey area of ~80 RRATs. We show that if PSR J2108+5055 were placed at the edge of the Local Group of galaxies, the Five-Hundred-Meter Aperture Spherical Telescope (FAST) would detect its brightest pulse. The second brightest pulse has a peak flux density of 1.5+/-0.7 kJy, and FAST would detect such a pulse out to ~1 Mpc, covering the inner Local Group, including Andromeda.

Recent DESI DR2 observations have renewed interest in extensions of the $\Lambda$CDM cosmological model, particularly through indications of a time-varying dark energy equation of state. In this work, we investigate whether such deviations may also involve interactions within the dark sector. We consider an interacting dark energy scenario in which the dark matter density evolves as $\rho_{\rm dm}\propto a^{-3+\delta}$, with the constant $\delta$ quantifying the interaction strength, and allow the dark energy equation of state to be either constant but different from $-1$, or dynamically evolving through the CPL parametrization. The models are constrained using Planck CMB data, DESI DR2 BAO measurements, and three Type Ia supernova compilations: PantheonPlus, Union3, and DES-Dovekie. For the constant equation-of-state case, the inclusion of DESI and supernova data leads to a preference for a small negative interaction parameter, with a significance above $2\sigma$. When dynamical dark energy is allowed, the evidence for interaction becomes weak, while the data favor a quintessence-like evolving dark energy component. In both scenarios, Bayesian model comparison still favors $\Lambda$CDM. Our results show that the inferred role of dark-sector interactions depends strongly on the nature of dark energy, highlighting the importance of jointly testing dark energy dynamics and interactions in the DESI era.

Context. R Coronae Borealis (RCB) stars are hydrogen-deficient and carbon-rich (HdC) supergiants that likely originate from a merger of two white dwarfs. As such they are not expected to be members of close binary systems, in agreement with current observations. Aims. We test the above hypothesis and search for eclipses of RCB stars that could be due to a companion. Methods. We exploited old archival photometric data from Harvard photographic plates (DASCH) as well as modern-era photometric surveys such as Gaia, ASAS-SN, and ZTF to search for eclipses, and conducted follow-up spectroscopic observations. Results. We found one RCB star (WISE J005128.08+645651.7) that features two 3.5 yr long dimming events. We argue that these may be eclipses separated by an orbital period of 66 yr and that the star is an eclipsing binary candidate. The next eclipse is predicted for 2080. If confirmed, this would be the first RCB star in a binary system, the second-longest orbital period of any known eclipsing binary, and the longest duration of an eclipse. Light curve simulations support the idea that the eclipses could be caused by a dusty disk around the companion. The star shows unusually large pulsations. Our follow-up spectroscopy indicates a stellar temperature of about 11000 K, the presence of sharp emission lines (including forbidden lines) during the normal brightness, and extremely strong diffuse interstellar bands. The star is likely surrounded by a hotter chromosphere and an extended low-density nebula. We calculated atmosphere models of the star. Models feature a density inversion and the derived chemical composition is compatible with some cool extreme He stars (EHe). The star constitutes a bridge between the classical RCB, hot RCB, and EHe stars. Conclusions. Future research is encouraged to confirm or refute the binary nature of the star and search for other RCB binaries.

Radio pulsars can power their coherent radio emission through intermittent plasma discharges in magnetospheric gaps. The nonlinear coupling between particle acceleration, quantum electrodynamic (QED) processes, and electric field screening remains difficult to model self-consistently. We present an analytical concurrency model and first-principles particle-in-cell (PIC) simulations of polar cap discharges. The model predicts limit-cycle behavior where the electric field and plasma density oscillate with a frequency set by the local plasma frequency, modified by the plasma inertia. One-dimensional simulations with exact QED rates and realistic plasma parameters validate these predictions. The discharges generate pair multiplicities up to $M_\pm \sim 10^4$ and excite electric field oscillations with spectra consistent with pulsar radio microstructure. These results provide a theoretical framework connecting microphysical gap dynamics to macroscopic observables.