We introduce pinta, a pipeline for reducing the uGMRT raw pulsar timing data developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. Furthermore, the results of a calibration experiment carried out to determine the relative time offsets between different back-end modes and the correct interpretation of the observatory frequency settings at the uGMRT, which are crucial for performing precision pulsar timing, are also discussed.

The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.

We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to make possible large-scale simulations of magnetized plasmas in strong, dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate that GiRaFFE passes a large suite of both flat and curved-spacetime code tests passed by a number of other state-of-the-art GRFFE codes, and is thus ready for production-scale simulations of GRFFE phenomena of key interest to relativistic astrophysics.

The young star V1331 Cyg received previous attention because it is surrounded by an optical, arc-like reflection nebula. V1331 Cyg is commonly considered to be a candidate for an object that has undergone an FU-Ori (FUOR) the outbreak in the past. This in turn could lead to a time-varying appearance of the dusty arcs that may be revealed by multi-epoch imaging. In particular, a radial colour analysis of the dust arcs can then be attempted to check whether radial grain size distribution was modified by a previous FUOR wind. Second-epoch imaging of V1331 Cyg was obtained by us in 2009 using the Hubble Space Telescope (HST). By comparing this to archival HST data from 2000, we studied the time evolution of the circumstellar nebulae. After a point spread function subtraction using model point spread functions, we used customised routines to perform a proper motion analysis. The nebula expansion was first derived by deconvolving and correlating the two-epoch radial brightness profiles. Additional data from other facilities TLS, UKIDSS, SPITZER, and HERSCHEL were also incorporated to improve our understanding of the star in terms of environment, viewing angle, bipolar outflow length, and the FUOR phenomenon. The derived radial colour profiles do not indicate a spatial separation of the dust grain sizes. The HERSCHEL 160 micron images show for the time thermal emission from dust probably residing in the outer arc. By viewing V1331 Cyg almost pole-on, the length of the bipolar outflow exceeds previous estimates by far. The outer arc expansion timescale is consistent with the implantation time of the CO torus, which supports the hypothesis of an outburst that occurred a few thousand years ago. The azimuthal colour variation of the outer arc is probably due to changes of the scattering angle, imposed by a tilt or helical geometry of the dust configuration.

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.

For axially symmetric flow of dissipationless inhomogeneous fluid onto a non rotating astrophysical black hole under the influence of a generalized pseudo-Schwarzschild gravitational potential, we investigate the influence of the background flow configuration on determining the salient features of the corresponding acoustic geometry. The acoustic horizon for the aforementioned flow structure has been located and the corresponding acoustic surface gravity $\kappa$ as well as the associated analogue Hawking temperature $T_{\rm AH}$ has been calculated {\it analytically}. The dependence of $\kappa$ on the flow geometry as well as on the nature of the back ground black hole space time (manifested through the nature of the pseudo-Schwarzschild potential used) has been discussed. Dependence of the value of $\kappa$ on various initial boundary conditions governing the dynamic and the thermodynamic properties of the background fluid flow has also been studied.