Papers for Monday, Jun 24 2024

Using deep imaging from the CANDELS and HFF surveys, we present bulge+disc decompositions with GalfitM for $\sim$17,000 galaxies over $0.2 \leq z\leq 1.5$. We use various model parameters to select reliable samples of discs and bulges, and derive their stellar masses using an empirically calibrated relation between mass-to-light ratio and colour. Across our entire redshift range, we show that discs follow stellar mass-size relations that are consistent with those of star-forming galaxies, suggesting that discs primarily evolve via star formation. In contrast, the stellar mass-size relations of bulges are mass-independent. Our novel dataset further enables us to separate components into star-forming and quiescent based on their specific star formation rates. We find that both star-forming discs and star-forming bulges lie on stellar mass-size relations that are similar to those of star-forming galaxies, while quiescent discs are typically smaller than star-forming discs and lie on steeper relations, implying distinct evolutionary mechanisms. Similar to quiescent galaxies, quiescent bulges show a flattening in the stellar mass-size relation at $\sim$10$^{10}$M$_\odot$, below which they show little mass dependence. However, their best-fitting relations have lower normalisations, indicating that at a given mass, bulges are smaller than quiescent galaxies. Finally, we obtain rest-frame colours for individual components, showing that bulges typically have redder colours than discs, as expected. We visually derive UVJ criteria to separate star-forming and quiescent components and show that this separation agrees well with component colour. HFF bulge+disc decomposition catalogues used for these analyses are publicly released with this paper.

Context. The recent detection of nitrogen-enhanced, metal-poor galaxies at high redshift by the James Webb Space Telescope has sparked renewed interest in exploring the chemical evolution of carbon, nitrogen, and oxygen (the CNO elements) at early times, prompting fresh inquiries into their origins. Aims. The main goal of this paper is to shed light onto the early evolution of the main CNO isotopes in our Galaxy and in young distant systems, such as GN-z11 at z=10.6. Methods. To this aim, we incorporate a stochastic star-formation component into a chemical evolution model calibrated with high-quality Milky Way (MW) data, focusing on the contribution of Population III (Pop III) stars to the early chemical enrichment. Results. By comparing the model predictions with CNO abundance measurements from high-resolution spectroscopy of an homogeneous sample of Galactic halo stars, we first demonstrate that the scatter observed in the metallicity range -4.5 < [Fe/H] <-1.5 can be explained by pre-enrichment from Pop III stars that explode as supernovae (SNe) with different initial masses and energies. Then, by exploiting the chemical evolution model, we provide testable predictions for log(C/N), log(N/O), and log(C/O) vs. log(O/H)+12 in MW-like galaxies observed at different cosmic epochs/redshifts. Finally, by calibrating the chemical evolution model to replicate the observed properties of GN-z11, we provide an alternative interpretation of its log(N/O) abundance ratio, demonstrating that its high N content can be reproduced through enrichment from high-mass faint Pop III SNe. Conclusions. Stochastic chemical enrichment from primordial stars explains both the observed scatter in CNO abundances in MW halo stars and the exceptionally high N/O ratios in some distant galaxies. These findings emphasize the critical role of Pop III stars in shaping early chemical evolution.

When a structure displays dependence on distance and azimuthal angle from a center (for example the spiral arms of galaxies or the diffraction spikes of stars), projecting the pixels to polar coordinates greatly simplifies their study. This projection from one pixel grid to another is known as a "polar plot". For this purpose, a new option has been added to the GNU Astronomy Utilities (Gnuastro) in version 0.23 to "astscript-radial-profile" script, which we describe in this research note. The figures of this research note are reproducible with Maneage, on the Git commit 5d34243.

During reionization, intergalactic ionization fronts (I-fronts) are sources of Ly$\alpha$ line radiation produced by collisional excitation of hydrogen atoms within the fronts. In principle, detecting this emission could provide direct evidence for a reionizing intergalactic medium (IGM). In this paper, we use a suite of high-resolution one-dimensional radiative transfer simulations run on cosmological density fields to quantify the parameter space of I-front Ly$\alpha$ emission. We find that the Ly$\alpha$ production efficiency -- the ratio of emitted Ly$\alpha$ flux to incident ionizing flux driving the front -- depends mainly on the I-front speed and the spectral index of the ionizing radiation. IGM density fluctuations on scales smaller than the typical I-front width produce scatter in the efficiency, but they do not significantly boost its mean value. The Ly$\alpha$ flux emitted by an I-front is largest if 3 conditions are met simultaneously: (1) the incident ionizing flux is large; (2) the incident spectrum is hard, consisting of more energetic photons; (3) the I-front is traveling through a cosmological over-density, which causes it to propagate more slowly. We present a convenient parameterization of the efficiency in terms of I-front speed and incident spectral index. We make these results publicly available as an interpolation table and we provide a simple fitting function for a representative ionizing background spectrum. Our results can be applied as a sub-grid model for I-front Ly$\alpha$ emissions in reionization simulations with spatial and/or temporal resolutions too coarse to resolve I-front structure. In a companion paper, we use our results to explore the possibility of directly imaging Ly$\alpha$ emission around neutral islands during the last phases of reionization.

We present a stellar dynamical mass measurement of the supermassive black hole in the elliptical (E1) galaxy NGC 3258. Our findings are based on Integral Field Unit spectroscopy from the Multi Unit Spectroscopic Explorer observations in Narrow Field Mode with adaptive optics and in the MUSE Wide Field Mode, from which we extract kinematic information by fitting the Ca II and Mg $b$ triplets, respectively. Using axisymmetric, three-integral Schwarzschild orbit-library models, we fit the observed line-of-sight velocity distributions to infer the supermassive black hole mass, the $H$-band mass-to-light ratio, the asymptotic circular velocity, and the dark matter halo scale radius of the galaxy. We report a black hole mass of $(2.2 \pm 0.2)\times10^9 \ \rm M_{\scriptscriptstyle\odot}$ at an assumed distance of $31.9 \ \rm Mpc$. This value is in close agreement with a previous measurement from ALMA CO observations. The consistency between these two measurements provides strong support for both the gas dynamical and stellar dynamical methods.

Long troughs observed in the $z > 5.5$ Ly$\alpha$ and Ly$\beta$ forests are thought to be caused by the last remaining neutral patches during the end phases of reionization -- termed neutral islands. If this is true, then the longest troughs mark locations where we are most likely to observe the reionizing intergalactic medium (IGM). A key feature of the neutral islands is that they are bounded by ionization fronts (I-fronts) which emit Lyman series lines. In this paper, we explore the possibility of directly imaging the outline of neutral islands with a narrowband survey targeting Ly$\alpha$. In a companion paper, we quantified the intensity of I-front Ly$\alpha$ emissions during reionization and its dependence on the spectrum of incident ionizing radiation and I-front speed. Here we apply those results to reionization simulations to model the emissions from neutral islands. We find that neutral islands would appear as diffuse structures that are tens of comoving Mpc across, with surface brightnesses in the range $\approx 1 - 5\times 10^{-21}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. The islands are brighter if the spectrum of ionizing radiation driving the I-fronts is harder, and/or if the I-fronts are moving faster. We develop mock observations for current and futuristic observatories and find that, while extremely challenging, detecting neutral islands is potentially within reach of an ambitious observing program with wide-field narrowband imaging. Our results demonstrate the potentially high impact of low-surface brightness observations for studying reionization.

