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

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

Convolutional neural networks are increasingly used to measure source amplitudes in survey maps, often with claims of outperforming the matched filter. That filter is the best linear unbiased estimator (BLUE) for any noise of a given covariance, and minimum-variance unbiased outright when that noise is Gaussian and known, so an advantage requires a covariance that varies from image to image, or noise that is non-Gaussian. We introduce a single number, the advantage ceiling eta >= 1, that quantifies both: the Fisher information for the amplitude in units of the matched filter's, computable from noise-only simulations before any network is trained, and a bound on any marginally unbiased estimator's variance. We compute eta for a taxonomy of millimeter/submillimeter survey noise with a variational score-matching ladder whose rungs are estimator classes of increasing statistical order. With an instrumental white-noise floor, Gaussian noise gives eta = 1 to within +- 0.04; covariance mixtures give exact ceilings of 10.2 (spectral tilt, extended source), 2.1 (leaked components) and 1.1-1.8 otherwise, of which a 2000-parameter mixture matched filter attains 70-81%. Source confusion, whose non-Gaussian component is the noise itself, gives the largest advantage, eta >= 3.8 (extended) and >= 2.4 (compact), reachable only above the bispectrum rung and calibrated at ~ 0.9 against a known answer. Read as observing time, eta multiplies a survey's integration time wherever the noise integrates down, excluding confusion. Unmodelled numerical channels can manufacture spurious advantages of up to two orders of magnitude. A ResNet regressor is 6-15% less efficient than the matched filter on Gaussian noise once prior shrinkage is divided out, and realizes 1.8 of the 10.2 available on the spectral tilt; a companion paper tests such regressors against these ceilings.

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.

We develop an on-shell description of low-frequency gravitational waveforms from black-hole mergers beyond leading order. Treating the strongly coupled merger as effective hard $S$-matrix data, we organize its long-wavelength response using soft theorems and the KMOC formalism. At next-to-leading order, the quantum soft theorem contains logarithmic terms absent from the classical soft theorem. We show that these extra terms cancel in the full KMOC in-in observable between the one-loop radiative amplitude and the corresponding graviton cut, leaving precisely the classical logarithmic contributions associated with gravitational drag and early-time acceleration. We also identify the $1/ω$ corrections from remnant recoil and Christodoulou non-linear memory. These results reveal a hierarchy of merger information accessible at low frequency: logarithmic tails depend only on asymptotic hard data, recoil probes total radiated momentum, while non-linear memory probes the angular distribution of the emitted radiation.

ESA's astrometric space mission Gaia is uniquely positioned for microlensing studies thanks to its all-sky coverage and simultaneous astrometric and photometric observations. In Gaia Data Release 4 (DR4), astrometric time series will be published for the first time, providing a rich dataset for isolated stellar remnant searches. In anticipation of DR4, we prepared tools for fitting photometric and 1D astrometric measurements. We tested the recovery of true event parameters and population distributions with Gaia-like observations. We created a mock survey of microlensing events mirroring Gaia DR4's design. The events were generated using dedicated Galactic simulations to provide realistic expectations on stellar and dark lens yields. For selected events, we simulated observations following Gaia's scanning law and precision, using the astromet package. Such simulations are needed to interpret forthcoming data; they also enable testing inference tools with known ground truth. We used nested sampling to infer full posterior distributions and compare them with the true lens parameters. Our simulation predicts $322 \pm 70$ microlensing events with astrometric signals above Gaia precision thresholds; $209 \pm 57$ of them have detectable photometric signals. Among those, $82 \pm 9$ and $49 \pm 7$, respectively, have remnant lenses. To assess Gaia's capabilities in constraining remnant mass distributions, we modelled $49$ events from this sample. We recovered the true Einstein times and radii for most events. Microlensing parallax measurements are consistent with true values but uncertain, limiting mass determination. We recommend focusing on events with bright sources or anomalous fits to optimize dark lens searches. We make available the mock Gaia time-series dataset, as well as a toolkit repository for working with the data format, fitting models, and visualizing results.

Origin and composition of ultra-high energy cosmic rays are still uncertain, particularly the rare events exceeding the Greisen-Zatsepin-Kuzmin cutoff and whose apparent arrival directions point to the Local Void. This work is an in-depth investigation of the possibility that such cosmic rays contain low-mass magnetic monopoles predicted in several recent theoretical models. Using a custom extension to CRPropa, a state-of-the-art code for modeling cosmic ray propagation, monopoles are tracked in the galactic environment under realistic assumptions on mass, magnetic charge, and initial phase-space distributions. Our simulations demonstrate that relic MMs follow filamentary ``Galactic magnetic funnels'', producing anisotropic arrival directions at Earth. The expected arrival directions are concentrated in a small region of the sky, which has significant implications for the design and interpretation of future searches. A comparison of the simulated arrival directions with the Pierre Auger Observatory dataset shows that the null hypothesis - no monopole contribution - is statistically preferred for the majority of cases. Using a profile-likelihood analysis that incorporates both directional and energy information, we set 90% C.L. upper limits on the integral MM flux of $< 2.1\times10^{-23}$ $\mathrm{cm^{-2}s^{-1}sr^{-1}}$.

