Locally authored papers of the past 5 days

The orbits of small bodies in the outer solar system are particularly sensitive to gravitational perturbations, including stellar flybys. Stellar clusters, with low velocity dispersions and high number densities, can be the source of strong and frequent flybys. As a result, we can infer what properties of the solar birth environment would be incompatible with the structure of the outer solar system observed today. Here, we explore with $n-$body simulations the implications of the low inclinations ($i < 20^{\circ}$) of the distant sednoids (objects with perihelia $q > 40 \mathrm{\; AU}$ and semimajor axes $a > 400 \mathrm{\; AU}$) for the properties of the solar birth cluster. We find that the existence of these orbits, if they were in place in the Sun's birth cluster phase, would limit the product of the stellar number density and the Sun's residence time in the birth cluster to $\lesssim 5 \times 10^3 \mathrm{\; Myr \; pc^{-3}}$, as compared to the weaker limit $\lesssim 5 \times 10^4 \mathrm{\; Myr \; pc^{-3}}$ implied by the low inclinations of the cold classical Kuiper belt.

Observations of the large-scale structure (LSS) provide a powerful test of gravity on cosmological scales, but high-resolution N-body simulations of modified gravity (MG) are prohibitively expensive. We present MG-NECOLA, a convolutional neural network that enhances fast MG-PICOLA simulations to near-N-body fidelity at a fraction of the cost. MG-NECOLA reproduces QUIJOTE-MG N-body results in the power spectrum and bispectrum with better than 1% accuracy down to non-linear scales ($k \simeq 1~h~\mathrm{Mpc}^{-1}$), while reducing computational time by several orders of magnitude. Importantly, although trained only on $f(R)$ models with massless neutrinos, the network generalizes robustly to scenarios with massive neutrinos, preserving accuracy to within 5% at non-linear scales. This combination of precision and robustness establishes MG-NECOLA as a practical emulator for producing large ensembles of high-fidelity simulations, enabling efficient exploration of modified gravity and beyond-$\Lambda$CDM cosmologies in upcoming surveys.

Recent observations suggest that the incidence of near-resonant planets declines as planetary systems age, making young planetary systems key signposts of early dynamical evolution. Here we investigate the dynamical states of three of the youngest multi-transiting planetary systems: AU Mic (3-planet, $\sim$20-Myr-old), V1298 Tau (4-planet, $\sim$23-Myr-old), and TOI-2076 (4-planet, $\sim$200-Myr-old). We find that most planet pairs in these systems lie near resonance with circulating rather than librating resonant angles. As a result, they are more susceptible to dynamical chaos than systems that are either securely locked in resonance or far removed from it. Even modest eccentricities of 0.04 to 0.08 may drive them to instability on timescales of tens to hundreds of Myr. Moreover, the observed orbital architectures are vulnerable to eccentricity excitation through mechanisms such as divergent resonance crossing triggered by planetesimal scattering. The observed near-resonant state may represent a transitional phase between a librating resonant chains and a mature non-resonant planetary system. Finally, we briefly discuss mechanisms that could give rise to the observed near-resonant configurations, including overstable libration, disk turbulence, and receding disk inner edge.

Galaxy mergers can change the rate at which stars are formed. We can trace when these changes occur in simulations of galaxy mergers. However, for observed galaxies we do not know how the star-formation rate (SFR) evolves along the merger sequence as it is difficult to probe the time before or after coalescence. We aim to derive how SFR changes in observed mergers throughout the merger sequence, from a statistical perspective. Merger times were estimated for observed galaxy mergers in the Kilo Degree Survey (KiDS) using a convolutional neural network (CNN). The CNN was trained on mock KiDS images created using IllustrisTNG data. The SFRs were derived from spectral energy density fitting to KiDS and VIKINGs data. To determine the change in SFR for the merging galaxies, each merging galaxy was matched and compared to ten comparable non-merging galaxies; matching redshift, stellar mass, and local density. Mergers see an increase in SFR for galaxies from 300~Myr before the merger until coalescence, continuing until at least 200~Myr after the merger event. After this, there is a possibility that SFR activity in the mergers begins to decrease, but we need more data to better constrain our merger times and SFRs to confirm this. We find that more galaxies with larger stellar mass (M$_{\star}$) have greater SFR enhancement as they merge compared to lower M$_{\star}$ galaxies. There is no clear trend of changing SFR enhancement as local density changes, but the least dense environments have the least SFR enhancement. The increasing SFR enhancement is likely due to closer proximity of galaxies and the presence of more close passes as the time before merger approaches 0~Myr, with SFR slowing 200~Myr after the merger event.