Crescent-shaped asymmetries are common in millimetre observations of protoplanetary discs and are usually attributed to vortices or dust overdensities. However, they often appear on a single side of the major axis and roughly symmetric about the minor axis, suggesting a geometric origin. In this work, we interpret such asymmetries as emission from the exposed inner cavity walls of inclined discs and use them to characterise their vertical extent. Here we focus on the discs around CIDA 9 and RY Tau, first modelling their observations in visibility space with a simple geometric prescription for the walls, and then exploring more detailed radiative transfer models. Accounting for the wall emission yields significantly better residuals than purely axisymmetric models, and we estimate the dust scale height of these systems to be 0.4 au at 37 au for CIDA 9 and 0.2 au at 12 au for RY Tau. Finally, we identify crescent-shaped asymmetries in twelve discs, nine of which have constraints on their orientation - in all cases, the asymmetry appears on the far-side of the disc, lending support to the hypothesis that they are due to their inner rims. Modelling this effect in larger samples of discs will help to build a statistical view of their vertical structure.

The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution gratings spectrum indicate an inflow balanced by an outflow ($s \sim 1$) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following $1/r$. The data require that the abundance of Iron relative to solar is $Z_{Fe} < 0.32 Z_\odot$ (90\% credible interval), much lower than the $2~Z_\odot$ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra's unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic Center in X-rays.

Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulation runtime of 3500 years. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe K alpha Sgr A* spectrum from Chandra HETG-S, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 1.5" of Sgr A*. We compare our results to the Radiatively Inefficient Accretion Flow (RIAF) model fit to the HETG-S spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9" angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*'s accretion flow in X-rays.

We report for the first time the detection of thermal free-free emission from post-flare loops at 34GHz in images from the Nobeyama Radioheliograph (NoRH). We studied 8 loops, 7 of which were from regions with extremely strong coronal magnetic field reported by Fedenev et al. (2023). Loop emission was observed in a wide range of wavelength bands, up to soft X-rays, confirming their multi-temperature structure and was associated with noise storm emission in metric wavelengths. The comparison of the 17GHz emission with that at 34GHz, after a calibration correction of the latter, showed that the emission was optically thin at both frequencies. We describe the structure and evolution of the loops and we computed their density, obtaining values for the top of the loops between 1 and 6 x 10^10 cm^-3, noticeably varying from one loop to another and in the course of the evolution of the same loop system; these values have only a weak dependence on the assumed temperature, 2 x 10^6 K in our case, as we are in the optically thin regime. Our density values are above those reported from EUV observations, which go up to about 10^10 cm^-3. This difference could be due to the fact that different emitting regions are sampled in the two domains and/or due to the more accurate diagnostics in the radio range, which do not suffer from inherent uncertainties arising from abundances and non-LTE excitation/ionization equilibria. We also estimated the magnetic field in the loop tops to be in the range of 10 to 30G.

The ring-like images of the two supermassive black holes captured by the Event Horizon Telescope (EHT) provide powerful probes of the physics of accretion flows at horizon scales. Specifically, the brightness asymmetry in the images carries information about the angular velocity profile of the inner accretion flow and the inclination of the observer, owing to the Doppler boosts photons experience at their site of emission. In this paper, we develop a method for quantifying the brightness asymmetry of black-hole images in the Fourier domain, which can be measured directly from interferometric data. We apply this method to current EHT data and find that the image of Sagittarius A* (Sgr A*) has an unusually low degree of asymmetry that is even lower than that inferred for M87. We then use a covariant semi-analytic model to obtain constraints on the inclinations and velocity profiles of the inner accretion flow for Sgr A*. We find that the lack of significant brightness asymmetry forces the observer inclination to uncomfortably small values ($6-10^\circ$), if the plasma velocity follows Keplerian profiles. Alternatively, larger inclination angles can be accommodated if the plasma velocities are significantly sub-Keplerian and the black hole is not spinning rapidly.

JWST has recently revealed a large population of accreting black holes (BHs) in the early Universe. Even after accounting for possible systematic biases, the high-z $M_*-M_{\rm \rm bh}$ relation derived from these objects by Pacucci et al. (2023 P23 relation) is above the local scaling relation by $>3\sigma$. To understand the implications of potentially overmassive high-z BH populations, we study the BH growth at $z\sim4-7$ using the $[18~\mathrm{Mpc}]^3$ BRAHMA suite of cosmological simulations with systematic variations of heavy seed models that emulate direct collapse black hole (DCBH) formation. In our least restrictive seed model, we place $\sim10^5~M_{\odot}$ seeds in halos with sufficient dense and metal-poor gas. To model conditions for direct collapse, we impose additional criteria based on a minimum Lyman Werner flux (LW flux $=10~J_{21}$), maximum gas spin, and an environmental richness criterion. The high-z BH growth in our simulations is merger dominated, with a relatively small contribution from gas accretion. For the most restrictive simulation that includes all the above seeding criteria for DCBH formation, the high-z $M_*-M_{\rm bh}$ relation falls significantly below the P23 relation (by factor of $\sim10$ at $z\sim4$). Only by excluding the spin and environment based criteria, and by assuming $\lesssim750~\mathrm{Myr}$ delay times between host galaxy mergers and subsequent BH mergers, are we able to reproduce the P23 relation. Overall, our results suggest that if high-z BHs are indeed systematically overmassive, assembling them would require more efficient heavy seeding channels, higher initial seed masses, additional contributions from lighter seeds to BH mergers, and / or more efficient modes for BH accretion.

Context. The scale height of the spatial distribution of open clusters (OCs) in the Milky Way exhibits a well known increase with age which is usually interpreted as evidence for dynamical heating of the disc or of the disc having been thicker in the past. Aims. We address the increase of the scale height with age of the OC population from a different angle. We propose that the apparent thickening of the disc can be largely explained as a consequence of a stronger disruption of OCs near the Galactic plane by disc phenomena, namely encounters with giant molecular clouds (GMCs). Methods. We present a computational model that forms OCs with different initial masses and follows their orbits while subjecting them to different disruption mechanisms. To setup the model and infer its parameters, we use and analyse a Gaia-based OC catalogue (Dias et al. 2021). We investigate both the spatial and age distributions of the OC population and discuss the completeness of the sample. The simulation results are then compared to the observations. Results. Consistent with previous studies, the observations reveal that the SH of the spatial distribution of OCs increases with age. We find that it is very likely that the OC sample is incomplete even for the solar neighbourhood. The model simulations successfully reproduce the SH increase with age and the total number of OCs that survive with age up to 1 Gyr. For older OCs, the predicted SH from the model starts deviating from the observations, although remaining within the uncertainties of the observations. This can be related with effects of incompleteness and/or simplifications in the model. Conclusions. A selective disruption of OCs near the galactic plane through GMC encounters is able to explain the SH evolution of the OC population.

We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions.

We perform an extensive study of equation of state (EoS) models featuring a phase transition from hadronic to deconfined quark matter in neutron star merger simulations. We employ three different hadronic EoSs, a constant speed of sound parameterization for the quark phase and a Maxwell construction to generate a large sample of hybrid EoS models. We systematically vary the onset density and density jump of the phase transition as well as the quark matter stiffness and simulate binary neutron star mergers to infer how the properties of the phase transition affect the gravitational-wave signal. In total we simulate mergers with 245 different hybrid EoS models. In particular, we explore in which scenarios a phase transition would be detectable by a characteristically increased postmerger gravitational-wave frequency compared to an estimate from the inspiral signal assuming a purely hadronic EoS. We find that the density jump at the transition (latent heat) has the largest impact on the gravitational-wave frequencies, while the influence of the stiffness of quark matter is smaller. We quantify which range of phase transition properties would be compatible with a certain magnitude or absence of the gravitational-wave postmerger frequency shift. By means of these dependencies, a future detection will thus directly yield constraints on the allowed features of the hadron-quark phase transition.

