Source author record

Shahram Abbassi

Shahram Abbassi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2026arXiv

Long Term (250 years) Hydrodynamical Simulation of the Supermassive Black Hole Binary OJ287

With upcoming facilities capable of detecting photometric and gravitational wave signals from supermassive black hole (SMBH) binaries, studying their long-term accretion-driven variability is timely. OJ287 is a bright, nearby ($z=0.3$), and well-studied candidate for a SMBH binary. As such, it is an excellent case study for how binary dynamics could influence observed active galactic nucleus (AGN) photometric variability. We present 3D hydrodynamic simulations of OJ287, using the code PHANTOM. We simulate two mass ratios: (i) M$_1$ $=$ 1.835$\times$10$^{10}$ M$_\odot$ with M$_2$ $=$ 1.4$\times$10$^{8}$ M$_\odot$, (ii) M$_1\approx$ M$_2$ ($\sim10^{8}$ M$_\odot$) along and (iii) control of a single SMBH and accretion disc. We find that the simulation with masses 1.835$\times$10$^{10}$ M$_\odot$ and 1.4$\times$10$^{8}$ M$_\odot$ evolves consistently with the most currently accepted model of OJ287 as a precessing SMBH binary. The secondary's impacts with the disc result in the formation of spiral density waves and a corresponding $\sim$10-20% increases in the mass accretion rate of the primary SMBH. The impact timings and the mass accretion rate spikes show quasi-periodic variability as a result of the precession of the secondary's orbit with intervals between impacts ranging from $\sim$ 1 year to $\sim$ 10 years. In the near-equal mass case, the disc of the primary becomes tidally disrupted after $\sim$ 2 years. Consequently, the near-equal mass system with a period of 12 years is not a viable candidate for OJ287. This modeling provides insights into the potential signatures of SMBH binaries by both gravitational wave observatories and the Rubin Legacy Survey of Space and Time.

preprint2022arXiv

Global structure and dynamics of slowly rotating accretion flows

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 $α$. Finally, we study the effect of galaxy potential on slowly rotating accretion flows.

preprint2021arXiv

Gravitational instability of non-isothermal filamentary molecular clouds, in presence of external pressure

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.

preprint2020arXiv

Analytical Solutions of Radiative Transfer Equations in Accretion Discs with Finite Optical Depth

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.

preprint2020arXiv

Does modified gravity predict fast stellar bars in spiral galaxies?

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.

preprint2019arXiv

Thermal instability of thin accretion disks in the presence of wind and toroidal magnetic field

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.

preprint2019arXiv

Viscous Evolution of Magnetized Clumps: a Source for X-ray Flares in Gamma-ray Bursts

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.

preprint2017arXiv

Self-Gravity in Magnetized Neutrino-Dominated Accretion Discs

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.

preprint2016arXiv

Exact analytical solutions for ADAFs

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φ}$. Furthermore we assume that the value of viscosity coefficient $α$ varies with $θ$. 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.

preprint2016arXiv

Gravitational instability of filamentary molecular clouds, including ambipolar diffusion

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.

preprint2016arXiv

On the gravitational stability of the Maclaurin disk

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 $α$ and $μ_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 $(α,μ_0)$ have destabilizing effects on the disk. Interestingly $μ_0$ excites only the bar mode $m=2$ while $α$ affects all the modes. More specifically, when $α>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.

preprint2016arXiv

The influence of outflow in supercritical accretion flows

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.

preprint2015arXiv

Hydrodynamical wind on vertically self-gravitating ADAFs in the presence of toroidal magnetic field

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 $α-$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.

preprint2015arXiv

On the stability of a galactic disk in modified gravity

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 $Σ_0<Σ_{\text{crit}}(v_s,κ,α,μ_0)$, where $Σ_{\text{crit}}$ is a function of $v_s$ (sound speed), $κ$ (epicycle frequency) and $α$ and $μ_0$ are the free parameters of the theory. In the case of a stellar disk the radial velocity dispersion $σ_r$ appears in $Σ_{\text{crit}}$ instead of $v_s$. We find the exact form of the function $Σ_{\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.

preprint2015arXiv

The Effect of Large Scale Magnetic Field on Outflow in ADAFs: an Odd Symmetry Configuration

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_θ$, and even symmetry in $B_r$ and $B_ϕ$. Following previous authors, we also assume that the polar velocity $v_θ$ is nonzero. We focus on four parameters: $β_{r0}$, $β_{\phi0}$ (the plasma beta parameters for associated with magnetic field components at the equatorial plane), the magnetic resistivity $η_0$, and the density index $n=-d\lnρ/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 $β_{r0}$ ($\le10^{-2}$ for $n>1$) and larger values of $β_{\phi0}$, $η_0$ and $n$.

preprint2015arXiv

Vertical structure of Advection dominated Accretion Flows

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, θ, ϕ) $. 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 $ τ_{r ϕ} $ component of viscous stress tensor is assumed and we have set $ v_θ = 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 $ θ$ 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.

preprint2014arXiv

Dynamics of clumps embedded in a hot accretion flow with toroidal magnetic field

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.

preprint2014arXiv

Jeans analysis in modified gravity

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.

preprint2013arXiv

Structure of Advection-Dominated Accretion Disks with Outflows: Role of Toroidal Magnetic Field

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, θ, φ) $ 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_θ, v_φ) $ are present in our work. We reduce the equations to a set of differential equations about $ θ$ 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 $ θ_{0} $, always exist. The results represent that the stronger the magnetic field, the smaller the inclination angle, $ θ_{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 $ θ_0 $ where $ v_{r}(θ_{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}(θ) > 0 $, is placed between $ θ_{0} $ and surface of the disk, $ θ_{0} < θ< θ_{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.

preprint2013arXiv

The Effect of Toroidal Magnetic Field on Thickness of a Viscose-Resistive Hot Accreting Flow

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, $Δθ$, 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, $β_0$, showing the magnetic field strength at the equatorial plane, decreases $Δθ$. The other one, $η_0$ is the magnetic resistivity parameter and when it increases, $Δθ$ increases, too.

preprint2012arXiv

Viscous accretion of a polytropic self-gravitating disk in the presence of wind

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 $α$ and innovative $β$ 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 $β$-disk which indicate that in such disks gravitational instabilities can occur in various distances from the central accretor and so the $β$-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 $α$-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 $γ$. It also may globally decrease if one prescribe $β$-model for the viscosity. Moreover, in both viscosity models, surface density and mass accretion rate reduce as wind gets stronger or $γ$ increases.

preprint2006arXiv

Self-Similar Evolutionary Solutions of Self-Gravitating, Polytropic $β$-Viscous Disks

We carry out the effect of $β$-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 $β$-model will have different dynamical behavior compare the well known $α$-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 $β$-viscosity prescription which has been derived from rotating shear flow experiments ($ν=βr^2Ω$). 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 $β$-disks are determined numerically.We have compared our results with $α$-disks under the same initial conditions. It has been found that the accretion rate onto the central object for $β$-disks more than $α$-disks at least in the outer regions where $β$-disks are more efficient. Our results show that Toomre instability parameter is less than one everywhere on the $β$-disk which means that in such disks gravitational instabilities can be occurred, so the $β$-disk model can be a good candidate for the origin of planetary systems. Our results show that the $β$-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.