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

Learned surrogates for expensive inner solver blocks are a widely pursued route to faster multiphysics simulation. We report a controlled, end-to-end negative result and identify three structural barriers, none of them a deficiency of the network we trained. The testbed is the implicit Newton solve coupling energy-dependent neutrino radiation to matter in a general-relativistic radiation-hydrodynamics code, its most expensive physics routine per call. First, per-call cost and share of runtime are different quantities, and only the second bounds acceleration. An exclusive self-time profile puts the target block at 16.9% of critical-rank wall clock, capping any surrogate at ~1.2x by Amdahl's law. A surrogate 5.8x cheaper per call merely ties the solver, and the configuration stable enough to run without fallback reaches only parity. Second, offline accuracy cannot rank surrogates for deployment: across fourteen networks the pooled Spearman error-versus-survival correlation (rho = +0.73) is a between-family confound that vanishes under control (rho = -0.04). Third, a correct out-of-distribution gate cannot accelerate a loop that leaves its training distribution. We give the break-even deferral fraction in closed form: because the visited states sit 73x off the data manifold, the gate defers 96.8 to 99.7% of cells, almost invariant to surrogate quality. Including its own cost, the gated loop is a 0.94 to 0.96x slowdown. We further separate stability from fidelity: a never-crashing gated run accumulates a linear -19.9% density bias over 6000 steps. The error is a directed, ballistically accumulating bias, not the variance-driven divergence the autoregressive literature targets.

The large-scale dipole anisotropy observed by the Pierre Auger Observatory above $8$~EeV provides important clues about the origin of ultra-high-energy cosmic rays (UHECRs). In this work, we investigate whether gamma-ray active galactic nuclei (AGNs) from the Fermi-LAT Fourth AGN Catalog (4LAC) can explain the observed dipole and examine the contribution of nearby starburst galaxies (SBGs). Incorporating source fluxes, redshift distributions, GZK attenuation, and a mixed-composition framework, we demonstrate that single-source populations of AGNs and SBGs alone fail to reproduce the observed Auger dipole. To resolve this, we construct a hybrid AGN+SBG model. The best-fit solution for the $8$--$16$~EeV interval requires a $23$\% AGN contribution and $45$\% SBG contribution, which gives a dipole amplitude of $5.03$\% while maintaining the dipole direction of ($104.1^\circ$, $-25.1^\circ$) with angular separation of $7.10^\circ$. In the other energy intervals also, we find that the nearby SBGs increasingly dominate the anisotropic component while AGNs provide a subdominant contribution alongside a quasi-isotropic background.

Electron-capture supernovae (ECSNe) are a promising formation channel for low-mass neutron stars, but the minimum gravitational mass of the neutron stars they produce depends on both the progenitor core structure and the neutron-star matter equation of state (EOS). We compute the electron-capture (EC) threshold gravitational mass ($M^\star_{WD}$) of cold white-dwarf-like O--Ne--Mg cores with representative compositions and map the baryon number onto cold neutron-star configurations constructed from a Bayesian ensemble of unified crust--core EOSs. Although the EC threshold density is sensitive to the concentrations of the O--Ne--Mg mixture, the threshold baryon number of the O--Ne--Mg core varies only weakly, producing a narrow remnant-mass window. In the baseline case with no baryonic mass loss during the transition from the EC threshold mass O--Ne--Mg core to the remnant neutron star, standard ECSNe yield remnant neutron stars with gravitational masses of $1.24$--$1.265\, M_\odot$, with only a small EOS-induced spread. Small baryonic mass losses of $0.01$--$0.02\, M_\odot$ shift this range modestly downward, but the $1.174\, M_\odot$ companion of PSR J0453+1559 would require an extreme mass loss close to $0.10\, M_\odot$, which is not favored by current ECSN simulations. We find that the residual EOS dependence of the remnant mass is controlled mainly by the pressure around nuclear saturation density, while the corresponding tidal deformability remains sensitive to the remnant radius and compactness. Thus, low-mass double neutron star systems can in principle connect ECSN-like formation channels with gravitational-wave constraints on the EOS. Our results show that ECSNe naturally form low-mass neutron stars, but within a restricted mass range; the lightest observed neutron stars likely require low-mass iron-core collapse, ultra-stripped supernovae, or other nonstandard channels.

We report results from a search for an isotropic stochastic gravitational-wave background using data collected by the LIGO--Virgo--KAGRA Collaboration. The analysis uses data from the first observing run through April 1, 2025, during the fourth observing run. New frequency-domain cuts are implemented to address a class of non-stationary spectral noise features that were not effectively identified and mitigated by existing data-quality checks in past analyses. Consequently, previously analyzed data from the fourth observing run are re-processed with the updated cuts. We find no evidence for a stochastic background signal and place upper limits on the gravitational-wave energy density. In particular, for a background following a power law with spectral index 2/3 as predicted by inspiralling compact binaries, we find $Ω_\mathrm{GW}(25\,\mathrm{Hz}) \leq 2.0 \times 10^{-9}$, while scale-invariant backgrounds are constrained to $Ω_\mathrm{GW}(25\,\mathrm{Hz}) \leq 2.8 \times 10^{-9}$, both at the 95\% credible level for a log-uniform prior on $Ω_\mathrm{GW}$. Relative to the constraints from previous data recomputed with the new frequency-domain cuts, these limits improve by a factor of 1.4. We also update bounds on alternative gravity scenarios predicting non-standard polarization modes, and we verify that correlated magnetic noise sources remain below the sensitivity of this search. Combining these observational constraints with population models of compact binary coalescences informed by the latest gravitational-wave transient catalog, GWTC-5.0, we predict the amplitude of the compact binary background to be $Ω_\mathrm{CBC}(25\,\mathrm{Hz}) = 6.3^{+5.0}_{-2.2} \times 10^{-10}$ at the 90\% credible level.