Compact groups (CGs) of galaxies are extreme environments for morphological transformations and the cessation of star formation. Our objective is to understand the dynamics of CGs and how their surrounding environment impacts galaxy properties. We selected a sample of 340 CGs in the Stripe 82 region, totaling 1083 galaxies, and a control sample of 2281 field galaxies. We find that at least 27\% of our sample of CGs are part of major structures, i.e. non-isolated CGs. We find a bimodality in the effective radius ($R_e$)-Sérsic index ($n$) plane for all transition galaxies (those with $(u-r) > 2.3$ and $n<2.5$) in CGs. Additionally, transition galaxies in isolated CGs populate more densely the $R_e-n$ plane for $n < 1.75$. In contrast, transition galaxies in non-isolated CGs have smoothly increasing $n$ values, suggesting these galaxies have already suffered morphological transformation, and primarily contribute to the distribution of more compact galaxies in the $R_e-n$ plane for all transition galaxies in CGs. We also find significant differences in the specific star-formation rate (sSFR) distribution between the late-type galaxies (LTGs) ($(u-r)<2.3$ and $n< 2.5$) in non-isolated CGs and the same type of galaxies in the control sample, suggesting that the evolution of LTGs differs in non-isolated CGs. Early-type galaxies ($(u-r)>2.3$ and $n>2.5$) and transition galaxies in non-isolated CGs have lower sSFR values and a higher fraction of quenched galaxies, compared to those in isolated CGs. Based on our results, we propose an evolutionary scenario where the major structures in which the CGs are embedded accelerate the morphological transformations of their members. Our findings highlight the importance of considering the larger structures in which CGs may be located, when analysing the properties of their galaxy, as this can significantly affect the evolution of CGs and their galaxies.

In this work, we present a new mosaic created with the second release of LOFAR Two-Metre Sky Survey data (LoTSS-DR2), which probes polarised synchrotron emission in the high-latitude inner Galaxy. Our objective is to characterise the observed emission through multi-tracer analysis to better understand the volume and the structures that may be observed with LOFAR. Furthermore, we exploit Faraday depth as a unique tool to probe the diffuse magnetised structure in the local ISM. We produced a mosaic Faraday cube of LoTSS-DR2 data by applying a rotation measure synthesis algorithm. From the cube, we constructed Faraday moment maps to characterise the nature of spectra. Additionally, we quantified the linear depolarisation canals using the Rolling Hough transform and used them to search for alignment with other data sets. Utilising LoTSS-DR2 observations alongside complementary data sets including Planck polarisation data, HI emission maps, and starlight polarisation measurements, we estimated the distance to the Faraday structures. The Faraday cube reveals an ordered structure across two-thirds of the observed area, whose orientation aligns well with that of both the HI filaments and the magnetic field. We estimate the minimum distance to the Faraday structures to be between 40 and 80 pc, which puts them in the vicinity of the Local Bubble wall. The emission is organised in a large gradient in Faraday depth whose origin we associate with the curved wall of the Local Bubble. Comparing our data with a model of the Local Bubble wall, we conclude that we might be probing a contribution of the medium inside the Local Bubble cavity as well, corresponding to the complex of local interstellar clouds. Moreover, we propose a toy model incorporating an ionised front of finite thickness into the Local Bubble wall, as a curved, cold neutral shell alone is insufficient to produce the observed gradient.

The chemical makeup of our solar system is a reflection of Galactic chemical evolution in the local interstellar medium (ISM) over the past ~9 Ga before the formation of the solar system. Although the incorporated ISM dust was mostly destroyed during the solar system formation, a small fraction of the ISM dust, known as presolar grains, is preserved in pristine extraterrestrial materials and identified through their exotic isotopic compositions, pointing to their formation in gas outflows or explosions of ancient stars. Since their stellar birth at more than 4.6 Ga, presolar grains have borne witness to a diverse array of astrophysical and cosmochemical processes. In this chapter, I will review recent progress in utilizing the isotopic and structural compositions of presolar grains to constrain physical mixing processes and dust formation in stars, stellar nucleosynthesis, ISM processes, and the origin and evolution of the solar system.

We present results from the occultation of Betelgeuse by asteroid (319) Leona on December 12, 2023, observed using a 64x64 pixel Single-Photon Avalanche Diode (SPAD) array mounted on a 10-inch telescope at the AstroCamp Observatory in Nerpio, Southeast of Spain, just a few kilometers from the center of the occultation shadow path. This study highlights remarkable advancements in applying SPAD technology in astronomy. The SPAD array's asynchronous readout capacity and photon-counting timestamp mode enabled a temporal resolution of 1 microsecond in our light curve observations of Betelgeuse. Our data analysis addressed challenges inherent to SPAD arrays, such as optical cross-talk and afterpulses, which typically cause the photon statistics to deviate from a Poisson distribution. By adopting a generalized negative binomial distribution for photon statistics, we accurately describe the observational data. This method yielded an optical cross-talk estimation of 1.07% in our SPAD array and confirmed a negligible impact of spurious detected events due to afterpulses. The meticulous statistical examination of photon data underscores our SPAD-array's exceptional performance in conducting precise astronomical observations. The observations revealed a major decrease in Betelgeuse's intensity by 77.78% at the occultation's peak, allowing us to measure Betelgeuse's angular diameter at 57.26 mas in the SDSS g-band. This measurement, employing a simplified occultation model and considering the known properties of Leona, demonstrates the potential of SPAD technology for astronomy and sets a new standard for observing ultra-rapid transient celestial events, providing a valuable public dataset for the astronomical community.

Galaxy clusters are useful cosmological probes and interesting astrophysical laboratories. With growing cluster samples, a deeper understanding of the sample characteristics and improved control of systematics becomes more crucial. In this analysis we create a new and larger ACT-DR5-based thermal Sunyaev-Zel'dovich Effect (tSZE) selected galaxy cluster catalog with improved control over sample purity and completeness. We employ the red sequence based cluster redshift and confirmation tool MCMF together with optical imaging data from the Legacy Survey DR-10 and infrared data from the WISE satellite to systematicallyidentify true clusters from a new cluster candidate detection run on the ACT-DR5 dataset. The resulting ACT-DR5 MCMF sample contains 6,237 clusters with a residual contamination of 10.7%. This is an increase of 51% compared to the previous ACT-DR5 cluster catalog, making this catalog the largest tSZE-selected cluster catalog to date. The z>1 subsample contains 703 clusters, three times more than in the previous ACT-DR5 catalog. Matching the ACT-DR5 MCMF cluster catalog with a deeper tSZE sample from SPTpol 500d allows us to confirm the completeness and purity of the new ACT-DR5 MCMF sample. Cross-matching to the two largest X-ray selected cluster samples, the all-sky RASS MCMF and the half-sky eRASS1, confirms the sample purity of the RASS MCMF sample and in the case of eRASS1 reveals that 43% of the matched clusters are designated in eRASS1 as X-ray point sources rather than clusters. Cross-correlating the ACT-DR5 MCMF cluster catalog with ACT-DR6 lensing maps results in a 16.4\sigma detection of CMB lensing around the clusters, corresponding to the strongest signal found so far for a galaxy cluster sample. Repeating the measurement for the z>1 cluster subsample yields a significance of 4.3\sigma, which is the strongest CMB lensing detection in a z>1 cluster sample to date.

Nancy Grace Roman Space Telescope will revolutionize our understanding of the Galactic Bulge with its Galactic Bulge Time Domain survey. At the same time, Rubin Observatories's Legacy Survey of Space and Time (LSST) will monitor billions of stars in the Milky Way. The proposed Roman survey of the Galactic Plane, with its NIR passbands and exquisite spacial resolution, promises groundbreaking insights for a wide range of time-domain galactic astrophysics. In this white paper, we describe the scientific returns possible from the combination of the Roman Galactic Plane Survey with the data from LSST.