The orbital architectures of compact exoplanet systems record their complicated dynamical histories. Recent research supports the ``breaking-the-chains'' hypothesis, which proposes that compact systems typically form in chains of mean-motion resonances (MMRs) but subsequently break out on a $\sim 100$Myr timescale. We investigate a scenario for breaking the chains through intermittent flybys of planetesimals originating from a distant reservoir. Using $N$-body simulations and semi-analytical calculations, we characterize the disruption of MMRs through these flybys. We find a planetesimal reservoir of total mass $\gtrsim 0.04 M_{\oplus}$ is required to disrupt MMR chains, depending on the mass distribution and the typical number of flybys executed by each planetesimal. We verify that systems disrupted in this way are frequently unstable to close encounters within $\sim 100$Myr of the final flyby. This mechanism operates in systems with both a sufficiently massive reservoir and an efficient mechanism for planetesimal injection. Consequently, we predict an anti-correlation between resonant inner systems and dynamically active outer configurations.

Background: Axion-like particles (ALPs) are hypothetical particles that emerge in numerous theoretical extensions to the Standard Model. Their coupling to electromagnetic field implies that ALPs would mix with photons in the presence of external magnetic fields. As ALP phenomenology is governed by the mass and strength of its coupling, there is a subset of this parameter space in which this mixing would be expected to leave an imprint on the spectra of TeV gamma-ray sources. Data: In 2017, the VERITAS gamma-ray observatory recorded the second day of a dramatic flare of the radio galaxy NGC 1275, embedded at the center of the Perseus galaxy cluster. This serendipitous locale provides a spatially-extended magnetic field of strength O(10$\mu$G) through which escaping photons traverse, making it an excellent target to study ALPs. Methods: We analyze the VERITAS data of NGC 1275's 2017 flare with the gammapy analysis package. Extensive fitting and modeling are performed to ultimately conduct a likelihood analysis used to search for any evidence of a preference for ALPs and to explore the confidence with which constraints can be set. We adopt the CLs method for this study for its conservative approach to setting limits in regimes where the search has limited sensitivity. Results: No evidence for the existence of ALPs is found, and no combination of mass and coupling strength can be excluded at or above 95% confidence level. We provide a map showing the strength of our exclusions in the mass and coupling parameter space. The strongest exclusions are found in the mass range $2 \times 10^{-7}$eV $\lesssim m_a \lesssim 4 \times 10^{-7}$eV and at the coupling strength of $g_{a\gamma} \gtrsim 3 \times 10^{-11}$ GeV$^{-1}$ up to 80% confidence level, which are consistent with previous studies. Conclusions: We find the CLs method to be a trustworthy approach, and advocate for its...

We report the discovery and characterization of three transiting giant planets in the TIC118798035 system. The three planets were identified as transiting candidates from data of the TESS mission, and confirmed with ground-based photometric transit observations along with radial velocity variations obtained with FEROS, HARPS and ESPRESSO. The three planets present transit timing variations (TTVs). We performed a N-body orbital fitting to the TTVs and radial velocities finding that TIC118798035 b is as warm low-density Neptune with a mass of 0.0250$\pm$0.0023 $M_J$, a radius of 0.655$\pm$0.018 $R_J$, and an orbital period of 11.507 d; TIC118798035 c is a warm Saturn with a mass of 0.403$\pm$0.024 $M_J$, a radius of 0.973$\pm$0.023 $R_J$, and an orbital period of 22.564 d; and TIC118798035 d is a warm Jupiter with a mass of 0.773$\pm$0.052 $M_J$, a radius of 0.923$\pm$0.044 $R_J$, and an orbital period of 48.925 d. The bulk metallicities of the three planets don't fully follow the mass-metallicity correlation found for the giant planets of the solar system, which hints at a somewhat different formation history for the planets of the TIC118798035 system. TIC118798035 is the only system having more than two transiting planets larger than 0.5 $R_J$ with a precise orbital and physical characterization, amenable for future atmospheric studies.

