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

Binary population synthesis (BPS) codes are extremely useful tools for investigating both the end-to-end lives of binary stars as well as a myriad of astrophysical phenomena observed in the Universe. Given the widespread use of BPS in modern astrophysical research, it is imperative to systematically compare BPS codes across the spectrum of computational efficiency, flexibility, and physical realism to gauge their consistency and robustness. In this work, we perform BPS using three modern codes---the rapid code $\texttt{COSMIC}$, the hybrid code $\texttt{METISSE}$ integrated into $\texttt{COSMIC}$, and the detailed code $\texttt{POSYDON}$---on three single-metallicity populations of identical initial binaries, ensuring consistent choices in physical parameterizations where possible. Investigating the final population of merging binary black holes (BBHs) as a test case, we find stark differences in the properties, formation pathways, and progenitors across the three codes. In an initial population of one million binary stars at $0.01 Z_\odot$, each code results in $\sim 5,000-9,000$ BBHs that merge within a Hubble time. However, only one initial binary becomes a BBH merger in all three codes, and $\lesssim 14 \%$ of BBH progenitors consistently merge in two codes. Binaries that become BBH mergers in two codes often go through different evolutionary pathways and result in different final properties. In short, the codes are inconsistent in predicting BBH merger properties, even for identical initial binary systems. Our results highlight the need for systematic comparisons of BPS techniques, for a deeper understanding of physical and computational differences between BPS codes, and for caution in over-interpreting the results from any BPS code. (Abridged)

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.

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.

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.

Learned surrogates for expensive inner solver blocks are a widely pursued route to faster multiphysics simulation. We report a controlled, end-to-end negative result and identify three structural barriers, none of them a deficiency of the network we trained. The testbed is the implicit Newton solve coupling energy-dependent neutrino radiation to matter in a general-relativistic radiation-hydrodynamics code, its most expensive physics routine per call. First, per-call cost and share of runtime are different quantities, and only the second bounds acceleration. An exclusive self-time profile puts the target block at 16.9% of critical-rank wall clock, capping any surrogate at ~1.2x by Amdahl's law. A surrogate 5.8x cheaper per call merely ties the solver, and the configuration stable enough to run without fallback reaches only parity. Second, offline accuracy cannot rank surrogates for deployment: across fourteen networks the pooled Spearman error-versus-survival correlation (rho = +0.73) is a between-family confound that vanishes under control (rho = -0.04). Third, a correct out-of-distribution gate cannot accelerate a loop that leaves its training distribution. We give the break-even deferral fraction in closed form: because the visited states sit 73x off the data manifold, the gate defers 96.8 to 99.7% of cells, almost invariant to surrogate quality. Including its own cost, the gated loop is a 0.94 to 0.96x slowdown. We further separate stability from fidelity: a never-crashing gated run accumulates a linear -19.9% density bias over 6000 steps. The error is a directed, ballistically accumulating bias, not the variance-driven divergence the autoregressive literature targets.

The large-scale dipole anisotropy observed by the Pierre Auger Observatory above $8$~EeV provides important clues about the origin of ultra-high-energy cosmic rays (UHECRs). In this work, we investigate whether gamma-ray active galactic nuclei (AGNs) from the Fermi-LAT Fourth AGN Catalog (4LAC) can explain the observed dipole and examine the contribution of nearby starburst galaxies (SBGs). Incorporating source fluxes, redshift distributions, GZK attenuation, and a mixed-composition framework, we demonstrate that single-source populations of AGNs and SBGs alone fail to reproduce the observed Auger dipole. To resolve this, we construct a hybrid AGN+SBG model. The best-fit solution for the $8$--$16$~EeV interval requires a $23$\% AGN contribution and $45$\% SBG contribution, which gives a dipole amplitude of $5.03$\% while maintaining the dipole direction of ($104.1^\circ$, $-25.1^\circ$) with angular separation of $7.10^\circ$. In the other energy intervals also, we find that the nearby SBGs increasingly dominate the anisotropic component while AGNs provide a subdominant contribution alongside a quasi-isotropic background.

Electron-capture supernovae (ECSNe) are a promising formation channel for low-mass neutron stars, but the minimum gravitational mass of the neutron stars they produce depends on both the progenitor core structure and the neutron-star matter equation of state (EOS). We compute the electron-capture (EC) threshold gravitational mass ($M^\star_{WD}$) of cold white-dwarf-like O--Ne--Mg cores with representative compositions and map the baryon number onto cold neutron-star configurations constructed from a Bayesian ensemble of unified crust--core EOSs. Although the EC threshold density is sensitive to the concentrations of the O--Ne--Mg mixture, the threshold baryon number of the O--Ne--Mg core varies only weakly, producing a narrow remnant-mass window. In the baseline case with no baryonic mass loss during the transition from the EC threshold mass O--Ne--Mg core to the remnant neutron star, standard ECSNe yield remnant neutron stars with gravitational masses of $1.24$--$1.265\, M_\odot$, with only a small EOS-induced spread. Small baryonic mass losses of $0.01$--$0.02\, M_\odot$ shift this range modestly downward, but the $1.174\, M_\odot$ companion of PSR J0453+1559 would require an extreme mass loss close to $0.10\, M_\odot$, which is not favored by current ECSN simulations. We find that the residual EOS dependence of the remnant mass is controlled mainly by the pressure around nuclear saturation density, while the corresponding tidal deformability remains sensitive to the remnant radius and compactness. Thus, low-mass double neutron star systems can in principle connect ECSN-like formation channels with gravitational-wave constraints on the EOS. Our results show that ECSNe naturally form low-mass neutron stars, but within a restricted mass range; the lightest observed neutron stars likely require low-mass iron-core collapse, ultra-stripped supernovae, or other nonstandard channels.

We report results from a search for an isotropic stochastic gravitational-wave background using data collected by the LIGO--Virgo--KAGRA Collaboration. The analysis uses data from the first observing run through April 1, 2025, during the fourth observing run. New frequency-domain cuts are implemented to address a class of non-stationary spectral noise features that were not effectively identified and mitigated by existing data-quality checks in past analyses. Consequently, previously analyzed data from the fourth observing run are re-processed with the updated cuts. We find no evidence for a stochastic background signal and place upper limits on the gravitational-wave energy density. In particular, for a background following a power law with spectral index 2/3 as predicted by inspiralling compact binaries, we find $Ω_\mathrm{GW}(25\,\mathrm{Hz}) \leq 2.0 \times 10^{-9}$, while scale-invariant backgrounds are constrained to $Ω_\mathrm{GW}(25\,\mathrm{Hz}) \leq 2.8 \times 10^{-9}$, both at the 95\% credible level for a log-uniform prior on $Ω_\mathrm{GW}$. Relative to the constraints from previous data recomputed with the new frequency-domain cuts, these limits improve by a factor of 1.4. We also update bounds on alternative gravity scenarios predicting non-standard polarization modes, and we verify that correlated magnetic noise sources remain below the sensitivity of this search. Combining these observational constraints with population models of compact binary coalescences informed by the latest gravitational-wave transient catalog, GWTC-5.0, we predict the amplitude of the compact binary background to be $Ω_\mathrm{CBC}(25\,\mathrm{Hz}) = 6.3^{+5.0}_{-2.2} \times 10^{-10}$ at the 90\% credible level.

We investigate how cosmological environment regulates the evolution of primordial black hole (PBH) seeds in the early universe using a suite of hydrodynamical simulations. Unlike traditional seeding channels, PBHs provide early BH seeds without the need of special conditions for gas collapse/fragmentation that favor overdense regions, thus covering diverse large-scale environments. We find that PBHs follow distinct evolutionary pathways depending on the gas supply and halo assembly history. In underdense regions, limited inflow suppresses both accretion and star formation, producing faint, metal-poor systems, while in overdense environments, sustained gas inflow drives rapid BH growth and the formation of compact, centrally concentrated stellar components. These differences lead to large variations in BH-to-stellar mass ratio, metallicity, and morphology, collectively bracketing a range of possible evolutionary pathways for PBH-seeded systems. We show that compact, BH-dominated sources resembling recently observed Little Red Dots naturally arise as one phase within these evolutionary pathways before evolving into more extended galaxy--AGN systems. Our results suggest that both the initial seed properties and cosmological environment jointly shape early BH--galaxy co-evolution, while subsequent environmental regulation can erase the memory of the initial seeding channel.

In this work, we present results from general-relativistic radiative-transfer calculations of black hole accretion and jet launching to reproduce the observed properties of the M\,87 jet at sub-parsec scales. For the first time, in the simulations we have included a blob component at the launching region of the jet that allowed us to reproduce the quasi-simultaneous low-energy spectral energy distribution ($10^{10} \text{ Hz} \leq ν\leq 10^{16} \text{ Hz}$). Moreover, we obtain a better fit of the synchrotron emission in the self-absorption region ($ν\leq 10^{11} \text{ Hz}$), showing that the inclusion of a blob component in our radiative transfer calculations results in a flatter spectrum ($α_{22-86\ \text{GHz}}=0.08$). Additionally, our simulations improved the accuracy of the observed morphology of M\,87's jet, matching the data up to a distance of 0.6 mas from the core at 86 GHz. Our results reproduce the edge-brightened structure up to $\sim1.3$ mas, i.e., in a more extended region compared with previous studies, with a brighter southern edge. The obtained synthetic image of the jet reproduces additional observed knots along the southern edge.

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.

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.

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.

