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

Recent JWST observations have revealed a population of compact, optically red sources known as Little Red Dots (LRDs). A popular interpretation is that LRDs host massive black holes embedded in dense gaseous envelopes, yet the physical origin of such envelopes remains unclear. We propose that the optically thick outflow driven by supercritical accretion onto black holes may explain the envelope. Based on an analytic radiative hydrodynamic outflow model, we relate the outflow properties to the black hole mass $M$ and dimensionless accretion rate $\dot{m}$. With $M\sim10^5-10^7\,M_\odot$ and $\dot m \sim 1500-5000$, the outflow photosphere reaches luminosities of $10^{43}-10^{45}\,\rm erg\,s^{-1}$ and temperatures of $\sim 3000-6000$ K, effectively matching the red optical continua observed in LRDs. The presence of a thick scattering region beyond the photosphere is central to deciphering the distinctive properties of LRDs. For each object, if one derives $M$ and $\dot{m}$ from the observed luminosity and temperature, while accounting for both kinematic broadening and scattering effects, the model predicts an emission line FWHM consistent with observations within a factor of 1.5 for more than 80\% objects. Furthermore, with Cloudy simulations, we find that the partially ionized gas beyond the photosphere produces Balmer breaks broadly consistent with measurements.

Solar wind alpha particles exhibit preferential heating and acceleration relative to protons; however, their behavior in the vicinity of turbulent coherent structures remains less understood. We report the first evidence of localized alpha particle and proton heating within coherent structures identified using the Partial Variance of Increments (PVI) method, based on Parker Solar Probe (PSP) observations. Our results show that high-PVI events are associated with significant, species-dependent temperature enhancements: protons undergo a relative larger temperature increase than alpha particles. This preferential proton heating produces a localized decrease in the alpha-to-proton temperature ratio, indicating that the plasma is driven toward thermal equilibration between species. The heating is also anisotropic, being dominated by enhancements in the perpendicular temperature. These temperature-signatures coincide with a pronounced reduction in the normalized alpha-proton differential flow speed and a localized minimum in the Coulomb collision age, suggesting that the relaxation is affected primarily by collisionless kinetic effects. These findings provide new insight into the intermittent energy conversion and ion thermodynamics in the solar wind.

Recent high-resolution observations from the Solar Orbiter mission can help detect the indirect signatures of heating at the smallest scales. In this work we measure the properties and investigate the physical origin of propagating disturbances (PDs) in the intensity at the smallest resolvable scales with Solar Orbiter/EUI along coronal loops. We used two sequences of EUI/HRIEUV at high spatial (down to 125 km per pixel) and temporal resolutions (5 s of cadence). We placed slits along 13 active region (AR) coronal loops. We measured the plane-of-sky (PoS) velocities and the intensity perturbation damping along the slits of PDs. We also measured the periodicity of PDs in one slit by using a Fourier analysis. We report the detection of PDs that have high PoS velocities ranging between 500 and 2000 km/s (which we call "fast" PDs). They are only visible in the upper part of the coronal loops. The intensity increase associated with these fast PDs is on the order of 4\% to 8\% of the HRIEUV intensity, and we measured little to no damping of their intensity with distance. We also measured a peak in the Fourier power spectra at 2 min above the 95\% confidence limit that is associated with fast PDs. The fast PDs are detected in the same coronal loops as PDs with a lower PoS velocity (70 to 90 kms/s), which we refer to as "slow" PDs. Unlike the fast PDs, these slow PDs are only visible in the lower part of the coronal loop, and they show clear intensity damping. Slow PDs show properties consistent with slow magneto-acoustic modes or upflows. On the other hand, fast PDs cannot be explained by slow magneto-acoustic modes. Instead, they show properties consistent with fast flows induced by magnetic reconnection, current sheet generated by propagating transverse oscillations, and fast magnetohydrodynamics modes or Alfvén waves.

We present optical observations on the young type Ia supernova (SN Ia) SN 2021fxy obtained within a few days after the explosion, with a focus on its prominent high-velocity features (HVFs). It reached a $B$-band maximum of $M_{\rm max}(B) = -19.36\pm0.31$ mag, corresponding to a bolometric luminosity of $\sim 1.3\times10^{43}~\rm{erg~s^{-1}}$ with a synthesized $^{56}$Ni mass of $0.58\pm0.14$ M$_{\odot}$. The early spectra exhibit strong HVFs of intermediate-mass elements that are significantly detached from the photospheric components. In particular, the velocity of the Si II $\lambda6355$ HVFs follows a power-law evolution ($β\approx 0.1$), shallower than the expected photospheric velocity evolution expected for an assumed $n=10$ density profile ($β\approx 0.22$) under homologous expansion. This behavior is consistent with the HVFs forming in intrinsic ejecta structures at least partially decoupled from the bulk outer ejecta, providing a possible constraint on the explosion physics of SNe Ia.

