Celestial Sources Research Articles

Intensity-sensitive Quality Assessment of Extended Sources in Astronomical Images

Radio astronomy studies the Universe by observing the radio emissions of celestial bodies. Different methods can be used to recover the sky brightness distribution (SBD), which describes the distribution of celestial sources from recorded data, with the output dependent on the method used. Image quality assessment (IQA) indexes can be used to compare the differences between restored SBDs produced by different image reconstruction techniques to evaluate their effectiveness. However, reconstructed images (for the same SBD) can appear to be very similar, especially when observed by the human visual system (HVS). Hence, current structural similarity methods, inspired by the HVS, are not effective. In the past, we have proposed two methods to assess point-source images, where low amounts of concentrated information are present in larger regions of noise-like data. But for images that include extended source(s), the increase in complexity of the structure makes the IQA methods for point sources oversensitive because the important objects cannot be described by isolated point sources. Therefore, in this article we propose the augmented low-information similarity index (augLISI), an improved version of LISI, to assess images including extended source(s). Experiments have been carried out to illustrate how this new IQA method can help with the development and study of astronomical imaging techniques. Note that although we focus on radio astronomical images herein, these IQA methods are also applicable to other astronomical images and imaging techniques.

Planetary radio interferometry and Doppler experiment as an operational component of the Jupiter Icy Moons Explorer mission

AbstractWe present an overview of the operations and engineering interface for Planetary Radio Interferometry and Doppler Experiment (PRIDE) radio astronomy observations as a scientific component of the ESA’s Jupiter Icy Moons Explorer (JUICE) mission, as well as other prospective planetary and space science missions. The article discusses advanced scheduling and planning methods that make it possible to create observing schedules for observations of specific spacecraft in concurrence with observations of natural radio sources. In order to put this into practice and find suitable natural background calibrator sources for PRIDE of JUICE mission, we developed planning and scheduling software. The conventional scheduling software for natural celestial radio sources is not set up to include spacecraft as observation targets in the necessary control files. Therefore, difficulties already arise during observation planning. We report on the development of new and the adaptation of existing routines used in astrophysical and geodetic VLBI for satellite scheduling and planning. The analysis of the PRIDE science observations led to improved observational planning, and the mission’s scheduling methodologies were studied using a systems engineering approach. In addition, we highlighted the new procedures, like finding charts for selecting calibrator radio sources over a range of frequency bands and the outcomes of those strategies for science operation planning. A simulation of the flyby of Venus during the cruise phase of the JUICE spacecraft, based on the Tudat software, is also presented, resulting in a promising opportunity to test PRIDE techniques and evaluate the improvements that PRIDE observables can make to natural bodies’ ephemerides. The first $$\hbox {K}_{a}$$ K a -band (32 GHz) observations of the ESA’s BepiColombo by a radio telescope in the VLBI network, which employs a similar radio communications system as JUICE, were also demonstrated as a test case. The primary objective of these activities is to serve as a practice run for the upcoming operational PRIDE JUICE operations. We showcase the capabilities of the planning and scheduling software for other space missions.

High-temporal-resolution optical spectroscopic observations of the transitional millisecond pulsar PSR J1023+0038

Transitional millisecond pulsars (tMSPs) represent a dynamic category of celestial sources that establish a crucial connection between low-mass X-ray binaries and millisecond radio pulsars. These systems exhibit transitions from rotation-powered states to accretion-powered ones and vice versa, highlighting the tight evolutionary link expected by the so-called recycling scenario. In their active phase, these sources manifest two distinct emission modes named high and low, occasionally punctuated by sporadic flares. Here, we present high-time-resolution spectroscopic observations of the binary tMSP J1023+0038, in the sub-luminous disc state. This is the first short-timescale (sim 1 min) optical spectroscopic campaign ever conducted on a tMSP. The campaign was carried out over the night of June 10, 2021 using the Gran Telescopio Canarias. The optical continuum shows erratic variability, without clear evidence of high and low modes or of orbital modulation. Besides, the analysis of these high-temporal-cadence spectroscopic observations reveals, for the first time, evidence for a significant (up to a factor of $ 2$) variability in the emission line properties (equivalent width and full width half maximum) over a timescale of minutes. Intriguingly, the variability episodes observed in the optical continuum and in the emission line properties seem uncorrelated, making their origin unclear.

LSPE-Strip on-sky calibration strategy using bright celestial sources

