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

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.

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.

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.

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.

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.