Periodic density structures (PDS) are trains of plasma enhancements carried into the heliosphere, but their repeated timing does not by itself identify where the cadence is set. We test whether PDS selected independently in the outer corona can be traced back to a low-coronal clock, and whether that organization survives through the EUVI--COR1--COR2 observing chain. We analyze 12 STEREO-A/SECCHI events on 2008 January 11--14 along a nonradial path. EUVI 171 A at 1.10--1.20 R_sun shows an unusual concentration of power in the predeclared 80--130 minute band relative to filter-matched red-noise controls, while the ensemble is most strongly organized in upper COR1 at 2.5--3.0 R_sun. These signatures identify a low-coronal modulation candidate and an intermediate-height organization domain, but they do not form a unique phase-preserving clock extending into COR2. Instead, the more persistent observable is spatial order. Events 9 and 12 retain expansion-stable outward ordering. In polar r--theta maps, intersections of oppositely inclined, expansion-aware matched-filter ridge supports form repeated X-/diamond-like patterns. The morphology is compatible with stationary or quasi-stationary shock-cell processing, although the brightness diagnostics are insufficient to establish a stationary MHD shock branch. The observations therefore favor an intermittent source--gate--transfer picture: a low-coronal cadence may modulate plasma release, while event-dependent propagation progressively destroys exact phase coherence. The outward ordering of the density structures can survive, providing a more persistent signature of the source-to-wind transfer than phase locking itself.

We present phase-connected timing and polarization measurements for two long-period FAST-CRAFTS pulsars, PSRs J0000+6252 and J2131+3642, and extend single-pulse emission-state analysis to a five-source sample including PSRs J1903+1407, J1502+4653, and J2112+4058. FAST L-band timing baselines span 414 to 511 days; the two pulsars have spin periods of 1.11 to 1.55 s, period derivatives of $(3.99$ to $4.06)\times10^{-15}~{\rm s~s^{-1}}$, characteristic ages of 4.34 to 6.16 Myr, surface magnetic fields of $(2.15$ to $2.52)\times10^{12}$ G, and spin-down luminosities of $(4.21\times10^{31}$ to $1.16\times10^{32})~{\rm erg~s^{-1}}$. Their rotation measures are $70.2\pm26.7$ and $-44.6\pm7.6~{\rm rad~m^{-2}}$, with linear polarization fractions of 17.6\% to 27.2\%. PSR J2131+3642 shows a short monotonic PA segment permitting a formal rotating-vector-model fit, though limited longitude coverage leaves the geometric parameters poorly constrained; PSR J0000+6252 has too few PA points for such a fit. Gaussian mixture modeling (GMM) of single-pulse energy distributions identifies null, weak, and burst components in PSRs J0000+6252, J1903+1407, J1502+4653, and J2112+4058, while J2131+3642 shows only weak and burst states. We also identify bright single pulses (peak intensity $\geq10\times$ the integrated average profile): 42 in J0000+6252, four each in J1903+1407 and J2112+4058, and none in J2131+3642 or J1502+4653. These bright pulses occur within the main emission window with no evidence of periodic recurrence, consistent with sporadic enhancements of the normal radio-emission beam. For J2112+4058, we measure a scattering timescale $τ_{\rm sc}=5.84\pm0.18$ ms at $ν_{\rm ref}=1.25$ GHz. Together, these results highlight the diversity of magnetospheric variability among slowly rotating neutron stars.

The Solar wind is a continuous flow of charged particles emanating from the solar surface and governed by complex, interacting magnetohydrodynamic processes. Accurate specification of inner-boundary conditions is essential for heliospheric modeling and solar-wind prediction. In many practical applications, only a subset of interacting multi-field variables is directly available, but for a comprehensive view of solar wind prediction and downstream magnetohydrodynamic simulations, a more complete boundary state is required. In this work, we study the problem of learning the multi-field multi-scale solar magnetohydrodynamic state at 30 solar radii ($R_\odot$) using operator learning. Specifically, given the radial velocity and radial magnetic field, we aim to reconstruct the non-radial velocity and magnetic field components, radial and non-radial current density, thermodynamic density, and pressure components. This mapping is highly nonlinear, spatially coupled, and multi-scale, making it a challenging task for data-driven scientific machine learning. To address this problem, we employ a Local Neural Operator (LocalNO) that learns mappings between input and output function spaces while retaining locality and resolution-awareness. Unlike conventional regression models and autoencoder models, neural operators are better suited for learning structured field-to-field transformations arising from physical systems. The resulting predictions along with inputs are intended to serve as boundary condition variables for future inner-heliospheric modeling pipelines.

Cosmological first-order phase transitions (FOPTs) in the early Universe are an exciting prediction of many beyond the Standard Model scenarios. Owing to the stochastic nature of bubble nucleation, some causal patches may remain trapped in the false vacuum long after the surrounding regions have completed the transition and started to redshift. Under suitable conditions, these delayed patches can become sufficiently overdense to collapse into primordial black holes (PBHs). One therefore expects characteristic parameters, such as the PBH mass and PBH formation time to be sensitive to the underlying phase-transition dynamics. Here, we adopt a model-independent framework with a fixed mean bubble separation scale ($R_*$) to directly compare the impact of exponential and Gaussian nucleation profiles on PBH formation. We further incorporate the effects of friction on the bubble walls arising from interactions with the surrounding plasma. We find that friction significantly modifies the bubble dynamics, leading to significant changes in the PBH mass and formation time for both nucleation histories. Moreover, the stochastic gravitational-wave spectrum generated by the FOPT can provide a complementary probe of the underlying dynamics, allowing the effects of friction on the bubble wall evolution to be distinguished from the frictionless case.

We present the capabilities of the Fluorescence detector Array of Single-pixel Telescopes (FAST) observatory in the single telescope configuration as an important step towards the possible future large-field observatory for detecting ultra-high-energy cosmic rays. Reconstruction of main shower physics parameters are explored on noise-free simulated events using machine learning techniques in the challenging domain of a low-intensity transient signal. We find a very good correlation between the true and reconstructed energy of the shower even with the information from just the four photomultipliers of the single FAST telescope, and a reduced performance for the maximum of the shower development Xmax, using various architectures of artificial deep and convolutional neural networks, with a comparison to a benchmark gradient boost regression model. The resolution in the energy is found at the sub-percent level, while in Xmax it is~$5\%$. The relative difference between predicted and true values is under one percent for energy, while for Xmax it ranges from $-8\%$ to $+17\%$, which can be attributed to the limited information from the single FAST telescope configuration. The results constitute an important capabilities verification and a lesson learned with implications for established FAST prototypes as well as for more complex configurations of the FAST observatory under construction.

Background: Localized structure in the equilibrium squared sound speed can affect neutron-star observables, but modifying $dP/dε$ requires thermodynamically consistent reconstruction and assessment over the complete declared domain. Purpose: We quantify the response to positive and negative localized sound-speed changes around a unified BSk24 baseline while keeping the deformation geometry, admission criteria, reconstruction, and fixed-mass comparison protocol fixed. Method: A Gaussian profile with quintic smootherstep activation is added to $c_{\rm s}^2$. Proposals violating $P>0$ or $0<c_{\rm s}^2\leq1$ anywhere in the retained domain are rejected before reconstruction. Passing cases are reconstructed as effective cold one-fluid barotropes and evolved with the Tolman-Oppenheimer-Volkoff and tidal equations. Results: The $A=0$ case reproduces the undeformed control exactly. A localized change produces a persistent pressure offset, relocates the fixed-mass central state, and yields a mass-dependent, nonlinear, sign-asymmetric response. For $A=-0.8$ and $+0.8$, respectively, $ΔΛ/Λ_0=(-46.02\%,+35.33\%)$ at $1.4\,M_\odot$ and $(-63.52\%,+87.17\%)$ at $2.0\,M_\odot$. The reached stiffened feature remains core connected at the reporting endpoints, whereas the softened feature becomes an off-center shell at high mass. Conclusions: These phenomenological BSk24-anchored effective barotropes illustrate how a local sound-speed intervention propagates through thermodynamic reconstruction and stellar structure. They are not BSk24 predictions or evidence for a microscopic composition or phase transition.

Recent observations by the James Webb Space Telescope (JWST) have indicated the presence of water in the inner regions of the disc around the young protostar PDS 70, but the origin of this water remains unclear. Here we report the discovery of variable absorption lines of neutral sodium in archival High Accuracy Radial velocity Planet Searcher (HARPS) spectra of PDS 70. These lines vary strongly and stochastically on a daily basis, in amplitude, number, and radial velocity. Our measurements indicate that this gas tends to be optically thick, and partially cover the stellar disc, meaning this fast-moving gas is often spatially confined. We explore several hypotheses for the origin of these lines, and conclude that a likely source of the observed sodium is sublimation of planetesimals that transit the star on highly elliptical orbits, reminiscent of the exocomet phenomenon seen in other seen in other extrasolar systems. These sublimating exocomets could play a role in sourcing the previously observed water in the inner, terrestrial-planet forming region of the PDS 70 system.

