Locally authored papers of the past 5 days

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular-momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local $\alpha$ viscosity to zeroth order in aspect ratio ($H/R$), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal O(\alpha^{1/2})$ deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular-momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in $H/R$ and $\alpha^{1/2}$. The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different $m$. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with $m$-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

We revisit the potential of 21cm radiation fluctuations, during the epoch of reionization, in constraining the amplitude of local primordial non-Gaussianity (PNG) $f_{\rm NL}^{\rm loc}$. There generically exists an epoch at which the linear bias of the 21cm field crosses zero, independent of the precise astrophysics of reionization. This epoch implies the 21cm radiation is a natural "zero-bias tracer" in the sense of Castorina et al (2018). We identify new noise-like contributions which directly compete with the zero-bias effect, but which should be mitigated through sophisticated analysis techniques such as field-level reconstruction. These noise-like terms act to hinder the constraining power on local PNG of the brightness temperature fluctuations, making $\sigma (f_{\rm NL}^{\rm loc}) \leq 1$ unachievable even in simplified forecasts. We show that analyses which can reach the 'sampling noise' floor for this tracer and harness its full power can potentially unlock a 10-fold reduction in error bars, even in the presence of large-scale cuts from foregrounds. The potential of this epoch motivates searching for it in future 21cm surveys, along with developing analysis techniques that can reach the noise floor required for the zero-bias epoch to saturate Fisher information.

We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) and \textit{Planck}. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing bandpowers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of $A_\mathrm{lens}^\mathrm{recon} = 1.025 \pm 0.017$ with respect to the theory prediction from the best-fit CMB \textit{Planck}-ACT cosmology. The bandpowers from all three lensing datasets, analyzed jointly, yield a $1.6\%$ measurement of the parameter combination $S_8^\mathrm{CMBL} \equiv \sigma_8\,(\Omega_m/0.3)^{0.25} = 0.825^{+0.015}_{-0.013}$. Including Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to $\sigma_8 = 0.829 \pm 0.009$ (a $1.1\%$ determination). When combining with uncalibrated supernovae from \texttt{Pantheon+}, we present a $4\%$ sound-horizon-independent estimate of $H_0=66.4\pm2.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}} $. The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a $\Lambda$CDM model fit to the primary CMB data from \textit{Planck} and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the $\Lambda\mathrm{CDM}$ cosmological model, the combination of primary CMB, BAO and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

What happens when a black box (neural network) meets a black box (simulation of the Universe)? Recent work has shown that convolutional neural networks (CNNs) can infer cosmological parameters from the matter density field in the presence of complex baryonic processes. A key question that arises is, which parts of the cosmic web is the neural network obtaining information from? We shed light on the matter by identifying the Fourier scales, density scales, and morphological features of the cosmic web that CNNs pay most attention to. We find that CNNs extract cosmological information from both high and low density regions: overdense regions provide the most information per pixel, while underdense regions -- particularly deep voids and their surroundings -- contribute significantly due to their large spatial extent and coherent spatial features. Remarkably, we demonstrate that there is negligible degradation in cosmological constraining power after aggressive cutting in both maximum Fourier scale and density. Furthermore, we find similar results when considering both hydrodynamic and gravity-only simulations, implying that neural networks can marginalize over baryonic effects with minimal loss in cosmological constraining power. Our findings point to practical strategies for optimal and robust field-level cosmological inference in the presence of uncertainly modeled astrophysics.

The Cosmic Dawn Survey Pre-launch (PL) catalogues cover an effective 10.13 deg$^{2}$ area with uniform deep Spitzer/IRAC data ($m\sim25$ mag, 5$\sigma$), the largest area covered to these depths in the infrared. These data are used to gain new insight into the growth of stellar mass across cosmic history by characterising the evolution of the galaxy stellar mass function (GSMF) through $0.2 < z \leq 6.5$. The total volume (0.62 Gpc$^{3}$) represents a tenfold increase compared to previous works that have explored $z > 3$ and significantly reduces cosmic variance, yielding strong constraints on the abundance of massive galaxies. Results are generally consistent with the literature but now provide firm estimates of number density where only upper limits were previously available. Contrasting the GSMF with the dark matter halo mass function suggests that massive galaxies ($M \gtrsim10^{11}$ M$_{\odot}$) at $z > 3.5$ required integrated star-formation efficiencies of $M/(M_{\rm h}f_{\rm b}) \gtrsim$ 0.25--0.5, in excess of the commonly-held view of ``universal peak efficiency" from studies on the stellar-to-halo mass relation (SHMR). Such increased efficiencies imply an evolving peak in the SHMR at $z > 3.5$ which can be maintained if feedback mechanisms from active galactic nuclei and stellar processes are ineffective at early times. In addition, a significant fraction of the most massive quiescent galaxies are observed to be in place already by $z\sim 2.5$--3. The apparent lack in change of their number density by $z\sim 0.2$ is consistent with relatively little mass growth from mergers. Utilising the unique volume, evidence for an environmental dependence of the galaxy stellar mass function is found all the way through $z\sim 3.5$ for the first time, though a more careful characterisation of the density field is ultimately required for confirmation.

