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

The extragalactic background light (EBL), ranging from the infrared to the ultraviolet bands, is the second most intense photon field in the universe, surpassed only by the cosmic microwave background (CMB). It is primarily generated by starlight in galaxies, either directly or through absorption by dust and re-emission at longer wavelengths. Very high energy (VHE, $E$ > $100\,$GeV) photons can be absorbed via $γγ$ interactions with the EBL during their propagation across cosmological distances, providing an indirect method to probe the EBL by studying its impact on the energy spectra of distant VHE sources. This paper examines the robustness of EBL constraints derived from gamma-ray data, critically assessing the assumptions made in previous studies about the intrinsic source spectra, the uncertainties of the observations and the validity of the analysis tools. We find that earlier studies likely underestimated the uncertainties in the EBL intensity constraints, underscoring the need to account for systematic uncertainties comprehensively. By employing a Monte Carlo (MC) simulation and a plausible model for unknown systematic errors, we compute more realistic uncertainties. Additionally, we discuss possible alternatives to set EBL constraints, relaxing the assumptions on the intrinsic spectra of gamma-ray sources, with the goal of achieving more robust constraints.

Neutron stars are among the densest objects in the universe, and the mass measurements are crucial to the understanding of their equation of state. Pulsars in tight binaries can be weighed dynamically, but accurate mass estimates are often hampered by uncertainties in orbital inclinations. In this paper, we report on multi-epoch optical monitoring of the pulsar binary, PSR J1048+2339, using the 0.5-m RIFT telescope at Lulin Observatory. Our observations reveal frequent optical flares with eclipses caused by the companion. Using the eclipse duration with archival photometric data and orbital ephemeris in the literature, we constrain the inclination to 78.0+/-1.0 degrees and the pulsar mass to 1.79 (+0.09, -0.08) solar masses. This measurement is among the most precise for any redback system to date, demonstrating a novel way to improve the pulsar mass measures with small aperture optical telescopes. Furthermore, the optical eclipse duration is significantly shorter than that measured in M/GeV gamma-rays, possibly suggesting that a fraction of the eclipsed gamma-ray emission originates from the apex of the intrabinary shock. Based on these constraints, we estimate the magnetic field strength of the companion to be 30-170 G. This is consistent with the high field strength implied by the X-ray orbital modulation observed.

In a companion paper, we demonstrated that static axion field gradients sourced by dense nuclear matter (\emph{axion hair}) can dominate the near-field electrodynamics of old rotation-powered pulsars, leading to new constraints on light QCD axions and on CP-violating axion-nucleon interactions. This article provides the extended theoretical and numerical framework underlying those results. We begin by providing a detailed description of the sourcing of axion hair from dense nuclear matter, computing self-consistent field profiles for each interaction across the relevant parameter space. We then study the modification to the electrodynamics induced in the polar gap region by these axion gradients; this is done at the analytic level by studying the modification induced by axion field gradients on the effective discharge parameter (computed in the force-free limit of the split monopole magnetic field configuration, and looking at leading deviations from the force-free limit for dipolar field configurations), and numerically by developing dedicated 1D particle-in-cell simulations which capture the leading order dynamical behavior near the star. Our results demonstrate that axion hair serves to either enhance acceleration, or enhance screening, where the relevant effect changes between the northern and southern hemispheres of the star, and between the field lines which support out-flowing and return currents.

Ultra-fast outflows (UFOs) from Active Galactic Nuclei (AGN) are winds which reach velocities up to 0.3c, and show mass outflow rates of 0.01-1 $M_{\odot} yr^{-1}$ and kinetic energies of $10^{42-45}$ erg/s. Their extreme properties allow UFOs to carry substantial kinetic power, sufficient to drive AGN feedback by injecting energy and momentum into the interstellar medium, regulating the host galaxy evolution and the supermassive black hole growth. In this work, we present a detailed analysis of XMM-Newton X-ray spectra of the low-luminosity quasar Mrk205. We find evidence of multi-component ultra-fast outflows, traced by two components detected in the high-resolution data, showing velocities in the range of $\sim$0.07c to $\sim$0.1c. In addition, we marginally detected an ultra-fast outflow feature in the Fe K band. According to our analysis, Mrk 205 shows episodes of wind launches, rather than persistent winds. We further discuss whether the energy carried by the UFOs in Mrk 205 is capable of producing significant feedback on its host galaxy.

Thanks to the almost 400 gravitational wave events detected, we are currently able to grasp the fundamental features of black hole mass, spin, and distance distributions. However, such a fast increase in the precision of the measurements does not necessarily correspond to a better theoretical understanding of gravitational wave sources, especially in current scenarios where the number of free parameters is significantly larger than the number of inferred properties of the black hole population. In this work, we showcase how the landscape of theoretical models can be chipped away by complementary data-analysis strategies, in particular by studying the statistical properties of gravitational wave anisotropic distribution. Specifically, we show how gravitational wave clustering is sensitive to two unique features of each binary formation mechanism: the time-delay distribution and the properties of the binary hosts. First, we consider a model-agnostic scenario and show the impact that different time-delay distributions have on the gravitational wave bias. Then, we consider a realistic scenario where gravitational wave events are sourced either by isolated binary evolution or dynamical processes in globular clusters, and study how the gravitational wave bias is unique sensitive to the specific properties of the environment. In both scenarios, we show how the cross-correlation between galaxy and gravitational wave catalogs is able to distinguish between models with different time delays or with different binary sub-populations originated in specific formation channels.

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

Evidence of a gravitational wave background suggests the existence of a population of sub-parsec supermassive black hole binaries (SMBHBs), with characteristic angular separations on the order of 1-10 $μ$as. Spaceborne extensions of Very Long Baseline Interferometry (VLBI) and the next generation ground arrays introduce the possibility of directly imaging these SMBHBs. In this work, a binary Spectral Energy Distribution (SED) model is used to predict the detectability of SMBHBs with ground and spaceborne VLBI. We consider the Black Hole Explorer (BHEX), a proposal for a spaceborne VLBI mission, as our primary case study. We explore the detectable SMBHB parameter space and identify distinguishable binary signatures that could exist in the visibility domain. We find that for a flux-density-limited sample, ground array observations are more effective at detecting a wider region of the binary parameter space, with $M_\mathrm{tot} \gtrsim 10^9$ solar mass systems detectable out to redshift, $z=0.075$ and beyond. Conversely, inclusion of a spaceborne element such as BHEX, offering finer angular resolution ($\sim6\,μ$as) and sampling of the (u,v) plane not limited by Earth rotation synthesis, will provide significant benefits in constraining binary properties, resulting in improvements in characterisation of the separation and position angle of SMBHBs by a factor of $\sim4$. Near-future ground and/or spaceborne VLBI may achieve the first direct observation of a SMBHB, contributing significantly to multi-messenger studies of such systems with pulsar timing arrays and observations across the electromagnetic spectrum.

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

We present NRHJSur3dq8, a family of numerical-relativity (NR) surrogates for non-eccentric aligned-spin binary-black-hole waveforms inspired by an action-angle parameterization. The model provides both the dynamics: the Hamiltonian values, azimuthal action $J_φ(Φ)$, the orbital clock $τ(Φ)$ as functions of its phasing, and the radiative pieces: the waveforms, their near-analytic derivatives, and modeling errors associated with each. The surrogated variables are conditioned into pieces that vary slowly on the radiation-reaction timescale, represented on shared $hp$-adaptive Chebyshev elements, and regressed across the intrinsic parameters with a Gaussian-process model. A time-parameterized merger-ringdown set of elements is attached where the adiabatic description breaks down. Leave-one-out validation of the central NRHJSur3dq8_AA flavor against NR reaches a $(2,2)$ unweighted mismatch of $1.3^{+32}_{-1.2}\times10^{-6}$ and an all-mode mismatch of $2.3^{+59}_{-2.2}\times10^{-6}$ over the full inspiral-merger-ringdown span. Sky-averaged and SNR-weighted over $32$ orientations per simulation, the released model's full-IMR in-sample mismatch against the NR simulations used for training is $8.9\times10^{-7}$ to $2.8\times10^{-5}$ across total masses $40$-$120\,M_\odot$. Waveforms are hybridized using post-Newtonian waveforms. Beyond accuracy, NRHJSur3dq8 provides two capabilities uncommon in NR surrogates: nearly analytic derivatives of the strain with respect to the physical parameters, and a calibrated Gaussian-process predictive uncertainty of the model. A full non-hybridized inspiral-merger-ringdown waveform evaluates in $12$ ms on a single CPU call. Batched evaluation amortizes this to $0.39$ ms per waveform at a batch of $4096$, measured at a starting frequency of $20$ Hz.