The growth of galaxies in the early Universe is thought to be dominated by compact, intensely star-forming systems, yet the corresponding growth of supermassive black holes (SMBHs) within such environments remains poorly constrained. Green Pea galaxies are nearby analogs of rapidly assembling galaxies in the early Universe owing to their compact morphologies, intense star formation, low metallicities, and extreme ionization conditions. Although galaxy interactions are thought to trigger episodes of rapid SMBH growth, direct observations of simultaneous accretion onto multiple SMBHs in compact, intensely star-forming galaxies are lacking. Here we report the discovery of the first confirmed dual active galactic nucleus (AGN) in a Green Pea system, SDSS J162209.41+352107.5. Chandra imaging resolves two luminous hard X-ray sources separated by 8.4 kpc in projection, demonstrating simultaneous accretion onto two SMBHs. Follow-up Keck spectroscopy confirms that the two optical nuclei share a common redshift and independently exhibit broad Balmer emission and high-ionization AGN emission lines. Unlike most known dual AGN, which are typically found in massive mergers, J162209.41+352107.5 is a compact low-mass system analogous to galaxies thought to dominate early phases of galaxy assembly. These findings demonstrate that efficient growth of multiple SMBHs can occur in such environments and establish Green Pea galaxies as nearby laboratories for investigating the interplay between galaxy interactions, star formation, and black-hole growth under conditions analogous to those prevalent in the young Universe.

Gravitational waves (GWs) have proven to be a powerful probe of compact binary populations. In the millihertz frequency range accessible to Laser Interferometer Space Antenna LISA, binary white dwarfs (BWDs) are expected to constitute a dominant source, forming both individually resolvable signals and an unresolved Galactic background. In this work, we construct a Milky Way like population model and calculate the GW background from unresolved Galactic BWD in the LISA sensitivity range, with particular emphasis on exploring and constraining uncertainties in binary stellar evolution. We employ COMPAS binary population synthesis framework to generate synthetic populations of BWD in the Milky Way. Various physically motivated evolution prescriptions and initial model parameters are used to study diverse population of BWDs. From these populations, we construct the GW background and investigate the dependence of the background spectrum on the assumptions on binary analysis. We discuss the possibility of constraints on binary evolution that LISA GW observations may yield. We find that the shape and amplitude of the background are sensitive to key binary evolution parameters like common envelope evolution and mass-transfer efficiencies. Variations in these assumptions lead to measurable differences in the predicted background spectrum. Our results demonstrate that LISA observations of the unresolved BWD background have the potential to constrain binary evolution models. This highlights the importance of GW background modelling as a complementary tool for studying the formation and evolution of compact binaries in the Milky Way.

Simulation-based inference (SBI) has become a major tool for extracting cosmological information from weak-lensing (WL) surveys, particularly from non-Gaussian observables. We compare its two main paradigms: explicit likelihood inference (ELI), based on a Gaussian likelihood built from an emulator and covariance matrix, and likelihood-free inference (LFI), which learns the likelihood directly from simulations using neural density estimators. Using Gaussian random field mocks representative of the non-tomographic final Euclid data release, we analyse shear two-point correlation functions (shear-2PCFs), compressed with linear or non-linear methods, together with a fundamentally different map-level convolutional neural network (CNN) statistic, focusing on $Ω_{\rm m}$ and $S_8$. We deploy posterior calibration diagnostics developed for LFI, including the test of accuracy with random points (TARP), showing that ELI becomes strongly miscalibrated under emulation inaccuracies or likelihood non-Gaussianity, whereas LFI remains well calibrated. These effects drive substantial disagreement between ELI and LFI, which largely vanishes once addressed. We further show that the compression scheme can significantly degrade ELI while leaving LFI largely unaffected. Although shear-2PCFs should capture all the information in Gaussian fields, finite compression and non-Gaussian likelihoods cause ELI constraints to differ by up to a factor of two from those inferred with the CNN, while the discrepancy drops to $\approx 30\%$ for LFI, underscoring the robustness of the deep-learning probe. Overall, our results indicate that in our simple setup, which neglects systematic biases, LFI provides a more robust and better-calibrated framework, while highlighting accurate non-Gaussian likelihood modelling and posterior calibration diagnostics as essential for future ELI analyses.

Cosmic Domain Wall networks are among the most interesting sources of a stochastic Gravitational Wave (GW) background from the early Universe. We present a thorough analysis of their annihilation, with a focus on scenarios where the collapse is induced by a population bias, whereby one of two degenerate vacua is initially preferred over the other. Our state-of-the-art $3+1$ lattice field theory simulations in the expanding Universe reveal that the network decays around the temperature $T_\text{ann}\sim T_s\,\mathcal{B}_s^{0.8}$, where $\mathcal{B}_s$ quantifies the preference for one vacuum over the other at the onset of the scaling regime at the temperature $T_s$. Furthermore, we obtain the spectrum of GWs from such networks, and provide a detailed comparison with the alternative potential bias annihilation mechanism that relies on a small explicit symmetry breaking in the potential. En passant, we update results on the evolution of these networks and on their GWs, and clarify existing disagreements in the recent literature. Our results sharpen the phenomenological viability of spontaneously broken discrete symmetries, and provide GW spectra that Pulsar Timing Arrays (PTAs) and ground-based interferometers (LIGO-Virgo-KAGRA) can readily use in their searches for a cosmological GW background.

