Search

In this paper, the redshift evolution of the galactic bar properties, like the bar length, pattern speed, and bar fraction, has been investigated for simulated galaxies at stellar masses $M_*>10^{10}\, M_{\odot}$ in the cosmological magnetohydrodynamical simulation TNG50. We focus on the redshift evolution of the bar pattern speeds and \textit{the fast bar tension}. We show that the median value of the pattern speed of the bars increases as the redshift grows. On the other hand, although the median value of the bar length increases over time, the ratio between the corotation radius and the bar radius, namely the $\mathcal{R}=R_{\text{CR}}/R_{\text{bar}}$ parameter, increases as well. In other words, the corotation radius increases with a higher rate compared to the bar length. This directly means that galactic bars slow down with time, or equivalently as the redshift declines. We discuss the possible mechanisms that reduce the pattern speeds in TNG50. We demonstrate that while mergers can have a significant impact on a galaxy's pattern speed, they do not play a crucial role in the overall evolution of mean pattern speed within the redshift range $z\leq 1.0$. Furthermore, we show that the $\mathcal{R}$ parameter does not correlate with the gas fraction. Consequently, the existence of gas in TNG50 does not alleviate the fast bar tension. We show that the mean value of the pattern speed, computed for all the galaxies irrespective of their mass, at $z=1.0$ is $\Omega_p=70.98\pm 2.34$ km s$^{-1}$ kpc$^{-1}$ and reduces to $\Omega_p=33.65 \pm 1.07$ km s$^{-1}$ kpc$^{-1}$ at $z=0.0$. This is a direct prediction by TNG50 that bars at $z=1.0$ rotate faster by a factor of $\sim 2$ compared to bars at $z=0.0$., Comment: 13 pages, Accepted for publication in A&A

We present the structure of a low angular momentum accretion flows around rotating compact objects incorporating relativistic corrections up to the leading post-Newtonian order. To begin with, we formulate the governing post-Newtonian hydrodynamic equations for the mass and energy-momentum flux without imposing any symmetries. However, for the sake of simplicity, we consider the flow to be stationary, axisymmetric, and inviscid. Toward this, we adapt the polytropic equation of state (EoS) and analyze the vertically integrated accretion flow confined to the equatorial plane. It is shown that the spin-orbit effects manifest themselves in the accretion dynamics. In the present analysis, we focus on global transonic accretion solutions, where a subsonic flow enters far away from the compact object and gradually gains radial velocity as it moves inwards. Thus, the flow becomes supersonic after reaching a certain radius, known as the critical point. To better understand the transonic solutions and examine the effect of post-Newtonian corrections, we classify the post-Newtonian equations into semi-relativistic (SR), semi-Newtonian (SN), and non-relativistic (NR) limits and compare the accretion solutions and their corresponding flow variables. With these, we find that SR and SN flow are in good agreement all throughout, although they deviate largely from the NR ones. Interestingly, the density profile seems to follow the profile $\rho \propto r^{-3/2}$ in the post-Newtonian regime. The present study has the potential to connect Newtonian and GR descriptions of accretion dynamics., Comment: 42 pages, 9 figures, and 5 appendices, Accepted for publication in JCAP

Observational evidence and many numerical simulations show the existence of wind (i.e., uncollimated outflow) in accretion systems of the elliptical galaxy center. One of the primary aims of this study is to investigate the solutions of slowly rotating accretion flows around the supermassive black hole with outflow. This paper presents two distinct physical regions: supersonic and subsonic, that extend from the outer boundary to the black hole. In our numerical solution, the outer boundary is chosen beyond the Bondi radius. Due to strong gravity, we ignore outflow (i.e., $s = 0$) in the inner region (within $\sim 10 r_s$). The radial velocity of the flow at the outer region is significantly increased due to the presence of the outflow. Compared to previous works, the one accretion mode - namely, slowly rotating case, which corresponds to low accretion rates that have general wind output is carefully described, and the effect of galaxy potential and the feedback effects by the wind in this mode are taken into account. Since the power-law form of the mass accretion rate is mathematically compatible with our equations, we consider a radius-dependent mass accretion rate ($\dot{M}_{in} \propto r^s$), where $s$ is a free parameter and shows the intensity of outflow. There is an unknown mechanism for removing the mass, angular momentum, and energy by outflows in this study. The effects of the outflow appear well on the outer edge of the flow., Comment: Accepted for publication in Astrophysical Journal

