Papers for Tuesday, Jun 10 2025

The Asteroid Terrestrial-impact Last Alert System (ATLAS) is an all-sky optical sky survey with a cadence of 24 to 48 hours and the ATLAS Transient Server processes the alert stream to enable the discovery and follow-up of extra-galactic transients. The data from the ATLAS server can be accessed through a REST API, which has allowed the development of bots that need direct access to the data to help rank alerts and trigger followup observations of promising targets. Here we present the python client we have developed for the ATLAS API to help connect bots and scientists to our data.

Red giants undergo dramatic and complex structural transformations as they evolve. Angular momentum is transported between the core and envelope during this epoch, a poorly understood process. Here, we infer envelope and core rotation rates from Kepler observations of $\sim$1517 red giants. While many measurements are consistent with the existing studies, our investigation reveals systematic changes in the envelope-to-core rotation ratio and we report the discovery of anomalies such as clump stars with rapidly rotating cores, and red giants with envelopes rotating faster than their cores. We propose binary interactions as a possible mechanism by which some of these cores and envelopes are spun up. These results pose challenges to current theoretical expectations and can have major implications for compact remnants born from stellar cores.

We present extensive ultraviolet (UV), optical, and near-infrared (NIR) photometric and spectroscopic observations of the nearby hydrogen-poor superluminous supernova (SLSN-I) SN2024rmj at z = 0.1189. SN 2024rmj reached a peak absolute magnitude of Mg $\approx$ -21.9, placing it at the luminous end of the SLSN-I distribution. The light curve exhibits a pronounced pre-peak bump ($\approx$ 60 d before the main peak) and a post-peak bump ($\approx$ 55 d after the main peak). The bulk of the light curve is otherwise well fit by a magnetar spin-down model, with typical values (spin: $\approx$ 2.1 ms; magnetic field: $\approx$ 6 $\times$ 10$^{13}$ G; ejecta mass: $\approx$ 12 M$_\odot$). The optical spectra exhibit characteristic SLSN-I features and evolution, but with a relatively high velocity of $\approx$ 8,000 km s$^{-1}$ post-peak. Most significantly, we find a clear detection of helium in the NIR spectra at He I $\lambda$1.083 $\mu$m and $\lambda$2.058 $\mu$m, blueshifted by $\approx$ 15,000 km s$^{-1}$ (13 d before peak) and $\approx$ 13,000 km s$^{-1}$ (40 d after peak), indicating that helium is confined to the outermost ejecta; based on these NIR detections, we also identify likely contribution from He I $\lambda$5876 Å in the optical spectra on a similar range of timescales. This represents the most definitive detection of helium in a bright SLSN-I to date, and indicates that progenitors with a thin helium layer can still explode as SLSNe.

We analyze two dusty star-forming galaxies at $z=6.6$. These galaxies are selected from the ASPIRE survey, a JWST Cycle-1 medium and ALMA Cycle-9 large program targeting 25 quasars and their environments at $z\simeq6.5 - 6.8$. These galaxies are identified as companions to UV-luminous quasars and robustly detected in ALMA continuum and [C II] emission, yet they are extraordinarily faint at the NIRCam wavelengths (down to $>28.0$ AB mag in the F356W band). They are more obscured than galaxies like Arp220, and thus we refer to them as "NIRCam-dark" starburst galaxies (star formation rate $\simeq 80 - 250\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$). Such galaxies are typically missed by (sub)-millimeter blank-field surveys. From the star-formation history (SFH), we show that the NIRCam-dark galaxies are viable progenitors of massive quiescent galaxies at $z\gtrsim4$ and descendants of UV-luminous galaxies at $z>10$. Although it is hard to constrain their number density from a quasar survey, we conclude that NIRCam-dark galaxies can be as abundant as $n\sim10^{-5.5}$ Mpc$^{-3}$ assuming a light halo occupation model. If true, this would equal to $\sim$30% of the number densities of both the quiescent galaxies at $z\gtrsim4$ and UV-luminous galaxies at $z>10$. We further predict that analogs at $z\sim8$ should exist according to the SFH of early massive quiescent galaxies. However, they may fall below the current detection limits of wide JWST and ALMA surveys, thus remaining "JWST-dark". To fully trace the evolution of massive galaxies and dust-obscured cosmic star formation at $z\gtrsim8$, wide-field JWST/NIRCam imaging and slitless spectroscopic surveys of early protoclusters are essential.

We present a systematic search for past supernovae (SNe) and other historical optical transients at the positions of fast radio burst (FRB) sources to test FRB progenitor systems. Our sample comprises 83 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and its KKO Outrigger, along with 93 literature FRBs representing all known well-localized FRBs. We search for optical transients coincident in position and redshift with FRBs and find no significant associations within the 5-sigma FRB localization uncertainties except for a previously identified potential optical counterpart to FRB 20180916B. By constraining the timescale for SN ejecta to become transparent to FRB emission, we predict that it takes at least 6-10 years before the FRB emission can escape. From this, we infer that approximately 7% of matched optical transients, up to 30% of currently known SNe, and up to 40% of core-collapse SNe could have an observable FRB based on timescales alone. We derive the number of new, well-localized FRBs required to produce one FRB-SN match by chance, and find it will take ~ 22,700 FRBs to yield one chance association at the projected CHIME/FRB Outrigger detection rate. Looking forward, we demonstrate redshift overlap between SNe detected by the upcoming Vera C. Rubin Observatory and CHIME/FRB Outrigger FRBs, indicating the prospect of an increase in potential associations at redshift z < 1. Our framework is publicly available, flexible to a wide range of transient timescales and FRB localization sizes, and can be applied to any optical transient populations in future searches.

We present a new framework for massive star evolution that is no longer driven by Dutch or other mass-loss rate Recipes, but which take the physics of Gamma or L/M dependent mass loss consistently into account. We first discuss the hot-star mass-loss kink and the transition mass loss rate between optically thin and thick winds, before discussing vertical stellar evolution, mass evaporation, and the maximum black hole (BH) mass. We end with a suggestion that a recently uncovered red supergiant (RSG) kink might be related to similar underlying L/M physics as the hot-star kink.

MeerKAT observations of the recently discovered, extremely low mass galaxy, Pavo, have revealed a neutral gas (HI) reservoir that was undetected in archival HI single dish data. Based on the latest distance measurement ($2.16^{+0.08}_{-0.07}$ Mpc), we measure Pavo's HI mass as $\log M_\mathrm{HI}/\mathrm{M_\odot} = 5.79 \pm 0.05$, making it the lowest mass HI reservoir currently known in an isolated galaxy (with a robust distance measurement). Despite Pavo's extreme isolation, with no known neighbor within over 700 kpc, its HI reservoir is highly disturbed. It does not show clear signs of rotation and its center of mass is offset from the stellar body center by 320 pc, while its peak is offset by 82 pc (both in projection). Despite this disturbed morphology, Pavo still appears to be consistent with the HI size--mass relation, although it is not possible to accurately determine a suitable inclination correction. Such disturbed, offset and disorganized HI reservoirs are predicted by simulations of low-mass, star-forming dwarfs in which supernova-driven outflows efficiently disrupt the interstellar medium after a star formation event. It is likely that we are witnessing Pavo in precisely this period, tens to a few hundred Myr after a star formation episode, when internal feedback has disrupted its gas reservoir.

ESA's Ariel mission will be uniquely suited to performing population-level studies of exoplanets. Most of these studies consist of quantifying trends between an Ariel-measured quantity, y, and an a priori planetary property, x; for example, atmospheric metallicity as inferred from Ariel transit spectroscopy vs. planetary radius. We define the leverage of a survey with N targets as L = sqrt(N)stdev(x) and show that it quantitatively predicts the precision of population-level trends. The target selection challenge of Ariel can therefore be summarized as maximizing L along some axes of diversity for a given cumulative observing time. To this end, we consider different schemes to select the mission reference sample for a notional three year transit spectroscopy survey with Ariel. We divide the exoplanets in the mission candidate sample into logarithmic classes based on radius, equilibrium temperature, and host star temperature. We then construct a target list by cyclically choosing the easiest remaining target in each class. We find that in many cases the leverage is greatest for a single class: dividing planets into multiple classes increases the diversity of targets, but reduces their numbers. The leverage on a single axis of diversity can be increased by dividing that axis into many planet classes, but this sacrifices leverage along other axes of diversity. We conclude that a modest number of classes, possibly only one, should be defined when selecting Ariel targets. Lastly, we note that the statistical leverage of the Ariel transit survey would be significantly increased if current candidate planets were confirmed. This highlights the urgency of vetting and confirming the easiest transmission and emission spectroscopy targets in the Ariel mission candidate sample.

Millimeter continuum spiral arms have so far only been detected in a handful of protoplanetary disks, and thus we have a limited understanding of the circumstances in which they can form. In particular, substructures in small disks ($R\lessapprox 50$ au) have not been well-characterized in comparison with large disks. We present ALMA 1.3 mm continuum observations of the disk around the T Tauri star Haro 6-13 at a resolution of $\sim0.04''$ ($\sim5$ au). A pair of low-contrast spiral arms are detected at disk radii from $\sim10-35$ au. They can be approximated as Archimedean spirals with pitch angles ranging from $\sim10-30^\circ$. The low value of the disk-averaged spectral index between 1.3 and 3 mm ($\alpha=2.1$) and the high brightness temperatures suggest that the millimeter continuum is likely optically thick and thus may hide sufficient mass for the disk to become gravitationally unstable and form spiral arms. CO observations have shown that Haro 6-13 is surrounded by an envelope, raising the possibility that infall is facilitating spiral arm formation.

High-redshift ($z\sim1$) galaxy clusters are the domain where environmental quenching mechanisms are expected to emerge as important factors in the evolution of the quiescent galaxy population. Uncovering these initially subtle effects requires exploring multiple dependencies of quenching across the cluster environment, and through time. We analyse the stellar-mass functions (SMFs) of 17 galaxy clusters within the GOGREEN and GCLASS surveys between $0.8<z<1.5$, and with $\log{(M/{\rm{M_\odot}})}>9.5$. The data are fit simultaneously with a Bayesian model that allows the Schechter function parameters of the quiescent and star-forming populations to vary smoothly with cluster-centric radius and redshift. The model also fits the radial galaxy number density profile of each population, allowing the global quenched fraction to be parameterised as a function of redshift and cluster velocity dispersion. We find the star-forming SMF to not depend on radius or redshift. For the quiescent population however, there is $\sim2\sigma$ evidence for a radial dependence. Outside the cluster core ($R>0.3\,R_{\rm200}$), the quenched fraction above $\log{(M/{\rm{M_\odot}})}=9.5$ is $\sim40{\rm\;per\,cent}$, and the quiescent SMF is similar in shape to the star-forming field. In contrast, the cluster core has an elevated quenched fraction ($\sim70{\rm\;per\,cent}$), and a quiescent SMF similar in shape to the quiescent field population. We explore contributions of 'early mass-quenching' and mass-independent 'environmental-quenching' models in each of these radial regimes. The core is well-described primarily by early mass-quenching, which we interpret as accelerated quenching of massive galaxies in protoclusters, possibly through merger-driven feedback mechanisms. The non-core is better described through mass-independent, environmental-quenching of the infalling field population.

Spot-crossing transits offer a unique opportunity to probe spot properties such as temperature, size, and surface distribution. TOI-3884 is a rare system in which spot-crossing features are persistently observed during every transit. This is due to its unusual configuration: a polar-orbiting super-Neptune transits a pole-on mid-M dwarf, repeatedly crossing a polar spot. However, previous studies have reported discrepant values in key system parameters, such as stellar inclination and obliquity. To address this, we conducted multi-band, multi-epoch transit observations of TOI-3884b using the MuSCAT instrument series, along with photometric monitoring with the LCO 1m telescopes/Sinistro. We detected time-dependent variations in the spot-crossing signals, indicating that the spot is not exactly on the pole. From the monitoring data, we measured a stellar rotation period of $11.043_{-0.053}^{+0.054}$ days with a modulation amplitude of $\sim$5% in the r-band, consistent with the time variability in the spot-crossing features. Our analysis reconciles previous discrepancies and improves the constraints on the parameters of the system geometry ($i_\star = 40.1^{+2.0}_{-1.2}$ deg and $\lambda =139.0^{+9.0}_{-3.7}$ deg) and those of the spot properties (spot radius of $0.425_{-0.011}^{+0.018} R_\star$ and spot-photosphere temperature difference of $200_{-9}^{+11}$ K). These results provide a critical context for interpreting upcoming transmission spectroscopy of TOI-3884b, as well as yielding new insights into the magnetic activity and spin-orbit geometry of M dwarfs.

