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

Binary black holes (BBHs) have been theorised to form through various different channels, but current LIGO-Virgo-KAGRA (LVK) sky locations make the direct study of formation channels a challenge. Association of strongly lensed BBHs to their hosts has been shown to be possible with current detector networks at upgraded sensitivity and survey telescopes like Euclid, with localisations that can reach sub-galactic scales in the source plane. With Hubble Space Telescope-like high-resolution imaging, we report that strong-lensing localisation of BBHs allows us to confidently associate the BBH with specific structures within their host galaxies, specifically, galactic centres (GCs). We demonstrate that a BBH localised to the GC of the galaxy can reject light-tracing BBH populations ($p < 5\%$) in a single observation, and in most cases demonstrate high confidence rejection ($p < 1\%$), providing a conservative basis to exclude light-tracing spatial distributions associated with field-like or non-central formation channels. Such association therefore directly supports formation channels confined to GCs, such as active galactic nuclei and nuclear stellar clusters.

The accretion-induced collapse (AIC) of carbon-oxygen (CO) white dwarfs (WDs) accreting from helium-star (He-star) donors in ultracompact binaries represents a previously unexplored gravitational-wave (GW) source channel that differs fundamentally from the well-studied oxygen--neon (ONe) WD AIC pathway. We perform Monte Carlo population-synthesis calculations for <mml:math><mml:mrow><mml:mn>5</mml:mn><mml:mo>×</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math> post-AIC CO-WD + He-star binaries spanning three He-donor sub-classes: Class A (low-mass He donors, <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>He</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.35</mml:mn></mml:mrow></mml:math>--<mml:math><mml:mrow><mml:mn>0.70</mml:mn><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math>, AM CVn progenitors), Class B (medium He donors, <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>He</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.70</mml:mn><mml:mo>-</mml:mo><mml:mn>1.20</mml:mn><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math>, ultracompact X-ray binary progenitors) and Class C (massive He donors, <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>He</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1.20</mml:mn></mml:mrow></mml:math>--<mml:math><mml:mrow><mml:mn>2.50</mml:mn><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math>, double-detonation survivors). CO WD masses span <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>WD</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>1.05</mml:mn></mml:mrow></mml:math>--<mml:math><mml:mrow><mml:mn>1.35</mml:mn><mml:mo>,</mml:mo><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math> and initial separations <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>a</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.01</mml:mn><mml:mo>-</mml:mo><mml:mn>0.50</mml:mn><mml:mspace></mml:mspace><mml:mi>R</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math>, placing these systems directly in the LISA millihertz band at birth. Natal kicks are drawn from a Maxwellian with <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>σ</mml:mi></mml:mrow><mml:mrow><mml:mi>kick</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>5</mml:mn><mml:mspace></mml:mspace><mml:mi>kms</mml:mi></mml:mrow></mml:math>, reflecting the smaller mass loss (<mml:math><mml:mrow><mml:mi>∆</mml:mi><mml:mi>M</mml:mi><mml:mo>≲</mml:mo><mml:mn>0.01</mml:mn><mml:mo>,</mml:mo><mml:mi>M</mml:mi><mml:mo>⊙</mml:mo></mml:mrow></mml:math>) of CO WD AIC relative to the ONe WD channel. We find that approximately 97% of binaries remain bound following collapse, with more than half born directly within the 1--20 mHz LISA sensitivity band and most of the remainder evolving through it during their inspiral. This behaviour contrasts sharply with the ONe-WD AIC channel, where only a small fraction of systems become observable by LISA. The three donor subclasses produce distinct chirp-mass distributions that are well separated from both the Galactic DWD foreground and stellar-mass black-hole binaries, providing a unique observational signature. We estimate subclass merger-rate densities of <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>AMC</mml:mi></mml:mrow></mml:msub><mml:mo>≈</mml:mo><mml:mn>0.04</mml:mn></mml:mrow></mml:math>--0.35, <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>U</mml:mi><mml:mi>C</mml:mi><mml:mi>X</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub><mml:mo>≈</mml:mo><mml:mn>0.09</mml:mn></mml:mrow></mml:math>--0.72 and <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>He</mml:mi></mml:mrow></mml:msub><mml:mo>≈</mml:mo><mml:mn>0.18</mml:mn></mml:mrow></mml:math>--<mml:math><mml:mrow><mml:mn>1.40</mml:mn><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>Gpc</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>yr</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math>, corresponding to a total rate of <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mrow><mml:mi>tot</mml:mi></mml:mrow></mml:msub><mml:mo>≈</mml:mo><mml:mn>0.31</mml:mn></mml:mrow></mml:math>--<mml:math><mml:mrow><mml:mn>2.47</mml:mn><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>Gpc</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>yr</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math>. At a Galactic-centre distance of <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mi>L</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>8</mml:mn><mml:mspace></mml:mspace><mml:mi>kpc</mml:mi></mml:mrow></mml:math>, we predict approximately 5--20 individually resolvable LISA sources from this channel. These results identify ultracompact CO-WD + He-star AIC binaries as a promising and previously unrecognized population of millihertz GW sources whose detection would constrain CO WD ignition physics, common-envelope evolution, stable mass transfer, and He-star binary evolution.

Little Red Dots (LRDs) display singular properties unlike any known class of AGN or galaxies, motivating novel mass estimators for their central engines. Inspired by their similarities to stellar phenomena, here we interpret the LRD continuum as being produced by a pseudo-photosphere. We fit tailored stellar atmosphere models to host-subtracted LRD central engines ("black hole stars," BH*s) represented by stacks of $117$ objects. Typical BH* continuum spectra are well fit by models in a narrow range of temperatures ($T_{\rm eff}\approx4200-4800$ K), with bolometric luminosities $\approx10^{43-45}$ erg s$^{-1}$, implying pseudo-photospheric radii $\approx700-2000$ au. Based on these parameters, we explore four different approaches to deriving BH* masses: 1) using the surface gravity from atmosphere models; 2) appealing to the resemblance to super-Eddington phenomena; 3) approximating the escape velocity from the outflowing material; and 4) exploiting the lack of variability to bound the dynamical time. For the typical BH*, all of these methods yield remarkably consistent masses of $\approx10^{4-5}\,M_\odot$, implying a highly super-Eddington luminosity of $L_{\rm{bol}}/L_{\rm{Edd}}\sim5-50$. These mass estimates place BH*s within the scatter of the local scaling relation between black hole mass and host galaxy stellar mass, providing a self-consistent alternative to "overmassive" black holes that lie $2-3$ dex above it. Crucially, our derived masses are consistent with BH*s arising from single supermassive stars (SMSs), whose masses cannot exceed $\approx10^{5-6}\,M_\odot$ due to general relativistic instabilities. Furthermore, for our derived $L_{\rm bol}/L_{\rm Edd}$, the sharp cutoff of the LRD luminosity function matches the maximum theoretical mass of an SMS. With LRDs, we may therefore be directly observing the birth of heavy black hole seeds.

The present thesis is a work in progress on improving the simulation of integrated pulse profiles of pulsars. Pulsars, highly magnetized rotating neutron stars, serve as precise cosmic clocks useful for studying gravity and the interstellar medium. Although periodic pulses from many pulsars are observed and modeled at differ- ent radio frequencies, a robust realistic simulation of their profiles utilizing their parameters remain underdeveloped. The present study focuses on developing a physically consistent model to reproduce observed pulse shapes across multiple frequencies. This study will be useful in coherently understanding the various as- pects of the pulsar emission mechanism. These integrated pulse profiles are the time-averaged properties of pulsar emission and are obtained after averaging indi- vidual pulses over a few thousand rotations. Observed to be highly stable, these are indicative of the global properties of the pulsar magnetosphere. By refining simulation techniques and incorporating the propagation effects of the interstellar medium, this work will prove to be helpful in generating synthetic, yet realistic models. The results are expected to contribute to the interpretation of observa- tional data and to the development of improved timing and emission models for pulsars

Long-term pulse-profile variability has emerged as a potentially important source of uncertainty for high-precision pulsar timing, as changes in the integrated pulse shape can bias pulse times of arrival and affect the sensitivity of pulsar timing arrays to nanohertz gravitational waves. We present a systematic search for long-term pulse-profile variability in the six-year MeerKAT Pulsar Timing Array data, analysing observations of 84 millisecond pulsars using a two-dimensional Gaussian-process framework to identify coherent temporal and phase-dependent profile evolution. We detect pulse-profile variability in 18 pulsars. Using the frequency dependence of the observed profile evolution together with scattering simulations, we classify six pulsars as exhibiting chromatic variability consistent with interstellar scattering, while the remaining twelve display predominantly frequency-independent behaviour suggestive of intrinsic changes in the pulsar emission process. Correcting for the S/N-dependent detection bias, we infer that intrinsic profile variability is substantially more common than implied by the detected sample, with at least 40 per cent of millisecond pulsars exhibiting intrinsic profile variability at the 95 per cent credible level. Motivated by the association between pulse-profile variability and spin-down-rate ($\dotν$) switching in canonical pulsars, we search for corresponding $\dotν$ variations but find no statistically significant detections. Simulations indicate that the expected variations remain below the sensitivity of current datasets, although unresolved changes may contribute to weak achromatic timing noise. These results suggest that long-term pulse-profile variability may be common among millisecond pulsars and should be explicitly modelled in future high-precision pulsar-timing experiments.

Bow-shock pulsar wind nebulae (PWNe) are synchrotron sources formed when the outflow of supersonic pulsars is confined by the surrounding interstellar medium. These sources are dominated by freshly injected particles, thus providing a unique laboratory for studying particle acceleration in relativistic outflows. The Mouse is a prototypical bow-shock PWN and is bright in the radio and X-ray bands, enabling detailed multi-wavelength spectral modeling. Using data from 17 telescopes, the Mouse is detected from 118 MHz to 353 GHz, covering over 3 decades in frequency, and in X-rays from 0.5 to 30 keV. A clear spectral break is identified at 3.7 GHz. The radio spectrum exhibits a rising trend at low frequency, peaks at the break, and then declines. We find a change in spectral index $Δα=0.65\pm0.05$ across the break. We attribute this to an intrinsic cutoff or steepening in the injected particle distribution. This could indicate the characteristic energy of leptons leaving the pulsar magnetosphere, or energy dissipation after the particles cross the termination shock.

We present the results of the first systematic search for refractive scattering substructure in active galactic nuclei (AGN) using ground-based very long baseline interferometry (VLBI) data. We selected 1074 AGNs showing a steep frequency dependence of their angular sizes and compared the observed long-baseline visibility amplitudes with the predictions of the refractive substructure model for scatter-broadened images. Nine sources showed consistency with the expected level of refractive noise, including the quasar TXS 2005+403 (J2007+4029), which confirms previously reported results for this source. Remarkably, the majority of these AGNs are located in the Cygnus region -- one of the key ``hotspots'' of strong Galactic scattering. For five selected AGNs, we identified at least two epochs with sufficiently high data quality showing long-baseline detections, suggesting the presence of scattering substructure. We showed that refractive scattering substructure in AGN is a rare effect, requiring high sensitivity, sufficient source brightness, and lines of sight with particularly strong turbulence in the interstellar medium. Nevertheless, AGNs can serve as powerful probes of the turbulent interstellar medium. We discuss strategies for improving the sensitivity of VLBI observations and identify the most promising sources for reliable detection and investigation of refractive substructure in future observations.