We examine the effect of thermal conduction on the low-angular momentum hot accretion flow (HAF) around non-rotating black holes accreting mass at very low rate. While doing so, we adopt the conductive heat flux in the saturated form, and solve the set of dynamical equations corresponding to a steady, axisymmetric, viscous, advective accretion flow using numerical methods. We study the dynamical and thermodynamical properties of accreting matter in terms of the input parameters, namely energy ($\varepsilon_0$), angular momentum ($\ell_0$), viscosity parameter ($\alpha$), and saturation constant ($\Phi_{\rm s}$) regulating the effect of thermal conduction. We find that $\Phi_{\rm s}$ plays a pivotal role in deciding the transonic properties of the global accretion solutions. In general, when $\Phi_{\rm s}$ is increased, the critical point ($r_{\rm c}$) is receded away from the black hole, and flow variables are altered particularly in the outer part of the disc. To quantify the physically acceptable range of $\Phi_{\rm s}$, we compare the global transonic solutions with the self-similar solutions, and observe that the maximum saturation constant ($\Phi^{\rm max}_{\rm s}$) estimated from the global solutions exceeds the saturated thermal conduction limit ($\Phi_{\rm sc}$) derived from the self-similar formalism. Moreover, we calculate the correlation between $\alpha$ and $\Phi^{\rm max}_{\rm s}$ and find ample disagreement between global solutions and self-similar solutions. Further, using the global flow variables, we compute the Bernoulli parameter ($Be$) which remains positive all throughout the disc, although flow becomes loosely unbound for higher $\Phi_{\rm s}$. Finally, we indicate the relevance of this work in the astrophysical context in explaining the possibility of massloss/outflows from the unbound disc., Comment: 12 pages, 8 figures; accepted for publication in MNRAS

We investigate the dynamics of clumps that coexisted with/in advection-dominated accretion flows by considering thermal conductivity. Thermal conduction can be one of the effective factors in the energy transportation of ADAFs; hence it may indirectly affect the dynamics of clumps by means of a contact force between them and their host medium. We first study the ensemble of clumps by assuming them as collision-less particles and secondly we find the orbital motion of these clouds as individuals. For both parts, clumps are subject to the gravity of the central object and a drag force. The strong coupling between clumps and ADAF leads to equality between the average treatment of the clumps and the dynamics of their background. By employing the collision-less Boltzmann equation we calculate the velocity dispersion of the clumps which turns out approximately one order of magnitude higher than the ADAF. In fact, involving drag force in such a system causes the angular momentum of the clumps can be transported outwards by the ADAF, and hence the clouds eventually will be captured at the tidal radius. The results show that the presence of thermal conduction increases the root of the averaged radial velocity square and this in turn speeds up the process of capturing the clouds through the tidal force. In the end, we focus on a typical individual cloud, the spiral orbits appear thanks to only the toroidal component of friction force. The parametric study again proves that the operation of thermal conduction helps for decreasing the lifetime of clumps., Comment: Accepted for publication in MNRAS

We study the global solutions of slowly rotating accretion flows around the supermassive black hole in the nucleus of an elliptical galaxy. The velocity of accreted gas surrounding the black hole is initially subsonic and then falls onto the black hole supersonically, so accretion flow must be transonic. We numerically solve equations from the Bondi radius to near the black hole. The focus of our discussion will be on the properties of slightly rotating accretion flows in which radiative losses have been ignored. This study discusses how outer boundary conditions (the temperature and specific angular momentum at the outer boundary) influence accretion flow dynamics. We investigate two physically discontinuous regimes: The Bondi-like type accretion and the Disk-like type accretion. A Bondi-like accretion occurs when the specific angular momentum at the Bondi radius $ \ell_{B} $ is smaller than the specific angular momentum at the marginally stable orbit $ \ell_{ms} $. In comparison, a Disk-like accretion occurs when the specific angular momentum at the Bondi radius $ \ell_{B} $ is larger than the specific angular momentum of the marginally stable orbit $ \ell_{ms} $. We also keep the assumption of hydrostatic equilibrium and compare our results with the case in which it is not considered. According to this study, considering the assumption of hydrostatic equilibrium reduces the mass accretion rate. Additionally, we find our solution for different ranges of the viscosity parameter $\alpha$. Finally, we study the effect of galaxy potential on slowly rotating accretion flows., Comment: Accepted for publication in MNRAS

The present work proposes a new formalism for the inner regions of a neutrino-dominated accretion flows (NDAFs) by considering the self-gravity, where the neutrino opacity is high enough to make neutrinos trapped becoming a dominant factor in the transportation of energy and angular momentum over the magneto rotational instability. We investigate the possibility of gravitational instability and fragmentation to model the highly variable structure of the prompt emission in gamma-ray bursts (GRBs). The results lead us to introduce the gravitational instability, in these inner regions, as a source of a new viscosity which is of the same functional form as that of the $\beta$-prescription of viscosity. Such a consideration brings about fragmentation in the unstable inner disk. In addition, we find the consequent clumpy structure of this area capable to account for the temporal variability of GRB's light curve, especially for the lower choices of the parameter $\beta$, $\sim 10^{-5}$. Finally, we predict the formation of gravitational waves through the migration of fragments before being tidally disrupted. These waves appear to be detectable via a range of current and future detectors from LIGO to Cosmic Explorer., Comment: Accepted for publication in MNRAS

Filamentary molecular clouds are omnipresent in the cold interstellar medium. Observational evidences show that the non-isothermal equations of state describe the filaments properties better than the isothermal one. In this paper we use the logatropic and the polytropic equations of state to study the gravitational instability of the pressure-confined filaments in presence of a uniform axial magnetic field. To fully explore the parameter space we carry out very large surveys of stability analysis that cover filaments with different radii in various magnetic fields. Our results show that for all the equations of state the instability of thinner filaments is more sensitive to the magnetic field variations than the thicker ones. Moreover, for all the equations of state, an intermediate magnetic field can entirely stabilize the thinner filaments. Albeit for the thicker ones this effect is suppressed for the magnetic field stronger than B ' 70 micro G., Comment: Accepted for publication in MNRAS