We present deep optical and near-infrared observations of the host galaxies of 34 fast radio bursts (FRBs) detected by the Commensal Real-time ASKAP Fast Transient (CRAFT) survey on the Australian SKA Pathfinder (ASKAP) to compare the locations of FRBs relative to their host light distributions. Incorporating three additional FRBs from the literature, for a total of four repeating and 33 apparently non-repeating FRBs, we determine their projected galactocentric offsets and find a median of $ 4.2^{+5.7}_{-2.5}$ kpc ($1.0^{+1.5}_{-0.6}r_e$). We model their host surface brightness profiles and develop synthetic spatial distributions of their globular clusters based on host properties. We calculate the likelihood the observed location of each FRB is consistent with the smooth light of its host galaxy, residual (primarily spiral) substructure, or globular cluster distributions. The majority of FRBs favor locations within the disks of their galaxies, while only 7-13\% favor a globular cluster origin, primarily those with galactocentric offsets $\gtrsim3r_e$. At $z<0.15$, where spiral structure is apparent in 86\% of our sample of FRB hosts, we find $\approx 20-46\%$ of FRBs favor an association with spiral arms. Assuming FRBs derive from magnetars, our results support multiple formation channels with the majority of progenitors associated with massive stars and a minority formed through dynamical channels. However, the moderate fraction of FRBs associated with spiral structure indicates that high star formation efficiency of the youngest and most massive stars is not a predominant driver in the production of FRB progenitors.

Predicting nuclear masses is a longstanding challenge. One path forward is machine learning (ML) which trains on experimental data, but can suffer large errors when extrapolating toward neutron-rich species. In nature, such masses shape observables for the rapid neutron capture process (r-process), which in principle could inform ML models. Here we introduce a multi-objective optimization approach using the Pareto Front algorithm. We show that this technique, capable of identifying models which generate r-process abundances aligning with both Solar and stellar data, is a promising method to select ML models with reliable extrapolation power.

The cores of neutron stars (NSs) contain the densest matter in the universe. Rapid advancements in neutron-star observations allow unprecedented empirical access to cold, ultra-dense Quantum Chromodynamics (QCD) matter. The combination of these observations with theoretical calculations has revealed previously inaccessible features of the equation of state (EoS) and the QCD phase diagram. In this thesis, I demonstrate how perturbative-QCD calculations at asymptotically high baryon densities provide robust constraints on the EoS at neutron-star densities. The method for constraint propagation is based solely on thermodynamical causality, stability, and consistency of the EoS. By constructing a large ensemble of EoSs using Gaussian processes regression and incorporating it into a Bayesian inference of EoS, I demonstrate that the novel pQCD constraints go beyond those obtained from current astrophysical observations alone, forcing the EoS to soften at the maximum densities of stable neutron stars. This softening of the EoS can be interpreted as an indication of approximate conformal symmetry restoration, a sign of a first-order phase transition (FOPT), or potentially both. I show that the conformal symmetry restoration is consistent with the hypothesis of quark matter cores inside the most massive NSs. Although current astrophysical data and theoretical inputs cannot definitively distinguish between the two scenarios, they slightly favor the occurrence of a phase transition of some kind - whether a crossover to quark matter or a destabilizing FOPT - in the cores of the most massive neutron stars.

A significant number of Lyman-break galaxies (LBGs) with redshifts 3 < z < 5 are expected to be observed by the upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This will enable us to probe the universe at higher redshifts than is currently possible with cosmological galaxy clustering and weak lensing surveys. However, accurate inference of cosmological parameters requires precise knowledge of the redshift distributions of selected galaxies, where the number of faint objects expected from LSST alone will make spectroscopic based methods of determining these distributions extremely challenging. To overcome this difficulty, it may be possible to leverage the information in the large volume of photometric data alone to precisely infer these distributions. This could be facilitated using forward models, where in this paper we use stellar population synthesis (SPS) to estimate uncertainties on LBG redshift distributions for a 10 year LSST (LSSTY10) survey. We characterise some of the modelling uncertainties inherent to SPS by introducing a flexible parameterisation of the galaxy population prior, informed by observations of the galaxy stellar mass function (GSMF) and cosmic star formation density (CSFRD). These uncertainties are subsequently marginalised over and propagated to cosmological constraints in a Fisher forecast. Assuming a known dust attenuation model for LBGs, we forecast constraints on the sigma8 parameter comparable to Planck cosmic microwave background (CMB) constraints.

We perform the first direct consistency check between the recently measured proton spectrum at the knee by LHAASO and the collaboration's own high-precision mapping of Galactic diffuse gamma-ray emission. By modeling the hadronic gamma-ray production using the updated cosmic-ray spectra, gas templates and cross-section models, we show that the predicted gamma-ray flux robustly overshoots the LHAASO data in both inner and lateral Galactic regions. This persistent mismatch in both normalization and spectral shape challenges conventional scenarios linking the local cosmic-ray sea to Galactic gamma-ray emission, and calls for a revision of current cosmic ray models in the TeV-PeV sky.

We present a modification of the expected local primordial black hole (PBH) count $N(y)$, typically seen in the context of the early PBH binary merger rate as a term in the merger rate suppression. We utilize recent results in small-scale PBH clustering to formulate $N(y)$ in such a way that accounts for variations in the binary mass ratio $q$. We then examine how this change affects the projected constraints on PBH abundance from simulated Einstein Telescope (ET) and LISA mergers. Our results indicate that for broadly extended mass distributions, the merger suppression is greatly reduced for binaries with $q \gg 1$. This leads to an enhanced merger rate for binary distributions favoring a lighter average mass. This change is best reflected in the constraints derived from the stochastic gravitational wave background (SGWB), as the increased merger rate does not have as much of an impact on constraints from resolvable mergers. Our results imply that the assumption of PBH clustering is testable at abundances and mass ranges much lower than anticipated, but only definitively through measurements of the SGWB.

A random Gaussian density field contains a fixed amount of Fisher information on the amplitude of its power spectrum. For a given smoothing scale, however, that information is not evenly distributed throughout the smoothed field. We investigate which parts of the field contain the most information by smoothing and splitting the field into different levels of density (using the formalism of indicator functions), deriving analytic expressions for the information content of each density bin in the joint-probability distribution (given a distance separation). We find that the information peaks at moderately rare densities (where the number of smoothed survey cells $N\sim 100$). Counter-intuitively, we find that, for a finite survey volume, indicator function analysis can outperform conventional two-point statistics while using only a fraction of the total survey cells, and we explain why. In light of recent developments in marked statistics (such as the indicator power spectrum and density-split clustering), this result elucidates how to optimize sampling for effective extraction of cosmological information.

The Resolve instrument on the X-ray Imaging and Spectroscopy Mission (XRISM) uses a 36-pixel microcalorimeter designed to deliver high-resolution, non-dispersive X-ray spectroscopy. Although it is optimized for extended sources with low count rates, Resolve observations of bright point sources are still able to provide unique insights into the physics of these objects, as long as high count rate effects are addressed in the analysis. These effects include {the loss of exposure time for each pixel}, change on the energy scale, and change on the energy resolution. To investigate these effects under realistic observational conditions, we observed the bright X-ray source, the Crab Nebula, with XRISM at several offset positions with respect to the Resolve field of view and with continuous illumination from {$^{55}$Fe sources} on the filter wheel. For the spectral analysis, we excluded data where exposure time loss was too significant to ensure reliable spectral statistics. The energy scale at 6 keV shows a slight negative shift in the high-count-rate regime. The energy resolution at 6 keV worsens as the count rate in electrically neighboring pixels increases, but can be restored by applying a nearest-neighbor coincidence cut (``cross-talk cut''). We examined how these effects influence the observation of bright point sources, using GX 13+1 as a test case, and identified an eV-scale energy offset at 6 keV between the inner (brighter) and outer (fainter) pixels. Users who seek to analyze velocity structures on the order of tens of km~s$^{-1}$ should account for such high count rate effects. These findings will aid in the interpretation of Resolve data from bright sources and provide valuable considerations for designing and planning for future microcalorimeter missions.

The evolution of low-mass galaxies is shaped by both internal processes and environmental factors, yet the role of environment and mergers in regulating their growth and star formation rates remains poorly understood, especially in the low-density regime. This study aims to compare the evolutionary pathways and merger histories of star-forming and quenched galaxies in dense (cluster) and under-dense (void) environments, focusing on galaxies with stellar masses in the range $10^{8.5} \leq M_\star/M_\odot \leq 10^{10.5}$. It presents the first statistical analysis within this stellar mass range, explicitly distinguishing between mini, minor, and major mergers across varying environments. Using the high-resolution TNG300-1 simulation from the IllustrisTNG project, we classify galaxies as star-forming and quenched based on sSFR and UVJ criteria. We track their physical properties over the last $\sim 11.5$ Gyr ($z < 3$), follow their complete merger histories, distinguishing between major, minor, and mini mergers, and assess the statistical impact of these mergers on star formation, gas content, and the mass assembly of their supermassive black holes.

We explore the inflationary evolution and observational viability of the $\alpha$-attractor E-model in the framework of warm inflation, focusing on both weak and strong dissipative regimes, with a dissipation coefficient linear in temperature. In the strong regime, we account for the growth of inflaton fluctuations due to coupling with the radiation bath via two different forms for the dissipation enhancement function: one associated with plateau-like potentials, and another motivated by the warm little inflation scenario. Employing slow roll conditions, we analytically derive the expressions for the key inflationary observables, the spectral index $n_s$ and the tensor-to-scalar ratio $r$, in both dissipative regimes. The resulting theoretical trajectories on the $n_s$--$r$ plane are then juxtaposed with the contour plots obtained from Planck 2018 data in order to constrain the model parameter. Our analysis shows that the warm $\alpha$-attractor E-model remains compatible with observations in both dissipative regimes, with dissipation playing a crucial role in shifting the predictions and enlarging the viable parameter space, highlighting observational robustness of the model when extended to warm inflation.

In recent years the idea, first formulated many decades ago, that the Be phenomenon could be causally related to the duplicity of Be stars, has been repeatedly reconsidered from various perspectives. It is important, therefore, to have reliable information on Be stars, which are confirmed members of binary systems. This study is devoted to V1624 Cyg = 28 Cygni, which was recently identified as a binary with a compact secondary. By measuring the radial velocities (RVs) of the wings of the H alpha emission line and using archival data and published RVs from the International Ultraviolet Explorer, we demonstrate that the Be primary moves in the 359.26 d orbit found recently from interferometry. Our preliminary radial-velocity solution leads to binary masses of 5.6, and 0.66 solar masses. Moreover, we documented large and irregular spectral, brightness, and colour changes over a time interval of several decades to show that the object never completely lost its circumstellar matter.

Reverberation mapping (RM) can measure black hole accretion disc sizes and radial structure through observed time lags that should increase with wavelength as $\tau\propto\lambda^{4/3}$. Our 250-day RM campaign on NGC 7469 combines sub-day cadence 7-band photometry from the Las Cumbres Observatory robotic telescopes and weekly X-ray and UVOT data from Swift. By fitting these light curves, we measure the spectral energy distribution of the variable accretion disc and inter-band lags of just 1.5 days across the UV to the optical range. The disc SED is close to the expected $f_\nu\propto\nu^{1/3}$, and the lags are consistent with $\tau\propto \lambda^{4/3}$, but three times larger than expected. We consider several possible modifications to standard disc assumptions. First, for a $9\times10^6$ M$_\odot$ black hole and 2 possible spins $a^\star=(0,1)$, we fit the X-ray-UV-optical SED with a compact relativistic corona at height $H_x=(46,27)R_g$ irradiating a flat disc with accretion rate $\dot{m}_{Edd}\sim(0.23,0.24)$ inclined to the line of sight by $i<20^\circ$. To fit the lags as well as the SED, this model requires a low spin $a^\star=0$ and boosts disc color temperatures by a factor $f_{col}=1.8$, which shifts reprocessed light to shorter wavelengths. Our Bowl model with $f_{col}=1$ neglects relativity near the black hole but fits the UV-optical lags and SEDs using a flat disc with $\dot{m}_{Edd}<0.06$ and a steep outer rim at $R_{out}/c\sim5-10$ days with H/R<1%. This rim occurs near the $10^3$K dust sublimation temperature in the disc atmosphere, supporting models that invoke dust opacity to thicken the disc and launch failed radiatively-driven dusty outflows at the inner edge of the broad line region (BLR). Finally, the disc lags and SEDs exhibit a significant excess in the $u$ and $r$ bands, suggesting the Balmer continuum and H$\alpha$ emission, respectively, from the BLR.