Machine learning enables powerful cosmological inference but typically requires many high-fidelity simulations covering many cosmological models. Transfer learning offers a way to reduce the simulation cost by reusing knowledge across models. We show that pre-training on the standard model of cosmology, $\Lambda$CDM, and fine-tuning on various beyond-$\Lambda$CDM scenarios -- including massive neutrinos, modified gravity, and primordial non-Gaussianities -- can enable inference with significantly fewer beyond-$\Lambda$CDM simulations. However, we also show that negative transfer can occur when strong physical degeneracies exist between $\Lambda$CDM and beyond-$\Lambda$CDM parameters. We consider various transfer architectures, finding that including bottleneck structures provides the best performance. Our findings illustrate the opportunities and pitfalls of foundation-model approaches in physics: pre-training can accelerate inference, but may also hinder learning new physics.

Among the emerging excess of massive, bright galaxies at Cosmic Dawn $z \gtrsim 9$ seen by the James Webb Space Telescope, several exhibit spectral features associated with active galactic nuclei (AGN). These AGN candidates suggest that supermassive black holes (SMBHs) grow rapidly in the early Universe. In a series of numerical experiments, we investigate how SMBHs grow within and influence the most massive galaxies at Cosmic Dawn using cosmological hydrodynamic zoom-in simulations run with the adaptive mesh refinement code RAMSES. Our suite of simulations explore how super-Eddington accretion, seed mass, and the strength of feedback influence SMBH-galaxy co-evolution in the early Universe. We find that SMBH growth is sensitive to stellar feedback which generates a turbulent-multiphase interstellar medium (ISM) that stochastically starves the SMBH. In the absence of AGN feedback, we find that the SMBH is starved $\sim 50\%$ of the time after the onset of star formation in the galaxy. SMBH growth can become self-regulated by AGN feedback if the SMBH becomes massive enough, either by accretion or seeding, for its feedback to dominate the surrounding nuclear region. We find no evidence of galaxy-scale, AGN-driven quenching in the star formation rate (SFR) across all simulations in our suite.

The inner region of a subhalo's density distribution is particularly sensitive to dark matter microphysics, with alternative dark matter models leading to both cored and steeply-rising inner density profiles. This work investigates how the lensing signature and detectability of dark matter subhalos in mock HST-, Euclid-, and JWST-like strong lensing observations depends on the subhalo's radial density profile, especially with regards to the inner power-law slope, $\beta$. We demonstrate that the minimum-mass subhalo detectable along the Einstein ring of a system is strongly dependent on $\beta$. In particular, we show that subhalos with $\beta \sim 2.2$ can be detected down to masses over an order-of-magnitude lower than their Navarro-Frenk-White (NFW) counterparts with $\beta \sim 1$. Importantly, we find that the detectability of subhalos with steep inner profiles is minimally affected by increasing the complexity of the main lens galaxy's mass model. This is a unique characteristic of these subhalos, as those with NFW or shallower profiles become essentially undetectable when multipole perturbations are added to the lens model. The results of this work highlight how the underlying dark matter physics can significantly impact the expected number of subhalo detections from strong gravitational lensing observations. This is important for testing Cold Dark Matter against alternatives, such as Self-Interacting Dark Matter, which predict the existence of subhalos with diverse inner density profiles.

Sequential scientific data span many resolutions and domains, and unifying them into a common representation is a key step toward developing foundation models for the sciences. Astronomical spectra exemplify this challenge: massive surveys have collected millions of spectra across a wide range of wavelengths and resolutions, yet analyses remain fragmented across spectral domains (e.g., optical vs. infrared) and object types (e.g., stars vs. galaxies), limiting the ability to pool information across datasets. We present a deep learning model that jointly learns from heterogeneous spectra in a self-supervised manner. Our universal spectral tokenizer processes spectra from a variety of object types and resolutions directly on their native wavelength grids, producing intrinsically aligned, homogeneous, and physically meaningful representations that can be efficiently adapted to achieve competitive performance across a range of downstream tasks. For the first time, we demonstrate that a single model can unify spectral data across resolutions and domains, suggesting that our model can serve as a powerful building block for foundation models in astronomy -- and potentially extend to other scientific domains with heterogeneous sequential data, such as climate and healthcare.