Nearly three decades after the first SOHO/EIT observations, reported EUV-wave speeds still range from a few tens to more than 1000 km/s. This Perspective asks what those numbers actually trace. A measured image speed depends on passband and line-of-sight weighting, cadence, difference imaging, the selected crest or leading edge, propagation sector, fitted interval, and projection geometry. Same-event comparisons quantify this observation-operator effect. In events common to the Nitta and Muhr analyses, early or fastest-sector measurements were typically higher; five of six directly shared events differed by about 30%. For 19 May 2007, the 171 Angstrom peak speed was 475 +/- 47 km/s, compared with 238 +/- 20 km/s in 304 Angstrom, while the 171 Angstrom cadence was four times faster than the 195 Angstrom cadence. A sampling audit of the 1 April 2017 SWAP event shows that a repeated speed near 834 km/s is almost exactly one 91.6-Mm radial ring per 110-s image interval. In a reconstructed 3 April sector, AIA 171 Angstrom gives about 405 km/s, compared with the published SWAP mean of 484 km/s; an independently selected AIA 193 Angstrom ridge gives about 250 km/s under strongly shortened exposures. Such differences are not automatic evidence for different MHD modes. A same-event energy pilot separates compact-source heating, dimming/ejecta, and one weak-compression front segment. The current data do not justify a universal power law, but define a falsifiable question: is the fraction of released energy carried by the front scale invariant from compact quiet-Sun events to global waves and shocks? Establishing this dependence would show which eruption scales can contribute materially to coronal heating.

High-resolution UV spectroscopy serves a diversity of science cases, from small bodies to planets, stars, and galaxies, but is currently limited to the Hubble Space Telescope and bright targets. Major advances require increasing sensitivity by at least one order of magnitude. Here we present the UV science cases for PEGASUS (Planets, Earths, Galaxies, And Stars UV Spectrograph), a UV Échelle high-resolution spectrograph concept, with $R = λ/δλ\sim 100\,000$ (full range 10 000-140 000) and covering 90--400 nm, with a foreseen extension to at least 800 nm. PEGASUS is ideally suited for the Habitable Worlds Observatory (HWO), enabling transformative science across the UV/optical wavelength ranges. PEGASUS will be unique in high sensitivity (effective area) and high spectral resolution -- an uncharted territory -- as well as robustness, thanks to the simplicity of its design. Its UV science cases include: I) Formation and evolution of planets and their habitability: properties of exoplanets and atmospheres, protoplanetary disks, Solar System bodies; II) Stellar lives and deaths at their extremes: the first stars and the origin of the elements, compact and massive stars, Supernovae; III) Gas and metals in the baryon cycle of galaxies: the interstellar, circumgalactic, and intergalactic medium and their roles in galaxy growth. These are essential for the Astro Decadal 2020 Survey, Voyage 2050, and HWO. While this paper focuses on high-impact science enabled by UV high-resolution spectroscopy, PEGASUS will extend into the optical regime and lower spectral resolution, making it a multi-purpose, widely used, workhorse spectrograph for HWO.

Constraint-satisfying initial data for binary black holes require the repeated solution of a nonlinear elliptic boundary-value problem throughout the physical parameter space. We present a para\-metric, differentiable initial-data model that combines sparse-grid spectral interpolation with Newton refinement over a family of binary configurations. Interpolation is used only to provide warm starts for the underlying elliptic solver, rather than to replace it. A small number of Newton iterations refines each interpolated solution to a discrete, row-equilibrated constraint residual $\|R\|_\infty\le10^{-10}$, independent of the interpolation error, so that every returned initial-data set satisfies the discretized constraint to that tolerance. Our JAX implementation provides analytic derivatives with respect to the physical parameters, enabling gradient-based parameter targeting, optimization, and sensitivity analysis. We demonstrate the method for quasi-circular binary-black-hole initial data using conformally flat Bowen--York punctures by constructing four-dimensional aligned-spin and eight-dimensional general-spin models. The resulting initial data are validated against TwoPunctures, and the method is applied to certified targeting of the physical Arnowitt--Deser--Misner mass and angular momentum and to effective-potential eccentricity reduction using a dedicated separation--tangential-momentum model.

Globular clusters are efficient factories of recycled pulsars, but searches toward high-dispersion-measure (DM) clusters can be strongly limited near 1 GHz by dispersive smearing and interstellar scattering. We present the first MeerKAT S-band (nu ~ 2.4 GHz) pulsar survey of 14 globular clusters. High time- and frequency-resolution observations were searched using segmented acceleration and jerk techniques, followed by candidate folding and targeted folding with available timing ephemerides. We re-detected 39 known pulsars and discovered four new millisecond pulsars in Glimpse-C01: J1848-0129C, D, E, and F. Multi-epoch follow-up enabled preliminary Keplerian orbital fits for J1848-0129C and J1848-0129D. J1848-0129C is an eclipsing MSP in a ~5 d orbit, placing it among long-period eclipsing systems known as huntsman binaries, while J1848-0129D is in a ~3.4 d nearly circular orbit with a massive white-dwarf companion of about 1 solar mass. Two Glimpse-C01 pulsars show large DM offsets from the cluster average. Comparison with other Galactic globular clusters indicates that intracluster DM spreads tend to increase with foreground DM, implying that narrow DM search windows may be sub-optimal for high-DM clusters. From detections and non-detections, we infer a practical single-epoch MeerKAT S-band tied-array detectability scale of about 10-20 microJy. These results demonstrate the value of high-frequency searches for pulsars in strongly dispersed and scattered cluster environments.

We present the second paper in our program characterizing the impact of massive neutrinos on the multiscale cosmic web using global topology and persistent homology. Building on the methodology established in Paper I, based on discrete Morse theory, we analyze a subset of the QUIJOTE simulations to compute persistent diagrams, Betti curves, and additional topological statistics for both dark matter and halo density fields, across redshifts z = 0, 1, 2. A central result of our study is the first clear demonstration that apex points in persistent diagrams are especially sensitive to neutrino mass, with enhanced sensitivity for specific pairs of saddle points at high redshift. In addition, Betti curves from dark matter density fields broaden and flatten with increasing neutrino masses, exhibiting two characteristic density thresholds where Betti numbers remain invariant. These mass-dependent signatures are detectable at the few-percent level, even for M<SUB>ν</SUB> ~ 0.1 eV, providing a robust, physically grounded probe of massive neutrinos in the cosmic web. While traditional two-point statistics encode only pairwise correlations and cannot fully break parameter degeneracies, persistent homology captures higher-order, multiscale information that can lift these degeneracies. Moreover, its high sensitivity to the sum of neutrino masses makes it a promising complement to conventional analyses. Our results thus establish a solid foundation for forward-modeling or emulator-based approaches using persistent homology and environment-based statistics to constrain neutrino mass---potentially enabling direct detection---and additional cosmological parameters, with immediate relevance for ongoing and upcoming galaxy surveys, including DESI, Euclid, and Rubin-LSST.

We present ABCMB, a differentiable Einstein-Boltzmann solver for the cosmic microwave background (CMB). ABCMB is an analysis-ready code package capturing important effects to linear order in ΛCDM cosmology. It computes the CMB power spectrum and includes effects like lensing, polarization, massive neutrinos, and a state-of-the-art treatment of Big Bang Nucleosynthesis and recombination. ABCMB has sub-percent-level agreement with CLASS and can be run on a GPU with competitive, and sometimes even faster, run times, owing to ABCMB's favorable run time scaling with maximum multipole moment. It is refactored compared to previous codes and takes advantage of object-oriented programming to improve extensibility, meaning new physics can be added to it without the need for modifying source files. ABCMB provides accurate and stable gradients to the user, making Fisher analyses straightforward, and enabling the use of efficient gradient-based sampling methods.

Neutrinos and dark energy (DE) have entered a new era of investigation, as the latest DESI baryon acoustic oscillation measurements tighten the constraints on the neutrino mass and suggest that DE may be dynamical rather than a cosmological constant. In this work, we obtain a high-confidence measurement of the neutrino mass within a dynamical DE framework. A joint analysis of DESI DR2, cosmic microwave background, DESY5 supernova, and DESY1 weak lensing data yields a total neutrino mass of <inline-formula><mml:math><mml:mo>∑</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>ν</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mn>0.098</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.037</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>0.016</mml:mn></mml:mrow></mml:msubsup><mml:mspace></mml:mspace><mml:mtext>eV</mml:mtext></mml:math></inline-formula>, indicating a measurement for a non-zero, positive neutrino mass at the 2.7σ level within the w<SUB>0</SUB>w<SUB>a</SUB>CDM framework. This high-confidence measurement is driven mainly by these factors: (i) the DESI's preference for a dynamical DE with its equation of state evolving from w &lt; -1 at early times to w &gt; -1 at late times, thus leading to a larger neutrino mass; (ii) treating N<SUB>eff</SUB> as a free parameter together with the inclusion of weak lensing data, which likewise allows for an increased neutrino mass. In the future, even higher-confidence measurements of neutrino mass are expected with stronger preferences for dynamical DE in light of more complete DESI data releases.

We propose a single-field mechanism that dynamically prepares the initial conditions required for exponential plateau inflation. The inflaton is decomposed into a coarse-grained mode and a high-occupancy, approximately Gaussian band of relativistic nonzero modes. Gaussian--Hartree averaging of the potential produces an exact variance-dependent minimum, which places the coarse-grained field high on the plateau when the excitation variance is large. Under a translation-invariant microcanonical measure, the available phase-space volume is exponentially concentrated near this minimum, while an order-one homogeneous kinetic-energy fraction is exponentially atypical. During expansion, the preparation modes redshift as radiation, $σ^2\propto a^{-2}$ and $ρ_χ\propto a^{-4}$. Although the induced minimum rapidly moves toward the vacuum, the force acting on the coarse-grained mode is bounded by an exponentially small factor, yielding analytic upper limits on its velocity, displacement, and kinetic fraction. We derive a sufficient condition for the subsequent number of slow-roll e-folds and exhibit a benchmark with more than 60 e-folds. The mechanism also motivates correlated finite-onset signatures: large-scale power suppression with decaying oscillations and a variance-convolution bispectrum.

The McVittie metric describes a self-gravitating compact object embedded in an expanding universe. It inherits two singularities: a curvature singularity at $r \to 0$ and a cosmological singularity at $a(t) \to 0$. Motivated by the construction of black-bounce geometries, we regularize both of these singularities and propose a regularized McVittie metric. The central singularity is regularized by replacing $r \to \sqrt{r^{2}+b^{2}}$, leading to a geometry that interpolates between a cosmological black hole, a black bounce, and a traversable wormhole. Similarly, the cosmological singularity is regularized by introducing a non-vanishing scale factor $a \rightarrow \sqrt{a^{2} + a_b^{2}}$. The resulting spacetime is supported by an effective imperfect fluid with finite anisotropic stress. We investigate the formation of the apparent horizon, analyze the circular geodesic structure, and argue that the geometry can also be interpreted as a conformally evolving Morris-Thorne wormhole embedded in a regular cosmological background.

Phenomenologically similar sign-switching dark-energy histories need not have comparable scalar-field realizations. We reconstruct minimally coupled scalars with fixed kinetic sign for three prescribed histories: the error-function model (ECDM), the compact smooth-step model (SSCDM), and the ladder-like model (L$Λ$CDM). Continuous negative-to-positive density crossings select the phantom branch. ECDM yields a smooth on-shell potential at every finite redshift; its equation-of-state pole at the density zero is only a ratio singularity. Exact SSCDM has a regular trajectory but a $C^1$, non-$C^2$ endpoint potential, $V-V_e\propto|φ-φ_e|^{4/3}$. Its non-Lipschitz force permits delayed departures from frozen plateaus, so the reconstructed potential and plateau data do not uniquely generate the prescribed history. The exact Ladder requires distributional kinetic stress and has no ordinary classical realization in the adopted one-field action. In conditional synthetic comparisons, the sigmoid--Gaussian family ranks highest for ECDM and the generalized axion-like family for SSCDM, although the fitted axion exponents $n<1/2$ imply a divergent endpoint force. Representative regular forward solutions exhibit negative-to-positive scalar-density crossings, with the potential zero preceding the density zero. In a closure test, unretuned evolution of the top-ranked ECDM template tracks the target approximately, $\max|Δ\widetildeΩ_φ|\simeq0.16$. Closure fails structurally for compact SSCDM: no potential with a locally Lipschitz force can reproduce its exact finite-duration frozen plateaus from exactly frozen initial data. Field-map existence, invertibility, and endpoint regularity must be assessed before interpreting a phenomenological history as scalar dynamics. The phantom action is used only as a homogeneous effective proxy.