We report the radio continuum detection of well known Galactic Planetary Nebula (PN) NGC 5189, observed at 943 MHz during the Australian Square Kilometre Array Pathfinder (ASKAP) Evolutionary Map of the Universe (EMU) survey. Two detections of NGC 5189 have been made during the survey, of better resolution than previous radio surveys. Both measurements of the integrated flux density are consistent with each other, at $S_{\rm 943\,MHz} = 0.33\pm0.03$ Jy, and the spectral luminosity is $L_{\rm{943\,MHz}}$ = 8.89 $\times$ 10$^{13}$ W m$^{-2}$ Hz$^{-1}$. Using available flux density measurements for radio-detections of NGC 5189, we calculate a radio surface brightness at 1 GHz and measure $\Sigma_{\rm 1~GHz}$ = 6.0 $\times$ 10$^{-21}$ W m$^{-2}$ Hz$^{-1}$ sr$^{-1}$, which is in the expected range for Galactic PNe. We measure an apparent size of 3.4${'}$ $\times$ 2.2${'}$ corresponding to physical diameters of 1.48 pc $\times$ 0.96 pc, and combine available radio observations of NGC 5189 to estimate a spectral index of $\alpha$ = 0.12 $\pm$ 0.05. Hence, we agree with previous findings that NGC 5189 is a thermal (free-free) emitting nebula. Additional measurements of the optical depth ($\tau = 0.00246$) and electron density ($N_{e} = 138~cm^{-3}$) support our findings that NGC 5189 is optically thin at 943 MHz. Furthermore, the radio contours from the ASKAP-EMU image have been overlaid onto a Hubble Space Telescope (HST) Wide Field Camera 3 image, demonstrating that the radio morphology closely traces the optical. Notably, the contour alignment for the innermost region highlights the two envelopes of gas previously reported to be low-ionisation structures, which is considered a defining feature of post common-envelope PNe that surround a central Wolf-Rayet star.

In light of the latest Planck-Atacama Cosmology Telescope (P-ACT) joint results on the primordial scalar power spectrum, we show that the $R^2$ inflation model extended with a non-minimally coupled scalar field $\chi$--namely the $R^2$-Higgs model--can yield a larger spectral index $n_s$ and a small positive running $\alpha_s$ at cosmic microwave background (CMB) scales, which are consistent with the data. This is because the $\chi$ field contributes a blue-tilted component to the primordial power spectrum which both modifies the large-scale power and, as a result, significantly enhances power on small scales. The consequent enhancement be large enough to lead to the formation of primordial black holes (PBHs) of mass $\lesssim 10^{20}\mathrm{g}$, accounting for all of the dark matter. Furthermore, future observations of the small-scale power spectrum, CMB spectral distortions, and stochastic gravitational waves will provide decisive tests of this model and its predictions for PBHs.

He II $\lambda$1640 emission in galaxies indicates the presence of sources that produce extreme ionizing photons. Here, we report the discovery of a He II $\lambda$1640 emitting galaxy, GNHeII J1236+6215, at $z=$ 2.9803 in the GOODS-north field. We use photometry in 17 wavebands from near-UV to infrared to characterize the galaxy SED and combine Keck LRIS and JWST NIRSpec spectra to identify 15 emission lines including He II $\lambda$1640. We infer that the He$^+$ ionization in the galaxy could be driven by small pockets of young Population III stars or low-metallicity Very Massive Stars (VMSs) rather than AGN or metal-rich Wolf-Rayet stars. The galaxy has a highly ionized ISM ([OIII]5007/[OII]3727 = 7.28$\pm$0.11, [SIII]/[SII] = 1.97$\pm$0.48 and detected Ly$\alpha$, H$\alpha$, H$\beta$, H$\gamma$ lines), little reddening by dust (E(B$-$V) = 0.04$\pm$0.12), low metallicity (12 + log(O/H) = 7.85$\pm$0.22), and high star formation rate (SFR$_{\rm SED}$ = 12.2$\pm$2.0 M$_{\odot}$ yr$^{-1}$). In addition to these ISM conditions, we also notice a significant [SII] deficiency ([SII]6718,6732/H$\alpha$ = 0.08$\pm$0.02, $\Delta$[SII] = $-$0.12) which may indicate the presence of density-bounded optically thin H~II regions that combined with the low dust extinction favor leaking of ionizing Lyman continuum (LyC) photons. Our best-fit SED model also infers a high nebular ionization (log U = $-2.0$) and a low stellar mass M = 7.8$\pm3.1\times$10$^8$M$_{\odot}$. This discovery not only adds one important object to the known sample of high-redshift He~II emitters but also highlights a potential connection between He$^+$ ionization and favorable ISM conditions for the leakage of ionizing photons from galaxies.

The sub-millimetre dust opacity spectral index is a critical observable to constrain dust properties, such as the maximum grain size of an observed dust population. It has been widely measured at galactic scales and down to protoplanetary disks. However, because of observational and analytical challenges, quite a gap exists in measuring dust properties in the envelopes that feed newborn protostars and their disks. To fill this gap, we use sensitive dust continuum emission data at 1.2 and 3.1 mm from the ALMA FAUST Large Program and constrain the dust opacity millimetre spectral index around a sample of protostars. Our high-resolution data, along with a more refined methodology with respect to past efforts, allow us to disentangle disk and envelope contributions in the uv-plane, and thus measure spectral indices for the envelopes uncontaminated by the optically thick emission of the inner regions. First, we find that the young disks are small and optically thick. Secondly, we measure the dust opacity spectral index at envelope scales for n=11 sources: the beta of n=9 sources had never been constrained in the literature. We effectively double the number of sources for which the dust opacity spectral index beta has been measured at these scales. Third, combining the available literature measurements with our own (total n=18), we show how envelope spectral indices distribute between ISM-like and disk-like values, bridging the gap in the inferred dust evolution. Finally, we statistically confirm a significant correlation between beta and the mass of protostellar envelopes, previously suggested in the literature. Our findings indicate that the dust optical properties smoothly vary from the ISM, through envelopes and all the way down to disks. Multi-wavelength surveys are needed to further this study and make more general claims on dust evolution in its pathway from cloud to disks.

We calculate the power spectrum of cosmological perturbations originated from quantum vacuum fluctuations in bouncing scenarios proposed in Ref.~\cite{chamseddine2014cosmology} in the framework of mimetic cosmology. We show that all physically relevant models produce scale invariant spectral indices, and amplitudes compatible with observations provided that the bounce occurs at length scales $t_0$ inside the physically reasonable interval $10^5 l_p < t_0 < 10^9 l_p$. We also show that by slightly modifying the scalar field potential proposed in Ref.~\cite{chamseddine2014cosmology}, we can also obtain the observed red-tilted spectral index, with the same amplitude constraints. Hence, mimetic cosmology provides reasonable bouncing cosmological models without the need of any background quantum effect.

Current and upcoming large-scale structure surveys place stringent requirements on the mitigation of observational systematics in order to achieve their unprecedented constraining power. In this work, we investigate the potential use of the spherical Fourier-Bessel (SFB) power spectrum in controlling systematics, leveraging its capability of disentangling angular and radial scales. We first clarify how the discrete SFB basis describes radial scales via the index $n$ and demonstrate that the SFB power spectrum reduces to the clustering wedge $P(k,\mu)$ in the plane-parallel limit, enabling it to inherit results from past literature based on the clustering wedge. Crucially, the separation of angular and radial scales allows systematics to be localized in SFB space. In particular, systematics with broad and smooth radial distributions primarily concentrate in the $n=0$ modes corresponding to the largest radial scales. This localization behavior enables one to selectively remove only particular angular and radial modes contaminated by systematics. This is in contrast to standard 3D clustering analyses based on power spectrum multipoles, where systematic effects necessitate the removal of all modes below a given $k_{\rm min}$. Our findings advocate for adopting the SFB basis in 3D clustering analyses where systematics have become a limiting factor.

We investigate key methodologies of Bhatnagar-Gross-Krook subgrid modeling for neutrino fast flavor conversions (FFC) in core-collapse supernova based on spherically symmetric Boltzmann radiation hydrodynamics simulations. We first examine time integration methods (explicit, implicit, or semi-implicit) and time step control for the subgrid term, and then compare various approaches in the literature approximating FFCs in two aspects: (1) angular dependent survival probability of neutrinos versus simple equipartition condition with a certain baryon mass density threshold, and (2) 4-species treatment versus 3-species assumption ($\nu_x=\bar\nu_x$). We find that the equipartition condition is reasonable for out-going neutrinos, but large deviations emerge in the incoming neutrinos, that has an influence on matter profiles. We also find that the 3-species model, in which flavor conversions evolve towards erasing electron neutrino lepton number (ELN) crossings, behave differently from the 4-species models where heavy leptonic neutrino number (XLN) are appropriately treated in FFC subgrid modeling. In 4-species models, we commonly observe noticeable differences between $\nu_x$ and $\bar\nu_x$, highlighting the limitation in 3-species treatments to study impacts of flavor conversion on neutrino signals. Our result also suggests that FFC models yield lower neutrino heating rate and smaller shock radii compared to cases with no FFC, in agreement with earlier studies employing quantum kinetic neutrino transport. This work provides valuable information towards robust implementation of FFC subgrid model into classical transport, and serves as a pilot study for future multi-dimensional simulations.

We study the formation of chemical sequences in the stellar disc of Milky Way (MW)-mass galaxies in a full cosmological context with the Auriga simulations. We focus on the conditions giving rise to bi-modal $\alpha$-chemistry in the MW disc and the potential influence of mergers (e.g. Gaia-Enceladus, GSE). We find a wide diversity of chemical sequences, without correlation between the emergence of dichotomous $\alpha$-chemistry and GSE-like mergers. The transition between multiple $\alpha$-sequences is sequential, and is mediated by modulations in the star formation rate (SFR). In some cases, this can be caused by the starburst and subsequent quiescence induced by mergers. In others, by exhaustion or violent disruption of the gas disc. Realisations with singular sequences either lack significant modulations in their SFR, or form too late to have a significant high-$\alpha$ sequence. The metallicity overlap between the high-$\alpha$ and low-$\alpha$ sequences (as seen in the Solar neighbourhood of the MW) arises from accretion of metal-poor gas from the circum-galactic medium. This depends on gas disc thickness, with thinner discs losing their metal-poor extremities. Gas donation from singular gas-rich merger events are incapable of driving long-lived metal dilution ($\Delta\text{[Fe/H]} \gtrsim 0.3$), and we rule-out this scenario for the low-$\alpha$ sequence in the MW. Finally, the shape of $\alpha$-sequences in the [Fe/H] versus [Mg/Fe] plane is related to long-term SFR trends, with a downward slanted locus (as is observed in the low-$\alpha$ sequence of the MW) owing to a sustained or declining SFR.