The dynamical method provides an efficient way to discover post-common-envelope binaries (PCEB) with faint white dwarfs (WDs), thanks to the development of time-domain survey projects. We perform a comprehensive analysis of the PCEB system TMTS J15530469+4457458 (J1553), discovered by the Tsinghua University-Ma Huateng Telescopes for Survey, to explore its physical origin and evolutionary fate. This system is characterized by double-peaked Balmer emission lines, and a cross-correlation function is applied to derive its radial velocity (RV) from a series of phase-resolved Keck spectra. Analyses with the cross-correlation function suggest that this system is a single-lined spectroscopic binary and only one star is optically visible. Further analysis through Doppler tomography indicates that J1553 is a detached binary without an accretion disk. Under such a configuration, the simultaneous light-curve and RV fitting reveal that this system contains an unseen WD with mass $M_{\rm A}=0.56\pm 0.09\, M_{\odot}$, and an M4 dwarf with mass $M_{\rm B}=0.37\pm 0.02\,M_{\odot}$ and radius $R_{\rm B}=0.403^{+0.014}_{-0.015}\,R_{\odot}$. The extra prominent Balmer emission lines seen in the spectra can trace the motion of the WD, which are likely formed near the WD surface as a result of wind accretion. According to the MESA simulation, J1553 could have evolved from a binary consisting of a 2.0-4.0 ${M}_{\odot}$ zero-age-main-sequence star and an M dwarf with an initial orbital period $P_i\approx 201-476$ d, and the system has undergone a common-envelope (CE) phase. After about $3.3\times10^6$ yr, J1553 should evolve into a cataclysmic variable, with a transient state as a supersoft X-ray source at the beginning. J1553 is an excellent system for studying wind accretion, CE ejection physics, and binary evolution theory.

Long gamma-ray bursts (LGRBs), including their subclasses of low-luminosity GRBs (LL-GRBs) and X-ray flashes (XRFs) characterized by low spectral peak energies, are known to be associated with broad-lined Type Ic supernovae (SNe Ic-BL), which result from the core collapse of massive stars that lose their outer hydrogen and helium envelopes. However, the soft and weak end of the GRB/XRF population remains largely unexplored, due to the limited sensitivity to soft X-ray emission. Here we report the discovery of a fast X-ray transient, EP250108a, detected by the Einstein Probe (EP) in the soft X-ray band at redshift $z = 0.176$, which was followed up by extensive multiband observations. EP250108a shares similar X-ray luminosity as XRF\,060218, the prototype of XRFs, but it extends GRBs/XRFs down to the unprecedentedly soft and weak regimes, with its $E_{\rm peak} \lesssim 1.8\,\mathrm{keV}$ and $E_{\rm iso} \lesssim 10^{49}\, \mathrm{erg}$, respectively. Meanwhile, EP250108a is found to be associated with SN\,2025kg, one of the most luminous and possibly magnetar-powered SNe Ic-BL detected so far. Modeling of the well-sampled optical light curves favors a mildly relativistic outflow as the origin of this event. This discovery demonstrates that EP, with its unique capability, is opening a new observational window into the diverse outcomes of death of massive stars.

