Papers for Thursday, Jan 23 2025

Gaia data have revealed vertically asymmetric phase-space structures in the Milky Way (MW) disc, such as phase spirals, indicating vertical oscillations. These oscillations exhibit two distinct modes: the bending mode and the breathing mode, associated with one-arm and two-arm phase spirals, respectively. This study aims to explore the excitation mechanisms of the bending and breathing modes and their subsequent evolution in the MW disc, focusing on the interplay between direct perturbations from the Sagittarius dwarf galaxy and indirect contributions from tidally induced spiral arms. We perform high-resolution $N$-body simulations to model the interaction between an MW-like disc galaxy and a Sagittarius dwarf-like satellite. These simulations resolve fine phase-space structures, enabling analysis of the bending and breathing modes at both macroscopic (global bending and breathing waves) and microscopic (local phase spirals) scales. Our simulations demonstrate that the satellite's perturbation directly excites the bending mode and induces spiral arms in the galactic disc. These spiral arms excite the breathing mode, making it an indirect consequence of the satellite interaction. Initially, the bending mode dominates, but it rapidly decays due to horizontal mixing. In contrast, the breathing mode persists for a longer duration, sustained by the spiral arms, leading to a transition from a bending-dominated to a breathing-dominated state. This transition progresses faster in the inner galaxy than in the outer regions. The simulations reproduce the one-arm phase spiral observed in the solar neighbourhood and reveal two-arm phase spirals, particularly in the inner galaxy, associated with spiral arm-induced breathing modes. Our findings highlight the combined effects of direct satellite perturbations and indirect spiral arm dynamics in shaping the vertical structure of the MW disc.

We introduce the AIDA-TNG project, a suite of cosmological magnetohydrodynamic simulations that simultaneously model galaxy formation and different variations of the underlying dark matter model. We consider the standard cold dark matter model and five variations, including three warm dark matter scenarios and two self-interacting models with constant or velocity-dependent cross-section. In each model, we simulate two cosmological boxes of 51.7 and 110.7 Mpc on a side, with the same initial conditions as TNG50 and TNG100, and combine the variations in the physics of dark matter with the fiducial IllustrisTNG galaxy formation model. The AIDA-TNG runs are thus ideal for studying the simultaneous effect of baryons and alternative dark matter models on observable properties of galaxies and large-scale structures. We resolve haloes in the range between $10^{8}$ and $4\times10^{14}\,$M$_{\odot}$ and scales down to the nominal resolution of 570 pc in the highest resolution runs. This work presents the first results on statistical quantities such as the halo mass function and the matter power spectrum; we quantify the modification in the number of haloes and the power on scales smaller than 1 Mpc, due to the combination of baryonic and dark matter physics. Despite being calibrated on cold dark matter, we find that the TNG galaxy formation model can produce a realistic galaxy population in all scenarios. The stellar and gas mass fraction, stellar mass function, black hole mass as a function of stellar mass and star formation rate density are very similar in all dark matter models, with some deviations only in the most extreme warm dark matter model. Finally, we also quantify changes in halo structure due to warm and self-interacting dark matter, which appear in the density profiles, concentration-mass relation and galaxy sizes.

X-ray observations of the optically selected $z=6.025$ QSO CFHQS J164121+375520 (hereafter J1641) revealed that its flux dropped by a factor $\gtrsim7$ from 2018, when it was a bright and soft X-ray source, to 2021. Such a strong variability amplitude has not been observed before among $z>6$ QSOs, and the underlying physical mechanism was unclear. We carried out a new X-ray and rest-frame UV monitoring campaign of J1641 over 2022-2024. We detected J1641 with Chandra in the 2-7 keV band, while no significant emission is detected at softer X-ray energies, making J1641 an X-ray changing look QSO at $z>6$. Comparing with the 2018 epoch, the 0.5-2 keV flux dropped dramatically by a factor $>20$. We ascribe this behaviour to intervening, and still ongoing, obscuration by Compton-thick gas intercepting our line of sight between 2018 and 2021. The screening material could be an inner disk or a failed nuclear wind that increased their thickness. Another possibility is that we have witnessed an occultation event due to dust-free clouds located at sub-pc/pc scales, similar to those recently invoked to explain the remarkable X-ray weakness of AGN discovered by JWST. These interpretations are also consistent with the lack of strong variations of the QSO rest-frame UV lightcurve over the same period. Future monitoring of J1641 and the possible discovery of other X-ray changing look QSOs at $z>6$ will provide us with precious information about the physics of rapid supermassive black-hole growth at high redshift.

The observations of eclipsing binary systems are of great importance in astrophysics, as they allow direct measurements of fundamental stellar parameters. By analysing high-quality space-based observations with ground-based photometric data, it becomes possible to detect these fundamental parameters with greater precision using multicolour photometry. Here, we report the first photometric analysis results of the V517 Cam eclipsing binary system by combining the Transiting Exoplanet Survey Satellite (TESS) light curve and new CCD observations in BVRI filters, obtained with a 60 cm robotic telescope (T60) at the TÜBİTAK National Observatory. By means of photometric analyses, the masses and radii of the primary and secondary stars were carefully determined to be $M_{1}= 1.47\pm 0.06\,M_\odot$, $M_{2}= 0.79\pm0.05\,M_\odot$, and $R_{1}=1.43\pm 0.03\,R_\odot$, $R_{2}= 0.75\pm 0.04\,R_\odot$, respectively. Furthermore, the distance to V517 Cam was calculated to be $284\pm20$ pc. The overall age of the system is estimated to be around $63\pm15$ Myr. At this age, the primary component stands near the onset of its main-sequence evolution, near the ZAMS, whereas the secondary component remains in the pre-main-sequence evolutionary phase. To better understand the evolutionary status and nature of V517 Cam, the mass ratio and temperature values, obtained with relatively low sensitivity by photometric measurements, need to be confirmed by spectral analysis.

Wide Field Survey Telescope (WFST), with a powerful sky survey capability in the northern hemisphere, will play an important role in asteroid searching and monitoring. However, WFST is not a telescope dedicated to near-Earth asteroids (NEOs) searching. In order to improve the efficiency of finding NEOs on the premise of meeting the needs of other scientific research, we ran mock observations for WFST to study its search capability for NEOs. The NEO population model, the WFST detection model and site conditions are taken into account in our simulations. Based on the original scheduling scheme, we present two new schemes. Compared to the original scheme, the optimized scheme can improve the search capability of known and unknown NEOs by 100% and 50%. We also emphasized the importance of trailing loss and proposed an improved effective field of view model. In addition, it is predicted that adopting the clear-day ratio of 0.7 and the optimized scheme, for NEOs with absolute magnitude from 17 to 25, WFST can provide tracklets for about 1800 NEOs if their orbits are known, and in the case of blind search, more than 600 NEOs can be found by WFST. The new schemes provide valuable reference and suggestions for the WFST's regular survey strategy.

We present time-resolved Gran Telescopio Canarias optical spectroscopy and William Herschel Telescope $i$-band photometry of the X-ray transient SWIFT J1753.5-0127 in quiescence. The $i$-band light curve is dominated by flickering with an amplitude of $\sim 0.5$ mag and shows no evidence of the ellipsoidal modulation of the companion star. The telluric-corrected average spectrum, on the other hand, reveals the presence of weak (strongly veiled) TiO bands at $7055$ Ȧ and $7589$ Ȧ. We used them for a spectral classification, finding an M4-5 V companion star. However, as velocity shifts are not clearly detected in the individual spectra, we turned the analysis to the double-peaked H$\alpha$ emission line from the accretion disc. By exploiting the empirical correlations established for quiescent X-ray transients between the line morphology and fundamental binary parameters, we estimated the radial velocity semi-amplitude of the companion $K_2 = 820 \pm 36$ km s$^{-1}$, a mass ratio $q = 0.023 \pm 0.006$ and an inclination $i = 79 \pm 5$ deg. Moreover, an orbital period of $3.26 \pm 0.02$ h was measured from the modulation of the centroid velocities and the double-peak trough depth of the H$\alpha$ profile. These quantities yielded a mass function $f(M_1) = 7.8 \pm 1.0$ M$_\odot$ and black hole and companion star masses of $M_1 = 8.8 \pm 1.3$ M$_\odot$ and $M_2 = 0.20 \pm 0.06$ M$_\odot$, respectively. The companion star mass is in line with the spectral classification obtained from the relative depth of the TiO bands. Based on the mean quiescent magnitude ($i = 21.4 \pm 0.1$), orbital period, and interstellar extinction, we estimate the distance to the source to be $3.9 \pm 0.7$ kpc and a Galactic plane elevation of $0.8 \pm 0.2$ kpc, supporting the case for a large natal kick.

This document specifies the structure of the Siril Catalog Format, designed to support efficient storage and querying of nested HEALpixel based astronomical catalogs such as Gaia DR3. The format includes a fixed length header, an index for rapid look-up, and structured data records optimized for space and speed. The document also outlines recommended search strategies for utilizing the format effectively.

The kinematic parameters identified from high-velocity stars situated within 100 kpc are examined and analyzed. We included three high velocity programs comprising 591, 87, and 519 stars as a function of distances ranging from 0.10 to nearly 109 kpc. In this analysis, we will determine the spatial velocities (U, V, W) in galactic coordinates along with their velocity dispersion (sigma_1, sigma_2, sigma_3), the convergent point (Ao, Do), and therefore, the solar motion (S_sun).