We present results from the Pan-STARRS optical search for kilonovae without the aid of gravitational wave and gamma-ray burst triggers. The search was conducted from 26 October 2019 to 15 December 2022. During this time, we reported 29,740 transients observed by Pan-STARRS to the IAU Transient Name Server. Of these, 175 were Pan-STARRS credited discoveries that had a host galaxy within 200 Mpc and had discovery absolute magnitudes M > -16.5. A subset of 11 transients was plausibly identified as kilonova candidates by our kilonova prediction algorithm. Through a combination of historical forced photometry, extensive follow-up, and aggregating observations from multiple sky surveys, we eliminated all as kilonova candidates. Rapidly evolving outbursts from massive stars (likely to be Luminous Blue Variable eruptions) accounted for 55% of the subset's contaminating sources. We estimate the rate of such eruptions using the ATLAS 100 Mpc volume-limited survey data. As these outbursts appear to be significant contaminants in kilonova searches, we estimate contaminating numbers when searching gravitational wave skymaps produced by the LIGO-Virgo-Kagra science collaboration during the Rubin era. The Legacy Survey of Space and time, reaching limiting magnitudes of m = 25, could detect 2-6 massive stellar outbursts per 500 deg^2 within a 4-day observing window, within the skymaps and volumes typical for binary neutron star mergers projected for Ligo-Virgo-Kagra Observing run 5. We conclude that while they may be a contaminant, they can be photometrically identified.

Context. Our position within the Galactic plane and the dust obscuration make it challenging to retrieve a true picture of the Milky Way's morphology. While the Milky Way has been recognized as a barred spiral galaxy since the 1960s, there is still uncertainty about the exact number of spiral arms it contains. Currently, our understanding of the Galaxy is evolving thanks to the unprecedented detail provided by Gaia's parallactic distances. Aims. To shed light on the spatial distribution of red supergiants (RSGs) on the Disk and their uniformity of parameters across it, a census of Galactic RSGs detected by Gaia is needed. Methods. Candidate RSGs were extracted from the combined Gaia DR3 and 2MASS catalogs using color criteria and parallactic distances. The sample includes 335 stars that were not included in catalogs of previously known RSGs detected by Gaia DR3. Interstellar and circumstellar extinction values were estimated from the infrared bands. Spectral types were collected from Simbad or VIZIER databases and, for 135 candidates, were inferred from the Gaia DR3 BPRP spectra. Stellar luminosities were inferred using photometric measurements and the Gaia DR3 distances. Results. The analysis yielded a genuine sample of O-rich late-type stars, and the calculated luminosities confirm that the sample is mostly made of stars brighter than Mbol=-5 mag. This new sample represents a 40% increase in the number of highly probable RSGs compared to previous studies. When looking at the X and Y distribution on the Galactic plane, beside the populous Perseus associations of RSGs and the Sagittarius group of RSGs, a novel population of highly probable RSGs populating the more distant Scutum-Centaurus arm appears.

We present a first step toward field-level cosmological inference beyond the standard $\Lambda$CDM model, focusing on optimizing precision tests in the nonlinear regime of large-scale structure (LSS). As an illustrative case, we study the model-independent ``bootstrap'' coefficient of the second-order perturbation theory (PT) kernel for matter in real space, which we use as a proxy for new physics effects in the nonlinear sector. We discuss in details the ultraviolet (UV) cutoff dependence induced by discretizing fields on a grid, which requires proper renormalization to eliminate grid artifacts. We formulate a Wilsonian perturbative framework in which the evolution from a UV theory defined at a high cutoff $\Lambda_\text{uv}$ down to lower cutoffs is computed analytically, even beyond the validity of a derivative expansion. Within this framework, we develop an extended version of the GridSPT code incorporating the bootstrap parameterization and demonstrate how cutoff-independent predictions can be achieved through the inclusion of appropriate counterterms. We validate our approach at third- and fifth-order in PT, emphasizing the importance of higher-derivative contributions for unbiased parameter extraction. Our framework is readily extendable to biased tracers and redshift-space distortions.

We present a Bayesian analysis of cosmological parameter constraints from early- and late-universe observations, focusing on the matter density parameter ($\Omega_m$) and the amplitude of matter fluctuations ($\sigma_8$) within the $\Lambda$CDM framework. Using data from the Planck 2018 satellite mission and the Dark Energy Survey (DES) Year 3, we compute theoretical predictions for angular and matter power spectra via Boltzmann solvers and perform Markov Chain Monte Carlo (MCMC) sampling using the \texttt{emcee} Python package. Our key contribution is a direct and quantitative comparison of DES and Planck constraints, assessing their consistency using chi-squared analysis and Gaussian tension metrics. We find a statistically significant $6.46\sigma$ tension in $\Omega_m$ and a $2.68\sigma$ tension in $\sigma_8$ between the two datasets. These results provide fresh evidence of persistent discrepancies in cosmological parameter estimates and suggest that simple extensions to the $\Lambda$CDM model may be insufficient to fully reconcile early- and late-time observations, motivating the need for more complex theoretical models or refined treatment of systematics.

The evolution of SN 1993J is unlikely to be self-similar. Spatially resolved VLBI-observations show that the velocity of the outer rim of the radio emission region brakes at a few hundred days. The reason for this break remains largely unknown. It is argued here that it is due to the transition between an initial piston phase to a later phase, which is described by the standard model. The properties of the reverse shock are quite different for a piston phase as compared to the standard self-similar model. This affects the expected X-ray emission; for example, the reverse shock becomes transparent to X-ray emission much earlier in the piston phase. Furthermore, it is shown that the observed box-like emission line profiles of H_alpha and other optical lines are consistent with an origin from the transition region between the envelope and the core. It is also pointed out that identifying the observed, simultaneous breaks at approximately 3100 days in the radio and X-ray light curves with the reverse shock reaching the core, makes it possible to directly relate the mass-loss rate of the progenitor star to observables.

The profiles of broad emission lines in active galactic nuclei (AGNs) provide critical insights into the geometry and kinematics of the broad-line region (BLR), which in turn influence the uncertainties in estimating the masses of central supermassive black holes. In this study, we report the discovery of a low-luminosity AGN, SDSS J1609+4902, in which the H$\alpha$ line exhibits two distinct BLR components: a Gaussian component and an extremely broad double-peaked component. Follow-up observations conducted at the Lijiang Observatory in 2025 reveal that the line profile remains roughly unchanged, suggesting that this BLR structure may remain stable over a timescale of $\sim$10 years. We find that the size of the central Gaussian (Full Width at Half Maximum, FWHM$\sim 3000\,{\rm km\, s^{-1}}$) component is consistent with the classical reverberation mapping correlation. In contrast, the asymmetric double-peaked wing (FWHM$\sim 23,000\,{\rm km\, s^{-1}}$) likely originates from a disk-like BLR with an inner radius of approximately 70 gravitational radii. These results provide new constraints on the structure and dynamics of BLRs in AGNs and highlight the potential for long-term stability in such systems.

The local Seyfert 1.5 active galactic nucleus (AGN), NGC 6814, is known to exhibit complex variability, eclipses, and even changing-look behaviour. In this work, we utilize optical-to-X-ray data obtained over 10-years with the Neil Gehrels Swift Observatory to examine the short-term (i.e. daily) and long-term (yearly) variations in the spectral energy distribution (SED). This includes three epochs of high-cadence monitoring (2012, 2016, and 2022), as well as two single observations (2019 and 2021). Model-independent methods of examining the variability suggest that the three monitored epochs exhibit distinct behaviour. X-ray weakness in 2016 can be attributed to the previously studied eclipses, while similar behaviour in 2012 is associated with continuum changes and slight neutral absorption. The multi-epoch SED models are consistent with a black hole (log (MBH / Msun) ~7.6) that is accreting between 0.01-0.1 of the Eddington rate. While the corona (primary X-ray source) is compact, all epochs are better fit with an accretion disk inner radius that is much larger than the innermost stable circular orbit, implying the possibility of a non-standard accretion disk or central structure in NGC 6814.

HD 213985 is an eccentric binary system with a post-AGB primary and a jet-launching secondary star. We confirm that the system photometric variability is likely due to obscuration by the inner edge of a circumbinary disc, similar to RVb-type RV Tau stars. The system has shown an increase in the orbital variability amplitude in optical photometric bands, along with irregular changes in its shape that often started to appear skewed. Variability in the Na D lines suggests that this behaviour may be driven by interactions between the circumbinary disc and outflows through the L2 Lagrange point. Moreover, HD 213985 has exhibited episodes of short-term fluctuations whose appearance is not strictly related to the orbital phase. This variability is consistent with obscuration by transient dust structure leading to weather-like variability patterns.