While no signal originating at a black hole's event horizon can reach future null infinity, information about the horizon and its immediate vicinity can be encoded in asymptotic properties of waves emitted by matter or field perturbations falling toward a growing/forming horizon. Within a response-filtered framework, "direct wave" denotes source-sensitive plunge and remnant-formation information revealed by filtering the black-hole response. Unlike a stationary damped sinusoid with a fixed complex frequency, the direct wave follows the evolving source and has an evolving instantaneous complex frequency tied at late times to the remnant horizon angular velocity $Ω_H$ and surface gravity $κ_H$. Using rational filters, we remove quasinormal modes from numerical-relativity waveforms to study this evolution. Before numerical contamination, trajectories depend on remnant spin: the real frequency evolves toward $2Ω_H$ from above for lower spins ($χ_f\lesssim0.7$) and from below for higher spins ($χ_f\gtrsim0.7$), while the instantaneous decay rate increases toward $\sim2κ_H$ and, in some high-spin cases, beyond it. As in particle-plunge results, the horizon-controlled value is approached only at late times, as frame dragging controls near-horizon motion. For $χ_f\sim0.7$ remnants of non-precessing, comparable-mass binaries, the early-time real frequency is close to $2Ω_H$ because the binary orbital frequency transitions smoothly to the remnant horizon frequency. Finite-time deviations therefore carry information about merger/collapse dynamics rather than undermining the horizon connection. Full direct-wave evolution requires numerical-relativity calibration. We further show that approximate pole-zero pairing in the Kerr response motivates a spin-dependent minimal filter set that suppresses quasinormal-mode features while revealing source-trajectory information.

We present a detailed observational analysis of BLG723.31.34651, the fourth reported eclipsing member of the class of Roche lobe-filling hot subdwarfs, using time series spectroscopy and multi-band photometric observations from the Goodman High Throughput Spectrograph on the 4.1-m SOAR telescope. Our light curves show periodic photometric variability at $P=64~{\rm min}$ from the tidally-deformed hot subdwarf donor and deep eclipses of the donor by an opaque accretion disc. We report two light curve model solutions that agree with each other to within $2σ$, both suggesting that BLG723.31.34651 contains a typical hot subdwarf donor, a white dwarf accretor, and an irradiated opaque accretion disc. Our spectroscopic analysis identifies statistically significant phase-dependent variations in the effective temperature and surface gravity coinciding with the phase-dependent variations seen in our observed light curves. We attribute these variations to a combination of intrinsic variations in local atmospheric parameters due to tidal deformation of the donor and unmodeled spectroscopic flux contributions from an accretion disc.

Atmospheric photochemical models incorporating the impacts of stellar flares often assume a $\sim$9,000 K spectrum at ultraviolet-optical wavelengths. Recent multiwavelength observations, however, reveal a more complex picture with temperature measurements spanning 4,000-40,000 K, although the occurrence rates for flares with different temperatures remain unknown. Here, we model the evolution of a Proterozoic Earth-like world with 0.01 bar of O$_2$ under repeated flaring to identify the impact of flare effective temperatures. We explore four scenarios - two host star types (K2V and M2.5V) and two flare temperatures (9,000 K and 19,000 K) - selected to bound the potential parameter space. The hotter flares have a larger impact on O$_3$ photochemistry for both stellar types. M-star planetary atmospheres are more volatile and exhibit rapid changes in their O$_3$ production and destruction rates. Meanwhile, K-star planetary atmospheres are more stable and are only impacted by the hottest flares, proving advantageous for biosignature searches. We simulate 0.2-1.0 $μ$m reflected light spectra for all four scenarios, and find that 19,000 K flares can result in either production or destruction of O$_3$ depending on the host star spectral type increasing the 0.2 $μ$m feature by $\sim$2$\times$ for the K2V star but decreasing it by 50% for the M2.5V star. Future missions such as the EVE SMEX mission concept will provide robust flare temperature constraints for young FGKM stars, which will serve as inputs to improve photochemical models to inform future HWO observations.

We prove that black-hole isospectrality between linearly independent solutions follows whenever radiative metric and matter perturbations are reconstructed from complex master variables obeying a closed, complex-linear system with complex-linear boundary conditions. If a phase rotation maps one parity sector to the other, then the even- and odd-parity QNM spectra coincide, and the perturbations are parity isospectral. A broad class of chiral-aligned theories satisfies these criteria. As an application, we demonstrate that subextremal Kerr--de Sitter black holes are parity isospectral.

We investigate the dynamical instability and gravitational collapse of charged, spherically symmetric anisotropic matter configurations within $f(R,\mathcal{L}_{m})$ gravity. The analysis focuses on the combined effects of anisotropic pressure, perturbations, electric charge, and modified-gravity source terms on the stability of compact objects. A specific equation of state is adopted to relate the static and perturbed variables through the adiabatic index. Using a perturbation scheme, we derive the modified hydrostatic equilibrium and collapse equations and obtain instability constraints in both Newtonian and post-Newtonian regimes. The results show that the stability of the system is governed by the competition between inward gravitational attraction and outward pressure support. In addition, the dark source terms generated by the non-minimal matter-geometry coupling modify the collapse conditions and can enhance the stability of the charged fluid. These findings provide a useful framework for understanding the evolution and collapse of dense self-gravitating compact objects in non-minimally coupled gravity.

It is widely repeated in the popular literature and elsewhere that the singularity at the center of a black hole is a point. It is not true. Two observers who free-fall into a spherical black hole along two different angular trajectories at the same time $t$ do not encounter each other at the central singularity; rather, they lose causal contact with each other already well away from the singularity. Counterintuitively, in general relativity two points can be spatially close yet causally distant. The singularity is a surface, not a point. The story for rotating black holes is more complicated, but the same conclusion holds. For a rotating black hole, the singular surface almost certainly resides at its inner horizon, where even the tiniest classical or quantum perturbations ignite the exponential mass inflation instability, precipitating collapse to a spacelike singular surface. There are implications for quantum gravity. We argue that, whatever the ultimate theory of quantum gravity may be, the quantum states of a black hole probably reside at its effectively 2-dimensional singular surface, which coevolves unitarily with, and in thermodynamic equilibrium with, the hot atmosphere of trapped Hawking radiation that the black hole generates within its event horizon.

We investigate the inflationary dynamics of a spectator massive vector field with a non-minimal curvature coupling in de Sitter (dS) and quasi-de Sitter (qdS) backgrounds, focusing on the transverse and longitudinal sectors. By solving the mode equations with Bunch-Davies initial conditions, we compute the spectral energy densities and quantify the deviations between the exact dS approximation and the slow-roll-corrected qdS evolution. We find that the coupling parameter xi significantly affects the evolution of the effective mass, the mass-crossing condition, and the behavior of the longitudinal mode. In particular, increasing xi suppresses the vector fluctuations and reduces transient tachyonic effects. We further follow the evolution of the modes into the radiation-dominated era, showing how the inflationary initial conditions determine the subsequent spectral energy densities. Our results identify the parameter regime in which the dS approximation provides an accurate description and determine when slow-roll corrections become relevant.

Newtonian mechanics is usually written on a spacetime of fixed geometry. Here we build that geometry, introducing a propagator and a bilinear form only as the description demands. Within this reconstruction the Principle of Inertia emerges as a theorem rather than an independent postulate. Inertial motion is blind to the scale of the form and to the torsion the propagator defines, and that blindness is a geometric precursor of the \emph{universality of free fall}, before gravity enters. Newton's Second Law is a compatibility between two derivatives along a curve. The force exceeds the mass times the lowered acceleration by an exact term, and the law survives arbitrary torsion. Requiring force-free motion to be inertial along curves parametrised by arc length restricts the torsion to be totally antisymmetric. The remaining torsion is a three-form which we postulate the particle carries along its worldline. Its contraction with the velocity annihilates that velocity, so the supplementary condition of spinning-particle mechanics holds identically, and we propose that contraction as a classical spin. The algebra is familiar from relativistic spin hydrodynamics; the assignment, to the best of our knowledge, has not been proposed elsewhere. Since inertial motion resolves neither the mass nor the spin, the predictive content lies in the transport law we postulate, which turns the spin against a parallel-transported frame at a rate fixed by the acceleration alone and by no coupling constant. That rate becomes a laboratory number only once the arc-length parameter is identified with a clock, which we leave open. The dynamics asks non-degeneracy of the form, so this spacetime is not Newton--Cartan's Galilean one, and gravitation lies outside the construction. We prove every mathematical result, make three constitutive postulates and one physical conjecture.

Exceptional control of systematic effects is required in order to achieve unbiased cosmic shear two-point correlation function measurements with the next generation of galaxy imaging surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). One critical challenge is accurately modeling the point-spread function (PSF), as errors in PSF estimation can introduce spatially correlated biases in galaxy shape measurements. The LSST Science Pipelines, which will be used to process Rubin data, include an implementation of the PIFF (PSFs in the Full Field of View) package originally developed for, and demonstrated to perform well on, DES-Y3 and Y6 data. In this work we use semi-realistic image simulations, mimicking LSST observing conditions in $i$-band over 100 square degrees, to perform an end-to-end test of the Rubin PSF modeling pipeline. PSF model residuals are quantified using both second- and fourth-order moment parameters and performance is evaluated, for both LSST year 1 and year 10 depth, with a series of diagnostic tests. We find that the additive bias contribution from PSF modeling errors to the non-tomographic cosmic shear data vector is well below 30% of the cosmic shear uncertainty estimated from our analytic covariance matrix. Though our simulations exclude several known effects that may further challenge PSF modeling, these results demonstrate promising performance from PIFF on LSST-like data and provide an early benchmark for ongoing PSF validation efforts for LSST weak lensing analyses.