The recent 248-keV nuclear recoil candidate reported by LUX-ZEPLIN has renewed interest in endothermic inelastic dark matter, motivating us to examine its capture in primordial dark stars. Relative to the conventional elastic-capture picture in dark star evolution, endothermic capture adds two qualitative features. First, it opens only after the growing star crosses a compactness threshold, expressed as a kinematic radius $R_{\rm kin}\proptoμ_χM_\star/δ$ set by the stellar mass, the dark-matter--nucleus reduced mass and the mass splitting. The accreting growth carries the star through the threshold, and the capture rate turns on quadratically above it. Second, although newly captured particles generically begin on nonthermal bound orbits, endothermic kinematics can keep this normally transient population spatially extended, turning it into a persistent reservoir rather than an intermediate step toward a thermal core. Following complete chains of state-changing collisions, we find that the reservoir compacts sharply and then stalls, because a ground state particle below a compactness-dependent orbital energy has no allowed up-scatter anywhere in the star. A percent-level radius contraction reopens the relaxation. These kinematic results do not depend on whether the excited state decays promptly or is long-lived. As a result, inelastic capture does not replenish a thermal annihilation core. The captured population forms an evolving orbital distribution that sets up the co-evolutionary dynamics between the star and the dark matter in the core and the reservoir, which we develop in a companion paper.

XRISM has mapped gas velocities across the core of the Perseus cluster, separating the kinematic effects of mergers and AGN feedback. The physical properties of these motions remain unclear: are they a superposition of bulk flows, predominantly random/turbulent motions, or a mixture of both? Without resolving this question, constraints on the nonthermal pressure fraction and heating rate remain uncertain, as both assume predominantly random motions. Unlike emission line broadening, resonant scattering is most sensitive to small-scale, random motions rather than coherent bulk flows. Taking advantage of the extensive XRISM coverage of the Perseus cluster, we detect the full effects of resonant scattering on the Heα w line for the first time. This includes flux suppression in the cluster center, enhancement in the outer regions, and non-Gaussianity in the emission line. We employ radiative transfer simulations to constrain the amplitude of small-scale ICM velocities in the inner 60 kpc of Perseus, finding them to be consistent with the line broadening measurements within the uncertainties. This indicates the observed velocity dispersion is primarily due to small-scale random motions in the central Perseus regions rather than coherent bulk flows. We further explore potential anisotropy of these motions, showing that they are consistent with isotropic or radial motions rather than tangential ones. Longer XRISM observations are required to improve these anisotropy constraints. Finally, we explore azimuthal variations between the two complete radial arms observed by XRISM.

We discuss the interpretation of an observed running Hubble constant with the redshift, in terms of the possible underlying physical scenario. Assuming, like in the $Λ$CDM model, that the matter and dark energy components do not directly interact, we arrive at the conclusion that the physical content of the cosmological dynamics can always be represented with an evolutionary dark energy paradigm. We discuss both the case of the power-law running Hubble constant, and that of an optimized fitting model, both studied with the binned Pantheon Sample of the Type Ia Supernovae data. These two cases turn out to be associated with a phantom-like evolutionary dark energy.

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.

An eV-scale QCD axion comprising the observed dark matter (DM) abundance is expected to generate a photon line at infrared energies that would be observable or near-observable in data collected by the James Webb Space Telescope (JWST), as might more general axion-like particles (ALPs) over a broader range of masses and couplings. This has motivated a number of efforts to either forecast JWST sensitivities to a QCD axion or realize them through analyses of publicly available datasets. At present, no consensus exists; leading analyses disagree by as much as an order of magnitude in terms of axion-coupling sensitivity, implying an orders-of-magnitude discrepancy in raw flux density sensitivity, and consistency between the analyses and prior forecasts is unclear. We address these outstanding discrepancies with a bespoke data reduction and flexible nonparametric inference procedure that lead to well-controlled and robust limits on the decay of eV-scale axion DM, consistent with previously forecasted sensitivities. We further demonstrate that the strongest previously claimed sensitivities exceed those attainable by any analysis of the datasets from which they were derived. We exclude QCD axion DM for masses between $500\,\mathrm{meV}$ and $2.5\,\mathrm{eV}$ using NIRSpec data, while setting limits on ALP DM complementary to other astrophysical constraints at masses between $100\,\mathrm{meV}$ and $500\,\mathrm{meV}$. However, we find the sensitivities to be systematically limited, and therefore unlikely to be improved upon by ongoing data collection or re-analysis unless instrumental modeling and data reduction pipelines improve considerably.

Cosmic rays can significantly alter the abundances of certain species in the upper layers of planetary atmospheres, especially in terms of their biosignatures. To fully understand the extent of this effect, it is essential to accurately model the interactions of cosmic rays with planetary magnetic fields. We used the CosmicTransmutation code to study the effect of a planetary magnetic field on the flux of galactic cosmic rays and stellar energetic particles. We found that both particle sources are significantly affected by magnetic fields, even though the effects are different due to the varying energy ranges that characterize each source. The stellar energetic particle energy flux is significantly higher for an Earth-like planet with no magnetic field, but with a magnetic field of 30\,$μ$T or higher, the energy flux of the two sources becomes comparable. We find that the atmospheric pressures the cosmic rays reach are significantly lower than those found in previous studies using simpler models and that the effect of the atmospheric composition on this outcome is small. We found similar results using the model for the exoplanet K2-18b, even though the radius is significantly larger than Earth's. Because the two cosmic ray sources cover different energy ranges and are affected by the magnetic field in different ways as a result, both sources should be considered. Future studies should focus on combining an accurate modelling procedure of the interaction between cosmic rays and planetary magnetic fields with atmospheric chemistry models. This will require general circulation models to capture the latitudinal and longitudinal dependencies in full.

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) is expected to significantly advance the study of Milky Way stellar streams. In particular, the deep, precise photometry from LSST should greatly increase the statistical sensitivity to density fluctuations in stellar streams, which can be used to probe the small-scale distribution of dark matter. However, current forecasts generally neglect the impact of observational systematics that will be imprinted on stream density measurements. In this study, we develop a realistic forward-modeling framework to inject stellar streams into LSST-like observations including photometric uncertainties, survey depth variations, background contamination, and imperfect star-galaxy classification. We develop a likelihood-ratio analysis to assess the detectability of gaps in stellar streams in the presence of these observational systematics. In the presence of realistic survey systematics, we find that after four years of operations, LSST will be sensitive to density reductions of $\sim50\%$ for gaps with widths of $5$ deg in streams with surface brightness of $\sim33$ mag arcsec$^{-2}$. Relative to the ideal case, this corresponds to a degradation in gap depth sensitivity by a factor of $\sim5$ due to the combined impact of background contamination and observational systematics. Assuming a simplified analytical mapping between gap depth and dark matter subhalo properties, these estimates correspond to a minimum detectable subhalo mass of $\sim1\times10^7$ M$_\odot$. Observational effects shift this accessible mass scale upward by a factor of $\sim16$, with background contamination contributing a factor of $\sim5$ and survey systematics a further factor of $\sim3$, dominated by star-galaxy classification.

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.