The fourth observing run (O4) of Advanced LIGO, Virgo, and KAGRA presents unparalleled sensitivity, rendering unsupervised pipelines highly vulnerable to the non-stationary domain shift of the detectors' noise manifolds. We present DANTE V3, concluding a longitudinal investigation into unmodeled anomalies during early O4a. By expanding to a multi-scale geometric framework (0.5 s to 4.0 s), we amplify morphological sensitivity, extracting 10,372 unique candidates. To mitigate domain-shift artifacts, we introduce a block-bootstrap Domain Shift Defense (DSD) against a vector-quantized native background index. While 28.3% of candidates survive recalibration, global topological analysis reveals they lack discrete morphological cohesion. The survivors coalesce into a single macro-cluster without compact substructure; we show explicitly that morphological "diffusivity" comparisons used previously are confounded and retire them. Pristine background is even more monolithic (100% in one cluster vs 99.77% for survivors): this topology is a property of the embedding geometry, not of anomaly status. Executing a definitive physical environment monitoring (PEM) cross-correlation defense using an empirically calibrated null, one singleton is vetoed by a control-line coupling, while another survives as an uncatalogued instrumental outlier. Replacing the embedding-similarity cross-detector coincidence test with a physical normalized cross-correlation test, we find no coincident events among 8,749 candidates. We quote 90% frequentist Poisson upper limits on the rate of novel uncatalogued instrumental morphologies -- $R_{90} \le 5.83 \mathrm{yr}^{-1}$ (H1) and $R_{90} \le 5.63 \mathrm{yr}^{-1}$ (L1). This underscores the absolute necessity of native background recalibration and physical auxiliary vetoes in unsupervised gravitational-wave astronomy.

The parameterized post-Newtonian (PPN) approach is the state of the art formalism for performing theory independent tests of weak-field gravity, and for constraining possible deviations from Einstein's theory. Within this framework, global conservation laws are useful for the calculation of dynamics and for giving meaning to parameters. In this paper we extend the concept of semi-conservative and fully-conservative theories of gravity to include situations in which compact astrophysical bodies are modeled as masses that are sensitive to their local environment, as relevant for theories that violate the strong equivalence principle. We find that globally conserved quantities can still exist in the presence of such sensitivities, and find their explicit forms when they do. We identify new ways of writing the coefficients that enter into the PPN metric when a theory of gravity admits conserved quantities in the presence of a sensitive body, and demonstrate the applicability of our approach by comparing it to known results in scalar-tensor theories of gravity.

We perform a search for continuous gravitational waves from J0435+3233 using LIGO O4a public data. J0435+3233 is unique among millisecond pulsars as it exhibits an exceptionally large spin-down and marks the first pulsar observed to date with a spin-down larger than $10^{-12}$ Hz/s in the sub $10$ ms spin period range, making it a potentially strong source of continuous gravitational waves. We target signals at exactly twice the rotation frequency, a narrow band around this frequency, and also signals corresponding to r-modes. Our results are consistent with a non-detection. Our most stringent upper limit on the intrinsic gravitational wave amplitude is $h_0=\targetedULninetyfive$, at the 95\% confidence level, which is a factor of $\approx14$ times lower than the spin-down limit. With an estimated source distance of $1.2$ kpc our upper limit on the gravitational wave amplitude excludes source ellipticities greater than $\targetedellipULninetyfive$, making this the first source for which the spin-down upper limit is beaten by over an order of magnitude {\it{and}} the ellipticity is constrained to a physically interesting range, in the low $10^{-8}$ region.

This study investigates the dynamical behavior of particle's test around a spinning charged Einstein-Maxwell-dilaton (EMd) quantum corrected black hole (QCBH) by examining the effects of spinning parameter $a$, quantum correction parameter $b$, and charge parameter $Q$. Using EMd spacetime geometry, we analyze the effective potential and effective force to determine the stability and structure of the circular orbits in a strong gravitational field. Furthermore, we calculate the oscillation frequencies radial $(Ω_r)$, vertical $(Ω_θ)$, and angular $(Ω_φ)$, as well as the corresponding Periapsis and Lense-Thirring precession frequencies. The spin parameter primarily determines the spinning properties of spacetime, while the quantum correction and charge particle introduce additional biases to Kerr-Newman geometry, particularly near the event horizon. These corrections alter the position and stability of the circular orbits, and the spectrum effectively associates with particle's motion. The results provide a comprehensive description of the orbital dynamics of spinning EMd QCBH through future high-precision astrophysical observations.

Recent DESI results exhibit preference for a Null Energy Condition violating dynamical dark energy, with early phantom behaviour. We explore an alternative interpretation in which this preference arises from unconventional dark matter dynamics rather than from dynamical dark energy. We propose a dynamical dark matter (DDM) model, with a non-zero equation of state (EoS) that smoothly interpolates between early time and late time asymptotes across a transition scale factor $a_t$, and study its consequences against cosmological datasets including CMB, DESI DR2 BAO, SNeIa (PantheonPlus, Union3, and DESY5) and growth rate data. We find the early time EoS to be consistent with zero, while the present day value is negative at a significance ranging from $0.42σ$ to $3.02σ$ depending on the dataset combination. The strongest preference occurs from the combination of CMB, DESI, and DESY5 giving the present day EoS to be $-0.060^{+0.013}_{-0.028}$ and $a_t = 0.41^{+0.088}_{-0.13}$ at 68\% CL. This preference for a non-zero, late time DM EoS persists when growth rate data are included and across all three SNeIa compilations considered, while the matter density $Ω_m$ mildly shifts to higher values relative to $Λ$CDM. The model also predicts a lower $σ_8$ and $S_8$ than $Λ$CDM, in better agreement with weak-lensing data, while $H_0$ remains unchanged and in tension with local distance-ladder measurements. The DDM model is preferred over $Λ$CDM ($Δχ^2_{\rm MAP} = -14.093$, $Δ{\rm DIC} = -7.838$ for Planck+DESI+DESY5) but disfavored relative to the CPL parameterization of DE ($Δχ^2_{\rm MAP} = 6.755$, $Δ{\rm DIC} = 7.966$). This preference is consistent among other combination of datasets as well.