We investigate how cosmological environment regulates the evolution of primordial black hole (PBH) seeds in the early universe using a suite of hydrodynamical simulations. Unlike traditional seeding channels, PBHs provide early BH seeds without the need of special conditions for gas collapse/fragmentation that favor overdense regions, thus covering diverse large-scale environments. We find that PBHs follow distinct evolutionary pathways depending on the gas supply and halo assembly history. In underdense regions, limited inflow suppresses both accretion and star formation, producing faint, metal-poor systems, while in overdense environments, sustained gas inflow drives rapid BH growth and the formation of compact, centrally concentrated stellar components. These differences lead to large variations in BH-to-stellar mass ratio, metallicity, and morphology, collectively bracketing a range of possible evolutionary pathways for PBH-seeded systems. We show that compact, BH-dominated sources resembling recently observed Little Red Dots naturally arise as one phase within these evolutionary pathways before evolving into more extended galaxy--AGN systems. Our results suggest that both the initial seed properties and cosmological environment jointly shape early BH--galaxy co-evolution, while subsequent environmental regulation can erase the memory of the initial seeding channel.

Periodic density structures (PDS) are trains of plasma enhancements carried into the heliosphere, but their repeated timing does not by itself identify where the cadence is set. We test whether PDS selected independently in the outer corona can be traced back to a low-coronal clock, and whether that organization survives through the EUVI--COR1--COR2 observing chain. We analyze 12 STEREO-A/SECCHI events on 2008 January 11--14 along a nonradial path. EUVI 171 A at 1.10--1.20 R_sun shows an unusual concentration of power in the predeclared 80--130 minute band relative to filter-matched red-noise controls, while the ensemble is most strongly organized in upper COR1 at 2.5--3.0 R_sun. These signatures identify a low-coronal modulation candidate and an intermediate-height organization domain, but they do not form a unique phase-preserving clock extending into COR2. Instead, the more persistent observable is spatial order. Events 9 and 12 retain expansion-stable outward ordering. In polar r--theta maps, intersections of oppositely inclined, expansion-aware matched-filter ridge supports form repeated X-/diamond-like patterns. The morphology is compatible with stationary or quasi-stationary shock-cell processing, although the brightness diagnostics are insufficient to establish a stationary MHD shock branch. The observations therefore favor an intermittent source--gate--transfer picture: a low-coronal cadence may modulate plasma release, while event-dependent propagation progressively destroys exact phase coherence. The outward ordering of the density structures can survive, providing a more persistent signature of the source-to-wind transfer than phase locking itself.

We present phase-connected timing and polarization measurements for two long-period FAST-CRAFTS pulsars, PSRs J0000+6252 and J2131+3642, and extend single-pulse emission-state analysis to a five-source sample including PSRs J1903+1407, J1502+4653, and J2112+4058. FAST L-band timing baselines span 414 to 511 days; the two pulsars have spin periods of 1.11 to 1.55 s, period derivatives of $(3.99$ to $4.06)\times10^{-15}~{\rm s~s^{-1}}$, characteristic ages of 4.34 to 6.16 Myr, surface magnetic fields of $(2.15$ to $2.52)\times10^{12}$ G, and spin-down luminosities of $(4.21\times10^{31}$ to $1.16\times10^{32})~{\rm erg~s^{-1}}$. Their rotation measures are $70.2\pm26.7$ and $-44.6\pm7.6~{\rm rad~m^{-2}}$, with linear polarization fractions of 17.6\% to 27.2\%. PSR J2131+3642 shows a short monotonic PA segment permitting a formal rotating-vector-model fit, though limited longitude coverage leaves the geometric parameters poorly constrained; PSR J0000+6252 has too few PA points for such a fit. Gaussian mixture modeling (GMM) of single-pulse energy distributions identifies null, weak, and burst components in PSRs J0000+6252, J1903+1407, J1502+4653, and J2112+4058, while J2131+3642 shows only weak and burst states. We also identify bright single pulses (peak intensity $\geq10\times$ the integrated average profile): 42 in J0000+6252, four each in J1903+1407 and J2112+4058, and none in J2131+3642 or J1502+4653. These bright pulses occur within the main emission window with no evidence of periodic recurrence, consistent with sporadic enhancements of the normal radio-emission beam. For J2112+4058, we measure a scattering timescale $τ_{\rm sc}=5.84\pm0.18$ ms at $ν_{\rm ref}=1.25$ GHz. Together, these results highlight the diversity of magnetospheric variability among slowly rotating neutron stars.

The Solar wind is a continuous flow of charged particles emanating from the solar surface and governed by complex, interacting magnetohydrodynamic processes. Accurate specification of inner-boundary conditions is essential for heliospheric modeling and solar-wind prediction. In many practical applications, only a subset of interacting multi-field variables is directly available, but for a comprehensive view of solar wind prediction and downstream magnetohydrodynamic simulations, a more complete boundary state is required. In this work, we study the problem of learning the multi-field multi-scale solar magnetohydrodynamic state at 30 solar radii ($R_\odot$) using operator learning. Specifically, given the radial velocity and radial magnetic field, we aim to reconstruct the non-radial velocity and magnetic field components, radial and non-radial current density, thermodynamic density, and pressure components. This mapping is highly nonlinear, spatially coupled, and multi-scale, making it a challenging task for data-driven scientific machine learning. To address this problem, we employ a Local Neural Operator (LocalNO) that learns mappings between input and output function spaces while retaining locality and resolution-awareness. Unlike conventional regression models and autoencoder models, neural operators are better suited for learning structured field-to-field transformations arising from physical systems. The resulting predictions along with inputs are intended to serve as boundary condition variables for future inner-heliospheric modeling pipelines.