We consider 4-dimensional asymptotically flat vacuum spacetimes near future null infinity endowed with the most general allowable Carroll geometry. We show that the near-boundary radial expansion (to a certain order) can be organised in terms of connections that can be obtained by gauging the conformal Carroll algebra. The only non-vanishing curvatures in this gauging procedure are those that are associated with the special conformal generators and we will refer to these as the $K$-curvatures. The vanishing of these $K$-curvatures defines an asymptotic vacuum spacetime and we use this to construct a boundary (i.e. Carroll) covariant expression for the vacuum (soft) shear in terms of two boundary Carroll scalar fields. The $K$-curvatures transform in a hierarchical fashion into one another under Carroll boosts. This leads to a classification of 4 types of spacetimes: vacuum, strongly and weakly non-radiative, and radiative spacetimes. It is shown that the $K$-curvatures correspond to 5 of the 10 Weyl tensor components at leading order in their $1/r$ expansion. We furthermore observe that one of the $K$-curvatures is equal to the recently found Carroll boost anomaly. Finally, we show that the Bondi loss equations for the boundary energy-momentum-news complex can be cast into a form involving another energy-momentum tensor with vanishing energy flux that is traceless and whose non-conservation is entirely captured by the $K$-curvatures. The BMS currents can be obtained by contracting this latter energy-momentum tensor with a Carroll conformal Killing vector.

Violations of Bell's inequalities offer a definitive signal of non-classical (quantum) behavior in local deterministic systems. Yet, in many physical settings where quantum mechanics is expected to play an important role, one cannot construct a Bell-type test using the available observables. Cosmology offers one concrete example, where cosmic structure may have originated from quantum vacuum fluctuations yet all observations are effectively classical. Nevertheless, recent work suggests that quantum and classical time evolution may still be differentiated by the pattern of statistical correlations of these observables. In this paper, we explore and generalize the use of correlations as Bell-type test of the nature of time evolution. We show that quantum vacuum fluctuations of closed systems undergoing Hamiltonian evolution produce unique correlations that are not mimicked by classical Hamiltonian systems. Correlations in the quantum vacuum are generally time-independent and are controlled by the energy gap between the ground and exited states. Classical systems exhibit apparent poles at physical frequencies that do not arise in the quantum vacuum. As one approaches the poles, the evolution becomes dominated by a resonant Hamiltonian giving rise to time dependent correlations that eventually decay through dephasing. Finally, we show that this behavior is distinct from other signals of quantum evolution, including applications to quantum optics, quantum walks, and quantum search.

We investigate the influence of the Garfinkle-Horowitz-Strominger (GHS) dilaton black hole on different quantum resources of Dirac fields beyond the single-mode approximation. By employing the negativity to characterize quantum entanglement and the $l_1$-norm and the relative entropy of coherence to characterize quantum coherence, we demonstrate that these resources exhibit remarkably different responses to the gravitational field. Specifically, increasing the dilaton parameter continuously suppresses quantum entanglement, leaving only a finite residual amount in the strong-gravity regime, whereas quantum coherence is enhanced, indicating that the dilaton-induced spacetime affects nonlocal quantum correlations and local quantum superposition in fundamentally different ways. Furthermore, we show that an initially maximally entangled state does not always possess the largest negativity after propagating in the GHS dilaton spacetime; instead, under appropriate conditions, certain non-maximally entangled states can retain stronger entanglement than the maximally entangled one. These findings reveal the resource-dependent nature of gravitational effects in dilaton black hole backgrounds and provide new insights into the manipulation and protection of quantum resources for relativistic quantum information processing in curved spacetime.

We extend the Dark Matter-Baryon Separability Condition to show that the same framework developed for dark matter deficient galaxies naturally admits a conjugate branch describing preferential baryonic depletion. Using the recently discovered almost-dark galaxy TTT J1237327+143535 as a worked example, we derive a family of dynamical consistency relations parameterized by the unknown progenitor ratio $μ_i$, including thresholds for $σ_{\rm los}$, $M_{\rm dyn}$, the enclosed baryonic fraction and the dynamical mass to light ratio. We further connect the separability framework to the expected globular cluster population, providing an independent consistency test of the inferred halo mass. We also provide some bounds on the baryon ratio using the tidal properties of Virgo Cluster. Future measurements of stellar kinematics, gas content and globular clusters can therefore determine whether this galaxy occupies the positive separability branch and also test whether dark matter deficient and baryon depleted systems can be described within a common framework.

