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

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.

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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