Search

Several space missions are proposed or planned for the coming two decades dedicated or including mid- to high-resolution spectropolarimetry on a wide UV band. This includes the European instrument Pollux for the NASA HWO flagship mission, the NASA SMEX candidate Polstar, and the French nanosatellite demonstrator CASSTOR. We are developing UV polarimeters for these missions thanks to a R&D program funded by CNES. For the mid- and near-UV, i.e. above 120 nm, birefringent material (MgF2) can be used to produce a polarimeter. This is the baseline for Polstar, CASSTOR, and the MUV and NUV channels of Pollux. Prototypes have been built and tested with excellent results, and further tests are ongoing to fully characterize them. For the FUV channel of Pollux however, it is not possible to use this technology and we have instead studied a design based on mirrors only. We will present the various missions and instruments, their technical challenges, as well as the R&D work performed on UV polarimeters and the proposed design solutions., Comment: 4 pages, ICSO 2024

CARMENES stands for Calar Alto high-Resolution search for M dwarfs with Exoearths with Near-infrared and optical \'Echelle Spectrographs. CARMENES took six years from a concept to the start of operations, and a couple more years of initial data collection until the first science publication, but now is revolutionising our knowledge on exoplanets and their stars in our immediate vicinity. Here we describe what CARMENES is: (i) an ultra-stabilised two-channel spectrograph at an almost dedicated 3.5 m telescope in southern Spain that covers in high spectral resolution and without big gaps from 0.52 mum to 1.71 mum; (ii) a science project aimed at comprehensively searching for and studying planetary systems with nearby, bright, M-dwarf hosts, but that also investigates transiting planets around other stars; and (iii) the German-Spanish consortium that designed and built the instrument and that has operated it under guaranteed and legacy time observations., Comment: To be published in: Handbook of Exoplanets, 2nd edition, Hans Deeg and Juan Antonio Belmonte (eds. in chief), Springer International Publishing AG, part of Springer Nature. The corresponding images in full size can be downloaded at this link: https://carmenes.caha.es/ext/instrument/CARMENES_Handbook_Caballero_Figs.pdf (massive file)

We determined phase curves for 301,272 asteroids in the orange filter and 280,953 in the cyan filter from the latest ATLAS Solar System Catalog V2 (SSCAT-2). Among them, 3,345 and 492 asteroids in the orange and cyan filters, respectively, have uncertainties below 15%. Our simple model, which considers only the apparition effect, showed good consistency with more sophisticated methods requiring much less computational time. Database cross-matching allowed us to analyze G1 and G2 distributions according to taxonomy. We conducted two-dimensional Kolmogorov-Smirnov tests to investigate two distinct aspects: similarities in paired G1, G2 distributions across different taxa and wavelength dependency within the same taxa. When comparing different taxa, we couldn't reject the null hypothesis for 11% of the orange sample and 31% of the cyan sample, indicating more disparities in the orange filter. For wavelength dependency, paired distributions of G1, G2 (o) vs. G1, G2 (c) showed statistically significant differences across all complexes, except for the A class. Our analysis suggests that while phase coloring behaviors are observed without a clear preference for reddening or bluening at phase angles below 5 degrees, reddening predominates in the 10 - 30 degrees range. We also observed smaller uncertainties in G2 than in G1. Simulations showed that G2 is less sensitive to lack of data at small phase angles. This is related to the definition of the H,G1,G2 function, where G1 contributes more to the opposition effect and G2 the linear part of the phase curve. Our catalog-independent algorithms are adaptable to new data sets, including future LSST data., Comment: Accepted for publication (ICARUS)

Comets are the most pristine planetesimals left from the formation of the Solar System. They carry unique information on the materials and the physical processes which led to the presence of planets and moons. Many important questions about cometary physics, such as origin, constituents and mechanism of cometary activity, remain unanswered. The next perihelion of comet 1P/Halley, in 2061, is an excellent opportunity to revisit this object of outstanding scientific and cultural relevance. In 1986, during its latest approach to the Sun, several flyby targeted Halley's comet to observe its nucleus and shed light on its properties, origin, and evolution. However, due to its retrograde orbit and high ecliptic inclination, the quality of data was limited by the large relative velocity and short time spent by the spacecraft inside the coma of the comet. A rendezvous mission like ESA/Rosetta would overcome such limitations, but the trajectory design is extremely challenging due to the shortcomings of current propulsion technology. Given the considerable lead times of spacecraft development and the long duration of the interplanetary transfer required to reach the comet, it is imperative to start mission planning several decades in advance. This study presents a low-thrust rendezvous strategy to reach the comet before the phase of intense activity during the close approach to the Sun. The trajectory design combines a gravity-assist maneuver with electric propulsion arcs to maximize scientific payload mass while constraining transfer duration. A propulsive plane change maneuver would be prohibitive. To keep the propellant budget within reasonable limits, most of the plane change maneuver is achieved via either a Jupiter or a Saturn flyby. The interplanetary low-thrust gravity-assisted trajectory design strategy is described, followed by the presentation of multiple proof-of-concept solutions.

The spatial extension of active regions (ARs) of the Sun can vary from one case to the next. This is a problem when studying solar flares with Convolutional Neural Networks (CNNs) as they generally use input images of a fixed size. Different processes can be performed to retrieve a database with homogeneous-sized data, such as resizing. Unfortunately, key features can be lost or distorted during these processes. This can lead to a deterioration of the ability of CNNs to classify flares of different soft X-ray classes, especially those from ARs with complex structures. Our work aims to implement and test a CNN architecture that retains the full features of the original resolution of the input images. We compare the performance of two CNN architectures for solar flare prediction: the first is a traditional CNN with resized input whereas the other implements a spatial pyramid pooling (SPP) layer without any input resizing. Both are trained on the Spaceweather HMI Active Region Patch line-of-sight magnetogram database. We also study two cases of binary classification. In the first case, our model distinguishes ARs producing flares in less than 24h of class greater or equal to C1.0 from ARs producing flares in more than 24h or never; in the second case, it distinguishes ARs producing flares in less than 24h of class greater or equal to M1.0 from the other ARs. Our models implementing an SPP layer outperform the traditional CNN models when predicting flares greater or equal to C1.0 within 24h. However, their performances degrade sharply along the other models studied in this paper, when trained to classify images greater or equal to M1.0 flares. The degradation in SPP models when classifying only images greater or equal to M1.0 flares as positive may be attributed to its success in identifying features that appear in ARs a few hours before the flare, independently of their soft X-ray class., Comment: This paper will be published in the Astrophysics and Astronomy journal, Volume 695, article A65. See https://www.aanda.org/component/article?access=doi&doi=10.1051/0004-6361/202449671

Context. We present the COMetary dust TAIL Simulator (COMTAILS), a numerical Monte Carlo code to generate images of dust tail brightness from comets and active asteroids in the Solar System. Aims. We describe a numerical code, available to interested users, capable of generating simulated images of dust tail brightness for comparison with observations, to retrieve key dust parameters including size distribution, ejection velocities, and dust loss rates. An optional stellar field can be included in the background allowing for the assessment of stellar extinction by the tail, which can be compared with observational data. Methods. A Monte Carlo procedure was used to obtain the simulated images. The orbital parameters of the targets and their helio-centric positions and velocities were obtained from the JPL Horizons on-line ephemeris system. Results. Earlier versions of this code have been used to characterize the dust environment of various targets. In this study we present recent examples that demonstrate its ability to fit observed images of the long-period comets C/2024 N3 (NEOWISE) and C/2023 A3 (Tsuchinshan-ATLAS). The code is readily applicable to future targets that may be identified for the upcoming European Space Agency Comet Interceptor mission., Comment: Accepted by A&A