We present an emulator of the halo mass function (HMF) based on Gaussian-process regression, that exploits the evolution mapping framework. This framework splits cosmological parameters into shape parameters ($\mathbfΘ_\mathrm{s}$), which set the shape of the linear power spectrum, and evolution parameters ($\mathbfΘ_\mathrm{e}$), which only rescale its amplitude, with their net effect entirely captured by the clustering amplitude, $σ_{12}$. The emulation acts in two steps: first, an emulator predicts the overall shape of the HMF given $\mathbfΘ_\mathrm{s}$ and $σ_{12}$; then, a second emulator applies a small correction accounting for the residual dependence on the integrated growth history of density fluctuations. The emulator is trained on a suite of 100 simulations of flat $Λ$CDM cosmologies but, owing to evolution mapping, can be applied to predict the abundance of haloes in cosmologies with any extended space of $\mathbfΘ_\mathrm{e}$, such as curvature, quintessence and dynamical dark energy. We test the emulator against independent $Λ$CDM runs and on an additional validation suite of dynamical dark energy simulations, achieving per cent-level accuracy over a wide range of masses and redshifts, with the performance only slightly degrading at early times, where the training data points are fewer. Compared to other state-of-the-art HMF emulators, its accuracy is higher in nearly all tested cases. Our results demonstrate the broad applicability of the evolution mapping framework to build very accurate emulators of statistics of the density field. The emulator is publicly available as part of the Aletheia emulator suite.

Ultralight dark matter could be an axial vector field, a simple possibility which motivates new experiments. Couplings to fermions, through an axial vector current or a dark electric dipole moment operator, lead to enhanced spin torques compared to axion dark matter. An axial vector also has an analogue of the axion-photon coupling, though its consistent realization requires a photon mass. Under this "axial-photon" coupling, the effect of a background magnetic field is suppressed, and the strongest experimental probes involve polarimetry, electric fields in superconducting cavities, and the cosmic microwave background. Since a massive axial vector is dual to a massive two-form, our results also apply to "Kalb-Ramond" dark matter motivated by string theory.

In a broad class of inflationary models, the leading non-Gaussian bispectrum arises from one-loop rather than tree-level processes. However, realistic one-loop contributions remain largely unexplored for phenomenologically relevant masses beyond the simplest bubble diagrams. This is insufficient for observational purposes because, for generic masses and couplings, triangle contributions need not be suppressed relative to their bubble counterparts. We present a numerical method for evaluating scalar one-loop in-in diagrams contributing to the inflationary bispectrum at fully general external momenta and masses. First, the Witten-Feynman parameterization reduces de Sitter loop integrals to generalized Euler-Mellin integrals governed by Symanzik graph polynomials. Unfortunately, the exponents of the polynomials in the integrand become complex for sufficiently heavy masses, leading to sign problems when evaluating the integral using standard numerical techniques. To remedy this problem, we introduce the reduced Schwinger method, which evaluates a highly oscillatory subintegral analytically to obtain a Gauss hypergeometric kernel, leaving the remaining integrals to standard numerical quadrature. We validate the algorithm by reproducing known analytic results for the tree-level bispectrum in terms of ${}_3F_2$ functions, and then apply it to the one-loop bubble and triangle contributions at general kinematics. This yields the first direct numerical evaluation of the complete one-loop bispectrum, valid across the full kinematic range. Using these numerical results, we construct bispectrum templates and compare them with CMB data.

We study the electromagnetic radiation and wave-optical imaging of an oscillating electric dipole orbiting a Kerr black hole. We derive the effective 4-current associated with a pointlike oscillating electric dipole in curved spacetime, and use black hole perturbation theory to compute the resulting radiation field at future null infinity, from first principles. We then develop a wave-optical imaging framework for a moving electromagnetic source in curved spacetime. We obtain images of an orbiting electric dipole, displaying relativistic beaming, gravitational lensing, and Einstein rings. We study polarization-dependent scattering by comparing the images produced by spinning dipoles with opposite helicities, finding a displacement that roughly decreases with the inverse of the radiation frequency, as expected for a beyond-geometric-optics effect. Our results provide a first-principles benchmark for beyond-geometric-optics descriptions of electromagnetic radiation in Kerr spacetime.

The Sommerfeld enhancement (SE) from long-range interactions, and the related bound state formation rate, can be important non-perturbative inputs to dark matter (DM) annihilation signals. In the presence of singular interaction potentials, the conventional boundary conditions of the Schrödinger equation fail, obfuscating the computation of SE, its physical origin, and its relation to bound states in such potentials. In this work, we clarify the origin of SE in singular potentials in a two-fold manner: using the framework of velocity power counting in non-relativistic effective field theory (NREFT), and via position-space regularization of singular potentials at short distances to compute SE and bound states. We illustrate our findings through the case of pseudoscalar mediators interacting with massive Dirac DM. We find that when such a system arises from a UV-complete theory at weak coupling, no SE is generated. However, in the case of a derivatively coupled pseudoscalar, where the interaction is described by a higher-dimension operator in an effective theory, SE (and bound states) can occur for weak couplings if there is a large hierarchy between the cutoff scale of the effective theory and the dark matter mass. This SE is sensitive to the choice of the UV completion of the potential at short distances, but a non-negligible SE can persist even when the UV physics alone would not generate any SE; we elucidate the interplay of UV and IR physics in this case.