Interest in axions and axion-like particles has resurged, driven by their role in solving the strong CP problem and their appeal as dark matter candidates. The luminosity of the tip of the red giant branch (TRGB) offers a key observable for probing these particle properties. This work aims to improve existing bounds on the axion-electron coupling by adopting a differential observable that is less sensitive to distance, interstellar extinction, the zero-point of bolometric corrections (BCs), and other systematic uncertainties. We use the bolometric magnitude difference between the TRGB and the RGB bump (RGBB) as a distance-independent constraint, applied to three globular clusters of intermediate-to-high metallicity: NGC 104 (47 Tuc), NGC 362, and NGC 5904 (M5). Using photometric catalogs of RGB stars, we determine V- and I-band magnitudes of both features and convert them to bolometric values. After validating that our stellar models reproduce the observed RGBB luminosity, we perform a maximum likelihood analysis with Monte Carlo simulations to propagate uncertainties and derive new bounds on the coupling. A combined analysis of the three clusters yields a maximum likelihood at g13 = gae/10^-13 = 0.8 and a 95% C.L. upper limit of 1.49. Although slightly less stringent than recent bounds from larger multi-cluster samples, this limit is substantially more robust. Under reasonable mass-loss assumptions, g13 > 7.5 is ruled out, as it predicts the disappearance of the HB and AGB phases routinely observed in globular clusters, a limit more than an order of magnitude stronger than current direct experimental bounds such as XENONnT. We demonstrate the effectiveness of this differential, distance-independent method for constraining physics beyond the Standard Model. Applying it to a wider sample of globular clusters would further refine the constraint on gae.

Following important experimental developments, tests of the Diffuse Supernova Neutrino Background (DSNB) are entering a new phase that might lead to a discovery. This motivates the development of simple and broadly applicable tools for comparing theoretical predictions and interpreting experimental results. We show that, over the energy range relevant to current detectors, the DSNB is well described by an exponential spectrum, characterized by the differential flux at 20 MeV, $φ_{20}$, and a spectral energy scale, $E_0$. We fit the spectra of 24 representative models and demonstrate that the exponential description provides an excellent approximation. We then recast the recent Super-Kamiokande results in the $(φ_{20},E_0)$ plane, showing that model-specific best fits and upper bounds can be interpreted naturally within this framework. Our results indicate that, above realistic detection thresholds, and until the precision phase is reached, the observable DSNB is effectively a two-parameter phenomenon. We advocate that the $(φ_{20},E_0)$ parameterization be adopted as a standard, to be used routinely to report future theoretical predictions and experimental results.

We present the first XRISM spectroscopy of the core of Abell 1413 (z = 0.143), combining 104.7 ks of Resolve with archival XMM-Newton and Chandra spectra of a matched sky region. A1413 is hot, massive and relaxed but lacks a strong cool core (t_cool = 6.4 Gyr), a regime in which the low turbulence reported to date has been tested only once, in A3571. The Resolve spectrum requires a single velocity-broadened component with kT = 7.06 (+0.33, -0.43) keV and a line-of-sight velocity dispersion sigma_v = 170 (+25, -24) km/s (M_3D = 0.217). The bulk velocity relative to the brightest cluster galaxy, v_bulk = -24 +/- 26 (z_ICM) +/- 71 (z_BCG) km/s, is limited by the optical redshift rather than by Resolve, and disfavours coherent line-of-sight sloshing. The non-thermal pressure fraction, 2.5 +/- 0.7 per cent, implies a hydrostatic mass bias b < 2.5 per cent, comparable to the strong cool-core cluster A2029 (2.6 +/- 0.3 per cent) and, at matched sampling scale, if anything below it. With A3571 (1.3 +/- 0.3 per cent), to which A1413 rescales within 0.2 sigma, two weak cool cores now show that this quiescence does not require central cooling. The resolved Fe-K He-alpha complex gives Fe = 0.38 +/- 0.03 (stat) +/- 0.04 (sys) solar; Ni/Fe, separated from Fe K-beta for the first time, is solar, as is the pattern as a whole. The Ne, Mg and Si lines that set Ne/Fe, Mg/Fe and Si/Fe lie below 2 keV, where our instruments disagree by 10-20 per cent; we therefore do not regard the SNIa fraction as constrained.

Continuous-wave searches for rotating neutron stars usually cover wide parameter spaces and produce a large number of candidate signals. Stochastic sampling methods therefore play an important role in candidate follow-up and parameter estimation. We present the integration of the Bayesian inference library BILBY into the open-source continuous-wave analysis package PyFstat, enabling access to the broad range of stochastic samplers available through BILBY while preserving the existing PyFstat analysis infrastructure. The implementation is validated through injection studies in Gaussian noise for two representative single-stage follow-up configurations: candidates from a directed search and candidates from an all-sky search. Using the nested sampler DYNESTY, we obtain detection efficiencies consistent with the theoretical sensitivity predictions. In this setup, DYNESTY performs comparably to the existing PTEMCEE-based implementation of PyFstat for the directed case, and also provides an effective single-stage approach for the all-sky candidate follow-up considered here, where we did not identify a PTEMCEE setup with similar recovery performance. Percentile-percentile tests for DYNESTY further show well-calibrated credible intervals, demonstrating reliable parameter estimation. While a full exploration of DYNESTY settings, multi-stage setups, and other samplers is left for future work, these results demonstrate that this flexible and publicly available framework opens useful new possibilities for Bayesian continuous-wave analyses.

We investigate the dynamics of a neutral scalar particle in a charged traversable wormhole background, considering harmonic and modified harmonic oscillator couplings. By solving the Klein--Gordon equation, the radial problem is reduced to a generalized Heun equation, whose polynomial truncation yields the bound-state energy spectrum and quantized oscillator frequencies. Although the particle is electrically neutral, the energy spectrum depends indirectly on the wormhole charge through the geometric deformation of spacetime. We also analyze the dependence of the oscillator frequency on the angular momentum, throat parameter, electric charge, and the additional interaction parameter $V_{0}$. The appropriate limits recover the corresponding results for the uncharged wormhole and the standard harmonic oscillator.

Cluster weak-lensing analyses often compress two-dimensional lensing fields into azimuthally averaged radial profiles and evaluate nonlinear observables from the averaged convergence and shear. However, the ring average of a nonlinear observable generally differs from the mean-field prediction formed from the averaged fields. This difference, divided by the mean-field prediction, defines the fractional residual. For a complete ring, the difference begins at second order in angular fluctuations. For centered elliptical halos, rotational symmetry makes the shape contribution even in signed ellipticity, while the leading miscentering dipole is orthogonal to the shape quadrupole. We test these predictions using projected triaxial NFW halos at six mass-redshift grid points spanning $3\leq M_{200\mathrm c}/(10^{14}\,h^{-1}M_\odot)\leq20$ and $0.2\leq z_l\leq0.5$, for $z_s=1$. For the reference offset model, centering offsets follow a Rayleigh distribution with scale $0.05r_{200\mathrm c}$. For $0.366\leq R/r_{200\mathrm c}\leq1$, where at least $95\%$ of each population satisfies a conservative subcriticality criterion, we calculate the median of the fractional residuals across the retained halos for each population and radius. The largest population-median magnitudes are $0.90\%$ for reduced shear, $0.0093\%$ for inverse magnification, $0.64\%$ for magnification, and $0.18\%$ across the two magnification-bias cases $μ^{α-1}$ ($α=0.3, 1.4$). In the paired calculation at $z_s=2$, the maximum magnitude increases for every observable, consistent with weak-lensing power counting. Across both source planes and the common radial domain, every population median remains below $2\%$ in magnitude. This accuracy follows from first-order cancellation, rotational symmetry, and harmonic orthogonality, with further weak-lensing suppression of the remaining nonlinear terms.

During their evolution, stars follow a distinct path in luminosity--temperature space. Low-mass stars that have exhausted hydrogen in their core follow the so-called red-giant branch which is predominantly in the direction of increasing luminosity and decreasing surface temperature. The luminosity bump is a temporary decrease in luminosity against this otherwise increasing trend. The bump is present in both observations and stellar evolution models. However, with canonical physics included in the models, the bump in computed evolutionary tracks appears at higher luminosities than the observed bump. To understand why the bump appears later in the models than in observations, the physics of the bump needs to be unravelled. The end point of the bump is well-understood; however, the onset of the bump is still an enigma. Here, we report on the physical origin of the onset of the luminosity bump. We show that the difference in specific entropy at the mean molecular weight discontinuity decreases due to the discontinuity moving in. This decrease is attributed to the decrease in the ratio of the temperature to the pressure at smaller radii. Upon reaching a critical value the specific entropy difference at the mean molecular weight discontinuity is reduced sufficiently to reduce the specific entropy in the convective envelope. The latter is a key signature of the bump. Hence the evolution of the specific entropy at the mean molecular weight discontinuity provides a viable description for the onset of the bump.