Polycyclic aromatic hydrocarbons (PAHs) are key contributors to interstellar infrared emission bands (AIBs). However, current spectral databases for AIB analysis face limitations because of the neglect of anharmonicity in vibrations and of temperature effects, which result from the high computational cost of traditional quantum chemical calculations (QCCs). In this work, we present a machine-learning-based molecular dynamics (MLMD) approach that efficiently computes anharmonic infrared (IR) spectra while accounting for temperature effects. MLMD achieves predictive accuracy comparable to that of QCCs but with significantly reduced computational cost, scaling linearly with the number of atoms in the system. We apply MLMD to calculate the anharmonic spectra of 1\,704 PAHs in the NASA Ames PAH IR Spectroscopic Database with up to 216 carbon atoms, demonstrating its ability to handle large molecular systems. Our results suggest that MLMD can facilitate the creation of extensive molecular spectral databases, advancing AI-assisted analyses of astronomical IR spectra, particularly with upcoming data from the James Webb Space Telescope.

This document provides a user guide for reducing UVIT data using CCDLAB. While CCDLAB offers a straightforward data reduction work-flow, users may encounter certain challenges that require additional guidance. This guide provides instructions by addressing common issues related to key processing steps, including WCS solutions and VIS drift tracking., Comment: 4 pages, 8 Figures

The magnetic field morphology of relativistic jets can be studied with circular polarization (CP). Recent 3D relativistic magnetohydrodynamic (RMHD) simulations coupled with radiative transfer calculations make strong predictions about the level (and morphology) of the jet's CP emission. These simulations show that the sign of CP and the electric vector position angle (EVPA) are both sensitive to the jet's magnetic field morphology within the radio core. We probe this theory by exploring if the jet's radio core EVPA orientation is consistent with the observed sign of the core CP in deep full-track polarimetric observations. We aim to probe the nature of linear polarization and CP in the innermost regions of jets from a small sample of nine blazars. This sample includes sources that have exhibited: (i) positive CP, (ii) negative CP, or (iii) positive & negative CP simultaneously in the radio core region. Nine blazar sources were observed using the VLBA at both 15 GHz and 23 GHz. Our self-calibration relies on a physically based model applied in DoG-HiT resulting in more accurate gains. We consider compact Stokes V structures instead of assuming it to be zero, which is crucial given the significant non-zero CP fraction observed at long baselines. We observe robust, relatively high degrees of fractional circular polarization m_c=(0.32 +- 0.2)% at 15 GHz and m_c=(0.59 +-0.56)% at 23 GHz. We observe consistent polarized structure and EVPA orientation over time when comparing our analysis with archival MOJAVE data. Theoretical predictions indicate a clear favored toroidal magnetic field orientation within the extended jet emission of the reconstructed signal of the blazar 0149+218. The jet structures of 1127-145 and 0528+134, even in superresolution, exhibit characteristics aligned with helical or poloidal magnetic nature.

This work introduces an advanced technique optimized for detecting photons generated by charged particles, leveraging Skipper-CCD sensors. By analyzing background sources and detection efficiencies, the technique achieves strong agreement between experimental results and Cherenkov-based simulations. It also provides a robust framework for investigating secondary photon production in environments with high fluxes of ionizing particles, such as those anticipated in space-based astronomical instruments. These secondary photons present a critical challenge as background noise for next-generation single-photon resolving imagers used to study faint celestial objects. Furthermore, the method exhibits significant potential for broader applications, including exploring photon generation in various substrate materials and examining their transport through multiple interfaces.

Light pollution, a rapidly escalating anthropogenic phenomenon driven by the excessive and often inefficient use of artificial lighting, has profound implications for astronomy, ecology, and human health. This study presents the first comprehensive characterization of night sky quality in Colombia, focusing on sites of astronomical and ecological significance. The selected locations include the Astronomical Observatory of UTP, the Tatacoa Desert, the Bogot\'a Botanical Garden, and Cerro Guadalupe. Utilizing the Sky Quality Camera, we collected all-sky data to measure surface brightness and correlated color temperature of the night sky. Our findings reveal a significant loss of natural sky visibility in urban areas and demonstrate the detrimental effects of artificial lighting on critical astronomical sites such as La Tatacoa. This study provides a crucial foundation for future research and informs the development of public policies aimed at preserving the night sky., Comment: Accepted for publication in RMxAA, 41 pages, 24 figures, 6 tables

The level and uncertainty of the particle induced background in CCD detectors plays a crucial role for future X-ray instruments, such as the Wide Field Imager (WFI) onboard Athena. To mitigate the background systematic uncertainties, which will limit the Athena science goals, we aim to understand the relationship between the energetic charged particles interacting in the detector and satellite, and the instrumental science background to an unprecedented level. These particles produce easily identified "cosmic-ray tracks" along with less easily identified signals produced by secondary particles, e.g., X-rays generated by particle interactions with the instrument and indistinguishable from genuine sky X-rays. We utilize the Small Window Mode of the PN camera onboard XMM-Newton to understand the time, spatial and energy dependence of the various background components, particularly the particle induced background. While the distribution of particle events follows expected detector readout patterns, we find a particle track length distribution inconsistent with the simple, isotropic model. We also find that the detector mode-specific readout results in a shifted Cu fluorescent line. We illustrate that on long timescales the variability of the particle background correlates well with the solar cycle. This 20-year lightcurve, can be reproduced by a particle detector onboard Chandra, the HRC anti-coincidence shield. We conclude that the self-anti-coincidence method of removing X-ray-like events near detected particle tracks in the same frame can be optimized with the inclusion of additional information, such as the energy of the X-ray. The results presented here are relevant for any future pixelated X-ray imaging detector, and could allow the WFI to probe to truly faint X-ray surface brightness., Comment: 35 pages, 16 figures. Accepted for publication in JATIS

Blazars are interesting source candidates for astrophysical neutrino emission. Multi-messenger lepto-hadronic models based on proton-photon (p-gamma) interactions result in predictions for the neutrino spectra which are typically strongly peaked at high energies. In contrast, statistical analyses looking to associate blazars and neutrinos often assume a power-law spectral shape, putting the emphasis at lower energies. We aim to examine the impact of such spectral modelling assumptions on the associations of neutrinos with blazars. We use hierarchical_nu, a Bayesian framework for point source searches, and incorporate the theoretical predictions for neutrino spectra through a dedicated spectral model and priors on the relevant parameters. Our spectral model is based on the predictions presented in Rodrigues et al. (2024) for a selection of intermediate and high synchrotron peaked blazars that have been connected with IceCube neutrino alert events. We apply our model to the 10 years of publicly available IceCube data, focusing on the Northern hemisphere. Out of 29 source candidates, we find five sources, including TXS 0506+056, that have an association probability greater than 0.5 to at least one event. The p-gamma spectra typically lead to a lower overall number of associated events compared to the power-law case, but retain or even enhance strong associations to high-energy events. Our results demonstrate that including more information from theoretical predictions can allow for more interpretable source-neutrino connections., Comment: Submitted to ApJ