We present new CO(2-1) observations (resolution $\sim1" = 460$pc) of the Coma cluster jellyfish galaxy NGC 4858 obtained from the ALMA-JELLY large program. Analyzing this data alongside complimentary Subaru H$\alpha$ and HST (F600LP / F350LP) observations, we find numerous structural and kinematic features indicative of the effects from strong, inclined ram pressure, including an asymmetric inner gas tail. We estimate a highly-inclined disk-wind angle of $\phi_{DW} = 75^{+10}_{-27}$. By subtracting a simple circular velocity model, we find (1): gas clumps that are being accelerated by ram pressure, and (2): signatures of gas clumps that had been previously pushed out of the disk but are now falling inwards. We also discuss head-tail morphologies in star complexes within the stellar disk that appear to be RPS-influenced. Lastly, we compare this galaxy to state-of-the-art galaxy ``wind tunnel'' simulations. We find that this galaxy is one of the best nearby examples of strong and inclined ram pressure gas stripping, and of gas that is perturbed by ram pressure but not fully stripped and falls back. We emphasize the importance of torques due to ram pressure in highly-inclined interactions, which help drive gas inwards on the side rotating against the wind, contributing to the formation of asymmetric inner RPS tails.

Modeling observations of transiting exoplanets or close binary systems by comparing the observations with theoretical light curves requires precise knowledge of the distribution of specific intensities across the stellar disk. We aim to facilitate this type of research by providing extensive tabulations of limb-darkening coefficients for 11 frequently used near- and mid-infrared passbands on the NIRCam, NIRISS, and NIRSpec instruments installed on board the James Webb Space Telescope. The calculation of the limb-darkening coefficients was based on spherically symmetric atmosphere models from the PHOENIX series, with high spectral resolution (approximately $10^{6}$ wavelengths), and covering the wavelength range $0.1-6.0~\mu$m. The models were computed for solar composition, and a microturbulent velocity of 1.0 km s$^{-1}$. We adopted two of the more accurate parametrizations for the coefficients: the 4-term law, and the power-2 law. We applied the Levenberg-Marquardt least-squares minimization method, with a strategy to determine the critical value $\mu_{\rm crit}$ of the cosine of the viewing angle near the limb that is designed to improve numerical accuracy. The limb-darkening coefficients were derived based on a total of 306 atmosphere models covering an effective temperature range of $2400-7800$ K, and a $\log g$ interval between 3.0 and 5.5. We discuss the quality of the fits to the specific intensities provided by the power-2 and 4-term laws, as well as by the often used quadratic law. Based on a comparison, we recommend the use of the 4-term or power-2 laws, in that order of preference.

We report on XRISM/Resolve spectroscopy of the recurrent transient and well-known black hole candidate 4U 1630$-$472 during its 2024 outburst. The source was captured at the end of a disk-dominated high/soft state, at an Eddington fraction of $\lambda_\mathrm{Edd} \sim 0.05~(10 M_{\odot}/M_\mathrm{BH})$. A variable absorption spectrum with unprecedented complexity is revealed with the Resolve calorimeter. This marks one of the lowest Eddington fractions at which highly ionized absorption has been detected in an X-ray binary. The strongest lines are fully resolved, with He-like Fe XXV separated into resonance and intercombination components, and H-like Fe XXVI seen as a spin-orbit doublet. The depth of some absorption lines varied by almost an order of magnitude, far more than expected based on a 10% variation in apparent X-ray flux and ionization parameter. The velocity of some absorption components also changed significantly. Jointly modeling two flux segments with a consistent model including four photoionization zones, the spectrum can be described in terms of highly ionized but likely failed winds that sometimes show red-shifts, variable obscuration that may signal asymmetric structures in the middle and outer accretion disk, and a tentative very fast outflow ($v = 0.026-0.033c$). We discuss the impact of these findings on our understanding of accretion and winds in stellar-mass black holes, and potential consequences for future studies.

Both of gyrochronology and the rotation-activity relationship are standard techniques to determine the evolution this http URL mismatch in the definition of the evolutionary phases has so far raised many issues in physics and mathematics and hampered the understanding of how the internal dynamo processes affect the observable this http URL address this problem, we seek a unified scheme that shows a one-to-one mapping from gyrochronology to the rotation--activity this http URL obtained the chromospheric activity $R'_{\rm HK}$ of 6846 stars and their rotation periods to investigate the rotation--activity this http URL apply a three-interval model to fit the relationship and find that it is best fitted by three intervals in the range of Ro$<0.7$.We associate those intervals to the convective, gap and interface phases of this http URL,we suggest an additional epoch at late times of the I this http URL further use the three-interval models to fit the period--activity relationship in temperature bins and determine the duration of the transition phase as a function of effective this http URL comparing the critical temperature and period of the g-to-I transition with the slowly rotating sequence of 10 young open clusters whose ages range from 1 Myr to 2.5 Gyr, we conclude that our new model finds the pure I sequence without fast rotating outliers, which defines the zero-age I sequence (ZAIS).We propose that there is an ambiguous consensus on when the I sequence starts to this http URL ambiguity is from the visually convergent sequence of the color--period diagrams in open this http URL visually convergent sequence is younger than the ZAIS and is actually the pre-I sequence that can be associated with the stall of the spin-down. Our results unify the rotation--activity relationship and gyrochonology for the stellar evolution of low-mass stars, for which we coined the ``CgIW" scenario.

It is assumed that superheavy dark matter particles (SHDM, $\phi$) with $\mathcal{O}$(ZeV) mass may decay to relativistic milli-charged particles (MCPs, $\chi$) via a channel $\phi\to\chi\bar{\chi}$ or $\gamma\chi\bar{\chi}$. The downward-going MCPs passing through the atmosphere can be measured by the Pierre Auger observatory (Auger). The massless hidden photon model is taken for MCPs to interact with nuclei, so that we evaluated the numbers and fluxes of expected MCPs at Auger assuming 14 years of Auger data. The hybrid data of Auger was fitted with the flux of UHE MCPs exceeding GZK-cutoff energy. Then the corresponding upper limits on $\epsilon^2$ are calculated at 90\% C. L.. These results indicate that MCPs can be directly detected with $\epsilon^2\gtrsim 2\times10^{-6}$ at Auger, when $m_{\phi}=2\times10^{21}$ eV and $\tau_{\phi}=10^{27}$ s. And a new region of 10$^{9.6}$ < $m_{MCP}$ < 10$^{11.6}$ eV and $\epsilon$ > $2.12\times10^{-2}$ is ruled out in the $m_{MCP}$-$\epsilon$ plane with 14 years of Auger data. These results indicate potential existence of MCPs and SHDM in the Universe.

We report detection and analysis of the largest low-frequency (200 - 231.25 MHz) sample of Crab giant pulses (GPs) reported in the literature. In total about 95000 GPs were detected. The observations were performed in 2024/2025 with the EDA2, a prototype station of the SKA-Low telescope. The fluence distribution of GPs in the entire sample is very well characterised with a single power law (no flattening at higher fluences) N(F) $\propto$ F$^\alpha$, where $\alpha = -3.17\pm0.02$ for all GPs, and $\alpha_{MP} = -3.13\pm0.02$ and $\alpha_{IP} = -3.59\pm0.06$ for GPs at the phases of the main pulse and interpulse respectively. The index of the power law fluence distribution remained approximately constant over the observing period, but the normalisation of the distribution was strongly correlated with the scatter broadening time ($\tau$). As a result, the measured fluence distribution increased for lower ($\tau \approx$ 2 ms) and decreased for higher ($\tau \approx$ 5 ms) scatter broadening time $\tau$ causing the GP rate to vary between 3000 and 100 per hour respectively. The timescale of variations (weeks) indicates refractive scintillation as the root cause. We also observe a strong positive correlation between the scatter broadening time and dispersion measure. Our modelling favours the screen of the size $\sim10^{-5}$ pc and mean electron density $\sim 400$e$^{-}$cm$^{-3}$ located within 100 pc from the pulsar. The frequency scaling of the scattering broadening time ($\tau \propto \nu^{\beta}$) with $\beta \approx -3.6\pm0.1$ is in agreement with earlier measurements. Our results agree with the current views that GPs from extra-galactic Crab-like pulsars can be responsible for very weak repeating FRBs, but cannot explain the entire FRB population. Finally, these results demonstrate an enormous scientific potential of individual SKA-Low stations.

We present an analysis of the splashback radius ($R_{\text{sp}}$) and the associated splashback mass ($M_{\text{sp}}$) for a sample of galaxy clusters using SDSS spectroscopic data and mock simulations. $R_{\text{sp}}$ marks a physical boundary between the virialized core and the outer infall regions of clusters, providing a robust measure of cluster mass accretion history without being affected by pseudo-evolution. We model the cumulative galaxy number profile of clusters, testing different halo density models and considering the impact of cluster properties, such as center definitions, magnitude limits, galaxy colors, and field contamination, on the estimation of splashback features. Our results show that observed splashback radii, measured in projection (2D), are consistently smaller than predicted by dark matter simulations, with $R_\text{sp}/R_{200m} \approx 1$, supporting previous discrepancies in the literature. We also explore the relationship between $M_{\text{sp}}$ and $R_{\text{sp}}$, proposing a new scaling relation for future cosmological studies, as $R_{\text{sp}}$ is easily observable. Our findings indicate that splashback masses strongly correlate with radii, with a dispersion of $\approx 0.15$ dex, competitive with other mass-observable relations. However, the fitted relation diverges from the constant density expectations of galaxy clusters around $R_\text{sp}$. Additionally, the $M_{\text{sp}} \textendash R_{\text{sp}}$ relation shows significant redshift evolution, though the predominantly low-redshift range of our sample limits our ability to confirm this trend conclusively. The approach developed here may play a key role in cluster characterization and cosmology in the era of large galaxy surveys.

Astrophysical transient phenomena on sub-millisecond timescales, such as fast radio bursts and giant radio pulses from the Crab pulsar, have been primarily observed in radio wavebands. To investigate their origins, a photon detector with high sensitivity and high time resolution is required also in other wavelengths. Recently, we developed the Imager of MPPC-based Optical photoN counter from Yamagata (IMONY), an observation system utilizing a Geiger-mode avalanche photodiode (GAPD) as a sensor. The sensor consists of 64 pixels, each comprising a GAPD and a quenching resistor, with pixel sizes of 75, 100, 150, and 200\,$\mu$m. Each pixel signal is read out independently, enabling single-photon detection. After successfully observing the Crab pulsar using two Japanese telescopes, we upgraded the readout boards to achieve a more compact and stable system. The new system incorporates an analog application-specific integrated circuit (ASIC) developed at KEK for multi-purpose fast readout for silicon photomultipliers. This ASIC features a fast transimpedance amplifier and a comparator, independently processing 16 channels. A Global Navigation Satellite System receiver and a Field Programmable Gate Array (FPGA) provide timestamps for each detected photon with a resolution of 100 ns. The FPGA transmits the acquired data to a PC via Ethernet. This paper presents the details of the new system and the results of its initial evaluation.

Heating mechanism in the solar atmosphere (from chromosphere to corona) is one of the top-challenges in modern astronomy. The classic mechanisms can be divided into two categories: wave heating (W) and magnetic reconnection heating (X). Both of them still face some problems currently difficult to overcome. Recently, we proposed a new mechanism, called magnetic-gradient pumping heating (MGP, or P) which seems to overcome those difficulties, but still lacks sufficient observational evidence. Which one really explained the physics of hot corona exactly? How can observations be used to identify and verify the heating mechanism? Since different heating mechanism will generate non-thermal particles from different accelerations and experience different propagations, they will have different response on the broadband spectral radio observations. Among them, the non-thermal electrons from W mechanisms are closely related to shock-wave acceleration, and their radio response should be group of spike bursts with random distribution of drifting rates; the non-thermal electrons from X mechanisms are accelerated by reconnecting electric field with bidirectional flow, and their radio response should be type III pairs or spike pairs; P mechanism will produce energetic particle upflows, and their radio response should be unidirectional fiber bursts with moderate negative drifting rates. Therefore, the heating mechanism can be identified and verified from the the broadband dynamic spectral radio observations. Additionally, using high-resolution radioheliographs and spectral-imaging observations, the heating mechanisms in different regions can be identified and verified separately, thereby demonstrating the physical essence of hot corona.

Leveraging spectroscopic data from the Sloan Digital Sky Survey, we conduct a comprehensive analysis of the central stellar velocity dispersion, stellar population properties, star formation history, and gas-phase chemical abundances in AGC721966, a unique ultra-diffuse galaxy (UDG) harboring a pseudobulge. Our findings reveal that the pseudobulge formed in the early universe but underwent a recent episode of rejuvenated star formation. The system exhibits a mass-weighted (light-weighted) stellar population age of $\tau_{\star}\sim 7.4\pm2.5$ ($2.9\pm1.5$)~Gyr, a stellar metallicity of [M/H]$\sim -0.62\pm0.26$ ($-0.55\pm0.20$), an $\alpha$-element enhancement of [$\alpha$/Fe]$\sim 0.36\pm0.09$ ($0.37\pm0.07$), and a gas-phase oxygen abundance of \Oabund$\sim 8.15\pm0.03$. The central stellar velocity dispersion is measured as $\sigma_{\rm c}\sim 57.9\pm15.7$~km/s. These results provide robust evidence supporting the early halo-halo merging formation scenario proposed by \cite{Rong25}, while unequivocally ruling out the ``failed'' $L^{\star}$ formation model, at least for AGC721966. Furthermore, through systematic application of the baryonic Tully-Fisher relation, we establish that these pseudobulge-hosting UDGs are neither misidentified nuclear star cluster-bearing dwarf galaxies nor bulge-dominated massive galaxies, thereby affirming their distinct evolutionary pathway.

Waves are an integral part of the solar atmosphere, and their characteristics (e.g., dominant period, range of periods, power, and phase angle) change on a diverse spatio-temporal scale. It is well well-established observationally that the dominant periods of solar oscillations are 5-min and 3-min in the photosphere and chromosphere, respectively. This shows that the wave spectra and their dominant periods evolve between these two layers. We present observational results that demonstrate variations of the dominant period with heights in the photosphere and chromosphere. Six photospheric absorption lines and one chromospheric line are analyzed by using the IRIS data, and the Doppler velocity time series at seven different atmospheric heights are determined. The wavelet analysis is applied to these time series, and the resulting spectrum of wave periods and its dominant period are deduced at these heights, which gives height variations of the dominant period. The obtained data shows that the dominant period decreases with height, and that there are also changes in the range of wave periods within the spectrum. Numerical simulations of filtered wave spectra through the solar atmosphere are also performed, and the obtained results match the observational data.

Stellar systems consisting of three or more stars are not an uncommon occurrence in the Galaxy. Nearly 50% of solar-type wide binaries with separations >1000 au are actually higher-order multiples with one component being a close binary. Additionally, the higher-order multiplicity fraction appears to be correlated with the physical separation of the widest component. These facts have motivated some of our current theories behind how the widest stellar systems formed, which can have separations on the order of or larger than protostellar cores. However, it is unclear if the correlation between wide binary separation and higher-order multiplicity extends to low-mass binaries. We present initial results of an ongoing speckle imaging survey of nearby low-mass wide binaries. We find an overall higher-order multiplicity fraction for our sample of $42.0\% \pm 10.9\%$. If we include systems where Gaia indicates that a companion is likely present, this fraction increases to $62.0\% \pm 14.2\%$. This is consistent with previous results from both higher-mass stars and a previous result for low-mass wide binaries. However, we do not detect the expected increase in higher-order multiplicity fraction with separation, as was seen with previous studies. We briefly explore why higher-order multiplicity statistics could be different in low-mass stars, and what the significance might be for models of wide binary formation.

The gravitational lensing effect of gravitational waves (GWs) has been extensively discussed as a probe of small-mass dark matter halos, which can provide missing information about dark matter. We propose a multi-band observation of lensed GWs from a compact binary to observe both geometrical optics (GO) and wave optics (WO) effects from the same source. This method is expected to be advantageous in breaking parameter degeneracies between a GW source and a dark matter halo acting as a lens. We assume DECIGO or B-DECIGO as a space-based detector observing the early inspiral phase, and the ET as a ground-based detector observing the merger phase. We perform a Fisher analysis of multi-band detection for a source with masses $m_1 = 30 M_{\odot}, m_2 = 20 M_{\odot}$ at redshift $z = 1.5$, and a lens with mass $3 \times 10^{3} M_{\odot}$ at redshift $z = 1.0$. With this setup, the GO effect appears in the ET frequency band, and the WO effect in that of DECIGO. For the halo density profile, we adopt Singular Isothermal Sphere, Cored Isothermal Sphere (CIS), and Navarro-Frenk-White models. We find that multi-band observation resolves parameter degeneracies and significantly reduces errors in estimated parameters. For the CIS model, in particular, we show that, by combining ET and DECIGO observations, the lens mass error improves by about 71 $\%$ and 58 $\%$ compared to ET and DECIGO alone, respectively. Similarly, the impact parameter error is reduced by about 65 $\%$ and 70 $\%$, and the core size error by 34 $\%$ and 68 $\%$, respectively. From these results, we conclude that the multi-band observation of GWs from compact binaries improves the estimation of the lens object properties by breaking the parameter degeneracy.

We use two subsets of 2000 and 1000 \textsc{Quijote} simulations to build two power spectrum emulators, allowing for fast computations of the non-linear matter power spectrum. The first emulator works on scales $k \in [0.015,1.8]\,h/ \rm{Mpc^{-1}}$ in terms of seven cosmological parameters: the matter and baryon fraction of the energy density of the Universe $\Omega_m$ and $\Omega_b$, the reduced Hubble constant $h$, the scalar spectral index $n_s$, the amplitude of matter density fluctuations $\sigma_8$, the total neutrino mass $M_{\nu}$ and the dark energy equation of state parameter $w$. The power spectra can be directly determined at redshifts 0, 0.5, 1, 2 and 3, while for intermediate redshifts these can be interpolated. The second emulator is based on five cosmological parameters, $\Omega_m$, $h$, $n_s$, $\sigma_8$ and the amplitude of equilateral non-Gaussianity $f_{\rm NL}^{\rm equilateral}$, at redshifts 0, 0.503, 0.733, 0.997 for $k \in [0.015,1.8]\,h/ \rm{Mpc^{-1}}$. The emulators are built on machine learning techniques. In both cases we have investigated both neural networks and tree-based methods and we have shown that the best accuracy is obtained for a neural network with two hidden layers. Both emulators achieve a root-mean-squared relative error of less then 5\% for all the redshifts considered on the scales discussed.

Historically, luminous edge-brightened Fanaroff-Riley type II (FRII) radio galaxies have been associated with radiatively efficient cold accretion disks. In contrast, faint edge-darkened Fanaroff-Riley Type Is (FRIs) are sustained by hot, inefficient accretion flows. However, several sources deviate from this paradigm, showing FRII morphologies and radiatively inefficient accretion. Three hypotheses have been proposed to explain the observations: (1) The evolutionary scenario: initially strong accretor FRIIs switch to having a hot inefficient flow due to the depletion of available material. (2) The environment scenario: radio structures are mainly shaped by the environment.\ (3) The engine scenario: both radiatively efficient (FRII-HERGs) and inefficient flows (FRII-LERGs) can expel jets powerful enough to maintain collimation up to kiloparsec distances from the nucleus. We explored these scenarios by studying the jet properties of 3CR radio galaxies using very long baseline interferometry data from the Fundamental Radio Catalog to investigate the distinction between FRII-LERGs and FRII-HERGs. In particular, we analyzed the 8 GHz luminosity ($L_{8GHz}$) distribution across different optical-radio classes at milliarcsecond scales. Our results favor the engine scenario. The separation between different radio-optical classes is evident even on parsec scales. In particular, the separation between FRII classes with different accretion rates occurs near the central engine before the jets reach kiloparsec-scale distances.

To investigate the physical origin of red giants (RGs) and red supergiants (RSGs), we construct steady-state models of stellar envelopes by explicitly solving the stellar structure equations, with boundary conditions set at the stellar surface and the hydrogen burning shell. For comparison, we consider both polytropic and realistic models. Polytropic models adopt a polytropic equation of state (EOS) and neglect energy transport, while realistic models incorporate radiative and convective energy transport, along with tabulated EOS and opacities. Our steady-state solutions reproduce three key features relevant to the evolution toward the RG/RSG phase. First, the refined mirror principle of \cite{Ou2024} is reproduced: the stellar radius varies inversely with the radius of the envelope's inner boundary, defined by the burning shell's surface. This feature arises purely from hydrostatic equilibrium, as it appears in both polytropic and realistic models. Second, realistic models reveal an upper limit to envelope expansion, corresponding to an effective temperature of $\sim4,000\,{\rm K}$, which is characteristic of RG/RSG stars and consistent with the classical Hayashi limit. Approaching this limit, the envelope undergoes a structural transition marked by a significant change in the density profile and the formation of an extended convective zone. The location of this limit is governed by the sharp drop in H$^{-}$ opacity as the envelope cools below $\sim10,000\,{\rm K}$. Finally, our solutions show that even a small shift in the envelope's inner boundary can induce substantial envelope expansion throughout the yellow regime, naturally explaining the bifurcation of giants and supergiants into red and blue branches.

The halo model provides a powerful framework for interpreting galaxy clustering by linking the spatial distribution of dark matter haloes to the underlying matter distribution. A key assumption within the linear bias approximation of the halo model is that the halo bias between two halo populations is a separable function of the mass of each population. In this work, we test the validity of the separable bias approximation on quasi-linear scales using both simulations and observational data. Unlike previous studies, we explicitly disentangle the effects of scale dependence of the bias to assess the robustness of bias separability across a broad range of halo masses and redshifts. In particular, we define a separability function based on halo or galaxy cross-correlations to quantify deviations from halo bias separability, and measure it from N-body simulations. We find significant departures from separability on quasi-linear scales ($\sim$ 1-5 Mpc) at high redshifts ($z \geq 3$), leading to enhancements in halo cross-correlations by up to a factor of 2 -- or even higher. In contrast, deviations at low redshifts remain modest. Additionally, using high-redshift ($z \sim 3.6$) galaxy samples, we detect deviations from bias separability that closely align with simulation predictions. The breakdown of the separable bias approximation on quasi-linear scales at high redshifts underscore the importance to account for non-separability in models of the galaxy-halo connection in this regime. Furthermore, these results highlight the potential of high-redshift galaxy cross-correlations as a probe for improving the galaxy-halo connection from upcoming large-scale surveys.

Blue straggler stars (BSSs) are core hydrogen-burning stars routinely observed in resolved stellar systems. According to the theory of single star evolution they should not exist because, due to their large mass, they should have already evolved to stellar remnants. Hence, they are suspected to be the result of mass-enhancement processes, like mass-transfer activity in binary systems, mergers in binary or higher order systems, or direct stellar collisions, possibly mediated by dynamical interactions. Galactic globular clusters (GGCs) are templates of collisional stellar systems, where frequent gravitational interactions among stars occur on timescales shorter than their age. Thus, in these systems the number of BSSs originated by stellar collisions is expected to increase with the local stellar density and the collision rate. Here we analyse a sample of more than 3000 BSSs homogeneously observed in 48 GGCs with different structures, and we find that the number of BSSs normalized to the sampled luminosity anti-correlates (instead of correlating) with the central density, the collision rate, and the dynamical age of the parent cluster. We also find strikingly similar trends between these environmental properties and the cluster binary fraction, indicating that the latter decreases in high density/high collisional regions. We finally demonstrate that the correlations found between BSSs and the environmental parameters are explained by an underlying dependence of the BSS specific frequency on binary fraction. Once inserted in the context of the current knowledge of the BSS phenomenon, these correlations indicate that low-density regions (possibly because of a higher binary production/survival rate) are the natural habitat of both BSSs and binary systems, and the observed BSSs mostly have a binary-related origin mediated by the environmental conditions.

The main-sequence turnoff (MSTO) stars well preserve the chemical properties where they were born, making them ideal tracers for studying the stellar population. We perform a detailed chemo-dynamical analysis on moderately metal-poor ($-2.0<\mathrm{[Fe/H]}<-1.0$) MSTO stars to explore the early accretion history of the Milky Way. Our sample includes four stars observed with high-resolution spectroscopy using CFHT/ESPaDOnS and 163 nearby MSTO stars selected from the SAGA database with high-resolution results. Within the action-angle spaces, we identified Gaia-Sausage-Enceladus (GSE, 35), stars born in the Milky Way (in situ, 31), and other substructures (21). We find that both GSE and in-situ stars present a similar Li plateau around $A(\mathrm{Li)}\sim 2.17$. GSE shows a clear $\alpha$-knee feature in Mg at $\mathrm{[Fe/H]}\sim-1.60\pm 0.06$, while the $\alpha$-elements of in-situ stars remain nearly constant within the metallicity range. The iron-peak elements show little difference between GSE and in-situ stars except for Zn and Ni, which decrease in GSE at $\mathrm{[Fe/H]}>-1.6$, while they remain constant in in-situ stars. Among heavy elements, GSE shows overall enhancement in Eu, with [Ba/Eu] increasing with the metallicity, while this ratio remains almost constant for in-situ stars, suggesting the contribution of longer time-scale sources to the $s$-process in GSE. Moreover, for the first time, we present the $r$-process abundance pattern for an extremely $r$-process enhanced ($r$-II) GSE star, which appears consistent with the solar $r$-process pattern except for Pr. Further investigation of larger GSE samples using high-resolution spectra is required to explore the reason for the significantly higher Pr in the GSE r-II star.

Galikyan et al. (2025) reported a statistically significant change in galaxy spectral properties at redshift $z$ $\simeq$ 2.7 based on a Kolmogorov analysis of JWST spectroscopic data of galaxies. In this comment, we demonstrate that their result is critically driven by a single outlier in the dataset. This outlier arises from the use of a questionable redshift estimate for one spectrum. When the outlier is removed or the redshift is corrected, the claimed transition at $z$ $\simeq$ 2.7 disappears entirely. By independently reproducing their analysis, we demonstrate that the claimed feature is not a robust statistical signal, but an artefact of this anomalous data point.

Galaxies in the most underdense regions of the Universe, known as cosmic voids, exhibit astrophysical properties that suggest a distinct evolutionary path compared to galaxies in denser environments. Numerical simulations indicate that the assembly of void galaxies occurs later, leading to galaxies with younger stellar populations, low metallicities, and a high gas content in their halos, which provides the fuel to sustain elevated star formation activity. Our objective in this work is to test these numerical predictions with observational data by comparing galaxies in voids with galaxies in non-void environments. We used voids identified in SDSS data and galaxies from the MaNGA survey, which provides galaxies with integral field spectroscopy (IFS). We separated the galaxies into void and non-void samples, mimicked the magnitude distribution, and compared their integrated astrophysical properties as well as the metallicity and age profiles through a stacking technique, ETGs and LTGs separately. We find that void galaxies have younger and less metal-rich stellar populations. Regarding gas mass, we do not find differences across environments. When dividing galaxies into ETGs and LTGs, we observe that ETGs show negative gradients in both age and metallicity, with void galaxies consistently appearing younger and less metal-rich. For LTGs, age gradients are also negative, showing younger populations in void galaxies. However, we do not find statistically significant differences in stellar metallicity gradients between void and non-void environments. Our results show how the astrophysical properties of galaxies in voids differ from those of galaxies in the rest of the Universe. This suggests that the void environment plays a role in the evolution of its galaxies, delaying their assembly and growth.

Homogeneous and nearly-isotropic cosmological models are natural extensions of standard Friedmann cosmologies. Constraining their features is crucial, as any detection of their properties would impact our understanding of inflation and the cosmological principle. Since these models evolve as a set of non-interacting scalar, vector, and tensor modes on top of homogeneous and isotropic spacetimes, their imprints on cosmological observables, particularly the CMB, can be obtained using standard line-of-sight methods. This requires (1) that one resorts on Laplacian eigenmodes on spatially curved spaces and (2) that radial functions for these modes are analytically continued to accommodate complex (i.e., supercurvature) wavenumbers. We introduce two line-of-sight integrators implementing the evolution of the CMB anisotropies in these models: \texttt{AniLoS}, a user-friendly and easy to modify \texttt{Python} package, and \texttt{AniCLASS}, an advanced and efficient extension of the Boltzmann solver \texttt{CLASS}. We discuss possible initial conditions that could generate such fluctuations and provide illustrative examples using our codes. This work offers a pathway for leveraging diverse cosmological datasets to constrain superhorizon anisotropies of the late-time universe.

Deuterium fractionation, as studied through mono-deuterated methanol, is frequently used as a diagnostic tool to trace the physical conditions and chemical evolution of interstellar sources. This study investigates methanol deuteration in the disk around V883 Ori, utilising recent laboratory spectroscopic data for CH$_2$DOH and CH$_3$OD along with ALMA observations. The derived column densities for CH$_2$DOH and CH$_3$OD are (5.14$\pm$0.08) $\times $10$^{16}$ cm$^{-2}$ and (4.22$\pm$0.06) $\times$ 10$^{16}$ cm$^{-2}$, respectively. The analysis demonstrates the influence of spectroscopic data on determining molecular column density, excitation temperature, and, most importantly, the inferred D/H ratio. The D/H ratio for CH$_2$DOH is calculated to be (7.3$\pm$1.5) $\times$ 10$^{-3}$ after applying a statistical correction, whilst the D/H ratio for CH$_3$OD is (1.79$\pm$0.36) $\times$ 10$^{-2}$. The discovery of an unexpectedly low CH$_2$DOH/CH$_3$OD ratio (1.22$\pm$0.02) in V883 Ori, however, raises further questions about the synthesis and chemical processes involved in CH$_3$OD formation. Overall, this study underscores the importance of accurate spectroscopic data for studies of isotopic fractionation and provides new insights into methanol deuteration chemistry in star-forming regions. Future research, combining updated spectroscopy and chemical modelling, will help further constrain these processes across different masses and evolutionary stages.

This work presents comprehensive multi-frequency radio monitoring of the black hole low-mass X-ray binary Swift J1727.8$-$1613, which underwent its first recorded outburst after its discovery in August 2023. Through a considerable community effort, we have coalesced the data from multiple, distinct observing programs; the light curves include ${\sim} 10$ months and 197 epochs of monitoring from 7 radio facilities with observing frequencies ranging from (approximately) 0.3$-$230GHz. The primary purpose of this work is to provide the broader astronomical community with these light curves to assist with the interpretation of other observing campaigns, particularly non-radio observing frequencies. We discuss the phenomenological evolution of the source, which included: (i) multiple radio flares consistent with the launching of discrete jet ejections, the brightest of which reached $\sim$ 1 Jy; (ii) temporally evolving radio spectral indices ($\alpha$), reaching values steeper than expected for optically-thin synchrotron emission ($\alpha {<} -1$) and emission with significant radiative cooling ($\alpha < -1.5$). We have published a digital copy of the data and intend for this work to set a precedent for the community to continue releasing comprehensive radio light curves of future low-mass X-ray binary outbursts.

Spatially resolved polarization measurements of extended X-ray sources are expanding our understanding of the emission mechanisms and magnetic field properties involved. Such measurements have been possible in the past few years thanks to the Imaging X-ray Polarimetry Explorer (IXPE). However, the analysis of extended sources suffers a systematic effect known as polarization leakage, which artificially affects the measured polarization signal. To address this issue, we built a hybrid reconstruction algorithm, which combines machine learning and analytic techniques to improve the reconstruction of photoelectron tracks in the Gas Pixel Detector and to significantly mitigate polarization leakage. This work presents the first application of this hybrid method to experimental data, including both calibration lab measurements and IXPE observational data. We confirmed the reliable performance of the hybrid method for both cases. Additionally, we demonstrated the algorithm's effectiveness in reducing the polarization leakage effect through the analysis of the IXPE observation of the supernova remnant G21.5-0.9. By enabling more reliable polarization measurements, this method can potentially yield deeper insights into the magnetic field structures, particle acceleration processes, and emission mechanisms at work within extended X-ray sources.

The goal of this paper is to describe the science verification of Milky Way Mapper (MWM) ASPCAP data products published in Data Release 19 (DR19) of the fifth phase of the Sloan Digital Sky Survey (SDSS-V). We compare MWM ASPCAP atmospheric parameters T$_{\rm eff}$, log g, 24 abundances of 21 elements (carbon, nitrogen, and oxygen have multiple sources for deriving their abundance values) and their uncertainties determined from APOGEE spectrograph spectra with those of the literature and evaluate their accuracy and precision. We also test the zero-point calibration of the v$_{\rm rad}$ derived by the APOGEE Data Reduction Pipeline. This data release contains ASPCAP parameters for 964,989 stars, including all APOGEE-2 targets expanded with new observations of 336,511 stars from the Apache Point Observatory (APO) observed until 4 July 2023. Overall, the new T$_{\rm eff}$ values show excellent agreement with the IRFM scale, while the surface gravities exhibit slight systematic offsets compared to asteroseisimic gravities. The estimated precision of T$_{\rm eff}$ is between 50 and 70 K for giants and 70$-$100 K for dwarfs, while surface gravities are measured with a precision of 0.07$-$0.09 dex for giants. We achieve an estimated precision of 0.02$-$0.04 dex for multiple elements, including metallicity, $\alpha$, Mg, and Si, while the precision of at least 10 elements is better than 0.1 dex.

The ultraluminous infrared galaxy (ULIRG) IRAS 00183-7111 ($z=0.328$) is one of the three ULIRGs that are currently known to host an active galactic nucleus (AGN) with young radio jets. We present a detailed study of the link between galaxy merger, AGN ignition, radio jet expansion and galactic-scale molecular outflow in IRAS 00183-7111, using high-resolution Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the $^{12}$CO(1-0) and $^{12}$CO(3-2) lines and very deep $i$-band VLT Survey Telescope (VST) imaging. The latter allows us to put - for the first time - tight constraints on the assembly history of the system, confirming that it formed through a major merger between two gas-rich spirals, likely characterised by a prograde encounter and no older than 2 Gyr. The recent merger channelled about $(1.0\pm0.1)\times10^{10}$ M$_{\odot}$ of molecular gas in the central regions of the remnant, as traced by the CO detections. The spatial correlation between the CO distribution and the radio core suggests that this gas likely contributed to the ignition of the AGN and thus to the launch of the radio jets, which are no older than 2 Myr. Furthermore, by comparing the relative strength of the two CO transitions, we find extreme gas excitation (i.e. $T_{\rm ex}\gg50$ K) around the radio lobes, supporting the case for a jet-ISM interaction. A qualitative study of the CO kinematics also demonstrates that, despite the overall disturbed dynamical state with no clear signs of regular rotation, at least one non-rotational kinematic component can be identified and likely associated to an outflow with $v_{\rm out}\approx 556$ km s$^{-1}$ and $\dot{M_{\rm out}}\approx 110$ M$_{\odot}$ yr$^{-1}$.

Observations of stellar occultations have revealed that small bodies are capable of hosting ring systems. The trans-Neptunian object (TNO) Quaoar, is the host of an enigmatic ring system, with two rings located well-outside the Roche limit. To better understand these structures, we observed a stellar occultation by Quaoar and its rings using the James Webb Space Telescope's (JWST) NIRCam instrument. Our observations detect both known rings, although Q2R -- the inner known ring -- is not detected on both sides of Quaoar, showing that it has substantial azimuthal variations similar to Q1R -- the outer ring. We also fit a model of the ring radii and pole orientation of the ring system, which confirms that Quaoar's spin-orbit and Weywot's mean motion resonances (especially Weywot's 6:1) may play a role in the rings' confinement and stability. In addition to examination of Quaoar's ring system, we also use our observations to place upper limits on a putative CH$_4$ atmosphere around Quaoar, finding that no global atmosphere with surface pressure $>1$ nbar can exist (at 3$\sigma$ significance). The lack of atmosphere supports the hypothesis that atmospheric processes are not the source of Quaoar's recently discovered inventory of light hydrocarbons.