We present stacked thermal Sunyaev-Zel'dovich (tSZ) effect measurements for three samples of galaxy groups and clusters: those traced by the Dark Energy Spectroscopic Intstrument Data Release 1 (DESI DR1) luminous red galaxies (LRG) and the DESI DR1 Bright Galaxy Sample (BGS), and an eROMaPPer optically-selected sample from the DESI Legacy Imaging Survey. We use the latest Atacama Cosmology Telescope DR6 (ACT)+Planck component-separated internal linear combination (ILC) Compton-$y$ maps and ACT+Planck coadded 90, 150, and 220 GHz temperature maps to extract the tSZ signal within a $\sim2'$ disk aperture for sources binned by luminosity, richness, or mass. We measure the average tSZ signal with high statistical significance, with signal-to-noise ratios surpassing 38 for LRG, 27 for BGS, and 39 for the eROMaPPer sample using the 90 GHz ACT DR6+Planck map. We conduct a detailed study of systematics and foregrounds such as dust and cosmic infrared background (CIB) contamination, which remain a core challenge for tSZ analysis. For the LRG and BGS samples, we find that dust and radio source emission dominate the tSZ signal at scales near and below the disk aperture radius. Large-scale ($R>4'$) contamination from the CIB is less significant. We mitigate these contaminants to isolate the tSZ signal and use a combination of simulated and real measurements to develop Compton-$y-$optical depth ($\bar y-\bar τ$) scaling relations to infer optical depths, which are found to be in agreement with values measured using the pairwise kinematic SZ effect for the same tracer samples. The $\bar y-\bar τ$ scaling relation for the eROMaPPer sample is the first such relationship to be derived directly from SZ measurements.

The recent observations of the supermassive black holes M87* and Sagittarius A* by the Event Horizon Telescope collaboration (EHT) have shown the need for robust, efficient, and accessible numerical tools to model light propagation in strongly curved spacetimes. Comparing theoretical accretion models with observational data relies heavily on ray-tracing algorithms that calculate the trajectories of photons traveling from a source plasma to a distant observer. In this paper, we present TARTARUS (Tracer for Astrophysical Ray Trajectories Around Relativistic Ultra-compact Sources), a new, highly modular, and extensible Python-based framework for computing null geodesics around compact objects. By leveraging Just-In-Time (JIT) compilation via Numba, TARTARUS bridges the gap between Python's high-level accessibility and the execution speed of compiled languages like C++ or FORTRAN. The code supports multiple native solvers (including adaptive embedded Runge-Kutta pairs, a Bulirsch-Stoer extrapolator, and a Verlet scheme), event-detection mechanics for handling structural intersections (e.g. event horizons and accretion disks), and support for both analytical and numerically generated background metrics. We demonstrate the code's physical accuracy and computational efficiency through some tests, including the evaluation of the Hamiltonian constraint conservation and ray-tracing of the shadow and the Novikov-Thorne thin accretion disks around a Kerr black hole.

The LUX-ZEPLIN (LZ) experiment has recently reported a nuclear-recoil candidate at $E_R=248\pm23_{\rm stat}\pm23_{\rm sys}\,{\rm keV}$ in an exposure of $2.84\,{\rm tonne\,yr}$. At this unusually high recoil energy, the spin-independent (SI) xenon response is strongly suppressed, motivating dark matter interactions with a harder recoil spectrum. We show that elastic spin-dependent (SD) scattering through the nonrelativistic operator $\mathcal O_4=\mathbf S_{\rm DM}\cdot\mathbf S_N$ provides a simple realization of such a spectrum. We consider singlet-doublet Majorana dark matter, for which a diagonal $Z$ coupling generates the required SD interaction. The accompanying Higgs-mediated SI interaction would generically produce too many low-energy recoils, but can be suppressed along a Higgs blind spot while retaining a nonzero $Z$ coupling. We find a region of the Higgs blind spot that simultaneously predicts an $\mathcal O(1)$ event rate in the LZ high-recoil window and reproduces the observed thermal relic abundance. The same SD interaction leads to solar capture, allowing this region to be tested independently with solar-neutrino searches. Current IceCube limits already probe part of the LZ-compatible thermal parameter space, while a viable region remains. Additional xenon exposure and improved solar-neutrino searches can therefore provide complementary tests of this interpretation.

The form and relationship between the mass distributions of Wolf-Rayet stars and black holes are analyzed using a direct method for estimating the distribution function from a sample of mass values with individual distributions obtained from current observations. It is shown that, taking into account the mass loss of WR stars by the end of this evolutionary stage, the mass distribution of black holes is close to that of the CO-cores of WR stars before gravitational collapse, $M_\mathrm{BH}\simeq M_\mathrm{CO}$.

In-situ measurements in space plasmas indicate that the velocity distributions of charged particles are not in thermal equilibrium, deviating from a standard Maxwellian mainly due to anisotropies and suprathermal populations which enhance high-energy tails. Although the Standard \k{appa}-Distribution (SKD) is well established in modeling these non-equilibrium distributions, its application is often viewed controversially due to certain unphysical implications, in particular divergent velocity moments. To address these issues, the Regularized \k{appa}-Distribution (RKD) was introduced. Such advanced, in general anisotropic RKDs are invoked here for the first time to investigate oblique firehose instabilities, including those induced by the temperature anisotropy of electrons and protons, the dominant species in space plasmas. In weakly collisional plasmas, both of these instabilities are expected to play significant roles in the self-regulation of the macroscopic properties of space plasmas (e.g. the expanding solar wind) reported by observations. Is not yet possible to resort to a general dispersion tensor related to RKD plasmas (whose derivation is still a challenge), instead the Arbitrary Linear Plasma Solver (ALPS) is exploited here. The unstable solutions obtained for the already established Maxwellian & SKDs successfully validate the capability of ALPS. For RKDs, the instabilities confirm the stimulating effect of suprathermal populations, with lower \k{appa} values generally enhancing the firehose growth rates. For the oblique electron firehose instability, RKDs substantially modify the competition between periodic & aperiodic branches only at intermediate angles & can sustain significantly increased growth rates of especially aperiodic modes at highly oblique propagation beyond both the Maxwellian unstable regime & the one that is accessible with SKDs.

Supermassive black holes (SMBHs) in "little red dots" (LRDs) discovered the James Webb Space Telescope (JWST) may result from runaway mergers of primordial black holes (PBHs) in clusters---through long--short mode coupling on small scales in the early Universe. In this framework, we derive the SMBH mass function, together with the compactness and overmassive features of LRDs. We also estimate that the dense gas residing in PBH clusters is consistent with LRD observations. In addition, SMBHs formed from PBH clusters can help accelerate galaxy formation at high redshifts, thus alleviating tension with $Λ$CDM cosmology.

Variability in O-type supergiants is a well known phenomena although their origin is not fully understood. In very massive, luminous stars, strange mode pulsations have been suggested to play a role in the variability. Here, we present a preliminary theoretical study of O-type supergiant models corresponding to the observed parameters of the star HD 152249, which shows line profile variability. Non-adiabatic linear stability analysis with respect to radial perturbations is carried out for the considered models of this star. Instabilities with growth rates of the order of dynamical timescale are found, with characteristics suggestive of strange modes. Non-linear simulations for selected models indicate that the instabilities could lead to significant envelope inflation, finite amplitude pulsation and may contribute to enhanced mass loss.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

In the second half of 2017, as Solar Cycle 24 approached its minimum, a pronounced and rigidity-dependent depression in the galactic cosmic ray (GCR) proton flux was observed using AMS-02. This event was driven by a sequence of strong coronal mass ejections (CMEs). In this work, we model the modulation of GCR protons during this period using a steady-state Parker transport framework with parameters constrained by AMS-02 observations. Treating the event as a purely long-term modulation interval yields excellent agreement with the observed proton intensities, but implies an inferred rigidity dependence of particle mean free paths that is inconsistent with basic turbulence theory, corresponding to an effective inversion of the scattering regime. We argue that this behaviour arises from the ability of highly flexible transport parameterizations to reproduce integrated modulation signatures through correlated adjustments of diffusion coefficients during transient events. To address this inconsistency, we relax the commonly adopted assumption of a time-independent ratio of perpendicular to parallel diffusion coefficients and introduce a physically motivated reduction of this ratio during the CME-dominated interval. This modification preserves the quality of the spectral fits whilst restoring the expected rigidity ordering of mean free paths. Our results indicate that allowing the perpendicular-to-parallel diffusion ratio to vary in time provides an effective means of constraining transport solutions during the 2017 depression and may help to resolve similar inconsistencies reported in other GCR modulation studies.

We present deep X-ray observations of the nearby X-ray-luminous galaxy cluster Abell 1795 obtained with the Einstein Probe Follow-up X-ray Telescope (EP-FXT), with a total exposure time of ~480 ks. Exploiting the large field of view and low particle background of EP-FXT, we directly measured the radial temperature profile of A1795 out to R<SUB>200</SUB> with full azimuthal coverage, which was shown to increase in line with the radius within 6' and then gradually decline toward larger radii. The surface-brightness residual map and 2D thermodynamic maps reveal a clockwise spiral structure extending from the cluster core towards the south-east, which traces low-temperature and low-entropy gas, which is also consistent with sloshing-induced cold fronts. In the north-west, the surface-brightness-enhanced region exhibits an arc-like high-temperature feature and both the temperature-derived and density-derived Mach numbers support the assumption that it is a weak shock. These substructures can be explained by a binary merger scenario: the perturbing subcluster induces sloshing during its first passage past the primary core and its subsequent return passage through the ICM might drive the shock towards the north-west. Our results indicate that relaxed galaxy clusters such as A1795 might still retain signatures of dynamical activity.

Long Period Variables (LPVs) are Asymptotic Giant Branch stars with periods >100 days. Around one third to half of LPVs display a Long Secondary Period (LSP) that is 5-10 times the length of the main periods. The exact cause of LSPs are unknown, with one potential explanation lying in cross-harmonics caused by the interaction of two simultaneously occurring periods. In this paper, we use the Fourier analysis software VStar to examine the visual range photometric light curves obtained through the American Association of Variable Star Observers of a sample of 5 stars. We find that one star, V CVn, displays a secondary period over 20 times the length of its main period, inconsistent with traditional LSPs. To test if cross-harmonics can result in an LSP or VCVn's extra-long secondary period, we then construct evolutionary tracks and pulsation models of LPVs using Modules for Evolution in Stellar Astrophysics (MESA) and GYRE. After examining the interactions of the fundamental pulsation mode and the first 4 overtones with each other, we do not find evidence that LSPs or V CVn's extra-long secondary period could be caused by cross-harmonics.

We introduce a merger tree builder, CTree. Its key feature is the identification of progenitors using the evolutionary information encoded in object branches, rather than relying solely on the instantaneous particle membership in a single snapshot. In this sense, CTree can be described as a branch-wise merger tree builder. Specifically, CTree extracts particles associated with objects over a defined timescale and uses this information to compute the merit functions employed in the progenitor search across multiple snapshots. During this process, CTree requires only the identifiers of particles and objects. This ensures high compatibility with different halo and galaxy catalogues and also with pre-existing merger trees. To identify `core' particles that remain bound to the host halo over extended physical timescales, CTree employs a simple and efficient algorithm that does not rely on explicit radial or binding-energy information. The timescale adopted for core particle extraction and the number of snapshots used in the progenitor search can be user-configured. Basic validation tests using suites of dark matter-only simulations demonstrate that CTree produces stable and physically consistent halo merger trees, particularly in high-density regions or in simulations with high output cadence. The code is written in C++ and employs a hybrid message passing interface and OpenMP parallelization, making it well-suited for high-resolution cosmological simulations. The code is publicly available and actively maintained to improve performance, reliability, and compatibility.

We analyze high-cadence spectropolarimetric observations of solar plage obtained with the Visible Spectro-Polarimeter (ViSP) on the Daniel K. Inouye Solar Telescope (DKIST) to investigate low-chromospheric waves sampled by the Na I D1 5896 Angstrom and Ca II 8542 Angstrom lines. Observations were taken at an oblique viewing angle, providing sensitivity to transverse motions. Magnetic oscillations within plage elements exhibit significant line-of-sight (LOS) RMS amplitudes of approximately 9 G on average in the Na I D1 line, whereas the corresponding Ca II oscillations are only marginally above the noise level. Adaptive profile tracking provided line-core intensities and Doppler velocities, the weak-field approximation provided LOS magnetic fields, and non-LTE inversions provided mass densities and formation heights. Cross-spectral analysis shows velocity-intensity phase distributions coupled at -90 degrees, consistent with compressive oscillations. The velocity phases in Na I-Ca II indicate upward propagation and the phase speeds are predominantly sub-Alfvenic, consistent with slow magnetoacoustic waves. Although the slit geometry prevents a definitive identification of the tube eigenmode, the differing magnetic phase relationships suggest a height-dependent sausage-like behavior in the coupling between the magnetic and compressive perturbations. Additionally, our inferred wave fluxes are insufficient to heat active-region plage.

The evolution of galaxy cluster abundance is a powerful probe of cosmology, but its efficacy depends on our ability to understand and control for systematics in the selection and mass estimation of clusters. We study redMaPPer-selected clusters in the Cardinal simulation, which simulates the DES Y3 observational conditions. Specifically, we investigate the line-of-sight velocity distributions of cluster member galaxies in 12,272 simulated clusters. We find that a significant fraction of clusters, 35%, have velocity dispersions that are high compared to their richness with an offset of $\sim 400$ km s$^{-1}$ between the outliers and the main population. A similar population of velocity dispersion outliers was found in the DES Y3 data. Nearly all of the outliers (98%) in a low-richness subsample have non-Gaussian velocity distributions, and they are more likely than the main population to have significant velocity substructure. These clusters also tend to have higher redshifts than clusters with typical velocity dispersions, and they show a larger scatter between halo mass and richness. Overall, our results point to a significant contamination of redMaPPer-selected samples from line of sight structure, particularly at high redshifts and low richness. Comparison to the Buzzard simulations shows similar results.

We investigate extremal black holes with scalar hair in the generic Einstein-Euler-Heisenberg (EEH)-scalar theory with two scalar couplings to the Maxwell term. One is an exponential coupling with coupling constant $α$ and the other is its polynomial coupling. We first construct scalarized black holes on the cold (C)-horizon of EEH black holes described by mass $M$ and magnetic charge $P$ both at the linear level, through the existence curves $α_n(q)$ with $q=P/M$, and as fully backreaction solutions. As extremality ($q=q_e$) is approached, all existence curves accumulate, following $α_n-α_c\sim 1/\ln^2(q_e-q)$, at the critical branch $α_c$ fixed by the Breitenlohner-Freedman bound for the near-horizon (AdS$_2\times S^2$) throat. We obtain scalarized extremal black hole (SEBH) with constant secondary hair and it is recovered exactly from the entropy function approach working on the near-horizon throat. Its entropy is an attractor invariant, being independent of the asymptotic modulus, so the scalar hair remains secondary. Exploiting this constant scalar, we seek further SEBHs with scalar hair keeping its charge $Q_s$ for the zero asymptotic scalar ($φ_\infty=0$). In the $(q,α)$ plane, we observe that these extremal solutions occupy $q\ge q_e$. Hence, the existence curves $α_n(q)$ existing for $q\le q_e$ and the extremal branches bound the scalarized domain from opposite sides and they meet only at $q=q_e$.

High-redshift quasars ($z>5$) are crucial for studying early black hole growth and reionisation. SDSS Data Release 16 catalogues 655 such sources with multiple redshift estimators. We systematically examined this subsample and found 480 sources (73.2\%) where pipeline redshift (Z\_DR16Q) disagrees with the derived estimate (Z\_FIT) significantly. Visual inspection confirms that 125 of these (19.1\% of the total) are actually low-redshift active galactic nuclei. In 114 cases, broad Mg\,{\sc ii} with C\,{\sc iii} and C\,{\sc iv} was misidentified as Ly$α$; in 11 cases, H$α$ was misidentified as Ly$α$. This error is immediately revealed by the absence of Gunn--Peterson trough blueward of Ly$α$. Crucially, Z\_FIT recovers correct redshift in the affected sources where emission lines are prominent, indicating the algorithm is sound; the issue lies in community preference for Z\_DR16Q or Z\_SPEC. We urge caution against uncritical use of pipeline redshifts for $z>5$ DR16Q sources and recommend Z\_FIT as the preferred estimator.

Active region atmospheres host many oscillatory phenomena. The chromosphere is delimited by steep temperature gradients at the photosphere and transition region, where magnetoacoustic waves are trapped and can form standing oscillations within this resonant cavity. We aim to detect the signature of the resonant nodes of standing waves, which are expected to produce sudden jumps in the oscillatory phase and power dips. Spectroscopic temporal series of H$α$ in a pore were acquired with the Swedish Solar Telescope. The velocity and temperature fluctuations at multiple atmospheric heights were inferred from the analysis of the intensity at many spectral positions along the line wings. Wavelet analysis was employed to characterize the phase differences and the power at different heights. The phase shift between velocity and temperature shows a $\pm90^{\circ}$ value, which is consistent with standing oscillations. Robust evidence of the presence of a nodal layer in the temperature at around the height probed by the intensity at H$α\pm0.30$ Å is found, such as the detection of 180$^{\circ}$ jumps in the phase of the temperature oscillations and remarkable power dips at the same atmospheric layer. The exact height of the resonant nodes depends on the spatial location and time. We generally find a mixture of standing and propagating waves. This is consistent with a leaky resonator where waves are partially reflected at the transition region, while some of them are able to propagate into the corona. For the first time, we report the detection and characterization of resonant nodes in the solar chromosphere. This result provides strong observational support for the chromospheric resonant cavity model and paves the way for the development of new seismological techniques to investigate the structure of active region chromospheres.

We present a simple model of a higher dimensional spacetime in which the 4d submanifold is a FLRW metric, while the extra dimensions are compactified to a hypersphere. We calculate the Friedmann equations of this model finding that the extra dimensions interact in a non-trivial way with the 4d submanifold evolution. We prove that what is observed as Dark Energy could be explained only as an effect due to the Gaussian curvature and the expansion rate of the extra dimensions. We calculate the luminosity distance in this model and derive the acceleration parameter q. We find that Ωm and q are in some sense dressed by terms coming from the extra dimensions, thus a Universe filled only with matter (both luminous and dark) and with extra dimensions can explain the observational data. We find that the extra dimensions do not affect the CMB spectrum, but this is not a solution to the Hubble tension. We fit our model to SNe data and find evidence for the presence of 3-4 extra dimensions which are presently contracting.

Rotating black holes in dynamical Chern-Simons gravity are known nonperturbatively mainly through numerical solutions, which limits their direct use in strong-field applications. We construct a Kerr-referenced analytical parametrization of the near-general-relativistic numerical branch using a compact radial coordinate, low-order angular multipoles, and continued fractions. The resulting two-parameter family accurately reproduces the numerical metric throughout the validated weak-coupling and moderate-spin domain. Independent off-grid solutions confirm that the parametrization retains its accuracy away from the calibration set. We further assess the model through equatorial light rings, critical impact parameters, and full two-dimensional null-geodesic ray tracing without assuming separability. The resulting shadow observables closely match those obtained from the numerical geometries. This parametrization therefore provides a practical analytical metric for geodesic and strong-field calculations without relying on a slow-rotation expansion.

We investigate the Casimir effect of a complex scalar field on a two-dimensional sphere threaded by a fixed Wu--Yang magnetic monopole at the center. In this background the charged scalar field is a section of a nontrivial complex line bundle over $S^2$, reflecting the monopole's nontrivial topology. We solve the Klein--Gordon equation analytically and compute the Casimir energy using a generalized Abel--Plana formula for multivalued functions. We find that the monopole can reverse the sign of the Casimir pressure, turning an attractive force into a repulsive one when the monopole strength is sufficiently large. This behavior persists across a broad range of curvature couplings and suggests that Wu--Yang monopoles may provide a source of repulsive vacuum stress, with possible implications for gravitational and cosmological settings.

We present an iteration algorithm for vacuum and Einstein scalar-field equations in double-null gauge, which transform the non-linear PDE into systems of ODE. The numerical realization combines characteristic constraint solves, LGL spectral elements, pole-free spherical operators, Galerkin projection, and independent first-order residual and consistency checks.

The matter power spectrum is one of the primary statistical descriptors of the large-scale distribution of matter in the Universe and provides a powerful probe of cosmic structure formation. Measurements of cosmic microwave background (CMB) lensing offer an integrated view of the matter distribution over a wide range of redshifts, enabling the reconstruction of the underlying matter power spectrum. In this work, we reconstruct the reference linear matter power spectrum $ P_\text{lin}(k,0)$ from the baseline joint CMB lensing measurements of Planck PR4, ACT DR6, and SPT-3G using a covariance-weighted modified Richardson-Lucy(MRL) deconvolution algorithm. The reconstructed spectrum is found to be consistent with the fiducial linear prediction on large scales, while exhibiting a systematic enhancement for $k \gtrsim 0.1\,{\rm Mpc}^{-1}$, where nonlinear gravitational evolution becomes important. To investigate this behavior, we introduce a scale-dependent correction factor, $A(k)$, defined through $P(k)=A(k)\,P_{\rm nl}(k),$ where $P_{\rm nl}(k)$ is the fiducial nonlinear matter power spectrum obtained from 2LPT simulations. The reconstructed correction factor remains consistent with unity within $2σ$ confidence over the reconstructed range, indicating that the observed enhancement is well explained by the standard nonlinear evolution of the matter power spectrum. In addition, the reconstruction shows agreement with the fiducial BAO template around the BAO feature at $k\sim(0.04-0.06)\ {\rm Mpc}^{-1}$, indicating that some BAO-scale information survives the lensing projection.

Within finite-temperature effective potential framework, we study the effects of PBHs on vacuum decay, the associated gravitational wave spectrum, and the baryon asymmetry. For the PBHs considered, the Hawking temperatures are below the ambient temperature, and their dominant effect is taken to be the gravitational distortion of the Higgs field configuration. Relative to flat spacetime, the PBH background reduces the vacuum decay exponent, enhances the gravitational wave peak amplitude, and shifts the peak frequency downward. For PBH masses of roughly $10^{11}~{\rm g}-10^{13}~{\rm g}$, the sphaleron energy is reduced, while the small PBH abundance renders the associated washout negligible. Additionally, over roughly the same mass range, PBH-emitted BSM particle decays generate a baryon asymmetry too small to spoil the agreement between the baryon abundances inferred from BBN and CMB.

Asteroseismic estimates of the magnetic field strength in the radiative interior of red giant stars depend strongly on the internal stellar structure derived from models. Since red giant branch merger products have been shown to be able to possess a different core structure than single stars of the same mass, we investigate how a mass-gain history influences our estimates of an internal magnetic field strength. We construct stellar models with and without a mass-gain event after the onset of the RGB evolutionary phase with masses of $1.1\,M_\odot \le M \le 2\,M_\odot$. First, by assuming a weak magnetic field, we investigate the influence of a mass-gain event on the global sensitivity of the oscillation frequencies to the magnetic field. We find that mass-gain models can be several times more sensitive to the field than single-star models of identical total mass at masses above $1.6\,M_\odot$. Therefore, considering a mass-gain evolutionary history for merger candidates allows a significant correction to the magnetic field strength. In the presence of strong magnetic fields, we also show that the critical field strength needed to suppress mixed dipole modes is significantly lower if a mass-gain event has occurred (for stars with masses $M\gtrsim1.6\,M_\odot$). The massive end of the suppressed stars' distribution is therefore strongly favored by a merger origin. We conclude that properly constraining the stellar evolutionary history is essential when aiming to constrain internal field strength estimates from asteroseismic observations.

SDSS J0946+1006 is an attractive target for measuring cosmological parameters through lens modelling. It is the best-studied galaxy-scale strong gravitational lens with multiple sources whose redshift separations are well-suited to constraining the dark energy equation of state. However, multi-plane lens models with free cosmological parameters risk a multi-plane mass-sheet degeneracy, although this can be lifted by a non-lensing deflector density profile tracer. We simultaneously reconstruct near-infrared and near-ultraviolet HST imaging whilst including a velocity dispersion measurement from VLT-MUSE to constrain the foreground deflector. Imaging is reconstructed faithfully regardless of whether the inferred kinematics are realistic, though we find the kinematic constraint essential for shifting the preferred cosmology into a region not in significant tension with other dark energy probes. Deflector density profile perturbations, via substructure and multipolar halo shape deformations, have only a modest effect on inferred cosmology. Combining our fiducial model with Planck CMB data yields $w=-1.01^{+0.08}_{-0.13}$; or, combined with Pantheon SNe Ia, $w=-0.99^{+0.13}_{-0.15}$. We further show that this system proves a remarkably valuable complementary probe in the $w_{0}$-$w_{a}$ plane of an evolving dark energy model, and that the DESI BAO tension with $Λ$CDM seen when combined with other datasets is not reproduced when combined with this lens, yielding $(w_{0}, w_{a})=(-0.87^{+0.10}_{-0.11}, -0.22^{+0.28}_{-0.25})$. The system's third source, visible with MUSE, only weakly constrains $w$CDM but may strengthen $w_{0}w_{a}$CDM constraints, though further mass-model complexity along its line of sight is required. Overall, kinematics-informed multi-plane lens modelling is a robust route to competitive dark energy constraints, even with a single system.

We revisit gravitational decoherence of matter superpositions using the leading soft sector of the graviton $S$-matrix. Starting from Weinberg's soft theorem and the overlap of coherent radiation states, we derive, for branch configurations with the same total four-momentum, a gauge-invariant, nonnegative leading-soft decoherence exponent and evaluate the polarization sum and angular integrals exactly at this order. In the vacuum, coherence between branches with distinguishable soft radiation vanishes as the infrared cutoff is removed. The leading nonrelativistic term is purely quadrupolar as a consequence of energy--momentum conservation. We extend the analysis to single- and two-mode squeezed graviton states, including the two-mode case motivated by inflation. Squeezing can enhance or suppress decoherence depending on the relative phase, while phase averaging gives an enhancement governed by the graviton occupation. For a phase-averaged relic background with an approximately flat energy-density spectrum over the relevant soft frequency band, the sharp-cutoff identification of the infrared scale with the inverse observation time gives a contribution that grows with the gravitational-wave energy density and with the fourth power of the observation time. These results connect soft-graviton scattering, infrared quantum information, and decoherence in nonvacuum gravitational backgrounds.

We present $\href{https://github.com/akmeena766/LensFactory.jl}{\texttt{LensFactory.jl}}$, an efficient, open-source, general-purpose strong lens modeling package written in Julia. It supports single- and multi-plane lensing with a wide range of analytic lens profiles, and its modular design, built on Julia's multiple dispatch, makes it straightforward to add new profiles or inference methods, while native multi-threading parallelizes the computationally expensive steps, enabling full cluster-scale lens reconstructions within hours on a desktop-class machine. We validate the mass-reconstruction pipeline of $\href{https://github.com/akmeena766/LensFactory.jl}{\texttt{LensFactory.jl}}$ on two simulated galaxy clusters (Ares and Hera), recovering the convergence and magnification maps well where the lensing constraints lie, and on the JWST-observed cluster SMACS J0723.3-7327, obtaining a strong lens model (image-plane RMS of 0.46'') consistent with previously published results. In addition, we introduce singularity maps as a complementary diagnostic for comparing lens models, probing differences in the higher-order derivatives of the lensing potential, which remain largely unconstrained by standard observables such as image positions, time delays, and flux ratios. $\href{https://github.com/akmeena766/LensFactory.jl}{\texttt{LensFactory.jl}}$ is accompanied by a set of worked examples covering all core functionality.

This work analyzes the detection prospects for weakly interacting massive particles (WIMPs) in dark matter (DM) density spike around the supermassive black hole (SMBH) in Andromeda galaxy M31 in microwave band. WIMP annihilation produces relativistic electrons and positrons, which generate the synchrotron emission in strong magnetic field of the accretion flow medium at microwave frequencies. Properties of this emission, as well as their dependence on DM and medium parameters, were modeled. Atacama Large Millimeter Array (ALMA) sensitivity to this emission and underlying WIMP parameters was estimated. It was obtained that M31* in microwave band represents the unique target, where DM spike signal can be detected in case of realistic spike and halo density profiles. In such case ALMA may reach the thermal s-wave annihilating WIMPs with masses up to around 1 TeV with just 10--100 hours exposure. Optimistic parameter configurations may allow probing all possible WIMP masses up to 100 TeV. Related systematic uncertainties are discussed.

Transiting dust and debris around white dwarfs provide a direct view of planetary disruption. We report six previously non-variable white dwarfs that underwent single, deep obscuration events, including ZTF~J2327+0019, which fully disappeared for more than a year. These discoveries raise the number of confirmed debris-transit white dwarfs to 20. In addition, the systematic search of negative alerts from the Zwicky Transient Facility (ZTF), followed by photometric and spectroscopic observations, also identified three hot ($T\approx 55000$--$92000$ K) irregular variables as candidate debris hosts. The six confirmed events last weeks to years, reach transit depths of 25--100%, and show no recurrence over the 8-yr ZTF baseline. They are therefore much deeper and longer lived than the recurring, minute-scale transits typical of previously known systems. The observed reddening is consistent with dust extinction, and high-speed photometry reveals minute-scale structure in several systems, indicating rapidly moving material close to the white dwarf. ZTF~J2327+0019 also shows a simultaneous WISE infrared brightening, the first observed coincidence of infrared brightening and optical obscuration in a white dwarf. Two of the six hosts, one DAZ and one DBZ, are strongly metal enriched, with accretion rates among the highest known for their spectral types and minimum accreted masses of $3.5\times 10^{19}$ and $2\times 10^{24}$ g. We interpret the one-off transit events as collisions between planetesimal remnants and pre-existing circumstellar debris that trigger collisional cascades: short-timescale variability traces dust-producing fragments crossing the star, while slower changes trace the build-up and clearing of an optically thick dusty structure. In contrast, ZTF~J0436+5728 shows a coherent 5.910-h period, consistent with a rocky body shedding material near the tidal-disruption radius.

Can a universe undergo a regular Euclidean-Lorentzian signature transition when the gravitational coupling itself is dynamical? We address this question in scalar-tensor gravity with nonminimal coupling $F(φ)R$. Although the spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) sector can be mapped to the Einstein frame for $F>0$ and a nondegenerate scalar redefinition, not all transition properties are frame independent. For a finite, positive, and sufficiently regular conformal factor, the existence and transverse character of the type change are preserved, whereas the extrinsic geometry is not. We therefore use the Einstein frame as a solution-generating representation and impose total geodesy in the physical Jordan frame. We construct two exactly integrable classes of solutions. In the oscillator-ghost-oscillator branch, the scalar field is stationary at the transition and Jordan-frame total geodesy follows under suitable regularity assumptions on the conformal factor. In the critical exponential-potential branch, the scalar is generically nonstationary and total geodesy instead requires \begin{equation*} H_E\big|_Σ= \frac12 \left( \frac{d\ln F}{dψ} \right)_Σ\dotψ_Σ. \end{equation*} We also distinguish this condition from the stronger smoothness requirements of a Kossowski-Kriele-type transverse metric. Finally, while the canonical Einstein-frame scalar does not allow effective phantom evolution, the dynamical nonminimal coupling can generate a locally superaccelerating Jordan-frame regime near the transition when $\left(\frac{d\ln F}{dψ}\right)_Σ>0$ for the chosen oscillator branch. Thus signature change persists beyond Einstein gravity, with regularity and effective cosmological behavior remaining intrinsically frame sensitive.

Currently, the only second-order gravitational self-force calculations have been based on directly solving the perturbative Einstein equations in the Lorenz gauge. That method relied on the complete separability of the Einstein equations in a Schwarzschild background. In this paper, we present a new scheme based on the second-order Teukolsky equation. Crucially, this method promises to extend (reasonably straightforwardly) to the more realistic case of a Kerr background. Here we implement the scheme in the simplest setting of quasicircular orbits around a Schwarzschild black hole. In addition to working with the Teukolsky equation, our scheme incorporates several other advances over previous second-order self-force calculations: compactified hyperboloidal slicing, transformation to a Bondi-Sachs gauge, and a combination of spectral and variation-of-parameters methods. We also use these tools to re-examine the infrared divergences that arise in second-order Lorenz-gauge calculations, showing they are less pronounced in the Teukolsky case and completely eliminated in the Bondi-Sachs gauge. We conclude by calculating the asymptotic energy fluxes and benchmarking them against previous Lorenz-gauge calculations.

We investigate three flavor neutrino phenomena in an effective holonomy corrected Schwarzschild geometry. We derive the weak-deflection angle through second post-Minkowskian order, obtain the semiclassical phases along radial and nonradial trajectories, and formulate the two-image flavor probability by including magnifications, Fermat phases, and wave packet overlap. The cross path phase acquires a logarithmic holonomy contribution and can retain information about the absolute neutrino mass scale. We also characterize flavor mode correlations and trajectory-flavor entanglement through the balance among path predictability, interference visibility, and I-concurrence. Neutrino-antineutrino annihilation outside an effective neutrinosphere is examined, with holonomy entering the integrated power through the radial proper volume measure. For a solar mass lens and $10\,\mathrm{MeV}$ neutrinos, the numerical results show displaced oscillation fringes and a geometry dependent redistribution among the electron, muon, and tau channels. The electron-neutrino trajectory-flavor concurrence reaches approximately $0.23$, with an entropy near $0.10$ bit, while remaining weakly sensitive to the mass ordering in the selected configuration. By contrast, the annihilation power increases monotonically, reaching an increase of approximately $39.3\%$ relative to Schwarzschild spacetime (for a particular configuration, i.e., $λ=2$ and $M/R_ν=1/3$). Furthermore, at fixed source parameters and $M/R_ν=1/3$, an assumed maximum excess of $10\%$ in the integrated annihilation power relative to Schwarzschild yields the conditional bounds $a/(2M)\leq0.2879$ and $|λ|\leq0.6358$.

Gravitational-wave searches do more than identify individual sources - they provide a way to infer the underlying astrophysical populations that produce them. Currently, Pulsar Timing Arrays (PTAs) constrain the supermassive black-hole binary (SMBHB) merger-rate density through both the stochastic gravitational-wave background (SGWB) and searches for individual SMBHB signals. The latter implies converting source upper limits into rate upper limits using detection efficiencies estimated empirically. We instead develop a hierarchical Bayesian framework that places the SMBHB merger-rate density model directly inside the PTA likelihood. The unknown catalog of individual SMBHB merger signals is modeled as a Poisson point process on source-parameter space, so that the number of merger signals is inferred from the data rather than imposed by a detection threshold. We describe two computational routes: a catalog-marginal likelihood based on analyses with fixed numbers of candidate sources, and an explicit transdimensional sampling approach that jointly samples population hyperparameters, latent merger signals, and noise. We validate the method with a toy merger-only population model and then apply it to an astrophysical SMBHB merger-rate density model that jointly predicts the SGWB amplitude and the expected merger-catalog size. In simulations, an individual merger signal adds complementary information to the SGWB and can tighten constraints on the merger-rate population.

Searches for parity violation (PV) in the Universe's large-scale structure have recently been carried out using parity-odd modes of the galaxy 4-Point Correlation Function (4PCF). Here, we introduce a new estimator for a 3-Point Correlation Function (3PCF) that is sensitive to PV. The construction exploits the handedness of triangles formed on 2D spherical shells, where mirror-image configurations cannot be related by proper rotations within the surface. By dotting the oriented area vector with an externally-defined normal vector to the shell, we obtain sensitivity to PV. We present a discrete estimator for this statistic and also show how to correct for the survey geometry. This parity-odd 3PCF provides a complementary test of PV in galaxy surveys and a potential cross-check of PV 4PCF measurements.

Baryonic feedback redistributes gas around dark matter halos, suppressing the matter power spectrum at scales now probed by weak lensing surveys. X-ray observations directly trace this hot gas, and are one of the main probes of its distribution and properties. We present a forward-modelling framework, built on the baryonification model, linking the three-dimensional gas density and temperature profiles of groups and clusters to observed X-ray surface brightness and luminosity profiles on one side, and to matter power spectrum suppression on the other. We validate the model against independent three-dimensional density reconstructions from the literature, and examine our temperature and metallicity treatment in the group-scale regime. Applying this framework to the SZ-selected CHEX-MATE and X-ray-selected eFEDs samples, we measure gas fractions across the group-to-cluster mass range while accounting for X-ray selection effects, with the first published gas fractions based on CHEX-MATE data. Combining both samples, we derive a joint constraint on the hot gas fraction retained by groups and clusters as a function of mass and on the baryonic suppression of the matter power spectrum. We find $f_{\rm gas} = 0.029 \pm 0.006$ at $M_{500c} = 3\times 10^{13}M_\odot$, $f_{\rm gas} = 0.078 \pm 0.004$ at $M_{500c} = 3\times 10^{14}M_\odot$, and suppression of 6% at $k=1\,h/\rm Mpc$ and 23% at $k=5\,h/\rm Mpc$. Our findings are consistent with recent kinematic Sunyaev-Zel'dovich results, hinting at strong feedback. We also show that the $L_X$-$M$ relation is degenerate with feedback strength, and that different feedback scenarios produce distinct X-ray profile shapes that map onto the same $L_X$-$M$ point. This work is a first step toward extending the framework to forward-model diffuse X-ray emission at the map level for simulation-based inference in upcoming wide-area X-ray surveys such as eROSITA.

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

This commentary briefly reviews recent measurements from DAMPE and LHAASO on cosmic rays to explore why the spectra of different species change shape across the TeV-PeV energy range. DAMPE finds that several primary nuclei show spectral softening at a similar magnetic rigidity, while the combined proton and helium data from DAMPE and LHAASO reveal a more complex evolution at higher energies. By comparing these spectral features and composition changes, we discuss what they may tell us about the sources of Galactic cosmic rays and how these particles propagate through the Galaxy.

This work focuses on the analysis of particle dynamics and the thermodynamic properties associated with a particular branch of hairy black hole solutions. These solutions are obtained by using extended geometric deformation gravitational decoupling on the seed solution of the Schwarzschild solution. This resulting geometry satisfies the dominant energy condition and the condition $Q^{2}=2χ\ell M$. We analyze the motion of massive particles and photons in this geometry and the thermodynamics of this solution. In particular, for massless particles, we study the effective potential and determine the radii of the photon sphere and the shadow of the black hole. For massive particles, we calculate the specific energy and angular momentum of circular orbits, the ISCO, the radiative efficiency of an accretion disk, and representative particle trajectories. Finally, we investigate the thermodynamic properties of this hairy black hole geometry, including the Hawking temperature, the Bekenstein-Hawking entropy, and the heat capacity in the fixed-$(Q,\ell)$ ensemble.

We discuss the cosmological dynamics of general teleparallel quadratic gravity. Using a suitable set of variables, we write the cosmological field equations as a homogeneous dynamical system. We find that this system generically features projective fixed points, where the dynamics can be solved analytically, and exhibits finite time singularities as well as an asymptotically vanishing cosmological expansion.

Magnetohydrodynamical simulations and interferometric observations toward young stellar objects reveal the accretion of gas that may alter the physical properties and chemical composition of protostellar disks. Our goal is to constrain the physical history and properties of late infalling gas, and assess its potential to change the physical and chemical properties of the disk. We carried out a core collapse simulation of 5 M$_{\odot}$ of gas endowed with tracer and sink particles during $\sim 3.6\times 10^{5}$ yr. We analyzed the properties of the rotationally supported disk, and the origin and physical history of the gas from which it is initially formed. We also selected the tracer particles that describe the late infall to analyze their origin and physical history, and compare them to the information obtained from the disk. The final mass of the sink reaches 1.02 M$_\odot$, while the mass of the disk stays fairly constant around 0.02 M$_\odot$. Throughout its evolution, the disk radius increases up to $\sim 40$ au and its density is reduced by a factor of $\sim 8$. Filaments that channel gas from the envelope to the disk and sink appear as a result of the magnetic interchange instability. Tracer particles can identify late infalling gas, whose physical properties resemble that of observed streamers and are only accreted in the final snapshots of the simulation. The origin and physical history of the late infalling gas are different from those of the gas that initially forms the disk. Core collapse simulations predict the appearance of filaments that connect the envelope to the disk. The late accretion of gas through these channels has properties similar to those of streamers and brings gas to the disk with a different physical history that could translate into a different chemical composition.

We present a deconvolution-based framework for testing scale-dependent departures of the late-time matter power spectrum from a fiducial cosmological model using $3\times2$pt measurements. We introduce a free-form scale-dependent modulation $A(k)$ of the fiducial nonlinear matter power spectrum and construct the linear response of binned galaxy-clustering, galaxy-galaxy-lensing, and cosmic-shear spectra to the discretized modulation $A(k)$. The response is evaluated with full-sky, beyond-Limber kernels including density, redshift-space-distortion, gravitational-shear, and intrinsic-alignment contributions. We reconstruct $A(k)$ using a regularized modified Richardson-Lucy algorithm, with weak diffusion in $\ln k$ and selection of the minimum-$χ^2$ solution along the iteration history. Using Rubin/LSST Year 10-like synthetic data, we find that oscillatory modulations with amplitudes $\gtrsim1\%$ can be recovered over $0.1\lesssim k\lesssim0.5\,{\rm Mpc}^{-1}$, provided the oscillation frequency $f\lesssim10$ on $\log_{10}[k/(0.2\,{\rm Mpc}^{-1})]$. We further introduce a posterior-weighted consistency statistic calibrated with posterior-predictive null mocks, thereby accounting for cosmological and nuisance-parameter uncertainties without relying on Wilks' theorem. The null case is consistent with $A(k)=1$, while a $1\%$ oscillatory modulation is detected at $\sim 2.6σ$. These results demonstrate the potential of regularized deconvolution as a model-independent consistency test of the matter power spectrum in future $3\times2$pt surveys.

With microarcsecond precision, \textit{Gaia} has made the astrometric discovery of substellar companions feasible. However, follow-up observations are needed to check on orbital solutions and rule out astrophysical false positives. We validate and characterise a sample of 20 \textit{Gaia} Data Release 3 (DR3) astrometric candidates for substellar companions, with the dual goal of identifying false positives and deriving robust physical parameters for genuine companions. We performed high-resolution spectroscopy using FIES, NIRPS, and NEID at the Nordic Optical Telescope, the ESO 3.6m telescope, and the WIYN Telescope, respectively. Double-lined spectroscopic binaries were identified as false positives, and surviving systems were characterised through joint \textit{Gaia} astrometry and radial velocity modelling. Thirteen of the candidates proved to be near-twin binary stars for which the astrometric motion of the center of light is small enough to mimic that of a single star in response to a substellar companion. We confirm that six companions are brown dwarfs with masses in the range ${\sim}26$--$68\,M_\mathrm{Jup}$, two of them for the first time. In contrast to earlier studies of lower-mass \textit{Gaia} candidates, we find that the astrometric orbital solutions are generally reliable to within the stated uncertainties.

Most stars and planets form in massive star-forming regions, where disks are exposed to external far-ultraviolet (FUV) radiation from nearby O- and B-type stars. The combined effects of stellar irradiation and external FUV fields on terrestrial planet-forming regions (< 10 au) across stellar masses remain unclear. We investigate how internal UV and X-ray irradiation and external FUV fields affect mid-infrared (mid-IR) gas emission and the atmospheric carbon-to-oxygen (C/O) ratios inferred in T Tauri and Herbig Ae/Be disks. We compute disk structures with ProDiMo and synthetic spectra with FLiTs, convolved to a representative JWST/MIRI-MRS resolution (R$\sim$2680). (1) We present the Internal and external irRadIation of diSks (IRIS) grid: four model sets spanning stellar masses of 0.3-4.0 solar masses, including stellar X-ray flares and an external FUV field of 1e4 G0 (Habing units). (2) We predict increasing flux densities with stellar mass for key atomic and molecular mid-IR tracers. (3) External FUV irradiation enhances CH3+ and H2 emission, whereas FUV-induced disk truncation yields inner disk chemistry resembling that of disks irradiated only by their host stars. (4) Mid-IR H2O, CO2, and C2H2 line ratios imply carbon-rich compositions (C/O$\sim$1-10) in the warm emitting layers of T Tauri and Herbig Ae/Be disks, primarily reflecting stellar irradiation, with little sensitivity to external FUV irradiation. The IRIS grid provides a flexible framework for interpreting JWST and future Extremely Large Telescope (ELT) infrared disk observations across a broad range of stellar properties and irradiation conditions. Future models should include FUV-driven photoevaporative winds, X-ray radiative transfer, and time-dependent X-ray irradiation.

Motivated by the singularity theorems of Poincaré Gauge (PG) theory, we investigate extensions of the Holst quadratic model by introducing cubic order invariants constructed from the curvature and torsion tensors into the gravitational action. Such models are characterised by a kinetic structure that is governed by a pseudoscalar mode, whereas the remaining irreducible modes of torsion contribute through nonlinear interactions that can have important implications for the space-time geometry. In particular, in line with other well-known models of PG theory, the Birkhoff theorem does not hold in general, allowing for new exact static and spherically symmetric black hole solutions with dynamical torsion. Across the different torsion sectors, corresponding to the irreducible modes and parity components of the torsion field involved in the analysis, we find both singular and regular configurations. Among the singular solutions, we obtain Kiselev-like and Boulware-Deser-like geometries, as well as new geometries with distinct algebraic and Lambert $W$ metric corrections. In addition, we find regular black holes with both primary and secondary torsion hair, which evade the singularity theorems through violations of the causal convergence conditions induced by the nonlinear torsion interactions. Therefore, we show that Holst-type PG models can support a rich variety of black hole geometries, while providing explicit mechanisms for evading the singularity theorems of PG theory.

Uniform acceleration induces a thermal response to the Minkowski vacuum and can thereby modify quantum correlations. For measurement-induced nonlocality (MIN), a form of nonclassical correlation distinct from entanglement, previous field-mode analyses of bosonic fields found that it decreases with increasing acceleration and vanishes in the infinite-acceleration limit. This raises the question of whether the disappearance of MIN is a generic consequence of the Unruh effect or depends on the physical description of the accelerated quantum system. To address this question, we study two uniformly accelerating Unruh-DeWitt detectors interacting with a massless scalar field. We find that the response of MIN to the Unruh temperature depends sensitively on the initial detector state: it can decrease monotonically, vanish at an intermediate temperature and subsequently recover, or increase monotonically. Moreover, for a broad class of initial states, MIN approaches a nonzero value in the high-temperature limit. These results show that the suppression of MIN under strong acceleration is not universal. Instead, Unruh-induced detector dynamics gives rise to a state-dependent response in which measurement-induced nonlocality can be suppressed, restored, or enhanced.

In this article we provide a systematic reduction to static spherical symmetry for (classical) gauge theoretic models with structure group SU(2) x U(1), coupled to gravity. In particular, we derive the most general static spherically symmetric ansatz for the Einstein-Yang- Mills-Higgs-Dirac-Yukawa equation system. This model is reminiscent of the lepton sector of the Standard Model coupled to gravity. The resulting total system of ordinary differential equations is rather large, so that finding global solutions is non-trivial. For this reason, an open-source interactive solver GUI written in C++ is provided together with the article.

We establish spacetime counterparts of fill-in inequalities for spin initial data sets satisfying the dominant energy condition in dimensions three through seven, controlling the minimum and the integral of the norm of the spacetime mean curvature vector by intrinsic boundary data. The proof relies on a Reilly--MIT estimate for complete spin manifolds with compact boundary. As further consequences, under suitable nonnegativity assumptions, Gromov's hyperspherical-radius fill-in inequality and Bär's recent total mean-curvature estimate extend to complete, possibly noncompact spin manifolds with compact boundary.

This is a concise and compact introduction to a framework for spacetime, gravity, and fundamental physics based on the IKKT matrix model, focusing on conceptual and structural aspects. The approach rests on nontrivial vacua or matrix backgrounds, which serve as 3+1-dimensional noncommutative spacetime. We discuss a specific class of backgrounds describing expanding homogeneous and isotropic FLRW spacetime, which provide a solution of the one-loop effective action with finitely many dof per volume. Perturbations of the background give rise to fields propagating on spacetime governed by a weakly-coupled gauge theory, as well as a nonstandard description of (quantum) gravity leading to some IR modifications. Maximal supersymmetry of the model is essential to ensure effective locality.

Open clusters remain among the most powerful laboratories for calibrating how fundamental stellar properties such as rotation, magnetic activity, surface chemistry, and internal structure evolve over the pre-main-sequence to main-sequence and post-main-sequence stellar lifetime. Because cluster members share a common age, galactic environment, and initial composition, they provide empirical anchors upon which age-dating techniques such as gyrochronology, activity-age relations, and chemical clocks are built. This paper summarizes the Cool Stars 23 splinter session - Open Clusters as Laboratories for Cool Star Evolution, held on 15 June 2026 in Tokyo, Japan. The session comprised one invited review, ten contributed talks, and ten poster pop-up presentations, organized around three themes: the evolution of rotation, magnetic activity, and chemical abundances (Li-depletion, [C/N], [Y/Mg], and neutron-capture elements). We close with a summary of the open questions identified in the panel discussion and the observational and theoretical work that new surveys (e.g. Gaia DR4, 4MOST, WEAVE) and missions (e.g. PLATO, Roman) will enable over the coming years. Together, these efforts will help assess the current state of the field and shape its future directions.

Reconstructing the initial conditions (ICs) of the matter field for an observed volume gives a complete picture of that region's temporal evolution. IC reconstruction is common at low redshifts but largely unexplored during the epoch of reionization (EoR; $z\gtrsim 5$), partly due to difficulties modeling inhomogeneous cosmic radiation fields. Yet the EoR spans half the observable Universe, offering unmatched potential for astrophysics and cosmology. Here we quantify how well upcoming galaxy and 21cm observations can constrain ICs during the EoR. We develop $\tt TILING$ ($\bf T$omographic $\bf I$nference of $\bf L$inear $\bf I$Cs via $\bf N$etwork $\bf G$rafting): a hybrid machine learning pipeline that first produces a point estimate of the ICs, which then improves a score-based diffusion network for generating posterior samples. We train $\tt TILING$ on mock galaxy maps at varying UV magnitude limits, plus corresponding 21cm maps at varying noise levels and foreground "wedge" contamination. For fiducial survey choices, we achieve accurate IC reconstruction (posterior mean cross-correlation coefficients >0.8 and power spectrum errors under a few percent) down to $\text{k}\lesssim 0.2 \text{cMpc}^{-1}$. While 21cm interferometry aids recovery of IC power spectra, most constraining power, especially for IC Fourier phases, comes from galaxy maps, underscoring the need for complementary observations when interpreting the 21cm signal. We show how $\tt TILING$ can recover 21cm power spectrum modes excised by foreground contamination. Our framework can also: (i) guide follow-up observations of sub-volumes of interest; (ii) reconstruct galaxy evolution and reionization morphology for specific volumes; and (iii) isolate the contribution to reionization from the vast majority of galaxies unobservable by optical/IR telescopes like $\it JWST$. Our code is publicly available.

Weak gravitational lensing is a premier cosmological probe, but its small-scale statistical power is compromised by baryonic feedback. Higher-order statistics capture non-Gaussian information that the power spectrum misses, yet their sensitivity to feedback remains a concern for Stage IV surveys. We quantify how unmodeled feedback biases the cosmological parameters inferred from the angular power spectrum (PS), starlet peak counts, and the starlet $\ell_1$-norm, and we determine the scale cuts needed to remove that bias. We also test the Bernardeau-Nishimichi-Taruya (BNT) transform as a strategy for more precise scale cuts. Our analysis is based on the cosmoGRID V1 suite, which imprints feedback on dark matter convergence maps with a baryon correction model, and on simulation-based inference with neural posterior estimation, carried out across footprints ranging from Stage III-like to the full sky. We find that biases grow with survey area, exceeding $2σ$ for all three statistics at Stage IV-like footprints, and removing them costs a substantial fraction of the signal. Restricted to these ``baryon-safe'' scales, the starlet $\ell_1$-norm still reaches a figure of merit almost twice that of the PS. The BNT transform localizes the baryonic sensitivity to the lowest transformed redshift bin and improves the PS figure of merit by a factor of $\sim$1.4, while its linear mixing of shape noise inflates the contours of the map-based higher-order statistics. Higher-order statistics therefore deliver a substantial gain over the power spectrum with no baryonic modeling at all, and an even larger one as modeling improves.

The cosmic microwave background birefringence angle $β$ is exactly degenerate with a common instrumental polarization-angle offset. The Minami--Komatsu(MK) likelihood separates these quantities using Galactic foreground polarization, so its inferred detector angles merit a diagnostic that does not reuse the foreground-$EB$ model. We develop a birefringence-blind differential estimator from antisymmetric and symmetric cross-spectra and apply it for the first time as an explicit consistency test of a foreground-assisted birefringence analysis. For a map pair, the common rotation cancels identically and the estimator measures $Δα_{ij}=α_i-α_j$; combining all pairs in a gauge-fixed network reconstructs the observable detector-dependent calibration pattern. Applied to eight Planck NPIPE detector-set maps, the differential reconstruction agrees with the MK pattern, giving $χ^2=8.18$ for seven degrees of freedom (PTE=0.32), with all shared-data residuals below $1.8σ$. Because both reconstructions use the same maps, this is an independent estimator and modeling route rather than an independent-data confirmation. The agreement therefore strengthens confidence in the calibration foundation of the MK-inferred $β$ without providing an independent absolute measurement of it. A conditional common-mode reconstruction gives $β=0.37\pm0.12^\circ$. The differential test uses existing data and can serve as a validation layer for future $EB$-based birefringence analyses.

Weak gravitational lensing provides a powerful way to map cosmic structure, but most current measurements rely on galaxy shape distortions from deep imaging surveys and are affected by systematics such as intrinsic alignments, photometric-redshift uncertainties and shape-measurement biases. Here we present a spectroscopic galaxy-galaxy lensing magnification measurement using Fundamental-Plane (FP) size residuals, $δ_r\equivΔ\log_{10}R_\mathrm{e}$, from 3.26 million DESI luminous red galaxy (LRG) sources behind DESI Bright Galaxy Survey lenses. The FP-like relation predicts the intrinsic sizes of LRGs from lensing-invariant quantities, including velocity dispersion $σ_0$ and surface brightness $I_\mathrm{e}$, with a scatter of about 0.06-0.07 dex. Lensing magnifies LRG sizes, giving the direct convergence response $δ_r(κ)=κ/\ln 10$, but we show that the full lensing response is modified by magnification bias, because fitting $R_\mathrm{e}$ with $I_\mathrm{e}$ inevitably induces a magnitude dependence in $\barδ_r(m)$. After calibrating this response, we recover surface-density profiles with uncertainties comparable to those from individual Stage-III shear surveys using 5-20 million higher-redshift sources. The corresponding excess surface-density profiles agree with shear-based measurements. We further show that the estimator is robust to size-measurement uncertainties, with a convergence multiplicative bias only $\simeq -0.2$ times the size bias. FP lensing therefore provides a clean, spectroscopic and complementary probe of cosmic structure.

We investigate the non-equilibrium thermodynamics of a Swiss-cheese (SC) braneworld universe. By employing Hayward's unified first law and the Clausius relation at the apparent horizon, we first show that the standard equilibrium treatment of the effective SC braneworld fluid, together with the usual Hawking temperature, reproduces the Bekenstein-Hawking area law without any braneworld correction to the entropy. This demonstrates that the brane quadratic energy density corrections cannot, by themselves, generate a modified entropy area relation within the equilibrium framework. We therefore formulate a non-equilibrium thermodynamic description in which deviations from the Bekenstein-Hawking entropy are accompanied by an internal entropy production term. For a general horizon entropy functional, we derive the corresponding entropy production rate and then specialize to a logarithmic and inverse area quantum-corrected entropy. The resulting irreversible contribution is expressed in terms of the Hubble parameter and its derivative. We apply the formalism to a nonsingular matter-bounce solution in the SC braneworld. During the expanding phase, the total entropy production remains positive satisfying the GSLT. Moreover, the late-time behavior of the second derivative of the total entropy becomes negative, indicating a decelerating entropy growth and an asymptotic approach toward thermodynamic equilibrium. Far from the bounce, the equation of state parameter $w \approx -1$, indicating a dark energy-dominated regime. Significant departures from $w = -1$ are confined to the bounce vicinity, where high energy braneworld corrections drive the nonsingular transition. Thus, the model naturally transitions from a high energy bounce into an asymptotic dark energy phase during late-time expansion.