Identifying weak gravitational wave signals in noise and estimating the source properties require high-precision waveform templates. Numerical relativity (NR) simulations can provide the most accurate waveforms. However, it is challenging to compute waveform templates in high-dimensional parameter space using NR simulations due to high computational costs. In this work, we implement a novel waveform mapping method, which is an alternative approach to the existing analytical approximations, based on closed-form continuous-time neural networks. This machine-learning-based method greatly improves the efficiency of calculating waveform templates for arbitrary source parameters, such as the binary mass ratio and the spins of component black holes. Based on this method, we present \textit{BHP2NRMLSur}, a class of models (including nonspinning and spin-aligned ones) that maps point-particle black hole perturbation theory waveforms into NR and surrogate waveforms. The nonspinning model provides highly accurate waveforms that match the NR waveforms to the level of $\gtrsim 0.995$. The spin-aligned model reduces the required input parameters and hence improves the efficiency of the waveform generation -- it takes a factor of $\sim 50$ less time than existing NR surrogate models to generate $100,000$ waveforms, with a mismatch of $<0.01$ compared to the NR waveforms from the Simulating eXtreme Spacetimes collaboration.

The study of the atmosphere of exoplanets orbiting white dwarfs is a largely unexplored field. With WD\,0806-661\,b, we present the first deep dive into the atmospheric physics and chemistry of a cold exoplanet around a white dwarf. We observed WD 0806-661 b using JWST's Mid-InfraRed Instrument Low-Resolution Spectrometer (MIRI-LRS), covering the wavelength range from 5 -- 12~$\mu \rm{m}$, and the Imager, providing us with 12.8, 15, 18 and 21\,$\mu$m photometric measurements. We carried the data reduction of those datasets, tackling second-order effects to ensure a reliable retrieval analysis. Using the \textsc{TauREx} retrieval code, we inferred the pressure-temperature structure, atmospheric chemistry, mass, and radius of the planet. The spectrum of WD 0806-661 b is shaped by molecular absorption of water, ammonia, and methane, consistent with a cold Jupiter atmosphere, allowing us to retrieve their abundances. From the mixing ratio of water, ammonia and methane we derive $\rm{C/O} = 0.34 \pm 0.06$, $\rm{C/N} = 14.4 ^{+2.5}_{-1.8}$ and $\rm{N/O} = 0.023 \pm 0.004$ and the ratio of detected metals as proxy for metallicity. We also derive upper limits for the abundance of CO and $\rm{CO_2}$ ($1.2\cdot10^{-6} \rm{\,and\,} 1.6\cdot10^{-7}$ respectively), which were not detected by our retrieval models. While our interpretation of WD\,0806-661\,b's atmosphere is mostly consistent with our theoretical understanding, some results -- such as the lack of evidence for water clouds, an apparent increase in the mixing ratio of ammonia at low pressure, or the retrieved mass at odds with the supposed age -- remain surprising and require follow-up observational and theoretical studies to be confirmed., Comment: Accepted for publication in ApJL in March 2025. 9 pages and 4 figures

This paper introduces Firmamento, a new online platform designed for astronomical research, particularly for studying blazars and other multi-messenger emitters. Firmamento provides access to a wealth of astronomical data, including imaging, spectral, and timing information, along with tools for data analysis and machine learning. A key feature is the Error Region Counterpart Identifier (ERCI), which helps locate potential counterparts to gamma-ray and other high-energy sources pinpointing objects exhibiting blazar-like characteristics. Firmamento grants access to spectral data across the entire electromagnetic spectrum, offering well-populated SED generation and VHE detectability estimation. The platform also supports user-provided data, allowing researchers to incorporate their own source lists., Comment: 8 pages, 7 figure, proceedings of the 8th Heidelberg gamma-raysymposium, Milan. Submitted to MemSAIt

Numerical simulations of multidimensional astrophysical fluids present considerable challenges. However, the development of exascale computing has significantly enhanced computational capabilities, motivating the development of new codes that can take full advantage of these resources. In this article, we introduce HERACLES++, a new hydrodynamics code with high portability, optimized for exascale machines with different architectures and running efficiently both on CPUs and GPUs. The code is Eulerian and employs a Godunov finite-volume method to solve the hydrodynamics equations, which ensures accuracy in capturing shocks and discontinuities. It includes different Riemann solvers, equations of state, and gravity solvers. It works in Cartesian and spherical coordinates, either in 1-D, 2-D, or 3-D, and uses passive scalars to handle gases with several species. The code accepts a user-supplied heating or cooling term to treat a variety of astrophysical contexts. In addition to the usual series of benchmarking tests, we use HERACLES++ to simulate the propagation of a supernova shock in a red-supergiant star envelope, from minutes after core collapse until shock emergence. In 1-D, the results from HERACLES++ are in agreement with those of V1D for the same configuration. In 3-D, the Rayleigh-Taylor instability develops and modifies the 1-D picture by introducing density and composition fluctuations as well as turbulence. The focus on a wedge, rather than the full solid angle, and the ability of running HERACLES++ with a large number of GPUs allow for long-term simulations of 3-D supernova ejecta with a sub-degree resolution. Future development is to extend HERACLES++ to a radiation-hydrodynamics code.

Axions are well-motivated dark matter particles. Many experiments are looking for their experimental evidence. For haloscopes, the problem reduces to the identification of a peak above a noisy baseline. Its modeling, however, may problematic. State-of-the-art analysis rely on the Savitzky-Golay (SG) filtering, which is intrinsically affected by any possible over fluctuation, leading to biased results. In this paper we study the efficiency that different extensions of SG can provide in the peak reconstruction in a standard haloscope-like experiment. We show that, once the correlations among bins are taken into account, there is no appreciable difference. The standard SG remains the advisable choice because of its numerical efficiency., Comment: 17 pages, 5 figures

The next generation of gravitational-wave observatories will achieve unprecedented strain sensitivities with an expanded observing band. They will detect ${\cal O}(10^5)$ binary neutron star (BNS) mergers every year, the loudest of which will be in the band for $\approx 90$ minutes with signal-to-noise ratios $\approx 1500$. Current techniques will not be able to determine the astrophysical parameters of the loudest of next-gen BNS signals. We show that subtleties arising from the rotation of the Earth and the free-spectral range of gravitational-wave interferometers dramatically increases the complexity of next-gen BNS signals compared to the one-minute signals seen by LIGO--Virgo. Various compression methods currently relied upon to speed up the most expensive BNS calculations -- reduced-order quadrature, multi-banding, and relative binning -- will no longer be effective. We carry out reduced-order inference on a simulated next-gen BNS signal taking into account the Earth's rotation and the observatories' free-spectral range. We show that standard data compression techniques become impractical, and the full problem becomes computationally infeasible, when we include data below $\approx 16$Hz -- a part of the observing band that is critical for precise sky localisation. We discuss potential paths towards solving this complex problem.