We present optical, radio, and X-ray observations of EP250108a/SN 2025kg, a broad-line Type Ic supernova (SN Ic-BL) accompanying an Einstein Probe (EP) fast X-ray transient (FXT) at $z=0.176$. EP250108a/SN 2025kg possesses a double-peaked optical light curve and its spectrum transitions from a blue underlying continuum to a typical SN Ic-BL spectrum over time. We fit a radioactive decay model to the second peak of the optical light curve and find SN parameters that are consistent with the SNe Ic-BL population, while its X-ray and radio properties are consistent with those of low-luminosity GRB (LLGRB) 060218/SN 2006aj. We explore three scenarios to understand the system's multi-wavelength emission -- (a) SN ejecta interacting with an extended circumstellar medium (CSM), (b) the shocked cocoon of a collapsar-driven jet choked in its stellar envelope, and (c) the shocked cocoon of a collapsar-driven jet choked in an extended CSM. All three models can explain the optical light curve and are also consistent with the radio and X-ray observations. We favor models (a) and (c) because they self-consistently explain both the X-ray prompt emission and first optical peak, but we do not rule out model (b). From the properties of the first peak in models (a) and (c), we find evidence that EP250108a/SN 2025kg interacts with an extended CSM, and infer an envelope mass $M_{\rm e} \sim 0.1\,\rm M_\odot$ and radius $R_{\rm e} \sim 4 \times 10^{13}$ cm. EP250108a/SN 2025kg's multi-wavelength properties make it a close analog to LLGRB 060218/SN 2006aj, and highlight the power of early follow-up observations in mapping the environments of massive stars prior to core collapse.

Open clusters offer unique opportunities to study stellar dynamics and evolution under the influence of their internal gravity, the Milky Way's gravitational field, and the interactions with encounters. Using the Gaia DR3 data for a catalog of open clusters within 500 parsecs that exhibit tidal features reported by the literature, we apply a novel method based on 3D principal component analysis to select a ``golden sample'' of nearby open clusters with minimal line-of-sight distortions. This approach ensures a systematic comparison of 3D and 2D structural parameters for tidally perturbed clusters. The selected golden sample includes Blanco 1, Melotte 20, Melotte 22, NGC 2632, NGC 7092, NGC 1662, Roslund 6 and Melotte 111. We analyze these clusters by fitting both 2D and 3D King Profiles to their stellar density distributions. Our results reveal systematic discrepancies: most of the golden sample clusters exhibit larger 3D tidal radii compared to their 2D counterparts, demonstrating that the 2D projection effects bias the measured cluster size. Furthermore, the 3D density profiles show stronger deviations from King profiles at the tidal radii ($\Delta \rho_{\rm 3D} > \Delta \rho_{\rm 2D}$), highlighting enhanced sensitivity to tidal disturbances. Additionally, we investigate the spatial distribution of cluster members relative to their bulk motion in the Galactic plane. We find that some clusters exhibit tidal features oriented perpendicular to their direction of motion, which can be attributed to the fact that the current surveys only detect the curved inner regions of the tidal features. In conclusion, this work offers a golden sample of nearby open clusters that are most reliable for 3D structure analysis and underscores the necessity of 3D analysis in characterizing OC morphological asymmetries, determining cluster size, and identifying tidal features.

The black hole at the center of M87 is observed to flare regularly in the very high energy (VHE) band, with photon energies $\gtrsim 100$ GeV. The rapid variability, which can be as short as $2$ days in the VHE lightcurve constrains some of the flares to originate close to the black hole. Magnetic reconnection is a promising candidate for explaining the flares, where the VHE emission comes from background soft photons that Inverse Compton (IC) scatter off of high energy electron-positron pairs in the reconnecting current sheet. In this work, we ray trace photons from a current sheet near the black hole event horizon during a flux eruption in a magnetically arrested state in a general relativistic magnetohydrodynamics simulation. We incorporate beaming of the Compton up-scattered photons, based on results from radiative kinetic simulations of relativistic reconnection. We then construct VHE lightcurves that account for the dynamics of the current sheet and lensing from general-relativistic effects. We find that most of the flux originates in the inner $5$ gravitational radii, and beaming is essential to explain the observed flux from the strongest VHE flares. The ray traced lightcurves show features resulting from the changing volume of the reconnecting current sheet on timescales that can be consistent with observations. Furthermore, we find that the amount of beaming depends strongly on two effects: the current sheet inclination with respect to the observer and the anisotropy in the direction of motion of the accelerated particles.

Radio astronomy is part of radio science that developed rapidly in recent decades. In the research of radio astronomy, pulsars have always been an enduring popular research target. To find and observe more pulsars, large radio telescopes have been built all over the world. In this paper, we present our studies on pulsars in M15 and NGC 6517 with FAST, including monitoring pulsars in M15 and new pulsar discoveries in NGC 6517. All the previously known pulsars in M15 were detected without no new discoveries. Among them, M15C was still detectable by FAST, while it is assumed to fade out due to precession [1]. In NGC 6517, new pulsars were continues to be discovered and all of them are tend to be isolated pulsars. Currently, the number of pulsars in NGC 6517 is 17, much more than the predicted before [2].