Stellar bowshock nebulae are arcuate shock fronts formed by the interaction of radiation-driven stellar winds and the relative motion of the ambient interstellar material. Stellar bowshock nebulae provide a promising means to measure wind-driven mass loss, independent of other established methods. In this work, we characterize the stellar sources at the center of bowshock nebulae drawn from all-sky catalogs of 24 $\mu$m-selected nebulae. We obtain new, low-resolution blue optical spectra for \numstars~stars and measure stellar parameters temperature \teff, surface gravity \logg, and projected rotational broadening \vsini. We perform additional photometric analysis to measure stellar radius \rstar, luminosity \lstar, and visual-band extinction $A_V$. All but one of our targets are O and early B stars, with temperatures ranging from $T$=16.5--46.8~k\kelvin, gravities $\log g=$2.57--4.60, and \vsini~from $<$100--400~\kms. With the exception of rapid rotator $\zeta$ Oph, bowshock stars do not rotate at or near critical velocities. At least 60 of 103 (60\%) OB bowshock stars are binaries, consistent with the multiplicity fraction of other OB samples. The sample shows a runaway fraction of 23\%, with \numrunaway~stars having $v_{\text{2D}}\geq25$~\kms. Of the 19 runaways, at least 15 ($\geq$79\%) are binaries, favoring dynamical ejection over the binary supernova channel for producing runaways. We provide a comprehensive census of stellar parameters for bowshock stars, useful as a foundation for determining the mass-loss rates for OB-type stars -- one of the single most critical factors in stellar evolution governing the production of neutron stars and black holes.