Modern Imaging Atmospheric Cherenkov Telescopes (IACTs) generate a huge amount of data that must be classified automatically, ideally in real time. Currently, machine learning-based solutions are increasingly being used to solve classification problems. However, these classifiers require proper training data sets to work correctly. The problem with training neural networks on real IACT data is that these data need to be pre-labeled, whereas such labeling is difficult and its results are estimates. In addition, the distribution of incoming events is highly imbalanced. Firstly, there is an imbalance in the types of events, since the number of detected gamma quanta is significantly less than the number of protons. Secondly, the energy distribution of particles of the same type is also imbalanced, since high-energy particles are extremely rare. This imbalance results in poorly trained classifiers that, once trained, do not handle rare events correctly. Using only conventional Monte Carlo event simulation methods to solve this problem is possible, but extremely resource-intensive and time-consuming. To address this issue, we propose to perform data augmentation with artificially generated events of the desired type and energy using conditional generative adversarial networks (cGANs), distinguishing classes by energy values. In the paper, we describe a simple algorithm for generating balanced data sets using cGANs. Thus, the proposed neural network model produces both imbalanced data sets for physical analysis as well as balanced data sets suitable for training other neural networks., Comment: 19 pages, 10 figures, Proceedings of The 8th International Conference on Deep Learning in Computational Physics, June 19-21, 2024, Moscow, Russia

We present a novel deep generative model, named GenMDI, to improve the temporal resolution of line-of-sight (LOS) magnetograms of solar active regions (ARs) collected by the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO). Unlike previous studies that focus primarily on spatial super-resolution of MDI magnetograms, our approach can perform temporal super-resolution, which generates and inserts synthetic data between observed MDI magnetograms, thus providing finer temporal structure and enhanced details in the LOS data. The GenMDI model employs a conditional diffusion process, which synthesizes images by considering both preceding and subsequent magnetograms, ensuring that the generated images are not only of high-quality, but also temporally coherent with the surrounding data. Experimental results show that the GenMDI model performs better than the traditional linear interpolation method, especially in ARs with dynamic evolution in magnetic fields., Comment: 11 pages, 7 figures

The Exoplanet Transmission Spectroscopy Imager (ETSI) amalgamates a low resolution slitless prism spectrometer with custom multi-band filters to simultaneously image 15 spectral bandpasses between 430 nm and 975 nm with an average spectral resolution of $R = \lambda/\delta\lambda \sim 20$. ETSI requires only moderate telescope apertures ($\sim2$ m) and is capable of characterizing an exoplanet atmosphere in as little as a single transit, enabling selection of the most interesting targets for further characterization with other ground and space-based observatories and is also well suited to multi-band observations of other variable and transient objects. This enables a new technique, common-path multi-band imaging (CMI), used to observe transmission spectra of exoplanets transiting bright (V$<$14 magnitude) stars. ETSI is capable of near photon-limited observations, with a systematic noise floor on par with the Hubble Space Telescope and below the Earth's atmospheric amplitude scintillation noise limit. We report the as-built instrument optical and optomechanical design, detectors, control system, telescope hardware and software interfaces, and data reduction pipeline. A summary of ETSI's science capabilities and initial results are also included., Comment: 33 pages, 17 figures, published in JATIS

One of the most prolific methods of studying exoplanet atmospheres is transmission spectroscopy, which measures the difference between the depth of an exoplanet's transit signal at various wavelengths and attempts to correlate the depth changes to potential features in the exoplanet's atmosphere. Here we present reconnaissance observations of 21 exoplanet atmospheres measured with the Exoplanet Transmission Spectroscopy Imager (ETSI), a recently deployed spectro-photometer on the McDonald Observatory Otto Struve 2.1 m telescope. ETSI measurements are mostly free of systematics through the use of a novel observing technique called common-path multi-band imaging (CMI), which has been shown to achieve photometric color precision on-par with space-based observations (300ppm or 0.03%). This work also describes the various statistical tests performed on the data to evaluate the efficacy of the CMI method and the ETSI instrument in combination. We find that none of the 8 comparisons of exoplanet atmospheres measured with ETSI and other observatories (including the Hubble Space Telescope) provide evidence that the spectra are statistically dissimilar. These results suggest that ETSI can provide initial transmission spectroscopy observations for a fraction of the observational and monetary overhead previously required to detect an exoplanet's atmosphere. Ultimately these reconnaissance observations increase the number of planets with transmission spectroscopy measurements by ~10% and provide an immediate prioritization of 21 exoplanets for future follow-up with more precious observatories, such as the James Webb Space Telescope. The reconnaissance spectra are available through the Filtergraph visualization portal at the URL: https://filtergraph.com/etsi/., Comment: 28 pages, 8 figures, 11 tables, article published in AJ. ETSI data is available on request and spectra can be found at the URL https://filtergraph.com/etsi/

We discuss the possible synergies for cosmology between SKAO and ESO facilities, focusing on the combinations SKA-Mid with the Multi-Object Spectrograph Telescope (4MOST) instrument built for ESO's Visible and Infrared Survey Telescope for Astronomy (VISTA), and SKA-Low with ESO's Extremely Large Telescope (ELT) multi-object spectrograph MOSAIC. Combining multiple tracers allows for tackling systematics and lifting parameter degeneracies. It will play a crucial role in the pursuit of precision cosmology., Comment: 4 pages, 3 figures. Based on the discussion outcomes of the splinter session "EoR and cosmology" of the 2023 "Coordinated Surveys of the Southern Sky" workshop, Garching b. Munchen, February 27 - March 3, 2023

Binary population synthesis simulations allow detailed modelling of gravitational-wave sources from a variety of formation channels. These population models can be compared to the observed catalogue of merging binaries to infer the uncertain astrophysical input parameters describing binary formation and evolution, as well as relative rates between various formation pathways. However, it is computationally infeasible to run population synthesis simulations for all variations of uncertain input physics. We demonstrate the use of normalising flows to emulate population synthesis results and interpolate between astrophysical input parameters. Using current gravitational-wave observations of binary black holes, we use our trained normalising flows to infer branching ratios between multiple formation channels, and simultaneously infer common-envelope efficiency and natal spins across a continuous parameter range. Given our set of formation channel models, we infer the natal spin to be $0.04^{+0.04}_{-0.01}$, and the common-envelope efficiency to be $>3.7$ at 90% credibility, with the majority of underlying mergers coming from the common-envelope channel. Our framework allows us to measure population synthesis inputs where we do not have simulations, and better constrain the astrophysics underlying current gravitational-wave populations., Comment: 20 pages, 6 figures, 2 tables, 1 appendix. Submitted to Astrophysical Journal. Data release at https://doi.org/10.5281/zenodo.14967688

Searches for continuous gravitational waves from isolated compact objects and those in binary systems aim to detect non-axisymmetric, deformed neutron stars at particular locations in the Galaxy or all-sky. However, a large fraction of known pulsars have rotational frequencies that lie outside the audio frequency band, rendering current detectors insensitive to these pulsars. In this work, we show that DECIGO, a future space-based deci-hertz gravitational-wave interferometer, will be sensitive to severely deformed compact objects, e.g. hybrid stars, neutron stars, or magnetars. We estimate the number of possible compact objects that could be detected with such high deformations, both via their individual continuous gravitational-wave emission and the stochastic gravitational-wave background created by a superposition of gravitational waves from the $\sim 10^8$ compact objects in the Galaxy. Furthermore, we show that the existence of such compact objects could be probed across a wide parameter space at a fraction of the computational cost of current searches for isolated compact objects and those in binary systems. For known pulsars, we will be able to both beat the spin-down limit and probe the Brans-Dicke modified theory of gravity parameter $\zeta<1$ for approximately 85% of known pulsars with $f_{\rm gw}<10$ Hz, the latter of which is currently only possible for $O(10)$ pulsars. DECIGO will thus open a new window to probe highly deformed compact objects and over half of the known pulsars, both of which are currently inaccessible to ground-based detectors., Comment: 12 pages, 9 figures, comments are welcome!