Ultrahigh-energy cosmic rays are often characterized indirectly by analyzing the properties of secondary cosmic ray particles produced in the collisions with air nuclei, the particle number $N_\mu$ of muon and the depth of shower maximum $X_\mathrm{max}$ after air shower cascade are mostly studied to infer the energy and mass of the incident cosmic rays. While, researches have shown that there is a significant excess comparing the observed number of muons arriving at the ground from extensive air showers (EAS) with the simulations by using the existing cosmic ray hadronic interaction model. To explain this muon excess phenomenon, a new theoretical model, the gluon condensation model (GC model), is introduced in this paper and simulated by using the AIRES engine. We asumme that the GC effect appearing mainly in the first colliding of the cascade, leads to a significant increase in the strangeness production, accodingly, the production rate of kaons is improved appearantly. The model assumes that only pions and kaons are the new prodcutions in the hadron cascades. It is found that, considering the GC effect, the value of $n_K/n_\pi$ increases and more energy of the incident cosmic rays paticipate in hadron cascades, and then increase the number of mouns finally. This model provides a new theoretical possibility to explain the muon excess puzzle.

Cosmography can be used to constrain the kinematics of universe in a model-independent way. In this work, we attempted to combine the Pad$\rm \acute{e}$ approximations with the latest Pantheon+ sample for testing cosmological principle. Based on the Pad$\rm \acute{e}$ approximations, we first gave the cosmographic constraints on the different order polynomials including third-order (Pad$\rm \acute{e}$$_{(2,1)}$), fourth-order (Pad$\rm \acute{e}$$_{(2,2)}$) and fifth-order (Pad$\rm \acute{e}$$_{(3,2)}$). Based on the Pad$\rm \acute{e}$$_{(2,1)}$ ($j_{0}$ = 1) polynomial and hemisphere comparison (HC) method, we tested the cosmological principle and found the preferred directions of cosmic anisotropy, such as (l, b) = (304.6$^{\circ}$$_{-37.4}^{+51.4}$, $-$18.7$^{\circ}$$_{-20.3}^{+14.7}$) and (311.1$^{\circ}$$_{-8.4}^{+17.4}$, $-$17.53$^{\circ}$$_{-7.7}^{+7.8}$) for $q_{0}$ and $H_{0}$, respectively. These two directions are consistent with each other in $1\sigma$ confidence level, but the corresponding results of statistical isotropy analyses including Isotropy and Isotropy with real positions (RP) are quite different. The statistical significance of $H_{0}$ are stronger than that of $q_{0}$, i.e., 4.75$\sigma$ and 4.39$\sigma$ for the Isotropy and Isotropy with RP respectively. Reanalysis with fixed $q_{0} = -0.55$ (corresponds to $\Omega_{m}$ = 0.30) gives similar results. Overall, our model-independent results provide clear indications for a possible cosmic anisotropy, which must be taken seriously. Further test is needed to better understand this signal.

The nature of the interaction between active galactic nuclei (AGNs) and their host galaxies remains an unsolved question. Therefore, conducting an AGN census is valuable to AGN research. Nevertheless, a significant fraction of AGNs are obscured by their environment, which blocks UV and optical emissions due to the dusty torus surrounding the central supermassive black hole (SMBH). To overcome this challenge, mid-infrared (IR) surveys have emerged as a valuable tool for identifying obscured AGNs, as the obscured light is re-emitted in this range. With its high sensitivity, the James Webb Space Telescope (JWST) uncovered more fainter objects than previous telescopes. By applying the SED fitting, this work investigates AGN candidates in JWST Cosmic Evolution Early Release Science (CEERS) fields. We identified 42 candidates, 30 of them are classified as composites ($0.2\leq f_{\rm AGN, IR}< 0.5$), and 12 of them are AGNs ($f_{\rm AGN, IR}\geq 0.5$). We report the AGN luminosity contributions and AGN number fractions as a function of redshift and total infrared luminosity, showing that previously reported increasing relations are not apparent in our sample due to the sample size. We also extend the previous results on ultra-luminous infrared galaxies (ULIRGs, $L_{\rm TIR}\geq 10^{12} L_{\odot}$) to less luminous AGNs, highlighting the power of JWST.

From the LOFAR Two-meter Sky Survey second data release (LoTSS DR2) at 144 MHz, we identified a peculiar radio galaxy, J0011+3217. It has a large, one-sided diffuse secondary wing that stretches up to 0.85 Mpc (roughly 85\% of the primary lobe). The linear size of the primary lobe of the galaxy is 0.99 Mpc. This peculiar source is a giant radio galaxy with a misaligned primary lobe. An optical galaxy resides at 16 kpc (7 arcsec) linear distance from the host active galactic nucleus (AGN) of J0011+3217. J0011+3217 has a radio luminosity of $1.65\times 10^{26}$ W Hz$^{-1}$ at 144 MHz with a spectral index of $-0.80$ between 144 and 607 MHz. J0011+3217 is located 1.2 Mpc away from the center of the Abell 7 cluster. The Abell 7 cluster has a redshift of 0.104 and a mass ($M_{500}$) of 3.71 $\times 10^{14}$ M$_\odot$. The cluster is associated with strong X-ray emission. To study the X-ray emission around the cluster and from the region surrounding J0011+3217, we present an XMM-Newton image of J0011+3217, and we analyse the velocity structure and spatial distribution of galaxies in the cluster, showing that J0011+3217 inhabits an offset group of galaxies that are moving with respect to Abell 7. The off-axis distortion or bending of the primary lobe of J0011+3217 in the outer edges has a strong effect on the relative motion of the surrounding medium that causes the bending of the jets in the opposite direction of the cluster (like WAT sources). The current paper suggests that the morphology of J0011+3217 is influenced by ram pressure created due to the Abell 7 cluster, highlighting the complex interactions between the source and the surrounding cluster environment.

The Fred Young Submillimeter Telescope (FYST), on Cerro Chajnantor in the Atacama desert of Chile, will conduct wide-field and small deep-field surveys of the sky with more than 100,000 detectors on the Prime-Cam instrument. Kinetic inductance detectors (KIDs) were chosen as the primary sensor technology for their high density focal plane packing. Additionally, they benefit from low cost, ease of fabrication, and simplified cryogenic readout, which are all beneficial for successful deployment at scale. The cryogenic multiplexing complexity is pulled out of the cryostat and is instead pushed into the digital signal processing of the room temperature electronics. Using the Xilinx Radio Frequency System on a Chip (RFSoC), a highly multiplexed KID readout was developed for the first light Prime-Cam and commissioning Mod-Cam instruments. We report on the performance of the RFSoC-based readout with multiple detector arrays in various cryogenic setups. Specifically we demonstrate detector noise limited performance of the RFSoC-based readout under the expected optical loading conditions.

We explore the spatial distribution of organics on Ceres using the visible and near-infrared data collected by the Dawn mission. We employ a spectral mixture analysis (SMA) approach to map organic materials within the Ernutet crater at the highest available spatial resolution revealing a discontinuous, granular distribution and a possible correlation with an ancient crater on which Ernutet has been superimposed. The SMA technique also helps us identify 11 new areas as potential sites for organics. These regions are predominantly located within craters or along their walls, resembling the distribution pattern observed in Ernutet, which implies a possible geological link with materials exposed from beneath the surface. In one of these candidate regions situated in the Yalode quadrangle, we detected the characteristic 3.4-micron absorption band in the infrared spectrum, indicative of organics and carbonates. By combining the spatial resolution of the Framing Camera data with the spectral resolution of the Visual and Infrared Imaging Spectrometer using SMA, we investigated the distribution of the 3.4-micron band in this quadrangle. The absorption pattern correlates with the Yalode/Urvara smooth material unit, which formed after significant impacts on Ceres. The association of organic-rich materials with complex and multiple large-impact events supports for an endogenous origin for the organics on Ceres.