A growing number of phenomena and theoretical problems indicate that quantum gravity theory is necessary. In this paper, we investigate the evaporation of covariant quantum BHs for particles with different spins and compare the results with the Schwarzschild case. Our results show that the Hawking radiation and mass loss rate of covariant quantum BHs differ from those of Schwarzschild BHs and they depend on the spins of the emitted particles. Therefore, these results suggest that it may be insufficient to consider only BH evaporation in the massless scalar field case and may provide a possible way to test loop quantum gravity in the future.

Double red giant stars are very important for studies of the stability of mass transfer, common-envelope evolution, and the formation of double white dwarfs with short orbital periods. However, no double red giant system undergoing mass transfer has yet been found. We present the discovery of a close Algol-type binary system composed of two red giant stars. This is the first known semi-detached system observed during the very short phase when the accretor has expanded into a red giant just before entering the common envelope phase. The $H_α$ line suggests that the system has recently lost some material, which is now moving toward us. We present a consistent analysis of all the available spectroscopic and photometric observations of this system, constraining its orbital parameters and the fundamental properties of the components. Our findings are supported by a binary evolution model that successfully reproduces the currently observed parameters. The model suggests that the system will eventually merge into a single star.

Primordial magnetic fields (PMFs) provide a well-motivated source of cosmological perturbations through their anisotropic stress and may leave observable imprints in both the scalar and tensor sectors. In this work, we study scalar metric perturbations sourced by a PMF evolving through a finite post-inflationary epoch characterized by a constant equation-of-state (EoS) parameter $w$. Working in a gauge-invariant framework, we derive the sourced evolution equation for the Bardeen gravitational potential in a general constant-$w$ background. Our analysis indicates faster (slower) growth of PMF-sourced scalar perturbations on superhorizon scales for a stiffer $w>1/3$ (softer $w<1/3$) background compared to the marginal radiation-dominated scenario with $w=1/3$. We then derive the scalar-induced gravitational wave (SIGW) background induced at second order by the magnetically generated scalar perturbations. Our analysis indicates a narrow viable parametric region for a kination-like post-inflationary era and a nearly scale-invariant PMF that may give rise to a PMF-sourced SIGW signal dominating over the direct PMF-generated tensor background. Simultaneously, the induced tensor modes are found to be perturbatively small compared to the magnetic sector, which, in turn, exerts negligible backreaction on the background. Interestingly enough, for an inflationary scale $H_{\rm inf}\sim10^3-10^9$ GeV and duration of reheating in $e$-folds $N_{\rm reh}\sim8-10$, the dominant PMF-sourced SIGW signal may be detectable by next-generation terrestrial and space-based interferometric GW detectors spanning the mHz$-$kHz frequency range.

We compute the instantaneous angular momentum flux from nonspinning black hole binaries in Einstein-Maxwell-dilaton theory for generic noncircular orbits up to relative first post-Newtonian order. Working in the Einstein frame and using the direct integration of the relaxed field equations approach, we construct the scalar, electromagnetic, and tensor gravitational fields in the wave zone and express the required source multipole moments in the center-of-mass frame. In addition to the leading $1/R$ radiative fields, we retain the next-to-leading $1/R^2$ terms in the wave-zone fields, which describe finite-distance corrections and do not contribute to the flux at null infinity. We obtain separately the scalar, electromagnetic, and tensor contributions to the angular momentum flux. The scalar and electromagnetic channels begin with dipole radiation, while the tensor channel begins at quadrupole order. Our results recover the quasicircular balance relation, the known general-relativistic limit, and the expected scalar- and electromagnetic-dipole suppression limits. Together with the previously known energy flux, these results provide the dissipative information required for future studies of orbit-averaged eccentric evolution and waveform phasing in Einstein-Maxwell-dilaton theory.

High-performance computing is increasingly dominated by hardware acceleration using Graphics Processing Units (GPUs). To take advantage of this long-term trend, we implement a major overhaul of the parallelisation approach in the relativistic particle-mesh N-body code gevolution. The new version 2.0 of gevolution employs three layers of parallelisation: MPI for scalability on a distributed memory system, shared memory parallelisation on each MPI rank using OpenMP, and offloading all compute intensive tasks to GPUs using CUDA. The code also includes many new features that have been developed over the past years. We provide an overview of the code structure and show key performance benchmarks. The public release of gevolution 2.0 can be found at https://github.com/gevolution-code/gevolution-2.0. We also release a GPU-ready version of the LATfield2 library which provides the parallelisation backend, available at https://github.com/gevolution-code/LATfield2.

