Locally authored papers of the past 5 days

Today, the observable cosmos exhibits a remarkable degree of isotropy and plausibly began in a nearly isotropic initial state. The properties of the Lorentzian Chern-Simons-Kodama (CSK) functional can provide an understanding of this initial state. In gravity with a positive cosmological constant, the Chern-Simons-Kodama (CSK) wavefunctional is an exact, chiral solution of the quantum gravitational constraints. We suggest that the normalizability and other issues with this functional, if interpreted as a proper state of quantum gravity, instead suggest an embedding into a larger quantum gravitational completion, and recast the CSK functional as a gravitational sphaleron with observationally desirable properties. By perturbing around the dominant de Sitter saddle of the wavefunctional with appropriate quantum gravitational boundary conditions, we find that for a closed universe the system is dynamically driven to spatial isotropy, while all anisotropic modes acquire positive quadratic curvature and are Gaussian-suppressed. The decay of this sphaleron therefore proceeds along an isotropic channel, providing an intrinsic quantum-gravitational mechanism for dynamical isotropization. This isotropization effect is robust under the inclusion of a slow-roll inflaton, and no analogous isotropic sphaleron exists for spatially flat or hyperbolic geometries. Taken together, these results recast the Lorentzian CSK functional as a chiral sphaleron that naturally prepares an approximately isotropic de Sitter background for inflation. Beyond this phenomenological study, we further suggest that the CSK functional can be understood as a boundary functional for a class of anomaly-free objects, including a complexified generalization of the Hartle-Hawking state.

Fast surrogate models for expensive simulations are now essential across the sciences, yet they typically operate as black boxes. We present \texttt{GWAgent}, a large language model (LLM)-based workflow that constructs interpretable analytic surrogates directly from simulation data. Surrogate modeling is well suited to agentic workflows because candidate models can be quantitatively validated against ground-truth simulations at each iteration. As a demonstration, we build a surrogate for gravitational waveforms from eccentric binary black hole mergers. We show that providing the agent with a physics-informed domain ansatz substantially improves output model accuracy. The resulting analytic surrogate attains a median Advanced LIGO mismatch of $6.9\times10^{-4}$ together with an $\sim 8.4\times$ speedup in waveform evaluation, surpassing both symbolic regression and conventional machine learning baselines. Beyond producing an accurate model, the workflow identifies compact physical structure from the learned representation. As an astrophysical application, we use \texttt{GWAgent} to analyze the eccentricity of GW200129 and infer $e_{20\mathrm{Hz}}=0.099^{+0.063}_{-0.044}$. These results show that validation-constrained agentic workflows can produce accurate, fast, and interpretable surrogates for scientific simulations and inference.

Tensions often arise between different datasets in cosmology, and consistency tests can serve as a powerful tool for diagnosing potential issues. The density-shear Baryon Acoustic Oscillations (GI BAO) are the imprint of the BAO feature on the shear field induced by the large-scale tidal field. We highlight that GI BAO can provide a robust consistency check for the density BAO, shear measurement, and alignment model. Failure of this check hints at systematics in any of these parts. As an illustration, we present the first GI BAO measurement on photometric data, using the DES Y3 dataset. We find the GI BAO constraint on the BAO scale dilation parameter $\alpha $ to be $ 0.966 \pm 0.252 $ (1$\sigma$), in good agreement with the density BAO constraint, $ 0.966 \pm 0.037 $, thereby validating the density BAO, shear measurement, and the linear alignment model. Furthermore, we argue that combining the density BAO with the GI BAO yields results that are more resilient to systematic effects. Thanks to the massive data volumes of stage IV surveys, the GI BAO will play an even more prominent role as a consistency check.