Cosmological first-order phase transitions (FOPTs) in the early Universe are an exciting prediction of many beyond the Standard Model scenarios. Owing to the stochastic nature of bubble nucleation, some causal patches may remain trapped in the false vacuum long after the surrounding regions have completed the transition and started to redshift. Under suitable conditions, these delayed patches can become sufficiently overdense to collapse into primordial black holes (PBHs). One therefore expects characteristic parameters, such as the PBH mass and PBH formation time to be sensitive to the underlying phase-transition dynamics. Here, we adopt a model-independent framework with a fixed mean bubble separation scale ($R_*$) to directly compare the impact of exponential and Gaussian nucleation profiles on PBH formation. We further incorporate the effects of friction on the bubble walls arising from interactions with the surrounding plasma. We find that friction significantly modifies the bubble dynamics, leading to significant changes in the PBH mass and formation time for both nucleation histories. Moreover, the stochastic gravitational-wave spectrum generated by the FOPT can provide a complementary probe of the underlying dynamics, allowing the effects of friction on the bubble wall evolution to be distinguished from the frictionless case.

We present the capabilities of the Fluorescence detector Array of Single-pixel Telescopes (FAST) observatory in the single telescope configuration as an important step towards the possible future large-field observatory for detecting ultra-high-energy cosmic rays. Reconstruction of main shower physics parameters are explored on noise-free simulated events using machine learning techniques in the challenging domain of a low-intensity transient signal. We find a very good correlation between the true and reconstructed energy of the shower even with the information from just the four photomultipliers of the single FAST telescope, and a reduced performance for the maximum of the shower development Xmax, using various architectures of artificial deep and convolutional neural networks, with a comparison to a benchmark gradient boost regression model. The resolution in the energy is found at the sub-percent level, while in Xmax it is~$5\%$. The relative difference between predicted and true values is under one percent for energy, while for Xmax it ranges from $-8\%$ to $+17\%$, which can be attributed to the limited information from the single FAST telescope configuration. The results constitute an important capabilities verification and a lesson learned with implications for established FAST prototypes as well as for more complex configurations of the FAST observatory under construction.

Background: Localized structure in the equilibrium squared sound speed can affect neutron-star observables, but modifying $dP/dε$ requires thermodynamically consistent reconstruction and assessment over the complete declared domain. Purpose: We quantify the response to positive and negative localized sound-speed changes around a unified BSk24 baseline while keeping the deformation geometry, admission criteria, reconstruction, and fixed-mass comparison protocol fixed. Method: A Gaussian profile with quintic smootherstep activation is added to $c_{\rm s}^2$. Proposals violating $P>0$ or $0<c_{\rm s}^2\leq1$ anywhere in the retained domain are rejected before reconstruction. Passing cases are reconstructed as effective cold one-fluid barotropes and evolved with the Tolman-Oppenheimer-Volkoff and tidal equations. Results: The $A=0$ case reproduces the undeformed control exactly. A localized change produces a persistent pressure offset, relocates the fixed-mass central state, and yields a mass-dependent, nonlinear, sign-asymmetric response. For $A=-0.8$ and $+0.8$, respectively, $ΔΛ/Λ_0=(-46.02\%,+35.33\%)$ at $1.4\,M_\odot$ and $(-63.52\%,+87.17\%)$ at $2.0\,M_\odot$. The reached stiffened feature remains core connected at the reporting endpoints, whereas the softened feature becomes an off-center shell at high mass. Conclusions: These phenomenological BSk24-anchored effective barotropes illustrate how a local sound-speed intervention propagates through thermodynamic reconstruction and stellar structure. They are not BSk24 predictions or evidence for a microscopic composition or phase transition.

While no signal originating at a black hole's event horizon can reach future null infinity, information about the horizon and its immediate vicinity can be encoded in asymptotic properties of waves emitted by matter or field perturbations falling toward a growing/forming horizon. Within a response-filtered framework, "direct wave" denotes source-sensitive plunge and remnant-formation information revealed by filtering the black-hole response. Unlike a stationary damped sinusoid with a fixed complex frequency, the direct wave follows the evolving source and has an evolving instantaneous complex frequency tied at late times to the remnant horizon angular velocity $Ω_H$ and surface gravity $κ_H$. Using rational filters, we remove quasinormal modes from numerical-relativity waveforms to study this evolution. Before numerical contamination, trajectories depend on remnant spin: the real frequency evolves toward $2Ω_H$ from above for lower spins ($χ_f\lesssim0.7$) and from below for higher spins ($χ_f\gtrsim0.7$), while the instantaneous decay rate increases toward $\sim2κ_H$ and, in some high-spin cases, beyond it. As in particle-plunge results, the horizon-controlled value is approached only at late times, as frame dragging controls near-horizon motion. For $χ_f\sim0.7$ remnants of non-precessing, comparable-mass binaries, the early-time real frequency is close to $2Ω_H$ because the binary orbital frequency transitions smoothly to the remnant horizon frequency. Finite-time deviations therefore carry information about merger/collapse dynamics rather than undermining the horizon connection. Full direct-wave evolution requires numerical-relativity calibration. We further show that approximate pole-zero pairing in the Kerr response motivates a spin-dependent minimal filter set that suppresses quasinormal-mode features while revealing source-trajectory information.

We present a detailed observational analysis of BLG723.31.34651, the fourth reported eclipsing member of the class of Roche lobe-filling hot subdwarfs, using time series spectroscopy and multi-band photometric observations from the Goodman High Throughput Spectrograph on the 4.1-m SOAR telescope. Our light curves show periodic photometric variability at $P=64~{\rm min}$ from the tidally-deformed hot subdwarf donor and deep eclipses of the donor by an opaque accretion disc. We report two light curve model solutions that agree with each other to within $2σ$, both suggesting that BLG723.31.34651 contains a typical hot subdwarf donor, a white dwarf accretor, and an irradiated opaque accretion disc. Our spectroscopic analysis identifies statistically significant phase-dependent variations in the effective temperature and surface gravity coinciding with the phase-dependent variations seen in our observed light curves. We attribute these variations to a combination of intrinsic variations in local atmospheric parameters due to tidal deformation of the donor and unmodeled spectroscopic flux contributions from an accretion disc.