We present a spectral test of the quasi-star hypothesis for Little Red Dots (LRDs) whereby a black hole grows inside a stellar-like envelope. We use \texttt{Prospector} to fit host galaxies and \texttt{TLUSTY} photospheres to four LRDs that show strong molecular or atomic absorption. We approximate each object's electron-scattered Eddington luminosity ratio, $φ\equiv κ_{\rm es}σT_{\rm eff}^{4}/(gc)$, and we use \texttt{MESA-QUEST} to simulate their envelopes. All four are super-Eddington at $φ= 4.0$--$299$, although the largest value sits at the edge of our atmosphere grid and beyond our simulations. We derive envelope masses between $700$--$27{,}000\,M_\odot$, where quasi-star theory requires the black hole to be less than a third of that. GN-28074 falls five decades below its published virial mass estimate, alleviating the overmassive black hole problem. Their black holes double in mass every $\sim0.05$--$3.5$~Myr and can produce intermediate-mass black holes in $\lesssim30$~Myr. The two fitted components of the water-absorbing object WIDE-EGS-2974 share a single $φ$, consistent with our hypothesis that the absorption originates in an extended atmosphere. We then extend our measurements to 82 archival LRDs, finding super-Eddington photospheres throughout the population that launch winds at $1$--$3\times$ their own escape speeds, which implies that LRDs evolve from being massive, continuum-driven sources to more eruptive winds as they mature and shed their outer envelopes. We thus constrain the Eddington ratios and masses of LRDs, and, for the first time, show self-consistently that quasi-stars may be the central engines powering LRDs.

The relative abundances of elements in galaxies serve as fossil records of the physical conditions and processes by which they were forged. While the Big Bang produced only the lightest elements, subsequent stellar nucleosynthesis imprinted characteristic abundance patterns onto the surrounding gas, set initially by the temperatures reached inside stars and subsequently shaped by how the processed material was mixed and released. Globular clusters - dense, ancient groups of stars - provide a striking unique example. Some contain stars depleted in magnesium and enriched in aluminum, showing that they formed from gas exposed to exceptionally hot hydrogen burning. The stars responsible remain unknown. Little Red Dots may provide this missing engine. These compact, luminous objects formed at cosmic epochs similar to those associated with globular-cluster formation and are enshrouded by dense gas whose chemical composition can be measured with the James Webb Space Telescope. Here, using deep spectroscopy from the SPURS program, we show that this abundance pattern characterizes the LRD central engine: magnesium-depleted and aluminum-enhanced gas with a metallicity only 1% that of the Sun. This pattern is not produced by ordinary massive stars at these redshifts and cannot be mimicked by ionization, gas geometry or dust. Instead, it is reproduced by hot hydrogen burning in fully convective supermassive stars, with the measured abundances implying masses of at least 10,000 solar masses - approximately 100 times larger than any star observed in the present-day Universe. Little Red Dots may therefore reveal supermassive stars during their brief lives or in the immediate aftermath of their direct collapse, simultaneously identifying the long-sought source of globular cluster abundance anomalies and a formation pathway for massive black hole seeds.

These notes trace large-scale structure from primordial curvature perturbations generated by inflationary quantum fluctuations to galaxy power-spectrum multipoles. Three core lectures develop the linear matter power spectrum, spherical and anisotropic collapse, galaxy bias, redshift-space distortions, the Kaiser model, and multipole estimators with Gaussian covariance. The extension develops nonlinear bias and the one-loop effective field theory model used in full-shape analyses. Derivations are explicit; appendices collect longer calculations and solutions. The core lectures assume undergraduate-level cosmology; the extension assumes familiarity with perturbation theory.

We present an analytical, statistical solution to the binary-binary (2+2) outcome of the chaotic non-hierarchical four-body problem. The solution is based on the density-of-states formulation pioneered by J. J. Monaghan. The method skips the computationally expensive integration of the equations of motion, and instead samples the outcome from the chaotic phase-space, subject to conservation of energy, momentum, and angular momentum. From the joint distribution, we extract marginal distributions of several key parameters using Monte-Carlo integration, and numerically verify them by comparing to an identical ensemble of scattering experiments produced by the FEWBODY code. From the comparison, we identify a regime not represented by the density-of-states formulation: the hard binary regime, where one binary is much harder then the four-body energy scale, and the system acts as an effective three-body system. We hypothesize that at this regime the system's probability distribution spreads over an effective reduced three-body chaotic phase space. The process of transfer from four-body to three-body phase-space is still not understood, and represents the next natural extension of density-of-states methods.

We study the gravitational collapse of a shell, modeled as a coherent state of quantum matter. Without adding any physics beyond semiclassical gravity, we show that the universal gravitational coupling between matter and geometry induces an interaction between the shell and the negative-frequency modes present in the $|in \rangle$ vacuum state. We find that this interaction occurs unavoidably immediately after a trapping horizon forms, making it a universal feature of semiclassical gravitational collapse. The resulting transfer of energy from the collapsing matter to outgoing Hawking radiation acts as a \textit{gravitational quantum friction}, opposing the formation of trapped regions while simultaneously imprinting the information carried by the collapsing matter onto the emitted quanta. We discuss the implications of this mechanism for the black hole information-loss problem and speculate about the possible endpoints of collapse, potentially preventing horizon and singularity formation.

Matter fields localized around black holes can produce small secondary bumps in the effective potentials governing black-hole perturbations, potentially leading to quasinormal-mode spectral instabilities and late-time echoes. Considering physically motivated models of anisotropic fluids surrounding a spherically symmetric black hole, we show that the appearance of such a secondary bump is accompanied by either singular matter profiles or divergent characteristic speeds in the fluid. These features indicate that the corresponding matter configurations are unphysical. We further demonstrate that a secondary bump in the effective potential does not necessarily imply the existence of additional light rings in the spacetime.

Triangular loop signals appear frequently in cosmological collider models but their analytical forms remain unknown due to multiple technical challenges. Here, we determine for the first time the leading analytical cosmological collider signals generated by a massive scalar triangle loop in inflationary bispectrum and trispectrum. Combining a directional cutting rule with partial Mellin-Barnes representations, we obtain both the local and nonlocal trispectrum signals in a soft limit and derive the leading signal in the squeezed bispectrum. Remarkably, all these signals are reproduced by the corresponding bubble diagram with a single effective pinched coupling. Signals beyond leading order and cases with arbitrary internal masses, spins, and interactions can be generated within the same framework. Our approach also provides a systematic route to extracting analytical loop signals from general cosmological correlators.

Chemical tracers provide some of the strongest observational signatures of ongoing planet formation and localized dynamical perturbations in protoplanetary disks. In particular, sulfur-bearing molecules are predicted to be enhanced in regions of shock heating, ice sublimation, and gravitational instability. Here, we present high-angular-resolution ($\approx$0.$^{\prime\prime}$2) Atacama Large Millimeter/submillimeter Array observations of $^{12}$CO J=3-2, $^{13}$CO J=3-2, CS J=7-6, and SO J$_{\rm N}$=8$_8$-7$_7$ toward the large, edge-on Gomez's Hamburger ('GoHam'; IRAS 18059-3211) disk. We detect a narrow, one-sided arc of SO emission that peaks near a previously-identified gas over-density, suggesting localized heating around an early-stage giant protoplanet or disk fragment. The edge-on geometry of GoHam enables us to place this chemical signature in the broader context of the disk gas and dust structure. To do so, we map the vertical distribution of molecular gas relative to millimeter- and (sub)-micron-sized dust, identify a pronounced north-south continuum asymmetry, and detect non-Keplerian $^{12}$CO and $^{13}$CO emission indicative of a disk wind. We also derive a dynamical stellar mass of 2.2 $\pm$ 0.5 M$_{\odot}$ and a revised dust-extinction-map-based distance of 139 $\pm$ 24 pc, which places GoHam in the outskirts of the Scorpius-Centaurus association. Together, these observations reveal a highly dynamic disk in which localized sulfur chemistry may trace one of the earliest observable stages of wide-separation giant planet formation.

In this work, we investigate the physical properties and cosmic evolution of galaxies and AGN in a sample of 280 protoclusters identified in the TNG300 simulation, selected to end up in $z = 0$ clusters with $M_{200, \rm c} > 10^{14} \, {\rm M_{\odot}}$. Our goal is to provide the first statistical view of AGN activity enhancement in a uniformly defined sample of overdense environments, and to identify the physical mechanisms driving it. We identify protoclusters as the progenitors of present-day galaxy clusters through merger-tree reconstruction and compare their galaxy and AGN populations with a control sample of field galaxies across the redshift range $0 \leq z \leq 6$. We investigate galaxy and SMBH demographics, AGN fractions, bolometric luminosity functions, and the SMBH accretion rate density, consistently applying homogeneous selection criteria in all environments. We find that TNG300 protoclusters host systematically more massive galaxies and SMBHs than the field by up to $1 \, {\rm dex}$ at all redshifts, with signatures of accelerated galaxy evolution already visible at $z \sim 3-4$. The AGN fraction increases with stellar mass in both environments and, at fixed host-galaxy stellar mass, is broadly consistent between protoclusters and field galaxies, indicating no strong environmental triggering of SMBH accretion. However, when analysed as a function of redshift, protoclusters exhibit a significant enhancement of AGN activity (by a factor $> 2$), particularly at high luminosities and early cosmic times. We show that this enhancement primarily arises from differences in the stellar-mass distributions of galaxies in overdense regions, where massive systems assemble earlier than in the field. Consistently, protoclusters dominate the bright end of the AGN bolometric luminosity function and contribute up to $\sim 50\%$ of the total SMBH accretion rate density at $z \sim 6$.

N-body simulations are an essential tool for modelling the evolution of gravitating systems. Differentiable N-body simulations allow solving complicated reconstruction problems by connecting observations of evolved systems to simple informative priors on their initial conditions. While so far most reconstruction efforts focus on mildly non-linear large scale dynamics, the internal dynamics of galaxies and clusters should exhibit rich information about their formation. The exploration of reconstruction possibilities in this highly non-linear regime naturally requires highly performant differentiable N-body simulations with small scale forces. Here, we present JZ-FMM, a GPU-native differentiable implementation of the fast multipole method for N-body simulations. We show that the GPU oriented design of tree structure, dual tree traversal and multipole translation operators improves performance over established CPU-native and hybrid N-body codes by notably more than an order of magnitude. Further, we show that gradients can be evaluated efficiently and accurately so that JZ-FMM can be used as an essential building block in future field-level reconstruction efforts. We apply the new code to the reconstruction problem of a tidally stripped satellite to show that deeply non-linear problems can indeed be solved efficiently with simulation gradients, although appropriate care must be taken to navigate the complicated optimization landscape.

We present the identification of the co-spatial and co-temporal small scale jets -- sometimes also named as jetlets -- at the base of the coronal hole plume associated with chromospheric and transition region (TR) flows. We identified 19 jets with coordinated observations recorded on March 19, 2016 using the Solar Dynamics Observatory (SDO) and the Interface Region Imaging Spectrograph (IRIS). Our analysis of \ion{Si}{4} 1393.7~Å line suggests blueshifts mostly in the range of 10 to 33~km~s$^{-1}$, but few less than -10~km~s$^{-1}$, while analysis of co-spatial \ion{Mg}{2}~h~\&~k lines show both blueshifts and redshifts in the range of {--}9 to 17~km~s$^{-1}$. We observed highly asymmetric, enhanced intensity in \ion{Si}{4} line during these events. The upflows observed in transition region through \ion{Si}{4} is strongly correlated with chromospheric downflows observed in \ion{Mg}{2}. We interpret these results as a signature of interchange reconnections creating bi-directional flows -- giving rise to upflows (downflows) in transition region (chromosphere).

