A self-gravity module for the PLUTO code

We present a novel implementation of an iterative solver for the solution of the Poisson equation in the PLUTO code for astrophysical fluid dynamics. Our solver relies on a relaxation method in which convergence is sought as the steady-state solution of a parabolic equation, whose time-discretization is governed by the \textit{Runge-Kutta-Legendre} (RKL) method. Our findings indicate that the RKL-based Poisson solver, which is both fully parallel and rapidly convergent, has the potential to serve as a practical alternative to conventional iterative solvers such as the \textit{Gauss-Seidel} (GS) and \textit{successive over-relaxation} (SOR) methods. Additionally, it can mitigate some of the drawbacks of these traditional techniques. We incorporate our algorithm into a multigrid solver to provide a simple and efficient gravity solver that can be used to obtain the gravitational potentials in self-gravitational hydrodynamics. We test our implementation against a broad range of standard self-gravitating astrophysical problems designed to examine different aspects of the code. We demonstrate that the results match excellently with the analytical predictions (when available), and the findings of similar previous studies.Comment: Submitted to ApJS. Comments are welcom

A CRITICAL ANALYSIS ON THE BASIC CONCEPT OF PRAMEHA IN CONVENTIONAL PARLANCE

Systemic information about several diseases has been documented in Ayurvedic classical texts. Although conventional entity of all of those diseases are not well established. Understanding of such diseases in conventional parlance is essential for an evidence based approach of Ayurveda. Prameha is one of such disease that is most widely described in almost all classical Ayurvedic texts but not well established in conventional parlance. The disease Premeha, has been named on its major clinical signs Avila-Prabhuta-Mutra (Excess and contaminated urine). In ancient text compiled by Acharya Charaka, Acharya Sushruta, Acharya Vagbhatta and many others, we get detailed description about this disease. Meda Dusti is considered as a key pathological phenomenon behind the development of Prameha. A conventional entity of this disease is still now doubtful. This review aims at scanning of both Ayurvedic and conventional medical literatures as well as published research articles to explore the basic concept of Prameha in conventional parlance

Kinetic control of competing nuclei in a dimer lattice-gas model

Nucleation is a key step in the synthesis of new material from solution. Well-established lattice-gas models can be used to gain insight into the basic physics of nucleation pathways involving a single nucleus type. In many situations a solution is supersaturated with respect to more than one precipitating phase. This can generate a population of both stable and metastable nuclei on similar timescales and hence complex nucleation pathways involving competition between the two. In this study we introduce a lattice-gas model based on two types of interacting dimer representing particles in solution. Each type of dimer nucleates to a specific space-filling structure. Our model is tuned such that stable and metastable phases nucleate on a similar timescale. Either structure may nucleate first, with probability sensitive to the relative rate at which solute is replenished from their respective reservoirs. We calculate these nucleation rates via Forward-Flux Sampling and demonstrate how the resulting data can be used to infer the nucleation outcome and pathway. Possibilities include direct nucleation of the stable phase, domination of long-lived metastable crystallites, and pathways in which the stable phase nucleates only after multiple post-critical nuclei of the metastable phase have appeared

Hard core lattice gas with third next-nearest neighbor exclusion on triangular lattice : one or two phase transitions?

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density, with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Though a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show unambiguously that there is only a single phase transition in the system between a low-density fluid phase and a high-density ordered sublattice phase, and that a hexatic phase is absent. The transition is shown to be first order in nature and the critical parameters are determined accurately

Rejection-free cluster Wang-Landau algorithm for hard-core lattice gases

We introduce a rejection-free, flat histogram, cluster algorithm to determine the density of states of hard-core lattice gases. We show that the algorithm is able to efficiently sample low entropy states that are usually difficult to access, even when the excluded volume per particle is large. The algorithm is based on simultaneously evaporating all the particles in a strip and reoccupying these sites with a new appropriately chosen configuration. We implement the algorithm for the particular case of the hard-core lattice gas in which the first k next-nearest neighbors of a particle are excluded from being occupied. It is shown that the algorithm is able to reproduce the known results for k=1,2,3 both on the square and cubic lattices. We also show that, in comparison, the corresponding flat histogram algorithms with either local moves or unbiased cluster moves are less accurate and do not converge as the system size increases

Modelling observable signatures of jet-ISM interaction: thermal emission and gas kinematics

Relativistic jets are believed to have a substantial impact on the gas dynamics and evolution of the interstellar medium (ISM) of their host galaxies. In this paper, we aim to draw a link between the simulations and the observable signatures of jet-ISM interactions by analyzing the emission morphology and gas kinematics resulting from jet-induced shocks in simulated disc and spherical systems. We find that the jet-induced laterally expanding forward shock of the energy bubble sweeping through the ISM causes large-scale outflows, creating shocked emission and high-velocity dispersion in the entire nuclear regions ($\sim2$ kpcs) of their hosts. The jetted systems exhibit larger velocity widths (> 800 km/s), broader Position-Velocity maps and distorted symmetry in the disc's projected velocities than systems without a jet. We also investigate the above quantities at different inclination angles of the observer with respect to the galaxy. Jets inclined to the gas disc of its host are found to be confined for longer times, and consequently couple more strongly with the disc gas. This results in prominent shocked emission and high-velocity widths, not only along the jet's path, but also in the regions perpendicular to them. Strong interaction of the jet with a gas disc can also distort its morphology. However, after the jets escape their initial confinement, the jet-disc coupling is weakened, thereby lowering the shocked emission and velocity widths.Comment: Matches the Published versio

Data for Nucleation rate in the two dimensional Ising model in the presence of random impurities

Nucleation phenomena are ubiquitous in nature and the presence of impurities in every real and experimental system is unavoidable. Yet numerical studies of nucleation are nearly always conducted for entirely pure systems. We have studied the behaviour of the droplet free energy in two dimensional Ising model in the presence of randomly positioned static and dynamic impurities. We have shown that both the free energy barrier height and critical nucleus size monotonically decreases with increasing the impurity density for the static case. We have compared the nucleation rates obtained from the Classical Nucleation Theory and the Forward Flux Sampling method for different densities of the static impurities. The results show good agreement. In the case of dynamic impurities, we observe preferential occupancy of the impurities at the boundary positions of the nucleus when the temperature is low. This further boosts enhancement of the nucleation rate due to lowering of the effective interfacial free energy

Data for Kinetic control of competing nuclei in a dimer lattice-gas model

Nucleation is a key step in the synthesis of new material from solution. Well-established lattice-gas models can be used to gain insight into the basic physics of nucleation pathways involving a single nucleus type. In many situations a solution is supersaturated with respect to more than one precipitating phase. This can generate a population both stable and metastable nuclei on similar timescales and hence complex nucleation pathways involving competition between the two. In this study we introduce a lattice-gas model based on two types of interacting dimer representing particles in solution. Each type of dimer nucleates to a specific space-filling structure. Our model is tuned such that stable and metastable phases nucleate on a similar timescale. Either structure may nucleate first, with probability sensitive to dimer mobility. We calculate these nucleation rates via Forward-Flux Sampling and demonstrate how the resulting data can be used to infer the nucleation outcome and pathway. Possibilities include direct nucleation of the stable phase, domination of long-lived metastable crystallites, and pathways in which the stable phase nucleates only after multiple post-critical nuclei of the metastable phase have appeared