In this paper, we present the linear decomposition method (LDM), which we developed to detect and analyze pulsar profile variations and mode changing behaviour. We developed LDM utilizing the likelihood function approach assuming the Gaussian noise. The LDM projects pulse profiles onto significance-ordered orthonormal vector bases. We show that the method is similar to the principal component analysis (PCA), but LDM can handle more general situations. We use simulated dataset and data from the Kunming 40-m radio telescope to demonstrate the application of the LDM. We found that the LDM successfully identified mode changes for well-known mode-changing PSR B0329+54 and found a continuous pulse profile evolution for PSR B0355+54 . We also show that the LDM can be used to improve the timing precision for mode changing PSR B0329+54., Comment: 10 pages, 10 figures, accepted for publication in ApJ

A recent paper demonstrated the existence of a secular aberration drift term in stellar proper motions that arises when transforming an astrometric catalogue defined for an observer at rest with respect to the solar system barycentre to some other reference frame in which, for example, the observer is at rest with respect to the Galactic centre. Such a transformation requires an accurate and precise estimate of the velocity of the solar system barycentre. It was argued that the Gaia catalogue construction should account for this effect and also for the aberrational effect due to acceleration of the solar system barycentre. We argue that these two effects should not be accounted for in the construction of the Gaia astrometric catalogue. We briefly review the Gaia catalogue reference frame, the concepts of stellar aberration and secular aberration drift, and their observable consequences. The Gaia catalogue is (and should be) constructed in the Barycentric Celestial Reference System: the reference system with the origin at the solar system barycentre as defined by the underlying solar system ephemerides. We explain that the Gaia catalogue is consistent with the International Celestial Reference System despite the presence of proper motion terms due to the acceleration of the solar system barycentre. We also explain why transformation of the astrometry to a frame in which the observer is at rest with respect to the Galactic centre or distant universe is not needed for the interpretation of stellar kinematics, and that there are practical concerns with such a transformation. The estimation of the velocity and acceleration of the solar system barycentre, although important as a matter of scientific investigation, are not needed for the construction of the Gaia astrometric catalogue., Comment: 6 pages, 1 figure, accepted for publication in Astronomy & Astrophysics

The Compton camera is a sensitive imaging detector for soft gamma-rays. Compton Reconstruction can not only give imaging capability but also remove background events to achieve good sensitivity. However, the angular resolution is in principle limited to several degrees. In this paper, we propose a novel concept of Compton camera incorporating shadow effects. We consider multi-column Compton camera (MCCC), consisting of stacked Si pixel sensors. Each of columns is separated each other to create shadow effects. This design achieves an angular resolution of less than 1 degree within around 1 degree from the center of the field-of-view by just modifying a conventional Si-stacked Compton camera and keeping advantages (wide field-of-view and good sensitivity) of conventional Compton camera. Here we validated the concept of proposed Compton camera through Monte-Carlo simulation. MCCC with 1 m column height, 0.5 mm pixel size, 100 layers, and 10 columns for the 1-D direction can distinguish two sources separated by 0.1 degree with 0.6M Compton-reconstructed events., Comment: published in NIM-A

With the advent of GRAVITY+, the upgrade to the beam combiner GRAVITY at the Very Large Telescope Interferometer (VLTI), fainter and higher redshift active galactic nuclei (AGNs) are becoming observable, opening an unprecedented opportunity to further our understanding of the cosmic coevolution of supermassive black holes and their host galaxies. To identify an initial sample of high-redshift type~1 AGNs that can be observed with GRAVITY+, we have obtained spectroscopic data with NTT/SOFI of the most promising candidates. Our goal is to measure their broad line region (BLR) fluxes and assess their physical geometries by analysing the spectral profiles of their Balmer lines. We present 29 $z$ $\sim$ 2 targets with strong H$\alpha$ emission in the $K$-band. Their line profiles are strongly non-Gaussian, with a narrow core and broad wings. This can be explained as a combination of rotation and turbulence contributing to the total profile or two physically distinct inner and outer regions. We find small H$\alpha$ virial factors, which we attribute to the low full-width-half-maximum (FWHM)/$\sigma$ ratios of their non-Gaussian profiles, noting that this can lead to discrepancies in black hole masses derived from scaling relations. We also find two targets that show tentative evidence of BLRs dominated by radial motions. Lastly, we estimate the expected differential phase signals that will be seen with GRAVITY+, which will provide guidance for the observing strategy that will be adopted., Comment: Accepted for publication at A&A

We present a new approach to report in the Section 4 of BIPM Circular T daily values of the offset between UTC and the predictions of UTC broadcast by the GNSS, this quantity we name bUTC_GNSS. In this approach, the determination of UTC - bUTC_GNSS is based on data collected by several multi-GNSS stations in selected time laboratories worldwide. Test computations over a 7-month period from July 2022 to January 2023 show that the offset between UTC and bUTC_GNSS was between 30 and 50 ns for GLONASS, between 5 and 20 ns for BeiDou, and between -5 and +5 ns for GPS and Galileo. We derive the uncertainty on the reported values, which is 4.1 ns for BeiDou and GPS, 3.7 ns for Galileo and 6.6 ns for GLONASS and show that, over the test period, the reported values of UTC - bUTC_GNSS and the solutions obtained from each multi-GNSS station are all consistent within the 1-sigma uncertainties.

We describe the basic instrument detrending software for the NIRWALS spectrograph on the SALT telescope. Its basic purpose is to process multiple non-destructive reads of increasing exposure time into a final image reflecting the observed source intensity as expressed in counts (or electrons) per second, including its uncertainty. The output products of this pipeline can then be used to as input for follow-up data processing to apply wavelength solutions and extract fluxes for individual fibers. All pipeline code is implemented in python and can be obtained via common portals such as PyPI and GitHub. Additional information on how to use the code is available online at nirwals.readthedocs.io, Comment: online documentation is available at https://nirwals.readthedocs.io/