While foundation models have shown promise across a variety of fields, astronomy still lacks a unified framework for joint modeling across its highly diverse data modalities. In this paper, we present AION-1, a family of large-scale multimodal foundation models for astronomy. AION-1 integrates heterogeneous imaging, spectroscopic, and scalar data using a two-stage architecture: modality-specific tokenization followed by transformer-based masked modeling of cross-modal token sequences. The model is pretrained on five large-scale surveys: Legacy Survey, Hyper Suprime-Cam (HSC), Sloan Digital Sky Survey (SDSS), Dark Energy Spectroscopic Instrument (DESI), and Gaia. These span more than 200 million observations of stars, galaxies, and quasars. With a single frozen encoder, AION-1 achieves strong results on a broad suite of downstream tasks, including galaxy and stellar property estimation, galaxy morphology classification, similarity-based retrieval, galaxy image segmentation, and spectral super-resolution. We release AION-1 model variants ranging from 300 M to 3.1 B parameters. Beyond astronomy, AION-1 provides a scalable blueprint for multimodal scientific foundation models that can seamlessly integrate noisy, instrument-specific observations. All code, tokenizers, pretrained weights, and a lightweight evaluation suite are released under an open-source license.

Hydrodynamic simulations have proposed that stellar feedback and bursty star-formation can produce dark matter cores in low-mass galaxies. A key prediction is that feedback-driven gas outflow and inflow cycles can lead to ``breathing modes'' (rapid fluctuations in the global gravitational potential) which drive correlated variations in galaxy size, kinematics, and star-formation rate. In this paper, we test the dynamical effects of feedback-driven breathing modes using a sample of 103 star-forming low-mass galaxies with stellar masses between $7.9<\rm \log M_*/M_\odot<9.6$ and $0.02<z<0.19$. We measure ionized gas velocity dispersions from H$\alpha$ emission lines and compare them to mock observations from the FIRE-2 simulations. We compare gas velocity dispersions ($\rm \sigma_{gas}$), stellar masses, and specific star-formation rates (sSFR). We find a positive correlation between gas velocity dispersion residuals at fixed stellar masses ($\rm \Delta\sigma_{gas}$) and sSFR in both data and simulations. However, the relation is tighter in FIRE-2 compared to the data. FIRE-2 produces more low-sSFR galaxies compared to our observational sample, however, the sSFR distributions agree after limiting both samples to a minimum sSFR. A deeper and more complete photometric sample further indicates that observed low-mass galaxies could span the full range of sSFR predicted in the FIRE-2 simulations. Our results support the existence of short-timescale dynamical effects driven by gas outflow and inflow cycles in low-mass galaxies and motivate additional tests of the breathing mode.

Observations by JWST have confirmed the presence of supermassive black holes (BHs) at redshifts $z\gtrsim10$, lending support to scenarios in which BHs experience rapid growth through intense gas accretion. Here we investigate the growth of a BH embedded at the center of a quasi-star, a theoretically predicted object formed via direct collapse. In a quasi-star, the central BH accretes at a highly super-Eddington rate, while the excess energy is transported outward by convection and radiated at approximately the Eddington luminosity of the entire star. We employ the open-source stellar evolution code \texttt{MESA} to construct quasi-star models and follow the time-dependent growth of the central BH under different prescriptions for the accretion rate at the inner boundary $R_i$, and further considering the effect of winds. For the case $R_i=NR_{\rm B}$, where $N$ is a constant and $R_{\rm B}$ is the Bondi radius corresponding to the mass of the BH and the gas infalling onto it, our models terminate when the BH mass reaches a critical value $M_{\mathrm{crit}}(N)=c_{s,i}^3/(12\sqrt{N^3G^3\pi\rho_i})$ (where $c_{s,i}$ and $\rho_i$ are the sound speed and density at $R_i$, respectively), a limit we also derive analytically. Models that feature an inner convective region matched to an outer adiabatic envelope exhibit BH growth up to approximately $M_{\mathrm{BH}}/M_\star\simeq 0.33$, largely independent of the stellar mass $M_\star$ itself. This ratio is approximately preserved even in the presence of mass loss, as several properties of the model are independent of the quasi-star's total mass.