The main purpose of this paper is to obtain analytical solutions for radiative transfer equations related to the vertical structure of accretion discs with finite optical depth. In the non-gray atmosphere, we employ the optical-depth dependent Eddington factor to define the relationship between the mean intensity and radiation stress tensor. Analytical solutions are achieved for two cases: (i) radiative equilibrium, and (ii) a disc with uniform internal heating and both cases are assumed to be in local thermodynamical equilibrium (LTE), too. These solutions enable us to study the probable role of scattering and disc optical depth on the emergent intensity and other radiative quantities. Our results show that for the first case, the surface value of mean intensity with constant Eddington factor is three times larger than that with a variable factor. Moreover, scattering has no role in the vertical radiative structure of discs with the assumptions of the first case. On the other hand, for the second case, we encounter reductions in all radiative quantities as the photon destruction probability decreases (which is equivalent to increasing scattering). Furthermore, for both cases with total optical depth less than unity, the outward intensity towards the polar direction becomes less than that from the edges of disc which is contrary to limb-darkening. In the end, we apply our results to find the spectrum from accretion systems, based on two dynamical models. Consequently, we can see how the total optical depth varies with frequency and causes remarkable changes in the emergent spectra., Comment: Accepted for publication in MNRAS

The evolution of disk galaxies in modified gravity is studied by using high-resolution N-body simulations. More specifically, we use the weak field limit of two modified gravity theories, i.e., nonlocal gravity (NLG) and scalar-tensor-vector gravity known as MOG, and ignore the existence of dark matter halo. On the other hand, we construct the same models in the standard dark matter model and compare their dynamics with the galactic models in modified gravity. It turns out that there are serious differences between galactic models in these different viewpoints. For example, we explicitly show that the galactic models in modified gravity, host faster bars compared to the dark matter case. On the other hand, final stellar bars are weaker in modified gravity. These facts are not new and have already been reported in our previous simulations for exponential galactic models. Therefore, our main purpose in this study is to show that the above-mentioned differences, with emphasis on the speed of the bars, are independent of the initial density profile of the adopted disk/halo. To do so, we employ different profiles for the disk and halo and show that the results remain qualitatively independent of the initial galactic models. Moreover, a more accurate method has been used to quantify the kinematic properties of the stellar bar. Our results imply that contrary to the dark matter models, bars in modified gravity are fast rotators which never leave the fast-bar region until the end of the simulation., Comment: Accepted for publication in ApJ

X-ray flares can be accounted for a hint to the late time activity of Gamma-ray bursts (GRBs) central engines. Such a long term activity has been described through some models, one of which is the viscous evolution of the outer disc fragments that proposed by Perna et al. (2006), and developed quantitatively by Dall Osso et al. (2017). Here, we reconstruct Dall Osso et al. (2017) framework through taking both small and large scale effects of magnetic field into account. To consider the magnetic barrier as a possible mechanism that might govern the accretion process of each magnetized clump, we make a simple pattern in boundary condition through which this mechanism might happen. Regarding various model parameters, we probe for their influence and proceed some key analogies between our model predictions and previous phenomenological estimates, for two different choices of boundary conditions (with and without a magnetic barrier). Our model is remarkably capable of matching flare bolometric and X-ray light-curves, as well as reproducing their statistical properties, such as the ratios between rise and decay time, width parameter and peak time, and the power-law correlation between peak luminosity and peak time. Combining our results with the conclusions of previous studies, we are led to interpret magnetic barrier as a less probable mechanism that might control the evolution of these clumps, especially the later created (or viscously evolved) ones., Comment: Accepted for publication in ApJ

We study the local thermal stability of thin accretion disks. We present a full stability analysis in the presence of a magnetic field and more importantly wind. For wind, we use a general model suitable for adequately describing several kinds of winds. First, we explicitly show that the magnetic field, irrespective of the type of wind, has a stabilizing effect. This is also true when there is no wind. In this case, we confirm the other works already presented in the literature. However, our main objective is to investigate the local thermal stability of the disk in the presence of the wind. In this case, interestingly, the response of disk is directly related to the type of wind. In other words, in some cases, the wind can stabilize the disk. On the other hand, in some cases, it can destabilize the disk. We found that in some thin disk models where the magnetic pressure cannot explain the stability of the disk by including a typical contribution for magnetic pressure, the wind can provide a viable explanation for the thermal stability., Comment: Accepted for publication in ApJ

By taking into account photon absorption, we investigate the vertical structure of accretion flows with comparable radiation and gas pressures. We consider two separate energy equations for matter and radiation in the diffusion limit. In order to solve the set of radiation hydrodynamic equations in steady state and axisymmetric configuration, we employ self-similar technique in the radial direction. We need the reflection symmetry about the mid-plane to find gas density at the equator. For a typical solution, we assume that the gas pressure has 10-50\% portion of the total pressure. In this paper, since the radiation energy is involved directly, we are able to estimate how much energy of viscous heating is transported in the radial direction and advected towards the central object. Our results show that although the mass accretion rate does not approach the Eddington limit, the energy advection is rather high. Moreover, in a disc with greater accretion rate and less portion of gas pressure at the total pressure, more energy is advected to its center. In addition, as we expect the accretion flow becomes thicker with greater values of gas pressure. Based on Solberg-H\~A\v{z}iland conditions, we notice that the flow is convectively stable in all parts of such a disc.