We investigate the high-energy collision of two spinless bodies within the Self-Force Effective Field Theory, taking one of the two bodies to be massless and highly energetic, such that it generates a gravitational shockwave (GSW) background. The GSW is perturbed by a second body, either massless or massive, satisfying respectively $E_2/E_1\ll1$ or $m_2/E_1\ll1$. Although the Self-Force expansion for shockwaves is frame-dependent, we choose a frame in which this EFT manifests its full power, obtaining results valid to all orders in Newton's constant that make the underlying resummation patterns manifest. Specifically, we push the expansion to 1SF order and conjecture the all-order vanishing of the so-called recoil operators in a given gauge. We obtain the exact graviton propagator on the shockwave background, matching previous results, and use it to compute the on-shell waveform at 1SF order in momentum space, specialising to the massless-massless collision. The waveform is computed to all orders in the Post-Minkowskian expansion, and its resummed form proved to vanish at the leading order in the super-Planckian regime. We use this result to prove the UV finiteness of the angular spectrum of radiated energy in the region collinear to the shockwave with energy $E_1$.

We present spatially resolved stellar kinematics of the deflector galaxies in two time-delay lens systems, B1608+656 and SDSS J1206+4332, measured from JWST NIRSpec integral field spectroscopy. B1608+656 is a quadruply imaged quasar lensed by a pair of tidally interacting galaxies, the primary deflector (G1) and a close companion (G2), for which we quantified G2's tidal influence on G1 and define a criterion identifying the region where this perturbation is negligible and standard Jeans-equation modeling remains valid. SDSS J1206+4332 is a doubly imaged quasar whose extended host galaxy is itself quadruply imaged by a foreground massive galaxy (G0), accompanied by a fainter perturber (G1). For both systems, we report the first spectroscopic redshift and velocity dispersion measurements of the perturber galaxies. For G2 in B1608+656, we measure σ_v = 188 +/- 7 (stat.) +/- 2 (sys.) km/s at z = 0.6307, identical to the redshift of the main deflector G1. For G1 in SDSS J1206+4332, we find σ_v = 100 +/- 21 (stat.) +/- 2 (sys.) km/s at z = 0.7403, slightly lower than the primary deflector G0. The resolved velocity maps of both primary deflectors show clear rotation signatures and we visually classify them as fast or regular rotators. The data products will be made publicly available and will be combined with time delays, lens models, and line-of-sight convergence to measure cosmological parameters in the upcoming TDCOSMO 2026 milestone publication.

We study linear gravitational perturbations and the corresponding quasinormal-mode spectra of electrically and magnetically charged black holes in the presence of higher-order corrections to the Einstein-Hilbert action. In particular, we consider a four-derivative interaction that couples the Riemann tensor to the electromagnetic field strength in the unique (Horndeski) combination that preserves second-order equations of motion. We provide a comprehensive analysis of the stability conditions. We numerically confirm the absence of exponentially growing modes throughout the physically admissible region of parameter space where the theory remains free of pathologies. We then investigate in detail the properties of the quasinormal-mode spectra and identify several distinctive features induced by the higher-order interaction.

White dwarfs (WDs) provide an unparalleled opportunity for detecting and characterizing exoplanets because of the large planet-to-star size ratio. However, only one transiting WD exoplanet is confirmed to date, and occurrence rate constraints remain loose because the rapid transits and weak phase curve modulations of WD exoplanets require high-precision, high-cadence photometry. The Galactic Bulge Time Domain Survey (GBTDS) of the Nancy Grace Roman Space Telescope will provide exactly this by observing thousands of WDs at 12.1-minute cadence. The expected WD exoplanet yield of this survey is not yet quantified. Here, we predict the GBTDS WD exoplanet yield over its first three high-cadence seasons, through both transiting and non-transiting (phase curve) channels. We run Monte Carlo detection probability simulations over a grid of planet radius, orbital period, and host properties for a Besançon synthetic population of WDs. Jupiter-sized planets dominate the transit yield, with a period-averaged effective survey size of $\overline{N}_\mathrm{eff}\approx8$, requiring an underlying occurrence rate of $\overlineη_1\approx12\%$ for one expected detection, followed by $\overline{N}_\mathrm{eff}\approx2.6$ ($\overlineη_1\approx38\%$) for Neptune-sized and $\overline{N}_\mathrm{eff}\approx1.8$ ($\overlineη_1\approx55\%$) for sub-Neptune-sized planets. Source blending is the major observational challenge in the crowded GBTDS fields, removing $\approx95\%$ of the raw predicted yield. Phase curve provides marginal detections, adding at most $\approx4\%$ to the transit yield. A null detection by GBTDS would constrain $η_\mathrm{WD}$ to $< 28\%$ for Jupiter-sized planets at $95\%$ confidence. GBTDS will deliver the first WD exoplanet occurrence-rate constraints beyond the solar neighborhood and identify prime targets for future characterization.

Inspired by a recent ground-breaking space-based measurement of Newton's gravitational constant, we outline a proposal for a general layout of an optimized space-based approach that combines the advantages of a largely disturbance-free environment with lock-in techniques, null coordinates and the method of covariances. The method of covariances, and cross-covariances, is used for an optimized statistical analysis of data originating from two distinct experimental channels. One of two distinct measurement channels relies on an approximately constant acceleration in the gravitational field of static sources masses, while the second channel involves a lock-in measurement based on test and source masses with oscillating coordinates. Conceivable mission parameters are analyzed in terms of the optimum scale for the source masses and launch capabilities of currently available space launch systems. A candidate location for a largely disturbance-free position in geostationary orbit is identified.

We investigate axial gravitational perturbations and shadow formation in black-bounce geometries supported by anisotropic fluids within general relativity. Two distinct classes of solutions are analyzed: a symmetric Simpson--Visser-like black bounce and an asymmetric deformation of the Reissner--Nordström metric, which can feature either a bounded interior region or an unbounded wormhole-like extension. Our findings indicate that configurations with horizons, whether symmetric or asymmetric, exhibit single-barrier effective potentials and quasinormal frequencies nearly indistinguishable from those of standard black holes. In contrast, horizonless symmetric configurations may develop multiple potential barriers, giving rise to gravitational-wave echoes whose amplitude and separation depend sensitively on the bounce parameter $a$ and the density parameter $ρ_0$. Horizonless asymmetric solutions, on the other hand, exhibit regularized effective potentials that do not produce echoes and remain phenomenologically close to the Reissner--Nordström case. From the optical perspective, both classes of solutions with horizons exhibit a single relevant photon sphere and shadow profiles nearly indistinguishable from those of their black-hole counterparts, while horizonless symmetric configurations can sustain multiple photon rings. Horizonless asymmetric configurations may instead support either a single photon ring or none. Finally, we recover the standard eikonal correspondence between quasinormal modes and the instability of circular null geodesics for configurations with horizons. For horizonless geometries, however, the emergence of multiple potential barriers prevents a straightforward extension of this relation, although the strong similarity between the wave and optical sectors suggests that a meaningful phenomenological correspondence may still persist.

We perform field-level likelihood-free inference of the matter density parameter $Ω_m$ from simulated galaxy catalogs using machine learning models with differing inductive biases. Using hydrodynamic simulations from CAMELS, we examine how observable choice and architecture govern cosmological information extraction. We consider galaxy positions and line-of-sight peculiar velocities, separately and jointly, and compare permutation-invariant Deep Sets, implemented with either multilayer perceptrons (MLPs) or Kolmogorov-Arnold Networks (KANs), to graph neural networks (GNNs), which explicitly encode spatial relations. We test in-distribution and out-of-distribution (OOD) performance across simulations with different subgrid galaxy-formation prescriptions. Deep Sets infer $Ω_m$ from velocities alone with mean relative errors of approximately $18\%$ in-distribution and $\sim25\%$ OOD, with KANs and MLPs achieving comparable performance. In contrast, the same set-based approach does not yield useful $σ_8$ predictions in either in-distribution or cross-suite tests. Adding positions does not improve Deep Sets, while GNNs infer $Ω_m$ with mean relative errors of about $10\%$ in-distribution and $10$--$17\%$ OOD. These results indicate that peculiar velocities provide the dominant source of $Ω_m$ information for set-based models in this setting, while spatial information is most effectively used by architectures that explicitly encode galaxy-galaxy relations. Because the velocity inputs are exact simulated peculiar velocities, applications to survey data will require validation under realistic velocity-measurement noise, selection effects, and survey geometry.

We show that a ringdown waveform can admit a constant-coefficient quasinormal-mode (QNM) representation while its individual QNM-pole contributions remain driven. To establish this result, we construct a first-principles rational approximation to the Kerr Green's function from its QNM poles, physical residues, and horizon-frequency zeros. For particles plunging from the innermost stable circular orbit into Kerr black holes with spins 0.5-0.9, each QNM-pole contribution initially follows the source's instantaneous complex frequency and decouples only after the source decays faster than the corresponding mode. At late times, the source approaches a sum of damped oscillations at the third and higher horizon frequencies. The corresponding zeros in the Green's function cancel these source-frequency components in the coherent waveform, while the continuing drive generates oscillations at QNM frequencies. Thus, an apparently free QNM superposition can emerge before its individual contributions dynamically decouple from the source.

Planetary nebulae (PNe) that are physically associated with star clusters are rare but crucial for stellar evolution studies because they allow their initial and final stellar properties to be determined independently, contributing to constraints on the initial-to-final-mass relation (IFMR), which traces stellar mass loss. Although the identification of such pairs is essential, most reported cases have been identified incidentally rather than through systematic searches. Here, we attempt to uncover new such associations in our Galaxy employing literature PN central stars (CSs) and star cluster catalogues. Exploring relevant angular positions on the sky and kinematic Gaia DR3 data, we identify 38 candidate CS-cluster pairs. The components of eight of these pairs present consistent interstellar extinctions, and their CS properties are compatible with stellar evolution models. Two of these candidate pairs are already reported in the literature as true associations. We evaluate separately the remaining six pairs using literature PN observational data, assessing the likelihood of each pair being a true association. The initial and final masses of our best candidates are mostly consistent with the latest IFMR predictions and confirmed cluster-PNe points. An additional 46 CS-cluster pairs have components with compatible angular positions but lack kinematic CS data and are thus presented separately. More data are needed for a precise estimation of the physical parameters of most of our candidates and to confirm these associations, which will be addressed in future studies.

Context. The China Space-station Survey Telescope (CSST) will provide deep, wide-field epoch astrometry during its 10-year mission. Astrometric binary orbits constrain the masses of stellar and compact-object components. Orbital recovery depends on astrometric precision and temporal coverage. Combining CSST and Gaia data extends the baseline and improves binary detection. Aims. We evaluate CSST, Gaia, and joint astrometry for binary-candidate selection and 12-parameter (12p) orbit fitting at faint magnitudes ($g>17.8$). We also test how regular CSST cadences affect the yield of 12p fits satisfying our criteria. Methods. We constructed a mock catalog, simulated CSST and Gaia epoch astrometry, and fitted five-parameter (5p) single-star models to derive astrometric diagnostics, proper-motion anomaly features, and observational-sampling features. A four-stage histogram-based gradient-boosting classifier used these features to select candidates for 12p orbit fitting and assessment. Results. On the independent test set, the classifier reaches a precision of 0.802 and a recall of 0.181 among eligible true binaries. In the scenario-specific fitted samples, joint astrometry raises the fiducial fraction from 6.76% for Gaia alone to 10.46%; for fitted binaries with $P_{\rm true}>15{\rm yr}$, it rises from 2.37% to 6.78%. The current CSST schedule yields few fiducial fits, while idealized regular cadences increase the yield mainly at $g\lesssim21$. Conclusions. In the simulation, joint CSST and Gaia epoch astrometry yields higher fractions of fitted unresolved binaries satisfying the stated criteria than Gaia-only solution. A practical strategy is to select candidates from 5p diagnostics and astrometric anomalies, obtain more regular CSST follow-up observations, and then fit 12p orbital models and apply the selection criteria.