Double-peaked Lyman-$\alpha$ (Ly$\alpha$) profiles provide critical insights into gas kinematics and the distribution of neutral hydrogen (HI) from the interstellar to the intergalactic medium (ISM to IGM), and serve as valuable diagnostics of ionising Lyman continuum (LyC) photon escape. We present a study of the global and spatially resolved properties of double-peaked Ly$\alpha$ emitters (LAEs) based on VLT/MUSE data from the MAGPI survey. From a parent sample of 417 LAEs at z = 2.9 - 6.6 in the first 35 fields, we identify 108 double-peaked LAEs using an automated peak classification technique. We measure a double-peak fraction of $\sim37\%$ at $z < 4$, decreasing to $\sim14\%$ at $z > 4$, likely due to enhanced IGM attenuation. Approximately $17\%$ of the double-peaked LAEs are blue-dominated, suggesting gas inflows. The blue-to-total flux ratio exhibits a luminosity dependence: fainter lines generally show higher blue flux, though some luminous sources also show strong blue peaks. We find a narrowing of the red peak at $z > 4$, despite the presence of the blue peak, indicating intrinsic galaxy evolution rather than IGM attenuation. Several LAEs exhibit residual flux in the absorption trough, with normalised trough flux anticorrelating with peak separation, reflecting variations in HI column density. We further investigate spatially resolved properties of ten red-dominated LAEs with extended Ly$\alpha$ halos. Despite azimuthal variations, both the blue-to-total flux ratio and normalised trough flux density increase with radius, while peak separation decreases. The red peak asymmetry shows only minor radial changes. These trends are consistent with variations in shell outflow velocity and HI column density across the halos. Based on peak separation, red peak asymmetry, and residual trough flux, we identify five LAEs as strong LyC-leaker candidates.

The China Jinping Underground Laboratory, characterized by a vertical rock overburden of 2,400 m, provides an exceptionally effective shield against cosmic muons with energies below 3 TeV. The surviving high-energy muons, produced as part of extensive air showers, open a unique observational window into primary cosmic rays with energies ranging from tens of TeV up to the PeV scale and beyond. This distinctive feature also enables detailed studies of the earliest stages of shower development. Using 1,338.6 live days of data collected with a one-ton prototype detector for the Jinping Neutrino Experiment, we measured the underground muon flux originating from air showers. The results show discrepancies of about 40%, corresponding to a significance of more than 5.5$\sigma$, relative to predictions from several leading hadronic interaction models. We interpret these findings from two complementary perspectives: (i) by adopting the expected cosmic ray spectra, we constrain the modeling of the initial hadronic interactions in air showers; and (ii) by assuming specific hadronic interaction models, we infer the mass composition of cosmic rays, and our data favor a lighter component in the corresponding energy range. Our study demonstrates the potential of deep underground laboratories to provide new experimental insights into cosmic rays.

Interstellar neutral (ISN) atoms enable studies of the physical conditions in the local interstellar medium surrounding the heliosphere. ISN helium, which is the most abundant species at 1 au, is directly observed by space missions, such as Interstellar Boundary Explorer (IBEX). However, some of these atoms are ionized by solar ultraviolet radiation before reaching 1 au, producing pickup ions (PUIs). A recent analysis of IBEX data suggests that the helium photoionization rates predicted by models are underestimated by up to 40%. The Solar Wind Around Pluto (SWAP) instrument on board New Horizons enables the study of PUIs giving complementary insight into the other side of the ionization process. Our goal is to verify this increased helium ionization by determining the ionization rate of ISN helium in the heliosphere based on the SWAP observations of helium PUIs. For this purpose, we analyze SWAP data collected between 2012 and 2022, at distances 22 to 54 au from the Sun. We develop a new method for fitting model distribution functions to the observational data using the maximum likelihood method. Our approach accounts for the spacecraft's rotation and the SWAP response function, which depends on both energy and inflow direction. We estimate SWAP's efficiency for helium relative to that for hydrogen and determine the ISN helium ionization rate. We find that the photoionization rate obtained from the SWAP observations is 43% larger than the rates predicted by models, confirming the IBEX results.

