Astro arXiv | astro-ph.IM

H-FISTA: A hierarchical algorithm for phase retrieval with application to pulsar dynamic spectra by Stefan Osłowski et al. on Tuesday 22 November

A pulsar dynamic spectrum is an inline digital hologram of the interstellar

medium; it encodes information on the propagation paths by which signals have

travelled from source to telescope. To decode the hologram it is necessary to

"retrieve" the phases of the wavefield from intensity measurements, which

directly gauge only the field modulus, by imposing additional constraints on

the model. We present a new method for phase retrieval in the context of pulsar

spectroscopy. Our method makes use of the Fast Iterative Shrinkage Thresholding

Algorithm (FISTA) to obtain sparse models of the wavefield in a hierarchical

approach with progressively increasing depth. Once the tail of the noise

distribution is reached the hierarchy terminates with a final, unregularised

optimisation. The result is a fully dense model of the complex wavefield that

permits the discovery of faint signals by appropriate averaging. We illustrate

the performance of our method on synthetic test cases and on real data. Our

algorithm, which we call H-FISTA, is implemented in the Python programming

language and is freely available.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11129v1

The Science Performance of JWST as Characterized in Commissioning by Jane Rigby et al. on Tuesday 22 November

This paper characterizes the actual science performance of the James Webb

Space Telescope (JWST), as determined from the six month commissioning period.

We summarize the performance of the spacecraft, telescope, science instruments,

and ground system, with an emphasis on differences from pre-launch

expectations. Commissioning has made clear that JWST is fully capable of

achieving the discoveries for which it was built. Moreover, almost across the

board, the science performance of JWST is better than expected; in most cases,

JWST will go deeper faster than expected. The telescope and instrument suite

have demonstrated the sensitivity, stability, image quality, and spectral range

that are necessary to transform our understanding of the cosmos through

observations spanning from near-earth asteroids to the most distant galaxies.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.05632v3

Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses by Aymeric Galan et al. on Tuesday 22 November

Modeling the mass distribution of galaxy-scale strong gravitational lenses is

a task of increasing difficulty. The high-resolution and depth of imaging data

now available render simple analytical forms ineffective at capturing lens

structures spanning a large range in spatial scale, mass scale, and morphology.

In this work, we address the problem with a novel multiscale method based on

wavelets. We tested our method on simulated Hubble Space Telescope (HST)

imaging data of strong lenses containing the following different types of mass

substructures making them deviate from smooth models: (1) a localized small

dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a

nonlocalized population of subhalos along the line of sight, and (3)

galaxy-scale multipoles that break elliptical symmetry. We show that wavelets

are able to recover all of these structures accurately. This is made

technically possible by using gradient-informed optimization based on automatic

differentiation over thousands of parameters, which also allow us to sample the

posterior distributions of all model parameters simultaneously. By

construction, our method merges the two main modeling paradigms - analytical

and pixelated - with machine-learning optimization techniques into a single

modular framework. It is also well-suited for the fast modeling of large

samples of lenses. All methods presented here are publicly available in our new

Herculens package.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.05763v2

A tactile model of the night summer northern sky for the teaching of astronomy to the BVI by E. Pérez-Montero et al. on Tuesday 22 November

Astroaccesible is an outreach project hosted by the Instituto de

Astrof\'{\i}sica de Andaluc\'{\i}a - CSIC aimed at the teaching and

popularisation of the astronomy among all publics independently of their

capabilities and abilities, paying special attention to the collective of blind

and visually impaired (BVI). Among the different strategies and resources using

in our project, we have developed new 3D models representing in relief some of

the stars, constellations and deep sky objects that can be observed during

night from the Northern hemisphere in spring and summer. These models can be

used by BVI to transmit to them the spatial configuration of the sky during

night, but can be also used as an additional resource for all kind of publics

to complement their sensorial experience. We also describe additional resources

based on sounds that can also be employed to get deeper into this

multisensorial experience. Finally, we summarize some of the activities and the

context in which this new material has been used in the last 2 years.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11707v1

3D Detection and Characterisation of ALMA Sources through Deep Learning by Michele Delli Veneri et al. on Tuesday 22 November

We present a Deep-Learning (DL) pipeline developed for the detection and

characterization of astronomical sources within simulated Atacama Large

Millimeter/submillimeter Array (ALMA) data cubes. The pipeline is composed of

