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

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.

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.

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.

Primordial black holes in the planetary-mass range have attracted renewed interest; however, the search for gravitational waves from such binaries remains challenging due to their long-lived nature. In this work, we present, define, and validate a fully operational search pipeline developed to detect planetary-mass binaries during their inspiral phase. We use the Viterbi algorithm, a dynamic programming technique that recovers the most likely track based on a Hidden Markov Model. To enhance its performance, we introduce a novel time-frequency representation of the data and a candidate isolation procedure that separates signals from background noise. The evaluation of candidates is carried out using the two detection statistics, $n_σ$ and NMSE, which quantify the power significance and the consistency with the expected binary evolution. We then validate the search method using O3 LIGO Hanford data with a population of injected signals. The pipeline is able to recover most of the signals with a fixed false-alarm ratio of $3\%$, covering Galactic scales across most of the parameter space and reaching luminosity distances $ \gtrsim 100$ kpc in the most sensitive region. For each candidate, we also obtain an estimate of the system's chirp mass, whose accuracy remains high throughout the detectable range, enabling a rapid characterization of the system upon detection.

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.

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.

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.

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.

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.

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.

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.

We present a comprehensive multiband photometric and spectroscopic study of the G-type binary CW Cas whose parameters have not been well determined. Our double-lined spectroscopic radial velocity curve of this system yields a reliable mass ratio of $q = 1.88(9)$. By combining $BVR_{c}I_{c}$ bands, TESS light curves and radial velocity curves, we found that CW Cas is a W-subtype shallow contact binary with a fill-out factor of 15\%. The components have masses of $0.98(6)M_{\odot}$ and $0.52(4)M_{\odot}$, separated by $2.25(5)R_{\odot}$. A notable asymmetry in the maxima of the light curves was detected and explained by a dark spot located on the surface of at least one component. Comparison of light curves from different years revealed that these dark spot activities exhibit cyclic variations with an approximate period of 1250 days. Orbital period analysis via O-C diagram spanning 125 years shows a long-term decrease superimposed with periodic oscillation caused by the light-travel-time effect (LTTE) due to a third body. This tertiary component has an orbital period of $P_{3} = 99.4(6)$ years and a minimal mass of $M_{3}=0.91(1)M_{\odot}$. The absence of detectable signatures for this massive object in either spectroscopic or photometric datasets implies it must be a compact object such as a white dwarf or neutron star. Furthermore, a visual companion was identified based on Gaia DR3 astrometric data, suggesting that CW Cas is part of a hierarchical quadruple system. As such, CW Cas represents a valuable laboratory for probing 2+1+1 hierarchical multiple system hosting compact object.

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.

Cosmic strings are not yet confirmed, topological defects formed in the early Universe. They can bend light or gravitational waves, which causes an effect similar to the gravitational lensing. Our goal is to check whether cosmic string could be detected as a lenses of gravitational waves by the Einstein Telescope (ET). To do that the apparatus of wave optics had been applied. Firstly we explored the amplification factor strength and behaviour. Next the wave effects in SNRs and waveforms was examined for different inclinations of the source. Lastly we estimated the string tensions based on eight mergers from the \textsc{StarTrack} simulation, using the Bayesian Inference methods. The wave effects were easily to see in waveforms, SNR and characteristic strains. Also most of the events could be detected, based on their Signal to Noise Ratio values. Almost whole characteristic strain lies in the range of the ET. When it comes to mock observations, we had got estimated vale of a logarithm of the CS tension equal to $\bar{\log{Gμ}}=-9.78^{+0.44}_{-0.49} $, which were consistent with injected value equal to $-10$. At the end we conclude that CSs could be detected by the ET.

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.

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.