We study the cosmological core-halo relation in vector dark matter using three-component Schrödinger-Poisson simulations. Starting from cosmological vector-field initial conditions, which due to the evolution of the vector field during inflation are enhanced on small scales, we find that nonlinear evolution begins almost immediately following matter-radiation equality and produces compact self-gravitating Proca-star condensates at the centers of halos. After confirming the central condensates through their radial density profiles, we find the empirical relation \(\widetilde M_\star\propto \widetilde M_{\rm h}^{0.6403}\) between the Proca star mass and halo mass, although interestingly we find that we are only able to confirm Proca stars in $\mathcal{O}(10\%)$ of halos. This serves as important input for future studies of the abundance and merger rates of Proca stars in models of vector dark matter. We also examine the vector-field structure of the objects through global longitudinal and transverse polarization fractions and local spin density inside halos, which increases over cosmic time.

Exploring how matter cycles through the Galaxy, from the birth of stars in dense interstellar clouds to the ejection of matter and energy during a star's final stages-requires a multi-faceted approach. Radio observations are essential to reveal the intricate interactions at play. By studying our Galaxy in detail, we can use it as a model to better understand these processes in galaxies as well. The Square Kilometre Array offers unique capabilities in wide-field, high-sensitivity, high-resolution spectroscopy and precise astrometry to revolutionise Galactic astrophysics. This chapter presents an overview of the science addressed in detail in chapters pertinent to the SKAO Science Working Group \textit{Our Galaxy}.

Wide stellar binaries are sensitive dynamical probes of dark matter substructure. We generalize existing point-mass disruption constraints to spherical extended dark objects and apply the resulting framework to the population of wide-binary candidates identified in the ultra-faint dwarf galaxy Bootes I. We derive 95% confidence limits on the dark matter fraction as a function of the perturber mass, the radius enclosing 90% of the mass, and the density profile, considering Navarro-Frenk-White, r^-3/2, and r^-9/4 models. Sufficiently compact objects converge to the point-mass limit, whereas finite size suppresses binary disruption when the perturber becomes comparable to the relevant encounter and binary scales. The constraints remain sensitive to dark matter fractions well below unity over a broad region of the mass-radius plane and extend to objects substantially larger than those accessible to conventional microlensing searches. As an application, we map the r^-3/2 results onto ultracompact minihalos and derive an illustrative constraint on the primordial curvature power spectrum. Ultra-faint-dwarf wide binaries therefore offer a purely gravitational probe of the abundance of extended dark matter objects and of the small-scale primordial curvature power spectrum.

A reliable computation of the bubble wall velocity during a cosmological phase transition requires an adequate modeling of the non-equilibrium dynamics in the vicinity of this expanding bubble. This task can be made computationally faster by imposing an \emph{Ansatz} on the shape of the non-equilibrium particle distribution function, thus simplifying the collision terms and making the Boltzmann equation solvable in terms of some out-of-equilibrium fluctuations. Two different \emph{Ansätze} have prevailed in the recent literature: the so-called fluid \emph{Ansatz} and an expansion in a basis of Chebyshev polynomials, consolidated in the public code \texttt{WallGo}. In this work we show that the two approaches yield essentially the same wall velocity in the regime of reasonably mild phase transitions, $α\lesssim 0.01$. Interestingly, the agreement is excellent when only top-quark annihilation is considered, but a noticeable discrepancy appears once scattering processes are included. We also investigate the limitations of linearizing the Boltzmann equation when the fluid \emph{Ansatz} is applied to stronger phase transitions, showing that non-linear contributions induce significant shifts in the predicted terminal velocity as $α\to 1$, even though the non-linear contribution to the wall pressure remain quantitatively small compared to the equilibrium and linearized non-equilibrium parts. We discuss possible consequences of this result for both \emph{Ansätze}, while also highlighting the possible limitations of the WKB approach itself when applied to the regime of strong transitions. Since strong phase transitions are precisely the primary targets for future gravitational waves observatories, our study emphasizes that not only higher precision computations of $v_w$ in the semi-classical approach are required, but a treatment beyond the WKB approximation may be needed.

Massive cosmological correlators exhibit a rich hypergeometric structure already at tree level, reflecting the distorted propagation of particles in de Sitter spacetime. In this paper, we reveal that this apparent complexity conceals a remarkably simple underlying mathematical structure. Using the spectral representation, we compute arbitrary tree-level correlators of scalar fields with generic masses and show that they are constructed from fundamental building blocks belonging to the family of Lauricella generalised hypergeometric functions, glued together by spectral integrals. We develop a spectral gluing algorithm that evaluates these integrals through elementary graph combinatorics, yielding explicit series representations that resum the dependence on internal energies away from soft limits. This algorithm naturally generates solutions to the differential equations satisfied by massive correlators as expansions in the corresponding eigenfunctions. Acting with a set of graph annihilators, we uncover a new class of magical identities among generalised hypergeometric functions, revealing an unexpected simplification: once the dynamical propagators are stripped away, the remaining hypergeometric kinematic dependence collapses to rational functions. Our results expose a hidden simplicity in the rigid hypergeometric analytic structure dictated by graph combinatorics, and hint at an intrinsic geometric principle from which properties of massive correlators naturally emerge.