Atmospheric photochemical models incorporating the impacts of stellar flares often assume a $\sim$9,000 K spectrum at ultraviolet-optical wavelengths. Recent multiwavelength observations, however, reveal a more complex picture with temperature measurements spanning 4,000-40,000 K, although the occurrence rates for flares with different temperatures remain unknown. Here, we model the evolution of a Proterozoic Earth-like world with 0.01 bar of O$_2$ under repeated flaring to identify the impact of flare effective temperatures. We explore four scenarios - two host star types (K2V and M2.5V) and two flare temperatures (9,000 K and 19,000 K) - selected to bound the potential parameter space. The hotter flares have a larger impact on O$_3$ photochemistry for both stellar types. M-star planetary atmospheres are more volatile and exhibit rapid changes in their O$_3$ production and destruction rates. Meanwhile, K-star planetary atmospheres are more stable and are only impacted by the hottest flares, proving advantageous for biosignature searches. We simulate 0.2-1.0 $μ$m reflected light spectra for all four scenarios, and find that 19,000 K flares can result in either production or destruction of O$_3$ depending on the host star spectral type increasing the 0.2 $μ$m feature by $\sim$2$\times$ for the K2V star but decreasing it by 50% for the M2.5V star. Future missions such as the EVE SMEX mission concept will provide robust flare temperature constraints for young FGKM stars, which will serve as inputs to improve photochemical models to inform future HWO observations.

AU-scale structures in the interstellar medium have become relevant across several areas of astrophysics, from ISM microphysics to pulsar timing arrays and cosmic-ray (CR) propagation. Current sheets formed at magnetic-field reversals have drawn particular attention because they can deflect the dominant GeV CR population through large angles. Pulsar scintillation, an outstanding probe of plasma density, resolves these structures---but not their magnetic configuration. Conventional Faraday rotation measures (RM) are only marginally sensitive at AU scales and are further compromised by ionospheric systematics. We propose a method that isolates birefringence-induced polarization---the signature of magnetic reversal---by applying phase retrieval to polarized data. Pulsar scintillation further resolves ray paths separated by $\sim 1\,\rm{AU}$, enabling differential RM measurements that suppress the ionospheric common mode while preserving the AU-scale magnetic signal. We demonstrate the method on an archival observation of PSR B0834+06. Phase retrieval resolves the 1-ms scattering feature into two branches in Doppler-delay space. Under the corrugated current-sheet interpretation, these branches correspond to ray paths sampling opposite sides of an AU-scale current sheet. We measure a branch-to-branch RM difference of $(-9.3\pm3.2)\times 10^{-3}\text{ rad m}^{-2}$ at 2.9$σ$ across a 10-MHz band. The signal is consistent with a field reversal of $|Δ\langle B_\parallel \rangle| \simeq4.4\pm2.1\,μ\rm G$, sufficient to sustain a long-lived sheet. With only six independent subbands, we present this as a proof of concept rather than a definitive detection. If it proves scalable, its full implementation could constrain cosmic-ray transport and inform models of propagation-induced noise in pulsar timing arrays.

Recent observations have revealed gas-rich debris disks around intermediate-mass stars at ages of tens of Myr. The origin of this gas remains unclear: it may be primordial, retained from the protoplanetary phase, or secondary, released from volatile-rich solids. Secondary-origin models reproduce CO emission but often overpredict neutral carbon. Recent observations and disk-evolution models suggest that primordial gas may survive longer than previously assumed, motivating thermochemical tests of the primordial-remnant scenario. We test the previously unexplored possibility that primordial-origin disks satisfy the observational constraints on gas-rich debris disks. Specifically, we determine under what conditions a disk around a $2\,M_{\odot}$ star reproduces substantial CO, low CI/CO ratios, and weak HCO+ emission consistent with current non-detections. We post-processed 20-40 Myr structures from 1D disk-evolution models with Cloudy, varying irradiation geometry, dust-to-gas mass ratio (DTG), and cosmic-ray ionisation rate. In the dust-poor models (DTG $=10^{-4}$), CO remains optically thick around $R\sim100$~au. The models yield low disk-integrated CI/CO mass ratios. Our model produces CO radial intensities of the observed order of magnitude, but its CI-emitting region is more extended than observed. The standard CR model overproduces HCO+, whereas the weak CR model brings its predicted luminosity within current observational limits. These results demonstrate that a primordial origin remains chemically viable for CO-rich debris disks. The main remaining tension is the excessive radial extent of the CI emission, although it may reflect our simplified modelling. Further testing of the primordial-origin scenario will require multidimensional, self-consistent modelling, spatially resolved CI observations, and deeper searches for HCO+.

We provide a geometric foundation for the kinetic theory of photons. Although the induced metric $\hat h$ on the light cone bundle $Γ_0^+$ is degenerate, which causes the corresponding volume element to vanish, we can still use a method similar to the Hodge dual to construct a volume element $η_{Γ_0^+}$ for the light cone bundle. Based on this geometric structure $(Γ^+_0,η_{Γ_0^+},\hat h)$, we establish the fully covariant Boltzmann equation for photons. More importantly, the volume element and the induced metric are linked in a nontrivial way, which allows the physical distributions to be defined consistently. This yields the corresponding hydrodynamic quantities and their divergences, which take the same form as in the case of massive particles.