In this paper we describe the global on-sky calibration strategy of the LSPE-Strip instrument. Strip is a microwave telescope operating in the Q- and W-bands (central frequencies of 43 and 95 GHz respectively) from the Observatorio del Teide in Tenerife, with the goal to observe and characterise the polarised Galactic foreground emission, and complement the observations of the polarisation of the cosmic microwave background to be performed by the LSPE-SWIPE instrument and other similar experiments operating at higher frequencies to target the detection of the B-mode signal from the inflationary epoch of the Universe. Starting from basic assumptions on some of the instrument parameters (NET, 1/f noise knee frequency, beam properties, observing efficiency) we perform realistic simulations to study the level of accuracy that can be achieved through observations of bright celestial calibrators in the Strip footprint (sky fraction of 30%) on the determination and characterisation of the main instrument parameters: global and relative gain factors (in intensity and in polarisation), polarisation direction, polarisation efficiency, leakage from intensity to polarisation, beams, window functions and pointing model.

Mitigating the effect of source structure in geodetic VLBI by re-weighting observations using closure delays and baseline-to-jet orientation

An ideal target for geodetic very long baseline interferometry (VLBI) is a strong and point-like radio source. In reality, most celestial sources used in geodetic VLBI have spatial structure. This is as a major source of error in VLBI Global Observing System (VGOS) and also affects legacy S/X observations. Source structure causes a systematic delay, which can affect the geodetic estimates if not modelled or otherwise accounted for. In this work, we aim to mitigate its impact by extending the stochastic model used in the least-squares fitting of the VLBI group delays. We have developed a weighting scheme to re-weight the observations by parameterizing the source structure component in terms of closure delays and jet orientation relative to the observing baseline. It was implemented in the Vienna VLBI Software. To assess the performance of the extended stochastic model, we analysed the CONT17 legacy sessions and generated suitable reference solutions for comparison. The effects of re-weighting were evaluated with respect to the session fit statistics, source-wise residuals, and geodetic parameters. We find that this relatively simple noise model consistently improves the session fit by about 5% with moderate variation from session to session. The geodetic estimates are not affected to a significant level by this new weighting method. Source-wise we see improved post-fit residuals for 63 out of a total of 91 sources observed.

Future Perspectives for Gamma-ray Burst Detection from Space

Since their first discovery in the late 1960s, gamma-ray bursts have attracted an exponentially growing interest from the international community due to their central role in the most highly debated open questions of the modern research of astronomy, astrophysics, cosmology, and fundamental physics. These range from the intimate nuclear composition of high-density material within the core of ultra-dense neuron stars, to stellar evolution via the collapse of massive stars, the production and propagation of gravitational waves, as well as the exploration of the early universe by unveiling the first stars and galaxies (assessing also their evolution and cosmic re-ionization). GRBs in the past ∼50 years have stimulated the development of cutting-edge technological instruments for observations of high-energy celestial sources from space, leading to the launch and successful operations of many different scientific missions (several of them still in data-taking mode currently). In this review, we provide a brief description of the GRB-dedicated missions from space being designed and developed for the future. The list of these projects, not meant to be exhaustive, shall serve as a reference to interested readers to understand what is likely to come next to lead the further development of GRB research and the associated phenomenology.

The Imaging X-ray Polarimetry Explorer (IXPE) and New Directions for the Future

An observatory dedicated to X-ray polarimetry has been operational since 9 December 2021. The Imaging X-ray Polarimetry Explorer (IXPE), a collaboration between NASA and ASI, features three X-ray telescopes equipped with detectors sensitive to linear polarization set to 120°. This marks the first instance of a three-telescope SMEX mission. Upon reaching orbit, an extending boom was deployed, extending the optics and detector to a focal length of 4 m. IXPE targets each celestial source through dithering observations. This method is essential for supporting on-ground calibrations by averaging the detector’s response across a section of its sensitive plane. The spacecraft supplies power, enables attitude determination for subsequent on-ground attitude reconstruction, and issues control commands. After two years of observation, IXPE has detected significant linear polarization from nearly all classes of celestial sources emitting X-rays. This paper outlines the IXPE mission’s achievements after two years of operation in orbit. In addition, we report developments for future high-throughput X-ray optics that will have much smaller dead-times by using a new generation of Applied Specific Integrated Circuits (ASIC), and may provide 3D reconstruction of photo-electron tracks.

Scientific Highlights of the AGILE Gamma-ray Mission

