2025-01-17 | | Total: 9
In recent years, numerical simulations have become indispensable for addressing complex astrophysical problems. The MagnetoHydroDynamics (MHD) framework represents a key tool for investigating the dynamical evolution of astrophysical plasmas, which are described as a set of partial differential equations that enforce the conservation of mass, momentum, and energy, along with Maxwell's equation for the evolution of the electromagnetic fields. Due to the high nonlinearity of the MHD equations (regardless of their specifications, e.g., classical/relativistic or ideal/resistive), a general analytical solution is precluded, making the numerical approach crucial. Numerical simulations usually end up producing large sets of data files and their scientific analysis leans on dedicated software designed for data visualization. However, in order to encompass all of the code output features, specialized tools focusing on the numerical code may represent a more versatile and built-in tool. Here, we present PyPLUTO, a Python package tailored for efficient loading, manipulation, and visualization of outputs produced with the PLUTO code (Mignone et al., 2007; Mignone et al., 2012). PyPLUTO uses memory mapping to optimize data loading and provides general routines for data manipulation and visualization. PyPLUTO also supports the particle modules of the PLUTO code, enabling users to load and visualize particles, such as cosmic rays (Mignone et al., 2018), Lagrangian (Vaidya et al., 2018), or dust (Mignone et al., 2019) particles, from hybrid simulations. A dedicated Graphical User Interface (GUI) simplifies the generation of single-subplot figures, making PyPLUTO a powerful yet user-friendly toolkit for astrophysical data analysis.
Low-energy (<300 keV) protons entering the field of view of XMM-Newton are observed in the form of a sudden increase in the background level, the so-called soft proton flares, affecting up to 40% of the mission observing time. In-flight XMM-Newton's observations of soft protons represent a unique laboratory to validate and improve our understanding of their interaction with the mirror, optical filters, and X-ray instruments. At the same time, such models would link the observed background flares to the primary proton population encountered by the telescope, converting XMM-Newton into a monitor for soft protons. We built a Geant4 simulation of XMM-Newton, including a verified mass model of the X-ray mirror, the focal plane assembly, and the EPIC MOS and pn-CCDs. We encoded the energy redistribution and proton transmission efficiency into a redistribution matrix file (RMF) and an auxiliary response file (ARF). For the validation, three averaged soft proton spectra, one for each filter configuration, were extracted from a collection of 13 years of MOS observations of the focused non X-ray background and analysed with Xspec. The best-fit model is in agreement with the power-law distribution predicted from independent measurements for the XMM-Newton orbit, spent mostly in the magnetosheath and nearby regions. For the first time we are able to link detected soft proton flares with the proton radiation environment in the Earth's magnetosphere, while proving the validity of the simulation chain in predicting the background of future missions. Benefiting from this work and contributions from the Athena instrument consortia, we also present the response files for the Athena mission and updated estimates for its focused charged background.
Medium-timescale (minutes to hours) radio transients are a relatively unexplored population. The wide field-of-view and high instantaneous sensitivity of instruments such as MeerKAT provides an opportunity to probe this class of sources, using image-plane detection techniques. The previous letter in this series describes our project and associated TRON pipeline designed to mine archival MeerKAT data for transient and variable sources. In this letter, we report on a new transient, a radio flare, associated with Gaia DR3 6865945581361480448, a G type star, whose parallax places it at a distance of 1334 pc. Its duration and high degree of circular polarization suggests electron cyclotron maser instability as the mechanism, consistent with an RS CVn variable.
Medium-timescale (minutes to hours) radio transients are a relatively unexplored population. The wide field-of-view and high instantaneous sensitivity of instruments such as MeerKAT provides an opportunity to probe this class of sources, using image-plane detection techniques. We aim to systematically mine archival synthesis imaging data in order to search for medium-timescale transients and variables that are not detected by conventional long-track image synthesis techniques. We deploy a prototype blind transient and variable search pipeline named TRON. This processes calibrated visibility data, constructs high-time cadence images, performs a search for variability on multiple timescales, and extracts lightcurves for detected sources. As proof of concept, we apply it to three MeerKAT observations of globular clusters, known to host transient or variable sources. We detect a previously known eclipsing MSP suspected to be a `black widow' system, in the globular cluster Omega Centauri, with a light curve confirming the eclipsing nature of the emission. We detect a previously known `red back' eclipsing MSP in the globular cluster Terzan 5. Using observations of the globular cluster 47 Tucanae, we detect two known millisecond pulsars (MSPs), and one previously reported MSP candidate, with hints of eclipsing behaviour.
The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) survey experiment that currently consists of three 0.42m small-aperture telescopes (SATs) and one 6m large-aperture telescope (LAT), located at an elevation of 5200m in the Atacama Desert in Chile. At the LAT's focal plane, SO will install >62,000 transition-edge sensor detectors across 13 optics tubes (OTs) within the Large Aperture Telescope Receiver (LATR), the largest cryogenic camera ever built to observe the CMB. Here we report on the validation of the LATR in the laboratory and the subsequent dark testing and validation within the LAT. We show that the LATR meets cryogenic, optical, and detector specifications required for high-sensitivity measurements of the CMB. At the time of writing, the LATR is installed in the LAT with six OTs (corresponding to >31,000 detectors), and the LAT mirrors and remaining seven OTs are undergoing development.
We present Mantis Shrimp, a multi-survey deep learning model for photometric redshift estimation that fuses ultra-violet (GALEX), optical (PanSTARRS), and infrared (UnWISE) imagery. Machine learning is now an established approach for photometric redshift estimation, with generally acknowledged higher performance in areas with a high density of spectroscopically identified galaxies over template-based methods. Multiple works have shown that image-based convolutional neural networks can outperform tabular-based color/magnitude models. In comparison to tabular models, image models have additional design complexities: it is largely unknown how to fuse inputs from different instruments which have different resolutions or noise properties. The Mantis Shrimp model estimates the conditional density estimate of redshift using cutout images. The density estimates are well calibrated and the point estimates perform well in the distribution of available spectroscopically confirmed galaxies with (bias = 1e-2), scatter (NMAD = 2.44e-2) and catastrophic outlier rate ($\eta$=17.53$\%$). We find that early fusion approaches (e.g., resampling and stacking images from different instruments) match the performance of late fusion approaches (e.g., concatenating latent space representations), so that the design choice ultimately is left to the user. Finally, we study how the models learn to use information across bands, finding evidence that our models successfully incorporates information from all surveys. The applicability of our model to the analysis of large populations of galaxies is limited by the speed of downloading cutouts from external servers; however, our model could be useful in smaller studies such as generating priors over redshift for stellar population synthesis.
The James Webb Space Telescope, launched in 2021, is an infrared observatory of novel design: deployable, with active optics, fully open to space for radiative cooling and orbiting the Lagrange point no. 2. This article explains the rationale leading to this specific design and describes the various other architectures that were considered along the way: from a monolithic 10-meter telescope in geosynchronous orbit to a 6-meter one in High Earth Orbit, then a 16-meter observatory on the Moon, a 4- or 6-meter one in an elliptical heliocentric orbit, and a segmented 8-meter one passively cooled to 50 K at L2, which was finally descoped to 6.6 meters. It also addresses the optimization for scientific performance, the challenge of dealing with such an ultra-low operating temperature, cost issues, supporting technology, modifications made during final design and, finally, how the architecture performs on orbit.
Recently, a double-lined binary (SB2) classical Cepheid, OGLE-LMC-CEP-1347, was discovered, with the orbital period (P$_{\rm orb} = 59$ days) five times shorter than of any binary Cepheid known before. The expected mass of the Cepheid was below $3.5$ M$_\odot$, which, if confirmed, would also probe the uncharted territory. The system configuration also pointed to the Cepheid as a merger. We present a novel method for determining precise physical parameters of binary Cepheids using both theory and observations. This q-PED method combines the measured mass ratio (q), pulsation (P), and evolutionary (E) models, and the known distance (D) supplemented with multi-band photometry. Applying it, we determined the mass of the Cepheid of $3.42 \pm 0.09$ M$_\odot$, its radius of $13.65 \pm 0.27$ R$_\odot$, the companion mass of $1.89 \pm 0.04$ M$_\odot$ and radius of $12.5 \pm 0.62$ R$_\odot$. With the current configuration, the apparent evolutionary age difference of almost $1$ Gyr between the components strongly favors the Cepheid merger origin scenario. If so, the actual age of the Cepheid would be $1.1$ Gyr, on the edge of Population II stars, indicating a significant fraction of Cepheids may be much older than typically assumed. We also applied our method to an eclipsing binary Cepheid with accurately determined physical parameters, obtaining a close agreement, which confirmed our method's reliability.
Einstein Telescope (ET) is a proposed next-generation Gravitational Wave (GW) interferometer designed to detect a large number of astrophysical and cosmological sources with unprecedented sensitivity. A key target for ET is the detection of a stochastic gravitational-wave background (SGWB), a faint signal from unresolved GW sources. In its proposed triangular configuration, correlated Newtonian noise of seismic origin poses some challenges for the SGWB detection. We study the impact of correlated noise on the SGWB detection and relative parameter estimation for ET in the triangular configuration, comparing it to a 2L configuration with two separated L-shaped detectors. We perform a Bayesian analysis on simulated data, which shows that accurate reconstruction of the SGWB parameters and instrumental noise is achievable if the noise is properly modeled. We illustrate that neglecting correlated noise leads to significant biases in the parameter reconstruction. Our results show that while the 2L configuration provides slightly better parameter estimation precision, mainly due to its longer arm length, the triangular configuration remains competitive when accurate noise modeling is provided.