There exist extremely massive spiral galaxies in isolated environments, with stellar masses several times that of the Milky Way, yet their star formation rates (SFRs) are comparable to or even lower than that of the Milky Way. In this paper, we investigate the molecular gas properties of such galaxies to better understand the origin of their low SFRs. We present IRAM 30m CO observations of five extremely massive spirals from the CGM-MASS sample. We compare their star formation efficiencies (SFEs) with the Kennicutt-Schmidt relation and find that these massive spirals generally exhibit low efficiency in converting molecular gas into stars. We further compare their molecular gas masses with their atomic gas and stellar masses, and also include the CHANG-ES sample galaxies observed with the IRAM 30m telescope in a similar manner for comparison. Our sample galaxies show low efficiency in converting atomic to molecular gas and have lower molecular gas fractions, suggesting that their suppressed star formation stems from both limited gas supply and inefficient star formation. Considering potential cold gas sources in massive spirals, we argue that their current reservoirs likely originate from past starburst or merger events rather than ongoing accretion in present isolated environments. Finally, we examine the location of these galaxies on the baryonic Tully-Fisher relation, finding them baryon-deficient and deviating from the trend of lower-mass galaxies. This suggests either a significant undetected baryonic component or a flattening/turnover of the relation at the high-mass end, consistent with the stellar mass-halo mass relation.

Utilizing cosmological hydrodynamic simulations we show that there is a brief super-Eddington accretion phase in typical halos at high redshift, impervious to AGN self-regulation. However, once having attained a black hole mass of $10^4-10^5\msun$, AGN feedback process can self-regulate to guide the SMBHs to grow at a significantly slower, sub-Eddington rate. By redshift $z\sim 10$ the black hole mass with an initial super-Eddington jump-start is caught up by that in the case with a steady Eddington limited case. Thus a continuous Eddington limit case represents the fastest possible route to maximally grow SMBHs. To account for the observed $z=7-10$ quasars with supermassive black holes of billions of solar masses, our analysis establishes firmer ground for the need of seed masses of $10^4-10^5\msun$ that are not grown via an earlier super-Eddington phase.

We presented a detailed analysis of seven thermonuclear X-ray bursts from Terzan 5 X-3/Swift J174805.3-244637, detected by NICER during the source's 2023 outburst. Our analysis reveals a clear evolution of burst properties, identifying four non-photospheric radius expansion (non-PRE) bursts, one PRE candidate occurring in a mixed hydrogen/helium environment, and two powerful PRE bursts from pure helium ignition. The time-resolved burst spectra were well described by a model including a variable persistent emission component, quantified by a factor $f_a$, due to the Poynting-Robertson drag. The strength of this interaction scales with burst luminosity: the enhancement is absent ($f_a \approx 1$) in the faintest bursts, becomes modest ($f_a \approx 1.5-2$) for the more luminous non-PRE burst and the PRE candidate, and is very strong ($f_a \approx 6-8$) during the pure-helium PRE bursts. This observed transition from mixed-fuel to pure-helium burning as the local mass accretion rate dropped below $\sim$10% of the Eddington limit, $\dot{m}_{\rm Edd}$, aligns with theoretical predictions. We verified this scenario with two independent methods. First, at the known distance to Terzan 5, the touchdown luminosities of both the pure helium PRE bursts and the mixed-fuel PRE candidate are consistent with reaching their respective, composition-dependent Eddington limits on the same plausible, massive neutron star of $\sim 2 M_\odot$. Second, the observed recurrence times of the non-PRE bursts were consistent with predictions for mixed-fuel burning.