We study cosmological perturbations in multi-field inflation in which the scalar fields couple to the Gauss-Bonnet terms through a general coupling function $f(φ^a)$. We derive the complete quadratic action for the field perturbations in spatially flat gauge and find that, besides corrections to the kinetic, gradient and mass matrices, the Gauss-Bonnet coupling induces an antisymmetric velocity coupling that vanishes in Einstein gravity. Decomposing the field perturbations into curvature and isocurvature modes, we show that the effective mass of the curvature perturbation and the non-derivative mixing mass vanish identically, $\mathcal{M}_{\mathcal{R}}^2 = \mathcal{M}_{\rm mix}^2 =0$ which are justified by the Weinberg's adiabatic mode. We then compute the curvature power spectrum in two limits. First, when the isocurvature mode is heavy, it can be integrated over all scales. This yields an effective single-field theory of the curvature perturbation with a modified sound speed. Second, when the gradient of the coupling function is aligned with the background trajectory ($f_N =0$), the derivative mixings disappear and the superhorizon transfer of isocurvature into curvature perturbations is governed by the turn rate, which is modified by the Gauss-Bonnet corrections. We also present the tensor power spectrum and the tensor-to-scalar ratio at linear order in the coupling.

In the present paper, the cases of a uniformly accelerated reference frame and a Schwarzschild black hole in two and four dimensions are discussed focusing on formal aspects of the corresponding quantum scalar field theories. The main emphasis is made on building consistent quantum field theories and highlighting the problems that are not considered elsewhere. In particular, some pathologies in the case of a uniformly accelerated reference frame are discussed in detail: it is shown that some quantum states that are well defined in Minkowski spacetime possess an infinite energy when considered in Rindler spacetime. In the case of a Schwarzschild black hole in four dimensions, it is shown that the standard theory providing the Unruh vacuum state turns out to be inconsistent. The problem of consistency of quantum scalar field theory in Kruskal-Szekeres spacetime is also discussed.

The central challenge in microscopic quantum-gravity models is to derive geometric and gravitational observables from one refining quantum family without changing dynamics or fitting a new coefficient at each interface. We construct a controlled benchmark on a declared fixed physical Fourier band. For every spatial regulator, one finite-Hilbert interacting Hamiltonian defines a unitary completely positive process and a Euclidean Schwinger functional. Its spectrum and source derivatives determine the Lorentzian gap, connected four-response, mixed matter--geometry response, dynamic geometry kernel, induced Newton coefficient, and quantum-matter stress entering a semiclassical FLRW evolution. Two inequivalent spatial symbols obey analytic second- and fourth-order bounds and converge to the same limiting Hamiltonian. Duhamel and Weyl inequalities, a uniform gap, analytic perturbation theory, and a kernel-remainder estimate yield a common fixed-band regulator limit. Within the declared local covariant FLRW two-derivative sector, the metric source uniquely fixes a positive induced Newton coefficient rather than leaving it as an interface fit. Direct diagonalization independently yields a nonzero connected fourth response and a non-fitted dimensionless spectral combination. The excitation gap supplies a conserved stress tensor for a semiclassical Friedmann trajectory. This is a same-family fixed-band closure theorem, not an all-band QFT or completed quantum gravity. The two-mode band, finite oscillator cutoff, covariant two-derivative sector, and semiclassical Einstein equation are inputs. It does not construct a quantum metric Hilbert space, quantum BV measure, topology-change control, or autonomous sector selection. The result provides a reproducible benchmark and isolates the further theorems needed for nonperturbative closure.

Primordial non-Gaussianities are at their largest when the curvature perturbation is coupled to a massive field by a mixing too strong to be treated perturbatively. We show that this regime, naturally reached in multifield inflation, is solved analytically by the Laplace-space representation we recently introduced, in which every mode becomes a superposition of plane waves weighted by a single kernel. The bispectrum then follows as a rapidly convergent series of elementary functions, in perfect agreement with independent numerical results. We find that the cosmological collider signal at strong mixing is larger than widely believed, with only half of the usual Boltzmann suppression. Its amplitude further grows with the number of time derivatives of the curvature perturbation in the interaction, singling out the most promising observational targets.

The inner disk truncation radius, $R_T$, plays a crucial role in the regulation of star-disk interaction and the early evolution of star-disk-planet systems; however, measuring this parameter is observationally challenging. We present a new method for determining $R_T$ in young accreting systems that hinges on the color dependence of the accretion shock emission in multi-band time-domain surveys. Based on the accretion simulations of Robinson et al. (2017, 2021), we produce synthetic color-magnitude diagrams at near-UV and optical wavelengths that predict the loci of accreting stars as a function of $R_T$. We test these model predictions on young stars with interferometric $R_T$ estimates, finding very good agreement in our results. We apply this novel technique to a pilot survey of 26 classical T Tauri stars in Taurus and Upper Scorpius. We find a predominance of sources with small truncation radii, $R_T < 4\ R_\star$, and an overall distribution of $R_T$ that is statistically similar to that inferred from interferometric studies, while differing from those inferred from emission line modeling. Finally, we discuss the application of this technique to NASA's mission concept EVE, with the goal to provide simultaneous measurements of inner disk truncation radii, corotation radii and mass accretion rates for hundreds of young stars across the Galaxy. The unprecedented survey of inner disk properties that the mission would produce would enable the first stringent test of angular momentum evolution theories in young stars and reveal the impact of the inner disk conditions on early planet architectures.