We aim to detect low-frequency radio emission from exoplanetary systems, which can provide insights into planetary magnetic fields, star-planet interactions, stellar activity, and exo-space weather. The HD 189733 system, hosting a well-studied hot Jupiter, is a prime target for such searches. We conducted NenuFAR imaging observations in the 15-62 MHz range, in order to cover the entire orbital phase of HD 189733 b. Dynamic spectra were generated for the target and other sources in the field, followed by a transient search in the time-frequency plane. The data processing pipeline incorporated direction-dependent calibration and noise characterization to improve sensitivity. We also searched for periodic signals using Lomb-Scargle analysis. A highly circularly polarized radio burst was detected at 50 MHz with a flux density of 1.5 Jy and a significance of 6 sigma at the position of HD 189733. No counterpart was found in Stokes I, likely because the emission is embedded in confusion noise and remains below the detection threshold. The estimated minimum fractional circular polarization of 38% suggests a coherent emission process. A periodicity search revealed no weaker signals linked to the planet's orbital period, the star's rotational period, or the synodic period and harmonic period between them. The burst's properties are consistent with cyclotron maser instability (CMI) emission, but the origin is still ambiguous. The comparison with theoretical models suggests star-planet interaction or stellar activity as potential origins. However, alternative explanations such as contamination from other sources along the line of sight (e.g. the companion M dwarf) or noise fluctuation cannot be ruled out.

For three decades, adaptive optic surveys have revealed an excess of T Tauri binaries across a = 10-100 au in nearby star-forming regions compared to the field population of main-sequence (MS) stars. Such an excess requires that most stars are born in dense clusters and subjected to significant dynamical processing that disrupts such binaries across intermediate separations. However, we demonstrate that the apparent excess is due to an observational selection bias. Close binaries within a < 100 au clear out their dusty circumstellar disks on faster timescales compared to wide binaries and single stars. A magnitude-limited sample is therefore biased toward close binaries that have preferentially cleared out their obscuring disks. We re-examine the separation distribution of pre-MS binaries in low-density Taurus, moderately dense Upper Scorpius, and the extremely dense Orion Nebula Cluster (ONC). By limiting the samples to primary spectral type / mass instead of magnitude, the artificial excess across a = 10-100 au disappears in all three environments. Across wider separations a = 100-4,000 au, Taurus exhibits an excess of companions (mostly tertiaries), the ONC displays a deficit, and Upper Scorpius matches the field MS population. The field derives from an amalgam of all three environments, where Upper Scorpius corresponds to the average birth environment of solar-type stars. The total binary fraction within a < 10,000 au in Taurus is only 52% +/- 7%, substantially lower than the 100% inferred from the biased observations and only slightly higher than the field MS value of 45%. N-body interactions preferentially disrupt outer tertiaries with only marginal dynamical processing of the inner binaries, especially those within a < 100 au.