The circumstellar envelope of the carbon star CW Leo exhibited various unexpected changes in recent optical imaging observations. We have performed a follow-up observation using the Near-infrared Integral-Field Spectrograph (NIFS) equipped on the Gemini-North telescope. We report the near-infrared counterparts of a local brightness peak in the optical at the stellar position of CW Leo. On the other hand, a second peak detected at short wavelengths in the J band coincides with the brightest, bluest position in the optical images. The absorption features in the K band are minimized at a radius of 0.2 arcsec from the predicted stellar position. The reduction of the absorption depths likely indicates dilution of the absorption features by thermal emission of dust grains newly formed at such a radius and heated by radiation from the central star. The broad absorption feature at 1.53 um is significantly deeper than in template carbon stars, consistent with the presence of a substantial amount of circumstellar material around CW Leo. Its northeastern quadrant lacks circumstellar absorption features and scattered light in the near-infrared regime, which are possibly manifestations of its conical cavity in both gas and dust. In addition, a cross correlation of CO overtone bands indicates that the average expansion velocity of dust grains is smaller to the northern direction, likewise the velocity of transverse wind components derived using the differential proper motion of a circumstellar whirled pattern. The gradual brightening of CW Leo and the changes in its innermost circumstellar envelope need further continuous monitoring observations to properly understand its transitional phase toward the post-asymptotic-giant-branch stage.