The constancy of Newton's gravitational constant G is a fundamental prediction underlying general relativity and a critical assumption in theories of gravity beyond it. We report a high-precision test of a possible temporal variation of G in a strong gravitational field, using high-resolution ultraviolet spectra of Ni V absorption lines in the white dwarf G191-B2B obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS). By comparing the measured line centroids with laboratory wavelengths and isolating the gravitational-redshift contribution, we obtain dG/G = (-0.014 +/- 0.016) x 10^-15 yr^-1 at a strong-field potential phi approximately 10^4, corresponding to an absorption redshift z_abs approximately 8.47 x 10^-5. The result is consistent with a constant G and places a stringent constraint on its secular variation in the strong-field regime. Our measurement demonstrates the potential of white-dwarf stars as a sensitive probe of temporal variations in the gravitational coupling.

A black hole need not ring down in isolation, since a nearby companion can subject it to an external gravitational tide. We calculate how a quadrupole tide modifies the Schwarzschild black hole resonances. Linearising the exact tidally distorted Schwarzschild solution in the tidal amplitude $A$, we derive the diagonal odd- and even-parity perturbation equations on the deformed background and compute the formal first order displacements of analytically continued Schwarzschild resonances. Each diagonal first order shift factorises as $δω_{n\ell m}=A\,c_{\ell m}κ_{n\ell}$, giving one reduced complex coefficient for each $(n,\ell)$, expressed as a ratio of contour integrals, multiplied by a universal Zeeman-like pattern that splits the $2\ell+1$ azimuthal multiplet. For the $\ell=2$ and $\ell=3$ sectors analysed explicitly, the odd- and even-parity shifts coincide, so Schwarzschild isospectrality survives at first order in the tidal amplitude. Quadratic tidal forces, by contrast, can break this degeneracy. Moreover, for a physical companion, the astrophysically dominant $(n,\ell,m)=(0,2,2)$ mode oscillates faster and decays more slowly. Finally, we show that the eikonal shifts admit a geometric description in terms of the Penrose limit about the tidally deformed photon ring.

We study how discrete self-similarity (DSS) in the critical behavior of scalar-field collapse is imprinted on primordial-black-hole (PBH) mass functions. Using the DSS-modulated critical scaling law found in cosmological simulations, we propagate the near-threshold mass map into normalized PBH mass functions during a {kination} era. We compare Gaussian window function, $k$-space top-hat window function, and real-space top-hat window function, together with a real-space top-hat window function multiplied by a kination transfer function. We find that critical behavior in gravitational collapse produces an irreducible minimum width even for an infinitesimally narrow primordial spectrum. DSS then modulates this critical-scaling profile, generating approximately log-periodic features in mass. At fixed horizon mass, successive equal-phase points of the DSS-modulated critical mass map satisfy $Δ\ln M_{\rm PBH}=γP_{\rm ln}$. Consequently, in the narrow-spectrum limit, the corresponding structures in the final mass function are expected to satisfy $Δ\ln m\simeqγP_{\rm ln}$, or $m_{n+1}/m_n\simeq5.6$, for the fiducial Choptuik DSS parameters $γ$ and $P_{\rm ln}$. For broad primordial spectra, the convolution over horizon masses dephases the DSS pattern and progressively washes out the critical substructure. We normalize the mass functions so that we may primarily study the profile shape and the survival of DSS substructure, rather than the absolute PBH abundance. These results provide a bridge between cosmological DSS collapse simulations and PBH phenomenology including population observables.

A large suite of N-body/hydrodynamical simulations is used to demonstrate that the physical properties of the merging cluster El Gordo are well reproduced within a two-state self-interacting dark matter scenario. Our findings show that, in addition to elastic scattering, the presence of an inelastic, endothermic, up-scattering channel is essential to satisfy observational constraints derived from the measured weak lensing profiles of the two colliding clusters. We find that the observational properties of El Gordo are adequately matched by our self-interacting dark matter model within the following interval of parameters: an elastic cross-section in the range of ~ 4 -5 cm^2/gr, and an up-scattering channel bounded by an inelastic cross-section of approximately ~ 2 -4 cm^2/gr with a threshold velocity of ~1200 - 1600 kms, respectively. We argue that our conclusions depend critically on current uncertainties in the primary mass and the projected separation d_{DM}. In fact, our findings suggest that if error estimates in weak lensing measurements were significantly reduced - with future observations showing evidence of a low primary mass (<~ 10^{15} M_sun ) and a larger projected separation d_{DM} >~ 740 kpc - it would then be difficult for the proposed self-interacting scenario to be reconciled with observational constraints. We further suggest that future lensing surveys of massive, high-velocity merging clusters could be decisive in testing the proposed dark matter scenario.

Computing boson star families traditionally requires repeated solution of nonlinear eigenvalue boundary-value problems and careful numerical continuation through turning points. We develop a physics-informed neural network (PINN) that learns the map from the physical parameters and radial coordinate directly to the scalar and metric fields over an equilibrium solution manifold. Regularity and asymptotic boundary conditions are incorporated into the network output, while the training objective combines pointwise supervision, Einstein-Klein-Gordon residuals, and curve-level constraints on the Arnowitt-Deser-Misner mass and Noether charge. A trained model generates a complete configuration in a single forward pass. Across representative one-, two-, and three-branch families, the method reconstructs the mass-frequency spirals and conserved quantities, including configurations on inner branches that require delicate continuation in conventional solvers. These results establish physics-informed surrogate learning as a practical route to amortized exploration of nonlinear self-gravitating solution families.

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.

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.

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 $μ$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 $μ$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.

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.

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 study the thermal history of the primordial Universe in 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.

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.

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.

Superconducting cosmic strings (SCSs) are characterized by three parameters: the dimensionless string tension $Gμ$, 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μ$ 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μ)$ 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μ\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.

The canonical scenario for the formation of subdwarf B (sdB) stars involves the ejection of the progenitor's envelope near the tip of the Red Giant Branch (RGB), concurrent with the onset of core He-burning. While binary interactions are known to dominate sdB formation, the origin of apparently single sdB stars remains uncertain. We investigate the conditions under which a sdB progenitor can eject its envelope through the engulfment of a substellar companion, using an energy balance approach. We further estimate the fraction of sdB stars that could form via this channel. We simulate the orbital evolution of substellar companions during the subgiant and RGB phases to determine the onset of engulfment and calculate the energy released during inspiral within the stellar envelope. By comparing this energy with the envelope binding energy and exploring different efficiencies for its deposition, we identify the conditions required for envelope ejection. We then combine these results with the observed population of substellar companions to estimate the occurrence rate of sdB formation. The engulfment of substellar objects can lead to envelope ejection within a limited region of the mass-semi-major axis parameter space, whose extent depend on the host star's properties and energy transfer efficiency. The minimum companion mass required for ejection increases with decreasing efficiency and increasing stellar mass. As a result, the frequency of envelope ejection is strongly influenced by the distribution of substellar companions, particularly by the paucity of objects in the brown-dwarf desert. Low-mass brown dwarfs and massive planets are viable candidates for explaining the existence of apparently single sdB stars. Our results define the region of parameter space where this mechanism is energetically possible and provide a guide for future multidimensional hydrodynamical simulations.

We comment on an error in the charged quark star literature ( Phys. Rev. D, 102(3):034031, 2020 ) that resulted in erroneous solutions to the radial stability equations for charged interacting quark stars in ( Phys. Rev. D, 104(12):123007, 2021 )[2]. In this comment article we recalculate the fundamental radial eigenfrequencies for charged, interacting quark stars in general relativity and discuss any departure from previous results [2]. Most of the qualitative trends persist, but the magnitudes of the separation between the maximum mass and stability points on a given branch can change significantly. We also note that all branches in [2] that contained no stable solutions now have stable solutions in their corrected counterparts. The net effect is that in all of the studied cases charged quark stars have a greater range of stability than originally predicted. We also compare four recent observational constraints not included in previous studies, and find that these constraints are satisfied by a subset of the mass-radius relations we obtain.

AU-scale structures in the interstellar medium have become relevant across several areas of astrophysics, from ISM microphysics to pulsar timing arrays and cosmic-ray (CR) propagation. Current sheets formed at magnetic-field reversals have drawn particular attention because they can deflect the dominant GeV CR population through large angles. Pulsar scintillation, an outstanding probe of plasma density, resolves these structures---but not their magnetic configuration. Conventional Faraday rotation measures (RM) are only marginally sensitive at AU scales and are further compromised by ionospheric systematics. We propose a method that isolates birefringence-induced polarization---the signature of magnetic reversal---by applying phase retrieval to polarized data. Pulsar scintillation further resolves ray paths separated by $\sim 1\,\rm{AU}$, enabling differential RM measurements that suppress the ionospheric common mode while preserving the AU-scale magnetic signal. We demonstrate the method on an archival observation of PSR B0834+06. Phase retrieval resolves the 1-ms scattering feature into two branches in Doppler-delay space. Under the corrugated current-sheet interpretation, these branches correspond to ray paths sampling opposite sides of an AU-scale current sheet. We measure a branch-to-branch RM difference of $(-9.3\pm3.2)\times 10^{-3}\text{ rad m}^{-2}$ at 2.9$σ$ across a 10-MHz band. The signal is consistent with a field reversal of $|Δ\langle B_\parallel \rangle| \simeq4.4\pm2.1\,μ\rm G$, sufficient to sustain a long-lived sheet. With only six independent subbands, we present this as a proof of concept rather than a definitive detection. If it proves scalable, its full implementation could constrain cosmic-ray transport and inform models of propagation-induced noise in pulsar timing arrays.