Primordial black holes (PBHs) evaporate through Hawking radiation, emitting high-energy photons and ionising their surrounding environment. As contributors to the density of dark matter, Omega_C, if 10^-18 solar mass PBHs are present in the Universe they delay recombination and move the surface of last scattering of the cosmic microwave background (CMB). We perform recombination simulations using the software Recfast, and calculate the PBH fraction of dark matter required to resolve the Hubble tension. We find that nominally a cosmic PBH energy density of Omega_PBH ~ 10^-3 Omega_C would cause an 8.9% increase in the value of H_0, enough to entirely reduce the tension between the early- and late-time observations. This PBH fraction is modified by Gray Body Factors affecting Hawking radiation. Furthermore, also fitting the CMB leaves the Hubble tension not fully relieved with our present PBH prescription. Until relevant non-gravitational properties of the dominant dark matter species are ruled out, we suggest that the hypothesis that the ionisation history of the universe matches the thermal history of the standard LCDM cosmology is too precarious to hang the expansion rate on, and that it is better to measure H_0 locally at z <~ 1.

Massive bosonic fields can form self-gravitating solitonic structures, which for vector fields are known as Proca stars. For ultralight fields, these structures can describe the cores of dark matter haloes surrounding the supermassive black holes at the center of galaxies. It has been argued that future gravitational-wave detectors might be able to probe the properties of dark matter structures. However, most of the analyses considering ultralight dark matter have focused on massive scalar fields. In this work, we study how Proca stars respond to a perturbing small object inspiralling in their interior, in the Newtonian limit. We consider both Newtonian spherically symmetric Proca stars and the non-spherically symmetric ground-state solutions. We compute the total energy lost by the orbiting object and compare it to the case where the bosonic star is composed of a scalar field. Our results show that, at the Newtonian level, the energy lost by the object in a Proca star ground-state is similar to that in its scalar field counterpart, with relative differences of at most $\sim 20\%$, while the maximum orbital energy loss rate in the ground-state configuration can be one to two orders of magnitude larger than in the spherically symmetric Proca star, when a central parasitic black hole is present. Our work motivates the need to study extreme mass-ratio inspirals into Proca stars using a fully relativistic setup, where differences between the Proca and scalar boson star ground states are expected to become more significant.

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.

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.

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.

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.

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}.

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.

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

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.

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.