Extremely metal-poor (EMP; [Fe/H] < -3.0) stars preserve key information about the earliest stages of Galactic chemical evolution. Based on a training set of spectroscopic metallicities, we constructed catalogs of more than 120,000 EMP candidates from large-scale narrow-band photometric surveys. To validate this approach, we obtained high-resolution CFHT/ESPaDOnS spectra for 13 candidates selected from the SkyMapper-based catalog. The follow-up observation confirms that the narrow-band photometric selection is effective in identifying very metal-poor stars and retaining a substantial EMP fraction: four targets are confirmed as EMP stars, and all targets remain very metal-poor with [Fe/H] < -2.4. The photometric metallicities are systematically lower than the spectroscopic values by 0.44 dex, with a relatively small scatter of 0.19 dex, indicating a systematic offset in the photometric metallicity scale that could be reduced through improved calibration. Abundances for more than 20 chemical species are derived, leading to the identification of a new potassium-enhanced r-II star, J2001-1215. By combining chemical abundances with orbital properties, five stars are found to be dynamically consistent with known Galactic substructures, including Gaia--Sausage--Enceladus, Thamnos, and Sequoia. These stars provide high-resolution chemical measurements for the extremely metal-poor regime of these Galactic substructures. These results underscore the role of narrow-band photometry in efficiently selecting EMP candidates for targeted high-resolution spectroscopic follow-up, thereby enabling detailed chemical and dynamical studies of the early Milky Way.

We present an explicit three-parameter family of toric Einstein-Maxwell gravitational instantons. These are complete Riemannian manifolds with vanishing scalar curvature which satisfy the Riemannian Einstein-Maxwell equations. We argue that the metric can be extended to a space with two asymptotic ends diffeomorphic to $\mathbb{H}^2 \times S^2$ with a non-collapsing two-cycle (bolt) in the bulk. The global metric is produced by suitably extending an analytic continuation of a local family of Lorentzian black hole geometries constructed by Ovcharenko-Podolsky.

We investigate the low-frequency behavior of the amplification factor in gravitational lensing and explore how it encodes information about the density profile of the lensing object. We derive the low-frequency expansion of the amplification factor under the Born approximation for a broad class of projected density profiles. For spherically symmetric profiles that decay faster than any power law at large distances, we derive a systematic expansion of the amplification factor in powers of frequency, with the logarithmic dependence appearing only in the leading term, and show that each expansion coefficient is determined by a finite set of moments of the density profile. We then extend the analysis to profiles with power-law tails and demonstrate that such profiles induce additional non-analytic frequency dependences, including fractional powers and logarithmic terms, which directly reflect the asymptotic behavior of the density distribution. Furthermore, we investigate the effects of non-sphericity and show that contributions from the quadrupole moment appear only as higher-order corrections relative to the spherically symmetric component in the low-frequency regime. Finally, we investigate the validity of the Born approximation in the low-frequency expansion. We derive a criterion for the maximum order of the low-frequency expansion up to which the Born approximation remains dominant over the post-Born corrections.

We report the discovery of a nearby FU Ori-type outburst (FUor), PR Ori B, in the L1641 cluster of the Orion star-forming region. The high-amplitude variability was first identified in the NEOWISE (3-5 $μ$m) photometry of the unresolved PR Ori binary system. Long-term, resolved optical photometric monitoring demonstrates that PR Ori B is the driver of a $ΔG=5~$mag outburst, while PR Ori A has remained constant over the last 20 years. The near-IR spectrum of PR Ori B changes from a late K-type spectral type during quiescence to a viscously heated disk during outburst, including deep absorption in $^{12}$CO and H$_2$O bands. The optical spectrum also exhibits features that are commonly associated with FUors, including P Cygni profiles in Na I D lines and absorption in the Ca II infrared triplet. The luminosity of the outburst (L$_{acc}\sim$30--40 L$_\odot$) is similar to that commonly observed in FUors. The comparison of Spitzer/IRS and VLT/VISIR spectroscopy shows some evidence of silicate crystallisation during the outburst. PR Ori B is one of the closest and brightest FUors discovered over the last few years, only one magnitude fainter than the archetype of the class FU Ori. The proximity and brightness will allow for future high angular resolution observations to probe the physics of the inner disk and to evaluate changes in the disk due to the increased luminosity.

We introduce CROCODILE-SIDM, a framework for treating self-interacting dark matter (SIDM) with the $N$-body part of GADGET4-Osaka code, as part of CROCODILE simulation family. As a test case, we investigate the impact of SIDM on the tidal formation of dark matter-deficient galaxies (DMDGs) with velocity-dependent cross-section models. We demonstrate that our implementation reproduces the analytic scattering rate in isolated halos. Including dynamical friction self-consistently, we evolve a dwarf satellite with $M_*=2\times10^8\,\mathrm{M}_\odot$ in a $\sim10^{11}\,\mathrm{M}_\odot$ halo on a decaying orbit around a massive host, comparing CDM with four SIDM cross sections for two initial satellite density profiles: a cuspy Navarro--Frenk--White (NFW) profile and a cored Burkert profile. We find that self-interactions primarily regulate the amount of DM retained between pericentric passages, and that the sign of this effect depends on the initial profile: a larger cross section retains more DM for the Burkert initial condition but less DM for the NFW initial condition. We show that this opposing behavior reflects the direction of SIDM heat conduction, which is set by the halo's evolutionary state at infall. Core formation in the cuspy profile assists DM stripping, whereas tidally accelerated gravothermal contraction in the cored profile suppresses tidal mass loss. Consequently, SIDM can either assist or hinder DMDG formation, depending on the satellite's inner structure, making DMDGs a potential probe of SIDM cross section.

The characterization of Sun-like stars, including determining whether those stars are single or have stellar companions, has become increasingly important for missions seeking planets beyond the Solar System. With upcoming missions like the Habitable Worlds Observatory (HWO) on the horizon, vetting stars for binarity and confirming potential targets as viable is essential precursor science. In this paper, we present an automated algorithm that searches for double-lined spectroscopic binaries (SB2s) in high-resolution spectra. Using MOLUSC, we simulate realistic stellar populations to quantify the capability of algorithm for detecting stellar companions around nearby stars, with potential applications to HWO target stars and Gaia DR4 planet-host candidates. In addition, we explore how the detection rate is influenced by complementary observations like high-resolution adaptive optics imaging and time-series radial velocities. We find that our SB2 detection algorithm alone can rule out 46.1% of stellar-mass companions. High-resolution adaptive optics alone rules out 79.6% of stellar companions. Radial velocity surveys alone ruled out 98.4% and 98.7% for 5-year and 10-year surveys, respectively. When all three methods are combined, the probability of a typical nearby star having a stellar companion is <1%. The combination of the methods allows us to better understand and resolve target lists for future surveys.

The motion of a spinning test particle in a static, spherically symmetric spacetime with a cosmological constant is studied within the Mathisson-Papapetrou-Dixon formalism, working to linear order in the particle's spin. By exploiting the symmetries of the background, the radial and latitudinal dynamics are reduced to first-order equations governed by a spin-deformed quintic polynomial. A complete classification of the real root structure of this quintic is carried out, from which the allowed orbital types for spinning particles in the Schwarzschild-anti-de Sitter and Schwarzschild-de Sitter geometries follow directly. The parameter dependence of the root structure and of the stable circular orbit interval is analyzed. In particular, for values of the cosmological constant beyond the Stuchlík limit, where spinless particles admit no bound orbits, spinning particles can still remain bound. The radial motion is expressed piecewise in terms of Lauricella hypergeometric functions, with the simple positive real roots of the quintic acting as branch cuts. The latitudinal motion consists of small oscillations about an equatorial plane, and its phase, governed by the same quintic, admits an analogous analytic representation.

Globular clusters (GCs) are well known to host stellar populations characterized by light-element variations. A subset of Galactic GCs, beyond such 'canonical' populations, contains an additional group enriched in iron, s-process elements, and total C+N+O abundance (i.e., the anomalous stars). We combine multi-facility photometry with APOGEE spectroscopy to investigate the spatial distribution, chemical properties, and formation history of the stellar populations in M22, with particular focus on its anomalous component. We trace the canonical and anomalous populations, together with their subpopulations, from the cluster center out to eight half-mass radii. The anomalous population becomes increasingly dominant in the outermost regions, whereas no significant radial gradients are detected among its subpopulations. Our chemical analysis reveals light-element anticorrelations within both the canonical and anomalous components, although the latter are shifted toward higher C, N, and Al abundances. For the first time, we show that the Fe, s-process, and C+N+O enhancements among anomalous stars are not uniform but correlate with their light-element composition: the most chemically extreme anomalous stars are also the most Fe-, Ce-rich, and C+N+O-poor. We identify a distinct red overdensity on the horizontal branch, likely populated by the most He-poor stars, and tentatively associate the extreme horizontal branch with the most chemically enriched anomalous population. These observations are difficult to reconcile with M22 being a merger between two GCs. Instead, they qualitatively favor a self-enrichment scenario regulated by dilution, similar to that recently proposed for $ω$Centauri, with their close chemical correspondence suggesting that they experienced analogous formation histories.

We investigate the influence of current helicity on the turbulent magnetic diffusivity $β$ using three complementary derivations of the $α$ and $β$ coefficients, based on the large-scale magnetic field $\overline{\mathbf{B}}$, the turbulent velocity $\mathbf{u}$, and the turbulent magnetic field $\mathbf{b}$. Applying these coefficients to raw DNS data, we reconstruct $\overline{\mathbf{B}}$ and compare the results with the original simulations. In the kinematic regime all models agree well with the DNS data. In the nonlinear regime, however, $β_{\mathrm{vv-vw}}$ alone produces unbounded growth of $\overline{\mathbf{B}}$. Including the contribution from turbulent magnetic fields ($β_{\mathrm{bb+jb}}$) suppresses this unphysical growth and restores agreement with the DNS results. We find that kinetic helicity drives $β$ more negative, while current helicity shifts it back toward zero. Weighted combinations of the coefficients further show that the $β$ effect dominates the evolution of $\overline{\mathbf{B}}$ throughout, whereas the $α$ effect becomes important mainly for sustaining the field in the nonlinear regime. The corresponding IDL analysis scripts are provided to facilitate practical implementation of the theoretical models.

$γ$ Cas stars are Oe/Be stars that exhibit bright and hard X-ray emission. HD45314 belonged to this category, but lost its $γ$ Cas characteristics as its circumstellar disc started to fade away. The star was monitored in optical spectroscopy, optical photometry, and X-ray spectroscopy to follow its variability on timescales of hours to years. Time series of magnitudes, colours, line equivalent widths and violet over red peak intensity ratios were analysed with Fourier methods to uncover possible periodicities. For five years, HD45314 exhibited oscillations of its magnitude and line strengths on a timescale of 230 d. It then declined towards a nearly disc-free stage and is now slowly rebuilding its disc. During these phases, the violet over red peak intensity ratio exhibited a modulation on a timescale of about 1190 d. Over the whole campaign, prominent short-term photometric variations, notably at a frequency of 3.369 d$^{-1}$, were probably due to $β$ Cep-type pulsations. The amplitude of these pulsations was stronger during phases of overall brightness changes. In parallel, the X-ray spectrum switched to a low state at the onset of the oscillation phase and has remained in this low state since then. Some residual hard emission is nonetheless observed. HD45314's long-term optical variations bear resemblance with the predictions of smooth particle hydrodynamics calculations, although a close look does reveal some inconsistencies. Unlike some other $γ$ Cas stars, where the X-ray emission remained nearly constant despite major changes of the circumstellar disc, HD45314's X-ray properties changed dramatically as the disc begun to dissipate. These differences might arise from a wider orbital separation between the Oe star and its putative white dwarf companion or from the ablation of the disc by the radiation field of the Oe star.