Recent observations have revealed gas-rich debris disks around intermediate-mass stars at ages of tens of Myr. The origin of this gas remains unclear: it may be primordial, retained from the protoplanetary phase, or secondary, released from volatile-rich solids. Secondary-origin models reproduce CO emission but often overpredict neutral carbon. Recent observations and disk-evolution models suggest that primordial gas may survive longer than previously assumed, motivating thermochemical tests of the primordial-remnant scenario. We test the previously unexplored possibility that primordial-origin disks satisfy the observational constraints on gas-rich debris disks. Specifically, we determine under what conditions a disk around a $2\,M_{\odot}$ star reproduces substantial CO, low CI/CO ratios, and weak HCO+ emission consistent with current non-detections. We post-processed 20-40 Myr structures from 1D disk-evolution models with Cloudy, varying irradiation geometry, dust-to-gas mass ratio (DTG), and cosmic-ray ionisation rate. In the dust-poor models (DTG $=10^{-4}$), CO remains optically thick around $R\sim100$~au. The models yield low disk-integrated CI/CO mass ratios. Our model produces CO radial intensities of the observed order of magnitude, but its CI-emitting region is more extended than observed. The standard CR model overproduces HCO+, whereas the weak CR model brings its predicted luminosity within current observational limits. These results demonstrate that a primordial origin remains chemically viable for CO-rich debris disks. The main remaining tension is the excessive radial extent of the CI emission, although it may reflect our simplified modelling. Further testing of the primordial-origin scenario will require multidimensional, self-consistent modelling, spatially resolved CI observations, and deeper searches for HCO+.

We provide a geometric foundation for the kinetic theory of photons. Although the induced metric $\hat h$ on the light cone bundle $Γ_0^+$ is degenerate, which causes the corresponding volume element to vanish, we can still use a method similar to the Hodge dual to construct a volume element $η_{Γ_0^+}$ for the light cone bundle. Based on this geometric structure $(Γ^+_0,η_{Γ_0^+},\hat h)$, we establish the fully covariant Boltzmann equation for photons. More importantly, the volume element and the induced metric are linked in a nontrivial way, which allows the physical distributions to be defined consistently. This yields the corresponding hydrodynamic quantities and their divergences, which take the same form as in the case of massive particles.

We study cosmological perturbations in multi-field inflation in which the scalar fields couple to the Gauss-Bonnet terms through a general coupling function $f(φ^a)$. We derive the complete quadratic action for the field perturbations in spatially flat gauge and find that, besides corrections to the kinetic, gradient and mass matrices, the Gauss-Bonnet coupling induces an antisymmetric velocity coupling that vanishes in Einstein gravity. Decomposing the field perturbations into curvature and isocurvature modes, we show that the effective mass of the curvature perturbation and the non-derivative mixing mass vanish identically, $\mathcal{M}_{\mathcal{R}}^2 = \mathcal{M}_{\rm mix}^2 =0$ which are justified by the Weinberg's adiabatic mode. We then compute the curvature power spectrum in two limits. First, when the isocurvature mode is heavy, it can be integrated over all scales. This yields an effective single-field theory of the curvature perturbation with a modified sound speed. Second, when the gradient of the coupling function is aligned with the background trajectory ($f_N =0$), the derivative mixings disappear and the superhorizon transfer of isocurvature into curvature perturbations is governed by the turn rate, which is modified by the Gauss-Bonnet corrections. We also present the tensor power spectrum and the tensor-to-scalar ratio at linear order in the coupling.

In the present paper, the cases of a uniformly accelerated reference frame and a Schwarzschild black hole in two and four dimensions are discussed focusing on formal aspects of the corresponding quantum scalar field theories. The main emphasis is made on building consistent quantum field theories and highlighting the problems that are not considered elsewhere. In particular, some pathologies in the case of a uniformly accelerated reference frame are discussed in detail: it is shown that some quantum states that are well defined in Minkowski spacetime possess an infinite energy when considered in Rindler spacetime. In the case of a Schwarzschild black hole in four dimensions, it is shown that the standard theory providing the Unruh vacuum state turns out to be inconsistent. The problem of consistency of quantum scalar field theory in Kruskal-Szekeres spacetime is also discussed.

The central challenge in microscopic quantum-gravity models is to derive geometric and gravitational observables from one refining quantum family without changing dynamics or fitting a new coefficient at each interface. We construct a controlled benchmark on a declared fixed physical Fourier band. For every spatial regulator, one finite-Hilbert interacting Hamiltonian defines a unitary completely positive process and a Euclidean Schwinger functional. Its spectrum and source derivatives determine the Lorentzian gap, connected four-response, mixed matter--geometry response, dynamic geometry kernel, induced Newton coefficient, and quantum-matter stress entering a semiclassical FLRW evolution. Two inequivalent spatial symbols obey analytic second- and fourth-order bounds and converge to the same limiting Hamiltonian. Duhamel and Weyl inequalities, a uniform gap, analytic perturbation theory, and a kernel-remainder estimate yield a common fixed-band regulator limit. Within the declared local covariant FLRW two-derivative sector, the metric source uniquely fixes a positive induced Newton coefficient rather than leaving it as an interface fit. Direct diagonalization independently yields a nonzero connected fourth response and a non-fitted dimensionless spectral combination. The excitation gap supplies a conserved stress tensor for a semiclassical Friedmann trajectory. This is a same-family fixed-band closure theorem, not an all-band QFT or completed quantum gravity. The two-mode band, finite oscillator cutoff, covariant two-derivative sector, and semiclassical Einstein equation are inputs. It does not construct a quantum metric Hilbert space, quantum BV measure, topology-change control, or autonomous sector selection. The result provides a reproducible benchmark and isolates the further theorems needed for nonperturbative closure.

Primordial non-Gaussianities are at their largest when the curvature perturbation is coupled to a massive field by a mixing too strong to be treated perturbatively. We show that this regime, naturally reached in multifield inflation, is solved analytically by the Laplace-space representation we recently introduced, in which every mode becomes a superposition of plane waves weighted by a single kernel. The bispectrum then follows as a rapidly convergent series of elementary functions, in perfect agreement with independent numerical results. We find that the cosmological collider signal at strong mixing is larger than widely believed, with only half of the usual Boltzmann suppression. Its amplitude further grows with the number of time derivatives of the curvature perturbation in the interaction, singling out the most promising observational targets.

The inner disk truncation radius, $R_T$, plays a crucial role in the regulation of star-disk interaction and the early evolution of star-disk-planet systems; however, measuring this parameter is observationally challenging. We present a new method for determining $R_T$ in young accreting systems that hinges on the color dependence of the accretion shock emission in multi-band time-domain surveys. Based on the accretion simulations of Robinson et al. (2017, 2021), we produce synthetic color-magnitude diagrams at near-UV and optical wavelengths that predict the loci of accreting stars as a function of $R_T$. We test these model predictions on young stars with interferometric $R_T$ estimates, finding very good agreement in our results. We apply this novel technique to a pilot survey of 26 classical T Tauri stars in Taurus and Upper Scorpius. We find a predominance of sources with small truncation radii, $R_T < 4\ R_\star$, and an overall distribution of $R_T$ that is statistically similar to that inferred from interferometric studies, while differing from those inferred from emission line modeling. Finally, we discuss the application of this technique to NASA's mission concept EVE, with the goal to provide simultaneous measurements of inner disk truncation radii, corotation radii and mass accretion rates for hundreds of young stars across the Galaxy. The unprecedented survey of inner disk properties that the mission would produce would enable the first stringent test of angular momentum evolution theories in young stars and reveal the impact of the inner disk conditions on early planet architectures.

Light spectator fields during inflation can acquire superhorizon fluctuations that cross a potential barrier between positive and negative regions of their potential. Motivated by the Standard Model Higgs instability, in this work we study the subsequent evolution of patches where this occurs in the radiation era after inflation ends for a Higgs-like spectator field. We utilise fully nonlinear, spherically symmetric numerical relativity. Across the black hole forming configurations in our investigation we find a robust two-stage evolution. First, the central negative potential region reverses its expansion, becomes kinetic dominated, and forms a primordial black hole that hides the runaway core. The positive potential barrier that survives outside this first horizon then determines one of two late-time branches. In the subcritical branch the original apparent horizon grows smoothly and engulfs the remaining scalar structure. In the supercritical branch, however, the potential energy of the barrier dominates the local evolution. The result is a transient wormhole throat, a bifurcating trapping horizon, and an inflating child universe branch. In both branches the parent radiation dominated universe is ultimately left with an ordinary primordial black hole whose subsequent growth is governed by radiation accretion.

We calculate the dynamical tidal response coefficients (TRCs) of 2+1D analogue black and white holes generated by draining and fountaining bathtub flows, respectively. The parameter space is characterized by the frequency and azimuthal number of the Fourier modes and the rotation of the analogue black or white hole. In general, the TRC is a complex-valued function of these parameters. Its real and imaginary parts are defined as the tidal Love number (TLN) and the tidal dissipation coefficient (TDC), respectively. The TRC of the analogue black hole (ABH) exhibits several interesting features. At certain points in the parameter space, including the static case of vanishing perturbation frequency, the TLN exhibits logarithmic running, while the TDC vanishes. Unlike in Einstein gravity, fluid dynamics allows for the physical existence of an analogue white hole (AWH) event horizon. Time-independent, torque-free, barotropic, inviscid, axisymmetric fluid dynamical equations can yield a pair of transonic background flow solutions. The solutions in the pair, corresponding to an ABH-AWH pair, share the same angular momentum and Bernoulli constant but have opposite mass flow rates, all of which are conserved quantities. For such an ABH-AWH pair, the TRC of the AWH is the complex conjugate of the TRC of the ABH.

