Abstract. Modified gravity theories can introduce modifications to the Poisson equation in the Newtonian limit. As a result, we expect to see interesting features of these modifications inside stellar objects. White dwarf stars are one of the most well studied stars in stellar astrophysics. We explore the effect of modified gravity theories inside white dwarfs. We derive the modified stellar structure equations and solve them to study the mass-radius relationships for various modified gravity theories. We also constrain the parameter space of these theories from observations.

A subset of Type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies (∼1052 erg) that cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) that suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae form a rapidly rotating and strongly magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work, we use the jetted outflow produced in a 3D general relativistic magnetohydrodynamic CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of 6.231 × 1051−1.725 × 1052 erg, the engine time-scale in the range of 0.479−1.159 s and the engine half-opening angle in the range of ∼9°−19°. Using these as central engines, we perform 2D special relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl that demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles (∼11°) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl’s without GRBs may be associated with wider outflows.

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 ∼ 40-50% on 2304 nodes (13824 NVIDIA V100 GPUs) with respect to single-node performance on OLCF's supercomputer Summit.

We present a new method to determine the probability distribution of the 3-D shapes of galaxy clusters from the 2-D images using stereology. In contrast to the conventional approach of combining different data sets (such as X-rays, Sunyaev-Zeldovich effect and lensing) to fit a 3-D model of a galaxy cluster for each cluster, our method requires only a single data set, such as X-ray observations or Sunyaev-Zeldovich effect observations, consisting of sufficiently large number of clusters. Instead of reconstructing the 3-D shape of an individual object, we recover the probability distribution function (PDF) of the 3-D shapes of the observed galaxy clusters. The shape PDF is the relevant statistical quantity which can be compared with the theory and used to test the cosmological models. We apply this method to publicly available Chandra X-ray data of 89 well resolved galaxy clusters. Assuming ellipsoidal shapes, we find that our sample of galaxy clusters is a mixture of prolate and oblate shapes, with a preference for oblateness with the most probable ratio of principle axes 1.4 : 1.3 : 1. The ellipsoidal assumption is not essential to our approach and our method is directly applicable to non-ellipsoidal shapes. Our method is insensitive to the radial density and temperature profiles of the cluster. Our method is sensitive to the changes in shape of the X-ray emitting gas from inner to outer regions and we find evidence for variation in the 3-D shape of the X-ray emitting gas with distance from the centre.