Previous studies of the strong lens system SDSSJ0946+1006 have reported a dark matter subhalo with an unusually high central density, potentially challenging the standard cold dark matter (CDM) paradigm. However, these analyses assumed the subhalo to be completely dark, neglecting the possibility that it may host a faint galaxy. In this work, we revisit the lensing analysis of SDSSJ0946+1006, explicitly modelling the subhalo as a luminous satellite. Incorporating light from the perturber broadens the range of allowed subhalo properties, revealing solutions with significantly lower central densities that are consistent with CDM expectations. The inferred luminosity of the satellite also aligns with predictions from hydrodynamical simulations. While high-concentration subhaloes remain allowed, they are no longer statistically preferred. The luminous subhalo model yields a better fit to the data, while also offering a more plausible explanation that is in line with theoretical expectations. We validate our methodology using mock data, demonstrating that neglecting subhalo light can lead to inferred mass distributions that are artificially compact.

This paper investigates resonant neutrino self-interactions in cosmology by employing, for the first time, the effective field theory of large-scale structure to model their impact on the matter distribution up to mildly nonlinear scales. We explore a broad range of mediator masses in two main analyses: one combining BOSS Full Shape (FS) galaxy clustering with Big Bang Nucleosynthesis (BBN), and another combining FS with Planck Cosmic Microwave Background (CMB) data. Our results place the strongest cosmological constraints to date on the resonant self interactions when using FS+Planck data, reaching up to $g_\nu< 10.8 \times 10^{-14}$ at $95\%$ confidence for the 1 eV mediator. Notably, FS+BBN can constrain the interaction for the 10 eV mediator independently of CMB data, yielding $g_\nu < 7.33 \times 10^{-12}$ at $95\%$ confidence. Our results suggest that resonant neutrino self-interactions are unlikely to resolve existing cosmological tensions within the standard $\Lambda$CDM framework.

M-type stars are prime targets for exoplanet searches within their habitable zones (HZs). M-type stars also exhibit significant magnetic flaring activity, especially during the first billion years, which can potentially accelerate the evaporation of the hydrogen-helium envelopes of close-in planets. We use the time-dependent photo-ionization hydrodynamics code ATES to investigate the impact of flares on atmospheric escape, focusing on a sub-Neptune-sized planet orbiting an early M-type star at distances of 0.01, 0.1, and 0.18-0.36 AU--the inner and outer edges of the HZ. Stellar flaring is modeled as a 1 Gyr-long high-activity phase followed by a 4 Gyr-long low-activity phase, each characterized by a suitable flare frequency distribution. We find that flares have a modest impact (i.e., less than a factor of two) on increasing the cumulative atmospheric mass loss, with the greatest enhancement occurring when the planet is at its closest separation. However, the relative enhancement in mass loss between flaring and non-flaring cases is greater at larger separations. This trend arises because, as stellar irradiation fluctuates between quiescent levels and peak flares, the proportion of time the planet spends in energy-limited versus recombination-limited mass loss regimes depends on its orbital separation. Additionally, we demonstrate the existence of a characteristic flare energy, intermediate between the minimum and maximum values, that maximizes the fractional contribution to flare-driven mass loss. Our results indicate that the flaring activity of M-dwarfs does not significantly affect the atmospheric retention of close-in Neptune-sized planets, including within the HZ. The potential occurrence of rare super-flares, which observational campaigns might be biased against, does not alter our conclusions.

This paper explores the phenomenon of particle creation associated with cosmic strings in de Sitter spacetime, a model that represents the universe exponential expansion. We examine how the presence of cosmic strings in a de Sitter background affects particle production, focusing on the roles of string tension and angular deficits. Utilizing the Klein Gordon equation adapted to curved spacetime with cosmic string defects, we derive solutions expressed through hypergeometric functions to describe particle states. Our findings highlight how string properties influence particle creation rates and energy distributions. By analyzing both point-like and linear potentials near the string, we determine exact solutions, investigate asymptotic behaviors, and calculate particle creation probabilities using Bogoliubov transformations.

Neutrinos from blazars can originate from inelastic scatterings between protons within their jets and sub-GeV dark matter (DM) around them, explaining IceCube detections of neutrinos from TXS 0506+056 that are otherwise challenging for models of its jet. In this paper we calculate such DM-induced high-energy neutrinos, from TXS 0506+056 as well as from a stacked blazar sample, in the four cases where DM-quark interactions are mediated by a new massive vector, axial, scalar, and pseudoscalar particle. Intriguingly, we find that this mechanism can saturate the diffuse astrophysical neutrino flux observed by IceCube at high energies. Our mechanism will be tested by additional blazar observations and by various searches for sub-GeV DM.

We demonstrate that the scalaron, a scalar degree of freedom, emerging from the $f(R)$ theory of gravity, can account for the observed dark matter (DM) abundance if its mass is around the MeV scale, to ensure its cosmological stability. Focusing on two well-known $f(R)$ gravity models, we systematically show that if scalaron production proceeds via the freeze-in mechanism, the right relic abundance is satisfied over a very narrow window of reheating temperature $10^{14}\lesssim T_{\rm rh}\lesssim 10^{16}$ GeV. We delineate the viable parameter space of the $f(R)$ models consistent with the observed DM abundance, and highlight relevant experimental constraints from searches targeting DM decay signatures.

In this work, we present new methods implemented in the GstLAL offline gravitational wave search. These include a technique to reuse the matched filtering data products from a GstLAL online analysis, which hugely reduces the time and computational resources required to obtain offline results; a technique to combine these results with a separate search for heavier black hole mergers, enabling detections from a larger set of gravitational wave sources; changes to the likelihood ratio which increases the sensitivity of the analysis; and two separate changes to the background estimation, allowing more precise significance estimation of gravitational wave candidates. Some of these methods increase the sensitivity of the analysis, whereas others correct previous mis-estimations of sensitivity by eliminating false positives. These methods have been adopted for GstLAL's offline results during the fourth observing run of LIGO, Virgo, and KAGRA (O4). To test these new methods, we perform an offline analysis over one chunk of O3 data, lasting from May 12 19:36:42 UTC 2019 to May 21 14:45:08 UTC 2019, and compare it with previous GstLAL results over the same period of time. We show that cumulatively these methods afford around a 50% - 100% increase in sensitivity in the highest mass space, while simultaneously increasing the reliability of results, and making them more reusable and computationally cheaper.

We model a single black hole in equilibrium with a dark photon-cold dark matter environment. Representing the dark photon as a Proca field, we show that a Schwarzschild black hole grows vector-field "hair" when allowed to accrete from an infinite homogeneous bath of particles far from the horizon. We solve the Proca equation in linear perturbation theory, separating it using the vector spherical harmonics and Frolov-Krtouš-Kubizňák-Santos approaches for the odd-parity and even-parity sectors, respectively. In the "particle" dark matter regime, the field is purely infalling and exhibits a sharply peaked density profile, in concordance with the particle dark matter "spikes" studied in the literature. In the "wave" regime, the field exhibits standing waves, and the profile is smeared. We find a dark-matter density amplification upward of $10^7$ near the horizon. Though small for most black holes, we find the mass enclosed in the cloud can reach $\sim 1 \%$ of the black hole mass for large supermassive black holes. These black holes are also most susceptible to vector dark matter accretion, with mass accretion rates as large as $10 M_\odot/$yr.

We use a simple and straightforward method to derive the energy partition in magnetohydrodynamics (MHD) turbulence that was first studied by Lee and then more rigorously by Chandrasekhar. By investigating the energy equation we find that the turbulent viscous and ohmic dissipations are comparable to each other. Under the condition that turbulent viscosity and turbulent magnetic diffusivity are comparable, we deduce that the ratio of kinetic to magnetic energies depends on the ratio of the turbulent magnetic lengthscale to turbulent velocity lengthscale of the largest eddies. When the two largest lengthscales are comparable, the two energies are in equipartition.

We investigate the production of primordial black holes and gravitational waves in composite hybrid inflation. Starting from an effective chiral Lagrangian with a dilaton and pions, we identify inflation occurring due to the walking dynamics of the theory. A $\mathbb{Z}_2$ symmetry-breaking term in the pion sector induces a shift in the inflaton's trajectory, which leads to a tachyonic instability phase. Curvature perturbations grow exponentially, producing copious primordial black holes and a stochastic gravitational wave background. We show that the primordial black hole mass and the gravitational wave frequency are strongly restricted by the anomalous dimensions of the pion operators, with larger anomalous dimensions giving lighter primordial black holes and higher frequency gravitational waves. In both cases, the associated signatures lie within reach of future gravitational wave observatories.

We present a study of LLM integration in final-year undergraduate astronomy education, examining how students develop AI literacy through structured guidance and documentation requirements. We developed AstroTutor, a domain-specific astronomy tutoring system enhanced with curated arXiv content, and deployed it alongside general-purpose LLMs in the course. Students documented their AI usage through homework reflections and post-course surveys. We analyzed student evolution in AI interaction strategies and conducted experimental comparisons of LLM-assisted versus traditional grading methods. LLM grading showed strong correlation with human evaluation while providing more detailed and consistent feedback. We also piloted LLM-facilitated interview-based examinations as a scalable alternative to traditional assessments, demonstrating potential for individualized evaluation that addresses common testing limitations. Students experienced decreased rather than increased reliance on LLMs over the semester, developing critical evaluation skills and strategic tool selection. They evolved from basic assistance-seeking to verification workflows, with documentation requirements fostering metacognitive awareness. Students developed effective prompting strategies, contextual enrichment techniques, and cross-verification practices. Our findings suggest that structured LLM integration with transparency requirements and domain-specific tools can enhance astronomy education while building essential AI literacy skills. We provide implementation guidelines for educators and make our AstroTutor repository freely available.

We propose a novel method to detect axions which uniquely depends on the dimensionless axion-photon coupling $g_{\gamma}$, independent of the suppressive axion decay constant $f_a$. Using neutron star PSR B1919+21 data from the Five-hundred-meter Aperture Spherical Telescope, we derive the first constraint $|g_{\gamma}|<0.93$ at $1\sigma$ confidence level for ultra-light axions ($m_a < 10^{-11}$ eV).

In the nanohertz band, the spatial correlations in pulsar timing arrays (PTAs) produced by interfering gravitational waves (GWs) from multiple sources likely deviate from the traditional ones without interference under the assumption of an isotropic Gaussian ensemble. This work investigates the impact of such interference within the framework of massive gravity. Through simulations, we show that while the resulting correlation patterns can be described by Legendre expansions with coefficients that depend on the interference configuration, they remain predominantly quadrupolar (l = 2), with this feature becoming more pronounced as the graviton mass increases--reflecting both the tensorial polarizations and the modified GW dispersion. However, the interference introduces significant variability in the angular correlation, making it difficult to distinguish massive gravity from general relativity based on a single realization of the Universe. We conclude that beyond a fundamental constraint set by the PTA observation time, achieving a substantially tighter bound on the graviton mass is statistically challenging and observationally limited under realistic conditions.

Collisionless shocks are frequently analyzed using the magnetohydrodynamic formalism (MHD), even though the required collisionality hypothesis is not fulfilled. In a previous work \citep{BretJPP2018}, we presented a model of collisionless shock displaying an important departure from the expected MHD behavior, in the case of a strong flow aligned magnetic field. This model was non-relativistic. Here, it is extended to the relativistic regime, considering zero upstream pressure and upstream Lorentz factor $\gg 1$. The result agrees satisfactorily with Particle-in-Cell simulations and shows a similar, and important, departure from the MHD prediction. In the strong field regime, the density jump $r$, seen in the downstream frame, behaves like $r \sim 2 + 1/\gamma_{\mathrm{up}}$ while MHD predicts 4 ($\gamma_{\mathrm{up}}$ is the Lorentz factor of the upstream measured in the downstream frame). Only pair plasmas are considered.