Massive stars characteristically form in clustered environments. Characterising young massive 'protoclusters' is therefore crucial to constraining the mechanism(s) of massive star formation, and of the assembly of stellar clusters. We present 1.3 mm continuum results from the ALMA EGO-10 imaging survey, targeting ten Spitzer GLIMPSE Extended Green Objects (EGOs) - massive protostars with active outflows traced by extended 4.5 $μ$m emission. Our sensitive 1'.6$\times$1'.6 mosaics reveal rich protoclusters associated with all targets. With a mean spatial resolution 2200$\times$1600 AU, we identify 570 cores - between 13 and 135 per field. We quantify protocluster structure with the $Q$-parameter, finding structural diversity with 0.5 $\lesssim Q \lesssim$ 0.9. The sample is notable for the wealth of complementary high-resolution multiwavelength data available. Correlating our cores with these observations, we find only 2%, 5% and 4% of cores host 6.7 GHz CH$_3$OH masers, 22 GHz H$_2$O masers and cm-$λ$ continuum sources, respectively. The massive protostars traced by 6.7 GHz masers typically reside near protocluster centres (median offset 0.045 pc), and all at $d<$ 3 kpc are found in clustered locales, with $>$10 cores within 10,000 AU. Using VLA cm-$λ$ continuum observations, we construct a new evolutionary indicator: the ratio of protocluster cm-$λ$ continuum luminosity to the mass of the associated ATLASGAL clump ($L_\text{cm}/M_\text{AGAL}$). This ratio correlates positively with $Q$, with the correlation driven primarily by the cm-$λ$ continuum emission from MYSOs. This suggests dynamic protocluster structure, evolving from subclustered to centrally condensed, consistent with the global collapse in hierarchical, clump-fed models of massive star formation.

The Euler product of the Riemann zeta function is the partition function of a free bosonic gas whose mode energies are the logarithms of the primes. We show that this arithmetic gas, combined with the assumption that the entropy--geometry correspondence holds, leads to a quantum-corrected black-hole metric. The prime gas is a Hagedorn system. Its entropy is linear in the energy, which is exactly what an entropy linear in the horizon area requires, and the simple pole of the zeta function at $β=1$ fixes the coefficient of the logarithmic correction to the area law. The nontrivial zeros cannot play this role: their level density grows only logarithmically, and far too slowly to be extensive. Demanding that the reconstructed geometry reduce to Schwarzschild at large radius then fixes the map from arithmetic energy to horizon area and yields the closed-form metric $f(r)=1-2GMr/(r^{2}+\lz^{2})$ with $\lz^{2}=αG/π$. This describes a two-horizon black hole with Reissner--Nordström horizon structure but no Coulombic hair, a positive-energy anisotropic source obeying the null energy condition, a softened central singularity, a bounded Hawking temperature, and a cold extremal remnant that ends the evaporation. The same length scale follows independently from requiring that the first law hold exactly with the corrected entropy. The nontrivial zeros survive only as exponentially suppressed log-periodic ripples in the area, which suggests a physical interpretation of the Riemann hypothesis as the statement that arithmetic corrections to black-hole thermodynamics are as small as they can be.

Strong gravitational lensing contains key information about cosmic acceleration. Modern and next-generation galaxy imaging surveys are expected to provide high-quality data on $\mathcal{O}(10^5)$ galaxy-galaxy lensing systems. The plethora and complexity of the data are likely to present computational challenges for parameter inference methods for fitting high-dimensional likelihoods, which are often analytically intractable. Neural Ratio Estimation (NRE) efficiently computes individual likelihood ratios that can be combined into population-level posteriors. We use simulations to study the capacity of NRE to jointly predict the dark energy equation-of-state parameter $w$ and the total matter density $Ω_{m}$ from lensing images and companion spectroscopic information. We also introduce a post hoc posterior coverage calibration procedure that mitigates the model overconfidence that is typically found in neural density estimation applications. Our experiments show that the errors on both parameters decrease with increasing inference population sizes. In particular, for 100 lenses in a standard $Λ$CDM Universe, our calibrated NRE model achieves median fractional uncertainty of $22.8\%$ in $w$ and $2.9\%$ in $Ω_{m}$. This proof of concept demonstrates a potentially scalable approach for efficient cosmological parameter inference with large populations of galaxy-scale lenses observed in future surveys.

We construct a new family of static, electrically charged, cylindrically symmetric black-string solutions of four-dimensional Einstein--Yang--Mills theory with a negative cosmological constant, supported by a genuinely non-Abelian $SU(2)$ vortex field. The coexistence of electric and vortex sectors renders the standard single-function Lemos metric inconsistent with the Einstein equations, requiring a three-function metric compatible with the anisotropic Yang--Mills stress tensor. Regular near-horizon and asymptotic AdS expansions provide boundary data for numerical integration of the complete field equations. Two limiting sectors are recovered: vanishing electric horizon datum yields the neutral Lemos black string with zero Yang--Mills field strength, while switching off the vortex component gives the embedded $U(1)\subset SU(2)$ charged Lemos--Zanchin solution. Perturbations about the latter show that genuinely non-Abelian hair appears at linear order through the commutator field strength, whereas its backreaction on the electric profile and geometry begins at quadratic order. The horizon equations yield the analytic local non-degeneracy bound $R_1^{\rm crit}=s_0α\sqrt{3/(4πG)}$. The asymptotic geometry exhibits a quadratic correction absorbable into an effective Abelian charge and a cubic coefficient giving the first irreducible asymptotic signature of the non-Abelian hair. The entropy obeys the Bekenstein--Hawking area law, while the temperature difference from the embedded Abelian solution arises through the asymptotic lapse normalization. The fully backreacted numerical solutions are nonlinear realizations of this transverse non-Abelian deformation.

Recently we have demonstrated that black hole radiation can be smoothly extended to the zero-temperature limit. Here we analyse entropy variations at the zero temperature limit and evaluate the validity of the Generalized Second Law (GSL). Typically one verifies that the entropy gain in the exterior region compensates for the entropy loss of the black hole due to evaporation. However an anomalous feature emerges at the zero temperature limit. The fact that the emission is restricted to the modes $ω<mΩ$, renders the variation of the horizon area of a Kerr black hole positive definite. Simultaneously the von Neumann entropy of the exterior region increases both for bosonic and fermionic modes. The dual entropy gain at the zero temperature limit marks a robust conservation of the GSL. We demonstrate that the dual increase also applies to Kerr-Taub-NUT black holes whose entropy is neither directly proportional to the horizon area nor analytically integrable.

We present a spacetime that describes a rotating and charged black hole immersed in an external magnetic field. This exact solution to the Einstein-Maxwell equations has a surprisingly simple form which reduces to standard Kerr-Newman metric, and to the uniform Bertotti-Robinson magnetic universe. We show that the Kerr-BR solution presented in [Phys. Rev. Lett. {\bf 135}, 181401 (2025)] has a special charge, and we identify a genuine Kerr black hole without any charge, in the magnetic field. In the weak-field limit it gives the Wald solution. We analyze main physical properties of the new family, namely singularities, horizons, regularity of axes, charges, entropy and temperature.

We formulate a data-completeness test for extracting low-energy characteristic-cone running from flat-space on-shell renormalization-group residues. In four-dimensional, parity-even Einstein--scalar--Maxwell effective field theory through four derivatives and four external fields, Wilson directions are first placed on the physical EFT quotient and background propagation is represented by a traceless response on the full massless root sector. The response is identifiable from chosen amplitude coordinates precisely when the amplitude kernel lies in the response kernel. For the published opposite-helicity vector--scalar residue we find $β_μC_{γφ}^{(2)}=-55/(96π^2)$ and, after an order-reduced frame translation, the contribution $β_μ(r_T-r_γ)=-55/(48π^2M_{\rm P}^4)$. A finite operator registry gives explicit rank-two incompleteness tests: independent $F^4$ directions split photon polarizations on a magnetic background, while $\mathcal CFF$ produces an analytic splitting along a principal Weyl axis on a Petrov-I Einstein--scalar background. For the isolated scalar--photon contact, direct Schur reduction finds nonzero raw mixing but no first-order physical off-diagonal entry, by the Maxwell Ward identity. Thus the residue fixes one diagonal cone contribution but does not determine the term-complete characteristic response over the declared admissible coefficient domain.

We present a multi-sector $TESS$ photometric and spectral energy distribution analysis of the eclipsing binary CPD--82~291 (RX~J0853.1--8244). Eight $TESS$ sectors reveal well-defined eclipses and substantial, changing out-of-eclipse modulation. The asymmetries are reproduced by different surface-brightness configurations at different epochs, consistent with strong and evolving stellar activity, although the fitted circular regions are not unique maps of individual starspots. A detached-binary light-curve model yields an orbital inclination of $i\simeq82.1^\circ$ and a photometric mass ratio of $q_{\rm phot}\simeq0.464$. Without a spectroscopic radial-velocity orbit, the mass ratio and absolute parameters are model-dependent rather than dynamical measurements. Combining the light-curve geometry with the SED radius scale and \textit{Gaia} distance gives $R_1\simeq1.25\,R_\odot$ and $R_2\simeq3.70\,R_\odot$; the corresponding model-inferred masses remain substantially less secure than the dimensionless light-curve parameters. The stellar photospheres reproduce the observed SED from the optical region through WISE W4, with no statistically significant mid-infrared excess. The only positive long-wavelength deviation is a low-significance $60\,μ$m measurement, while the $90\,μ$m result is an upper limit. A cool circumbinary-dust model can reproduce this tentative signal, but such a component is not required by the photometry, and its inferred mass and radial extent are non-unique and strongly model dependent. Evolutionary tracks place the primary close to the zero-age main sequence and the secondary well above the main-sequence locus. However, the absence of radial velocities and spectroscopic youth indicators, together with $Gaia$ astrometry inconsistent with Chamaeleon~I membership, prevents a secure pre-main-sequence classification.