Light primordial black holes heat the surrounding plasma via Hawking radiation, forming localized hotspots whose temperature may far exceed that of the cosmological background. Previous studies of hotspot formation and cooling have treated the subsequent energy transport in flat spacetime, thereby neglecting the expansion of the Universe. We formulate the diffusion equation governing the hotspot evolution, in an expanding universe, and clarify the regime in which the formalism is valid. We find that hotspot formation is robust against cosmological expansion. We show that the critical distance scale, where Hubble expansion overtakes diffusion, coincides with the decoupling radius introduced in earlier work, and the temperature profile $T\propto r^{-7/11}$ essentially remains unchanged. However, the cooling stage is substantially modified. We find that the plateau temperature of a cooling hotspot initially undergoes a rapid drop and then follows $T_{\rm plt} \propto t^{-11/15}$, steeper than the flat-spacetime scaling $t^{-7/15}$. This scaling cannot be obtained by simply redshifting the flat-spacetime solution, because expansion also suppresses diffusive transport. As a consequence, all hotspots disappear within a finite time, as opposed to the flat-spacetime prediction of everlasting hotspots in part of the parameter space.

We present the first public data release of DDO51 band from the Stellar Abundances and Galactic Evolution Survey (SAGES), based on Nanshan One-meter Wide-field Telescope (NOWT) observations obtained between 2023 September and 2024 January. This release initiates the DDO51-band component of the survey, covering $\sim$ 2,500 deg$^2$ of the northern sky and including more than 10 million sources. The DDO51 filter is centered near the \ion{Mg}{1}~$b$ triplet and the adjacent MgH feature, offering sensitivity to stellar surface gravity. The data reduction pipeline incorporates an improved astrometric solution anchored to Gaia DR3 and a photometric calibration strategy tied to synthetic photometry from Gaia XP spectra. These procedures yield a point-source depth of $\sim$18.9 mag at S/N$\sim$10 and an internal photometric precision $\approx$6-7 mmag at the bright end. A preliminary color--color analysis using Gaia broadband photometry confirms the expected sensitivity of the DDO51 band to stellar surface gravity, demonstrating a clear photometric separation between dwarf and giant sequences for late-type stars. This dataset, when combined with existing SAGES photometry in other bands, provides a crucial tool for disentangling the substructures of the Milky Way. All data products from this release upon publication will be available.

Linear C4H and cyclic c-C3H2, as small unsaturated hydrocarbons, are the key precursors to complex organic molecules and are critical components of the interstellar medium. We present on-the-fly mapping observations of C4H 9-8 lines, c-C3H2 2-1, H13CO+ 1-0, and H42 toward a sample of 22 massive star-forming regions using the IRAM 30m telescope. Our aim is to further explore the evolution of these carbon-chain molecules by combining observational results obtained in cold cores. We employed H13CO+ 1-0 and H42 as tracers to probe the positions of molecular cloud cores and ionised hydrogen regions (HII regions), respectively. One chemical model in particular, which includes gas, dust grain surface, and icy mantle phases for C4H and c-C3H2 molecules, was used to make comparisons with observed abundances. From mapping observations targeting 31 regions across 22 sources, C4H 9-8 (J = 19/2-17/2) and C4H 9-8 (J = 17/2-15/2) were detected in only 17 regions, while H13CO+ 1-0 and c-C3H2 2-1 were successfully detected in all 31 regions. We find that the emission of C4H 9-8 and c-C3H2 2-1 is concentrated at the edges of H42 emission regions. The C4H/H13CO+ and c-C3H2/H13CO+ relative abundance ratios range from 0.17 to 1.77 and 1.42 to 6.69, respectively, with a median C4H/c-C3H2 ratio of 0.13. By combining the observational results of cold cores, we find that C4H/H13CO+ and c-C3H2/H13CO+ ratios show a strong decreasing trend as molecular cores evolve. The decreasing trends in C4H/H13CO+ and c-C3H2/H13CO+ ratios imply that small unsaturated hydrocarbons can be consumed and converted into other organic molecules during the evolution of molecular cores. The spatial concentration of C4H and c-C3H2 emission at the edges of H42 regions further supports their role as precursors in the chemical pathways that lead to complex organic molecules in the interstellar medium.