six DL models: a Convolutional Autoencoder for source detection within the

spatial domain of the integrated data cubes, a Recurrent Neural Network (RNN)

for denoising and peak detection within the frequency domain, and four Residual

Neural Networks (ResNets) for source characterization. The combination of

spatial and frequency information improves completeness while decreasing

spurious signal detection. To train and test the pipeline, we developed a

simulation algorithm able to generate realistic ALMA observations, i.e. both

sky model and dirty cubes. The algorithm simulates always a central source

surrounded by fainter ones scattered within the cube. Some sources were

spatially superimposed in order to test the pipeline deblending capabilities.

The detection performances of the pipeline were compared to those of other

methods and significant improvements in performances were achieved. Source

morphologies are detected with subpixel accuracies obtaining mean residual

errors of $10^{-3}$ pixel ($0.1$ mas) and $10^{-1}$ mJy/beam on positions and

flux estimations, respectively. Projection angles and flux densities are also

recovered within $10\%$ of the true values for $80\%$ and $73\%$ of all sources

in the test set, respectively. While our pipeline is fine-tuned for ALMA data,

the technique is applicable to other interferometric observatories, as SKA,

LOFAR, VLBI, and VLTI.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11462v1

Daksha: On Alert for High Energy Transients by Varun Bhalerao et al. on Tuesday 22 November

We present Daksha, a proposed high energy transients mission for the study of

electromagnetic counterparts of gravitational wave sources, and gamma ray

bursts. Daksha will comprise of two satellites in low earth equatorial orbits,

on opposite sides of earth. Each satellite will carry three types of detectors

to cover the entire sky in an energy range from 1 keV to >1 MeV. Any transients

detected on-board will be announced publicly within minutes of discovery. All

photon data will be downloaded in ground station passes to obtain source

positions, spectra, and light curves. In addition, Daksha will address a wide

range of science cases including monitoring X-ray pulsars, studies of

magnetars, solar flares, searches for fast radio burst counterparts, routine

monitoring of bright persistent high energy sources, terrestrial gamma-ray

flashes, and probing primordial black hole abundances through lensing. In this

paper, we discuss the technical capabilities of Daksha, while the detailed

science case is discussed in a separate paper.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.12055v1

Line Emission Mapper LEM : Probing the physics of cosmic ecosystems by Ralph Kraft et al. on Tuesday 22 November

The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will

answer the outstanding questions of the Universe's structure formation. It will

also provide transformative new observing capabilities for every area of

astrophysics, and to heliophysics and planetary physics as well. LEM's main

goal is a comprehensive look at the physics of galaxy formation, including

stellar and black-hole feedback and flows of baryonic matter into and out of

galaxies. These processes are best studied in X-rays, and emission-line mapping

is the pressing need in this area. LEM will use a large microcalorimeter

array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft

X-ray line emission from objects that constitute galactic ecosystems. These

include supernova remnants, star-forming regions, superbubbles, galactic

outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs

in other galaxies), the Circumgalactic Medium in the Milky Way and other

galaxies, and the Intergalactic Medium at the outskirts and beyond the confines

of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV

band will make it possible to disentangle the faintest emission lines in those

objects from the bright Milky Way foreground, providing groundbreaking

measurements of the physics of these plasmas, from temperatures, densities,

chemical composition to gas dynamics. While LEM's main focus is on galaxy

formation, it will provide transformative capability for all classes of

astrophysical objects, from the Earth's magnetosphere, planets and comets to

the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling

gas in galaxy clusters. In addition to pointed observations, LEM will perform a

shallow all-sky survey that will dramatically expand the discovery space.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.09827v2

Low-Energy Electron-Track Imaging for a Liquid Argon Time-Projection-Chamber Telescope Concept using Probabilistic Deep Learning by M. Buuck et al. on Monday 21 November

The GammaTPC is an MeV-scale single-phase liquid argon

time-projection-chamber gamma-ray telescope concept with a novel dual-scale

pixel-based charge-readout system. It promises to enable a significant

improvement in sensitivity to MeV-scale gamma-rays over previous telescopes.

The novel pixel-based charge readout allows for imaging of the tracks of

electrons scattered by Compton interactions of incident gamma-rays. The two

primary contributors to the accuracy of a Compton telescope in reconstructing