Recent DESI DR2 observations have renewed interest in extensions of the $Λ$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 $ρ_{\rm dm}\propto a^{-3+δ}$, with the constant $δ$ 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σ$. 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 $Λ$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.

The Stone-von Neumann theorem makes canonical quantization unique only after one imposes strong continuity of the Weyl translation groups. We show that this continuity assumption cannot be certified by any finite-resolution interrogation of the Weyl relations. A finite protocol probes only finitely many phase-space displacements and bounded expectation values, thereby determining a protocol algebra \(C_S\) generated by a finitely generated subgroup of Weyl translations, but not the regularity of an extension to the full Weyl algebra. Using the polymer representation as a canonical non-regular representation, we prove that for every regular Weyl representation, every normal state, and every finite family of bounded protocol observables, a polymer state reproduces the same expectation values to arbitrary prescribed accuracy. Thus finite agreement with Schrödinger quantum mechanics does not operationally certify the regularity hypothesis entering uniqueness. The result does not assert physical equivalence between regular and polymer quantizations, it identifies where additional input, such as Hamiltonian dynamics, energy regularity, semiclassicality, or a continuum limit, is required to distinguish or select representations. In this precise sense, non-regular polymer kinematics are not excluded by finite Weyl data alone, a point of direct relevance to quantum-gravity motivated quantizations.

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 $Λ$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.

We propose a thermodynamic description of the genuine Kerr-Bertotti-Robinson (Kerr-BR) black hole spacetimes. Contrary to previous studies, such solutions no longer have electric charge. Consequently, the thermodynamics is governed by simple relations. Namely, upon introducing a suitable normalization factor for the timelike Killing vector, these spacetimes are characterized by the standard Hawking temperature, Bekenstein's entropy, Komar mass and angular momentum, and horizon angular velocity. The corresponding mass is shown to obey the standard Christodoulou-Ruffini formula, along with the standard first law and Smarr relations. In fact, by a proper rescaling of coordinates and solution parameters, and up to a non-trivial parameter space constraint, we formally recover the thermodynamic quantities of the vacuum Kerr solution. This is used to bring the genuine Kerr-BR spacetime into a Kerr-like form.

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.

We study supersymmetric near-horizon geometries of gauged D=5 supergravity coupled to vector multiplets, on the branch where the canonical rotational Killing vector $\tilde V$ of the horizon section $\mathcal{S}$ is non-vanishing. No rotational symmetry is assumed. On a compact, connected $\mathcal{S}$ without boundary, a second rotational Killing field, independent of $\tilde V$, always exists on the open set where $Z$ and $W$ are linearly independent, and extends to an isometry of all of $\mathcal{S}$, unconditionally if the moduli are constant and under a checkable extension hypothesis otherwise. With $K=Q_{IJ}C^IC^J$ vanishing exactly in the minimal theory, $K\equiv0$ recovers the known result of Grover, Gutowski, Papadopoulos and Sabra. When $K>0$, $P=Z+W$ vanishes at every critical point of $\|η_-\|^2$, and two sub-branches occur, both realised. On the first, $P\equiv0$, the moduli are constant and the isometry follows by homogeneity, recovering Kunduri and Lucietti local geometries as a consequence of supersymmetry and compactness rather than an ansatz. On the second, $P\not\equiv0$, the moduli genuinely vary; closure of the first-order system ($dP=0$) turns the Killing equation for a vector orthogonal to $P$ into a flat rank-two system transported along $P$, solved explicitly by $U=\|η_-\|^2(h_YN-h_NY)$. Both sub-branches are non-empty on compact $\mathcal{S}$: known $U(1)^3$ black holes populate the first, and an explicit compact horizon from the non-static Kunduri-Lucietti family, with $P$ vanishing at two points, populates the second and verifies the extension hypothesis numerically. We also prove $K$ is identically zero or everywhere positive, exclude a hidden Killing-tensor symmetry, and leave the $S^1\times S^2$ black-ring window $\tfrac13Φ^2<χ^2K\leq\tfrac43Φ^2$ unexcluded on either sub-branch.

This work represents an initial step towards a systematic framework for identifying hidden symmetries in algebraically general (Petrov type I) spacetimes. Motivated by this, we focus on Weyl's class, which provides a simple yet a class of algebraically general solutions of static, axisymmetric geometries. Employing the canonical forms of rank-2 Killing tensors, we systematically derive the constraints that must be satisfied for these spacetimes to admit genuine hidden symmetries. Our analysis yields an analytical characterization of the obstructions to integrability within subclasses of Weyl's class, both in vacuum and electrovacuum. In particular, we find that no two-Killing-vector member of Weyl's class, in vacuum or electrovacuum, admits an irreducible Killing tensor capable of supplying the additional integral of motion required for complete integrability. Members of Weyl's class admitting a third additional Killing-vector isometry (i.e. Levi-Civita) are already completely integrable through their manifest symmetries alone admitting only reducible Killing tensors. This result provides an analytic complement to previous analytical and numerical investigations, and establishes a systematic connection between the absence of hidden symmetries and the breakdown of complete integrability. Beyond the specific classification obtained here, these results demonstrate the potential of our approach regarding the canonical forms of Killing tensor as a systematic tool for probing hidden symmetries in Petrov type I spacetimes.

In particle physics the virtual N$^3$LO contribution to gluon-fusion Higgs production in the Standard Model requires four-loop multi-scale QCD amplitudes that in dimensional regularization are dominated by singular integral reduction and difficult master-integral evaluations. In this note we show that a nonlocal UV completion makes each fixed-order four-loop Feynman diagram fully convergent in four dimensions while keeping the analytic structure needed for Lehmann-Symanzik-Zimmermann (LSZ) reduction formula. For gluon-fusion $gg\!\to\!H$ we derive a closed form Schwinger-parameter representation where all loop-momentum integrals are performed analytically and thus leaving finite parameter integrals viable for direct numerical evaluation without needing the usual UV subtractions. This method heps to isolate the collider-hard contribution as a convergent four-loop form factor. The strong coupling then is defined by a finite normalization condition in the entire-function regularization scheme and it is related to the $\overline{\mathrm{MS}}$ coupling by finite matching where after this matching, the regulator-dependent corrections decouple in the local limit $E_M\!\to\!\infty$.

We study the one-loop relation between an effective-field-theory extension of Poincaré gauge theories of gravity with propagating torsion and its unimodular counterpart. In the Einstein representation, we consider two representative quadratic torsion subsectors, namely totally antisymmetric (axial) torsion and the vector component of hook-antisymmetric torsion. On maximally symmetric metric backgrounds with vanishing background torsion, the diffeomorphism-invariant and genuine unimodular theories yield identical one-loop determinants in both sectors, reproducing the same local logarithmic divergences. We further examine flat backgrounds with homogeneous axial or vector torsion, where metric-torsion mixing is present. In each case, the torsion-dependent one-loop effective action in unimodular gauge and in the genuine unimodular formulation is controlled by the same physical determinant. Thus, for the actions and backgrounds considered, the unimodular constraint does not modify local one-loop torsion dynamics. We also delineate the limitations of this equivalence and the extensions needed for generic torsionful curved backgrounds.

We investigate the Joule-Thomson expansion of a non-commutative Schwarzschild-AdS black hole, whose space-time geometry is described by the Lie-algebraic $κ$-deformed space-time, and show that the $κ$-deformation parameter induces Joule-Thomson expansion in an uncharged-AdS black hole. We find the ratio of minimum inversion temperature to critical temperature $T_i^{(min)}/T_c \simeq 0.493$, is independent of the deformation parameter, and is in close agreement with the value for Reissner-Nordstrom black hole. Furthermore, we show the existence of a minimum black hole mass (depending on the $κ$-deformation parameter) below which the Joule-Thomson expansion vanishes.

The study presents a confirmation of six quadruples with two sets of eclipses that have the 2+2 architecture. These so-called doubly eclipsing systems still present a quite rare group of stars. We collected all available photometric data and carried out a detailed analysis of them. In addition to the precise TESS photometry used to model the light curves for both inner eclipsing binaries, photometric survey data were also used, and more than 100 nights of our own dedicated observations were carried out. These were mainly used for the detection of the long-term evolution of orbital periods. Thanks to these data, we were able to derive the mutual orbits of the inner doubles via eclipse timing variations. The systems studied are: CzeV1254 (periods 0.715348 d + 0.362593 d, mutual period 15.6 yr); ASASSN-V J070838.27-171952.9 (4.300490 d + 3.767235 d, 4 yr); ASASSN-V J091951.17-593306.9 (2.668674 d + 14.342218 d, 0.95 yr); WISE J100820.07-731554.0 (0.368315 d + 7.945339 d, 3.8 yr); ASASSN-V J143536.01-721459.4 (7.353842 d + 6.368567 d, 27.4 yr); and Corot 310284765 (1.875418 d + 2.371126 d, 22.6 yr). Some of the binaries show significantly eccentric orbits. We also estimated their apsidal motion periods. With our six new systems, the number of confirmed orbits of 2+2 quadruples has been increased to 64 in total.