The identification of physically associated kiloparsec-scale quasar pairs is important for understanding galaxy evolution, the growth of supermassive black holes, and their co-evolution with host galaxies. However, their rarity and the high contamination from stellar superpositions and projected alignments require efficient pre-selection methods. We develop a machine-learning framework to produce photometric-redshift point estimates and redshift probability density functions for quasars, with the main goal of identifying high-probability quasar pair candidates in the MGQPC catalogue. We construct two large spectroscopically confirmed quasar samples with multi-wavelength photometry, based on SDSS and DESI Legacy Imaging Surveys data. CatBoost is used for point-estimate photometric-redshift regression, and FlexZBoost is used for full redshift-PDF estimation. The workflow achieves robust performance, with a normalised median absolute deviation of 0.036 and an outlier fraction of 5.6% on the test sample. Applying the trained model to the MGQPC catalogue, we identify 185 high-probability quasar pair candidates based on photometric-redshift consistency. Among them, 20 systems have been subsequently confirmed as genuine physical pairs by independent spectroscopic observations. The resulting MGQPC photometric-redshift catalogue provides a useful resource for future spectroscopic follow-up of quasar pairs and dual supermassive black holes.

Nonlinear Force-free Field (NLFFF) models are widely used to investigate coronal magnetic field structure in solar active regions, but methods to validate them remain limited. Here, we use Gaussian separation, recently applied to solar vector magnetogram data, to assess the accuracy of NLFFF models constructed with two methods: optimization and the current-field iteration (CFIT) implementation of the Grad-Rubin method. Gaussian separation partitions the photospheric vector magnetic field into three components associated with currents flowing below, above, and passing through the photosphere, respectively. Comparing the photospheric field components due to coronal currents in an NLFFF model with those in the original vector magnetogram data provides a check on the accuracy of the model's coronal currents. We consider NLFFF models constructed for the active region AR 11429. The photospheric signatures of coronal currents in both the models and the vector magnetogram data indicate currents flowing above and parallel to central, sheared polarity inversion lines (PILs), consistent with other recent studies. We find that while both models reproduce the coronal current signatures along the upper section of the main PIL, the CFIT model significantly alters the signature of a flux rope along the lower section of the PIL, including shifting its positive-polarity footpoint. These differences arise from modifications to the vector magnetogram boundary data when solving the NLFFF equations, and from the assumptions underlying the models. We propose Gaussian separation as a useful tool to validate coronal magnetic field models, in addition to existing methods.

Intensity ratios of aromatic emission features are widely used to diagnose the size and ionization state of polycyclic aromatic hydrocarbons (PAHs) in astronomical environments. However, PAHs are known to typically carry aliphatic side chains, a structural feature that may compromise the reliability of traditional diagnostic methods. This study systematically investigates the effects of aliphatic components on the aromatic emission properties of PAHs. Based on theoretical data from the NASA Ames PAH IR Spectroscopic Database, we compare the emission behavior of purely aromatic PAHs with those containing aliphatic substituents, revealing that aliphatic functionalization may modify the intensity ratio of the 11.2 $\mu$m band relative to the 7.7 $\mu$m and 3.3 $\mu$m bands. This potentially leads to misidentification of their ionization state if molecular structural effects are neglected. Further analysis indicates that the impact of aliphatic components on diagnostic band ratios strongly depends on PAH size: small PAHs exhibit significant emission ratio shifts, deviating from traditional size/ionization trends, while larger PAHs are minimally affected. Despite these shifts, the classic $(I_{11.2/7.7})$ versus $(I_{11.2/3.3})$ diagnostic grid remains largely applicable to mixed aromatic-aliphatic PAHs, although some systematic calibration may be needed. Our findings emphasize the necessity for caution when interpreting PAH band ratios in aliphatic-rich environments, as variations in PAH molecular composition may distort inferences about physical conditions.

Recent advances in astronomical observations have ushered in an era of remarkable discoveries. We now probe the Universe through multi-messenger signals, image the sky with unprecedented depth and resolution, and investigate individual sources using powerful large-aperture telescopes. Yet, a critical gap persists: the lack of wide-field, highly multiplexed spectroscopic capabilities needed to fully exploit the wealth of imaging data from current and upcoming surveys. In this review, we trace the historical development of large optical telescopes and spectroscopic surveys, assess the capabilities of ongoing and near-future facilities, and motivate the need for next-generation Stage-V spectroscopic experiments. As a representative example, we present the MUltiplexed Survey Telescope (MUST), the first Stage-V spectroscopic facility currently under construction. MUST is a 6.5-meter telescope designed to obtain optical spectra for over 20,000 targets simultaneously within a $\sim$5 deg$^2$ field, using a modular focal plane populated with 6.2-mm pitch fiber-positioning robots. Over an 8-year survey in the 2030s, MUST aims to build the most comprehensive 3D spectroscopic map of the Universe to date, measuring redshifts for over 100 million galaxies and quasars and opening new windows into cosmology, Galactic structure, and time-domain astrophysics.