We implement a new semi-analytical approach to investigate radially self-similar solutions for the steady-state advection-dominated accretion flows (ADAFs). We employ the usual $\alpha$-prescription for the viscosity and all the components of the energy-momentum tensor are considered. In this case, in the spherical coordinate, the problem reduces to a set of eighth-order, nonlinear differential equations with respect to the latitudinal angle $\theta$. Using the Fourier expansions for all the flow quantities, we convert the governing differential equations to a large set of nonlinear algebraic equations for the Fourier coefficients. Using the Newton-Raphson method we solve the algebraic equations and ADAF properties are investigated over a wide range of the model parameters. We also show that the implemented series are truly convergent. The main advantage of our numerical method is that it does not suffer from the usual technical restrictions which may arise for solving ADAF differential equations near to the polar axis. In order to check the reliability of the approach, we recover some widely studied solutions. We also introduce a new varying $\alpha$ viscosity model and new outflow and inflow solutions for ADAFs are presented using Fourier expansion series., Comment: Accepted for publication in ApJ. Comments welcome

Recent observations of the filamentary molecular clouds show that their properties deviate from the isothermal equation of state. Theoretical investigations proposed that the logatropic and the polytropic equations of state with negative indexes can provide a better description for these filamentary structures. Here, we aim to compare the effects of these softer non-isothermal equation of states with their isothermal counterpart on the global gravitational instability of a filamentary molecular cloud. By incorporating the ambipolar diffusion, we use the non-ideal magnetohydrodynamics framework for a filament that is threaded by a uniform axial magnetic field. We perturb the fluid and obtain the dispersion relation both for the logatropic and polytropic equations of state by taking the effects of magnetic field and ambipolar diffusion into account. Our results suggest that, in absence of the magnetic field, a softer equation of state makes the system more prone to gravitational instability. We also observed that a moderate magnetic field is able to enhance the stability of the filament in a way that is sensitive to the equation of state in general. However, when the magnetic field is strong, this effect is suppressed and all the equations of state have almost the same stability properties. Moreover, we find that for all the considered equations of state, the ambipolar diffusion has destabilizing effects on the filament., Comment: Accepted for publication in MNRAS, 11 pages, 6 figures

We solved the set of two-dimensional magnetohydrodynamic (MHD) equations for optically thin black hole accretion flows incorporating toroidal component of magnetic field. Following global and local MHD simulations of black hole accretion disks, the magnetic field inside the disk is decomposed into a large scale field and a fluctuating field. The effects of the fluctuating magnetic field in transferring the angular momentum and dissipating the energy are described through the usual $ \alpha $ description. We solved the MHD equations by assuming steady state and radially self-similar approximation in $ r-\theta $ plane of spherical coordinate system. We found that as the amount of magnetic field at the equatorial plane increases, the heating by the viscosity decreases. In addition, the maximum amount of the heating by the viscous dissipation is produced at the mid-plane of the disk, while that of the heating by the magnetic field dissipation is produced at the surface of the disk. Our main conclusion is that in terms of the no-outflow solution, thermal equilibrium still exists for the strong magnetic filed at the equatorial plane of the disk., Comment: Accepted for publication in The Astrophysical Journal

We investigate the properties of a hot accretion flow bathed in a poloidal magnetic field. We consider an axisymmetric viscous resistive flow in the steady state configuration. We assume the dominant mechanism of energy dissipation is due to turbulence viscosity and magnetic diffusivity. A certain fraction of that energy can be advected towards the central compact object. We employ self-similar method in the radial direction to find a system of ODEs with just one varible, $\theta$ in the spherical coordinates. For the existence and maintaining of a purely poloidal magnetic in a rotating thick disk, we find the necessary condition is a constant value of angular velocity along a magnetic field line. We obtain an analytical solution for the poloidal magnetic flux. We explore possible changes in the vertical structure of the disk under the influences of even symmetric and asymmetric magnetic fields. Our results reveal that a polar magnetic field with even symmetry about the equatorial plane makes the disk vertically thin. Moreover, the accretion rate decreases when we consider a strong magnetic field. Finally, we notice hot magnetized accretion flows can be fully advected even in a slim shape., Comment: Accepted for publication in Astrophysical Journal

What is conducted in the present work is the study of self-gravity effects on vertical structure of magnetized neutrino-dominated accretion disc (termed NDAF and considered as a central engine for Gamma-ray bursts (GRBs)). Some of the disc physical time scales such as viscous, cooling and diffusion time scales, that are supposed to play a pivotal role in the late time evolutions of the disc, have been studied. It is of our interest to investigate the possibility of the X-ray flares occurrence, observed in late time GRB's extended emission through the "magnetic barrier" and "fragmentation" processes in our model. The results lead us to interpret self-gravity as an amplifier for Blandford-Payne luminosity (BP power) and the generated magnetic field, but a suppressor for neutrino luminosity and magnetic barrier process via highlighting fragmentation mechanism in the outer disc, especially for the higher mass accretion rates., Comment: 11 pages, 7 figures, accepted at ApJ