The growth of primordial black holes through the absorption of the surrounding plasma has recently been shown to exhibit a critical behavior during radiation domination. We extend this analysis to the reheating era and derive analytical solutions for general reheating histories. We show that reheating modifies the critical condition for runaway absorption, making it dependent on both the reheating dynamics and the black-hole formation time. We identify two distinct regimes: runaway growth occurring during reheating, or being triggered after the onset of radiation domination by the mass accumulated during reheating. More generally, we derive a simple composition law describing how independent mass- growth mechanisms combine across successive cosmological eras. Applying it to radiation absorption and inflaton accretion, we obtain analytical results in excellent agreement with the full numerical evolution.

We investigate the internal structure of two galaxies from the LEWIS sample: the ultra-diffuse galaxy UDG-1 and the extended dwarf LSB-6. Both show coherent stellar rotation combined with a non-negligible fraction of random motion, with no signs of ongoing disturbance. This is the first attempt to constrain the dynamics and dark matter physics of rotation-supported UDGs through integral-field stellar kinematics. We model the galaxies as stellar spheroids with typical dwarf-like thickness embedded in spherical dark matter halos, deriving two-dimensional velocity fields that we compare with the observations in a Bayesian framework. Besides cold dark matter, we test fuzzy, self-interacting, and non-minimally coupled dark matter. The data prefer a cuspy halo in UDG-1 and a cuspless one in LSB-6. All the alternative models remain viable, although current data cannot conclusively discriminate among them. For fuzzy dark matter we obtain mutually consistent boson masses, $m_α= 4.7^{+1.8}_{-1.1} \times 10^{-22}$ eV for UDG-1 and $m_α= 6.3^{+0.4}_{-0.4} \times 10^{-22}$ eV for LSB-6. For self-interacting dark matter, LSB-6 yields a robust constraint on the velocity-weighted cross section, $\langle σv\rangle/m = 14.3^{+2.0}_{-1.9}$ cm$^2$ km g$^{-1}$ s$^{-1}$. For non-minimally coupled dark matter we derive upper bounds on the coupling length, implying marginal deviations from $Λ$CDM. Regardless of the model, both galaxies exhibit a dark matter content comparable with halos of typical dwarf galaxies with similar stellar masses. Together with stellar population analyses and globular cluster information, these results support a scenario in which UDG-1 and LSB-6 originate from puffed-up dwarfs. This work is the first in a series exploiting kinematic data of extremely low-surface-brightness galaxies to test fundamental physics, from dark matter to modified theories of gravity.

Baryonic feedback is a dominant source of systematic uncertainty for upcoming weak-lensing surveys, but current tools for modeling its effect rely on spherical approximations and density profiles calibrated almost entirely on two-point statistics. We introduce BIND (Baryonic INpainting with Deep learning), a conditional flow-matching model that learns a field-level mapping from dark-matter-only halos to their hydrodynamical counterparts. BIND is trained on halos from the 1024 paired hydrodynamical and dark-matter-only simulations of the CAMELS $50\,h^{-1}\,\mathrm{Mpc}$ SB35 suite and samples dark matter, gas, and stellar mass fields over redshift across the full 35-dimensional $Λ$CDM and IllustrisTNG galaxy formation parameter space. BIND recovers dark matter, gas, and stellar masses at the percent level, reproduces azimuthally averaged profiles to $\lesssim10\%$ at all radii, and matches halo shape distributions with high fidelity. The learned parameter dependence captures the rank correlations between the generated fields and the subgrid parameters, and the field-level response to individual parameter variations is recovered in both sign and morphology. Halo mass is never supplied as conditioning, yet the baryon fraction, stellar-to-halo mass relation, inter-component scaling relations, and the joint covariance of their residuals are all reproduced. We finally show that, applied halo-by-halo to a $(50\,h^{-1}\,\mathrm{Mpc})^3$ $N$-body volume with $512^3$ particles, BIND reproduces the projected matter power spectrum suppression to the accuracy ceiling set by pasting in the hydrodynamical halos themselves, in minutes on one GPU. We release the trained BIND models and all generated halos as open-source tools. A companion paper extends BIND to thermodynamic fields and non-Gaussian weak-lensing statistics.

Recent multiwavelength observations of galaxy group and cluster gas suggest stronger baryonic feedback than our best-calibrated hydrodynamical simulations produce, while modeling this feedback remains a primary source of uncertainty in Stage-IV weak-lensing (WL) analyses. We present a suite of ray-traced maps generated with BIND (Baryonic INpainting with Deep learning), a conditional flow-matching model that paints baryonic mass and gas thermodynamics onto the halos of dark-matter-only simulations, and which was developed in a companion paper. Applied to IllustrisTNG300-Dark and ray-traced, we generate convergence, optical depth, and Compton-$y$ maps at five source redshifts, each with $1000$ pseudo-independent realizations. We build lightcones across a 256-node Sobol sequence spanning the thirty-dimensional IllustrisTNG galaxy formation prior, with individual parameter variations and at the fiducial model. In validation, the maps match those built from the IllustrisTNG300 halos to within LSST-Y10-like precision for a range of WL statistics. Across the prior, the response to feedback exceeds Stage-IV statistical precision by more than an order of magnitude on small scales, and different statistics respond to different model sectors: galactic winds control the WL power spectrum and gas auto- and cross-spectra, while the stellar initial mass function slope and AGN parameters shape the morphological statistics (PDF, peaks, minima, and Minkowski functionals). Finally, we find that the response of statistics can be compressed into seven halo properties which linearly predict a range of WL and SZ statistics. We publicly release the maps, statistics, and model tables.