Constraining local primordial non-Gaussianity through scale-dependent bias requires an accurate determination of the PNG bias parameter $b_φ$, which is completely degenerate with $f_\mathrm{NL}$ in the galaxy power spectrum. The standard universality relation predicts $b_φ$ from the linear bias, but breaks down for realistic galaxy populations due to PNG assembly bias. The time evolution estimator offers a promising alternative, inferring $b_φ$ from the redshift evolution of tracer number counts. We test it against the separate universe estimator using the IllustrisTNG and SIMBA suites of CAMELS. Our main contribution is to demonstrate that the estimator can be extended to realistic observer-frame photometric selections, provided correction factors derivable from the galaxy SEDs are applied to account for luminosity-distance dimming and $K$-corrections. We further show that physical evolution not directly tracing the growth of matter fluctuations can also be corrected under assumptions about the star formation history: magnitude cuts require correcting for passive stellar fading, after which agreement with the separate universe estimator is restored. This establishes the robustness of the time evolution estimator of PNG bias with realistic survey selections, paving the way to future applications with real data.

Future interferometers may detect a gravitational-wave (GW) signal from a cosmological first-order phase transition. Reconstructing the underlying particle-physics model from such a signal requires theoretical control over the map from microphysics to the spectrum. For classically scale-invariant extensions of the Standard Model, which generically predict strongly supercooled transitions and strong GW signals, this map depends sensitively on the treatment of quantum and thermal corrections to the nucleation rate. Taking the classically conformal ${\rm U}(1)_X$ model as representative of this class, we scan its parameter space and compare two resummation schemes. The first is a high-temperature effective field theory, matched at two-loop level and including next-to-leading-order corrections to the bounce action, with the nucleation-rate prefactor given by the full one-loop functional determinants. The second is a commonly employed daisy-resummed effective potential, with the prefactor estimated on dimensional grounds. Reconstructing the fundamental model parameters through a Fisher-matrix analysis of injected GW signals at LISA, we find that the daisy-resummation scheme is strongly disfavored, as its theoretical error dominates over the reconstruction uncertainty.

In this paper, we investigate the shadow of the generalized Vaidya black hole described by the Husain solution. Assuming a barotropic equation of state for the additional matter sector, we first identify the self-similar class admitting a homothetic Killing vector and transform the geometry to a conformally static form. This allows us to derive the effective potential for null geodesics in the conformally static representative, determine the corresponding homothetic photon surface, and then transform the result back to the original generalized Vaidya coordinates. We show that, within the weak-energy-condition branches, the parameter $α$ controlling the equation of state determines whether the shadow is enlarged or reduced relative to the Vaidya case: the branch $0\leqα<1/2$ increases the photon-sphere and shadow radii, whereas the branch $1/2<α\leq 1$ decreases them. We then turn to the genuinely time-dependent Husain metric and derive quasistatic slow-evolution equations for the instantaneous photon surface and critical impact parameter. In particular, we show that the growth or contraction of the shadow is governed by the local effective influx in the photon region rather than by the signs of $\dot M$ and $\dot N$ separately. For the charge-like branch, this yields a transparent criterion for when accretion enlarges the shadow and when rapid charge growth instead causes it to shrink.

We build an ultralight dark matter model with repulsive self-interactions, starting from a 5D action that only includes a $U(1)$ gauge field and several massive and charged free fermions. The dark matter scalar field corresponds to the fifth component of the gauge field, after compactification to 4D. Although the single-fermion case only leads to attractive self-interactions, two fermions can already give rise to repulsive self-interactions, thanks to a Scherk-Schwarz twist around the fifth dimension. We check the observational and theoretical self-consistency of this scenario, from the inflation stage to the current time. We find that large ranges of model parameters are allowed. For scalar masses below $10^{-12}$ eV the self-interactions are negligible and the model behaves as fuzzy dark matter. For higher masses the repulsive self-interactions govern the formation of solitons of astrophysical size. Larger sizes require a large initial misalignment or small masses in the fuzzy dark matter regime.

The intracluster light (ICL) has been proposed as a luminous tracer of the dark matter (DM) halo shape in galaxy clusters, with recent observational studies and $Λ$CDM simulations suggesting close alignment between the two components. We test whether this correspondence holds in a self-interacting dark matter (SIDM) cosmology using matched CDM and SIDM ($σ_\mathrm{DM}/m = 0.1$ and $0.5\;\mathrm{cm^{2}\,g^{-1}}$) hydrodynamical cluster simulations from TheThreeHundred project. We measure 3D axis ratios and major-axis orientations as a function of radius for both the DM halo and the stellar (BCG+ICL) distribution, and quantify the spatial correspondence using the Weighted Overlap Coefficient. In the SIDM simulations, DM haloes become significantly rounder in their inner regions with increasing DM cross-section, while the stellar distribution retains the same elongated morphology in both CDM and SIDM. This decoupling is robust across two baryonic physics models and two resolution levels, and is detected at $>5σ$ for $σ_\mathrm{DM}/m = 0.5\;\mathrm{cm^{2}\,g^{-1}}$. It arises because the stellar shape is governed by the accretion geometry of tidally stripped satellites along the large-scale structure, which is identical in CDM and SIDM, whereas, in SIDM, halo shapes respond additionally to local self-interactions. The ICL-DM shape correspondence demonstrated in CDM therefore breaks down in an SIDM universe.