This paper presents a high-precision gravitational redshift test using the China Space Station (CSS) Laser Time Transfer (CLT) system. We develop a comprehensive observation equation based on a c^{-3} order relativistic model for space-ground clock comparison. While the CSS optical clock system is currently in the orbital debugging phase, our simulation using actual CSS orbit data achieves a gravitational redshift verification precision of (1.8 \pm 47)*10^{-7} -- approximately one order of magnitude improvement over previous experiments. Our work represents the first application of laser-based time transfer for gravitational redshift verification at such precision, and the first use of the CSS CLT link for testing this fundamental aspect of General Relativity. Unlike microwave-based methods, our laser approach avoids ionospheric effects and first-order Doppler shifts. Residual analysis identifies tropospheric delay variations and atmospheric turbulence as the primary remaining uncertainty contributors. The achieved precision enables gravitational potential difference measurements with 0.1 m^2/s^2 precision -- offering new capabilities for both fundamental physics investigations and geodetic applications including intercontinental height transfer. This work establishes a new benchmark for high-precision tests of relativistic physics and demonstrates the transformative potential of space-based optical time transfer.

Galaxy formation scenarios can be interpreted through galaxy morphology and the level of rotational versus pressure support, quantified through the ratio of a galaxy's rotation speed to its velocity dispersion: $V/\sigma$. Observational studies of dwarf galaxies find that $V/\sigma$ does not strongly depend on environment, and may weakly depend on galaxy mass, which could shift our understanding of how dwarf galaxies form. We utilize the Marvelous Massive Dwarfs suite to examine whether $V/\sigma$ depends on mass in simulations, and understand how this varies for different baryonic components of the galaxy: HI gas, young stars ($<$ 1 Gyr) and old stars ($>$ 1 Gyr). We use a simulation sample of 67 isolated dwarf galaxies with M$_\star=10^6-10^9$ M$_\odot$ and produce line-of-sight maps for rotation speed and dispersion for different viewing angles of each galaxy. We find that $V/\sigma$ increases with mass, and that HI gas and young stars are more rotation-supported ($V/\sigma\approx 1-13$) while old stars are more dispersion-supported ($V/\sigma\approx 0.2-5$). This result is consistent with the scenario where young stars are born from dynamically cold gas in the interstellar medium and undergo dynamical heating over time. We quantify the effects of spatial resolution in observational determinations of $V/\sigma$ and find that existing observations using old stars may underestimate the intrinsic $V/\sigma$. We find a correlation between $V/\sigma_\mathrm{HI,global}$ and HI line profile shape that is qualitatively similar to previous simulation results, but we find higher $V/\sigma_\mathrm{HI,global}$ compared to prior work which found values $\lesssim 2$ for most galaxies in this mass range. Our results motivate future work to examine $V/\sigma$ and dwarf galaxy formation with different kinematic tracers of the galaxy.

Stellar streams trace the gravitational potential of their host galaxies and offer a direct probe of dark matter halo geometry. Cosmological simulations predict that halo shapes depend on both baryonic physics and the nature of dark matter, yet observational constraints on halo flattening and orientation remain limited, especially for individual galaxies. We present Potamides, which utilizes the curvature of extragalactic stellar streams to derive constraints on halo shapes. We apply Potamides to 15 stellar streams from the Stellar Stream Legacy Survey to infer the projected axis ratios and orientation of their host halos. We find that some streams in our sample exclude large regions of halo flattenings and halo orientations. Systems with edge-on wrapping loops or sharp turning points yield the strongest constraints, whereas great circle-like streams remain largely uninformative. All streams in our sample support a spherical halo for a given flattening direction. These results demonstrate that stream morphology can provide halo shape constraints for individual external galaxies. With upcoming surveys (such as Euclid, Rubin, Roman, and ARRAKIHS) expected to discover large numbers of stellar streams, this curvature-based technique will enable rapid statistical tests of dark matter and baryonic physics through the shapes and alignments of halos and disks across cosmic time.