In this work, we extend the thermodynamic analyses of the neutral and specially charged Kerr-Bertotti-Robinson black holes to the general Kerr-Newman-Bertotti-Robinson family, in which the electric and external-field parameters are independent. Using covariant surface charges and canonical integrability methods, we determine the total angular momentum and the canonical mass. The angular momentum follows analytically from the combined gravitational and electromagnetic surface charges. However, the infinitesimal charge associated with coordinate time translations is not integrable in solution space, so the mass must be associated with a more general symmetry generator. Imposing the canonical integrability conditions on this generator, together with the Kerr-Newman mass as the zero-field boundary condition, selects the Christodoulou-Ruffini mass and determines the associated thermodynamic potentials. The first law and Smarr formula take the same form as the ones in usual Kerr-Newman case, and the two previously studied Kerr-Bertotti-Robinson cases are recovered as special limits.

In space-based gravitational-wave detectors, laser phase locking transfers the phase of weak received light to a local laser and maintains the interspacecraft heterodyne frequencies within the phasemeter bandwidth. Fluctuations of the onboard ultrastable oscillators introduce clock noise into these heterodyne phase measurements, requiring additional calibration in time-delay interferometry (TDI). Existing formulations based on six independently formed one-way measurements do not explicitly describe the common-reference clock structure of a master--slave phase-locking configuration. In this work, we formulate clock-noise propagation and calibration directly for this configuration. Starting from carrier and sideband readouts, we derive the clock-noise couplings and construct an ordered sideband-calibration rule. The rule retains noncommuting delays and applies to two-branch TDI combinations satisfying propagation and clock-coefficient closure. For mutually independent onboard clocks with identical fractional-frequency noise spectra, we compare matched realizations with and without phase locking. In the adopted frozen equal-arm model, 45 second-generation combinations with up to 16 links considered here have smaller clock-noise residuals in the phase-locking realization throughout 0.1--10~mHz, before final clock calibration. Numerical simulations of two representative 16-link combinations reproduce the analytical spectra and verify the predicted reduction and calibration. These results show that the phase-locking realizations reduce the clock-noise contribution before final calibration, while the cancellation of master-laser phase noise still relies on TDI.

Context. Scattering polarization in the Sr I 4607 Å line is a key diagnostic for small-scale, unresolved magnetic fields in the quiet solar photosphere, inaccessible to Zeeman-based techniques. The complete frequency redistribution (CRD) limit is commonly used to model this line, but its validity and impact on magnetic sensitivity have not been systematically investigated. This is timely given new facilities observing this line with unprecedented accuracy, and future synoptic programs targeting long-term quiet-Sun magnetism. Aims. We assess the CRD limit against the general partial frequency redistribution (PRD) description of scattering for the Sr I 4607 Å polarization, focusing on its magnetic sensitivity via the Hanle effect. Methods. We solved the non-LTE radiative transfer problem for polarized radiation in a semi-empirical 1D plane-parallel static atmosphere, with CRD and PRD scattering, across magnetic field configurations relevant for Hanle diagnostics. Results. For small-scale unresolved fields, where the field acts only via Hanle depolarization of Q/I, CRD and PRD emergent profiles and Hanle sensitivity are essentially identical. For deterministic fields, CRD remains accurate when polarization signals exceed about 0.3%; for weaker signals, PRD effects can appreciably affect the emergent profiles. Conclusions. CRD is adequate for modeling Sr I 4607 Å scattering polarization in the vast majority of observationally relevant cases; PRD effects matter only for weak polarization signals.

New time series photometry of the short period ($P<0.2$~d) binary stars ZTF~J032906.36+070408.3 (ZTF~0329) and ATO~J090.7238+00.5996 (ATL~0602) is presented. The primary (hotter) stars in both systems are white dwarfs, while the cooler stars are M dwarfs. Archival photometry and new and archival spectra are used to derive the temperatures of the components in the two binaries. Both are single-lined binaries. Velocity curves are fitted to radial velocity measurements of the red dwarfs and used to place lower limits on the masses of the white dwarfs. New and archival time series photometry show evolution in the shapes of the light curves of both stars. Spectra of ZTF~0329 show an excess of radiation in the $R_C$ band which is ascribed to electron cyclotron radiation, suggesting that the star is a low accretion rate polar. Satisfactory model fits to the multi-filter light curves of ATL~0602 require the presence of cool spots on the red component of the binary. The latter star exhibits a variety of velocities associated with Balmer and Ca~II emission lines within individual spectra.

We present a detailed stellar evolution study of BD$-18^{\circ}$5550, an extremely metal-poor red giant star in the Galactic halo. Such stars are among the oldest known objects and preserve chemical signatures of the earliest stages of Galactic history. BD$-18^{\circ}$5550 exhibits $α$-enhancement, $\rm[α/Fe]=+0.4$, typical of Population II halo stars, with enhanced abundances of oxygen, magnesium, and other $α$-elements relative to iron. Using the \texttt{MESA} stellar evolution code, we perform the first direct determination of its mass, age, radius, and luminosity by simultaneously fitting eight observed stellar properties and multiple chemical abundance ratios. Our analysis yields a current mass of $\rm M=0.76,M_{\odot}$, radius of $R=38.4,R_{\odot}$, luminosity of $\log L=2.8,L_{\odot}$, and an age of $t=11.87$ Gyr, with a 7-fold improvement in age precision. The model reproduces all observed constraints, supporting our treatment of stellar physics, nuclear reactions, and internal mixing. We further investigate the predicted surface chemical composition up to iron and compare it with observations. The surface abundances of elements up to iron remain essentially unchanged throughout the evolution, while nitrogen shows gradual enrichment, indicating the onset of mixing processes at intermediate RGB stages. The preservation of the other surface abundances provides a valuable window into the chemical conditions of the early Galactic environment. Our results establish BD$-18^{\circ}$5550 as a benchmark star for constraining stellar evolution and early Galactic chemical evolution. Its derived age places it among the oldest stars in the Milky Way, providing important constraints on the epoch of early star formation and the initial assembly and evolution of our Galaxy.

Can propagating gravitational waves serve as operational probes of the quantum nature of gravity? We address this question by developing a unified theoretical framework that combines spacetime geometry, quantum information, and gravitational-wave physics. Starting from the geodesic deviation equation in linearized General Relativity, we derive the effective detector Hamiltonian directly from spacetime geometry and construct the complete quantum dynamics for detector subsystems interacting with both classical and quantized propagating gravitational-wave fields. This unified formulation enables a direct comparison between classical and quantum descriptions of gravitational radiation within the same physical framework. We demonstrate that classical gravitational-wave backgrounds can induce mixedness in the detector state but cannot generate genuine quantum correlations between the detector subsystems. In contrast, quantized gravitational waves coherently mediate gravity-induced entanglement, quantum coherence, quantum memory, and nonclassical correlations, providing clear operational signatures of the quantum nature of propagating gravitational radiation. We further discuss how mesoscopic quantum mechanical oscillators offer a promising route towards experimentally probing these effects. Our results establish a geometric and quantum-information-based framework for exploring quantum gravity through propagating gravitational waves.

The cosmic web, consisting of an intricate network of voids, walls, filaments, and nodes, encodes key information about structure formation and the cosmological parameters that govern it. In the era of high-precision cosmology, large ensembles of numerical simulations are required to analyse next-generation galaxy surveys, motivating the use of generative models to circumvent their high computational cost. While they have shown promise in emulating high-fidelity cosmic web simulations, their ability to capture distinct cosmic web environments remains largely unexplored. For this study, we trained a diffusion model on the Quijote N-body simulation suite to investigate the semantic information learnt by its self-attention maps. Using statistical estimators such as the Dice coefficient and cross-power spectra, we quantified the correspondence between attention maps and cosmic web environments defined by the T-Web classifier. We find that attention maps of varying spatial resolutions across different layers capture overdense and underdense structures in distinct ways, exhibiting strong positive correlations and anti-correlations with both the overall matter distribution and individual cosmic web environments. Moreover, the diffusion model predominantly encodes cosmological information at intermediate-to-large spatial scales, indicating that attention maps primarily capture globally coherent structures. Our results show that, beyond accurately reproducing two-point statistics, diffusion models learn a multi-scale representation of the cosmic web through self-attention, including non-Gaussian information.

The millihertz gravitational-wave sky will be dominated by $\sim10^{4}$--$10^{7}$ resolved double white dwarf binaries in the Milky Way, whose three-dimensional spatial distribution traces the Galaxy's structural parameters. We present \textsc{Galena}, a hierarchical Bayesian pipeline that recovers these parameters from a double white dwarf catalogue through an inhomogeneous Poisson-process likelihood. Each source's $3\times3$ Galactic position covariance is computed from the waveform Fisher information matrix and propagated analytically through a $3\times5$ Jacobian; the resulting measurement-error convolution is factored into a pre-computed sparse weight matrix, rendering the nested-sampling inference tractable. Applied to GBSIEVER-reported LISA Data Challenge sources ($N=2151$ after quality cuts), \textsc{Galena} recovers the disk scale length $R_d=2175^{+51}_{-52}$~pc and scale height $z_d=282^{+7}_{-7}$~pc, consistent with a canonical thin disk, together with the bulge fraction $A=0.187^{+0.011}_{-0.012}$ and bulge scale radius $R_b=773^{+26}_{-25}$~pc; the bulge fraction, now recovered much closer to the literature value of $0.25$ than in our earlier exact-position fit, is measured at $\sim6\%$ statistical precision. An independent particle-swarm optimisation of the same likelihood reproduces these values to within a fraction of a percent, supporting the attribution of the improvement to the Fisher-matrix error propagation rather than to the sampler. The search-stage astrophysical prior improves per-source distance estimates but leaves the hierarchical inference essentially unchanged.

In this work, we model the long-term impact of recurrent stellar flares on the atmospheres of metal-rich gaseous exoplanets. Using synthetic flare spectra from a fiducial flare model integrated with a photochemical kinetics code, we track the changes in atmospheric composition with time. We further analyze the spectral variability by feeding these abundance profiles into a radiative transfer code at various time steps. Our simulations showed variability and persistent changes in key atmospheric species, such as CH4, CO2, and SO2, when compared to their quiescent state. Extreme flare events cause rapid depletion of molecules in the upper atmosphere and a temporary disappearance of spectral features, especially the SO2 feature at 7-8 microns, which shifted by about 75 ppm. Many species did not fully return to their quiescent state after flares, resulting in lasting changes in abundance, especially for SO2 and CO2, key species when inferring the atmospheric metallicty. We also explored the cumulative effects of recurrent flares, showing that species like H2O and CH4 followed a decreasing abundance trend, with half-lives of around 28 to 31 years. These results indicate that flare activity plays a significant role in shaping both the short- and long-term atmospheric composition and spectral features of giant gaseous exoplanets orbiting M-dwarf stars, underscoring the need to account for stellar activity when characterising such atmospheres. Our findings also highlight that the atmosphere of the modelled planet is not static, suggesting that a probabilistic approach to atmospheric abundances may be more appropriate than static retrievals, particularly for gaseous planets orbiting active stars.

We establish a minimal nonrelativistic realization of scalarization in a self-gravitating bosonic condensate coupled to a scalar response field. A local effective-mass shift and nonlinear saturation generate a multibranch structure with two distinct transition routes. In the linearly stable regime, a finite perturbation drives a first-order transition with type-I logarithmic scaling near threshold. Beyond the linear onset, small perturbations grow tachyonically, and the response time controls both the growth and the subsequent breathing dynamics. These results identify a common mechanism for scalarization across relativistic compact objects and nonrelativistic condensates, and point toward laboratory analogues in coherent media.

