We present the first results of a project called SAGAN, dedicated solely to the studies of relatively rare megaparsec-scale radio galaxies in the Universe, called the Giant Radio Galaxies (GRGs). We have identified 162 new GRGs primarily from the NVSS with sizes ranging from ~0.71 Mpc to ~2.82 Mpc in the redshift range of ~0.03 - 0.95. As part of the project SAGAN, we have created a database of all known GRGs-`GRG-catalogue' from literature (including our new sample) of 820 sources. For the first time, we present the multi-wavelength properties of the largest sample of GRGs, providing new insights about their nature. Our results firmly establish that the distributions of radio spectral index and the black hole mass of GRGs do not differ from the corresponding distributions of normal sized radio galaxies (RGs). However, GRGs have lower Eddington ratio than RGs. Using the mid-infrared data, we have classified GRGs in terms of their accretion mode: either high-power radiatively-efficient, high-excitation state or a radiatively-inefficient low-excitation state. This enables us to compare key physical properties of their AGN like the black hole mass, spin, Eddington ratio, jet kinetic power, total radio power, magnetic field and size. We find that GRGs in high excitation state statistically have larger sizes, radio power, jet kinetic power and Eddington ratio than those in low excitation state. Our analysis reveals a strong correlation between black hole's accretion efficiency and jet kinetic power, thus suggesting a disk-jet coupling. Our environmental study reveals that ~10% of all GRGs may reside at the centres of galaxy clusters, in a denser galactic environment while majority seem to reside in sparse environment. The probability of finding the brightest cluster galaxy (BCG) as GRG is quite low and even lower for high mass clusters. Abridged.

The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust particles which overcomes the stiffness introduced by the linear coupling term. Besides being time-reversible and globally second-order accurate in time, the exponential integrator provides energy errors which are always bounded and it remains stable in the limit of arbitrarily small particle stopping times yielding the correct asymptotic solution. Such properties make this method preferable to the more widely used semi-implicit or fully implicit schemes at a very modest increase in computational cost. Coupling between particles and grid quantities is achieved through particle deposition and field-weighting techniques borrowed from Particle-In-Cell simulation methods. In this respect, we derive new weight factors in curvilinear coordinates that are more accurate than traditional volume- or area-weighting. A comprehensive suite of numerical benchmarks is presented to assess the accuracy and robustness of the algorithm in Cartesian, cylindrical and spherical coordinates. Particular attention is devoted to the streaming instability which is analyzed in both local and global disk models. The module is part of the PLUTO code for astrophysical gas-dynamics and it is mainly intended for the numerical modeling of protoplanetary disks in which solid and gas interact via aerodynamic drag.

In this paper, we explore the behavior of optical polarization during the multi-wavelength outburst of the blazar Mrk 421 detected in February 2010. We use optical polarization measurements in the wavelength range 500-700 nm from SPOL observations available between January 1, 2010 and March 31, 2010 (MJD 55197-55286) including the period of multi-wavelength flaring activity detected from the source around February 16-17, 2010 (MJD 55243-55244). We also use near simultaneous optical and radio flux measurements from SPOL in V and R bands and OVRO at 15 GHz respectively. We find that the emissions in the optical and radio bands do not show any significant change in the source activity unlike at X-ray and $\gamma$--ray energies during the outburst. The optical and radio flux measurements are found to be consistent with the long term quiescent state emission of the source. Moreover, the linear polarization in the wavelength range 500-700 nm decreases to a minimum value of 1.6$\%$ during the X-ray and $\gamma$--ray outburst which is significantly lower than the long term average value of $\sim$ 4.2$\%$. The angle of polarization varies between 114$^\circ$-163$^\circ$ with a preferred average value of $\sim$ 137$^\circ$ during this period. We estimate the degree of polarization intrinsic to the jet taking into account the host galaxy contamination in R band and compare this with the theoretical synchrotron polarization estimated for a power law distribution of relativistic electrons gyrating in an emission region filled with ordered and chaotic magnetic fields. The intrinsic linear polarization estimated for different epochs during the above period is found to be consistent with the theoretical synchrotron polarization produced by the relativistic electrons with power law spectral index $\sim$ 2.2.

The outflowing magnetized wind from a host star shapes planetary and exoplanetary magnetospheres dictating the extent of its impact. We carry out three-dimensional (3D) compressible magnetohydrodynamic (MHD) simulations of the interactions between magnetized stellar winds and planetary magnetospheres corresponding to a far-out star-planet system, with and without planetary dipole obliquity. We identify the pathways that lead to the formation of a dynamical steady-state magnetosphere and find that magnetic reconnection plays a fundamental role in the process. The magnetic energy density is found to be greater on the night-side than that on the day-side and the magnetotail is comparatively more dynamic. Magnetotail reconnection events are seen to associated with stellar wind plasma injection into the inner magnetosphere. We further study magnetospheres with extreme tilt angles keeping in perspective the examples of Uranus and Neptune. High dipole obliquities may also manifest due to polarity excursions during planetary field reversals. We find that global magnetospheric reconnection sites change for large planetary dipole obliquity and more complex current sheet structures are generated. We discuss the implications of these findings for injection of interplanetary species and energetic particles into the inner magnetosphere, auroral activity and magnetospheric radio emission. This study is relevant for exploring star planet interactions in the solar and extra-solar systems.

We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through an Eulerian grid where the (relativistic) MHD equations are solved concurrently. Lagrangian particles follow fluid streamlines and represent ensembles of (real) relativistic particles with a finite energy distribution. The spectral distribution of each particle is updated in time by solving the relativistic cosmic ray transport equation based on local fluid conditions. This enables us to account for a number of physical processes, such as adiabatic expansion, synchrotron and inverse Compton emission. An accurate semi-analytically numerical scheme that combines the method of characteristics with a Lagrangian discretization in the energy coordinate is described. In presence of (relativistic) magnetized shocks, a novel approach to consistently model particle energization due to diffusive shock acceleration has been presented. Our approach relies on a refined shock-detection algorithm and updates the particle energy distribution based on the shock compression ratio, magnetic field orientation and amount of (parameterized) turbulence. The evolved distribution from each \textit{Lagrangian} particle is further used to produce observational signatures like emission maps and polarization signals accounting for proper relativistic corrections. We further demonstrate the validity of this hybrid framework using standard numerical benchmarks and evaluate the applicability of such a tool to study high energy emission from extra-galactic jets.

We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged particles, or cosmic rays (CR). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques. The module can be used to describe either test particles motion in the fluid electromagnetic field or to solve the fully coupled MHD-PIC system of equations with particle backreaction on the fluid as originally introduced by \cite{Bai_etal.2015}. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm's law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided the Larmor gyration scale is properly resolved. When applicable, this formulation avoids to resolve microscopic scales offering a substantial computational saving with respect to PIC simulations. We present a fully-conservative formulation which is second-order accurate in time and space and extends to either Runge-Kutta (RK) or corner-transport-upwind (CTU) time-stepping schemes (for the fluid) while a standard Boris integrator is employed for the particles. For highly-energetic relativistic CRs and in order to overcome the time step restriction a novel sub-cycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration and test particle acceleration in reconnecting layers are discussed.

There is considerable observational evidence of implosion of magnetic loop systems inside solar coronal active regions following high energy events like solar flares. In this work, we propose that such collapse can be modeled in three dimensions quite accurately within the framework of ideal magnetohydrodynamics. We furthermore argue that the dynamics of loop implosion is only sensitive to the transmitted disturbance of one or more of the system variables, e.g. velocity generated at the event site. This indicates that to understand loop implosion, it is sensible to leave the event site out of the simulated active region. Towards our goal, a velocity pulse is introduced to model the transmitted disturbance generated at the event site. Magnetic field lines inside our simulated active region are traced in real time, and it is demonstrated that the subsequent dynamics of the simulated loops closely resemble observed imploding loops. Our work highlights the role of plasma $\beta$ in regards to the rigidity of the loop systems and how that might affect the imploding loops' dynamics. Compressible magnetohydrodynamic modes such as kink and sausage are also shown to be generated during such processes, in accordance with observations.

An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat parabolic terms arising from the diffusion operator. An explicit treatment of parabolic terms will considerably reduce the simulation time step due to its dependence on the square of the grid resolution ($\Delta x$) for stability. Although an implicit scheme relaxes the constraint on stability, it is difficult to distribute efficiently on a parallel architecture. Treating parabolic terms with accelerated super-time stepping (STS) methods has been discussed in literature but these methods suffer from poor accuracy (first order in time) and also have difficult-to-choose tuneable stability parameters. In this work we highlight a second order (in time) Runge Kutta Legendre (RKL) scheme (first described by Meyer et. al. 2012) that is robust, fast and accurate in treating parabolic terms alongside the hyperbolic conversation laws. We demonstrate its superiority over the first order super time stepping schemes with standard tests and astrophysical applications. We also show that explicit conduction is particularly robust in handling saturated thermal conduction. Parallel scaling of explicit conduction using RKL scheme is demonstrated up to more than $10^4$ processors.

We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and not in others and we derive a simple estimate for the shock strength required for this to occur. These results are relevant to dense cores and Bok globules in star forming regions subjected to shocks produced by stellar feedback.

Magnetic reconnection is a plasma phenomenon where a topological rearrangement of magnetic field lines with opposite polarity results in dissipation of magnetic energy into heat, kinetic energy and particle acceleration. Such a phenomenon is considered as an efficient mechanism for energy release in laboratory and astrophysical plasmas. An important question is how to make the process fast enough to account for observed explosive energy releases. The classical model for steady state magnetic reconnection predicts reconnection times scaling as $S^{1/2}$ (where $S$ is the Lundquist number) and yields times scales several order of magnitude larger than the observed ones. Earlier two-dimensional MHD simulations showed that for large Lundquist number the reconnection time becomes independent of $S$ ("fast reconnection" regime) due to the presence of the secondary tearing instability that takes place for $S \gtrsim 1 \times 10^4$. We report on our 3D MHD simulations of magnetic reconnection in a magnetically confined cylindrical plasma column under either a pressure balanced or a force-free equilibrium and compare the results with 2D simulations of a circular current sheet. We find that the 3D instabilities acting on these configurations result in a fragmentation of the initial current sheet in small filaments, leading to enhanced dissipation rate that becomes independent of the Lundquist number already at $S \simeq 1\times 10^3$.

Recent high-resolution high-sensitivity observations of protostellar jets have shown many to possess deviations to their trajectories. HH 211 is one such example where sub-mm observations with the SMA have revealed a clear reflection-symmetric wiggle. The most likely explanation is that the HH 211 jet source could be moving as part of a protobinary system. Here we test this assumption by simulating HH 211 through 3D hydrodynamic jet propagation simulations using the PLUTO code with a molecular chemistry and cooling module, and initial conditions based on an analytical model derived from SMA observations. Our results show the reflection-symmetric wiggle can be recreated through the assumption of a jet source perturbed by binary motion at its base, and that a regular sinusoidal velocity variation in the jet beam can be close to matching the observed knot pattern. However, a more complex model with either additional heating from the protostar, or a shorter period velocity pulsation may be required to account for enhanced emission near the source, and weaker knot emission downstream. Position velocity diagrams along the pulsed jet beam show a complex structure with detectable signatures of knots and show caution must be exercised when interpreting radial velocity profiles through observations. Finally, we make predictions for future HH 211 observations with ALMA.

An Equation of State (\textit{EoS}) closes the set of fluid equations. Although an ideal EoS with a constant \textit{adiabatic index} $\Gamma$ is the preferred choice due to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Here, we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation as well as temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation- equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes.We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice-versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS.The first one employs root-finder methods and it is best suited for EoS in analytical form. The second one leans on lookup table and interpolation and results in a more computationally efficient approach although care must be taken to ensure thermodynamic consistency. A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is $500-10^4$ K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although using lookup table methods can significantly reduce the overhead by a factor $3\sim 4$.

We investigated a semi-analytic and numerical model to study the geometrically thin and optically thick accretion disk around Maclaurin spheroid (MS). The main interest is in the inner region of the so called {\alpha}-disk, {\alpha} being the viscosity parameter. Analytical calculations are done assuming radiation pressure and gas pressure dominated for close to Eddington mass accretion rate and $\dot{M}\lesssim 0.1\dot{M_{Edd}}$ respectively. We found that the change in eccentricity of MS gives a change at high frequency region in the emitted spectra. We found that disk parameters are dependent on eccentricity of MS. Our semi-analytic results show that qualitatively an increase in eccentricity of MS has same behavior as decrease in mass accretion rate. Numerical work has been carried out to see the viscous time evolution of the accretion disk around MS. In numerical model we showed that if the eccentricity of the object is high the matter will diffuse slowly during its viscous evolution. This gives a clue that how spin-up or spin-down can change the time evolution of the accretion disk using a simple Newtonian approach. The change in spectra can be used to determine the eccentricity of MS and thus period of the MS.

An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even for relatively weak shocks. Self-gravity plays only a minor role in determining the highest density that is reached, but it does play a role in the subsequent evolution. The densities at comparable times differ by a factor of a few for shocks initially propagating perpendicularly or obliquely to the magnetic field in the ambient medium.

Very long baseline interferometry (VLBI) observations of SiO masers in Orion Source I has enabled for the first time to resolve the outflow from a high-mass protostar in the launch and collimation region. Therefore, Source I provides a unique laboratory to study mass-loss and mass-accretion in a high-mass protostar. We numerically simulate the dynamics of the disk-wind inside 100 AU from Source I. This enables us to investigate the balance of different forces (gravitational, magnetic, thermal) regulating gas dynamics in massive star formation. In this work, we adopt magnetohydrodynamic (MHD) disk-wind models to explain the observed properties of the disk-wind from Orion Source I. The central source is assumed to be a binary composed of two 10\,$\msun$ stars in a circular orbit with an orbital separation of 7 AU. High resolution ideal MHD wind launching simulations (which prescribe disk as a boundary) are performed using the PLUTO code. The simulations are allowed to run until a steady state is obtained. MHD driven disk-wind provides a consistent model for the wide-angle flow from Source I probed by SiO masers, reproducing the bipolar morphology, the velocity amplitude and rotational profile, the physical conditions, and the magnetic field strength.

We investigate the launching of jets and outflows from magnetically diffusive accretion disks. Using the PLUTO code we solve the time-dependent resistive MHD equations taking into account the disk and jet evolution simultaneously. The main question we address is which kind of disks do launch jets and which kind of disks do not? In particular, we study how the magnitude and distribution of the (turbulent) magnetic diffusivity affect mass loading and jet acceleration. We have applied a turbulent magnetic diffusivity based on \alpha-prescription, but have also investigate examples where the scale height of diffusivity is larger than that of the disk gas pressure. We further investigate how the ejection efficiency is governed by the magnetic field strength. Our simulations last for up to 5000 dynamical time scales corresponding to 900 orbital periods of the inner disk. As a general result we observe a continuous and robust outflow launched from the inner part of the disk, expanding into a collimated jet of super fast magneto-sonic speed. For long time scales the disk internal dynamics changes, as due to outflow ejection and disk accretion the disk mass decreases. For magneto-centrifugally driven jets we find that for i) less diffusive disks, ii) a stronger magnetic field, iii) a low poloidal diffusivity, or a iv) lower numerical diffusivity (resolution), the mass loading of the outflow is increased - resulting in more powerful jets with high mass flux. For weak magnetization the (weak) outflow is driven by the magnetic pressure gradient. We further investigate the jet asymptotic velocity and the jet rotational velocity in respect of the different launching scenarios. We find a lower degree of jet collimation than previous studies, most probably due to our revised outflow boundary condition.

Observations indicate that outflows from massive young stars are more collimated during their early evolution compared to later stages. Our paper investigates various physical processes that impacts the outflow dynamics, i.e. its acceleration and collimation. We perform axisymmetric MHD simulations particularly considering the radiation pressure exerted by the star and the disk. We have modified the PLUTO code to include radiative forces in the line-driving approximation. We launch the outflow from the innermost disk region (r < 50 AU) by magneto-centrifugal acceleration. In order to disentangle MHD effects from radiative forces, we start the simulation in pure MHD, and later switch on the radiation force. We perform a parameter study considering different stellar masses (thus luminosity), magnetic flux, and line-force strength. For our reference simulation - assuming a 30 Msun star, we find substantial de-collimation of 35 % due to radiation forces. The opening angle increases from 20 deg to 32 deg for stellar masses from 20 Msun to 60 Msun. A small change in the line-force parameter 'alpha' from 0.60 to 0.55 changes the opening angle by ~ 8 deg. We find that it is mainly the stellar radiation which affects the jet dynamics. Unless the disk extends very close to the star, its pressure is too small to have much impact. Essentially, our parameter runs with different stellar mass can be understood as a proxy for the time evolution of the star-outflow system. Thus, we have shown that when the stellar mass (thus luminosity) increases (with age), the outflows become less collimated.

The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launching surface of the disk up to 6000^2 Schwarzschild radii allowing to reach highly relativistic Lorentz factors. The Poynting dominated disk wind develops into a jet with Lorentz factors of 8 and is collimated to 1 degree. In addition to the disk jet, we evolve a thermally driven spine jet, emanating from a hypothetical black hole corona. Solving the linearly polarized synchrotron radiation transport within the jet, we derive VLBI radio and (sub-) mm diagnostics such as core shift, polarization structure, intensity maps, spectra and Faraday rotation measure (RM), directly from the Stokes parameters. We also investigate depolarization and the detectability of a lambda^2-law RM depending on beam resolution and observing frequency. We find non-monotonic intrinsic RM profiles which could be detected at a resolution of 100 Schwarzschild radii. In our collimating jet geometry, the strict bi-modality in polarization direction (as predicted by Pariev et al.) can be circumvented. Due to relativistic aberration, asymmetries in the polarization vectors across the jet can hint to the spin direction of the central engine.

We investigate the structure of accretion disks around massive protostar applying steady state models of thin disks. The thin disk equations are solved with proper opacities for dust and gas taking into account the huge temperature variation along the disk. We explore a wide parameter range concerning stellar mass, accretion rate, and viscosity parameter \alpha . The most essential finding is a very high temperature of the inner disk. For e.g. a 10 M_sun protostar and an accretion rate of 10^-4 M_sun/yr, the disk midplane temperature may reach almost 10^5 K. The disk luminosity in this case is about 10^4 L_sun and, thus, potentially higher than that of a massive protostar. We motivate our disk model with similarly hot disks around compact stars. We calculate a dust sublimation radius by turbulent disk self-heating of more than 10AU, a radius, which is 3 times larger than caused by stellar irradiation. We discuss implications of this result on the flashlight effect and the consequences for the radiation pressure of the central star. In difference to disks around low mass protostars our models suggest rather high values for the disk turbulence parameter \alpha close to unity. However, disk stability to fragmentation due to thermal effects and gravitational instability would require a lower \alpha value. For \alpha = 0.1 we find stable disks out to 80AU. Essentially, our model allows to compare the outer disk to some of the observed massive protostellar disk sources, and from that, extrapolate on the disk structure close to the star which is yet impossible to observe.