an incident gamma-ray's original direction are its energy and position

resolution. In this work, we focus on using deep learning to optimize the

reconstruction of the initial position and direction of electrons scattered in

Compton interactions, including using probabilistic models to estimate

predictive uncertainty. We show that the deep learning models are able to

predict locations of Compton scatters of MeV-scale gamma-rays from simulated

pixel-based data to better than 0.6 mm RMS error, and are sensitive to the

initial direction of the scattered electron. We compare and contrast different

deep learning uncertainty estimation algorithms for reconstruction

applications. Additionally, we show that event-by-event estimates of the

uncertainty of the locations of the Compton scatters can be used to select

those events that were reconstructed most accurately, leading to improvement in

locating the origin of gamma-ray sources on the sky.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.07805v2

GRaM-X: A new GPU-accelerated dynamical spacetime GRMHD code for Exascale computing with the Einstein Toolkit by Swapnil Shankar et al. on Monday 21 November

We present GRaM-X (General Relativistic accelerated Magnetohydrodynamics on

AMReX), a new GPU-accelerated dynamical-spacetime general relativistic

magnetohydrodynamics (GRMHD) code which extends the GRMHD capability of

Einstein Toolkit to GPU-based exascale systems. GRaM-X supports 3D adaptive

mesh refinement (AMR) on GPUs via a new AMR driver for the Einstein Toolkit

called CarpetX which in turn leverages AMReX, an AMR library developed for use

by the United States DOE's Exascale Computing Project (ECP). We use the Z4c

formalism to evolve the equations of GR and the Valencia formulation to evolve

the equations of GRMHD. GRaM-X supports both analytic as well as tabulated

equations of state. We implement TVD and WENO reconstruction methods as well as

the HLLE Riemann solver. We test the accuracy of the code using a range of

tests on static spacetime, e.g. 1D MHD shocktubes, the 2D magnetic rotor and a

cylindrical explosion, as well as on dynamical spacetimes, i.e. the

oscillations of a 3D TOV star. We find excellent agreement with analytic

results and results of other codes reported in literature. We also perform

scaling tests and find that GRaM-X shows a weak scaling efficiency of $\sim

40-50\%$ on 2304 nodes (13824 NVIDIA V100 GPUs) with respect to single-node

performance on OLCF's supercomputer Summit.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2210.17509v2

StarDICE I: sensor calibration bench and absolute photometric calibration of a Sony IMX411 sensor by Marc Betoule et al. on Monday 21 November

The Hubble diagram of type-Ia supernovae (SNe-Ia) provides cosmological

constraints on the nature of dark energy with an accuracy limited by the flux

calibration of currently available spectrophotometric standards. The StarDICE

experiment aims at establishing a 5-stage metrology chain from NIST photodiodes

to stars, with a targeted accuracy of \SI{1}{mmag} in $griz$ colors. We present

the first two stages, resulting in the calibration transfer from NIST

photodiodes to a demonstration \SI{150}{Mpixel} CMOS sensor (Sony IMX411ALR as

implemented in the QHY411M camera by QHYCCD). As a side-product, we provide

full characterization of this camera. A fully automated spectrophotometric

bench is built to perform the calibration transfer. The sensor readout

electronics is studied using thousands of flat-field images from which we

derive stability, high resolution photon transfer curves and estimates of the

individual pixel gain. The sensor quantum efficiency is then measured relative

to a NIST-calibrated photodiode. Flat-field scans at 16 different wavelengths

are used to build maps of the sensor response. We demonstrate statistical

uncertainty on quantum efficiency below \SI{0.001}{e^-/\gamma} between

\SI{387}{nm} and \SI{950}{nm}. Systematic uncertainties in the bench optics are

controlled at the level of \SI{1e-3}{e^-/\gamma}. Uncertainty in the overall

normalization of the QE curve is 1\%. Regarding the camera we demonstrate

stability in steady state conditions at the level of \SI{32.5}{ppm}.

Homogeneity in the response is below \SI{1}{\percent} RMS across the entire

sensor area. Quantum efficiency stays above \SI{50}{\percent} in most of the

visible range, peaking well above \SI{80}{\percent} between \SI{440}{nm} and

\SI{570}{nm}. Differential non-linearities at the level of \SI{1}{\percent} are

detected. A simple 2-parameter model is proposed to mitigate the effect.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.04913v2