The variability of helium abundance in the solar corona and the solar wind is an important signature of solar activity, solar cycle, solar wind sources, as well as coronal heating processes. Motivated by recently reported remote sensing UV imaging observations by Helium Resonance Scattering in the Corona and Heliosphere (HERSCHEL) payload sounding rocket of helium abundance in inner corona on 14-Sep-2009 near solar minimum, we present the results of the first three-dimensional three-fluid (electrons, protons, alpha particles) model of tilted coronal streamer belt and slow solar wind that illustrate the various processes leading to helium abundance differentiation and variability. We find good qualitative agreement between the three-fluid model and the coronal helium abundances variability reported from UV observations of streamers, providing insight on the effects of the physical processes, such as heating, gravitational settling and inter-species Coulomb friction in the out-flowing solar wind that produce the observed features. The study impacts our understanding of the origins of the slow solar wind.

Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW) during the inspiral of a binary NS merger. To this aim, it is important to estimate the theoretical uncertainties on the EoS, one of which is the possible bias coming from an inconsistent treatment of the low-density region; that is, the use of a so called non-unified NS crust. We propose a numerical tool allowing the user to consistently match a nuclear-physics informed crust to an arbitrary high-density EoS describing the core of the star. We introduce an inversion procedure of the EoS close to saturation density that allows users to extract nuclear-matter parameters and extend the EoS to lower densities in a consistent way. For the treatment of inhomogeneous matter in the crust, a standard approach based on the compressible liquid-drop (CLD) model approach was used in our work. A Bayesian analysis using a parametric agnostic EoS representation in the high-density region is also presented in order to quantify the uncertainties induced by an inconsistent treatment of the crust. We show that the use of a fixed, realistic-but-inconsistent model for the crust causes small but avoidable errors in the estimation of global NS properties and leads to an underestimation of the uncertainties in the inference of NS properties. Our results highlight the importance of employing a consistent EoS in inference schemes. The numerical tool that we developed to reconstruct such a thermodynamically consistent EoS, CUTER, has been tested and validated for use by the astrophysical community.

There are serious prospects that the next $\gamma$-ray space mission could use a telescope based on silicon pixel detectors. I characterize the potential of such active targets for polarimetry with gamma-ray conversions to pairs and find it excellent both in terms of selection efficiency and of effective polarization asymmetry.

On the second of November 2019 the black hole X-ray binary MAXI J0637$-$430 went into outburst, at the start of which it was observed in a thermal ``disc-dominated'' state. High photon energy (extending above 10 keV) observations taken by the NuSTAR telescope reveal that this thermal spectrum can not be fit by conventional two-component (disc plus corona) approaches which ignore disc emission sourced from within the plunging region of the black hole's spacetime. Instead, these models require a third ``additional'' thermal component to reproduce the data. Using new disc solutions which extend classical models into the plunging region we show that this ``additional'' thermal emission can be explained self-consistently with photons emitted from the accretion flow at radii within the innermost stable circular orbit of the black hole. This represents the second low mass X-ray binary, after MAXI J1820+070, with a detection of plunging region emission, suggesting that signatures of this highly relativistic region may well be widespread but not previously widely appreciated. To allow for a detection of the plunging region, the black hole in MAXI J0637$-$430 must be at most moderately spinning, and we constrain the spin to be $a_\bullet < 0.86$ at 99.9$\%$ confidence. We finish by discussing the observational requirements for the robust detection of this region.

Heavy nuclei can be synthetized or entrained in gamma-ray bursts (GRBs) with implications on the high-energy neutrino emission. By means of a Monte-Carlo algorithm, we model nuclear cascades and investigate their impact on the neutrino production considering kinetic dominated jets (in the internal shock model, including a dissipative photosphere) as well as Poynting flux dominated jets (for a jet model invoking internal-collision-induced magnetic reconnection and turbulence, ICMART). We find that the ICMART model allows for efficient nuclear cascades leading to an overall larger neutrino fluence than in the other two jet models. The survival of nuclei and inefficient nuclear cascades lead to an overall reduction of the neutrino fluence up to one order of magnitude. However, if nuclei are disintegrated, the neutrino fluence may be comparable to the one emitted from a jet loaded with protons. Exploring the parameter space of jet properties, we conclude that the composition and the bulk Lorentz factor have significant impact on the efficiency of nuclear cascades as well as the spectral shape of the expected neutrino fluence. On the other hand, the neutrino spectral distribution is less sensitive to the power-law index of the accelerated population of protons or heavier nuclei. For what concerns the diffuse emission of neutrinos from GRBs, we find that the uncertainty due to the jet composition can be at most comparable to the one related to the GRB cosmological rate.

Aims. The origin of the hard X-ray emission in the Be/X-ray binary system X Persei has long been debated as its atypical 'two-hump' spectrum can be modelled in multiple ways. The main debate focuses on the the high-energy hump, which is fit as either a cyclotron resonance scatter frequency (CRSF) or inverse Comptonization due to bulk Comptonization. Methods. Using INTEGRAL/JEM-X and ISGRI data, we studied the temporal and spectral variability in the 3-250 keV energy range during observations over ~15 years. A NuSTAR observation was also included in a joint spectral fit with the INTEGRAL spectrum. Results. We find that the joint spectrum can be described well by a low-energy component due to thermal Comptonization and a high-energy component due to bulk Comptonization, a CRSF, or a cyclotron emission line. The three models begin to diverge above ~120 keV, where statistics are low. Conclusions. We compare our results with observations of other Be/X-ray binaries that show similar 'two-hump' spectra while in a low-luminosity state. As the sources are in a low accretion state, the bulk Comptonization process is likely inefficient, and thus not an explanation for the high-energy component. The broad CRSF (27 +/- 2 keV) in X Persei suggests that the high-energy emission is not due to a CRSF. Thus, the high-energy component is potentially due to cyclotron emission, though other scenarios are not definitively excluded.

There is something about inverse-square laws that makes them fail at extremities. The inverse-square law, which is used commonly in studying the irradiance on exoplanets fails at the extreme case of planets closer than 0.01 AU. Therefore, in order to correctly predict possible climate states of such systems, we need a new model which accurately calculates angles subtended by various surface elements and integrate the generalized equation. A direct consequence of such a model is the shift of the terminator about 20 degrees beyond the poles. The irradiance at the poles is about 100-200 kW/m2 higher than what is predicted by the inverse-square law. This work therefore becomes crucial because it underscores the need to modify the current GCM models. The error in the numerical integration values of irradiance is less than one percent making our estimate very reliable.

The Keck Planet Imager and Characterizer (KPIC), a series of upgrades to the Keck II Adaptive Optics System and Instrument Suite, aims to demonstrate high-resolution spectroscopy of faint exoplanets that are spatially resolved from their host stars. In this paper, we measure KPIC's sensitivity to companions as a function of separation (i.e., the contrast curve) using on-sky data collected over four years of operation. We show that KPIC is able to reach contrasts of $1.3 \times 10^{-4}$ at 90 mas and $9.2 \times 10^{-6}$ at 420 mas separation from the star, and that KPIC can reach planet-level sensitivities at angular separations within the inner working angle of coronagraphic instruments such as GPI and SPHERE. KPIC is also able to achieve more extreme contrasts than other medium-/high-resolution spectrographs that are not as optimized for high-contrast performance. We decompose the KPIC performance budget into individual noise terms and discuss limiting factors. The fringing that results from combining a high-contrast imaging system with a high-resolution spectrograph is identified as an important source of systematic noise. After mitigation and correction, KPIC is able to reach within a factor of 2 of the photon noise limit at separations < 200 mas. At large separations, KPIC is limited by the background noise performance of NIRSPEC.

