Related papers: Mahakala: a Python-based Modular Ray-tracing and R…
X-ray propagation-based imaging techniques are well-established at synchrotron radiation and laboratory sources. However, most reconstruction algorithms for such image modalities, also known as phase retrieval algorithms, have been…
A new CAD program, PolyCAD, designed for X-ray photon tracing in polycapillary optics is described. To understand the PolyCAD code and its results, the theoretical basis of X-ray transmission by a single cylindrical channel (monocapillary)…
Trajectory data represent a trace of an object that changes its position in space over time. This kind of data is complex to handle and analyze, since it is generally produced in huge quantities, often prone to errors generated by the…
The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on accuracy and precision of numerical models of plasma…
We present {\tt radpol} - a numerical scheme for integrating multifrequency polarized radiative transfer equations along rays propagating in a curved spacetime. The scheme includes radiative processes such as synchrotron emission,…
We describe Rabacus, a Python package for calculating the transfer of hydrogen ionizing radiation in simplified geometries relevant to astronomy and cosmology. We present example solutions for three specific cases: 1) a semi-infinite slab…
Goupil is a software library designed for the Monte Carlo transport of low-energy gamma-rays, such as those emitted from radioactive isotopes. The library is distributed as a Python module. It implements a dedicated backward sampling…
When modeling astronomical objects throughout the universe, it is important to correctly treat the limitations of the data, for instance finite resolution and sensitivity. In order to simulate these effects, and to make radiative transfer…
Synchrotron processes, the radiative processes associated with the interaction of energetic charged particles with magnetic field, are of interest in many areas in astronomy, from the interstellar medium to extreme environments near compact…
We present MCRaT, a Monte Carlo Radiation Transfer code for self-consistently computing the light curves and spectra of the photospheric emission from relativistic, unmagnetized jets. We apply MCRaT to a relativistic hydrodynamic simulation…
A new numerical method is introduced to study the problem of time evolution of generic non-linear dynamical systems in four-dimensional spacetimes. It is assumed that the time level surfaces are foliated by a one-parameter family of…
This paper presents \textsc{hankel}, a pure-python code for solving Hankel-type integrals and transforms. Such transforms are common in the physical sciences, especially appearing as the radial solution to angularly symmetric Fourier…
The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable…
Eikonal, or ray tracing, methods are commonly used to estimate the propagation of radio frequency fields in plasmas. While the information gained from the rays is quite useful, an approximate solution for the fields would also be valuable,…
Retrieving images transmitted through multi-mode fibers is of growing interest, thanks to their ability to confine and transport light efficiently in a compact system. Here, we demonstrate machine-learning-based decoding of large-scale…
We present \code{phaser}, an open-source Python package that provides a unified interface to both conventional and gradient descent-based ptychographic algorithms. Features such as mixed-state probe, probe position correction, and…
High-$p_T$ jets are an important tool for characterizing the quark-gluon plasma (QGP) created in heavy-ion collisions. However, a precise understanding of the jet-medium interaction is still lacking, and the development of more…
In the context of track fitting problems by a Kalman filter, the appropriate functional forms of the elements of the random process noise matrix are derived for tracking through thick layers of dense materials and magnetic field. This work…
A concise and highly performant branchless formulation of a Joseph-type interpolating ray-casting algorithm for the computation of X-ray projections is presented. It efficiently utilizes the hardware resources of modern graphics processing…
Infinitesimal electric and magnetic dipoles are widely used as an equivalent radiating source model. In this paper, an improved method for dipole extraction from magnitude-only electromagnetic-field data based on genetic algorithm and…