Related papers: Mahakala: a Python-based Modular Ray-tracing and R…
We present RAYTHEIA, a high-performance reverse ray-tracing algorithm designed to efficiently solve three-dimensional direction-dependent equations in astronomical simulations. The algorithm uses a dual-grid framework in which the native…
We present a PDE-based approach for the multidimensional extrapolation of smooth scalar quantities across interfaces with kinks and regions of high curvature. Unlike the commonly used method of [2] in which normal derivatives are…
This paper presents OCTOPUS, a relativistic ray-tracing algorithm developed within a Fortran-based, OpenMP-accelerated framework, designed for asymptotically flat, spherically symmetric curved spacetimes. The code efficiently and accurately…
A new Markov Chain Monte Carlo method for simulating the dynamics of molecular systems characterized by hard-core interactions is introduced. In contrast to traditional Kinetic Monte Carlo approaches, where the state of the system is…
We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code CIGALE (Code Investigating GALaxy emission). The original CIGALE builds a grid of…
X-ray speckles have been used for a wide variety of experiments, ranging from imaging (and tomography), wavefront sensing, spatial coherence measurements all the way to x-ray photon correlation spectroscopy (XPCS) and ptychography. In the…
We present a new code, named COCAL - Compact Object CALculator, for the computation of equilibriums and quasiequilibrium initial data sets of single or binary compact objects of all kinds. In the cocal code, those solutions are calculated…
Aims. We present an implementation of an algorithm for 3D time-dependent Monte Carlo radiative transfer. It allows one to simulate temperature distributions as well as images and spectral energy distributions of the scattered light and…
We propose that high frequency quasi-periodic oscillations (HFQPOs) can be produced from randomly-formed X-ray bursts (flashes) by plasma interior to the ergosphere of a rapidly-rotating black hole. We show by direct computation of their…
We present GANDALF, a JAX-based spectral solver for Kinetic Reduced MHD (KRMHD) turbulence designed to lower infrastructure barriers to plasma turbulence research. Existing production codes require specialized HPC infrastructure and…
The differential evolution indicators, recently introduced for imaging mechanical evolution in highly scattering solids, is examined in a laboratory setting with the focus on spatiotemporal tracking of an advancing damage zone in an elastic…
We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We…
Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with…
The JILA Multidimensional Optical Nonlinear SpecTRometer (JILA-MONSTR) is a robust, ultra-stable platform consisting nested and folded Michelson interferometers that can be actively phase stabilized. This platform generates a square of…
Scattering can rapidly degrade our ability to form an optical image, to the point where only speckle-like patterns can be measured. Truly non-invasive imaging through a strongly scattering obstacle is difficult, and usually reliant on a…
We present a numerical method for studying the normal modes of accretion flows around black holes. In this first paper, we focus on two-dimensional, viscous, hydrodynamic disks, for which the linear modes have been calculated analytically…
Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the…
We introduce microJAX, the first fully differentiable implementation of the image-centered ray-shooting (ICRS) algorithm for gravitational microlensing. Built on JAX and its XLA just-in-time compiler, microJAX exploits GPU parallelism while…
A new high time resolution observing mode for the Murchison Widefield Array (MWA) is described, enabling full polarimetric observations with up to 30.72 MHz of bandwidth and a time resolution of ~0.8 $\mu$s. This mode makes use of a…
We present $\texttt{PyBird-JAX}$, a differentiable, $\texttt{JAX}$-based implementation of $\texttt{PyBird}$, using internal neural network emulators to accelerate computationally costly operations for rapid large-scale structure (LSS)…