The ground-based technique for imaging atmospheric Cherenkov telescopes became a rapidly developing and powerful branch of science. Thanks to this technique, over 250 very high-energy gamma-ray sources of galactic and extragalactic origin have been discovered. Many fundamental questions of astrophysics, astro-particle physics, the physics of cosmic rays and cosmology are the focus of this technique. In the past 33 years since the discovery of the first gamma-ray source, the Crab Nebula, the discipline has made remarkable progress. Today, the technology boasts highly sensitive telescopes capable of detecting a point source 100 times fainter than the standard candle, the Crab Nebula, in 25 hours of observation. Further developments in this technology led to the Cherenkov Telescope Array (CTA), the next-generation large instrument. The sensitivity of CTA will be several times higher than that of the current best instruments. This article presents a brief history of ground-based very high energy gamma-ray astrophysics.

Edge-on galaxies have many important applications in galactic astrophysics, but they can be difficult to identify in vast amounts of astronomical data. To facilitate the search for them, we have developed a deep learning algorithm designed to identify and extract edge-on galaxies from astronomical images. We utilised a sample of edge-on spiral galaxies from the Galaxy Zoo database, retrieving the corresponding images from the Sloan Digital Sky Survey (SDSS). Our dataset comprises approximately $16,000$ galaxies, which we used to train the YOLOv5 algorithm for detection purposes. To isolate galaxies from their backgrounds, we trained the SCSS-Net neural network to generate segmentation masks. As a result, our algorithm detected $8,000$ edge-on galaxies, for which we compiled a catalogue including their parameters obtained from the SDSS database. We describe the basic properties of our sample, finding that most galaxies have redshifts $0.02<z<0.10$, have low values of $b/a$ and are mostly red, which is expected from edge-on galaxies and is consistent with our training sample, as well as other literature. The cutouts of the detected galaxies can be used for future studies and the algorithm can be applied to data from future surveys as well.

Palaeomagnetic evidence shows that the behaviour of the geodynamo has changed during geological times. Variations in the heat flux at the core-mantle boundary (CMB) due to mantle convection could be responsible. Previous studies, based on unrealistically viscous dynamo simulations, have shown that large-scale CMB heat flux heterogeneities impact the magnetic dipole stability. To better understand how they affect the geodynamo, we used several simulations, ranging from standard numerical dynamos to more extreme parameters, including strong-field cases and turbulent cases. We show that heterogeneities with realistic amplitudes can favour a multipolar dynamo by either forcing equatorially antisymmetric zonal flows or eastward zonal flows. Strong-field dynamo models are found to be less sensitive, due to significant westward flows. We also find that the dipolar fraction of the magnetic field is best captured by $M^*=M\ E_{\eta}\ \dfrac{l_c}{\pi}$ where $M$ is the magnetic to kinetic energy ratio, $E_{\eta}$ is the magnetic Ekman number, and $l_c$ is the dominant harmonic degree of the flow, with multipolar dynamos found at lower $M^*$. $M^*$ estimated for the Earth's core is consistent with a reversing dipolar magnetic field. Within the range of $M^*$ susceptible to reversals, breaking the equatorial symmetry or forcing eastward zonal flows in our simulations consistently triggers reversals or a transition towards multipolar dynamos. Our results support that time variations of heat-flux heterogeneities driven by mantle convection through Earth's history are capable of inducing the significant variations in the reversal frequency observed in the palaeomagnetic record.

The high luminosity of massive, early-type stars drives strong stellar winds through line scattering of the stars continuum radiation. Their momenta contribute substantially to the dynamics and energetics of the ambient interstellar medium in galaxies. The detailed multi-wavelength study of massive O-type and Wolf-Rayet binaries is essential to explore the hydrodynamics of the shocks formed in the stellar outflows and wind structure. Further, deep analysis of some of the interesting phenomena like particle acceleration and dust formation associated with hot stars winds provides a global view of stellar outflows. In this context, a few massive binaries have been explored using photometric and spectroscopic measurements in different wavebands. This paper highlights important insights gained from investigating massive binaries with several ground and space-based facilities.

The Dragonfly Spectral Line Mapper is a mosaic telescope comprising 120 Canon telephoto lenses, based on the design of the Dragonfly Telephoto Array. With a wide field of view, and the addition of the "Dragonfly Filter-Tilter" instrumentation holding ultra narrow bandpass filters in front of each lens, the Dragonfly Spectral Line mapper is optimized for ultra low surface brightness imaging of visible wavelength line emission. The Dragonfly Spectral Line Mapper was constructed and commissioned in four phases from March 2022 to November 2023. During this time, four individual mounts of 30 lenses each were constructed and commissioned. The commissioning of the telescope included the deployment of the "Dragonfly StarChaser" which carries out image stabilization corrections in the telephoto lens, to enable hour-long exposures to be taken. In addition, we introduced new instrumentation such as a film to cover the optics to keep the filters clean. Here we describe the updated design of the complete 120-lens array, and the implementation of the instrumentation described above. Additionally, we present updated characterization of the cameras and filter transmission for the full array. Finally, we reflect on the construction and commissioning process of the complete 120-lens array Dragonfly Spectral Line Mapper, and remark on the feasibility of a larger 1000-lens array.

LHS 1140 b is the second-closest temperate transiting planet to the Earth with an equilibrium temperature low enough to support surface liquid water. At 1.730$\pm$0.025 R$_\oplus$, LHS 1140 b falls within the radius valley separating H$_2$-rich mini-Neptunes from rocky super-Earths. Recent mass and radius revisions indicate a bulk density significantly lower than expected for an Earth-like rocky interior, suggesting that LHS 1140 b could either be a mini-Neptune with a small envelope of hydrogen ($\sim$0.1% by mass) or a water world (9--19% water by mass). Atmospheric characterization through transmission spectroscopy can readily discern between these two scenarios. Here, we present two JWST/NIRISS transit observations of LHS 1140 b, one of which captures a serendipitous transit of LHS 1140 c. The combined transmission spectrum of LHS 1140 b shows a telltale spectral signature of unocculted faculae (5.8 $\sigma$), covering $\sim$20% of the visible stellar surface. Besides faculae, our spectral retrieval analysis reveals tentative evidence of residual spectral features, best-fit by Rayleigh scattering from an N$_2$-dominated atmosphere (2.3 $\sigma$), irrespective of the consideration of atmospheric hazes. We also show through Global Climate Models (GCM) that H$_2$-rich atmospheres of various compositions (100$\times$, 300$\times$, 1000$\times$solar metallicity) are ruled out to $>$10 $\sigma$. The GCM calculations predict that water clouds form below the transit photosphere, limiting their impact on transmission data. Our observations suggest that LHS 1140 b is either airless or, more likely, surrounded by an atmosphere with a high mean molecular weight. Our tentative evidence of an N$_2$-rich atmosphere provides strong motivation for future transmission spectroscopy observations of LHS 1140 b.

The Be star zeta Tau was recently reported to be a gamma Cas analog; that is, it displays an atypical (bright and hard) X-ray emission. The origin of these X-rays remains debated.The first X-ray observations indicated a very large absorption of the hot plasma component (N_H~ 10^{23}/cm^2). This is most probably related to the edge-on configuration of the zeta Tau disk. If the X-ray emission arises close to the companion, an orbital modulation of the absorption could be detected as the disk comes in and out of the line of sight. New XMM-Newton data were obtained to characterize the high-energy properties of zeta Tau in more detail. They are complemented by previous Chandra and SRG/eROSITA observations as well as by optical spectroscopy and TESS photometry. The high-quality XMM-Newton data reveal the presence of a faint soft X-ray emission, which appears in line with that recorded for non-gamma Cas Be stars. In addition, zeta Tau exhibits significant short-term variability at all energies, with larger amplitudes at lower frequencies (``red noise''), as is found in X-ray data of other gamma Cas stars. Transient variability (softness dip, low-frequency signal) may also be detected at some epochs. In addition, between X-ray exposures, large variations in the spectra are detected in the 1.5-4.keV energy band. They are due to large changes in absorption toward the hottest (9keV) plasma. These changes are not correlated with either the orbital phase or the depth of the shell absorption of the Halpha line. These observed properties are examined in the light of proposed gamma Cas models.

