Search

At optical/ultraviolet energies, blazars display an underlying thermal (unpolarized) contribution from the accretion disc, torus and line emitting regions, diluting the polarized emission from the jet-component. Optical polarimetry can be used to disentangle the thermal and non-thermal components, and place constraints on the particle populations and acceleration mechanisms responsible for the non-thermal emission. We present the results of a linear optical spectropolarimetric observing campaign of 18 blazars (6 BLLs and 12 FSRQs) undertaken with the Southern African Large Telescope between 2016 and 2022. This was done to observe these systems during flaring states, as well as long term monitoring of PKS1510-089, AP Lib and PKS 1034-293. The observations traced the frequency dependence of the degree and angle of polarization, as well as changes in the spectral line strengths. We investigated possible correlations between the polarization and other observed characteristics for the sources. While an indication of correlation was found between the frequency dependence and the average level of polarization for some sources, a correlation was not found for the population as a whole. These results highlight that continuous observations and in-depth modelling of polarization and its frequency dependence is required to obtain a more holistic view of TeV blazars., Comment: 17 pages, 21 figures, 3 tables, accepted for publication in MNRAS

The evolution of the spectral energy distribution during flares constrains models of particle acceleration in blazar jets. The archetypical blazar BL Lac provided a unique opportunity to study spectral variations during an extended strong flaring episode from 2020-2021. During its brightest $\gamma$-ray state, the observed flux (0.1-300 GeV) reached up to $2.15\,\times\,10^{-5}\,\rm{ph\,cm^{-2}\,s^{-1}}$, with sub-hour scale variability. The synchrotron hump extended into the X-ray regime showing a minute-scale flare with an associated peak shift of inverse-Compton hump in gamma-rays. In shock acceleration models, a high Doppler factor value $>$100 is required to explain the observed rapid variability, change of state, and $\gamma$-ray peak shift. Assuming particle acceleration in mini-jets produced by magnetic reconnection during flares, on the other hand, alleviates the constraint on required bulk Doppler factor. In such jet-in-jet models, observed spectral shift to higher energies (towards TeV regime) and simultaneous rapid variability arises from the accidental alignment of a magnetic plasmoid with the direction of the line of sight. We infer a magnetic field of $\sim0.6\,\rm{G}$ in a reconnection region located at the edge of BLR ($\sim0.02\,\rm{pc}$). The scenario is further supported by log-normal flux distribution arising from merging of plasmoids in reconnection region., Comment: 6 pages, 3 figures, 1 table, Accepted for publication in MNRAS-L

Radio relics are diffuse radio sources in galaxy clusters that are associated with merger shock waves. Detailed observations of radio relics in total intensity and in polarisation show complex structures on kiloparsec scales. The relation between the observed features and the underlying morphology of the magnetic field is not clear. Using three-dimensional magneto-hydrodynamical simulations, we study the polarised emission produced by a shock wave that propagates through a turbulent medium that resembles the intracluster medium. We model the polarised synchrotron emission on the basis of diffusive shock acceleration of cosmic-ray electrons. We find that the synchrotron emission produced in a shocked turbulent medium can reproduce some of the observed features in radio relics. Shock compression can give rise to a high polarisation fraction at the shock front and a partial alignment of the polarisation $E$-vectors with the shock normal. Our work confirms that radio relics can also be formed in an environment with a tangled magnetic field. We also discuss the effect of Faraday Rotation intrinsic to the source, and how our results depend on the angular resolution of observations., Comment: Accepted for publication in MNRAS, 22 pages, 23 figures

High-resolution radio observations of cluster radio relics often show complex spatial and spectral features. However, it is not clear what these features reveal about the underlying magnetic field properties. We performed three-dimensional magneto-hydrodynamical simulations of merger shock waves propagating through a magnetised, turbulent intracluster medium. Our model includes the diffusive shock acceleration of cosmic-ray electrons, their spatial advection and energy losses at run-time. With this set-up we can investigate the relation between radio substructure and pre-shock plasma conditions in the host cluster. We find that upstream turbulence plays a major role in shaping the properties of radio relics produced downstream. Within the assumption of diffusive shock acceleration, we can reproduce the observed discrepancy between the X-ray derived Mach number of shocks, and the Mach number inferred from radio spectra. Our simulated spectral index maps and profiles across the radio relic also suggest that the standard deviation of the upstream magnetic field must be relatively small ($\sigma_B\leq 1 \, \mu$G) in order to reproduce observations and therefore, radio relics can potentially constrain the distribution of magnetic fields in galaxy clusters outskirts., Comment: Accepted for publication in MNRAS, 22 pages, 25 figures

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., Comment: 22 pages, 13 figures

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., Comment: 4 Figures, 3 Tables, Accepted for publication in Astrophysics and Space Science Journal, Abstract shortened

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., Comment: 27 pages, 16 figures. Accepted for publication in ApJ Supplement series

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., Comment: 9 pages 9 figures

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$., Comment: 17 pages, 10 figures, Accepted for publication in A&A

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., Comment: 11 pages, 10 figures, accepted by MNRAS

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., Comment: 5 pages, 3 figures, accepted for publication in MNRAS

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., Comment: 5 pages, 4 figures, Accepted for publication in MNRAS Letters

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., Comment: 23 pages, 8 figures, accepted by ApJ, For full version see http://www.mpia.de/homes/vaidya/