We present GW Explorer: A Beginner's Guide, an outreach curriculum designed to introduce high school students to gravitational-wave (GW) astrophysics through interactive Python Jupyter notebooks. Most existing GW resources target beginning audiences and advanced students, leaving a gap at the pre-college level that we directly address. The curriculum integrates foundational physics with hands-on computation implemented through both self-directed and workshop-based instructional formats. In the self-directed format, students completed the curriculum independently on cloud-based platforms such as Google Colab. During the workshop format, students worked through the same activities under the guidance of University of Nevada, Las Vegas graduate student mentors. Topics span gravity, spacetime, GW sources, interferometric detection, and data analysis. An implementation in local high school classrooms informed the content and pacing, and survey results demonstrate gains in conceptual understanding and coding confidence. GW Explorer offers a scalable, open-access framework for authentic astrophysics research in the high school classroom.

We investigate a quantum massive scalar field in the interior of a charged and spherically-symmetric (Reissner-Nordström) black hole. We examine the behaviour for varying values of the black hole charge, field mass and coupling constant when the field is in two quantum states: Hartle-Hawking (representing a black hole in thermal equilibrium) and Unruh (representing a black hole evaporating via the emission of Hawking radiation). We show that the vacuum polarization as well as the angular components of the quantum stress-energy tensor diverge on the Cauchy horizon, in stark contrast to what happens for massless fields. We also calculate the energy fluxes in Eddington-Finkelstein coordinates $\{u,v\}$. We show that these fluxes on the Cauchy horizon do not generically vanish. This implies, in particular, that in regular, Kruskal coordinates $\{U,V\}$ the ingoing flux diverges like $V^{-2}$ on the Cauchy horizon (where $V=0$). This divergence suggests that its backreaction via the semiclassical Einstein equations would yield a strong singularity, as opposed to its weaker, classical counterpart. Interestingly, there are exceptions, in which the energy fluxes in Eddington-Finkelstein coordinates vanish: (i) in the extremal limit (where the black hole is maximally charged); (ii) certain fine-tuned regions of parameter space, where the fluxes change sign.

Bimetric theory is an extension of general relativity that perturbatively describes one massless and one massive graviton. It admits observationally viable cosmologies, passes local tests of gravity, and provides candidates for dynamical dark energy and spin-2 dark matter. Scattering amplitudes provide a complementary probe of its consistency, e.g. through analyticity, unitarity, and causality. While the tree-level $2\to2$ amplitudes of the closely related theory of massive gravity have previously been computed and analysed, the richer amplitude structure of bimetric theory has mostly been studied for selected processes and helicity sectors. We present the complete set of tree-level $2\to2$ scattering amplitudes of ghost-free bimetric theory around proportional Minkowski backgrounds. We expand the action through quartic order in the mass eigenstates and use a computer-algebra workflow to compute all of the $199$ symmetry-inequivalent helicity amplitudes. We find a clear hierarchy: amplitudes with a single massive external state vanish, those with exactly two are independent of the nonlinear bimetric parameters, and nonlinear ghost-free parameter dependence first appears with three massive external states. Taking the massive-gravity limit yields the complete set of tree-level $2\to2$ massive gravity amplitudes, which agree exactly with previous results. Finally, all amplitudes grow with energy at most as $E^6$, corresponding to the characteristic $Λ_3$ strong-coupling scale, and this maximal growth cannot be eliminated by any nontrivial choice of the theory parameters.

We examine an on-disk solar quiet region of enhanced magnetic network using Hbeta images from Daniel K. Inouye Solar Telescopes (DKISTs) Visible Broadband Imager (VBI) and line-of-sight magnetograms from DKISTs Visible Spectro-Polarimeter (ViSP). We also compare the ViSP magnetograms with co-aligned co-temporal line-of-sight magnetograms from the Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI). We find:(i) Two types of chromospheric jet-like features: Hbeta spicules, and small-scale jets. Both are rooted near edges of magnetic network lanes. Several sit close to ViSP-detected tiny islands of either minority-polarity flux or dips in majority-polarity flux. Hence, those several plausibly stem from mixed-polarity magnetic flux, and are caused by chromospheric magnetic reconnection. For H\b{eta} spicules, the average width, length, lifetimes, and speeds are 420$\pm$400 km, 2600$\pm$1600 km, 4.3$\pm$0.25 min, and 12$\pm$4.7 kms. For small-scale jets, those are 385$\pm$100 km, 620$\pm$40 km, 5$\pm$3min, and 4.4$\pm$1.5 kms. (ii) H$β$ microflashes sit in evidently unipolar flux. Their widths, lengths, lifetimes, and speeds are 180$\pm$65 km, 365$\pm$100 km, 4.3$\pm$2.6 min, and 3.1$\pm$0.2 kms. (iii) At the base of a coronal plume, due to its higher spatial resolution than HMI, ViSP shows both some very strong (>800 G) majority-polarity flux not shown by HMI and some sub-arc-second minority-polarity-flux inclusions not shown by HMI. These sub-arcsecond chromospheric transients carry sufficient energy (10^{24}erg) to transiently heat the local chromosphere and corona and potentially contribute to solar-wind acceleration. They may represent the small-scale end of a continuum of magnetic-reconnection-driven activity, highlighting the importance of DKIST's high-resolution magnetic-field measurements for understanding small-scale chromospheric dynamics.