The Fourier Domain Acceleration Search (FDAS) and Fourier Domain Jerk Search (FDJS) are proven matched filtering techniques for detecting binary pulsar signatures in time-domain radio astronomy datasets. Next generation radio telescopes such as the SPOTLIGHT project at the GMRT produce data at rates that mandate real-time processing, as storage of the entire captured dataset for subsequent offline processing is infeasible. The computational demands of FDAS and FDJS make them challenging to implement in real-time detection pipelines, requiring costly high performance computing facilities. To address this we propose Pulscan, an unmatched filtering approach which achieves order-of-magnitude improvements in runtime performance compared to FDAS whilst being able to detect both accelerated and some jerked binary pulsars. We profile the sensitivity of Pulscan using a distribution (N = 10,955) of synthetic binary pulsars and compare its performance with FDAS and FDJS. Our implementation of Pulscan includes an OpenMP version for multicore CPU acceleration, a version for heterogeneous CPU/GPU environments such as NVIDIA Grace Hopper, and a fully optimized NVIDIA GPU implementation for integration into an AstroAccelerate pipeline, which will be deployed in the SPOTLIGHT project at the GMRT. Our results demonstrate that unmatched filtering in Pulscan can serve as an efficient data reduction step, prioritizing datasets for further analysis and focusing human and subsequent computational resources on likely binary pulsar signatures.

The stochastic nature of star formation and photon propagation in high-redshift galaxies can result in sizable galaxy-to-galaxy scatter in their properties. Ignoring this scatter by assuming mean quantities can bias estimates of their emissivity and corresponding observables. We construct a flexible, semi-empirical model, sampling scatter around the following mean relations: (i) the conditional halo mass function (CHMF); (ii) the stellar-to-halo mass relation (SHMR); (iii) galaxy star formation main sequence (SFMS); (iv) fundamental metallicity relation (FMR); (v) conditional intrinsic luminosity; and (vi) photon escape fraction. In our fiducial model, ignoring scatter in these galaxy properties overestimates the duration of the EoR, delaying its completion by up to $\Delta z$ ~ 2. We quantify the relative importance of each of the above sources of scatter in determining the ionizing, soft-band X-ray and Lyman Werner (LW) emissivities as a function of scale and redshift. We find that scatter around the SFMS is important for all bands, especially at the highest redshifts where the emissivity is dominated by the faintest, most "bursty" galaxies. Ignoring this scatter would underestimate the mean emissivity and its standard deviation computed over 5 cMpc regions by factors of up to $\sim$2-10 at $5< z < 15$. Scatter around the X-ray luminosity to star formation rate relation is important for determining X-ray emissivity, accounting for roughly half of its mean and standard deviation. The importance of scatter in the ionizing escape fraction depends on its functional form, while scatter around the SHMR contributes at the level of ~10-20%. Although scatter does flatten the UV luminosity functions, shifting the bright end by 1-2 magnitudes, the level of scatter in our fiducial model is insufficient to fully explain recent estimates from JWST photometry (consistent with previous studies).

The $\mathrm{H}_2$ ionisation rate in the Central Molecular Zone, located in the Galactic Centre region, is estimated to be $\zeta\sim2\times10^{-14}~\mathrm{s}^{-1}$, based on observations of H$_3^+$ lines. This value is 2-3 orders of magnitude larger than that measured anywhere else in the Galaxy. A high cosmic-ray density has been invoked to explain the unusually high ionisation rate. However, this excess is not seen in the $\gamma$-ray emission from this region, which is produced by high-energy cosmic rays. Therefore, an excess is expected only in the low-energy cosmic-ray spectrum. Here, we derive constraints on this hypothetical low-energy component in the cosmic-ray spectra and we question its plausibility. To do so, we solve numerically the cosmic-ray transport equation in the Central Molecular Zone, considering spatial diffusion, advection in the Galactic wind, reacceleration in the ambient turbulence, and energy losses due to interactions with matter and radiation in the interstellar medium. We derive stationary solutions under the assumption that cosmic rays are continuously injected by a source located in the Galactic Centre. The high-energy component in the cosmic-ray spectrum is then fitted to available $\gamma$-ray and radio data, and a steep low-energy component is added to the cosmic-ray spectrum to explain the large ionisation rates. We find that injection spectra of $p^{-7}$ for protons below $p_{enh,p}c\simeq780~\mathrm{MeV}$ and $p^{-5.2}$ for electrons below $p_{enh,e}c=1.5~\mathrm{GeV}$ are needed to reach the observed ionisation rates. This corresponds to an extremely large cosmic-ray power of the order $\sim10^{40-41}~\mathrm{erg}\,\mathrm{s}^{-1}$ injected at the Galactic Centre. We conclude that cosmic rays alone can not explain the high ionisation rates in the Galactic Centre region.

We have selected candidate massive stars in the galaxy IC 342 based on archival images from the Hubble Space Telescope and images from the 2 m telescope at the Natioal Astronomical Observatory Rozhen, Bulgaria. Spectral observations of 24 out of 27 selected stars are carried out with the 6 m BTA telescope at the SAO RAS and with the 3.5 m Apache Point Observatory telescope (USA) as part of the program for searching bright massive stars in galaxies outside the Local Group. Our analysis reveals that 12 objects have spectra lacking prominent features, except for the emission lines of the surrounding nebulae and are identified as single supergiants of classes O9 to F5 or spatially unresolved young compact clusters. One source with an absorption spectrum probably belongs to our Galaxy. The spectra of seven other objects show features typical of Wolf-Rayet stars or compact clusters containing Wolf-Rayet stars. Another source is a compact supernova remnant. Two other objects are tentatively classified as cold LBV candidates, and one object is classified as a B[e]-supergiant candidate.

The magnetic properties of the umbra-penumbra (UP) boundary of sunspots and the boundary of pores at various evolutionary stages have been characterised using datasets from different instruments. We aim to study the differences between the intensity and vector magnetic field properties derived from Hinode/SP and SDO/HMI observations of a decaying sunspot. We analysed the sunspot embedded in AR 12797 during six days in 30 SP/Hinode scans and 704 HMI/SDO for both regular maps and HMI_dcon. We studied the correlation of the magnetic properties and continuum intensity in the datasets within the spot, and at the UP boundary. We examined the decaying process using the full temporal resolution of the HMI_dcon maps. We find a good correspondence between the magnetic properties in the SP and HMI_dcon maps, but the continuum intensity of the spots in the SP maps is 0.04I_qs brighter than in the HMI_dcon maps. The influence of scattered light in the HMI maps makes it the least ideal dataset for studying the boundary of spots without a penumbra. The properties at the UP boundary evolve slowly during the sunspot decay stage, while the penumbra still provides some stability. In contrast, they respond more abruptly to areal changes in the naked-spot stage. During the sunspot decay, we find linear decay in the area and in the magnetic flux. Moreover, the umbra shows two characteristic decaying processes: a slow decay during the first three days, and a sudden fast decay during the final dissipation of the penumbra. We find indications of a 3.5h lag between the dissipation of the vertical fields in the umbral region and the photometric decay of the umbral area. The differences found in the continuum intensity and in the vertical component of the magnetic field, Bver, between the analysed datasets explain the discrepancies among the Bver values found at the boundaries of umbrae in previous studies.