The recent observations of rippled structures on the surface of the Orion molecular cloud (Bern\'{e} et al. 2010), have been attributed to the Kelvin-Helmholtz (KH) instability. The wavelike structures which have mainly seen near star-forming regions taking place at the interface between the hot diffuse gas, which is ionized by massive stars, and the cold dense molecular clouds. The radiation pressure of massive stars and stellar clusters is one of the important issues that has been considered frequently in the dynamics of clouds. Here, we investigate the influence of radiation pressure, from well-known Trapezium cluster in the Orion nebula, on the evolution of KH instability. The stability of the interface between HII region and molecular clouds in the presence of the radiation pressure, has been studied using the linear perturbation analysis for the certain range of the wavelengths. The linear analysis show that consideration of the radiation pressure intensifies the growth rate of KH modes and consequently decreases the e-fold time-scale of the instability. On the other hand the domain of the instability is extended and includes the more wavelengths, consisting of smaller ones rather than the case when the effect of the radiation pressure is not considered. Our results shows that for $\lambda_{\rm KH}>0.15\rm pc$, the growth rate of KH instability dose not depend to the radiation pressure. Based on our results, the radiation pressure is a triggering mechanism in development of the KH instability and subsequently formation of turbulent sub-structures in the molecular clouds near massive stars. The role of magnetic fields in the presence of the radiation pressure is also investigated and it is resulted that the magnetic field suppresses the effects induced by the radiation pressure., Comment: Accepted for publication in MNRAS

We examine the effects of a global magnetic field and outflow on radiatively inefficient accretion flow (RIAF) in the presence of magnetic resistivity. We find a self-similar solutions for the height integrated equations that govern the behavior of the flow. We use the mixing length mechanism for studying the convection parameter. We adopt a radius dependent mass accretion rate as $\dot{M}=\dot{M}_{out}{(\frac{r}{r_{out}})^{s}}$ with $s> 0$ to investigate the influence of outflow on the structure of inflow where $s$ is a constant and indication the effect of wind. Also, we have studied the radiation spectrum and temperature of CDAFs. The thermal bermsstrahlung emission as a radiation mechanism is taken into account for calculating the spectra emitted by the CDAFs. The energy that powers bremsstrahlung emission at large radii is provided by convective transport from small radii and viscous and resistivity dissipation. Our results indicate that the disc rotates slower and accretes faster, it becomes hotter and thicker for stronger wind. By increasing all component of magnetic field, the disc rotates faster and accretes slower while it becomes hotter and thicker. We show that the outflow parameter and all component of magnetic field have the same effects on the luminosity of the disc. We compare the dynamical structure of the disc in two different solutions (with and without resistivity parameter). We show that only the radial infall velocity and the surface density could changed by resistivity parameter obviously. Increasing the effect of wind increases the disc's temperature and luminosity of the disc. The effect of magnetic field is similar to the effect of wind in the disc's temperature and luminosity of the disc, but the influence of resistivity on the observational properties is not evident.

We present the effects of ordered large scale magnetic field on the structure of supercritical accretion flow in the presence of outflow. In the cylindrical coordinates ($r, \varphi,z$), We write the 1.5 dimensional, the steady state ($\frac{\partial}{\partial t}=0$) and axisymmetric ($\frac{\partial}{\partial \varphi}=0$) inflow-outflow equations by using the self similar solutions. Also a model for radiation pressure supported accretion flow threaded by both toroidal and vertical components of magnetic field has been formulated. For studying the outflows, we adopt a radius dependent mass accretion rate as $\dot{M}=\dot{M}_{out}{(\frac{r}{r_{out}})^{s+\frac{1}{2}}}$ with $s=\frac{1}{2}$. Also by following the previous works, we have considered the interchange of mass, radial and angular momentum and the energy between inflow and outflow. we have found numerically that two components of magnetic field have the opposite effects on the thickness of the disc and similar effects on the radial and angular velocities of the flow. We have found that the existence of the toroidal component of magnetic field will lead to increasing of radial and azimuthal velocities as well as the relative thickness of the disc, which is increased . Moreover, the thickness of the disc will decrease when the vertical component of magnetic field becomes important in magnetized flow. The solutions indicated that the mass inflow rate, the specific energy of outflow affect strongly on the advection parameter.We have shown that by increasing the two components of magnetic field, the temperature of the accretion flow will decrease significantly. On the other hand we have shown the bolometric luminosity of the slim discs for high values of $\dot{m} (\dot{m}>>1)$ is not sensitive to mass accretion rate and is kept constant ($L \approx 10 L_{E}$)., Comment: Accepted for publication in MNRAS

This paper explores the effects of large-scale magnetic fields in hot accretion flows for asymmetric configurations with respect to the equatorial plane. The solutions that we have found show that the large-scale asymmetric magnetic field can significantly affect the dynamics of the flow and also cause notable outflows in the outer parts. Previously, we treated a viscous resistive accreting disc in the presence of an odd symmetric $\textbf{B-}$ field about the equatorial plane. Now we extend our earlier work by taking into account another configuration of large-scale magnetic field which is no longer symmetric. We provide asymmetric field structures with small deviations from even and odd symmetric $\textbf{B-}$field. Our results show that the disc's dynamics and appearance become different above and below the equatorial plane. The set of solutions also predicts that even a small deviation in a symmetric field causes the disc to compress on one side and expand on the other. In some cases, our solution represents a very strong outflow from just one side of the disc. Therefore, the solution may potentially explain the origin of one-sided jets in radio galaxies., Comment: accepted for publication in MNRAS

The Jeans analysis is studied in the first post-Newtonian limit. In other words, the relativistic effects on the local gravitational instability are considered for systems where characteristic velocity of the system and corresponding gravitational field are higher than what permitted in Newtonian limit. The dispersion relation for propagation of small perturbations is found in the post-Newtonian approximation using two different techniques. A new Jeans mass is derived and compared to the standard Jeans mass. In this limit, the relativistic effects make the new Jeans mass to be smaller than the Newtonian Jeans mass. Furthermore, the fractional difference between these two masses increases when temperature/pressure of the system increases. Interestingly, in this limit pressure can help the gravitational instability instead of preventing it. Finally the results are applied to high temperature astrophysical systems and the possibility of local fragmentations in some relativistic systems is investigated., Comment: 14 pages, 5 figures, accepted for publication in ApJ

The gravitational instability of a filamentary molecular cloud in non-ideal magnetohydrodynamics is investigated. The filament is assumed to be in hydrostatic equilibrium. We add the effect of ambipolar diffusion to the filament which is threaded by an initial uniform axial magnetic field along its axis. We write down the fluid equations in cylindrical coordinates and perform linear perturbation analysis. We integrate the resultant differential equations and then derive the numerical dispersion relation. We find that, a more efficient ambipolar diffusion leads to an enhancement of the growth of the most unstable mode, and to increase of the fragmentation scale of the filament., Comment: Accepted for publication in MNRAS, 10 pages, 7 figures

We obtain two-dimensional exact analytic solutions for the structure of the hot accretion flows without wind. We assume that the only non-zero component of the stress tensor is $T_{r\varphi}$. Furthermore we assume that the value of viscosity coefficient $\alpha$ varies with $\theta$. We find radially self-similar solutions and compare them with the numerical and the analytical solutions already studied in the literature. The no-wind solution obtained in this paper may be applied to the nuclei of some cool-core clusters., Comment: to appear in MNRAS

We study the global gravitational stability of a gaseous self-gravitating Maclaurin disk in the absence of a halo. Further, we replace Newtonian gravity with the specific Modified gravity theory known as MOG in the relevant literature. MOG is an alternative theory for addressing the dark matter problem without invoking exotic dark matter particles, and possesses two free parameters $\alpha$ and $\mu_0$ in the weak field limit. We derive the equilibrium gravitational potential of the Maclaurin disk in MOG and develop a semi-analytic method for studying the response of the disk to linear non-axysymmetric perturbations. The eigenvalue spectrum of the normal modes of the disk is obtained and its physical meaning has been explored. We show that Maclaurin disks are less stable in MOG than in Newtonian gravity. In fact both parameters $(\alpha,\mu_0)$ have destabilizing effects on the disk. Interestingly $\mu_0$ excites only the bar mode $m=2$ while $\alpha$ affects all the modes. More specifically, when $\alpha>1$, the bar mode is strongly unstable and unlike in Newtonian gravity can not be avoided, at least in the weak field limit, with increasing the pressure support of the disk., Comment: to appear in ApJ

We solve the radiation-hydrodynamic (RHD) equations of supercritical accretion flows in the presence of radiation force and outflow by using self similar solutions. Compare with the pioneer works, in this paper we consider power-law function for mass inflow rate as $ \dot{M} \propto r^{s} $. We found that $ s = 1 $ when the radiative cooling term is included in the energy equation. Correspondingly, the effective temperature profile with respect to the radius was obtained as $ T_{\text{eff}} \varpropto r^{-1/2} $. In addition, we investigated the influence of the outflow on the dynamics of the accretion flow. We also calculated the continuum spectrum emitted from the disk surface as well as the bolometric luminosity of the accretion flow. Furthermore, our results show that the advection parameter, $ f $, strongly depends on mass inflow rate., Comment: Accepted for publication in Astrophysical Journal

We construct self-similar inflow-outflow solutions for a hot viscous-resistive accretion flow with large scale magnetic fields that have odd symmetry with respect to the equatorial plane in $B_\theta$, and even symmetry in $B_r$ and $B_\phi$. Following previous authors, we also assume that the polar velocity $v_\theta$ is nonzero. We focus on four parameters: $\beta_{r0}$, $\beta_{\phi0}$ (the plasma beta parameters for associated with magnetic field components at the equatorial plane), the magnetic resistivity $\eta_0$, and the density index $n=-d\ln\rho/d\ln r$. The resulting flow solutions are divided into two parts consisting of an inflow region with a negative radial velocity ($v_r<0$) and an outflow region with $v_r>0$. Our results show that stronger outflows emerge for smaller $\beta_{r0}$ ($\le10^{-2}$ for $n>1$) and larger values of $\beta_{\phi0}$, $\eta_0$ and $n$., Comment: 10 pages, 9 figures, Accepted for publication MNRAS

We present the effect of a hydrodynamical wind on the structure and the surface temperature of a vertically self-gravitating magnetized ADAFs using self-similar solutions. Also a model for an axisymmetric, steady-state, vertically self-gravitating hot accretion flow threaded by a toroidal magnetic field has been formulated. The model is based on $\alpha-$prescription for turbulence viscosity. It is found that the thickness and radial velocity of the disc are reduced significantly as wind gets stronger. In particular, the solutions indicated that the wind and advection have the same effects on the structure of the disc. We also find that in the optically thin ADAF becomes hotter by including the wind parameter and the self-gravity parameter., Comment: 7 pages, 6 Figures, Accepted for publication MNRAS

We solve the set of hydrodynamic (HD) equations for optically thin Advection Dominated Accretion Flows (ADAFs) by assuming radially self-similar in spherical coordinate system $ (r, \theta, \phi) $. The disk is considered to be steady state and axi-symmetric. We define the boundary conditions at the pole and the equator of the disk and to avoid singularity at the rotation axis, the disk is taken to be symmetric with respect to this axis. Moreover, only the $ \tau_{r \phi} $ component of viscous stress tensor is assumed and we have set $ v_{\theta} = 0 $. The main purpose of this study is to investigate the variation of dynamical quantities of the flow in the vertical direction by finding an analytical solution. As a consequence, we found that the advection parameter, $ f^{adv} $, varies along the $ \theta $ direction and reaches to its maximum near the rotation axis. Our results also show that, in terms of no-outflow solution, thermal equilibrium still exists and consequently advection cooling can balance viscous heating., Comment: Accepted for publication to ApJ

In the framework of metric $f(R)$ gravity, we find the dispersion relation for the propagation of tightly wound spiral density waves in the surface of rotating, self-gravitating disks. Also, new Toomre-like stability criteria for differentially rotating disks has been derived for both fluid and stellar disks., Comment: To appear in Astrophysics and Space Science

We find the dispersion relation for tightly wound spiral density waves in the surface of rotating, self-gravitating disks in the framework of Modified Gravity (MOG). Also, the Toomre-like stability criterion for differentially rotating disks has been derived for both fluid and stellar disks. More specifically, the stability criterion can be expressed in terms of a matter density threshold over which the instability occurs. In other words the local stability criterion can be written as $\Sigma_0<\Sigma_{\text{crit}}(v_s,\kappa,\alpha,\mu_0)$, where $\Sigma_{\text{crit}}$ is a function of $v_s$ (sound speed), $\kappa$ (epicycle frequency) and $\alpha$ and $\mu_0$ are the free parameters of the theory. In the case of a stellar disk the radial velocity dispersion $\sigma_r$ appears in $\Sigma_{\text{crit}}$ instead of $v_s$. We find the exact form of the function $\Sigma_{\text{crit}}$ for both stellar and fluid self-gravitating disks. Also, we use a sub-sample of THINGS catalog of spiral galaxies in order to compare the local stability criteria. In this perspective, we have compared MOG with Newtonian gravity and investigated the possible and detectable differences between these theories., Comment: 9 pages, 5 figures, 1 table

MOdified Gravity (MOG) is a covariant modification of Einstein's general relativity. This theory is one of the current alternatives to dark matter models. We describe dynamics of collisionless self-gravitating systems in the context of MOG. By studying the weak field approximation of this theory, we derive the equations governing the dynamics of the self-gravitating systems. More specifically, we consider the Jeans instability for self-gravitating fluid and stellar systems, and derive new Jeans mass limit $\tilde{M}_J$ and wave-number $\tilde{k}_J$. Furthermore, considering the gravitational instability in star forming regions, we show that MOG has not a significant difference with general relativity on this astrophysical scale. However, at larger scales such as intergalactic space MOG may lead to different galaxy and structure formation processes.

Dynamics of clumps within a magnetized advection dominated accretion flow is investigated by solving the collisionless Boltzmann equation and considering the drag force due to the relative velocity between the clumps and the gas. Toroidal component of the magnetic field is assumed to be dominant. Dynamical properties of the hot gaseous component such as the radial and the rotational velocities are affected by the magnetic effects, and so, drag force varies depending on the strength of magnetic field and the velocity dispersion of the clumps is then modified significantly. We show that when magnetic pressure is less than the gas pressure, the root of the averaged radial velocity square of the clumps decreases at the inner parts of the hot flow and increases slightly at its outer edge., Comment: 6 pages, Accepted for publication in MNRAS

By taking into account the effect of toroidal magnetic field and its correspond heating, we determine the thickness of advection-dominated accretion flows. We consider an axisymetric, rotating, steady viscous-resistive, magnetized accretion flow under an advection dominated stage. The dominant mechanism of energy dissipation is assumed to be turbulence viscosity and magnetic diffusivity. We adopt a self-similar assumption in the radial direction to obtain the dynamical quantities, that is, radial, azimuthal, sound and Alfv\' en velocities. Our results show the vertical component of magnetic force acts in the opposite direction of gravity and compresses the disc, thus compared with the non-magnetic case, in general the disc half-thikness, $\Delta\theta$, significantly is reduced. On the other hand, two parameters, appearing due to action of magnetic field and reaction of the flow, affect the disc thickness. The first one, $\beta_0$, showing the magnetic field strength at the equatorial plane, decreases $\Delta\theta$. The other one, $\eta_0$ is the magnetic resistivity parameter and when it increases, $\Delta\theta$ increases, too., Comment: To appear in MNRAS

The main aim of this paper is studying the effect of toroidal magnetic field on the structure of Advection-Dominated Accretion Flows (ADAF) in the presence of the turbulence viscosity and diffusivity due to viscosity and magnetic field respectively. We use self-similar assumption in radial direction to solve the magnetohydrodynamic (MHD) equations for hot accretion disk. We use spherical coordinate $ (r, \theta, \varphi) $ to solve our equation. The toroidal component of magnetic field is considered and all three components of the velocity field $ \mathbf{v}\equiv (v_{r}, v_{\theta}, v_{\varphi}) $ are present in our work. We reduce the equations to a set of differential equations about $ \theta $ and apply the symmetric boundary condition at the equatorial plane of the disk. Our results indicate that the outflow region, where the redial velocity becomes positive in a certain inclination angle $ \theta_{0} $, always exist. The results represent that the stronger the magnetic field, the smaller the inclination angle, $ \theta_{0} $ becomes. It means that a magnetized disk is thinner compared to a non-magnetized disk. According to the work by \citealt{jw}, we can define three regions. The first one is called inflow region, which starts from the disk midplane to a certain inclination $ \theta_0 $ where $ v_{r}(\theta_{0}) = 0 $. In this region, the velocity has a negative value and the accretion material moves toward the central object. The outflow region, where $ v_{r}(\theta) > 0 $, is placed between $ \theta_{0} $ and surface of the disk, $ \theta_{0} < \theta < \theta_{s} $. In this area, the accretion flow moves away from the central object. The third region, which is located between the surface of the disk and the polar axis, is called wind region., Comment: To appear in MNRAS

We study the dynamics of domain walls in a double-field model in which the U(1) symmetry is broken both spontaneously and explicitly. The global U(1) symmetry of the system is restored when the symmetry breaking parameter $\epsilon$ is set to zero. Two pairs of degenerate kinks exist in the model with are related to each other by a $Z_2$ transformation. We first calculate the single domain wall solutions and then investigate collision processes. These include simple scattering, pair annihilation, pair capture, and other interesting processes. The possibility of the domain wall being punctured by a string is also investigated., Comment: 15 pages, 7 figures

Given the fact that accretion discs are associated with their parent molecular cloud, we studied its effects as a constraint on the outer boundary of the viscous-resistive polytropic self-gravitating accretion flows subject to the mass and angular momentum loss., Comment: Proceedings of IAUS 290 "Feeding Compact Objects: Accretion on All Scales", C. M. Zhang, T. Belloni, M. Mendez & S. N. Zhang (eds.)

Self-similar and semi-analytical solutions are found for the height-averaged equations govern the dynamical behavior of a polytropic, self-gravitating disk under the effects of winds, around the nascent object. In order to describe time evolution of the system, we adopt a radius dependent mass loss rate, then highlight its importance on both the traditional $\alpha$ and innovative $\beta$ models of viscosity prescription. In agreement with some other studies, our solutions represent that Toomre parameter is less than one in most regions on the $\beta$-disk which indicate that in such disks gravitational instabilities can occur in various distances from the central accretor and so the $\beta$-disk model might provide a good explanation of how the planetary systems form. The purpose of the present work is twofold. First, examining the structure of disk with wind in comparison to no-wind solution; and second, to see if the adopted viscosity prescription affects significantly the dynamical behavior of the disk-wind system. We also considered the temperature distribution in our disk by a polytropic condition. The solutions imply that, under our boundary conditions, the radial velocity is larger for $\alpha$-disks and increases as wind becomes stronger in both viscosity models. Also, we noticed that the disk thickness increases by amplifying the wind or adopting larger values for polytropic exponent $\gamma$. It also may globally decrease if one prescribe $\beta$-model for the viscosity. Moreover, in both viscosity models, surface density and mass accretion rate reduce as wind gets stronger or $\gamma$ increases., Comment: Accepted for publication in Astrophysical Journal

We carry out the effect of $\beta$-prescription for viscosity which introduced by Duschel et al. 2000 & Hure, Richard & Zhan 2001, in a standard self-gravitating thin disks. We were predicted in a self-gravitating thin disk the $\beta$-model will have different dynamical behavior compare the well known $\alpha$-prescriptions. We used self-similar methods for solving the integrated equations which govern the dynamical behavior of the thin disks. We present the results of self-similar solutions of the time evolution of axisymmetric, polytropic, self-gravitating viscous disks around a new born central object. We apply a $\beta$-viscosity prescription which has been derived from rotating shear flow experiments ($\nu=\beta r^2\Omega$). Using reduced equations in a slow accretion limit, we demonstrate inside-out self-similar solutions after core formation in the center. Some physical quantities for $\beta$-disks are determined numerically.We have compared our results with $\alpha$-disks under the same initial conditions. It has been found that the accretion rate onto the central object for $\beta$-disks more than $\alpha$-disks at least in the outer regions where $\beta$-disks are more efficient. Our results show that Toomre instability parameter is less than one everywhere on the $\beta$-disk which means that in such disks gravitational instabilities can be occurred, so the $\beta$-disk model can be a good candidate for the origin of planetary systems. Our results show that the $\beta$-disks will decouple in the outer part of the disk where the self-gravity plays an important role which is in agreement with Duschl predictions., Comment: 8 pages, 6 figures, accepted by A&A