As they orbit within the Galaxy, stars swim through a vast population of interstellar objects (ISOs). In this paper, we use N-body simulations to show that a fraction of ISOs within $\sim$1 pc of the Sun (its tidal radius) may be captured during the flyby of another star -- a mechanism that requires no planets. Capture is most efficient when the impulse imparted by the flyby is comparable to the escape speed at the widest stable orbit, which is roughly 0.1 km/s for the Sun. ISO capture is dominated by the few highest-impulse stellar flybys, typically involving relatively slow encounters with massive stars. Most ISOs are captured in the outer parts of the Oort cloud, with semimajor axes greater than $\sim$50,000 au. Using Monte Carlo simulations, we show that the Sun underwent only a small number of ISO-capturing flybys in its history (median [mean] of 1 [1.7]). Using the {Ō}tautahi-Oxford population model, we estimate that a few times $\sim$$10^{4}$ `Oumuamua-sized ISOs were likely captured by the Sun. This only represents a $\sim$$10^{-8\pm1}$ contribution to the total Oort cloud population, yet it contains roughly as many present-day captured ISOs as Jupiter-assisted capture provides. Given that flybys are unavoidable in the Galactic field, most stars should, at any given time, host have sparse Oort clouds populated with ISOs captured during stellar flybys. Massive stars are both the main drivers of capture when they fly by a given star, and more efficient at capturing ISOs around themselves than low-mass stars.

We constrain the class of modified cosmologies derived in Paper I [1] from a generalized mass-to-horizon relation (MHR) that enforces thermodynamic consistency between the Cai-Kim horizon temperature and generalized horizon entropies. The modified Friedmann equations depend on an entropy exponent $m$, an MHR coupling parameter $γ$, and an entanglement-correction amplitude $f_B$, with standard $Λ$CDM recovered in the appropriate limit. Using Pantheon$+$/SH0ES Type~Ia supernovae, cosmic chronometers, DESI DR2 baryon acoustic oscillations, and Planck 2018 CMB distance priors, we constrain eight physically motivated sub-cases via Markov chain Monte Carlo and compare models through the Bayesian log-evidence. The entropy exponent is tightly bounded, $|m-1|\lesssim O(10^{-4})$ when the MHR coupling is fixed ($γ=1$), relaxing to $O(10^{-3})$ along the $m$-$γ$ degeneracy, excluding any macroscopically significant departure from standard horizon thermodynamics. Freeing the MHR coupling parameter or the entanglement amplitude raises the inferred Hubble constant to $h\simeq0.70$-$0.71$, reducing the CMB-SH0ES tension from ${\sim}4σ$ to ${\sim}1.2$-$2.6σ$, but no scenario fully resolves it within a flat universe. The Bayesian log-evidence nevertheless disfavors every extension relative to $Λ$CDM in all dataset combinations ($-16\lesssimΔ\ln Z \lesssim-1$): the improved fits obtained when the SH0ES calibration is included reflect an absorption of the Hubble tension by the additional parameters rather than genuine evidence for modified horizon entropy.

We analyze the vacuum dynamics of Gravity from Entropy, including its algebraically constrained $G$-field formulation. Evaluating the curvature traces over zero-, one-, and two-form sectors, we show that the complete Minkowski Hessian is exactly that of the quadratic-gravity action $A R+B R_{μν}R^{μν}$, with $A=3β/\ell_{\rm P}^{4}$ and $B=5β^{2}/(2\ell_{\rm P}^{4})$. For diagonalizable curvature blocks, the same action reduces to a sum over eigenvalue logarithms and reproduces these coefficients exactly. A strict diagonal-curvature restriction on the perturbations is instead only a reduced subsector and excludes non-diagonalizable type-N wave curvatures. Linearizing the $G$-field equations and subsequently imposing the algebraic vacuum constraint reproduces the same reduced metric equation and covariant Minkowski Hessian. The spectrum contains the massless graviton, a scalar with $m_{0}^{2}=3/(5β)$, and an opposite-residue spin-2 branch with $m_{2}^{2}=-6/(5β)=-2m_{0}^{2}$. For the foundational choice $β>0$, conventional Einstein normalization therefore implies a tachyonic spin-2 instability. We also show that every four-dimensional Ricci-flat metric solves the local bulk equations through quadratic curvature order, while square-zero Ricci-flat pp-waves are exact local vacuum solutions of the analytic metric-only logarithmic branch. On the isolated massless transverse-traceless eigenspace, the quadratic translation current has the standard general-relativistic normalization.

The Nancy Grace Roman Space Telescope's Wide Field Instrument (WFI) will observe hundreds of thousands of bright stars across its Core Community Surveys, particularly in the dense stellar fields of the Galactic Bulge Time Domain Survey (GBTDS). Sources brighter than ~17th magnitude will saturate WFI detector pixels in typical survey exposures, with the brightest stars deeply saturating large pixel regions and potentially producing persistence signals that may impact subsequent observations. Prior detector characterization did not explore the regime of deep point source saturation. To address this gap, we conducted a bright star saturation test during WFI's second Thermal Vacuum test campaign (TVAC2) at BAE Space & Mission Systems in Boulder, CO. Using the Stimulus of Ray Cones (SORC) telescope simulator, we projected nine in-focus point sources through the F146 filter onto two Sensor Chip Assemblies (SCAs), with fluxes tuned to approximate stellar magnitudes ranging from ~4 to ~18 in ~170 s exposures. We present analyses of the saturation response and persistence properties of these detectors. We find that the saturated region of a ~4 mag source grows to ~150 pixels in diameter after ~170 s of illumination, compared to ~15 pixels for a ~12 mag source. Pixels adjacent to the expanding saturation front exhibit pronounced non-linear behavior consistent with charge leakage from saturated neighbors. For persistence, we find that the median signal in the first post-illumination dark exposure is broadly consistent across source magnitudes spanning ~4 to ~17, and that persistence decays to detector background levels ($\lesssim$0.05 e$^{-}$ s$^{-1}$) within approximately 20 minutes, consistent with flat field persistence measurements from the same TVAC2 campaign. These pre-flight characterization results inform community understanding of WFI detector response to prepare for Roman science.

We investigate the nonradial quasinormal-mode spectrum of static charged Ellis--Bronnikov wormholes in Einstein--Maxwell theory minimally coupled to a phantom scalar field. The background solutions are known in closed form and comprise three classes: subcritical, critical and supercritical, which all approach the extremal Reissner--Nordström geometry at the boundary of their domain of existence. We derive the linear perturbation equations for axial and polar sectors, including the coupled gravitational, electromagnetic and phantom-scalar degrees of freedom, and compute the corresponding spectra by means of a Chebyshev spectral method. The uncharged limit reproduces the known Ellis--Bronnikov spectrum and exhibits the expected electromagnetic isospectrality. For charged configurations we track the axial and polar branches across the three families of solutions and identify the effect of the charge on the damping times and oscillation frequencies. In particular, we find that charge can substantially reduce damping rates as the extremal Reissner--Nordström limit is approached. We also uncover a nonradial polar instability, most clearly visible in the fundamental $l=2$ branch for sufficiently large wormhole masses. This instability is distinct from the familiar radial Ellis--Bronnikov instability and shows that the nonradial sector imposes additional constraints on the dynamical viability of charged wormholes.

We compute, for the first time, the Hawking emission spectrum of massive vector (Proca) fields by spinning Kerr black holes, determining the associated greybody factors and the resulting mass and spin loss functions. We show, in particular, that the scalar (longitudinal) polarization of the Proca field has a spectrum approaching that of a free scalar field in the massless limit (in which it becomes a pure gauge mode), although we find substantial differences for finite mass. The contribution of the two vector (transverse) polarization modes coincides, as expected, with the one obtained by Page for the Maxwell field in the massless limit. The black hole's evaporation rate is dominated by the scalar mode for slowly spinning black holes and by the two vector modes as the black hole approaches extremality. As for other fields, we find that Proca Hawking emission is Boltzmann-suppressed for Hawking temperatures $T_H\lesssim |μ-Ω_H|$, where $μ$ is the field mass and $Ω_H$ is the angular velocity of the black hole's horizon. This implies that highly spinning black holes can efficiently emit massive vector fields at temperatures parametrically below the field's mass. Finally, we also find that superradiant emission is more pronounced for massive vector fields, with a maximum amplification factor of $\simeq 7\%$ (compared to $\simeq 4\%$ for massless photons).

Recent DESI results are in tension with the constant dark energy density predicted by the $Λ$CDM model. If dark energy is associated with the vacuum energy of a scalar field in a metastable state, it will undergo a first-order phase transition through the nucleation of bubbles containing a reduced dark energy density. In this paper, we explore the consequences of this model, where dark energy density varies in both space and time. We model a single bubble spacetime using the Israel junction conditions and derive each of the usual distance measures in this inhomogeneous cosmology. We find that the model predictions of Alcock-Paczynski distortion have features that align strikingly well with the DESI measurements if the dark energy phase transition occurred at roughly a redshift of 1.4 and if the bubble of lower dark energy density has roughly 10$\%$ less dark energy density than the outer cosmology. Despite this feature, we find that the dark energy bubble is heavily constrained by the CMB which excludes the region of parameter space that reproduces the DESI BAO measurements. Still, the peculiar features in the distance measurements of the dark energy bubble cosmology serve as a useful toy model to motivate and inform future work in the currently poorly explored area of spatially varying dynamical dark energy.

The cosmological collider program probes ultra-high-energy physics by searching for subtle oscillatory signatures induced by inflationary particle exchange. Due to the non-linear symmetries usually assumed in inflation, these signals do not appear in isolation; moreover, their amplitudes are bounded by perturbativity and unitarity. To comprehensively probe cosmological collider physics, we must jointly analyze the full multi-field inflationary Lagrangian: in this work, we conduct such a study, probing single, double, and triple scalar field exchange and self-interactions across a wide range of masses (in both the complementary and principal series), mixings, and sound-speeds. Using modern theoretical tools (including the cosmological bootstrap and the cosmological flow), we construct a vast library of tree-level primordial bispectra: combining these with recent measurements of the inflationary shape function from Planck and modern sampling techniques, we perform Bayesian exploration of the eight-parameter multi-field likelihood. Most previous studies of cosmological collider physics assume weak mixing, such that the coupling between the scalar field and the Goldstone mode can be treated perturbatively: in our companion study, we demonstrate that this assumption is incompatible with current datasets except at the smallest masses. By solving for the inflationary bispectra numerically using CosmoFlow, we perform the first analysis of the strongly-mixed collider, demonstrating that the Planck constraints are much tighter than the theoretical bounds, though the oscillatory contributions are heavily suppressed. Across the full multi-field landscape, we find no evidence for new physics with a maximal $χ^2$ improvement of $5.3$.

We investigate how localized quantum corrections to the inflaton potential affect preheating dynamics and the resulting stochastic gravitational wave (GW) spectrum. When these corrections sufficiently suppress the quadratic term of the potential near its minimum, inflaton self-resonance can produce a peaked GW spectrum. On the other hand, we find that a smooth and enhanced spectrum can appear if the quadratic term acquires a negative coefficient. As a concrete realization, we analyze the $α$-attractor T-model with a one-loop Coleman--Weinberg correction induced by a heavy scalar and compute the resulting GW spectra using lattice simulations. The GW signals lie in the ultra-high-frequency regime at frequencies above the kHz range. These results suggest that GW signals from preheating may probe quantum corrections to the inflaton potential, thereby providing indirect information about the underlying UV physics.

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

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.

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 map the observational predictions of cosmological stasis in the inflationary gravitational wave background onto the sensitivity bands of current and planned gravitational wave detectors. Using the closed-form piecewise spectral template derived in the companion paper, we generate detectability maps for four stasis scenarios: canonical, dynamical scalar, vacuum-energy/matter, and vacuum-energy/radiation across the frequency bands probed by NANOGrav, SKA, LISA, DECIGO, BBO, the Einstein Telescope, and Cosmic Explorer. For scenarios in which the spectrum is suppressed, $w_s < 1/3$, the stasis feature is detectable by BBO in the region of $(w_s,ΔN)$ parameter space in which $w_s\gtrsim 0.2$ for tensor-to-scalar ratios close to the Planck upper limit, r = 0.036. For scenarios in which the spectrum is enhanced, $w_s > 1/3$, the stasis feature is detectable by BBO across the entire $(w_s,ΔN)$ parameter space for tensor-to-scalar ratios of $O(0.01)$. We characterize the Standard Model (SM) $g_*$ fine structure of the IGWB, showing that SM phase transitions introduce spectral steps of $\approx 20\%$ (electroweak, at $\sim 2.6\times10^{-6}$~Hz) and $\approx 53\%$ (QCD, at $\sim 3.6\times 10^{-9}$~Hz). For stasis scenarios with end-of-stasis temperatures below the QCD scale these steps fall inside the stasis band and constitute additional spectral features that complement the primary signature. Finally, we model the finite-width end-of-stasis transition phenomenologically, demonstrating that the spectral break at $f_{end}$ is smoothed over a log-frequency window $ΔN_\mathrm{trans}\times 3(1+w_s)/4$, and that the consistency relation $C^2=C^2(α)$ remains testable provided $ΔN_\mathrm{stasis}\gg ΔN_\mathrm{trans}$, a condition easily satisfied for all scenarios of phenomenological interest.

GUEST is a space mission concept whose central objective is the detection of gravitational waves (GWs) in the microhertz band -- a physics-rich frequency window that no other present or planned detector can reach at a significant level. The concept is simple: two dense, passive spheres, covered with cube-corner retroreflectors, deployed in {highly eccentric} Earth orbits ($e \gtrsim 0.7$, period $P \gtrsim 33$ h), tracked continuously by the global network of satellite laser-ranging stations over a minimum observation time of 10 years, with an expected total duration of 30 years. The orbits themselves act as resonant detectors of the oscillating gravitational perturbations, with the microhertz sensitivity emerging from the selected orbital parameters. From the same data stream, GUEST delivers a programme of fundamental and applied science that cuts across particle physics, gravitational-wave astronomy, cosmology, astrophysics, and geodesy: the first coherent search for GWs from supermassive black-hole binaries in the $μ$Hz band, the exploration of primordial GW backgrounds in the unexplored energy-scale gap between pulsar-timing arrays and LISA, a dedicated probe of ultra-light dark matter in a parameter region untouched by any other experiment, a new way to search for ultra-light bosons, order-of-magnitude-improved tests of new gravitational interactions at astronomical ranges, and a step change in the absolute determination of $GM_\oplus$ that underpins the Global Geodetic Observing System and future navigation and Earth-observation missions. This white paper presents the motivation, scientific reach, and mission concept of GUEST.

Probing inelastic dark matter at large mass splittings requires heavy target nuclei, an extended recoil-energy range, and the high-velocity tail of the dark-matter distribution. We exploit these features with the RES-NOVA prototype detector, featuring a PbWO4 cryogenic calorimeter, produced from archaeological Pb and operated at the deep-underground laboratory of Gran Sasso of INFN (Italy), analyzing a 32.4 g day exposure over 2.5 keV - 1 MeV under both the Standard Halo Model (SHM) and a Large Magellanic Cloud (LMC)-motivated velocity distribution. We extend direct-detection constraints beyond the 330 keV reach of established technologies (e.g. Xe-based TPCs), probing splittings up to 510 (780) keV in the SHM (LMC) benchmark, while future exposures will probe new regions of the parameter space.

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.

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.

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.

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.

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.

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.

Quantum gravity impacts Standard Model fields through effective interactions generated by renormalization. These appear as Standard Model Effective Field Theory (SMEFT) operators, where the effects of new physics are parametrized by the values of the higher-order SMEFT coefficients. As a step towards predicting these SMEFT coefficients from the asymptotically safe Standard Model with quantum gravity, we investigate the flow of a representative set of four-fermion operators under gravitational fluctuations. We strengthen the evidence for the near-perturbative nature of asymptotic safety by finding that these dimension-six-interactions remain irrelevant, resulting in the prediction of specific values for the dimensionless ratios of SMEFT coefficients in the IR. We also find a mechanism that can make a specific subset of these operators relevant. This subset is determined by symmetry considerations. On this basis, we provide a categorization of dimension-six-SMEFT interactions into two categories. Interactions in the first category are predicted to be non-zero, but Planck-scale suppressed. Interactions in the second category are predicted to be zero at small gravitational coupling, but may become free parameters of the theory at larger gravitational couplings.

We evaluate the Baryon Acoustic Oscillation (BAO) signal in the bispectrum as a tool to detect and characterize the relative velocity effect. We extend the existing framework by presenting an updated model for the redshift-space tree-level bispectrum that comprehensively incorporates all relative velocity terms. We introduce a novel, unbiased technique to extract the isotropic BAO dilation parameter ($α_{\rm iso}$) solely from the bispectrum monopole. Validated against N-body simulations, this template-based extraction successfully recovers the acoustic scale and enhances the statistical constraining power by $\sim30\%$, when analyzed in tandem with the pre-reconstruction power spectrum, offering a powerful complement to standard post-reconstruction pipelines. We quantify how individual relative velocity components distort both two- and three-point statistics. We find that these effects induce distinct systematic shifts in the extracted $α_{\rm iso}$ between the two probes. We find systematic discrepancies of up to $2\%$ for $b_{v^2}=\pm0.05$ and up to $20\%$ for $b_{δ^{bc}}\le-2$. This differences demonstrate that a direct comparison of independent power spectrum and bispectrum BAO measurements can break parameter degeneracies and isolate the amplitude of these biases. Finally, we provide a concrete prescription to detect and constrain the three associated relative velocity bias parameters, showing that a joint analysis is highly sensitive to the $b_{v^2}$ and $b_{δ^{bc}}$ amplitudes. This establishes the bispectrum BAO as a robust cosmological probe for current and next-generation galaxy surveys, serving both as a cross-check for standard analyses and a crucial diagnostic tool against systematic biases.