Decaying pulsations have been simultaneously detected in the low-energy X-rays of solar/stellar flares, which are supposed to be associated with standing slow magnetoacoustic or kink-mode waves. The physical mechanism behind rapidly decaying remains unknown. We present the detection of quasi-periodic pulsations (QPPs) with rapidly decaying in high-energy emissions produced in two major flares on 10 January and 14 May 2024. Using empirical mode decomposition, decaying QPPs are identified in hard X-ray and microwave emissions during the flare impulsive phase, suggesting a process of oscillatory magnetic reconnection. The quasi-periods and decay times are determined by a damped harmonic function, which are approximately 177$\pm$8 s (249$\pm$25 s) and 118$\pm$4 s (124$\pm$5 s), respectively. The restructured X-ray images reveal double footpoints connected by hot flare loops. Their phase speeds are estimated to about 400 km/s and 670 km/s, both below the local sound speed in high-temperature plasmas, indicating the presence of slow-mode waves in hot flare loops. We perform coronal diagnostics based on standing slow-mode waves and derive key physical parameters, including the polytropic index, the thermal ratio, viscous ratio and radiation ratio, which are consistent with previous results. Our observations support that the decaying QPPs are triggered by oscillatory magnetic reconnection that is modulated by standing slow magnetoacoustic waves, with their rapid decay attributable to a co-effect of viscous damping and localized magnetic reconnection rate.

Detecting coherent radio bursts from nearby M dwarfs provides opportunities for exploring their magnetic activity and interaction with orbiting exoplanets. However, it remains uncertain if the emission is related to flare-like activity similar to the Sun or magnetospheric process akin to magnetized planets. Using observations (1.0 - 1.5 GHz) taken by the Five-hundred-meter Aperture Spherical radio Telescope, we found a type of millisecond-scale radio bursts with exceptionally high frequency drift rates ($\sim 8\;\rm{GHz\;s^{-1}}$) from an active M dwarf, AD Leo. The ultrafast drift rates point to a source region with a notably low magnetic scale height ($<0.15\; r_\star$, $r_\star$ as the stellar radius), a feature not expected in a commonly assumed dipole-like global field but highly possible in localized strong-field structures, i.e. starspots. Our findings suggest that a concentrated magnetic field above starspots could be responsible for some of the most intense radio bursts from M dwarfs, supporting a solar-like electron acceleration mechanism.

We compare the performance of the flat-sky approximation and Limber approximation for the clustering analysis of the photometric galaxy catalogue of Euclid. We study a 6 bin configuration representing the first data release (DR1) and a 13 bin configuration representative of the third and final data release (DR3). We find that the Limber approximation is sufficiently accurate for the analysis of the wide bins of DR1. Contrarily, the 13 bins of DR3 cannot be modelled accurately with the Limber approximation. Instead, the flat-sky approximation is accurate to below $5\%$ in recovering the angular power spectra of galaxy number counts in both cases and can be used to simplify the computation of the full power spectrum in harmonic space for the data analysis of DR3.

We present the formulation, algorithm and numerical tests of the magnetohydrodynamic-particle-in-cell (MHD-PIC) method with particles treated under the guiding center approximation, which we term the MHD-gPIC method, and it is implemented in the Athena++ MHD code. The new MHD-gPIC model consists of thermal (cold) fluid and high-energy particles whose dynamics are integrated through guiding center equations including drift motion, with carefully evaluated source terms as particle backreaction. The code is validated with a series of tests, and it is expected to be primarily applicable to study particle acceleration and transport in systems where gyro-resonance is considered insignificant. We also present preliminary studies of particle acceleration during non-relativistic magnetic reconnection.

Baryonic feedback remains one of the largest uncertainties in cosmological hydrodynamical simulations, with different prescriptions producing divergent predictions for the fraction of gas expelled from halos, the radial extent of the gas expulsion and the impact on large scale matter clustering. We present the first systematic study of the kinetic Sunyaev-Zel'dovich (kSZ) effect across a wide range of simulations (FLAMINGO, ANTILLES, BAHAMAS, SIMBA, FABLE and their variants), and compare them directly to DESI Year 1 + ACT kSZ measurements. We ensure a like-for-like comparison with observations by developing a robust methodology that accounts for the halo mass selection using galaxy-galaxy lensing, cosmic variance, miscentering and satellites, establishing the kSZ effect as a new benchmark for the simulations. We find that fiducial feedback models are disfavoured by >3 sigma, while simulations with more powerful AGN feedback within the FLAMINGO and BAHAMAS suites, as well as SIMBA, reproduce the observed kSZ signal within <2 sigma. We use the ANTILLES simulation suite to demonstrate that the amplitude of the kSZ effect is a strong predictor of matter power spectrum suppression, competitive with baryon fraction metrics. These results establish the kSZ as a critical probe for evaluating feedback physics and for advancing the fidelity of cosmological simulations.