We present a ``cyclic zoom'' method to capture the dynamics of accretion flows onto black holes across a vast range of spatial and temporal scales in general relativistic magnetohydrodynamic (GRMHD) simulations. In this method, we cyclically zoom out (derefine) and zoom in (refine) the simulation domain while using a central mask region containing a careful treatment of the coarsened fluid variables to preserve the small-scale physics, particularly the magnetic field dynamics. The method can accelerate GRMHD simulations by $\gtrsim 10^5$ times for problems with large scale separation. We demonstrate the validity of the technique using a series of tests, including spherically symmetric Bondi accretion, the Blandford-Znajek monopole, magnetized turbulent Bondi accretion, accretion of a magnetized rotating torus, and the long-term evolution of an accreting torus about both Schwarzschild and Kerr black holes. As applications, we simulate Bondi and rotating torus accretion onto black holes from galactic scales, covering an extremely large dynamic range. In Bondi accretion, the accretion rate is suppressed relative to the Bondi rate by $\sim(10r_\mathrm{g}/r_\mathrm{B})^{1/2}$ with a feedback power of $\sim 0.01 \dot{M} c^2$ for vanishing spin, and $\sim 0.1 \dot{M} c^2$ for spin $a\approx0.9$. In the long-term evolution of a rotating torus, the accretion rate decreases with time as $\dot{M}\propto t^{-2}$ on timescales much longer than the viscous timescale, demonstrating that our method can capture not only quasi-steady problems but also secular evolution. Our new method likewise holds significant promise for applications to many other problems that need to cover vast spatial and temporal scales.

Using 15 years of data from the Fermi Large Area Telescope (Fermi-LAT), we performed a comprehensive analysis on the gamma-ray binary HESS J0632+057. Its spectrum in 0.1-300 GeV band is well described by a power law model with an index of $2.40\pm0.16$, leading to an energy flux of (5.5$\pm$1.6$)\times$ 10$^{-12}$ erg cm$^{-2}$ s$^{-1}$. The GeV Spectral Energy Distribution (SED) of HESS J0632+057 hints for a spectral turn-over between $\sim$10-100 GeV. Orbital analysis reveals a flux enhancement during the phase range of 0.2-0.4, consistent with the X-ray and TeV light curves, indicating an origin of a common particle population. We carried out six deep radio observations on HESS J0632+057 with the Five-hundred-meter Aperture Spherical Telescope (FAST), evenly distributed across its orbit, reaching a detection sensitivity of 2$\mu$Jy. However, no radio pulsation was detected within these observations. The absence of radio pulsation may be attributed to the dense stellar wind environment of HESS J0632+057.

We identified known Trans-Neptunian Objects (TNOs) and Centaurs in the complete Dark Energy Survey (DES) year six catalog (DES Y6) through the Sky Body Tracker (SkyBoT) tool. We classified our dataset of 144 objects into a widely used 4-class taxonomic system of TNOs. No such previous classification was available in the literature for most of these objects. From absolute magnitudes and average albedos, an estimation of the diameters of all these objects is obtained. Correlations involving colours, orbital parameters, dynamical classes and sizes are also discussed. In particular, our largest reddest object has a diameter of $390^{+68}_{-53}$ km and our largest cold classical, $255^{+19}_{-17}$ km. Also, a weak correlation between colour and inclination is found within the population of resonant TNOs in addition to weak correlations between colour and phase slope in different bands.

International and U.S. strategies and protocols have identified the need to develop rapid-response spacecraft reconnaissance capabilities as a priority to advance planetary defense readiness. A space-based reconnaissance response is recommended for potential impactors as small as 50 m, making these small objects the most likely to trigger a space-based response and the ones that drive the reconnaissance capabilities needed. Even following the successful completion of the NEO Surveyor mission and Rubin Observatory survey efforts, roughly half of the 50-m near-Earth object (NEO) population will remain undiscovered. As a result, 50-m impactors may not be found with long warning times, and a rapid-response flyby mission may be the only reconnaissance possible. To develop a robust flyby reconnaissance capability for planetary defense, four major requirements are defined for a demonstration mission. 1. Enable a flyby of greater than 90 percent of the potential asteroid threat population. 2. Demonstrate the flyby reconnaissance for a 50 m NEO. 3. Obtain the information needed to determine if and where it would impact the Earth. 4. Determine key properties of the asteroid to inform decision makers. As commonly noted in the planetary defense community, in planetary defense, you do not pick the asteroid, the asteroid picks you. Thus, a planetary defense flyby reconnaissance demonstration mission is not about just flying by an asteroid, but rather it is about developing a robust capability for the objects that are most likely to require a short-warning-time, space-based response.