We present Gravity.jl, a new software for the modeling of gravitational lens systems. Gravity.jl is written in the Julia programming language, and is designed to be fast, accurate, and flexible. It can be used to model gravitational lens systems composed of multiple lensing planes, and to perform Bayesian inference on the lens model parameters. In this paper we present the theoretical and statistical ideas behind the code, and we describe its main features. In this first paper of the series, we focus on the modeling of point-like and small extended sources, for which we can linearize the lens equation. We show a practical use of Gravity.jl on a galaxy-scale lens, and we compare the results with those obtained with other codes. We also show how Gravity.jl can be used to perform Bayesian inference on cosmological parameters.

We analyse photometric data of nine supernovae (SNe) in filters V, R and I obtained during observational campaigns at the OAUNI site in 2016, 2017 and 2023. The calibrated magnitudes of the observed SNe were compared with their respective light curves available in the literature to study their evolution after their maximum brightness. In some cases, the supernova color-color diagnostic diagram was used to determine our observation date and correctly locate our magnitudes on the light curves. For this purpose, the use of supernova light curve templates, as well as reference supernovae, was also helpful. This work allowed us to verify the feasibility of performing precision astronomical photometry at the OAUNI.

The Dragonfly Spectral Line Mapper (DSLM) is a semi-autonomous, distributed-aperture based telescope design, featuring a modular setup of 120 Canon telephoto lenses, and equal numbers of ultra-narrowband filters, detectors, and other peripherals. Here we introduce the observatory software stack for this highly-distributed system. Its core is the Dragonfly Communication Protocol (DCP), a pure-Python hardware communication framework for standardized hardware interaction. On top of this are 120 REST-ful FastAPI web servers, hosted on Raspberry Pis attached to each unit, orchestrating command translation to the hardware and providing diagnostic feedback to a central control system running the global instrument control software. We discuss key features of this software suite, including docker containerization for environment management, class composition as a flexible framework for array commands, and a state machine algorithm which controls the telescope during autonomous observations.

The far-ultraviolet (FUV) flare activity of low-mass stars has become a focus in our understanding of the exoplanet atmospheres and how they evolve. However, direct detection of FUV flares and measurements of their energies and rates are limited by the need for space-based observations. The difficulty of obtaining such observations may push some works to use widely available optical data to calibrate multi-wavelength spectral models that describe UV and optical flare emission. These models either use single temperature blackbody curves to describe this emission, or combine a blackbody curve with archival spectra. These calibrated models would then be used to predict the FUV flare rates of low-mass stars of interest. To aid these works, we used TESS optical photometry and archival HST FUV spectroscopy to test the FUV predictions of literature flare models. We tested models for partially (M0-M2) and fully convective (M4-M5) stars, 40 Myr and field age stars, and optically quiet stars. We calculated FUV energy correction factors that can be used to bring the FUV predictions of tested models in line with observations. A flare model combining optical and NUV blackbody emission with FUV emission based on HST observations provided the best estimate of FUV flare activity, where others underestimated the emission at all ages, masses and activity levels, by up to a factor of 104 for combined FUV continuum and line emission and greater for individual emission lines. We also confirmed previous findings that showed optically quiet low-mass stars exhibit regular FUV flares.

The calibration of future wide field adaptive optics (WFAO) systems requires knowledge of the geometry of the system, in particular the alignment parameters between the sub-apertures of the wavefront sensors (WFS), pupil and deformable mirror (DM) actuator grid. Without this knowledge, closed-loop operation is not possible and the registration must be identified with an error significantly smaller than the sub-aperture size to achieve the nominal performance of the adaptive optics system. Furthermore, poor accuracy in this estimation will not only affect performance, but could also prevent the closed loop from being stable. Identification is not trivial because in a WFAO system several elements can move with respect to each other, more than in a SCAO system. For example, the pairing of the sub-aperture and the actuator grating on a DM conjugated to an altitude different from 0 can depend on the size of the pupil on the WFS, the exact conjugation of the DM, the position of the guide star and the field rotation. This is the same for each WFS/DM pair. SPRINT, System Parameters Recurrent INvasive Tracking, is a strategy for monitoring and compensating for DM/WFS mis-registrations and has been developed in the context of single conjugate adaptive optics (SCAO) systems for the ESO Extremely Large Telescope (ELT). In this work, we apply SPRINT in the context of WFAO systems with multiple WFSs and DMs, investigating the best approach for such systems, considering a simultaneous identification of all parameters or subsequent steps working on one DM at a time.

The population of strong lensing galaxies is a sub-set of intermediate-redshift massive galaxies, whose population-level properties are not yet well understood. In the near future, thousands of multiply imaged systems are expected to be discovered by wide-field surveys like Rubin Observatory's Legacy Survey of Space and Time (LSST) and Euclid. With the soon-to-be robust population of quadruply lensed quasars, or quads, in mind, we introduce a novel technique to elucidate the empirical distribution of the galaxy population properties. Our re-imagining of the prevailing strong lensing analysis does not fit mass models to individual lenses, but instead starts with parametric models of many galaxy populations, which include generally ignored mass distribution complexities and exclude external shear for now. We construct many mock galaxy populations with different properties and obtain populations of quads from each of them. The mock `observed' population of quads is then compared to those from the mocks using a model-free analysis based on a 3D sub-space of directly observable quad image properties. The distance between two quad populations in the space of image properties is measured by a metric $\eta$, and the distance between their parent galaxy populations in the space of galaxy properties is measured by $\zeta$. We find a well defined relation between $\eta$ and $\zeta$. The discovered relation between the space of image properties and the space of galaxy properties allows for the observed galaxy population properties to be estimated from the properties of their quads, which will be conducted in a future paper.

In this work we provide a data analysis of the BOSS galaxy clustering data with the recently proposed FreePower method, which adopts as parameters the power spectrum band-powers instead of a specific cosmological parametrization. It relies on the Alcock-Paczyński effect and redshift-space distortions, and makes use of one-loop perturbation theory for biased tracers. In this way, we obtain for the first time constraints on the linear growth rate, on the Hubble parameter, as well as on various bias functions, that are independent of a model for the power spectrum shape and thus of both the early and late-time cosmological modelling. We find at $z_{\rm eff}=0.38$, $f=0.60^{+0.11}_{-0.09}$, $H/H_0=1.052^{+0.060}_{-0.041}$, and at $z_{\rm eff}=0.61$, $f=0.69^{+0.14}_{-0.13}$, $H/H_0=1.36\pm0.10$. The low-$z$ $H/H_0$ result is at over 2$\sigma$ tension with Planck 2018 $\Lambda$CDM results. The precision in these results demonstrates the capabilities of the FreePower method. We also get constraints on the bias parameters which are in agreement with constraints from previous BOSS analyses, which serves as a cross-check of our pipeline.

Future low-noise cosmic microwave background (CMB) lensing measurements from e.g., CMB-S4 will be polarization dominated, rather than temperature dominated. In this new regime, statistically optimal lensing reconstructions outperform the standard quadratic estimator, but their sensitivity to extragalactic polarized foregrounds has not been quantified. Using realistic simulations of polarized radio and infrared point sources, we show for the first time that optimal Bayesian lensing from a CMB-S4-like experiment is insensitive to the expected level of polarized extragalactic foregrounds after masking, as long as an accurate foreground power spectrum is included in the analysis. For more futuristic experiments where these foregrounds could cause a detectable bias, we propose a new method to jointly fit for lensing and the Poisson foregrounds, generalizing the bias hardening from the standard quadratic estimator to Bayesian lensing.