The γ-ray sky above a few tens of megaelectronvolts (MeV) reveals some of the most powerful and energetic phenomena of our Universe. The Astrorivelatore Gamma ad Immagini LEggero (AGILE) Gamma-ray Mission was launched in 2007 with the aim of observing celestial sources by means of three instruments covering a wide range of energies, from hard X-rays up to 30 GeV. Thanks to its wide field of view, AGILE set to observe and detect emission from pulsars, pulsar wind nebulae, gamma-ray bursts, active galactic nuclei, fast radio bursts, terrestrial gamma-ray flashes, and the electromagnetic counterparts of neutrinos and gravitational waves. In particular, the fast on-ground processing and analysis chain allowed the AGILE team to promptly respond to transient events, and activate or participate in multiwavelength observing campaigns. Eventually, after 17 years of operations, the AGILE Italian scientific satellite re-entered the atmosphere on 14 February 2024, ending its intense activity as a hunter of some of the most energetic cosmic sources in the Universe that emit X and γ-rays. We will review the most relevant AGILE results to date and their impact on the advancements of theoretical models.

A Hierarchical Method for Locating the Interferometric Fringes of Celestial Sources in the Visibility Data

In source detection in the Tianlai project, locating the interferometric fringe in visibility data accurately will influence downstream tasks drastically, such as physical parameter estimation and weak source exploration. Considering that traditional locating methods are time-consuming and supervised methods require a great quantity of expensive labeled data, in this paper, we first investigate characteristics of interferometric fringes in the simulation and real scenario separately, and integrate an almost parameter-free unsupervised clustering method and seeding filling or eraser algorithm to propose a hierarchical plug and play method to improve location accuracy. Then, we apply our method to locate single and multiple sources’ interferometric fringes in simulation data. Next, we apply our method to real data taken from the Tianlai radio telescope array. Finally, we compare with unsupervised methods that are state of the art. These results show that our method has robustness in different scenarios and can improve location measurement accuracy effectively.

Science with the ASTRI Mini-Array: From Experiment to Open Observatory

Although celestial sources emitting in the few tens of GeV up to a few TeV are being investigated by imaging atmospheric Čerenkov telescope arrays such as H.E.S.S., MAGIC, and VERITAS, at higher energies, up to PeV, more suitable instrumentation is required to detect ultra-high-energy photons, such as extensive air shower arrays, as HAWC, LHAASO, Tibet AS-γ. The Italian National Institute for Astrophysics has recently become the leader of an international project, the ASTRI Mini-Array, with the aim of installing and operating an array of nine dual-mirror Čerenkov telescopes at the Observatorio del Teide in Spain starting in 2025. The ASTRI Mini-Array is expected to span a wide range of energies (1–200 TeV), with a large field of view (about 10 degrees) and an angular and energy resolution of ∼3 arcmin and ∼10 %, respectively. The first four years of operations will be dedicated to the exploitation of Core Science, with a small and selected number of pointings with the goal of addressing some of the fundamental questions on the origin of cosmic rays, cosmology, and fundamental physics, the time-domain astrophysics and non γ-ray studies (e.g., stellar intensity interferometry and direct measurements of cosmic rays). Subsequently, four more years will be dedicated to Observatory Science, open to the scientific community through the submission of observational proposals selected on a competitive basis. In this paper, I will review the Core Science topics and provide examples of possible Observatory Science cases, taking into account the synergies with current and upcoming observational facilities.

Low Complexity Radio Frequency Interference Mitigation for Radio Astronomy Using Large Antenna Array

With the ongoing growth in radio communications, there is an increased contamination of radio astronomical source data, which hinders the study of celestial radio sources. In many cases, fast mitigation of strong radio frequency interference (RFI) is valuable for studying short lived radio transients so that the astronomers can perform detailed observations of celestial radio sources. The standard method to manually excise contaminated blocks in time and frequency makes the removed data useless for radio astronomy analyses. This motivates the need for better RFI mitigation techniques for array of size M antennas. Although many solutions for mitigating strong RFI improves the quality of the final celestial source signal, many standard approaches require all the eigenvalues of the spatial covariance matrix ([Formula: see text]) of the received signal, which has [Formula: see text] computation complexity for removing RFI of size d where [Formula: see text]. In this work, we investigate two approaches for RFI mitigation, (1) the computationally efficient Lanczos method based on the Quadratic Mean to Arithmetic Mean (QMAM) approach using information from previously-collected data under similar radio-sky-conditions, and (2) an approach using a celestial source as a reference for RFI mitigation. QMAM uses the Lanczos method for finding the Rayleigh–Ritz values of the covariance matrix [Formula: see text], thus, reducing the computational complexity of the overall approach to [Formula: see text]. Our numerical results, using data from the radio observatory Long Wavelength Array (LWA-1), demonstrate the effectiveness of both proposed approaches to remove strong RFI, with the QMAM-based approach still being computationally efficient.

A Multimodal Transfer Learning Method for Classifying Images of Celestial Point Sources