We investigate the impact of self-interactions on the Rényi and entanglement entropies of a scalar field on $(2+1)$-dimensional spacetimes, whose spatial sections are modeled by fuzzy spaces, specifically the fuzzy sphere and the fuzzy disc. We compute the first-order perturbative correction induced by a $λφ^4$ interaction using the Green's function approach. In contrast to the free theory, where the entanglement entropy is dominated by degrees of freedom near the entangling boundary and obeys an area law, we find that the interaction correction has an extensive bulk contribution, receiving significant contributions from degrees of freedom throughout the fuzzy space. For the fuzzy sphere, the correction exhibits strong infrared sensitivity associated with the zero mode. We isolate and resolve this zero-mode IR divergence by projecting out the zero mode, thereby obtaining a physically meaningful quantity. In the commutative continuum limit, the interaction correction has the same degree of UV divergence as the free entropy but does not obey a pure area law. Furthermore, we analyze the Moyal plane limit, where the interaction correction exhibits a distinct IR divergence. We discuss the physical origin of these extensive bulk features and examine their possible connection to the celebrated UV/IR mixing phenomenon in noncommutative quantum field theories.

Motion of plasma in pulsar magnetospheres essentially follows bent magnetic field lines which may result in radiation beams sweeping past the observer's sightline. We show that under certain conditions such a passing beamlet can produce two orthogonal polarization modes (OPMs) solely from averaging of radiative contributions, without involving birefringence. This requires that radiation from different stages of the passage is contributed incoherently, hence that the timescale for the coherence is much shorter than the typical sampling time. The model produces two OPMs of comparable amount, with a wide range of (single-pulse) polarized fractions, which is consistent with general observed pulsar properties. Under the specific conditions required (essentially incoherent summation of signal at appropriate time scale) the averaging effects alone can thus lead to the appearance of two orthogonal polarization tracks with OPM jumps.

Minimal Geometric Deformation (MGD) applied to a static Schwarzschild black hole seed generates a single decoupler function $h(r)$, obtained by solving the $θ$-sector field equations together with an equation of state. Once $h(r)$ is fixed, the resulting one-parameter family is controlled by the coupling strength $k$ through $F(r;k)=1+k\,h(r)$. We show that, whenever the deformation develops a simple outermost root that crosses the seed horizon, the same fixed decoupler leads to two mutually exclusive branches associated with different global completions: on one side of the critical coupling the deformed metric preserves the seed horizon as a black hole, whereas on the other side the root $r_*>2M$ lies in the exterior and cannot be interpreted as an interior modification of the black hole geometry. We prove that this root forces a loss of Lorentzian signature on the interval $(2M,r_*)$, so that no smooth extension of the exterior metric through the seed horizon $r=2M$ exists once $r_*$ lies outside it. Within the static, spherically symmetric class considered here, the corresponding smooth Lorentzian completion is a two-ended wormhole obtained by excising $(2M,r_*)$ and doubling the region $r\geq r_*$ across the minimal sphere $\mathcal T=\{r=r_*\}$. No topology change of any single spacetime is claimed or required: $k>k_c$ and $k<k_c$ simply correspond to two different, non-diffeomorphic manifolds, and Lemma~1 below shows that the metric itself dictates which of the two is the admissible completion for a given $k$. We compute the second homology group of both completions explicitly, $H_2(Σ_{\rm BH},\mathcal H)=0$ for the black hole exterior relative to its horizon and $H_2(Σ_{\rm WH})\cong\mathbb Z$ for the completed wormhole manifold, giving a discrete invariant that distinguishes the two branches.

The magnetic fields that make up the internal structure of coronal mass ejections (CMEs) are thought to be organised in a flux rope configuration, consisting of twisted magnetic fields that wind about a central axis. Remote-sensing observations of CMEs show a wide range of morphologies, dynamics, and evolution, including rotation, non-uniform expansion, deflection, and interaction with the ambient solar wind. In-situ measurements, however, typically consist of a single 1D spacecraft trajectory through a large 3D structure (or few at best), restricting our ability to determine how CME magnetic structure varies in space and time. In this work, we analyse the magnetic configuration of two idealised CMEs during their early evolution in the range 1-30 $R_{\odot}$ using the Magnetohydrodynamic Algorithm Outside a Sphere code. The initial flux ropes differ only in their orientation with respect to the global dipole and erupt from a bipolar active region below the streamer belt in a simplified coronal environment. We deploy a fleet of synthetic spacecraft throughout the CMEs' paths at different combinations of latitudes, longitudes, and heliocentric distances. We identify and examine flux rope signatures in the synthetic in-situ profiles to characterise radial variations as well as latitudinal/longitudinal ones. We find that, even in a simplified CME erupting under solar minimum-like conditions, the sampling location significantly affects the global structure that would be deduced from common flux rope reconstruction and analysis techniques used for in-situ measurements. Additionally, we consider possible in-situ signatures that may serve as indicators of the source region's topology.

We constrain a damped harmonic oscillator (DHO) dark-energy equation of state using the full cosmic microwave background (CMB) likelihoods in combination with DESI BAO and three distinct Type Ia supernova compilations: Pantheon+, DES-Dovekie, and Union3. The equation of state obeys a second-order damped oscillator equation in number of $e$-folds, so that its frequency $f$, damping rate $b$, and equilibrium value $w_{\rm m}$ fully specify the late-time dynamics. The model exhibits oscillatory behavior only at low redshifts, around the equilibrium value $w=-1$, with distinct characteristics for the different supernova compilations: an underdamped solution for DES-Dovekie and Union3, and an overdamped solution for Pantheon+. At higher redshifts, the model closely mimics $Λ$CDM and deviates significantly only at $z<0.6$, with the magnitude of the deviation depending on the supernova compilation. We further identify a region of the $(f,b)$ parameter space, corresponding to rapid variation of the equation of state at low redshift, in which the perturbation equations become numerically stiff and cannot be integrated with a canonical dark-energy sound speed i.e., $c_s^2 =1$. We show that reducing the rest-frame sound speed removes this obstruction while leaving the observables unchanged at the $10^{-3}$ level, and therefore treat it as a numerical prescription rather than a physical modification of the model. The model yields $H_0 = 67.53^{+1.22}_{-1.18}$ km/s/Mpc for Pantheon+, $H_0 = 69.08^{+1.23}_{-1.16}$ km/s/Mpc for DES-Dovekie, and $H_0 = 70.67^{+1.89}_{-1.87}$ km/s/Mpc for Union3. The present-day equation-of-state parameter is constrained to $w_0 = -0.521^{+0.891}_{-0.414}$, $-3.01^{+1.14}_{-1.21}$, and $-3.16^{+1.11}_{-1.16}$ for Pantheon+, DES-Dovekie, and Union3, respectively.

There is growing evidence for cosmic-ray acceleration in massive stellar clusters. Furthermore, extended gamma-ray emission suggests that particle transport in the vicinity of their sources is influenced by physical processes markedly different from large-scale diffusion in the Milky Way. We characterize extended gamma-ray emission in the direction of the Westerlund 1 stellar cluster and Kes 41 supernova remnant using > 16 years of data from the Fermi Large Area Telescope (LAT) at energies > 0.8 GeV. We test whether clusters of gamma-ray sources not associated to multiwavelength counterparts are better described by extended emission components. We report the detection of three new extended emission components with soft spectra towards regions of high gas column density in the Galactic plane. One extended component associated with the natal cloud of Kes 41 is statistically preferred over the point source previously reported towards the supernova remnant shell. The other two extended components overlap with neutral gas within ~100 pc from the edge of the Westerlund 1 superbubble. The extended emission may be explained either by mismodeled gas in the interstellar background model or by the local injection of particles. Under the latter hypothesis, explaining the component associated with Kes 41 requires converting $\lesssim 5\%$ of the supernova remnant energy into accelerated particles, while accounting for the two components near Westerlund 1 requires converting $10^{-4}$ of the cluster wind mechanical power into gamma rays. Nevertheless, our results strengthen the evidence for gamma-ray emission structures arising at intermediate spatial scales between isolated objects and the large-scale diffuse emission from the interstellar medium. This could explain a part of the soft unassociated Galactic sources detected by the Fermi LAT. (Abridged)

In the nanohertz band, spatial correlations between pulsar timing residuals provide a key observable for characterizing stochastic gravitational-wave backgrounds and probing their polarization content. In torsionless Palatini spacetime, nonmetricity can generate two additional shear modes, referred to as the shear-$x$ and shear-$y$ modes. In this work, we investigate the pulsar timing response produced by the shear-induced motions of the emitting pulsar and the receiving Earth. Assuming that the electromagnetic field is minimally coupled to the physical metric and that the Earth and pulsar possess non-negligible effective hypermomentum responses, we derive the single-pulsar redshift, the frequency-domain two-point correlation function, and the corresponding spatial correlation. For a stationary and isotropic stochastic shear background, the normalized overlap reduction function (ORF) for distinct pulsars reduces in the short-wavelength limit to the pure dipolar form $Γ_{ab}^{\mathrm{sh}}(ζ)=\cosζ$, where $ζ$ is the angular separation between the two pulsars. The same dipolar correlation can also be produced by isotropic Solar System ephemeris errors, leading to a spatial degeneracy between the two signals. This degeneracy highlights the importance of information beyond the angular correlation for identifying Palatini shear signatures in PTA data.

We present the connected even-parity galaxy 4-Point Correlation Function (4PCF) of DESI Year 1 (Y1) Luminous Red Galaxies (LRGs). This is one of the first measurements of the connected 4PCF on data, and shows a clear detection at $\sim$12 to 17$σ$ in the full-sky analysis. We test the robustness of the signal by varying three factors: hemisphere (north vs. south), redshift range (either the full range $0.4 < z < 1.1$ or the higher-number-density interval $0.4 < z < 0.8$), and sample completeness (the full sample vs. only the most complete regions within the number density-motivated cut). Finally, we cross-correlate different patches that are spatially well-separated, and thus should have largely independent noise. This approach, at leading order, is able to remove mismatch between the covariance matrix of the data and that of the mocks; however, in certain cases at the cost of sensitivity. We find $\sim$15$σ$ evidence for an even-parity 4PCF in this approach. Overall, the studies of the 4PCF presented here probe sensitively any mismatch between mock catalogs and data, and also open up the prospect of fitting a model, constraining cosmological parameters and galaxy biases, and searching for Baryon Acoustic Oscillation (BAO) features.

We investigate whether quantum fluctuations of a circular probe string develop a quantum-chaotic response while traversing a wormhole throat. The classical circular-string embedding is periodic and radially stable, but its two physical transverse polarizations experience time-dependent tidal potentials. Expanding the world-sheet action to quadratic order, we canonically quantize these modes and construct out-of-time-ordered correlator(OTOC) amplitudes from their unequal-time commutators. For the Ellis--Bronnikov wormhole, both polarizations exhibit finite intervals of approximately exponential OTOC growth associated with the first throat passage. The corresponding dimensionless rate measured with respect to physical time is positive over the parameter range studied and generally decreases as the probe energy is increased relative to the throat scale. In the global-monopole extension, increasing the solid-angle deficit narrows the band of locally amplifiable modes and suppresses the extracted rates; a sufficiently strong defect can nearly quench the radial signal, while the angular channel retains a polarization-dependent non-monotonic structure when the energy-to-throat-scale ratio is small. These quantities characterize finite-time dynamical sensitivity in the Gaussian fluctuation sector and should not be identified with asymptotic many-body chaos or a thermodynamic phase transition. Although the numerical analysis uses two representative wormhole geometries, the construction depends only on covariant world-sheet fluctuations and real-time commutators. It therefore provides a transferable, non-holographic framework for applying quantum-chaos diagnostics directly to quantum probes in curved spacetimes.

Manchak and Weatherall (2018) formulate the black hole information paradox as a failure of predictability in black hole evaporation spacetimes, diagnosed by non-global hyperbolicity. I offer a strategy for resolving this paradox. I argue that failures of predictability in black hole evaporation are not well diagnosed by non-global hyperbolicity. I then consider two weakenings of global hyperbolicity: prediction and retrodiction friendliness, the failure of which could ground a new paradox. However, deidealized black hole evaporation models can be prediction and retrodiction friendly. Therefore, the information paradox cannot be based upon failures of global hyperbolicity, nor either retrodiction or prediction unfriendliness.