We present ultra-precise fundamental parameters for the detached eclipsing binary UX Men based on 31 TESS sectors ($\sim$882 d) and radial velocities from HARPS, FEROS, CORALIE, and UCLES. The light curves were detrended for spot-induced variability and modelled independently with JKTEBOP and ALLESFITTER, while radial velocities were extracted via TODCOR and fitted with a Keplerian orbit. Spectral disentangling yielded individual component spectra analysed with iSpec to determine $T_{\rm eff}$ and $\mathrm{[M/H]}$. Multi-band photometry was used to determine the fluxes, flux ratios and ultimately to measure photometric $T_{\rm eff}$. We obtain the masses $M_{\rm A}$ = $1.2300 \pm{0.0012}$ M$_{\odot}$, $M_{\rm B}$ = $1.1946\pm{0.0012}$ M$_{\odot}$ and radii $R_{\rm A}$ = $1.3605\pm{0.0031}$ R$_{\odot}$, $R_{\rm B}$ = $1.2801\pm{0.0069}$ R$_{\odot}$ with sub-percent precision. Uncertainties were found to be dominated by systematic effects coming from (weak) stellar activity. $T_{\rm eff,A}$ = $6302\pm{77}$, $T_{\rm eff,B}$ = $6292\pm65$ \& Metallicity estimates are $[M/H]_{\rm A} = 0.02\pm{0.13}$ dex and $[M/H]_{\rm B}= 0.07\pm{0.10}$ dex. Isochrone fitting yields a consistent age of $τ$ = $2.75\pm{0.13}$ Gyr. The system exhibits a $\sim$ 673-day periodicity in eclipse-timing variations, likely due to stellar activity, with no evidence for a third body. Results are independent of the modelling code used, but without proper treatment of activity-induced systematics the uncertainties would be underestimated. UX Men is a high-precision benchmark binary providing stringent tests of stellar models and a robust calibration source for PLATO. Expanding the sample of well-characterised detached eclipsing binaries will strengthen the empirical foundations for stellar physics and the calibration of future photometric missions

We carry out a search for the relative baryon-cold dark matter (CDM) density and velocity modes in the DESI full-shape (FS) data, and investigate whether DESI's preference for evolving dark energy can be affected by these relative perturbations, which are a known contaminant of the baryon acoustic oscillations. For the first time, we use a holistic effective field theory (EFT) treatment of the power spectrum and bispectrum both for the standard terms (without relative perturbations), modeled in EFT to one-loop order, and for the relative perturbations. For the latter we implement an efficient computation of the one-loop power spectrum that allows us to scan over different cosmologies. We obtain robust constraints on the relative density and velocity bias parameters $b_{δ_{bc}}$, and $b_{v^2_{bc}}$ for all the DESI DR1 tracers, finding e.g. $b_{δ_{bc}}=5 \pm 6$, $b_{v^2_{bc}}=-0.031\pm 0.023$ (at 68% CL) for the LRG2 sample in a DESI FS-only analysis within $Λ$CDM, supplemented by CMB priors on the baryon density and spectral tilt. Assuming instead a $w_0w_a$CDM cosmological model, and combining DESI FS with the CMB primary anisotropies, CMB lensing, and Pantheon+ supernovae, our LRG2 constraints read $b_{δ_{bc}}=2 \pm 5$, $b_{v^2_{bc}}=-0.032\pm 0.022$. Adding the relative velocity divergence biases we find $b_{δ_{bc}}=-17_{-17}^{+19}$, $b_{θ_{bc}}=-760_{-680}^{+700}$, $b_{v^2_{bc}}=-0.032\pm 0.024$ for the same sample. The constraints on the dark energy equation of state parameters are virtually unchanged when marginalizing over the baryon-CDM modes. We conclude that DESI's preference for evolving dark energy is robust with respect to the baryon-CDM relative perturbations.

We develop a largely automatic semi-analytic method for calculating weakly damped quasi-bound states and trapped modes supported by a local minimum of an effective potential. The real part of the frequency is obtained from an arbitrarily high-order local WKB expansion and Padé resummation, while the exponentially small imaginary part is estimated using the Gamow approximation for tunnelling through one or two potential barriers. Because the local quantization requires only derivatives of the potential at its minimum, the method applies readily to non-rational effective potentials and to different compact-object geometries. We provide a Mathematica notebook implementing the procedure. Comparisons with continued-fraction and direct numerical results for massive scalar fields around Schwarzschild and Kerr black holes and for axial trapped modes of a uniform-density star show that the method yields accurate real frequencies and useful estimates of decay rates outside the superradiant regime.

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.

We derive an entropy formula for recently proposed tensor network models which prepare diffeomorphism invariant states of topological quantum field theories with non-compact and/or continuous gauge groups. We show that our entropy formula generalizes the notion of ``topological entanglement entropy'' to incorporate the infinite number of particle-like excitations in such theories. When our networks are endowed with gauge group $\mathrm{SL}(2,\mathbb{R})$, we can interpret them as models of three-dimensional gravity with small Newton's constant and possibly non-invertible metrics.

We associate a gravitational Hilbert space $\mathbf{H}_σ$ to any closed compact $(d-1)$-manifold $σ$ with real metric. A quantum state $\mathcal{J}(σ)$ is a $d$-manifold bounded by $σ$ and equipped with elliptic data. An inner product is defined by gluing states pairwise across $σ$ and evaluated by viewing the resulting closed $d$-manifold as a boundary condition on the gravitational path integral (GPI) over $(d+1)$-manifolds. If $σ$ is nonempty and the GPI is dominated by a single $(d+1)$-manifold $M$ in the $G_N\to 0$ limit, then $M$ contains a Lorentzian CRT fixed-point set, providing $\mathcal{J}(σ)$ with a classical spacetime interpretation. Conversely, given a finite Lorentzian domain with edge $σ$, a state $\mathcal{J}(σ)$ may be associated to it by deforming its initial data off the real Lorentzian section and retaining only elliptic data. This establishes a broad correspondence between non-asymptotic spacetime regions and quantum states. Assuming that $\mathbf{H}_σ$ factorizes over connected components of $σ$, our framework admits operators and partial traces. This allows us to explore the information-theoretic structure of the states we define. As an example, we construct a family of states by deforming partial Cauchy slices $Σ$ that straddle a two-sided black hole; $σ$ consists of two spheres. We construct the reduced state on one sphere and find that its Rényi entropies are positive, monotonic, and sensitive to all aspects of $Σ$ and its complex deformation. The von Neumann entropy, however, is controlled only by the maximin surface in the causal domain of $Σ$, independently of other parameters, so long as the complex deformation does not vanish. Our proposal may thus explain the efficacy of tensor network toy models of holography while transcending their limitations.