Using $10,\!080^3$ grid simulations, we analyze scale-dependent alignment in driven, compressible, no net-flux magnetohydrodynamic turbulence. The plasma self-organizes into localized, strongly aligned regions. Alignment spans all primitive variables and their curls. Contrary to incompressible theory, velocity-magnetic alignment scales as $\theta(\lambda) \sim \lambda^{1/8}$, where $\lambda$ is the scale, suggesting a distinct three-dimensional eddy anisotropy and a much higher critical transition scale toward a reconnection-mediated cascade.

High-energy gamma rays have been detected in the region of LHAASO~J2108+5157 by the Fermi--LAT, HAWC and LHAASO-KM2A observatories. Cygnus~OB2 in Cygnus--X has been confirmed as the first strong stellar cluster PeVatron in our Galaxy. Thus, the star--forming regions Kronberger~80 and Kronberger~82, located in the field of LHAASO~J2108+5157, are analyzed to evaluate their stellar population and potential as associated PeVatron candidates. A distance of 10~kpc is adopted for Kronberger~80, while $\sim$1.6~kpc is estimated for Kronberger~82. Based on stellar densities, we report that their cluster radii are 2.5$\arcmin$ and 2.0$\arcmin$, while IR photometry reveals poor stellar content in massive O-type stars in both cases. From optical data, the estimation of cluster ages are 5--12.6~Myr and $\lesssim$ 5~Myr, respectively. We conclude that, in contrast to the stellar content of Cygnus~OB2, it is unlikely that Kronberger~80 and Kronberger~82 are PeVatrons associated with LHAASO~J2108+5157. The presence of a PeVatron in this region remains a mystery, but we confirm that the two Kronberger regions are star--forming regions undergoing formation rather than evolution.

The evolution of large-scale structure, galaxies and the intergalactic medium (IGM) during the Epoch of Reionization (EoR) can be probed by upcoming Line Intensity Mapping (LIM) experiments, which sample in redshift and direction without needing to resolve individual galaxies. We predict the intensity and sources of hydrogen H$\alpha$ emission, dominated by radiative recombination following ionization by UV from the same massive stars that caused reionization, down to redshift 4.6, using the largest fully-coupled, radiation-hydro simulation of galaxy formation and reionization to date, Cosmic Dawn (CoDa) III. We compute the mean intensity and Voxel Intensity Distribution (VID) vs. redshift, including the relative contributions of galaxies and IGM. This will provide mock data to guide and interpret LIM experiments such as NASA's SPHEREx and proposed Cosmic Dawn Intensity Mapper (CDIM).

Eclipsing binary systems (EBs), as foundational objects in stellar astrophysics, have garnered significant attention in recent years. These systems exhibit periodic decreases in light intensity when one star obscures the other from the observer's perspective, producing characteristic light curves (LCs). With the advent of the Transiting Exoplanet Survey Satellite (TESS), a vast repository of stellar LCs has become available, offering unprecedented opportunities for discovering new EBs. To efficiently identify such systems, we propose a novel method that combines LC data and generalized Lomb-Scargle periodograms (GLS) data to classify EBs. At the core of this method is CNN Attention LSTM Net (CALNet), a hybrid deep learning model integrating Convolutional Neural Networks (CNNs), Long Short-Term Memory (LSTM) networks, and an Attention Mechanism based on the Convolutional Block Attention Module (CBAM). We collected 4,225 EB samples, utilizing their 2-minute cadence LCs for model training and validation. CALNet achieved a recall rate of 99.1%, demonstrating its robustness and effectiveness. Applying it to TESS 2-minute LCs from Sectors 1 to 74, we identified 9,351 new EBs after manual visual inspection, significantly expanding the known sample size. This work highlights the potential of advanced deep-learning techniques in large-scale astronomical surveys and provides a valuable resource for further studies on EBs.