A large fraction of celestial objects exhibit point shapes in CCD images, such as stars and QSOs, which contain less information due to their few pixels. Point source classification based solely on image data may lead to low accuracy. To address this challenge, this paper proposes a Multi-modal Transfer Learning-based classification method for celestial objects with point shape images. Considering that spectral data possess rich features and that there is a correlation between spectral data and image data, the proposed approach fully utilizes the knowledge gained from celestial spectral data and transfers it to the original image-based classification, enhancing the accuracy of classifying stars and QSOs. Initially, a one-dimensional residual network is employed to extract a 128-dimensional spectral feature vector from the original 3700-dimensional spectral data. This spectral feature vector captures important features of the celestial object. The Generative Adversarial Network is then utilized to generate a simulated spectral vector of 128 dimensions, which corresponds to the celestial object image. By generating simulated spectral vectors, data from two modals (spectral and image) for the same celestial object are available, enriching the input features of the model. In the upcoming multimodal classification model, we only require the images of celestial objects along with their corresponding simulated spectral data, and we no longer need real spectral data. With the assistance of spectral data, the proposed method alleviates the above disadvantages of the original image-based classification method. Remarkably, our method has improved the F1-score from 0.93 to 0.9777, while reducing the error rate in classification by 40%. These enhancements significantly increase the classification accuracy of stars and QSOs, providing strong support for the classification of celestial point sources.

A sign in space: An interdisciplinary exploration of the potential reception of an extraterrestrial signal

A Sign in Space is an interdisciplinary project by media artist Daniela de Paulis, in collaboration with the Green Bank Observatory, the National Institute for Astrophysics, the SETI Institute and the European Space Agency.The project consists in transmitting a simulated extraterrestrial message as part of a live performance, using an ESA spacecraft as celestial source. The objective of the project is to involve the worldwide search for extraterrestrial intelligence (SETI) community, professionals from different fields and the broader public in the reception, decoding and interpretation of the message. This process will require global cooperation, bridging a conversation around the topics of SETI, space research and society, across multiple cultures and fields of expertise.

Digging a little deeper: characterizing three new extreme ULX candidates

ABSTRACT A prime motivation for compiling catalogues of any celestial X-ray source is to increase our numbers of rare subclasses. In this work, we take a recent multimission catalogue of ultraluminous X-ray sources (ULXs) and look for hitherto poorly-studied ULX candidates that are luminous ($L_{\rm X} \ge 10^{40}\,{\rm erg\,s}^{-1}$), bright ($f_{\rm X} \ge 5 \times 10^{-13}\,{\rm erg\,cm}^{-2}\,{\rm s}^{-1}$), and have archival XMM–Newton data. We speculate that this luminosity regime may be ideal for identifying new pulsating ULXs (PULXs), given that the majority of known PULXs reach similar high luminosities. We find three sources that match our criteria and study them using archival data. We find 4XMM J165251.5−591503 to possess a bright and variable Galactic optical/IR counterpart, and so conclude it is very likely to be a foreground interloper. 4XMM J091948.8−121429 does appear to be an excellent ULX candidate associated with the dwarf irregular galaxy PGC 26378, but has only one detection to date with low data quality. The best data set belongs to 4XMM J112054.3+531040 which we find to be a moderately variable, spectrally hard (Γ ≈ 1.4) X-ray source located in a spiral arm of NGC 3631. Its spectral hardness is similar to known PULXs, but no pulsations are detected by accelerated pulsation searches in the available data. We discuss whether other missions provide objects for similar studies and compare this method to others suggested for identifying good PULX candidates.

Likelihood of the sky

Catalogs with celestial positions and uncertainties are among the primary data products of observational astronomy. Consequently, the (cross-)matching problem, namely, the question of how many (and which) pairs of positions pertain to the same celestial source, remains a central part of this scientific exploitation. We address the matching problem for two catalogs from a purely geometric point of view by adopting the concept of point processes and nearest-neighbor distributions. Suitable expressions for the nearest-neighbor distributions were derived and used to model the distribution of spatial offsets in the sky. Thus, we are able to estimate a total number of matching pairs, along with individual matching probabilities between sources in any two catalogs of positions and uncertainties. We demonstrate the workings of the model using mock data and apply it to cross-matching the X-ray sources of the eROSITA Final Equatorial-Depth survey with opticalGaiacounterparts.

Design of Extensible Structured Interferometric Array Utilizing the “Coarray” Concept