Light spectator fields during inflation can acquire superhorizon fluctuations that cross a potential barrier between positive and negative regions of their potential. Motivated by the Standard Model Higgs instability, in this work we study the subsequent evolution of patches where this occurs in the radiation era after inflation ends for a Higgs-like spectator field. We utilise fully nonlinear, spherically symmetric numerical relativity. Across the black hole forming configurations in our investigation we find a robust two-stage evolution. First, the central negative potential region reverses its expansion, becomes kinetic dominated, and forms a primordial black hole that hides the runaway core. The positive potential barrier that survives outside this first horizon then determines one of two late-time branches. In the subcritical branch the original apparent horizon grows smoothly and engulfs the remaining scalar structure. In the supercritical branch, however, the potential energy of the barrier dominates the local evolution. The result is a transient wormhole throat, a bifurcating trapping horizon, and an inflating child universe branch. In both branches the parent radiation dominated universe is ultimately left with an ordinary primordial black hole whose subsequent growth is governed by radiation accretion.

We calculate the dynamical tidal response coefficients (TRCs) of 2+1D analogue black and white holes generated by draining and fountaining bathtub flows, respectively. The parameter space is characterized by the frequency and azimuthal number of the Fourier modes and the rotation of the analogue black or white hole. In general, the TRC is a complex-valued function of these parameters. Its real and imaginary parts are defined as the tidal Love number (TLN) and the tidal dissipation coefficient (TDC), respectively. The TRC of the analogue black hole (ABH) exhibits several interesting features. At certain points in the parameter space, including the static case of vanishing perturbation frequency, the TLN exhibits logarithmic running, while the TDC vanishes. Unlike in Einstein gravity, fluid dynamics allows for the physical existence of an analogue white hole (AWH) event horizon. Time-independent, torque-free, barotropic, inviscid, axisymmetric fluid dynamical equations can yield a pair of transonic background flow solutions. The solutions in the pair, corresponding to an ABH-AWH pair, share the same angular momentum and Bernoulli constant but have opposite mass flow rates, all of which are conserved quantities. For such an ABH-AWH pair, the TRC of the AWH is the complex conjugate of the TRC of the ABH.

Massive stars characteristically form in clustered environments. Characterising young massive 'protoclusters' is therefore crucial to constraining the mechanism(s) of massive star formation, and of the assembly of stellar clusters. We present 1.3 mm continuum results from the ALMA EGO-10 imaging survey, targeting ten Spitzer GLIMPSE Extended Green Objects (EGOs) - massive protostars with active outflows traced by extended 4.5 $μ$m emission. Our sensitive 1'.6$\times$1'.6 mosaics reveal rich protoclusters associated with all targets. With a mean spatial resolution 2200$\times$1600 AU, we identify 570 cores - between 13 and 135 per field. We quantify protocluster structure with the $Q$-parameter, finding structural diversity with 0.5 $\lesssim Q \lesssim$ 0.9. The sample is notable for the wealth of complementary high-resolution multiwavelength data available. Correlating our cores with these observations, we find only 2%, 5% and 4% of cores host 6.7 GHz CH$_3$OH masers, 22 GHz H$_2$O masers and cm-$λ$ continuum sources, respectively. The massive protostars traced by 6.7 GHz masers typically reside near protocluster centres (median offset 0.045 pc), and all at $d<$ 3 kpc are found in clustered locales, with $>$10 cores within 10,000 AU. Using VLA cm-$λ$ continuum observations, we construct a new evolutionary indicator: the ratio of protocluster cm-$λ$ continuum luminosity to the mass of the associated ATLASGAL clump ($L_\text{cm}/M_\text{AGAL}$). This ratio correlates positively with $Q$, with the correlation driven primarily by the cm-$λ$ continuum emission from MYSOs. This suggests dynamic protocluster structure, evolving from subclustered to centrally condensed, consistent with the global collapse in hierarchical, clump-fed models of massive star formation.

The Euler product of the Riemann zeta function is the partition function of a free bosonic gas whose mode energies are the logarithms of the primes. We show that this arithmetic gas, combined with the assumption that the entropy--geometry correspondence holds, leads to a quantum-corrected black-hole metric. The prime gas is a Hagedorn system. Its entropy is linear in the energy, which is exactly what an entropy linear in the horizon area requires, and the simple pole of the zeta function at $β=1$ fixes the coefficient of the logarithmic correction to the area law. The nontrivial zeros cannot play this role: their level density grows only logarithmically, and far too slowly to be extensive. Demanding that the reconstructed geometry reduce to Schwarzschild at large radius then fixes the map from arithmetic energy to horizon area and yields the closed-form metric $f(r)=1-2GMr/(r^{2}+\lz^{2})$ with $\lz^{2}=αG/π$. This describes a two-horizon black hole with Reissner--Nordström horizon structure but no Coulombic hair, a positive-energy anisotropic source obeying the null energy condition, a softened central singularity, a bounded Hawking temperature, and a cold extremal remnant that ends the evaporation. The same length scale follows independently from requiring that the first law hold exactly with the corrected entropy. The nontrivial zeros survive only as exponentially suppressed log-periodic ripples in the area, which suggests a physical interpretation of the Riemann hypothesis as the statement that arithmetic corrections to black-hole thermodynamics are as small as they can be.

Strong gravitational lensing contains key information about cosmic acceleration. Modern and next-generation galaxy imaging surveys are expected to provide high-quality data on $\mathcal{O}(10^5)$ galaxy-galaxy lensing systems. The plethora and complexity of the data are likely to present computational challenges for parameter inference methods for fitting high-dimensional likelihoods, which are often analytically intractable. Neural Ratio Estimation (NRE) efficiently computes individual likelihood ratios that can be combined into population-level posteriors. We use simulations to study the capacity of NRE to jointly predict the dark energy equation-of-state parameter $w$ and the total matter density $Ω_{m}$ from lensing images and companion spectroscopic information. We also introduce a post hoc posterior coverage calibration procedure that mitigates the model overconfidence that is typically found in neural density estimation applications. Our experiments show that the errors on both parameters decrease with increasing inference population sizes. In particular, for 100 lenses in a standard $Λ$CDM Universe, our calibrated NRE model achieves median fractional uncertainty of $22.8\%$ in $w$ and $2.9\%$ in $Ω_{m}$. This proof of concept demonstrates a potentially scalable approach for efficient cosmological parameter inference with large populations of galaxy-scale lenses observed in future surveys.