The Solar, Heliospheric and Ionospheric (SHI) Physics Science Working Group of the Square Kilometre Array Observatory (SKAO) addresses the full chain of plasma processes linking the solar corona to the terrestrial environment. This overview chapter synthesises 16 topical contributions to Advancing Astrophysics with the SKA-II}, spanning the quiet and active solar atmosphere, eruptive phenomena, heliospheric turbulence and solar-wind diagnostics, ionospheric science, stellar-solar connections, and the observational frameworks required to deliver these science goals. The primary focus is on the capabilities of Array Assembly 4 (AA4), the design baseline for both SKA-LOW (50-350,MHz) and SKA-MID (0.35-15.4,GHz), which together provide continuous spectral coverage, sub-arcsecond angular resolution, full-Stokes polarimetry, and sensitivity gains of an order of magnitude over existing facilities. From resolving fine-scale coronal heating events to mapping coronal mass ejection magnetic fields and characterising multi-scale heliospheric turbulence, SKAO will deliver transformative advances in solar and space-weather science. We frame these contributions as a single end-to-end Sun-to-Earth system. We identify cross-cutting themes and gaps not fully addressed by individual chapters and outline the staged roadmap from early operations through to the full AA4 capability.

In general relativity, cosmic acceleration is timelike defocusing of the comoving congruence and requires a negative total active gravitational mass density, $\mathcal{M}_{\rm tot}=ρ_{\rm tot}+3p_{\rm tot}<0$. The criterion $w\equiv p/ρ<-1/3$ diagnoses sector repulsion only for $ρ>0$: the inequality reverses for $ρ<0$, and ratio variables are ill-defined at $ρ=0$ even when the stress-energy tensor is finite. For sign-changing effective dark energy (DE), as in $Λ_{\rm s}$CDM-type histories, we instead use the signed density $ρ_{\rm de}$ and two branch-independent combinations. The regular null energy condition (NEC) boundary $\mathcal{I}_{\rm de}=ρ_{\rm de}+p_{\rm de}=0$ replaces the phantom divide $w_{\rm de}=-1$, while $\mathcal{M}_{\rm de}=ρ_{\rm de}+3p_{\rm de}<0$ governs sector-level Raychaudhuri repulsion. For a separately conserved DE sector with a smooth negative-to-positive density crossing of finite odd order $n$ at $z_\dagger$, we prove that $\mathcal{I}_{\rm de}$ and $\mathcal{M}_{\rm de}$ are negative in a punctured neighborhood and non-positive at the crossing, while $w_{\rm de}$ develops a kinematic pole with universal residue $n(1+z_\dagger)/3$. If $\mathcal{M}_{\rm de}>0$ at some sufficiently high redshift, continuity requires at least one repulsion boundary $z_{\rm rep}>z_\dagger$: the sector is already repulsive while $ρ_{\rm de}<0$. We derive the exact range of $ρ_{\rm de}'(z_\dagger)$ for acceleration at the crossing with the total NEC satisfied. Under the stated single-impulse and stationary-point assumptions, the deceleration parameter has one or three sign-changing zeros. A smooth $Λ_{\rm s}$CDM profile, an exponential infrared $f(T)$ model, and the minimal phantom brane illustrate the results. These results motivate organizing late-time inference around $(ρ,p,\mathcal{I},\mathcal{M})$ rather than around $w$ alone.

The ultrahigh-sensitivity polarization imaging of the cosmic microwave background (CMB) is a treasure trove for new physics. Searching for a predicted polarization angle rotation known as the cosmic birefringence effect is a pursuit of the next-generation CMB experiments, and this new polarization signal would be important for testing fundamental physics theories. However, such a delicate rotation effect can be confused by different mechanisms including polarized Galactic foregrounds and instrumental calibration effects. Also, the cosmological rotation effects may arise from both the background evolution and spatial fluctuations of axion-like particles, leaving complicated imprints on CMB polarization fluctuations. In this work, we establish a comprehensive modeling of cosmological rotation effects, the instrumental miscalibration of polarization angles and complex foreground polarization with higher-order fluctuations, and perform a Bayesian analysis to infer both isotropic and anisotropic rotation angles simultaneously for foreground-mitigated polarization datasets from submillimeter to millimeter wavelengths. We find that such a Bayesian formalism can effectively extract the rotation effects for the next-generation CMB experiments and investigate the impact of the known sources including Galactic foregrounds and lensing on the inferred cosmological signals. The method in this work is well suited for future polarization data analyses, and the inferred rotation effects may shed light on the nature of the axion-like particles.