Hot Jupiters with orbital periods shorter than a few days have probably been affected by tidal orbital migration. We develop an analytical framework for constraining tidal migration from the present-day hot Jupiter period distribution, taking into account the uncertain rate and period distribution of hot Jupiters produced by mechanisms such as high-eccentricity migration or disk-driven migration. Assuming the tidal migration timescale is proportional to $P^{χ_τ}$, solutions with $χ_τ\simeq 3, 1.7,$ and 5.6 are all compatible with the present-day period distribution. The $χ_τ\simeq 3$ solution is consistent with equilibrium tides with suppression of dissipation at short periods, and implies that newly circularized hot Jupiters have periods concentrated near $3-4$ days, as predicted in some high-eccentricity migration models. The $χ_τ\simeq 1.7$ solution is also compatible with the $3-4$ day peak but has no clear counterpart in existing tidal theories and is more finely tuned. The $χ_τ\simeq 5.6$ solution is compatible with enhanced short-period equilibrium tidal dissipation or weakly nonlinear gravity-wave dissipation, but requires circularization at unexpectedly short periods. Thus, we find the model with $χ_τ\simeq 3$ most appealing. Transit timing of individual systems and observational constraints on the rate of hot Jupiter engulfment provide additional constraints, which are presently inconclusive but should improve with future data. Improved measurements of the occurrence of short-period planets as a function of planet mass and system age could also help to sharpen the constraints on tidal migration.

As parameters are varied, the set of saddle points contributing to a (path) integral may change discontinuously, leading to a corresponding change in the asymptotic expansion of the integral. This behavior is known as the Stokes phenomenon. We explore this phenomenon in the context of the analytic continuation problem relating the AdS/CFT and dS/CFT correspondences. In this paper, we study these correspondences for three-dimensional pure gravity and two-dimensional Liouville theory, using independent calculations in bulk minisuperspace and in the boundary Liouville zero-mode. In the bulk, the dS contour selects a single saddle and yields the tunneling wave function. Upon continuation to AdS, the contour instead selects an infinite family of saddles. The boundary calculation independently reproduces the same Stokes structure at leading semiclassical order, providing a nontrivial holographic consistency check. Our construction also offers a contour prescription for the conformal factor problem in Euclidean AdS$_3$ quantum gravity within minisuperspace.

We investigate the fate of asymptotic particle states in supersymmetric quantum field theories on de Sitter backgrounds, distinguishing algebraic field content from localized LSZ-like excitations. We study how long-wavelength gravitational fluctuations may suppress pole residues and redistribute spectral weight from isolated contributions toward continuum sectors. The analysis is formulated within a Kallen-Lehmann spectral framework and is intended as a structural infrared mechanism rather than a non-perturbative proof of de Sitter quantum gravity. The paper develops three connected levels. First, we analyze residue suppression in a quasi-de Sitter setting and its possible transmission through Yukawa and gravitational interactions. Second, we formulate the spectral decomposition as a map between the full field-theoretic spectrum and the observable particle sector. Third, we examine a conditional cosmological realization of redistributed spectral weight, including possible effects on expansion, structure growth, and lensing. This cosmological layer is a downstream phenomenological test and is not used as evidence that pole loss has occurred. In the current MCMC analysis we find S8 approximately 0.803 and, when the smooth fraction is varied, fS8 = 0.0025 +/- 0.0017, consistent with the benchmark value 0.00262. These results test the internal consistency of the proposed realization but do not establish the dark-sector assignment or the underlying infrared mechanism.

The scatter in elemental abundances among stars of similar age and metallicity reflects chemical inhomogeneity in their birth environments, making abundance scatter a key observable for chemical-tagging studies of Galactic formation and evolution. Using a large sample of subgiant stars with precise ages and elemental abundances derived from LAMOST low-resolution spectra, we investigate the intrinsic chemical abundance scatter of the low-$α$ thin disk near the solar neighborhood ($7 < R < 10$ kpc). We model the abundance ratio [X/Fe], for each of the 19 elements of concern, as a function of age, [Fe/H], and [Mg/Fe], and deduce the intrinsic dispersions with a forward modelling technique. Our results confirm previous findings that the intrinsic scatters are small, typically $\lesssim0.05$~dex, for light elements (C, Al), $α$-elements (O, Mg, Si, Ca, Ti), and iron-peak elements (Mn, Ni). A dedicated analysis of M67 yields similarly small scatter values for these elements, implying limited discriminatory power from light-element abundances alone. In contrast, neutron-capture elements exhibit substantially larger scatters, typically $\gtrsim$0.1 dex, which are significantly larger than those of M67 member stars ($\sim$0.07~dex). In particular, our analysis suggests that the abundance variations of individual neutron-capture elements cannot be explained by a single tracer such as [Ba/Fe]. These findings clarify the utility of neutron-capture elements for chemical tagging and highlight the potential of low-resolution spectroscopy in such studies.

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.