We study the one-loop relation between an effective-field-theory extension of Poincaré gauge theories of gravity with propagating torsion and its unimodular counterpart. In the Einstein representation, we consider two representative quadratic torsion subsectors, namely totally antisymmetric (axial) torsion and the vector component of hook-antisymmetric torsion. On maximally symmetric metric backgrounds with vanishing background torsion, the diffeomorphism-invariant and genuine unimodular theories yield identical one-loop determinants in both sectors, reproducing the same local logarithmic divergences. We further examine flat backgrounds with homogeneous axial or vector torsion, where metric-torsion mixing is present. In each case, the torsion-dependent one-loop effective action in unimodular gauge and in the genuine unimodular formulation is controlled by the same physical determinant. Thus, for the actions and backgrounds considered, the unimodular constraint does not modify local one-loop torsion dynamics. We also delineate the limitations of this equivalence and the extensions needed for generic torsionful curved backgrounds.

We investigate the Joule-Thomson expansion of a non-commutative Schwarzschild-AdS black hole, whose space-time geometry is described by the Lie-algebraic $κ$-deformed space-time, and show that the $κ$-deformation parameter induces Joule-Thomson expansion in an uncharged-AdS black hole. We find the ratio of minimum inversion temperature to critical temperature $T_i^{(min)}/T_c \simeq 0.493$, is independent of the deformation parameter, and is in close agreement with the value for Reissner-Nordstrom black hole. Furthermore, we show the existence of a minimum black hole mass (depending on the $κ$-deformation parameter) below which the Joule-Thomson expansion vanishes.

At a deeper level the principle of least action is interpreted as the law of least entropy increase consistent with Prigogine's principle of minimum entropy production, and the implications of quasistatic information quantization rule (Proc. R. Soc. A 480: 20240024), are explored for the conjectured second law of quantum complexity. It is thus shown that the conjectured second law of complexity is derivable from the information quantization rule such that long after heat-death the quantum state complexity evolves as $C(t)=C_{\rm max}\exp(-1/t)$, increasing with time to a saturation value $C_{\rm max}$ that is exponential in the equilibrium entropy $S_{\rm max}$. For the out-of-equilibrium circumstances, however, the quantum complexity can decrease with the time, asymptotically tending to a minimum determined by the distance from equilibrium.

Gravitational-decoupling provides a systematic framework for constructing black hole geometries from known (seed) solutions (such as the Schwarzschild), whose exterior properties closely resemble those of the seed solution while still retaining potentially distinguishable features. Using this approach, we study a class of static, spherically symmetric hairy black holes and examine how such deformations affect their horizon structure, photon spheres, shadow scales, and scalar quasinormal modes (QNMs). We first determine the physically controlled regions of the parameter space by imposing the relevant horizon and energy condition requirements, while retaining configurations outside these domains only for exploratory comparison. To isolate geometric effects, we compare solutions at a common horizon radius, while accounting explicitly for the distinction between the mass parameter entering the seed geometry and the ADM mass measured asymptotically. The resulting ADM normalized shadow scales are used for an illustrative comparison with the characteristic angular scales of M87* and Sgr A*. We also compute massless scalar QNMs using the sixth-order Wentzel-Kramers-Brillouin approximation and the improved Asymptotic Iteration Method, finding close agreement over the modes considered. Our results show that fixing the horizon radius does not uniquely determine either the optical or perturbative properties of the spacetime, and that mass normalization can significantly alter the interpretation of geometric trends. Horizon structure, null geodesic observables, and scalar QNMs therefore provide complementary probes of gravitationally-decoupled black hole geometries.

The Sun is one of the fundamental benchmarks in astronomy, and there is ever growing interest in precise and accurate determinations of its chemical composition. We present the current state-of-the-art in spectroscopic determinations of the solar elemental abundances, that now routinely employ three-dimensional (3D) radiative-hydrodynamic simulations of the solar photosphere and post-processing in non-local thermodynamic equilibrium (non-LTE). We critically review the recent literature and present recommended present-day and protosolar abundances of the 83 long-lived elements, and discuss their implications vis-à-vis primitive meteorites, helioseismology, neutrino fluxes, and the Solar Modelling Problem.

As a natural extension of the single-field fuzzy dark matter (FDM) model, multifield FDM has attracted increasing attention in recent years. This scenario is motivated both by the axiverse scenario predicted by string theory and by the possibility that multifield FDM may provide a better match to astrophysical observations than its single-field counterpart. In this work, we perform high-resolution numerical simulations to systematically investigate the dynamics of wave structures in multifield FDM halos. In particular, we study the oscillatory and stochastic motions of the central core, the evolution and statistical properties of granules, and the resulting dynamical heating of embedded stellar systems. We find that the frequency spectra of the core density oscillations develop multiple peaks and shift toward higher frequencies relative to the single-field case. The centers of different field components undergo nearly synchronized random walks, while subdominant components exhibit larger random-walk amplitudes. We further show that the suppression of granule density fluctuations with increasing number of fields is largely insensitive to the fractional abundance of each component over a broad parameter range. Moreover, using self-consistent simulations, we find that the dynamical heating induced by granules is progressively suppressed as the number of fields increases. However, once the contribution from the central core is taken into account, this trend would become much less pronounced.