Context. Radiation plays a significant role in solar and astrophysical environments as it may constitute a sizeable fraction of the energy density, momentum flux, and the total pressure. Modelling the dynamic interaction between radiation and magnetized plasmas in such environments is an intricate and computationally costly task. Aims. The goal of this work is to demonstrate the capabilities of the open-source parallel, block-adaptive computational framework MPI-AMRVAC, in solving equations of radiation-magnetohydrodynamics (RMHD), and to present benchmark test cases relevant for radiation-dominated magnetized plasmas. Methods. The existing magnetohydrodynamics (MHD) and flux-limited diffusion (FLD) radiative-hydrodynamics physics modules are combined to solve the equations of radiation-magnetohydrodynamics (RMHD) on block-adaptive finite volume Cartesian meshes in any dimensionality. Results. We introduce and validate several benchmark test cases such as steady radiative MHD shocks, radiation-damped linear MHD waves, radiation-modified Riemann problems and a multi-dimensional radiative magnetoconvection case. We recall the basic governing Rankine-Hugoniot relations for shocks and the dispersion relation for linear MHD waves in the presence of optically thick radiation fields where the diffusion limit is reached. The RMHD system allows for 8 linear wave types, where the classical 7-wave MHD picture (entropy and three wave pairs for slow, Alfven and fast) is augmented with a radiative diffusion mode. Conclusions. The MPI-AMRVAC code now has the capability to perform multidimensional RMHD simulations with mesh adaptation making it well-suited for larger scientific applications to study magnetized matter-radiation interactions in solar and stellar interiors and atmospheres., Comment: Accepted to Astronomy and Astrophysics journal

In order to reduce light pollution we have to reduce the overall amount of artificial light emissions. This is a consequence of the basic rules governing the propagation of light in the terrestrial environment. In this work I revisit the physical laws causing that (i) "all" artificial light emitted outdoors is lost or pollutant, (ii) the negative effects of light pollution depend monotonically on the local concentration of photons of anthropogenic origin, and (iii) this concentration depends linearly on the total artificial light emissions, weighted by the light pollution propagation functions. Setting total emission limits becomes necessary in order to ensure that the negative effects of light pollution do not surpass red-lines of unacceptable degradation of the natural night. Once these red-lines are socially agreed, monitoring compliance becomes a relevant task. Several complementary methods for assessing total light emissions are being used nowadays, including public inventories of installed lights, direct radiance measurements from ground or low Earth orbit satellites, and scattered radiance measurements using ground based detectors (night sky brightness monitoring). While updated administrative inventories could in principle be trusted, independent verification is a must. The required measurements pose in turn significant challenges: we discuss here how the variability of the terrestrial atmosphere sets a lower limit on the minimum emission changes that can be reliably detected by measurements, and propose some ways to improve this performance., Comment: 18 pages, 4 figures. This preprint is an extended version of the talk with the same title presented at the 3rd International Conference on Environmental and Astronomical Light Pollution, EALPO 2024, Niepo{\l}omice, Poland. September 6th-8th 2024

Distinguishing galaxies as either fast or slow rotators plays a vital role in understanding the processes behind galaxy formation and evolution. Standard techniques, which are based on the $\lambda_R$-spin parameter obtained from stellar kinematics, frequently face difficulties to classify fast and slow rotators accurately. These challenges arise particularly in cases where galaxies have complex interaction histories or exhibit significant morphological diversity. In this paper, we evaluate the performance of a Convolutional Neural Network (CNN) on classifying galaxy rotation kinematics based on stellar kinematic maps from the SAMI survey. Our results show that the optimal CNN architecture achieves an accuracy and precision of approximately 91% and 95% on the test dataset, respectively. Subsequently, we apply our trained model to classify previously unknown rotator galaxies for which traditional statistical tools have been unable to determine whether they exhibit fast or slow rotation, such as certain irregular galaxies or those in dense clusters. We also used Integrated Gradients (IG) to reveal the crucial kinematic features that influenced the CNN's classifications. This research highlights the power of CNNs to improve our comprehension of galaxy dynamics and emphasizes their potential to contribute to upcoming large-scale Integral Field Spectrograph (IFS) surveys., Comment: 18 pages, 10 figures. Accepted to the Special Issue "Universe: Feature Papers 2024: Galaxies and Clusters" of MDPI Galaxies

For space-borne gravitational wave detectors,such as LISA and TianQin ,the disturbance caused by the coupling of test masses and the external magnetic fields is one of the main sources of the residual acceleration noise. Although the detection frequency band is from 0.1 mHz to 1 Hz, magnetic fields with frequencies higher than 1 Hz can still contribute to the noise through down conversion effect. Therefore, it is necessary to measure the AC magnetic susceptibility or magnetic response of the test mass at higher frequency for the evaluation of the magnetic noise. In this work, we propose a magnetic field response measurement method by directly probing the induced magnetic field of the test mass placed in a spatially uniform magnetic field. The frequency can be measured up to 1500 Hz, satisfying the requirement of space-borne gravitational wave detection., Comment: 18 pages,6 figures

The Solar Ultraviolet Imaging Telescope (SUIT) is one of the seven payloads on board Aditya-L1 mission of the Indian Space Research Organization (ISRO). SUIT provides full and partial disk images of the Sun in the 200-400 nm wavelength range. This would help us probe the solar atmosphere at different heights and understand the mass and energy transfer process between its layers. For the first time, SUIT will also help us measure spatially resolved solar spectral irradiance at this wavelength band, which is significant for studying the sun-climate relationships. To perform these studies, it is necessary to photometrically calibrate the payload and validate the spectral coverage of the various bandpasses. We perform the photometric calibration and spectral validation of 8 bandpasses using light of known intensity and spectral coverage. For photometric calibration, the telescope throughput is modeled using sun-as-a-star spectrum from SOLSTICE and SOLSPEC. The modeled throughput is compared with in-lab measurements taken with light of known intensity. The ratio of measured photoelectrons gathered with the modeled prediction agree within 20%. For spectral validation, readings are taken across the transmission spectrum of each filter, keeping adjacent readings independent of each other. The relative intensity measured at each wavelength is seen to trace the modeled telescope bandpass for that filter. These tests could not be performed for filters with bandpasses operating below 250 nm (NB01, BB01 and BB02), primarily due to heavy atmospheric attenuation in these wavelengths leading to decreased SNR of the data. The experimentally measured results agree closely with the modeled values, validating SUIT's optical performance and presenting the reliability of the developed throughput model., Comment: 21 Pages, 8 Figures, 4 Tables

The R2D2 Deep Neural Network (DNN) series was recently introduced for image formation in radio interferometry. It can be understood as a learned version of CLEAN, whose minor cycles are substituted with DNNs. We revisit R2D2 on the grounds of series convergence, training methodology, and DNN architecture, improving its robustness in terms of generalisability beyond training conditions, capability to deliver high data fidelity, and epistemic uncertainty. Firstly, while still focusing on telescope-specific training, we enhance the learning process by randomising Fourier sampling integration times, incorporating multi-scan multi-noise configurations, and varying imaging settings, including pixel resolution and visibility-weighting scheme. Secondly, we introduce a convergence criterion whereby the reconstruction process stops when the data residual is compatible with noise, rather than simply using all available DNNs. This not only increases the reconstruction efficiency by reducing its computational cost, but also refines training by pruning out the data/image pairs for which optimal data fidelity is reached before training the next DNN. Thirdly, we substitute R2D2's early U-Net DNN with a novel architecture (U-WDSR) combining U-Net and WDSR, which leverages wide activation, dense connections, weight normalisation, and low-rank convolution to improve feature reuse and reconstruction precision. As previously, R2D2 was trained for monochromatic intensity imaging with the Very Large Array (VLA) at fixed $512 \times 512$ image size. Simulations on a wide range of inverse problems and a case study on real data reveal that the new R2D2 model consistently outperforms its earlier version in image reconstruction quality, data fidelity, and epistemic uncertainty., Comment: 17 pages, 6 figures

