More Like this

1. Stingray : A Modern Python Library for Spectral Timing

   https://doi.org/10.3847/1538-4357/ab258d  

   Huppenkothen, Daniela; Bachetti, Matteo; Stevens, Abigail L.; Migliari, Simone; Balm, Paul; Hammad, Omar; Khan, Usman Mahmood; Mishra, Himanshu; Rashid, Haroon; Sharma, Swapnil; et al (August 2019, The Astrophysical Journal)
2. p † q : a tool for prototyping many-body methods for quantum chemistry

   https://doi.org/10.1080/00268976.2021.1954709  

   Rubin, Nicholas C.; DePrince, A. Eugene (November 2021, Molecular Physics)
3. PyQMC : An all-Python real-space quantum Monte Carlo module in PySCF

   https://doi.org/10.1063/5.0139024  

   Wheeler, William A.; Pathak, Shivesh; Kleiner, Kevin G.; Yuan, Shunyue; Rodrigues, João N. B.; Lorsung, Cooper; Krongchon, Kittithat; Chang, Yueqing; Zhou, Yiqing; Busemeyer, Brian; et al (March 2023, The Journal of Chemical Physics)

   We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for a simple comparison between QMC calculations and other many-body wave function techniques, as well as access to high accuracy trial wave functions. 

   more » « less
4. Mahakala : A Python -based Modular Ray-tracing and Radiative Transfer Algorithm for Curved Spacetimes

   https://doi.org/10.3847/1538-4357/adc104  

   Sharma, Aniket; Medeiros, Lia; Wong, George N; Chan, Chi-kwan; Halevi, Goni; Mullen, Patrick D; Stone, James M (May 2025, The Astrophysical Journal)

   Abstract We introduceMahakala, aPython-based, modular, radiative ray-tracing code for curved spacetimes. We employ Google’sJAXframework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr spacetimes.JAXalso enablesMahakalato run in parallel on both CPUs and GPUs.Mahakalanatively uses the Cartesian Kerr–Schild coordinate system, which avoids numerical issues caused by the pole in spherical coordinate systems. We demonstrateMahakala’s capabilities by simulating 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature ofMahakalaallows us to quantitatively explore how different regions of the flow influence different image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole and peaks near the photon orbit. We also quantify the relative contribution of the disk, forward jet, and counterjet to 1.3 mm images. 

   more » « less
5. Radiation hydrodynamics modelling of kilonovae with SNEC

   https://doi.org/10.1093/mnras/stac399  

   Wu, Zhenyu; Ricigliano, Giacomo; Kashyap, Rahul; Perego, Albino; Radice, David (February 2022, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work. 

   more » « less