The optimum placement of receiving telescope antennas is a central topic for designing radio interferometric arrays, and this determines the performance of the obtained information. A variety of arrays are designed for different purposes, and they perform poorly in scalability. In this paper, we consider a subclass of structured sparse arrays, namely nested arrays, and examine the important role of “coarray” in interferometric synthesis imaging, which is utilized to design nested array configurations for a complete uniform Fourier plane coverage in both supersynthesis and instantaneous modes. Both nested arrays and the theory of the coarray have rich research achievements, and we apply them to astronomy to design arrays with good scalability and imaging performance. Simulated celestial source image retrieval performance validates the effectiveness of nested interferometric arrays.

QUIJOTE scientific results – IV. A northern sky survey in intensity and polarization at 10–20 GHz with the multifrequency instrument

ABSTRACT We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17, and 19 GHz, and covering approximately 29 000 deg2, including most of the northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE multifrequency instrument (MFI), and have angular resolutions of around 1°, and sensitivities in polarization within the range 35–40 µK per 1° beam, being a factor ∼2–4 worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterize the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5 per cent. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A, and 3C274), which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11 GHz and WMAP 23 GHz map peaks at β = −3.09 with a standard deviation of 0.14. The measured BB/EE ratio at scales of ℓ = 80 is 0.26 ± 0.07 for a Galactic cut |b| > 10°. We find a positive TE correlation for 11 GHz at large angular scales (ℓ ≲ 50), while the EB and TB signals are consistent with zero in the multipole range 30 ≲ ℓ ≲ 150. The maps discussed in this paper are publicly available.

Accounting for systematic uncertainties in the Imaging X-ray Polarimetry Explorer (IXPE) detector response

Launched on December 9, 2021, the Imaging X-ray Polarimetry Explorer (IXPE) is the first imaging polarimeter ever flown, providing sensitivity in the 2–8 keV range, and during the 2-year initial phase of the mission will sample tens of X-ray sources among different source classes. While most of the measurements will be statistics-limited, for some of the brightest objects observed and long integration times, the systematic uncertainties in the detector response (primarily the effective area, the modulation factor and the absolute energy scale) will be important. In this contribution, we describe a framework to propagate on high-level observables (e.g.: spectro-polarimetric fit parameters) the systematic uncertainties connected with the response of the detector, that we estimate from relevant ground calibrations and from observations of celestial point sources.

Semi-supervised classification and clustering analysis for variable stars

ABSTRACT The immense amount of time series data produced by astronomical surveys has called for the use of machine learning algorithms to discover and classify several million celestial sources. In the case of variable stars, supervised learning approaches have become commonplace. However, this needs a considerable collection of expert-labelled light curves to achieve adequate performance, which is costly to construct. To solve this problem, we introduce two approaches. First, a semi-supervised hierarchical method, which requires substantially less trained data than supervised methods. Second, a clustering analysis procedure that finds groups that may correspond to classes or subclasses of variable stars. Both methods are primarily supported by dimensionality reduction of the data for visualization and to avoid the curse of dimensionality. We tested our methods with catalogues collected from the Optical Gravitational Lensing Experiment (OGLE), the Catalina Sky Survey (CSS), and the Gaia survey. The semi-supervised method reaches a performance of around 90 per cent for all of our three selected catalogues of variable stars using only $5{{\ \rm per\ cent}}$ of the data in the training. This method is suitable for classifying the main classes of variable stars when there is only a small amount of training data. Our clustering analysis confirms that most of the clusters found have a purity over 90 per cent with respect to classes and 80 per cent with respect to subclasses, suggesting that this type of analysis can be used in large-scale variability surveys as an initial step to identify which classes or subclasses of variable stars are present in the data and/or to build training sets, among many other possible applications.

Weak celestial source fringes detection based on channel attention shrinkage networks and cluster-based anchor boxes generation algorithm

Detecting weak celestial source signals from massive radio data is a very challenging task because the radiation received by radio telescope is very weak and prone to disturbances. In order to detect these weak signals, we propose a two-stage object detection method that performs more finely in computer vision tasks. The novelty of the proposed method is to combine traditional soft thresholding denoising methods with attention mechanisms in deep neural networks. We propose a channel attention shrinkage network as the backbone of the object detection model to extract the features of weak signals from celestial sources by removing noise-related information. Moreover, targeting the characteristics of celestial source fringes in phase images, we propose a cluster-based anchor boxes generation algorithm to improve the accuracy of fringes position detection. We also introduce the CIoU loss function to improve the performance of the model because of the large aspect ratio of the celestial source fringes in the phase image. We generate simulated celestial source fringes data based on the parameters of the observation system to train our model and conduct experiments to evaluate the performance of the proposed algorithm. Our model obtains satisfactory detection accuracy and accurate for the location of celestial source fringes.