The Compton Spectrometer and Imager (COSI) is a Compton telescope designed to survey the 0.2-5 MeV sky, consisting of a compact array of cross-strip germanium detectors. As part of its development, in 2016 COSI had a successful 46 day flight on board NASA's Super Pressure Balloon platform. This was a precursor to the COSI Small Explorer (COSI-SMEX) satellite mission that will launch in 2027 into a equatorial low Earth (530 km) orbit. The observation of MeV gamma-rays is dominated by background radiation, especially due to the activation of the detector materials induced by cosmic-ray interactions. Thus, background simulation and identification are crucial for the data analysis. Because the COSI-SMEX detectors will be similar to the ones used for the balloon flight, the balloon measurements provide an important tool for testing and cross-checking our background simulations for the upcoming space mission. In this work we perform Monte Carlo simulations of the background emission from the 2016 COSI balloon flight. Including a phenomenological shape correction, we obtain an agreement with the data at the 10-20% level for energies between 0.1-1.6 MeV, and we successfully reproduce most of the activation lines induced by cosmic ray interactions., Comment: submitted to ApJ the 3 March 2025

Galactic chemical abundances provide crucial insights into fundamental galactic parameters, such as the high-mass slope of the initial mass function (IMF) and the normalization of Type Ia supernova (SN Ia) rates. Constraining these parameters is essential for advancing our understanding of stellar feedback, metal enrichment, and galaxy formation processes. However, traditional Bayesian inference techniques, such as Hamiltonian Monte Carlo (HMC), are computationally prohibitive when applied to large datasets of modern stellar surveys. We leverage simulation-based-inference (SBI) as a scalable, robust, and efficient method for constraining galactic parameters from stellar chemical abundances and demonstrate its the advantages over HMC in terms of speed, scalability, and robustness against model misspecifications. We combine a Galactic Chemical Evolution (GCE) model, CHEMPY, with a neural network emulator and a Neural Posterior Estimator (NPE) to train our SBI pipeline. Mock datasets are generated using CHEMPY, including scenarios with mismatched nucleosynthetic yields, with additional tests conducted on data from a simulated Milky Way-like galaxy. SBI results are benchmarked against HMC-based inference, focusing on computational performance, accuracy, and resilience to systematic discrepancies. SBI achieves a $\sim75,600\times$ speed-up compared to HMC, reducing inference runtime from $\gtrsim42$ hours to mere seconds for thousands of stars. Inference on $1,000$ stars yields precise estimates for the IMF slope ($\alpha_{\rm IMF} = -2.298 \pm 0.002$) and SN Ia normalization ($\log_{10}(N_{\rm Ia}) = -2.885 \pm 0.003$), deviating less than 0.05% from the ground truth. SBI also demonstrates similar robustness to model misspecification than HMC, recovering accurate parameters even with alternate yield tables or data from a cosmological simulation. (shortened...), Comment: submitted to A&A, comments welcome, all source code to reproduce this work can be found on GitHub under url: https://github.com/TobiBu/sbi-chempy

LOFAR is a low-frequency array distributed across several European countries. Each LOFAR station contains thousands of antennas and associated electronics, making monitoring and thorough testing of those components essential to ensuring station reliability. This paper discusses various anomalies that may arise in LOFAR antennas, tile elements, modems, and summators. We also introduce two diagnostic pipelines designed to detect these anomalies: a real-time station monitoring system and an offline stationtest system. These pipelines provide valuable insights into the operational status of each antenna, issuing alerts to minimize observational disruptions while maximizing station uptime, reliability, and sensitivity. By enhancing the efficiency and stability of LOFAR stations, they also serve as a foundation for future large-scale arrays like SKA-Low. The experience gained from their development and deployment will contribute to the construction and maintenance of SKA-Low, improving monitoring and diagnostic capabilities for large-scale antenna networks. Ultimately, these systems play a crucial role in ensuring continuous observations and maintaining data integrity.

Point spread functions (PSFs) describe the distribution of light for a pure point source in an astronomical image due to the instrument optics. For deconvolution, as for point source photometry and for source removal, it is key to have an accurate PSF for a particular image. Space-based telescopes can then pose a challenge as their PSFs are informed by their complex construction, and the myriad of pointings and rotations used to capture deep images. These telescopes also capture the highest resolution images of astronomical sources, resolving stars around even relatively distant galaxies. Proper co-addition of PSFs at a specific source position for space-based imaging is then both critical and challenging. This code, spike, generates model PSFs and runs them through the same processing pipeline used to derive deep, co-added images, providing correctly co-added and resampled PSFs for images from the Hubble Space Telescope, the James Webb Space Telescope, and the Nancy Grace Roman Space Telescope., Comment: Submitted to JOSS; code available at https://github.com/avapolzin/spike

The quest for primordial $B$-mode polarization signatures in the Cosmic Microwave Background (CMB) is a major goal of contemporary cosmology. Detecting these signatures would confirm primordial gravitational waves and allow precise determination of the tensor-to-scalar ratio, $r$, which is crucial for distinguishing between inflationary models. This requires high-precision, full-sky observations from space-based platforms to avoid atmospheric fluctuations and windows. Since $B$-mode signatures are much fainter than CMB temperature anisotropies, their detection requires well-designed calibration and mitigation strategies. Next-generation space observatories such as LiteBIRD, which use Half-Wave Plate (HWP) modulation technology, represent a significant advance. This technology allows single-detector observations, eliminating the need for differential detection and its associated systematic complexities. The study begins with the optimization of scan strategy parameters for missions with HWP to improve in-flight calibration, suppress systematic effects, and develop robust null-test methods. In addition, we present a method for calculating the impact of systematic effects on the estimation of $r$ using the SBM software framework, which performs fast mapping in spin space. This method allows the decomposition and elimination of systematic effects in terms of spins, and we demonstrate the elimination of several systematic effects., Comment: Ph.D. thesis. Abstracts are abridged for arXiv submission

In the era of time-domain, multi-messenger astronomy, the detection of transient events on the high-energy electromagnetic sky has become more important than ever. The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a dedicated mission to monitor gamma-ray transients, launched in December, 2020. A real-time on-board trigger and location software, using the traditional signal-to-noise ratio (SNR) method for blind search, is constrained to relatively bright signals due to the limitations in on-board computing resources and the need for real-time search. In this work, we developed a ground-based pipeline for GECAM to search for various transients, especially for weak bursts missed by on-board software. This pipeline includes both automatic and manual mode, offering options for blind search and targeted search. The targeted search is specifically designed to search for interesting weak bursts, such as gravitational wave-associated gamma-ray bursts (GRBs). From the ground search of the data in the first year, GECAM has been triggered by 54 GRBs and other transients, including soft gamma-ray repeaters, X-ray binaries, solar flares, terrestrial gamma-ray flashes. We report the properties of each type of triggers,such as trigger time and light curves. With this search pipeline and assuming a soft Band spectrum, the GRB detection sensitivity of GECAM is increased to about 1.1E-08 erg cm-2 s-1 (10 keV - 1000 keV, burst duration of 20 s). These results demonstrate that the GECAM ground search system (both blind search and targeted search) is a versatile pipeline to recover true astrophysical signals which were too weak to be found in the on-board search., Comment: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy (SCPMA)