Motivated by the old idea of treating gravity as a gauge theory, we review the gauge-theoretic construction of Conformal Gravity and of Noncommutative (Fuzzy) Gravity. We then outline how these gauge formulations of gravity can be unified with internal interactions described by an SO(10) Grand Unified Theory, based on the critical observation that the dimension of the tangent-space symmetry group of a curved manifold need not coincide with the spacetime dimension.

The deflection angles of rays diverge logarithmically near outside and inside of a photon sphere and they could be detected by near-future space observations on black hole shadows. We show that the deflection angles of the rays near a photon sphere at the throat of an example among some commonly used types of wormholes or one of the simple wormholes, which is often studied by the Event Horizon Telescope collaborations, does diverge not logarithmically but in power. This is because it has a marginally unstable photon sphere at a $Z_2$-symmetric wormhole throat and the wormhole has not a common nor simple but a special property. We investigate numerically gravitational lensing of rays with the deflection angles diverging in power outside and inside of a marginally unstable photon sphere at $Z_2$-symmetric wormhole throat in strong deflection limits. We apply our numerical method to the commonly-used simple wormhole spacetime, a Damour-Solodukhin wormhole spacetime, and a Simpson-Visser black-bounce spacetime. Our numerical results imply universal property of the deflection angle by the marginally unstable photon sphere on $Z_2$-symmetric wormhole throat in the strong deflection limits as similar to the deflection angles of ray bent by a marginally unstable photon sphere off a throat in the strong deflection limits which have a different power. We also correct the coefficients of power-divergent terms of deflection angles of the rays bent nearly outside of the marginally unstable photon sphere at the $Z_2$-symmetric wormhole throat in a semi-analytic approach investigated by the author and others.

We prove the black hole uniqueness conjecture in the axially symmetric, stationary, vacuum setting: there is no regular asymptotically flat configuration with more than one horizon component when every degenerate component has nonzero angular momentum. More precisely, it is shown that in any such multi-horizon configuration, the logarithmic angle defect along every finite axis rod is strictly negative, and hence the interaction force along such axis rods is always attractive. The proof is based on a refined asymptotic analysis of the associated singular harmonic maps, in part from the companion paper [16], together with a global bound for the Weyl conformal factor arising from a certain scalar curvature related differential inequality. Furthermore, as a corollary of the extremal version of this uniqueness result, we obtain the conjectured mass-angular momentum inequality for multiple dynamical black holes using the main theorem of [17].

I prove that a collapsing charged scalar field with scalar charge $\mathfrak{e}$ that forms an exactly extremal Reissner--Nordström black hole à la Kehle--Unger with radius $r_+$ must satisfy $|\mathfrak{e}|r_+ > \frac{1}{3}$. In contrast, I also show that no such bound exists for subextremal collapse: a scalar field with an arbitrary $\mathfrak{e} \ne 0$ can form a subextremal black hole with an arbitrary (subextremal) charge-to-mass ratio. Complementing the extremal bound, I prove that a Schwarzschild black hole can become extremal if $|\mathfrak{e}|r_+ > \sqrt{\frac{3 + \sqrt{33}}{3}}$. The fact that $|\mathfrak{e}|r_+$ can be taken to be of order unity could be considered evidence that the third law of black hole mechanics in vacuum is false.

In multi-fluid, matter-type gravity, the Bianchi identity constrains the divergence of the total effective stress tensor but does not, by itself, determine the currents assigned to its constituent sectors. Those currents are selected only after the off-shell matter action and equations, or an additional phenomenological closure, have been specified. We formulate this distinction and examine a spatially flat two-fluid background inspired by scale-independent EMSG. The case study combines an algebraic perfect-fluid prescription with a vanishing contracted-Hessian contribution, together with separate conservation of the total conventional and modification sectors. Because neither assumption is derived here from a concrete off-shell fluid action, the resulting system is an effective background closure rather than a microscopic EMSG model. For unequal EoS and every $α\ne0$, the closure yields a Barrow-Clifton system whose transfer coefficients depend on $(w_1,w_2)$ but not on $α$. Hence the nonzero-$α$ family is singular: its $α\to0$ limit does not recover the uncoupled GR conservation laws. We solve the two density eigenmodes and derive the associated modified Liénard equation for $H$. We also prove that the discriminant governing rank loss of the density-reconstruction map is strictly positive for every finite $w_1\ne w_2$. Thus each genuinely quadratic case has two distinct real rank-degenerate couplings. Exact vacuum and stiff-fluid families illustrate modal cancellation. When today's conventional densities are positive, the intervals on which both remain positive generally terminate at finite endpoints in the parameter ranges analyzed explicitly. These results provide a consistency diagnostic for separating action-level predictions from closure artifacts in multi-fluid, matter-type gravity; they do not establish an observationally viable EMSG model.

Light spectator fields provide a compelling inflationary origin for scalar dark matter, but their primordial fluctuations often generate excessive isocurvature perturbations on large scales. We develop a novel, exactly solvable filtering mechanism based on a tower of $δ-$function contributions to the spectator mass, which we call kicks, whose dense limit converges to a smooth Coulomb-like profile. The filter produces a isocurvature spectrum during inflation that is strongly suppressed on large scales, rises with momenta as $k^3$ in the infrared, and, after a rapid turnover, smoothly approaches a constant plateau at small scales. This construction is motivated by a discrete translational invariance, admits an EFT regime in which radiative corrections can be systematically controlled, and can find a concrete realization in terms of an inflationary massive-clock model. We embed the filter in a minimal dark-matter scenario in which a spectator condensate oscillates about the minimum of a quadratic potential and behaves as pressureless matter. In this framework, the spectator dark-matter mass is directly related to the turnover scale of the inflationary isocurvature spectrum. We determine the resulting dark-matter abundance and analyze the relevant isocurvature and backreaction constraints. The most restrictive realization, in which the spectator condensate is generated during inflation, favors a relatively low inflationary scale.