Planet formation begins with dust growth and planetesimal formation within protoplanetary disks surrounding young stars. To understand these processes, it is essential to estimate dust grain sizes from disk observations. In this study, we develop a new method to constrain grain size based on the estimation of cooling timescales. Our approach applies to transitional disks that possess an inclined inner disk casting shadows on the outer disk, whose temperature variations serve as a tracer of dust properties. By constructing a three-dimensional model of the disk surface using near-infrared scattering light images and comparing it with submillimeter dust continuum maps, we estimate the spatial offset between the irradiated and shadowed regions to derive the cooling timescale. We then build an analytic model that calculates the cooling timescale at the dust thermal emission height with an assumed turbulent diffusion intensity to infer the dust surface density and dust grain size. Applying this method to the protoplanetary disk around HD~142527, we find that the disk's northern shadowed region cools on a timescale of a few percent of the orbital period and that the maximum grain size consistent with the observations is approximately 0.1-1 mm. We also find that the conditions required for the vertical shear instability, which needs a short cooling timescale, are satisfied, allowing turbulence with an intensity consistent with near-infrared observations. This study demonstrates that estimating cooling timescales is an effective tool for constraining dust grain size. Our approach can be generally applied to other transition disks with inner-disk-induced shadows.

Stellar age is a fundamental quantity for Galactic archaeology, but reliable age estimation for large stellar samples remains challenging. In this work, we develop an uncertainty aware NGBoost framework for stellar age estimation using Gaia XP-derived atmospheric parameters and chemical abundances. Different from the standard NGBoost model, we modify the loss function by incorporating the uncertainties of the training age labels. We further use a Monte Carlo strategy to quantify the influence of input-feature uncertainties on the predicted ages. The resulting model provides age estimates together with uncertainty estimates. Applying this framework to Gaia XP stars, we construct a stellar age catalog containing 15,175,107 stars.

AI Phe is an eclipsing binary star with an orbital period of 24.6 days for which the surface gravity and effective temperature are known from direct measurements to very high precision and accuracy. We have obtained high-quality spectroscopy of the K0IV star during the total eclipse of the F7V companion, and also obtained spectra with a very high signal-to-noise ratio for this star and its F7V companion using the spectral disentangling technique. We have used these spectra to measure the abundances of iron and magnesium for both stars. We compare the values of [Fe/H] and [Mg/H] for the F7V star and the K0IV star to stars in M67, an open cluster of similar age and metallicity to AI Phe. We find that our [Fe/H] and [Mg/H] measurements clearly show the signature of elemental diffusion in the F7V star. This suggests that AI Phe can be used to test models of single stars that include diffusion and mixing of elements.

A compact form of the stress tensor $\langle T^{ab}T^{cd}\rangle$ correlation function of a quantum scalar field of arbitrary mass $m$ and curvature coupling $ξ$ is presented, with particular emphasis on decoupling in the $m\to\infty$ limit and infrared origin of the conformal anomaly as $m \to 0$. The Ward Identity (WI) of covariant conservation is verified for the full second order metric variation of the one-loop effective action, of which $\langle T^{ab}T^{cd}\rangle$ is part, including local contact terms. This WI is satisfied by two tensors, one spin-2 (traceless) and the second spin-0 (traceful), each multiplied by a Lorentz invariant form factor computed in $n$-dimensional regularization. Minimal subtraction of pole terms at $n =4$ produces form factors that fail to satisfy decoupling for general $ξ\neq 1/6$, but can be simply amended to do so. Particular interest attaches to the $ξ=1/6$ case, in which the spin-0 form factor at $n =4$ is completely finite, requires no UV regularization or subtractions, and satisfies decoupling directly. Its imaginary part defines a spectral function that obeys a UV finite sum rule for any $m$, while its real part unambiguously determines the $\square R$ in the massless limit. Its $m^2$ dependence describes a Wilsonian renormalization group flow to an effective field theory limit at large distances. At $m= 0$ the spin-0 spectral function becomes a Dirac $δ$-function at zero energy, demonstrating the existence of a massless scalar Goldstone collective excitation due to the conformal anomaly in the low energy effective theory of gravity.

We re-investigate the open cluster UBC 63 using the Gaia DR3 data and show that, rather than being a single cluster as previously classified, it is a compelling candidate for a double cluster undergoing an unbound flyby interaction. A GMM decomposition performed in the 5D astrometric space reveals the two statistically distinct components, namely UBC 63A (98 members, Age = 21 $\pm$ 4 Myr) and UBC 63B (148 members, Age = 562 $\pm$ 43 Myr). A significant age difference of $Δ\mathrm{Age} = 541 \pm 43$ Myr between the clusters, rules out coeval formation. Their 3D separation of $60 \pm 29$ pc at the birth-epoch of the younger cluster, indicates that the clusters might have originated from the same molecular cloud complex. At present, the system exhibits a 3D separation of $26 \pm 8$ pc, with a relative velocity of $3.60 \pm 1.80$ km s$^{-1}$. Orbital integrations and \textit{N}-body simulations of the pair suggest that the systems had a close encounter, reaching a separation of $7 \pm 2$ pc only $\sim$~6 Myr ago and predict a rapid divergence to a separation of $491 \pm 213$ pc within the next $\sim$100 Myr. The low escape velocity ($V_{\rm esc} = 0.51 \pm 0.12$ km s$^{-1}$) of the system compared to the relative 3D velocity indicates that they are gravitationally unbound. Their low tidal factors, elongated structures and populations extending beyond the Jacobi radii may reflect a strong transient tidal interaction between the clusters.