Large spectroscopic and astrometric surveys have revealed complex wave-like features in the Milky Way disk, suggesting that its kinematic and chemical structures are shaped by time-dependent perturbations. Recent studies have reported oscillatory patterns in the Rg-Vphi-VR space, hinting at a possible structural transition in the outer disk. We aim to characterise the transition between the inner and outer Galactic thin disk and to investigate whether radial corrugations can provide a plausible physical interpretation of the observed features. We analysed two large stellar samples from LAMOST DR8 and Gaia DR3, combining spatial, kinematic, and chemical diagnostics. A simplified corrugation model consisting of two radial waves propagating in opposite directions was constructed and fitted to the observed VR pattern. We further validated the model using N-body simulations. Both LAMOST and Gaia samples reproduce the previously reported wave-like pattern in the Rg-Vphi-VR plane. We identify a clear transition between the inner and outer disks via the variations in rotational velocity and metallicities. The corrugation model naturally reproduces the periodic variation of VR with galactocentric radius, and the superposition of the inward and outward propagating modes gives rise to a comparable oscillatory pattern in both observations and simulations. Our modelling suggests that radial corrugations can provide a plausible interpretation of the observed kinematic signatures. The results highlight the complex, multi-perturber nature of the Galactic disk and motivate further investigation with upcoming surveys.

The 3D mass distributions of galaxy clusters are generally triaxial, a geometry that is difficult to constrain from projected observations. In this work, we measure the projected halo shapes of clusters from their weak lensing signatures using the triaxiality functionality in the Cluster Lensing Mass Modeling software, a tool developed by the Dark Energy Science Collaboration to analyze data from NSF-DOE Rubin Observatory's Legacy Survey of Space and Time (LSST). We measure ensemble halo ellipticity on the plane of the sky via axis-aligned stacking and multipole expansion of the weak lensing data. We study a precursor dataset -- the redMaPPer cluster catalog, the metacalibration shape catalog, and the Directional Neighborhood Fitting photometric redshift catalog from the Dark Energy Survey Year 3 public data release. We select clusters that have a high centering probability (>90%) of the identified central galaxy, and use the satellite galaxy distribution to determine the major-axis orientation for stacking. We extend the analysis to the second order of ellipticity in the monopole and quadrupole measurement. The projected ellipticity of the cluster sample is found to be $0.310^{+0.017}_{-0.016}$ (axis ratio $0.527^{+0.018}_{-0.019}$). The projected cluster ellipticity shows no statistically significant dependence on mass and redshift. We further verify the accuracy of the cluster shape measurement using mock catalogs. This analysis is applicable to datasets from upcoming wide-area cosmic surveys such as LSST, Euclid, and the Roman Space Telescope, where larger sample sizes will lead to tighter constraints on the cluster ellipticities.

Cosmic magnetic fields are typically inhomogeneous and often highly tangled due to large-scale plasma flows, turbulence, and instabilities. If the variations in the magnetic field occur on scales that are large compared to the gyro-radius of the plasma electrons, the electrons are primarily confined to gyro-centre trajectories along the field lines. Therefore, in-situ electron measurements help us map out the connectivity of the magnetic field in space plasmas. Gyro-centre drifts, wave-particle interactions, trapping, and cross-field diffusion are processes related to field inhomogeneities and fluctuations; they have the potential to modify or even disrupt the transport of electrons along field lines. We introduce the basic principles of electron transport in tangled magnetic fields and review the creation of tangled fields through turbulence and instabilities as well as the modulation of parallel electron transport through kinetic instabilities. We then describe trapping and de-trapping effects in inhomogeneous magnetic fields, as well as electron diffusion and energisation across the magnetic field. The transport of electrons in tangled fields results from a complex interplay of plasma processes that occur on a broad range of scales. A combination of in-situ plasma measurements, remote-sensing plasma observations, and plasma theory and simulations is required to resolve this contemporary challenge to the fields of heliophysics and astrophysics.