Strong lensing of gravitational wave (GW) sources allows the observer to see the GW source from different lines-of-sight (LOS) through the corresponding images, which provides a way for constraining the relative proper motion of the GW source. This is possible as the GW signals received from each image will have slightly different projected velocity components, from which one can `Doppler-Triangulate' for the GW source velocity vector. The difference in projected velocity between the different images can be observationally inferred through pairwise GW phase measurements that accumulate over the time-of-observation. In this paper we study lensed eccentric GW sources and explore how the observable GW phase shift between images evolve as a function of time, eccentricity, lens- and binary parameters. Next generation GW observatories, including the Einstein Telescope and Cosmic Explorer, will see $\sim $hundreds/year of lensed GW sources, where a significant fraction of these are expected to be eccentric. We discuss the expected unique observables for such eccentric lensed GW sources, and the relation to their observable relative linear motion, which otherwise is exceedingly difficult to constrain in general.

In this work, we propose a novel ensemble of recurrent neural networks (RNNs) that considers the multiband and non-uniform cadence without having to compute complex features. Our proposed model consists of an ensemble of RNNs, which do not require the entire light curve to perform inference, making the inference process simpler. The ensemble is able to adapt to varying numbers of bands, tested on three real light curve datasets, namely Gaia, Pan-STARRS1, and ZTF, to demonstrate its potential for generalization. We also show the capabilities of deep learning to perform not only classification, but also regression of physical parameters such as effective temperature and radius. Our ensemble model demonstrates superior performance in scenarios with fewer observations, thus providing potential for early classification of sources from facilities such as Vera C. Rubin Observatory's LSST. The results underline the model's effectiveness and flexibility, making it a promising tool for future astronomical surveys. Our research has shown that a multitask learning approach can enrich the embeddings obtained by the models, making them instrumental to solve additional tasks, such as determining the orbital parameters of binary systems or estimating parameters for object types beyond periodic ones.

Tidally induced spin-up of stripped helium stars in short-period (<\,1\,d) binaries with black holes has been presented as one of the possible mechanisms to reproduce the high-spin tail of the black hole spin distribution derived from gravitational wave (GW) merger observations. At such short periods, a fraction of the strong stellar wind from the stripped helium stars may be accreted by the black holes, and its gravitational potential energy may be released as observable radiation in the X-ray regime. We estimate the X-ray luminosity and its observability from black holes in orbit with stripped helium stars, which evolve into binary black hole or black hole+neutron star binaries that merge within Hubble time. We post-process recent advancements for estimating X-ray luminosities (via wind accretion onto stellar mass black holes) into two rapid population synthesis codes, BSE and StarTrack. We derive lower limits on the X-ray luminosity distribution from populations of stripped helium+black hole binaries at four metallicities (0.01, 0.1, 0.5, 1 $Z_{\odot}$) and two mass transfer stability criteria. We find that a large fraction (0.1-0.5) of stripped-helium stars that get spun up by tides also transfer enough wind matter onto the black hole to produce X-ray luminosities above $10^{35}$\,erg\,s$^{-1}$, up to $\sim10^{39}$\,erg\,s$^{-1}$. Such binaries should be observable as X-ray bright systems at 0.1\,$Z_{\odot}$, 0.5\,$Z_{\odot}$ and $Z_{\odot}$, representative of Sextans A, the Large Magellanic Cloud (LMC) and the Solar neighbourhood, respectively. The formation efficiency of these systems increases with decreasing metallicity. However, accounting for the local star formation rates, our population synthesis predicts $\sim$2 and $\sim$1 such observable binaries in the Milky Way and LMC, respectively, that will produce a binary compact object merger within a Hubble time. (Abridged)

TeV-emitting BL Lac type blazars represent the extreme end of the blazar population. They are characterized by relatively weak jets and radiatively inefficient accretion disks. Particles accelerated in these jets experience fewer radiative losses, allowing them to reach energies beyond the TeV scale and produce TeV gamma-ray emission. The study of TeV blazars is constrained by the limited number of known sources in this category. Currently, only 56 high synchrotron-peaked BL Lacs have been detected at energies above 0.1\,TeV. Searches for TeV emission from BL~Lacs typically target sources with bright X-ray emission and a synchrotron peak at or above 1\,keV. The recently released eRASS1 catalog by the eROSITA collaboration, which covers half of the sky, represents the deepest X-ray survey in the soft X-ray band to date. Utilizing the eROSITA survey, combined with infrared data from WISE and archival radio observations, we have identified 135 TeV-emitting blazar candidates. Our search introduces selection criteria based on the radio to infrared that remove quasar-like objects that have similar infrared spectra and X-ray fluxes as TeV-emitting BL~Lacs. In our search, we find 26 objects that had not been detected in the ROSAT X-ray survey and 20 that have not been previously associated with blazars. The candidates resulting from our search are suitable for follow-up observations with currently operating imaging atmospheric Cherenkov telescopes, as well as future facilities like the CTAO Observatory.

Gravitationally lensed supernovae (glSNe) are a powerful tool for exploring the realms of astronomy and cosmology. Time-delay measurements and lens modeling of glSNe can provide a robust and independent method for constraining the expansion rate of the universe. The study of unresolved glSNe light curves presents a unique opportunity for utilizing small telescopes to investigate these systems. In this work, we investigate diverse observational strategies for the initial detection of glSNe using the 7-Dimensional Telescope (7DT), a multi-telescope system composed of twenty 50-cm telescopes. We implement different observing strategies on a subset of 5807 strong lensing systems and candidates identified within the Dark Energy Camera Legacy Survey (DECaLS), as reported in various publications. Our simulations under ideal observing conditions indicate the maximum expected annual detection rates for various glSNe types (Type Ia and core-collapse (CC)) using the 7DT target observing mode in the $r$-band at a depth of 22.04 mag, as follows: 8.25 events for type Ia, 3.06 for type Ic, 1.04 for type IIb, 0.66 for type IIL, 0.89 for type IIn, 5.26 for type IIP, and 1.44 for type Ib. Furthermore, in the case of medium-band filter observations (m6000) at a depth of 20.61 in the Wide-Field Survey (WFS) program, the predicted detection rate for glSNe Ia is 3.30 $yr^{-1}$. Given targeted follow-up observations of these initially detected systems with more powerful telescopes, we can apply a model-independent approach to forecast the ability to measure $H_{0}$ using a Gaussian process from Type Ia Supernovae (SNe Ia) data and time-delay distance information derived from glSNe systems, which include both Ia and CC types. We forecast that the expected detection rate of glSNe systems can achieve a $2.4\%$ precision in estimating the Hubble constant.

In this work, we report for the first time two repeated quasi-periodic oscillations (QPOs) in the light curve of the Flat Spectrum Radio Quasar (FSRQ) S5 1044+71. This source was observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple sectors. We used the generalized Lomb-Scargle periodogram method and weighted wavelet Z-transform method to search for significant periodic signals. The main results are as follows: We found QPOs of $\sim$ 7.0 days (persisted for 4 cycles, with a significance of $\sim3.5\sigma$) and $\sim$ 7.3 days (persisted for 5 cycles, with a significance of $\sim3.8\sigma$) in the light curves of Sector 47 and EP1, respectively. Considering range of error, we consider them to be the same. We discussed two likely models of these rapid quasi-periodic variations: One comes from the jet and the other from the accretion disk. For the first one, we consider kink instability of the jet as a plausible explanation. Second, the QPO is probable to come from the main hot spots in the accretion disk, which is located approximately within the innermost stable circular orbit allowed by general relativity. Based on this model, we estimate the mass of the black hole in S5 1044+71 to be $3.49 \times 10^9 M_{\odot}$.

The Ground-based Wide-Angle Cameras array (GWAC) necessitates the integration of over 100 hardware devices, more than 100 servers, and upwards of 2500 software modules, all synchronized within a 3-second imaging cycle. However, the complexity of real-time and high concurrency processing of big data have historically resulted in a substantial failure rate, with estimated observation efficiency of less than 50% in 2023. To address these challenges, we developed a monitoring system aimed at enhancing fault diagnosis efficiency. This paper details the system's architecture, data collection methods, and the design philosophy of monitoring views. After a year of practical fault diagnostics, the system has demonstrated the ability to identify and localize faults within minutes, achieving fault localization speeds nearly ten times faster than traditional methods. Additionally, the system's design exhibits high generalizability, making them applicable to other telescope array systems.

We investigate the presence and kinematics of NV absorption proximate to high redshift quasars selected upon the presence of strong $H_{2}$ and HI absorption at the quasar redshift. Our spectroscopic observations with X-shooter at the VLT reveal a 70% detection rate of NV (9 of 13 quasars with 2.5 < z < 3.3), remarkably higher than the 10% detection rate in intervening DLA systems and the 30% rate observed within a few thousand km/s of the source in the general quasar population. While many NV components lie within the velocity range of the neutral gas, the kinematic profiles of high-ionization species appear decoupled from those of low-ionization species, with the former extending over much larger velocity ranges, particularly towards bluer velocities. We also observe significant variations in the NV/SiIV, which we attribute to varying ionization conditions, with a velocity-dependent trend: blueshifted NV components systematically exhibit higher ionization parameters compared to those near the quasar's systemic redshift. Furthermore, the most redshifted systems relative to the quasar show no evidence of NV absorption. The results suggest that proximate $H_{2}$ absorption systems select critical stages of quasar evolution, during which the quasar remains embedded in a rich molecular environment. Redshifted systems trace infalling gas, potentially associated with mergers, preceding the onset of outflows. Such outflows may reach or even carry out neutral and molecular this http URL latter stage would correspond to proximate $H_{2}$ systems located around or blueshifted relative to the quasar's systemic z. Finally, the only case in our sample featuring highly blueshifted neutral gas shows no evidence of an association with the this http URL findings highlight the need to account for the ionization state when defining a velocity threshold to distinguish quasar-associated systems from intervening.

The amplitude, polarization, frequency spectrum and energy fluence carried by the electric field at a given measurement position are the key parameters for retrieving information from radio signals generated by extensive air showers. Accurate reconstruction of the electric field from the signals recorded by the antennas is therefore essential for the radio detection technique. Conventional reconstruction methods primarily focus on electric field reconstruction for antennas with two horizontal polarizations. In this paper, we introduce an analytical least-squares ($\chi^2$) reconstruction method that operates with both two and three polarizations, providing the reconstructed electric field directly at each antenna. This solution has been verified for simple and realistic antenna responses, with a particular focus on inclined air showers. Our method achieves an accuracy better than 4\% in determining the Hilbert peak amplitude of an electric field and better than 6\% accuracy, with minimal bias, when estimating energy fluence. Additionally, this method is reliable for almost all arrival directions, and the direction dependence has been investigated. This work also demonstrates that incorporating vertically polarized antennas enhances the precision of reconstruction, leading to a more accurate and reliable electric field estimation for inclined air showers. Consequently, the method enhances our ability to extract information about cosmic rays from the detected signals in current and future experiments.

Reconstructing the dark matter density, velocity, and tidal (MTV) fields from galaxy surveys is essential for advancing our understanding of the large-scale structure of the Universe. In this work, we present a machine learning-based framework using a UNet convolutional neural network to reconstruct the MTV fields from mock samples of the DESI bright galaxy survey within the redshift range $0.1 < z < 0.4$. Our approach accounts for realistic observational effects, including geometric selection, flux-limited data, and redshift space distortion (RSD) effects, thereby improving the fidelity of the reconstructed fields. Testing on mock galaxy catalogs generated from the Jiutian N-body simulation, our method achieves significant accuracy level. The reconstructed density field exhibits strong consistency with the true field, effectively eliminating most RSD effects and achieving a cross-correlation power spectrum coefficient greater than 0.985 on scales with $k < 0.1 \, h \, \mathrm{Mpc}^{-1}$. The velocity field reconstruction accurately captures large-scale coherent flows and small-scale turbulent features, exhibiting slopes of the grid-to-grid relationship close to unity and a scatter below $\sim$100 $\mathrm{km} \, \mathrm{s}^{-1}$. Additionally, the tidal field is reconstructed without bias, successfully recovering the features of the large-scale cosmic web, including clusters, filaments, sheets, and voids. Our results confirm that the proposed framework effectively captures the large-scale distribution and dynamics of dark matter while addressing key systematic challenges. These advancements provide a reliable and robust tool for analyzing current and future galaxy surveys, paving the way for new insights into cosmic structure formation and evolution.

The spectra of galactic cosmic rays (GCRs) contain crucial information about their origin and propagation through the interstellar medium. When GCRs reach Earth, they are significantly influenced by the solar wind and the heliospheric magnetic field, a phenomenon known as solar modulation. This effect introduces time-dependent variations in GCR fluxes. The AMS-02 experiment has released time-dependent flux data for protons, electrons, and positrons, revealing clear correlations with solar modulation. Studies suggest that cosmic rays with the same charge, such as protons and helium nuclei, exhibit similar/same solar modulation parameters. In this work, we derive the LIS for protons and positrons under the assumption of a common solar modulation potential, using data from Voyager 1 and a 7-year average from AMS-02. Similarly, the LIS for antiprotons and electrons is derived by assuming they are governed by a separate solar modulation potential. We demonstrate that the time-dependent fluxes of positrons and protons can be accurately modeled using the same set of solar modulation parameters within a modified force-field approximation framework. Based on this, we predict the time-dependent fluxes of antiprotons using the corresponding electron flux data.

Decaying dark matter (DDM) affects the evolution of cosmic structure relative to standard cold, collisionless, stable dark matter (CDM). We introduce a new semi-analytic model for the effects of two-body DDM on halo structure and subhalo populations. In this scenario, cold parent dark matter particles decay into less massive daughter particles plus dark radiation with a lifetime comparable to the age of the universe. Our DDM model is implemented in the open-source software $\texttt{Galacticus}$ and accounts for heating (due to the velocity kicks imparted on daughter particles) and mass loss (due to the parent--daughter mass splitting). We show that decays flatten and reduce the amplitude of halos' inner density profiles. These effects make DDM subhalos susceptible to tidal disruption, which we show yields a mass-dependent suppression of the subhalo mass function relative to CDM. Our predictions for DDM density profiles, velocity dispersion profiles, and subhalo populations are consistent with results from isolated and cosmological DDM N-body simulations. Thus, our model enables efficient and accurate exploration of DDM parameter space and will be useful for deriving constraints from upcoming small-scale structure observations.

We present a theoretical model to investigate a two-planet pair that undergoes convergent type I migration in a gaseous protoplanetary disk and traps into the first-order mean motion resonance. Our study identifies the conditions for resonant capture and explores the subsequent dynamical stability of the system. We derive the analytical criteria for planets with an arbitrary mass ratio and validate through numerical N-body simulations. Slow migration and weak eccentricity damping are required for resonant capture, the latter of which has not received enough attention in the literature. Once capture into resonance, their following stability can be classified into three regimes: stable trap, overstable trap and escape. Notably, resonant capture remains stable when the inner planet significantly outweighs the outer one. However, the subsequent evolution can be diverse when the mass of the inner planet is lower or comparable to that of the outer planet. Stability weakens as the relative strength between migration and eccentricity damping increases. The key results can be comprehensively demonstrated in a $\tau_{a}$-$\tau_{a}/\tau_{\rm e}$ plot, where $\tau_{a}$ and $\tau_{\rm e}$ are orbital decay and eccentricity damping timescales, respectively.

To study the formation of star clusters and their properties in a dwarf-dwarf merging galaxy, we have performed a numerical simulation of a dwarf-dwarf galaxy merger by using the Tree+GRAPE $N$-body/SPH code ASURA. In our simulation, 13 young star clusters are formed during the merger process. We show that our simulated star clusters can be divided into two types: with and without [Fe/H] abundance variations. The former is created by a seed star cluster (the first-generation stars) formed in compressed gas. These stars contaminate the surrounding gas by Type II supernovae (SNe). At that time, the energy injection is insufficient to induce an outflow of the surrounding gas. After that, the contaminated gas falls into the seed, thereby forming a new generation of stars from the contaminated gas. We also show that most star clusters are formed in the galactic central region after the second encounter and fall into the galactic center due to dynamical friction within several hundred Myr. As a result, close encounters and mergers between the clusters take place. Although the clusters with shallower gravitational potential are tidally disrupted by these close encounters, others survive and finally merge at the center of the merged dwarf galaxies to create a nuclear star cluster. Therefore, the nuclear star cluster is comprised of various stellar components in [Fe/H] abundance and age. We discuss our work in the context of observations and demonstrate the diagnostic power of high-resolution simulations in the context of star cluster formation.

We investigate the angular power spectrum ($C_\ell)$ and angular correlation function ($w(\theta)$) of galaxy number density field in the presence of the local-type primordial non-Gaussianity (PNG), explicitly accounting for the integral constraint in an all-sky survey. We show that the PNG signature in $C_{\ell}$ is confined to low multipoles in the linear regime, whereas its signature in $w(\theta)$ extends across a wide range of angular scales, including those below the nonlinear scale. Therefore, the equivalence between $C_\ell$ and $w(\theta)$ can be violated when scale cuts of multipoles or angular scales -- for example, to mitigate systematic effects -- are applied in the analysis. Assuming samples of photometric galaxies divided into multiple redshift bins in the range $0<z<7$, we forecast the precision of constraining the PNG parameter ($f_{\rm NL}$) from the hypothetical measurements of $C_\ell$ or $w(\theta)$ assuming different scale cuts in the multipoles or angular scales, respectively. Our results imply that the PNG information can be extracted from $w(\theta)$ on relatively small angular scales such as $\lesssim 10$ degree for a high-redshift galaxy sample or from $w(\theta)$ measured in a survey with partial area coverage.

The wide survey of the Chinese Space Station Telescope (CSST) will observe a large field of 17,500 $\text{deg}^2$. The GU, GV, and GI grism observations of CSST will cover a wavelength range from 2550 to 10000Å at a resolution of $R\sim 200$ and a depth of about 22 AB magnitude for the continuum. In this paper, we present a pipeline to identify quasars and measure their physical properties with the CSST mock data. We simulate the raw images and extract the one-dimensional grism spectra for quasars, galaxies, and stars with the r-band magnitudes of $18<\text{m}_{\text{r}}<22$ using the CSST Cycle 6 simulation code. Using a convolution neural network, we separate quasars from stars and galaxies. We measure the redshifts by identifying the strong emission lines of quasars. We also fit the 1D slitless spectra with QSOFITMORE to estimate the black hole masses and Eddington ratios. Our results show that the CSST slitless spectroscopy can effectively separate quasars with redshifts $z=0-5$ from other types of objects with an accuracy of 99\%. Among those successfully classified quasars, 90\% of them could have precise redshift measurements with $\sigma_{\mathrm{NMAD}}=0.002$. The scatters of black hole masses and Eddington ratios from the spectral fittings are 0.13 and 0.15 dex, respectively. The metallicity diagnosis line ratios have a scatter of 0.1-0.2 dex. Our results show that the CSST slitless spectroscopy survey has the potential to discover about 0.9 million new quasars and provide important contributions to AGN science and cosmology.

We have been developing silicon-on-insulator (SOI) pixel detectors with a pinned depleted diode (PDD) structure, named "XRPIX", for X-ray astronomy. In our previous study, we successfully optimized the design of the PDD structure, achieving both the suppression of large leakage current and satisfactory X-ray spectroscopic performance. Here, we report a detailed study on the X-ray spectroscopic performance of the XRPIX with the optimized PDD structure. The data were obtained at $-60^\circ\mathrm{C}$ with the "event-driven readout mode", in which only a triggering pixel and its surroundings are read out. The energy resolutions in full width at half maximum at 6.4 keV are $178\pm1$ eV and $291\pm1$ eV for single-pixel and all-pixel event spectra, respectively. The all-pixel events include charge-sharing pixel events as well as the single-pixel events. These values are the best achieved in the history of our development. We argue that the gain non-linearity in the low energy side due to excessive charge injection to the charge-sensitive amplifier is a major factor to limit the current spectroscopic performance. Optimization of the amount of the charge injection is expected to lead to further improvement in the spectroscopic performance of XRPIX, especially for the all-pixel event spectrum.

We present a study of the radio and optical properties of the high-frequency peaker (HFP) blazar PKS 1614+051 at $z=3.21$ based on the data covering the time period of 1997-2024. The radio data are represented by the almost instantaneous 1-22 GHz measurements from the SAO RAS RATAN-600 radio telescope, the 5 and 8 GHz data from the IAA RAS RT-32 telescopes, and the 37 GHz data from the RT-22 telescope of CrAO RAS. The optical measurements in the $R$ band were collected with the SAO RAS 1-m Zeiss-1000 and 0.5-m AS-500/2 telescopes and the ZTF archive data. We have found low overall variability indices (0.1-0.2) and a spectral peak around 4.6 GHz, which is stable during the long-term period of monitoring. An analysis of the radio light curves reveals significant time delays (0.6 to 6.4 years) between the radio frequencies along with variability timescales ranging from 0.2 to 1.8 years in the source's rest frame, which is similar to the blazars at lower redshifts. Spectral modeling suggests the presence of both synchrotron-self absorption (SSA) and free-free absorption (FFA) processes. Based on the SSA model, we provide estimates of the magnetic field strength, which peaks at $\sim\!30$ mG. A spectroscopic study with the BTA SCORPIO-I spectrograph has found signs of the regular motion of a neutral hydrogen envelope around the blazar center, which confirms the presence of enough gaseous matter to form an external FFA screen. The results highlight the importance of multiwavelength and long-term monitoring to understand the physical mechanisms driving the variability in high-redshift blazars.

This paper presents an analysis of the predicted optical-to-X-ray spectral index ($\alpha_{\rm ox}$) within the context of ultra-luminous X-ray sources (ULXs) associated with stellar mass black holes and neutron stars. We use the population synthesis code COSMIC to simulate the evolution of binary systems and investigate the relationship between UV and X-ray emission during the ULX phase, namely the $\alpha_{\rm ox}$ relation. The study investigates the impact of metallicity on $\alpha_{\rm ox}$ values. Notably, it predicts a significant anti-correlation between $\alpha_{\rm ox}$ and UV luminosity ($L_{\rm UV}$), consistent with observations, with the slope of this relationship varying with metallicity for BH-ULXs. The NS-ULX population shows a relatively consistent slope around $-0.33$ across metallicities, with minor variations. The number of ULXs decreases with increasing metallicity, consistent with observational data, and the X-ray luminosity function shows a slight variation in its slope with metallicity, exhibiting a relative excess of high-luminosity ULXs at lower metallicities. Inclusion of beaming effect in the analysis shows a significant impact on the XLF and $\alpha_{\rm ox}$, particularly at high accretion rates, where the emission is focused into narrower cones. Furthermore, the study finds that UV emission in ULXs is predominantly disk-dominated, which is the likely origin of the $\alpha_{\rm ox}$ relation, with the percentage of disk-dominated ULXs increasing as metallicity rises.

During its annual conference in 2024, the French Society of Astronomy \&amp; Astrophysics (SF2A) hosted, for the fourth time, a special session dedicated to discussing the environmental transition within the French astrophysics research community. This year had a special context: both the CNRS-INSU and the CNES were preparing their scientific perspectives for the period 2025-2030 in the field of Astronomy-Astrophysics (AA). In this proceeding, we first describe the main actions undertaken by the Commission Transition Environnementale. Then, we summarize the discussions held during the half-day workshop, which brought together about 100 participants, and point to forthcoming proceeding, reports and other related resources. A key message is that the French A\&amp;A community is now fully aware that astronomical activities simply cannot thrive indefinitely in the current situation, and seems now eager to seize the opportunity of developing our profession towards a better social and environmental impact.

Supernova remnants (SNRs) are among the most important sources of non-thermal X-rays in the sky and likely contributors to Galactic cosmic rays and represent ideal targets to showcase the capabilities of the Imaging X-ray Polarimetry Explorer (IXPE) in performing spatially-resolved X-ray polarimetry. For the first time, we can determine the turbulence level (through the measurement of polarization degree) and the orientation (through polarization direction) of the magnetic field near the shocks, where particle acceleration occurs. IXPE reported so far the results of the observations of four SNRs: Cas A, Tycho, SN 1006 and RX J1713.7-3946. These objects exhibit a wide range of characteristics, including dynamical age, spectral composition of emission, and progenitor type. Aptly, they revealed significantly different results among them in terms of magnetic field properties and morphology, providing unexpected insights and shedding light on the particle acceleration mechanisms in astrophysical shocks.

The masses of Population III stars are largely unconstrained since no simulations exist that take all relevant primordial star formation physics into account. We perform the first suite of radiation magnetohydrodynamics (RMHD) simulations of Population III star formation, with the POPSICLE project. Compared to control simulations that only include magnetic fields (MHD), protostellar ionizing and dissociating feedback, or neither, the RMHD simulation best resembles the MHD simulation during the earliest stages of collapse and star formation. In $5000\,\rm{yrs}$, the mass of the most massive star is $65\,\rm{M_{\odot}}$ in the RMHD simulation, compared to $120\,\rm{M_{\odot}}$ in simulations without magnetic fields. This difference arises because magnetic fields act against gravity, suppress mass transport, and reduce compressional heating. The maximum stellar mass of Population III stars is thus already limited by magnetic fields, even before accretion rates drop to allow significant protostellar radiative feedback. Following classical main sequence stellar evolution with MESA reveals that it is difficult to create Population III stars with masses larger than $600\,\rm{M_{\odot}}$ in typical dark matter minihaloes at $z \gtrsim 20$, with maximum stellar masses $\sim 100\,\rm{M_{\odot}}$ more likely due to expected negative feedback from both magnetic fields and stellar radiation. This work lays the first step in building a full physics-informed mass function of Population III stars.

Context. Solar filament oscillations have been observed for many years, but recent advances in telescope capabilities now enable daily monitoring of these periodic motions, offering valuable insights into the structure of filaments. A systematic study of filament oscillations over the solar cycle can shed light on the evolution of the prominences. Until now, only manual techniques have been used to analyze these oscillations. Aims. This work serves as a proof of concept, aiming to demonstrate the effectiveness of Convolutional Neural Networks (CNNs). These networks automatically detect filament oscillations by applying power-spectrum analysis to H$\alpha$ data from the GONG telescope network. Methods. The proposed technique studies periodic fluctuations of every pixel of the H$\alpha$ data cubes. Using the Lomb Scargle periodogram, we computed the power spectral density (PSD) of the dataset. The background noise is well fitted to a combination of red and white noise. Using Bayesian statistics and Markov chains -- Monte Carlo (MCMC) algorithms we can fit the spectra and determine the confidence threshold of a given percentage to search for real oscillations. We built two CNN models to obtain the same results as the MCMC approach. Results. We applied the CNN models to some observations reported in the literature, proving its reliability in detecting the same events as the classical methods. A day with events not previously reported has been studied to check for the model capabilities outside of a controlled dataset where we can check with previous reports. Conclusions. CNNs proved to be a useful tool to study solar filament oscillations, using spectral techniques. Computing times have been reduced significantly while getting results similar enough to the classical methods. This is a relevant step towards the automatic detection of filament oscillations.

Light-travel-time delays provide one of the most powerful ways of learning about the structure and kinematics of active galactic nuclei (AGNs). Estimating delays from observations of AGN variability presents statistical challenges because the time series are almost invariably irregularly sampled. Correct assessment of errors in lag estimates is important for evaluating results. The most widely used method of determining phase lags has been the interpolated cross-correlation function method introduced by Gaskell and Sparke (1986, GS). It is shown here that the widely used modified smooth bootstrap method of Peterson et al. (1998) significantly overestimates the error in lags derived using the GS, especially for poorly sampled light curves. The remarkably high accuracy claimed for lags obtained by the JAVELIN method of Zu et al. (2011) (more than an order of magnitude improvement for 25% of cases) is spurious and the accuracy no better than the GS method. A related method proposed by Li et al. (2013) suffers from similar but less serious problems. A slightly modified version of the analytic formula of Gaskell and Peterson (1987) gives an easy and accurate way of estimating lag errors for well-sampled, high-quality data. The width of the continuum autocorrelation function is shown to be proportional to the square root of the luminosity. This is useful in planning observing campaigns. The discrete correlation function method is less powerful than the GS method and for many typical situations gives errors in the lag twice as large as those of the GS method.

The subsurface convective zones (CZs) of massive stars significantly influences many of their key characteristics. Previous studies have paid little attention to the impact of rotation on the subsurface convective zone (CZ), so we aim to investigate the evolution of this zone in rapidly rotating massive stars. We use the Modules for Experiments in Stellar Astrophysics (MESA) to simulate the subsurface CZs of massive stars during the main sequence phase. We establish stellar models with initial masses ranging from 5 $M_{\odot}$ to 120 $M_{\odot}$, incorporating four metallicities (Z = 0.02, 0.006, 0.002, and 0.0001) and three rotational velocities ($\mit\omega/\omega_{\text {crit}}$ = 0, $\mit\omega/\omega_{\text {crit}}$ = 0.50, and $\mit\omega/\omega_{\text {crit}}$ = 0.75). We find that rapid rotation leads to an expansion of the subsurface CZ, increases convective velocities, and promotes the development of this zone. Additionally, subsurface CZs can also emerge in stars with lower metallicities. Comparing our models with observations of massive stars in the Galaxy, the Large Magellanic Cloud, and the Small Magellanic Cloud, we find that rotating models better encompass the observed samples. Rotation significantly influences the evolution of the subsurface CZ in massive stars. By comparing with the observed microturbulence on the surfaces of OB stars, we propose that the subsurface CZs may be one of the sources of microturbulence.

We study the impact of bar-induced non-circular motions on the derivation of galactic rotation curves (RCs) using hydrodynamic simulations and observational data from the PHANGS-ALMA survey. We confirm that non-circular motions induced by a bar can significantly bias RCs derived from the conventional tilted-ring method, consistent with previous findings. The shape of the derived RC depends on the position angle difference ($\Delta \phi$) between the major axes of the bar and the disk in the face-on plane. For $\left|\Delta \phi\right|\lesssim40^\circ$, non-circular motions produce a bar-induced "dip" feature (rise-drop-rise pattern) in the derived RC, which shows higher velocities near the nuclear ring and lower velocities in the bar region compared to the true RC (${\mathrm{RC_{true}}}$). We demonstrate convincingly that such dip features are very common in the PHANGS-ALMA barred galaxies sample. Hydrodynamical simulations reveal that the "dip" feature is caused by the perpendicular orientation of the gas flows in the nuclear ring and the bar; at low $\left|\Delta \phi\right|$ streamlines in the nuclear ring tend to enhance $V_\mathrm{los}$, while those in the bar tend to suppress $V_\mathrm{los}$. We use a simple {\misaell} model to qualitatively explain the general trend of RCs from the tilted-ring method (${\mathrm{RC_{tilted}}}$) and the discrepancy between ${\mathrm{RC_{tilted}}}$ and ${\mathrm{RC_{true}}}$. Furthermore, we propose a straightforward method to implement a first-order correction to the RC derived from the tilted-ring method. Our study is the first to systematically discuss the bar-induced "dip" feature in the RCs of barred galaxies combining both simulations and observations.

Massive Star Clusters (SCs) have been proposed as important CR sources, with the potential of explaining the high-energy end of the Galactic cosmic-ray (CR) spectrum, that Supernova Remnants (SNRs) seem unable to account for. Thanks to fast mass losses due to the collective stellar winds, the environment around SCs is potentially suitable for particle acceleration up to PeV energies and the energetics is enough to account for a large fraction of the Galactic CRs, if the system is efficient enough. A handful of star clusters have been detected in gamma-rays confirming the idea that particle acceleration is taking place in this environment. However, contamination by other sources often makes it difficult to constrain the contribution arising from SCs only. Here we present a new analysis of Fermi-LAT data collected towards a few massive young star clusters. The young age (< 3 Myr) of the clusters guarantees that no SN has exploded in the region, allowing us to determine the power contributed by the stellar component alone, and to quantify the contribution of this type of sources to the bulk of CRs. Moreover, we will present a recent statistical investigation that quantifies the degree of correlation between gamma-ray sources and these astrophysical objects and briefly discuss the observational prospect for ASTRI and CTAO.

The Cryogenic AntiCoincidence Detector (CryoAC) is a key element of the X-ray Integral Field Unit (X-IFU) on board the future ATHENA X-ray observatory. It is a TES-based detector designed to reduce the particle background of the instrument, thereby increasing its sensitivity. The detector design is driven by an end-to-end simulator which includes the electro-thermal modelling of the detector and the dynamics of its readout chain. Here, we present the measurements carried out on the last CryoAC single pixel prototype, namely DM127, in order to evaluate the critical thermal parameters of the detector and consequently to tune and validate the CryoAC end-to-end simulator.

We fit X-ray reverberation models to Rossi X-ray Timing Explorer data from the X-ray binary Cygnus X-1 in its hard state to yield estimates for the black hole mass and the distance to the system. The rapid variability observed in the X-ray signal from accreting black holes provides a powerful diagnostic to indirectly map the ultra-compact region in the vicinity of the black hole horizon. X-ray reverberation mapping exploits the light crossing delay between X-rays that reach us directly from the hard X-ray emitting 'corona', and those that first reflect off the accretion disc. Here we build upon a previous reverberation mass measurement of Cygnus X-1 that used the RELTRANS software package. Our new analysis enhances signal to noise with an improved treatment of the statistics, and implements new RELTRANS models that are sensitive to distance. The reduced uncertainties uncover evidence of mass accretion rate variability in the inner region of the disc that propagates towards the corona. We fit two different distance-sensitive models, and both return reasonable values of distance and mass within a factor 2 of the accepted values. The models both employ a point-like 'lamppost' corona and differ only in their treatment of the angular emissivity of the corona. The two models return different mass and distance estimates to one another, indicating that future reverberation models that include an extended corona geometry can be used to constrain the shape of the corona if the known mass and distance are utilised via Bayesian priors.

ExoSim 2 is the next generation of the Exoplanet Observation Simulator (ExoSim) tailored for spectro-photometric observations of transiting exoplanets from space, ground, and sub-orbital platforms. This software is a complete rewrite implemented in Python 3, embracing object-oriented design principles, which allow users to replace each component with their functions when required. ExoSim 2 is publicly available on GitHub, serving as a valuable resource for the scientific community. ExoSim 2 employs a modular architecture using Task classes, encapsulating simulation algorithms and functions. This flexible design facilitates the extensibility and adaptability of ExoSim 2 to diverse instrument configurations to address the evolving needs of the scientific community. Data management within ExoSim 2 is handled by the Signal class, which represents a structured data cube incorporating time, space, and spectral dimensions. The code execution in ExoSim 2 follows a three-step workflow: the creation of focal planes, the production of Sub-Exposure blocks, and the generation of non-destructive reads (NDRs). Each step can be executed independently, optimizing time and computational resources. ExoSim 2 has been extensively validated against other tools like ArielRad and has demonstrated consistency in estimating photon conversion efficiency, saturation time, and signal generation. The simulator has also been validated independently for instantaneous read-out and jitter simulation, and for astronomical signal representation. In conclusion, ExoSim 2 offers a robust and flexible tool for exoplanet observation simulation, capable of adapting to diverse instrument configurations and evolving scientific needs. Its design principles and validation results underscore its potential as a valuable resource in the field of exoplanet research.

Early energy injection leaves an imprint on the observed blackbody spectrum of the CMB, allowing us to study the thermal history of the Universe. For small energy release, the distortion can be efficiently computed using the quasi-exact Green's function method. For pre-recombination injections, the Green's function has already been studied previously. Here we reconsider the pre- and post-recombination periods, showcasing both the spectral distortion intensity and the relative temperature difference, which encrypt precious information about physical processes such as free-free interactions and thermal decoupling. We present the associated distortion visibility function, investigating the impact of various physical effects. We then study improvements to the so-called frequency hierarchy (FH) treatment, a method that was developed for the modelling of anisotropic distortions, which like the average distortion signals encode valuable cosmological information. Specifically, the FH treatment has shortcomings even in the $\mu$ era, that in principle should be easy to overcome. In this paper, we introduce a new approach to reduce the mismatch, concluding with a redefinition of the $\mu$ spectral shape using CosmoTherm. This solution takes into account double Compton and Bremsstrahlung effects in the low tail, which can be included in the FH. This opens the path towards a refined modeling of spectral distortion anisotropies.

Context. Classical chemically peculiar stars exhibit atmospheres that are often structured by the effects of atomic diffusion. As a result of these elemental diffusion and horizontal abundance variations, photospheric temperature varies at a given height in the atmosphere. This may lead to horizontal flows in the photosphere. In addition, the suppression of such flows by magnetic field can alter the elemental transport processes. Aims. Using a simplified model of such a structured atmosphere and 2D MHD simulations of a typical He-rich star, we examine atmospheric flows in these chemically peculiar stars which often are strongly magnetic. Methods. We use Zeus-MP which is a Fortran 90 based publicly available parallel finite element modular code. Results. We find that for non-magnetic stars of spectral type BA, atmospheric flow related to horizontal temperature gradient can reach 1.0 km/s yielding mixing timescale of order of tens of days. For the magnetic counterparts, the flow speeds are an order of magnitude lower allowing for stratification of chemical elements. Conclusions. Magnetic field can influence the dynamics in atmospheres significantly. Strong horizontal magnetic field inhibits flow in the vertical direction, while strong vertical magnetic field can suppress horizontal atmospheric flow preventing elemental mixing.

We study the influence of hyperons in binary neutron star (NS) mergers considering a total of 14 temperature dependent equations of state (EoSs) models which include hyperonic degrees of freedom and partly delta resonances. Thermally produced hyperons induce a higher heat capacity and a lower thermal index, i.e. a reduced thermal pressure for a given amount of thermal energy, compared to purely nucleonic models. We run a large set of relativistic hydrodynamics simulations of NS mergers to explore the impact on observables of these events. In symmetric binaries, we describe a characteristic increase of the dominant postmerger gravitational-wave (GW) frequency by a few per cent, which is specifically linked to the occurrence of hyperons and can thus be potentially used as a discriminator between purely nucleonic and hyperonic systems. We corroborate that this effect occurs similarly for asymmetric binaries and becomes more prominent with increasing total binary mass. Hyperonic models tend to stick out in relations between the dominant postmerger GW frequency and the tidal deformability of massive stars providing a signature to identify the presence of hyperons. Distinct secondary postmerger GW spectral features are differently affected by the presence of hyperons, in the sense that one feature exhibits a characteristic frequency shift due to the specific thermal properties of hyperonic EoSs while the other does not. The dynamical mass ejection of mergers is tentatively enhanced for hyperonic models in comparison to nucleonic EoSs which yield roughly the same stellar properties of cold NSs. This may serve useful to identify exotic degrees of freedom in these systems by kilonova observations. Also, for hyperonic EoSs the threshold mass for prompt black hole formation is reduced by about 0.05 $M_{\odot}$ in comparison to nucleonic systems with the same stellar parameters of cold NSs.

For minor planet observations to be archived and used by the scientific community, observations of individual objects must be linked together into groups called tracklets. This linking is nontrivial, as linking software must find real tracklets in noisy data within a reasonable amount of time. We describe FindPOTATOs, a linking software written in Python that can assemble minor planet tracklets. With the appropriate parameters, FindPOTATOs assembles tracklets for a variety of objects, including close-approaching near-Earth objects and trans-Neptunian objects. FindPOTATOs is ideal for processing data sets taken at small observatories, processing archival data, or finding tracklets of minor planets on unusual orbits. This paper describes the code structure, usage, and validation.

The G-ring arc of Saturn, confined by the 7:6 corotation eccentric resonance with Mimas, is primarily composed of micrometric particles. These particles, significantly influenced by the solar radiation pressure, are subject to rapid depletion. This study investigates a mechanism for dust replenishment in the arc, specifically analyzing collisions between macroscopic bodies and the satellite Aegaeon. Utilizing N-body and Smoothed particle hydrodynamics (SPH) simulations, we assess the dust generation from these impacts, with a focus on the most likely collision parameters derived from the N-body simulations. Our findings indicate that, while collisions among macroscopic bodies are inefficient for dust production, impacts involving Aegaeon are substantially more effective, providing conservative lower limits for dust generation. This mechanism, in conjunction with the natural decay processes and continuous dust generation from impacts, potentially keeps the arc population over thousands of years with a possible variation in brightness.

First-principle modeling of dense hydrogen is crucial in materials and planetary sciences. Despite its apparent simplicity, predicting the ionic and electronic structure of hydrogen is a formidable challenge, and it is connected with the insulator-to-metal transition, a century-old problem in condensed matter. Accurate simulations of liquid hydrogen are also essential for modeling gas giant planets. Here we perform an exhaustive study of the equation of state of hydrogen using Density Functional Theory and quantum Monte Carlo simulations. We find that the pressure predicted by Density Functional Theory may vary qualitatively when using different functionals. The predictive power of first-principle simulations is restored by validating each functional against higher-level wavefunction theories, represented by computationally intensive variational and diffusion Monte Carlo calculations. Our simulations provide evidence that hydrogen is denser at planetary conditions, compared to currently used equations of state. For Jupiter, this implies a lower bulk metallicity (i.e., a smaller mass of heavy elements). Our results further amplify the inconsistency between Jupiter's atmospheric metallicity measured by the Galileo probe and the envelope metallicity inferred from interior models.

BayesEoR is a GPU-accelerated, MPI-compatible Python package for estimating the power spectrum of redshifted 21-cm emission from interferometric observations of the Epoch of Reionization (EoR). Utilizing a Bayesian framework, BayesEoR jointly fits for the 21-cm EoR power spectrum and a "foreground" model, referring to bright, contaminating emission between us and the cosmological signal, and forward models the instrument with which these signals are observed. To perform the sampling, we use MultiNest [arXiv:1402.0004], which calculates the Bayesian evidence as part of the analysis. Thus, BayesEoR can also be used as a tool for model selection [see e.g. arXiv:1701.03384].

AGILE (Astrorivelatore Gamma ad Immagini LEggero) has been a unique and hugely successful mission of Italian Space Agency (ASI), built and operated with the programmatic and technical support of the National Institute for Astrophysics (INAF), the National Institute for Nuclear Physics (INFN), several universities and Italian industries. During almost 17 years of observations in orbit (from April 23, 2007, to January 18, 2024), AGILE contributed to high-energy astrophysics and terrestrial physics with many discoveries and detections. Two co-aligned X- and gamma-ray detectors, a silicon-strip-based tracker, a wide field of view gamma-ray imager and the fast-reaction ground segment were the AGILE innovative solutions with respect to the previous generation of gamma-ray satellites. With the AGILE's re-entry, the in-orbit operational phase ends, but a new phase of scientific work on the satellite legacy data archive opens: AGILE may still hold future surprises.

The Harvard College Observatory was the preeminent astronomical data center of the early 20th century: it gathered and archived an enormous collection of glass photographic plates that became, and remains, the largest in the world. For nearly twenty years DASCH (Digital Access to a Sky Century @ Harvard) actively digitized this library using a one-of-a kind plate scanner. In early 2024, after 470,000 scans, the DASCH project finished. Now, this unique analog dataset can be integrated into 21st-century, digital analyses. The key DASCH data products include ~200 TB of plate images, ~16 TB of calibrated light curves, and a variety of supporting metadata and calibration outputs. Virtually every part of the sky is covered by thousands of DASCH images with a time baseline spanning more than 100 years; most stars brighter than B ~ 15 have hundreds or thousands of detections. DASCH Data Release 7, issued in late 2024, represents the culmination of the DASCH scanning project.

Population III stars are possible precursors to early massive and supermassive black holes (BHs). The presence of soft UV Lyman Werner (LW) background radiation can suppress Population III star formation in minihalos and allow them to form in pristine atomic cooling halos. In the absence of molecular hydrogen ($\rm H_2$) cooling, atomic-cooling halos enable rapid collapse with suppressed fragmentation. High background LW fluxes from preceding star-formation have been proposed to dissociate $\rm H_2$. This flux can be supplemented by LW radiation from one or more Population III star(s) in the same halo, reducing the necessary background level. Here we consider atomic-cooling halos in which multiple protostellar cores form close to one another nearly simultaneously. We assess whether the first star's LW radiation can dissociate nearby $\rm H_2$, enabling the prompt formation of a second, supermassive star (SMS) from warm, atomically-cooled gas. We use a set of hydrodynamical simulations with the code ENZO, with identical LW backgrounds centered on a halo with two adjacent collapsing gas clumps. When an additional large local LW flux is introduced, we observe immediate reductions in both the accretion rates and the stellar masses that form within these clumps. While the LW flux reduces the $\text{H}_2$ fraction and increases the gas temperature, the halo core's potential well is too shallow to promptly heat the gas to $\gtrsim$ 1000 K and increase the accretion rate onto the second protostar. We conclude that internal LW feedback inside atomic-cooling halos is unlikely to facilitate the formation of SMSs or massive BH seeds.

It is unclear which fraction of cosmic rays above an energy of 1 PeV is accelerated by the observed Galactic PeVatron population. These sources' gamma-ray data is typically degenerate between hadronic and leptonic emission scenarios, which hinders their straightforward association with the cosmic ray population. Here, we aimed to distinguish between leptonic and hadronic particle acceleration scenarios for the PeVatron candidate HESS J1646-458, associated with the star cluster Westerlund 1 (Wd 1). We first studied the diffuse X-ray emission from Wd 1 to better understand if its origin is of thermal or nonthermal nature. In addition, we searched for X-ray synchrotron emission from the associated PeVatron candidate HESS J1646-458 to put new constraints on the magnetic field strength and the leptonic particle population of this source. We used data from eROSITA on board the SRG orbital platform to spectrally analyze the diffuse emission from Wd 1 and HESS J1646-458. For Wd 1, we compared a purely thermal model and a model with a thermal and a nonthermal component. Next, we analyzed the spectra of four annuli around Wd 1 which coincide with HESS J1646-458 to search for synchrotron radiation. We find that eROSITA data cannot distinguish between thermal and nonthermal source scenarios for the diffuse emission from Wd 1 itself. In the case of HESS J1646-458, we find no evidence of synchrotron emission. We estimated an upper confidence bound of the synchrotron flux up to 40' around Wd 1 of 1.9e-3 keV-1 cm-2 s-1, which we used to study the spectral energy distribution of the source. From this, we obtained an upper 1 sigma bound on the magnetic field strength of HESS J1646-458 of 7 uG. This is compatible with a previous estimate in the literature for a fully leptonic source scenario. A purely leptonic emission scenario, a hadronic, and a hybrid one are compatible with our results.

This paper presents the development and characterization of a fiber laser hydrophone designed for deep-sea applications, with a focus on detecting neutrino interactions via their acoustic signatures. The hydrophone design includes a static pressure compensation mechanism, ensuring reliable operation at depths exceeding 1 km. The performance of the hydrophone was evaluated through laboratory tests and experiments in an anechoic basin, where its transfer function was measured before and after a 140-bar pressure cycle. The results show that the hydrophone maintains its sensitivity, with resonance peaks identified in both low- and high-frequency ranges. The hydrophone's sensitivity to acoustic signals was also compared to ambient sea state noise levels, demonstrating compatibility with the lowest noise conditions.

We present a new technique to constrain the gravitational potential of a galaxy from the observed stellar mass surface density alone under a number of assumptions. It uses the classical Eddington Inversion Method to compute the phase-space distribution function (DF) needed for the stars to reside in a given gravitational potential. In essence, each potential defines a set of density profiles, and it is the expansion of the observed profile in this database that provides the DF. If the required DF becomes negative then the potential is inconsistent with the observed stars and can be discarded. It is particularly well-suited for analyzing low-mass low surface brightness galaxies, where photometric but not spectroscopic data can be obtained. The recently discovered low surface brightness galaxy 'Nube' was used to showcase its application. For the observed Nube's stellar core to be reproduced with non-negative DF, cuspy NFW (Navarro, Frenk, and White) potentials are highly disfavored compared with potentials having cores (Schuster-Plummer or rho-230). The method assumes the stellar system to have spherical symmetry and isotropic velocity distribution, however, we discuss simple extensions that relax the need for isotropy and may help to drop the spherical symmetry assumption.

New observations obtained with JWST of the proto-stellar HH~270 jet and the "deflected" HH 110 system, show that HH 110 has a morphology of a series of distorted working surfaces. These working surfaces appear to be "deflected versions" of the heads of the incident, HH 270 jet. We compute a series of 3D numerical simulations, in which we explore the possible parameters of a shearing environment that could give origin to the deflection of HH 270 into the HH 110 flow. We find that we need a quite high sideways velocity for the streaming environment (of ~30km/s) in order to produce the complex structure observed in HH 110. This high velocity would be possible in an environment which has been strongly perturbed by the passage of other outflows.

Variability studies offer a compelling glimpse into black hole dynamics, and NICER's (Neutron Star Interior Composition Explorer) remarkable temporal resolution propels us even further. NICER observations of an Active Galactic Nucleus (AGN), NGC 4051, have charted the geometry of the emission region of the central supermassive black hole. Our investigation of X-ray variability in NGC 4051 has detected extreme variations spanning a factor of 40 to 50 over a mere 10 to 12 hours. For the first time, we have constrained the X-ray Power Spectral Density (PSD) of the source to 0.1 Hz, corresponding to a temporal frequency of 10,000 Hz in a galactic X-ray binary (GXRB) with a mass of 10 M_{\odot}. No extra high-frequency break/bend or any quasi-periodic oscillations are found. Through detailed analysis of energy-dependent PSDs, we found that the PSD normalization, the high-frequency PSD slope as well as the bending frequency remains consistent across all energies within the 0.3-3 keV band, revealing the presence of a constant temperature corona. These significant findings impose critical constraints on current models of X-ray emission and variability in AGN.

Super-resolution techniques have the potential to reduce the computational cost of cosmological and astrophysical simulations. This can be achieved by enabling traditional simulation methods to run at lower resolution and then efficiently computing high-resolution data corresponding to the simulated low-resolution data. In this work, we investigate the application of a Wasserstein Generative Adversarial Network (WGAN) to increase the particle resolution of dark-matter-only simulations, reproducing and building on prior results. Our WGAN models successfully generate high-resolution data with summary statistics, including the power spectrum and halo mass function, that closely match those of true high-resolution simulations. We also identify a limitation of the WGAN model in the form of smeared features in the generated high-resolution snapshots, particularly in the shapes of dark-matter halos.

We report results on the internal dynamical evolution of old star clusters located in the outer regions of the Small Magellanic Cloud (SMC). Because the SMC has been imprinted with evidence of tidal interaction with the Large Magellanic Cloud (LMC), we investigated at what extend such an interaction has produced extra tidal structures or excess of stars beyond the clusters' tidal radii. For that purpose, we used the Survey of the Magellanic Stellar History (SMASH) DR2 data sets to build number density radial profiles of suitable star clusters, and derived their structural and internal dynamics parameters. The analysed stellar density profiles do not show any evidence of tidal effects caused by the LMC. On the contrary, the Jacobi volume of the selected SMC star clusters would seem underfilled, with a clear trend toward a smaller percentage of underfilled volume as their deprojected distance to the SMC centre increases. Moreover, the internal dynamical evolution of SMC star clusters would seem to be influenced by the SMC gravitational field, being star clusters located closer to the SMC centre in a more advanced evolutionary stage. We compared the internal dynamical evolution of SMC old star clusters with those of LMC and Milky Way globular clusters, and found that Milky Way globular clusters have dynamical evolutionary paths similar to LMC/SMC old star clusters located closer to their respective galaxy's centres. Finally, we speculate with the possibility that globular clusters belonging to Magellanic Clouds like-mass galaxies have lived a couple of times their median relaxation times.

Recent observations from JWST have revealed an abundant population of active galactic nuclei (AGN) and so-called ``Little Red Dots'' (LRDs) at $2\lesssim z \lesssim 11$, many of which are characterized by V-shaped UV-to-optical continua with turnovers around the Balmer limit. The physical nature of these LRDs is unclear, and it remains debated whether the peculiar spectral shape originates from AGN, compact galaxies, or both. We present the analysis of new NIRSpec-IFU data from the BlackTHUNDER JWST Large Programme and archival NIRSpec-MSA data of a lensed LRD at $z=7.04$. The spectra confirm the presence of a smooth Balmer break and a broad H$\beta$ tracing the Broad Line Region (BLR) of an AGN. The small velocity dispersion of the H$\beta$ narrow component indicates a small dynamical mass of the host galaxy of $M_{\rm dyn}<4 \times 10^8~M_{\odot}$, which implies that the stellar population cannot contribute more than 10% to the optical continuum. We show that the Balmer break can be well described by an AGN continuum absorbed by very dense ($n_{\rm H}\sim 10^{10}~{\rm cm^{-3}}$) and nearly dust-free gas along our line-of-sight (possibly gas in the BLR or its surrounding). The same gas is expected to produce H$\beta$ absorption, at a level consistent with a tentative detection ($3\sigma$) in the high-resolution spectrum. Such a non-stellar origin of the Balmer break may apply to other LRDs, and would alleviate the issue of extremely high stellar mass surface densities inferred in the case of a stellar interpretation of the Balmer break. We note that this is a rare case of a black hole that is overmassive relative to both the host galaxy stellar and dynamical masses. We finally report indications of variability and the first attempt of AGN reverberation mapping at such an early epoch.

With the evolution of radio astronomy, related education and training, the demand for scalable, efficient, and remote systems in data acquisition, storage, and analysis has significantly increased. Addressing this need, we have developed a web interface for a log-periodic dipole antenna (LPDA) array integral to the SKA Test activities at the Gauribidanur Radio Observatory (77.428 E, 13.603 N). This interface, employing Python-based technologies such as Streamlit and PyVISA, along with SCPI commands, offers a seamless and user-friendly experience. Our solution introduces a unique data acquisition approach, employing SCPI through Python to communicate with the setup's Data Acquisition System (DAS). The web interface, accessible remotely via a secure WLAN network or VPN, facilitates user-initiated observations and comprehensive logging and offers advanced features like manual RFI masking, transit plotting, and fringe plot analysis. Additionally, it acts as a data hub, allowing for the remote downloading of observational data. These capabilities significantly enhance the user's ability to conduct detailed post-observation data analysis. The effectiveness of this interface is further demonstrated through a successful solar transit observation, validating its utility and accuracy in real-world astronomical applications. The applications of this web tool are expandable and can be tailored according to the Observatory's Goals and Instrumentation as well as for the growing radio astronomy instrumentation and observing facilities coming up at various educational institutions.

The current knowledge in cosmology deals with open problems whose solutions are still under investigation. The main issue is the so-called Hubble constant ($H_0$) tension, namely, the $4-6 \sigma$ discrepancy between the local value of $H_0$ obtained with Cepheids+Supernovae Ia (SNe Ia) and the cosmological one estimated from the observations of the Cosmic Microwave Background (CMB). For the investigation of this problem, probes that span all over the redshift $z$ ranges are needed. Cepheids are local objects, SNe Ia reached up to $z=2.9$, and CMB is observed at $z=1100$. In this context, the use of probes at intermediate redshift $z>3$ is auspicious for casting more light on modern cosmology. The Gamma-ray Bursts (GRBs) are particularly suitable for this task, given their observability up to $z=9.4$. The use of GRBs as standardizable candles requires the use of tight and reliable astrophysical correlations and the presence of gaps in the GRB time-domain data represents an obstacle in this sense. In this work, we propose to improve the precision of the lightcurve (LC) parameters through a reconstruction process performed with the functional forms of GRB LCs and the Gaussian Processes (GP). The filling of gaps in the GRB LCs through these processes shows an improvement up to $41.5\%$ on the precision of the LC parameters fitting, which lead to a reduced scatter in the astrophysical correlations and, thus, in the estimation of cosmological parameters.

The Equation of State of Quantum Chromodynamics with $N_f=3$ flavours is determined non-perturbatively in the range of temperatures between $3$ and $165$ GeV with a precision of about $0.5$-$1.0\%$. The calculation is carried out by numerical simulations of lattice gauge theory discretized à la Wilson with shifted boundary conditions in the compact direction. At each given temperature the entropy density is computed at several lattice spacings in order to extrapolate the results to the continuum limit. Taken at face value, data point straight to the Stefan-Boltzmann value by following a linear behaviour in the strong coupling constant squared. They are also compatible with the known perturbative formula supplemented by higher order terms in the coupling constant, a parameterization which describes well our data together with those present in the literature down to $500$ MeV.

According to recent lore, supercooled first-order phase transitions (FOPTs) cannot explain the evidence for a stochastic gravitational wave background obtained by pulsar timing arrays. We demonstrate that supercooled FOPTs in Majoron-like U(1)' models with a conformal dark sector easily explain the nHz signal at NANOGrav in conjunction with the observed pattern of active neutrino masses and mixing. We find that models with a U(1)' charge of the Higgs doublet close to that of the B-L model are favored.

The recent unveiling of the images of Sgr A* and M87* has significantly advanced our understanding of gravitational physics. In this study, we derive a class of Kerr-Taub-NUT metrics in the presence of a scalar field (KTNS). Treating these metrics as models for supermassive objects, we constrain the parameters using shadow size estimates done by observations of M87* and Sgr A* from the Event Horizon Telescope (EHT). Comparing the obtained results with M87* data, we show an upper limit on the NUT charge $n$ such that the constraint on the shadow deviation from circularity ($ \Delta C $) will be fulfilled for $ n<0.5 $, and this allowed range changes with a variation in other parameters. Additionally, our findings reveal that fast-rotating KTNS metrics are better candidates for supermassive M87* than slowly rotating ones. We continue our study by estimating parameters using Keck and VLTI observations of Sgr A* and find that the constraint on the fraction deviation $ \delta $ is maintained within a certain range of the NUT charge such that the Keck bound is satisfied for $ n<0.41 $. In contrast, the VLTI bound can be fulfilled for $ n>0.34 $. Finally, we investigate weak gravitational lensing using the Gauss-Bonnet theorem and illustrate that all model parameters increase the deflection angle, causing light rays to deviate more significantly near fast-rotating KTNS objects.

The RES-NOVA project is an experimental initiative aimed at detecting neutrinos from the next galactic supernova using PbWO$_{4}$ cryogenic detectors, operated at low temperatures in a low-background environment. By utilizing archaeological lead (Pb) as the target material, RES-NOVA leverages its high radiopurity, large nuclear mass, and the natural abundance of $^{207}$Pb, making it well-suited for exploring both spin-independent and spin-dependent Dark Matter (DM) interactions via nuclear scattering. This work presents a background model developed for the RES-NOVA technology demonstrator and evaluates its implications for Dark Matter detection. Detailed calculations of nuclear matrix elements, combined with the unique properties of archaeological Pb, demonstrate RES-NOVA's potential as a complementary tool to existing direct detection experiments for studying Dark Matter interactions. The experiment will conduct DM searches over a broad mass range spanning 4 orders of magnitude, from sub-GeV/$c^2$ to TeV/$c^2$. In the most optimistic scenario, RES-NOVA is expected to probe DM-nucleon cross-sections down to 1$\times 10^{-43}$ cm$^2$ and 2$\times 10^{-46}$ cm$^2$ for candidates with masses of 2 GeV/$c^2$ and 20 GeV/$c^2$, respectively.

Neutron stars (NS's) with their strong magnetic fields and hot dense cores could be powerful probes of axions, a classic benchmark of feebly-coupled new particles, through abundant production of axions with the axion-nucleon coupling and subsequent conversion into X-rays due to the axion-photon coupling. In this article, we point out that the pulsation structures in both the intensity and polarization of X-rays from NS's could provide us additional information about axions and their couplings. We develop new analytical formalisms of pulsation-polarization structure applicable to a wide range of NS's in the axion scenario and argue that they hold in complicated astrophysical environments. As a case study, we apply our formalism to a representative X-ray Dim Isolated Neutron Star, RX J1856.6-3754, with an unexpected hard X-ray excess which might be axion-induced. We show with an updated fit that the axion explanation is compatible with both the intensity and pulsation data available. Yet the preferred parameter space is close to being excluded by other astrophysical constraints. With an achievable reduction of the uncertainties in the pulsation data, we could potentially draw a definite conclusion on the axion-induced X-rays.

We present a class of relativistic fluid models for cold and dense matter with bulk viscosity, whose equilibrium equation of state is polytropic. These models reduce to Israel-Stewart theory for small values of the viscous stress $\Pi$. However, when $\Pi$ becomes comparable to the equilibrium pressure $P$, the evolution equations "adjust" to prevent the onset of far-from-equilibrium pathologies that would otherwise plague Israel-Stewart. Specifically, the equations of motion remain symmetric hyperbolic and causal at all times along any continuously differentiable flow, and across the whole thermodynamic state space. This means that, no matter how fast the fluid expands or contracts, the hydrodynamic equations are always well-behaved (away from singularities). The second law of thermodynamics is enforced exactly. Near equilibrium, these models can accommodate an arbitrarily complicated dependence of the bulk viscosity coefficient $\zeta$ on both density and temperature.

The nature of Dark Matter (DM) remains mysterious despite the substantial evidence from astrophysical and cosmological observations. While the majority of DM in our universe is non-relativistic, collisionless and its equation of state (EoS) is approximately pressureless $p\simeq 0$, DM becomes relativistic near the massive black holes in galactic center. Yet its EoS is seldom discussed in the relativistic regime. Here we initially explore the possible equation of state for DM in the vicinity of Schwarzschild black holes. We work in a spherical and quasi-static background spacetime, and describe DM as a perfect fluid in equilibrium. Through numerically solving the TOV equations with physical boundary conditions, we show that DM can have static profiles near black holes and its pressure should be negative in order to support the viable density profiles $\rho$. We illustrate with two simple general equations of state, namely the power law $p \propto \rho^\gamma$ and the radius-dependent $p \propto r\cdot \rho$, and compare them with the observations of the Milky Way. Our findings provide insights into the model-building of DM, which should incorporate the possibility of negative pressure in the relativistic regime around black holes.

Astrophysical neutrinos with energy in the TeV-PeV range traverse megaparsecs (Mpc) to gigaparsecs (Gpc) scale distances before they reach the Earth. Tiny physics effects that get accumulated over these large propagation paths during their journey may become observable at the detector. If there is some new interaction between neutrinos and the background matter, that can potentially affect the propagation of the astrophysical neutrinos. One such possible case is the flavor-dependent long-range interaction of neutrinos, which can affect the standard neutrino flavor transition, modifying the flavor composition of the astrophysical neutrinos at Earth. Using the present-day and future projection of the flavor-composition measurements of IceCube and IceCube-Gen2 along with the present and future measurement of the oscillation parameters, we explore the sensitivity of these experiments to probe long-range neutrino interaction with matter.

GRTresna is a multigrid solver designed to solve the constraint equations for the initial data required in numerical relativity simulations. In particular it is focussed on scenarios with fundamental fields around black holes and inhomogeneous cosmological spacetimes. The following overview has been prepared as part of the submission of the code to the Journal of Open Source Software. The code is based on the formalism in Aurrekoetxea, Clough \& Lim arXiv:2207.03125 and can be found at this https URL