Understanding the ages of stars is crucial for unraveling the formation history and evolution of our Galaxy. Traditional methods for estimating stellar ages from spectroscopic data often struggle with providing appropriate uncertainty estimations and are severely constrained by the parameter space. In this work, we introduce a new approach using normalizing flows, a type of deep generative model, to estimate stellar ages for evolved stars with improved accuracy and robust uncertainty characterization. The model is trained on stellar masses for evolved stars derived from asteroseismology and predicts the relationship between the carbon and nitrogen abundances of a given star and its age. Unlike standard neural network techniques, normalizing flows enable the recovery of full likelihood distributions for individual stellar ages, offering a richer and more informative perspective on uncertainties. Our method yields age estimations for 378,720 evolved stars and achieves a typical absolute age uncertainty of approximately 2 Gyr. By intrinsically accounting for the coverage and density of the training data, our model ensures that the resulting uncertainties reflect both the inherent noise in the data and the completeness of the sampled parameter space. Applying this method to data from the SDSS-V Milky Way Mapper, we have produced the largest stellar age catalog for evolved stars to date., Comment: Submitted to AJ. Age data to be released with SDSS DR19

We propose a Roman Space Telescope survey to investigate fundamental properties of the distant solar system in the region of the Kuiper Belt where object characteristics and the size distribution are inaccessible from any other telescope. Our pointing is coincident with the search space accessible to NASA's New Horizons spacecraft meaning, that a discovered object sufficiently near the orbit of New Horizons would potentially be investigated by a close flyby. In addition, numerous objects expected to be discovered by this search can be observed in the distance by New Horizons allowing their surface properties and satellite systems to both be probed. As designed, this survey will discover and determine orbits for as many as 900 Kuiper Belt objects (KBOs), providing a unique opportunity for ground-breaking Kuiper Belt science. It will simultaneously: (1) Probe and characterize the deep Kuiper Belt by identifying objects as small as a few km and taking our understanding of the size distribution to a new level. This has implications for understanding the the standard model (the Streaming Instability) of KBO formation and elucidating crater formation physics on these icy bodies. (2) Open KBO rotation studies, in particular of those objects with long rotation periods,(3) Discover and characterize KBO binaries at large distances, important because their duplicity offers information about object densities at these distant locations from the Sun. (4) Shed light on the cratering history of KBOs and improving the dating of the surfaces of Arrokoth, Pluto and Charon in addition to helping to place the 32 distant KBOs New Horizons has observed in context. This project also has synergies with transiting exoplanet studies due to the stellar density of our search fields. Coupled with our timing requirements it is sensitive to discovery of hot Jupiters and hot Neptunes., Comment: 12 Pages, 6 Figures, 1 Table. Submitted to Roman Galactic Plane Survey White paper call

We present the design and science case for a new array of radio antennas to be located at the Pierre Auger Observatory. Six stations of three SKALA antennas each will be deployed around a single water-Cherenkov surface detector triggering the radio readout. The planned antenna layout will allow for the detection of cosmic rays above a few tens of PeV and reach full efficiency for vertical air showers at several hundred PeV in primary energy. The array will thus be a pathfinder to demonstrate that fully-efficient radio detection in combination with the underground muon detectors already present at the location is possible. This will enable combined studies of the muon component and the depth of shower maximum, relevant for hadronic interaction models studies and more accurate determination of the cosmic-ray mass composition in the energy range of the Galactic-to-extragalactic transition., Comment: Presented at the 7th International Symposium on Ultra High Energy Cosmic Rays (UHECR2024)

We present brutus, an open source Python package for quickly deriving stellar properties, distances, and reddenings to stars based on grids of stellar models constrained by photometric and astrometric data. We outline the statistical framework for deriving these quantities, its implementation, and various Galactic priors over the 3-D distribution of stars, stellar properties, and dust extinction (including $R_V$ variation). We establish a procedure to empirically calibrate MIST v1.2 isochrones by using open clusters to derive corrections to the effective temperatures and radii of the isochrones, which reduces systematic errors on the lower main sequence. We also describe and apply a method to estimate photometric offsets between stellar models and observed data using nearby, low-reddening field stars. We perform a series of tests on mock and real data to examine parameter recovery with MIST under different modeling assumptions, illustrating that brutus is able to recover distances and other stellar properties using optical to near-infrared photometry and astrometry. The code is publicly available at https://github.com/joshspeagle/brutus., Comment: 33 pages, 17 figures main text (+23 pages, +10 figures in appendices). To be re-submitted to AAS Journals. Package documentation to be added in coming weeks (as of 2025-03-03)

When a gravitational wave signal encounters a massive object, such as a galaxy or galaxy cluster, it undergoes strong gravitational lensing, producing multiple copies of the original signal. These strongly lensed signals exhibit identical waveform morphology in the frequency domain, allowing analysis without the need for complex lens models. However, stellar fields and dark matter substructures within the galactic lens introduce microlensing effects that alter individual signal morphologies. Identifying these microlensing signatures is computationally challenging within Bayesian frameworks. In this study, we propose a novel residual test to efficiently search for microlensing signatures by leveraging the fact that current Bayesian inference pipelines are optimized solely for the strong lensing hypothesis. Using cross-correlation techniques, we investigate the microlensing-induced deviations from the strong hypothesis, which are imprinted in the residuals. Our simulated signals from realistic microlensing populations reveal that while most events exhibit small mismatches, a fraction exhibit significant deviations. We find that 31% (55%) and 45% (65%) of microlensed events with mismatch > 0.03 and > 0.1, respectively, can be discerned with O4 (O5) detector sensitivities, which demonstrates that high-mismatch events are more likely to be identified as microlensed. Including all events from a realistic population, 12% (21.5%) are identifiable with O4 (O5) sensitivity using our approach., Comment: 13 pages, 6 figures, 2 tables

Recently, the R2D2 paradigm, standing for ''Residual-to-Residual DNN series for high-Dynamic-range imaging'', was introduced for image formation in Radio Interferometry (RI) as a learned version of the traditional algorithm CLEAN. The first incarnations of R2D2 are limited to planar imaging on small fields of view, failing to meet the spherical-imaging requirement of modern telescopes observing wide fields. To address this limitation, we propose the spherical-imaging extension S-R2D2. Firstly, as R2D2, S-R2D2 encapsulates its minor cycles in existing 2D-Euclidean deep neural network (DNN) architectures, but adapts its iterative scheme to incorporate the wide-field measurement model mapping a spherical image to visibility data. We implemented this model as the composition of an efficient Fourier-based interpolator mapping the spherical image onto the equatorial plane, with the standard RI operator mapping the equatorial-plane image to visibility data. Importantly, the interpolation step must inevitably be performed at a lower-than-optimal resolution on the plane, to meet the high-resolution requirement on the sphere of wide-field imaging while preserving scalability. Therefore, secondly, we design S-R2D2's DNN training loss to jointly learn to correct the interpolation approximations and identify residual image structures on the sphere, ensuring consistency with the spherical ground truth using the adjoint plane-to-sphere interpolator. Finally, we demonstrate through simulations S-R2D2's capability to perform fast and accurate reconstructions of spherical monochromatic intensity images, across high-resolution, high-dynamic-range settings., Comment: 16 pages, 13 figures