Baryonic feedback from active galactic nuclei (AGN) is often invoked as a major source of suppression in the matter power spectrum, with implications for precision cosmology and the $S_8$ tension. We present Astrid-DMO, the dark matter-only counterpart to the large-volume Astrid hydrodynamical simulation, and measure baryonic effects through $P_{\rm hydro}(k)/P_{\rm DMO}(k)$. We find no significant suppression at $z=0$ and mild suppression at $z=0.2$, weaker than in other state-of-the-art simulations. Using controlled small-volume runs, we identify a key driver of this discrepancy: the treatment of black hole (BH) dynamics. The widely used BH repositioning scheme artificially enhances BH mergers and boosts kinetic AGN feedback (e.g., by a factor of $2$ at $z=1.5$), leading to overly strong suppression. By contrast, a more physical dynamical friction model reduces feedback efficiency and weakens clustering suppression. Consequently, reconciling large-scale structure measurements with cosmic microwave background (CMB)-inferred $\Lambda$CDM cosmology with AGN feedback becomes more challenging. Although strengthening AGN feedback can increase suppression, in our model this induces tensions with the observed galaxy stellar mass and AGN luminosity functions. These results motivate considering either new non-baryonic physics that suppresses late-time matter clustering, or novel mechanisms that can efficiently eject gas from halos without compromising other galaxy properties.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino observatory designed to detect the impulsive radio waves produced by relativistic particle cascades in the Antarctic glacial ice. Using a significantly enhanced simulation pipeline, which adds data-driven detector simulations and fully incorporates secondary particle production, we calculate the trigger-level acceptance of the entire array. We compare the resulting trigger-level sensitivity to constraints on the UHE neutrino flux from other detectors. Given its exposure from 2013 to 2023, we find that ARA achieves a world-leading sensitivity above about $10^{19}$ eV, depending on the details of the event selection used in a search. Moreover, we find that up to 13 neutrinos are predicted to have been observed in this period at trigger-level, assuming the most optimistic neutrino flux models. We show that observations of secondary particles account for up to 30\% of the total acceptance starting at $10^{19}$ eV, and we explore the potential signatures and implications of both multi-pulse (from direct and refracted pulses and/or from secondary particle interactions) and multi-station events. Finally, we comment on the implications of this study for the design of next-generation UHE neutrino experiments, in particular IceCube-Gen2 Radio.

We run and analyze a suite of high-redshift zoom-in cosmological simulations with varying supernova feedback and supermassive black hole (SMBH) accretion prescriptions to study the joint evolution of stellar and SMBH mass in high-redshift galaxies down to $z=10$. The simulations reproduce the observed high-$z$ $M_{\mathrm{BH}}/M_{\star}$ relation if super-Eddington accretion is allowed prior to the final self-regulated phase. To extend the evolution to lower redshift, we model subsequent black hole and host growth using analytic halo assembly histories combined with a redshift-dependent effective Eddington duty cycle, $f_{\rm duty}=0.0004(1+z)^3$, calibrated to observations at $z\le6$, with conservative uncertainties at higher redshift. Within this framework, $M_{\mathrm{BH}}/M_{\star}$ exhibits a broad peak at $z\sim7$--10, reaching a few percent up to $\sim30\%$, followed by a steady, approximately power-law decline toward $z=0$. The model predicts $M_{\mathrm{BH}}/M_{\star}\sim(0.002,0.003,0.006,0.016,0.071,0.156)$ at $z=(0,1,2,3,5,10)$, consistent with available observations. This evolution is driven by rapid SMBH growth at high redshift, with effective mass e-folding times shorter than those of stellar mass, while at later times galaxy growth dominates, leading to the decline in $M_{\mathrm{BH}}/M_{\star}$. These results demonstrate that the emergence of a high-redshift peak and subsequent decline is robust despite uncertainties in the duty-cycle normalization.