At a deeper level the principle of least action is interpreted as the law of least entropy increase consistent with Prigogine's principle of minimum entropy production, and the implications of quasistatic information quantization rule (Proc. R. Soc. A 480: 20240024), are explored for the conjectured second law of quantum complexity. It is thus shown that the conjectured second law of complexity is derivable from the information quantization rule such that long after heat-death the quantum state complexity evolves as $C(t)=C_{\rm max}\exp(-1/t)$, increasing with time to a saturation value $C_{\rm max}$ that is exponential in the equilibrium entropy $S_{\rm max}$. For the out-of-equilibrium circumstances, however, the quantum complexity can decrease with the time, asymptotically tending to a minimum determined by the distance from equilibrium.

Gravitational-decoupling provides a systematic framework for constructing black hole geometries from known (seed) solutions (such as the Schwarzschild), whose exterior properties closely resemble those of the seed solution while still retaining potentially distinguishable features. Using this approach, we study a class of static, spherically symmetric hairy black holes and examine how such deformations affect their horizon structure, photon spheres, shadow scales, and scalar quasinormal modes (QNMs). We first determine the physically controlled regions of the parameter space by imposing the relevant horizon and energy condition requirements, while retaining configurations outside these domains only for exploratory comparison. To isolate geometric effects, we compare solutions at a common horizon radius, while accounting explicitly for the distinction between the mass parameter entering the seed geometry and the ADM mass measured asymptotically. The resulting ADM normalized shadow scales are used for an illustrative comparison with the characteristic angular scales of M87* and Sgr A*. We also compute massless scalar QNMs using the sixth-order Wentzel-Kramers-Brillouin approximation and the improved Asymptotic Iteration Method, finding close agreement over the modes considered. Our results show that fixing the horizon radius does not uniquely determine either the optical or perturbative properties of the spacetime, and that mass normalization can significantly alter the interpretation of geometric trends. Horizon structure, null geodesic observables, and scalar QNMs therefore provide complementary probes of gravitationally-decoupled black hole geometries.

The Sun is one of the fundamental benchmarks in astronomy, and there is ever growing interest in precise and accurate determinations of its chemical composition. We present the current state-of-the-art in spectroscopic determinations of the solar elemental abundances, that now routinely employ three-dimensional (3D) radiative-hydrodynamic simulations of the solar photosphere and post-processing in non-local thermodynamic equilibrium (non-LTE). We critically review the recent literature and present recommended present-day and protosolar abundances of the 83 long-lived elements, and discuss their implications vis-à-vis primitive meteorites, helioseismology, neutrino fluxes, and the Solar Modelling Problem.

As a natural extension of the single-field fuzzy dark matter (FDM) model, multifield FDM has attracted increasing attention in recent years. This scenario is motivated both by the axiverse scenario predicted by string theory and by the possibility that multifield FDM may provide a better match to astrophysical observations than its single-field counterpart. In this work, we perform high-resolution numerical simulations to systematically investigate the dynamics of wave structures in multifield FDM halos. In particular, we study the oscillatory and stochastic motions of the central core, the evolution and statistical properties of granules, and the resulting dynamical heating of embedded stellar systems. We find that the frequency spectra of the core density oscillations develop multiple peaks and shift toward higher frequencies relative to the single-field case. The centers of different field components undergo nearly synchronized random walks, while subdominant components exhibit larger random-walk amplitudes. We further show that the suppression of granule density fluctuations with increasing number of fields is largely insensitive to the fractional abundance of each component over a broad parameter range. Moreover, using self-consistent simulations, we find that the dynamical heating induced by granules is progressively suppressed as the number of fields increases. However, once the contribution from the central core is taken into account, this trend would become much less pronounced.

Motivated by the old idea of treating gravity as a gauge theory, we review the gauge-theoretic construction of Conformal Gravity and of Noncommutative (Fuzzy) Gravity. We then outline how these gauge formulations of gravity can be unified with internal interactions described by an SO(10) Grand Unified Theory, based on the critical observation that the dimension of the tangent-space symmetry group of a curved manifold need not coincide with the spacetime dimension.

The deflection angles of rays diverge logarithmically near outside and inside of a photon sphere and they could be detected by near-future space observations on black hole shadows. We show that the deflection angles of the rays near a photon sphere at the throat of an example among some commonly used types of wormholes or one of the simple wormholes, which is often studied by the Event Horizon Telescope collaborations, does diverge not logarithmically but in power. This is because it has a marginally unstable photon sphere at a $Z_2$-symmetric wormhole throat and the wormhole has not a common nor simple but a special property. We investigate numerically gravitational lensing of rays with the deflection angles diverging in power outside and inside of a marginally unstable photon sphere at $Z_2$-symmetric wormhole throat in strong deflection limits. We apply our numerical method to the commonly-used simple wormhole spacetime, a Damour-Solodukhin wormhole spacetime, and a Simpson-Visser black-bounce spacetime. Our numerical results imply universal property of the deflection angle by the marginally unstable photon sphere on $Z_2$-symmetric wormhole throat in the strong deflection limits as similar to the deflection angles of ray bent by a marginally unstable photon sphere off a throat in the strong deflection limits which have a different power. We also correct the coefficients of power-divergent terms of deflection angles of the rays bent nearly outside of the marginally unstable photon sphere at the $Z_2$-symmetric wormhole throat in a semi-analytic approach investigated by the author and others.

We prove the black hole uniqueness conjecture in the axially symmetric, stationary, vacuum setting: there is no regular asymptotically flat configuration with more than one horizon component when every degenerate component has nonzero angular momentum. More precisely, it is shown that in any such multi-horizon configuration, the logarithmic angle defect along every finite axis rod is strictly negative, and hence the interaction force along such axis rods is always attractive. The proof is based on a refined asymptotic analysis of the associated singular harmonic maps, in part from the companion paper [16], together with a global bound for the Weyl conformal factor arising from a certain scalar curvature related differential inequality. Furthermore, as a corollary of the extremal version of this uniqueness result, we obtain the conjectured mass-angular momentum inequality for multiple dynamical black holes using the main theorem of [17].

I prove that a collapsing charged scalar field with scalar charge $\mathfrak{e}$ that forms an exactly extremal Reissner--Nordström black hole à la Kehle--Unger with radius $r_+$ must satisfy $|\mathfrak{e}|r_+ > \frac{1}{3}$. In contrast, I also show that no such bound exists for subextremal collapse: a scalar field with an arbitrary $\mathfrak{e} \ne 0$ can form a subextremal black hole with an arbitrary (subextremal) charge-to-mass ratio. Complementing the extremal bound, I prove that a Schwarzschild black hole can become extremal if $|\mathfrak{e}|r_+ > \sqrt{\frac{3 + \sqrt{33}}{3}}$. The fact that $|\mathfrak{e}|r_+$ can be taken to be of order unity could be considered evidence that the third law of black hole mechanics in vacuum is false.

In multi-fluid, matter-type gravity, the Bianchi identity constrains the divergence of the total effective stress tensor but does not, by itself, determine the currents assigned to its constituent sectors. Those currents are selected only after the off-shell matter action and equations, or an additional phenomenological closure, have been specified. We formulate this distinction and examine a spatially flat two-fluid background inspired by scale-independent EMSG. The case study combines an algebraic perfect-fluid prescription with a vanishing contracted-Hessian contribution, together with separate conservation of the total conventional and modification sectors. Because neither assumption is derived here from a concrete off-shell fluid action, the resulting system is an effective background closure rather than a microscopic EMSG model. For unequal EoS and every $α\ne0$, the closure yields a Barrow-Clifton system whose transfer coefficients depend on $(w_1,w_2)$ but not on $α$. Hence the nonzero-$α$ family is singular: its $α\to0$ limit does not recover the uncoupled GR conservation laws. We solve the two density eigenmodes and derive the associated modified Liénard equation for $H$. We also prove that the discriminant governing rank loss of the density-reconstruction map is strictly positive for every finite $w_1\ne w_2$. Thus each genuinely quadratic case has two distinct real rank-degenerate couplings. Exact vacuum and stiff-fluid families illustrate modal cancellation. When today's conventional densities are positive, the intervals on which both remain positive generally terminate at finite endpoints in the parameter ranges analyzed explicitly. These results provide a consistency diagnostic for separating action-level predictions from closure artifacts in multi-fluid, matter-type gravity; they do not establish an observationally viable EMSG model.

Light spectator fields provide a compelling inflationary origin for scalar dark matter, but their primordial fluctuations often generate excessive isocurvature perturbations on large scales. We develop a novel, exactly solvable filtering mechanism based on a tower of $δ-$function contributions to the spectator mass, which we call kicks, whose dense limit converges to a smooth Coulomb-like profile. The filter produces a isocurvature spectrum during inflation that is strongly suppressed on large scales, rises with momenta as $k^3$ in the infrared, and, after a rapid turnover, smoothly approaches a constant plateau at small scales. This construction is motivated by a discrete translational invariance, admits an EFT regime in which radiative corrections can be systematically controlled, and can find a concrete realization in terms of an inflationary massive-clock model. We embed the filter in a minimal dark-matter scenario in which a spectator condensate oscillates about the minimum of a quadratic potential and behaves as pressureless matter. In this framework, the spectator dark-matter mass is directly related to the turnover scale of the inflationary isocurvature spectrum. We determine the resulting dark-matter abundance and analyze the relevant isocurvature and backreaction constraints. The most restrictive realization, in which the spectator condensate is generated during inflation, favors a relatively low inflationary scale.