The results of detailed long-slit spectroscopy of the bipolar planetary nebula \ESO\ are presented. The observations were obtained with the facility spectrograph RSS of the SALT telescope. The total measured extent of the nebula reaches $\approx$55$^{\prime\prime}$, and the kinematics along the slit directly confirms its bipolar structure. The inclination of the nebula to the line of sight is estimated as $i\approx 45^{\circ}$, and the kinematic ages of two morphological structures are measured: the expanding ring and the polar outflow in the direction perpendicular to the ring. The age of the former is $\approx$1.5$\times$10$^4$~yr and that of the latter is $\approx$(3--4)$\times$10$^4$~yr, i.e. the polar outflow is 2--3 times older than the central ring. The ratio of these ages does not depend on the adopted distance. The weighted mean heliocentric velocity of \ESO\ is V$_{\rm hel}$=74.7$\pm$1.6~\kms. Chemical abundances are obtained for 19 regions along the major axis of the nebula, over an extent of 40$^{\prime\prime}$, i.e. for different morphological parts of \ESO. The abundances were derived with the T$_{\rm e}$--method, where the electron temperature was calculated directly using the measured faint auroral lines [\ion{O}{iii}]~$λ$4363~Å and [\ion{N}{ii}]~$λ$5755~Å. The abundances of O, N, Ne, S, Ar, Cl and He were determined. Despite a significant stratification, the measured abundances along the slit can be considered nearly constant.

We construct a holographic dark-energy model using Kaniadakis entropy with the particle horizon as the infrared cutoff, consistently modifying both the holographic density and the Friedmann background. In standard Einstein gravity, noninteracting particle-horizon holographic dark energy (HDE) does not produce the sufficiently negative pressure required for late-time accelerated expansion. We show that the Kaniadakis deformation changes this behavior while the particle horizon is retained. For the representative parameter choice $K=1.90\times10^{-36}$, with $Ω_{\rm DE0}=0.7$ and $c^2=0.64$, the deceleration parameter changes sign at $z\simeq0.585$ in the noninteracting case. Including the interaction strengths $b^2=0.03$ and $0.06$ changes the transition only slightly, giving $z\simeq0.586$ and $0.589$, respectively. Thus, for the parameter set considered here, the interaction does not generate the acceleration; rather, the transition is already present in the noninteracting Kaniadakis model and the interaction produces only a small shift in its timing. We derive the autonomous evolution equations, the effective dark-energy equation of state, and the deceleration parameter. The adiabatic squared sound speed is also examined as a diagnostic of the effective-fluid stability, while statefinder variables are used to characterize deviations from $Λ$CDM. Finally, we verify that the $K\to0$ limit continuously recovers standard Einstein-gravity particle-horizon HDE, for which the noninteracting model remains decelerating.

We develop a Plebański--Demiański-adapted parametrization of the general Ovcharenko--Podolský solution and construct new families of generalized black holes in four-dimensional Einstein--Maxwell theory. These Petrov type~D spacetimes possess non-null, fully non-aligned electromagnetic fields. A restriction imposed in the previous parametrization is not required by the field equations in shifted coordinates. Removing it retains an independent charge parameter $q$, which may be real or purely imaginary. For $q^2\geq0$, both the metric and electromagnetic field admit a smooth aligned limit to the charged Plebański--Demiański solution with $Λ=0$ and $e^2+g^2=q^2$. We impose angular conditions in the Griffiths--Podolský form allowing Killing horizon cross sections with spherical topology. The resulting eight-parameter class extends the seven-parameter class of Ovcharenko and Podolský [arXiv:2508.04850]. The remaining normalization condition is generically quartic, defining four algebraic branches. We give explicit metrics and electromagnetic fields for the non-twisting and $l=0$ twisting families and analyze their black hole and acceleration horizons, extremality, conicity, and further limits. We determine the effective-NUT-free branches of the latter family and their degenerate cases. Using the genuine Griffiths--Podolský form, we establish an explicit correspondence between the non-accelerating effective-NUT-free subclass and the Kerr--Newman--Bertotti--Robinson family and express its electromagnetic flux charges in our parameters. The $q^2<0$ sector includes the genuine uncharged Kerr--Bertotti--Robinson solution recently identified by Ovcharenko and Podolský [arXiv:2608.21672].

Adopting pseudo-random noise (PRN) codes for inter-satellite ranging, the resulting data-encoded PRN signal induces excessive laser phase measurement noise under the currently adopted integrated scheme in space-based gravitational wave (GW) detection systems. To suppress the noise, a scheme of inter-satellite integrated laser link is presented in this Letter, in which the data-encoded PRN signal for inter-satellite absolute distance measurement and communication is phase-modulated onto the ultra-stable oscillator (USO) clock signal before phase modulation onto the laser carrier. Theoretical analysis and simulation results show that in the proposed scheme the noise is significantly suppressed compared with the currently-adopted one, while the noise stemming from the data-encoded PRN signal on clock noise transfer is limited within the requirement.

This work concerns a procedure for removing gauge dependence from graviton corrections to the effective field equations for a massless, minimally coupled scalar on de Sitter background. The procedure involves combining diagrams in the same way as the scattering amplitude for the $t$-channel exchange of the massless scalar between two massive particles, but {\it without} taking the asymptotic limits which are problematic in cosmology. Implementing this at 1-loop on flat space background requires five classes of diagrams, in addition to the naive exchange, and the combination of those diagrams does eliminate dependence on the graviton gauge fixing functional. When those same diagrams are generalized to de Sitter background, their combination cancels the most important (``nonlocal'') secular gauge dependence, but it leaves (``local'') secular gauge dependence that is associated with external legs. Local secular gauge dependence can be removed by adding two additional classes of diagrams, and making secular corrections to the external wavefunctions. The purpose of this paper is to include these additional corrections on de Sitter background in order to derive fully gauge independent effective field equations at 1-loop order.

Quantification of energy transport and dissipation in weakly collisional heliospheric plasmas that are far from local thermodynamic equilibrium is an outstanding scientific problem. A central challenge is determining how non-Maxwellian velocity-space structure affects damping and emission of coherent ion-scale waves, especially compared to simplified analytical models for background plasma velocity distributions. In this work, we study the damping and emission of parallel-propagating proton cyclotron waves for two models of proton velocity distributions measured by the SPAN-I instrument on board Parker Solar Probe during an extended storm of waves with left-hand polarization in the solar wind at a heliocentric distance of 30.1 solar radii. Using the measured velocity distribution rather than a two-component bi-Maxwellian model predicts instabilities consistent with the observed coherent waves. For intervals in which both models predict net damping, the observed VDF model yields weaker damping in 90\% of cases, with a reduction in the integrated heating rate of 0.44 relative to the bi-Maxwellian model. These results suggest that simplified analytical velocity distribution models may overestimate cyclotron damping and underestimate wave emission in the near-Sun solar wind.

Point particles are generically ill defined in fully nonlinear general relativity, although they are naturally well defined in linear perturbation theory. In a previous paper, two of us showed how they can be rigorously defined from matched asymptotic expansions at second order in perturbation theory on a generic background spacetime. Here we show how the resulting field equations can be recast in a more practical, more fully ``skeletonized'' form. We also show that this form of the field equations enables simple derivations of equations of motion. Finally, we show how calculations in second-order perturbation theory can bypass the use of singular ``punctures'' and make rigorous use of off-the-shelf regularization methods such as Hadamard regularization, opening new avenues to second-order self-force calculations in black hole binaries.

Recent hints of dynamical dark energy, including a possible crossing of the phantom divide, have motivated the search for viable theoretical explanations. Scalar fields non-minimally coupled to gravity with a coupling of the form $F(\varphi)=1+α\varphi^2$ have emerged as promising candidates, providing an improved fit to the data compared to $Λ$CDM and even the $w_0w_a$CDM parametrization. However, in the absence of a screening mechanism, the regions of parameter space examined exhaustively so far (with $α<0$) violate stringent local constraints on the strength and time variation of gravity and enhance the matter growth rate at low redshifts, thereby compromising the viability of these solutions. In this work, we investigate a largely unexplored region of parameter space (with $α>0$) that realizes the crossing while remaining consistent with local constraints. This setup, which is subject to some fine-tuning, respects the BBN constraints and leads only to deviations of the effective gravitational coupling from Newton's constant below $0.3%$ in the past. We review the background and linear perturbation dynamics of the model, deriving analytical results that provide physical insight into its behavior, and update the constraints on the model using the latest cosmological data, including the DES-Dovekie supernova sample. We perform for the first time a full-fledged Monte Carlo analysis with $α>0$, studying the impact of the variation of $G$ on the supernovae absolute luminosity, and compare the results with those obtained for $α<0$, using in both cases a linear potential. In the context of the model with non-minimal coupling and $α<0$, the standard model is excluded at the $2.50σ$ CL, whereas for $α>0$, the exclusion significance decreases to $\sim 2σ$ CL, indicating a moderate preference for these models over $Λ$CDM. [abridged]

For the first time, we establish a direct connection between inflationary dynamics and the subsequent magnetohydrodynamic evolution. Starting from vacuum fluctuations, we self-consistently evolve the coupled inflaton-plasma system through inflation, reheating, and into the radiation-dominated era. We find an inflaton-driven inverse cascade followed by freely decaying turbulence, with kinetic energy dominating over magnetic energy and driving dynamo amplification. The inflaton-driven turbulence causes steeper magnetic energy decay than helical turbulence and should therefore be accounted for subsequent evolution of primordial fields and their present-day observables.

We provide a unified analysis of scalar-field quantization and detector response in relativistic rotating frames, comparing rigid and Trocheris-Takeno (TT) descriptions in unbounded Minkowski spacetime. We distinguish three questions that are often conflated: (i) how a fixed quantum state is represented in rotating coordinates, (ii) whether rotating time evolution selects a global ground state, and (iii) how a rotating detector responds to the vacuum. Coordinate-transformed rigid and TT modes span the same positive-frequency subspace as inertial modes: their Bogoliubov beta coefficients vanish and their mode sums reproduce the Minkowski Wightman function. Neither rotating time flow selects a global scalar ground state, although for different reasons: the rigid generator is not globally timelike and is unbounded below, whereas TT time translation is not a Killing symmetry. Nevertheless, a circular Unruh-DeWitt detector in the Minkowski vacuum has a stationary, nonzero response. We relate this response to negative corotating frequencies, evaluate it after an explicit Hadamard subtraction, and show that it is nonthermal. Thus vanishing Bogoliubov mixing, the absence of a symmetry-selected rotating vacuum, and detector excitation are mutually consistent.

Substructures appear to be a common trait of all extended protoplanetary discs. They are found notably in discs still surrounded by ambient, possibly infalling material. In this study, we revisit the observational signatures of gravitational instability (GI) in the dust and gas emission of protoplanetary discs. We focus on discs undergoing a short-lived episode of late infall that triggers spirals via GI, in order to investigate the long-term dynamics of millimetre-sized dust. We performed 2D hydrodynamical simulations modelling gas and millimetre-sized dust in a self-gravitating disc, with a simplified treatment of stellar irradiation and infall. Results were post-processed by 3D radiative transfer calculations. During infall, GI sets in and the disc develops a gravitoturbulent state characterised by the continuous formation of spirals. Shortly after infall ends, the millimetre-sized dust remains well coupled to the gas, and the dust emission exhibits large-scale spirals in both near-infrared scattered light and continuum emission. Even long after infall has ceased, two-armed spirals are persistently excited by GI in specific regions of the disc. Their dissipation entails the formation of multiple long-lasting pressure maxima, which can be viewed as vestiges of shocks induced by the spirals. They form persistent dust rings that are more or less axisymmetrical. Consequently, once late infall ends, a massive disc can gradually evolve from a disc with spirals in the continuum into one featuring multiple, potentially lopsided bright rings in the continuum. In near-infrared scattered light, the disc initially displays multiple spirals during infall or shortly after it ceases, before ultimately exhibiting multiple rings. The persistent excitation of spirals implies, however, that the residual line-of-sight velocity, derived from line emission, shows large-scale spirals across the disc.