Source author record

Martin E. Pessah

Martin E. Pessah 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

18works
6topics
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

18 published item(s)

preprint2022arXiv

Finite-time response of dynamo mean-field effects in magnetorotational turbulence

Accretion disc turbulence along with its effect on large-scale magnetic fields plays an important role in understanding disc evolution in general, and the launching of astrophysical jets in particular. Motivated by enabling a comprehensive sub-grid description for global long-term simulations of accretions discs, we aim to further characterize the transport coefficients emerging in local simulations of magnetorotational disc turbulence. For the current investigation, we leverage a time-dependent version of the test-field method, which is sensitive to the turbulent electromotive force (EMF) generated as a response to a set of pulsating background fields. We obtain Fourier spectra of the transport coefficients as a function of oscillation frequency. These are well approximated by a simple response function, describing a finite-time build-up of the EMF as a result of a time-variable mean magnetic field. For intermediate timescales (i.e., slightly above the orbital frequency), we observe a significant phase lag of the EMF compared to the causing field. Augmented with our previous result on a non-local closure relation in space, and incorporated into a suitable mean-field description that we briefly sketch out here, the new framework will allow to drop the restrictive assumption of scale separation.

preprint2022arXiv

From Pericenter and Back: Full Debris Stream Evolution in Tidal Disruption Events

When a star passes too close to a supermassive black hole, it gets disrupted by strong tidal forces. The stellar debris then evolves into an elongated stream of gas that partly falls back towards the black hole. We present an analytical model describing for the first time the full stream evolution during such a tidal disruption event (TDE). Our framework consists in dividing the stream into different sections of elliptical geometry, whose properties are independently evolved in their co-moving frame under the tidal, pressure, and self-gravity forces. Through an explicit treatment of the tidal force and the inclusion of the gas angular momentum, we can accurately follow the stream evolution near pericenter. Our model evolves the longitudinal stream stretching and both transverse widths simultaneously. For the latter, we identify two regimes depending on whether the dynamics is entirely dominated by the tidal force (ballistic regime) or additionally influenced by pressure and self-gravity (hydrostatic regime). We find that the stream undergoes transverse collapses both shortly after the stellar disruption and upon its return near the black hole, at specific locations determined by the regime of evolution considered. The stream evolution predicted by our model can be used to determine the subsequent interactions experienced by this gas that are at the origin of most of the electromagnetic emission from TDEs. Our results suggest that the accretion disk may be fed at a rate that differs from the standard fallback rate, which would provide novel observational signatures dependent on black hole spin.

preprint2020arXiv

Local Simulations of Instabilities Driven by Composition Gradients in the ICM

The distribution of Helium in the intracluster medium (ICM) permeating galaxy clusters is not well constrained due to the very high plasma temperature. Therefore, the plasma is often assumed to be homogeneous. A non-uniform Helium distribution can however lead to biases when measuring key cluster parameters. This has motivated one-dimensional models that evolve the ICM composition assuming that the effects of magnetic fields can be parameterized or ignored. Such models for non-isothermal clusters show that Helium can sediment in the cluster core leading to a peak in concentration offset from the cluster center. The resulting profiles have recently been shown to be linearly unstable when the weakly-collisional character of the magnetized plasma is considered. In this paper, we present a modified version of the MHD code Athena, which makes it possible to evolve a weakly-collisional plasma subject to a gravitational field and stratified in both temperature and composition. We thoroughly test our implementation and confirm excellent agreement against several analytical results. In order to isolate the effects of composition, in this initial study we focus our attention on isothermal plasmas. We show that plasma instabilities, feeding off gradients in composition, can induce turbulent mixing and saturate by re-arranging magnetic field lines and alleviating the composition gradient. Composition profiles that increase with radius lead to instabilities that saturate by driving the average magnetic field inclination to roughly $45^{\circ}$. We speculate that this effect may alleviate the core insulation observed in homogeneous settings, with potential consequences for the associated cooling flow problem.

preprint2020arXiv

On Helium Mixing in Quasi-global Simulations of the Intracluster Medium

The assumption of a spatially uniform helium distribution in the intracluster medium can lead to biases in the estimates of key cluster parameters if composition gradients are present. The helium concentration profile in galaxy clusters is unfortunately not directly observable. Current models addressing the putative sedimentation are one-dimensional and parametrize the presence of magnetic fields in a crude way, ignoring the weakly-collisional, magnetized nature of the medium. When these effects are considered, a wide variety of instabilities can play an important role in the plasma dynamics. In a series of recent papers, we have developed the local, linear theory of these instabilities and addressed their non-linear development with a modified version of Athena. Here, we extend our study by developing a quasi-global approach that we use to simulate the mixing of helium as induced by generalizations of the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI), which feed off thermal and composition gradients. In the inner region of the ICM, mixing can occur on few Gyrs, after which the average magnetic field inclination angle is $\sim 45^{\circ}$ resulting in an averaged Spitzer parameter higher by about 20 % than the value obtained in homogeneous simulations. In the cluster outskirts the instabilities are rather inefficient, due to the shallow gradients. This suggests that compositions gradients in cluster cores might be shallower than one-dimensional models predict. More quantitative statements demand more refined models that can incorporate the physics driving the sedimentation process and simultaneously account for the weakly-collisional nature of the plasma.

preprint2016arXiv

Stability Analysis of Convection in the Intracluster Medium

We use the machinery usually employed for studying the onset of Rayleigh--Bénard convection in hydro- and magnetohydro-dynamic settings to address the onset of convection induced by the magnetothermal instability and the heat-flux-buoyancy-driven-instability in the weakly-collisional magnetized plasma permeating the intracluster medium. Since most of the related numerical simulations consider the plasma being bounded between two `plates' on which boundary conditions are specified, our strategy provides a framework that could enable a more direct connection between analytical and numerical studies. We derive the conditions for the onset of these instabilities considering the effects of induced magnetic tension resulting from a finite plasma beta. We provide expressions for the Rayleigh number in terms of the wave vector associated with a given mode, which allow us to characterize the modes that are first to become unstable. For both the heat-flux-buoyancy-driven-instability and the magnetothermal instability, oscillatory marginal stable states are possible.

preprint2016arXiv

Sustained Turbulence in Differentially Rotating Magnetized Fluids at Low Magnetic Prandtl Number

We show for the first time that sustained turbulence is possible at low magnetic Prandtl number for Keplerian flows with no mean magnetic flux. Our results indicate that increasing the vertical domain size is equivalent to increasing the dynamical range between the energy injection scale and the dissipative scale. This has important implications for a large variety of differentially rotating systems with low magnetic Prandtl number such as protostellar disks and laboratory experiments.

preprint2015arXiv

Characterizing the mean-field dynamo in turbulent accretion disks

The formation and evolution of a wide class of astrophysical objects is governed by turbulent, magnetized accretion disks. Understanding their secular dynamics is of primary importance. Apart from enabling mass accretion via the transport of angular momentum, the turbulence affects the long-term evolution of the embedded magnetic flux, which in turn regulates the efficiency of the transport. In this paper, we take a comprehensive next step towards an effective mean-field model for turbulent astrophysical disks by systematically studying the key properties of magnetorotational turbulence in vertically-stratified, isothermal shearing boxes. This allows us to infer emergent properties of the ensuing chaotic flow as a function of the shear parameter as well as the amount of net-vertical flux. Using the test-field method, we furthermore characterize the mean-field dynamo coefficients that describe the long-term evolution of large-scale fields. We simultaneously infer the vertical shape and the spectral scale dependence of these closure coefficients, with the latter describing non-local contributions to the turbulent electromotive force. Based on this, we infer a scale-separation ratio of about ten for the large-scale dynamo. We finally determine scaling properties of the mean-field dynamo coefficients. The relevant component of the dynamo α effect is found to scale linearly with the shear rate, as is the corresponding turbulent diffusion, η. Together, these scalings allow us to predict, in a quantitative manner, the cycle period of the well-known butterfly diagram. This lends new support to the importance of the αΩ mechanism in determining the evolution of large-scale magnetic fields in turbulent accretion disks.

preprint2015arXiv

On the axisymmetric stability of stratified and magnetized accretion disks

We conduct a comprehensive axisymmetric, local linear mode analysis of a stratified, differentially rotating disk permeated by a toroidal magnetic field which could provide significant pressure support. In the adiabatic limit, we derive a new stability criteria that differs from the one obtained for weak magnetic fields with a poloidal component and reduces continuously to the hydrodynamic Solberg-Høiland criteria. Three fundamental unstable modes are found in the locally isothermal limit. They comprise of overstable: (i) acoustic oscillations, (ii) radial epicyclic (acoustic-inertial) oscillations and (iii) vertical epicyclic (or vertical shear) oscillations. All three modes are present for finite ranges of cooling times but they are each quickly quenched past respective cut-off times. The acoustic and acoustic-inertial overstable modes are driven by the background temperature gradient. When vertical structure is excluded, we find that the radial epicyclic modes appear as a nearly degenerate pair. One of these is the aforementioned acoustic-inertial mode and the other has been previously identified in a slightly different guise as the convective overstability. Inclusion of vertical structure leads to the development of overstable oscillations destabilized by vertical shear but also has the effect of suppressing the radial epicyclic modes. Although our study does not explicitly account for non-ideal effects, we argue that it may still shed light into the dynamics of protoplanetary disk regions where a strong toroidal field generates as a result of the Hall-shear instability.

preprint2015arXiv

On Vertically Global, Horizontally Local Models for Astrophysical Disks

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

preprint2015arXiv

Plasma Instabilities in the Context of Current Helium Sedimentation Models: Dynamical Implications for the ICM in Galaxy Clusters

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable, to different kinds of instabilities depending on the magnetic field orientation, at all radii. The fastest growing modes are usually related to generalizations of the Magnetothermal Instability (MTI) and the Heat-flux-driven Buoyancy Instability (HBI) which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow fast compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

preprint2014arXiv

On the Anisotropic Nature of MRI-Driven Turbulence in Astrophysical Disks

The magnetorotational instability is thought to play an important role in enabling accretion in sufficiently ionized astrophysical disks. The rate at which MRI-driven turbulence transports angular momentum is related to both the strength of the amplitudes of the fluctuations on various scales and the degree of anisotropy of the underlying turbulence. This has motivated several studies of the distribution of turbulent power in spectral space. In this paper, we investigate the anisotropic nature of MRI-driven turbulence using a pseudo-spectral code and introduce novel ways to robustly characterize the underlying turbulence. We show that the general flow properties vary in a quasi-periodic way on timescales comparable to 10 inverse angular frequencies motivating the temporal analysis of its anisotropy. We introduce a 3D tensor invariant analysis to quantify and classify the evolution of the anisotropic turbulent flow. This analysis shows a continuous high level of anisotropy, with brief sporadic transitions towards two- and three-component isotropic turbulent flow. This temporal-dependent anisotropy renders standard shell-average, especially when used simultaneously with long temporal averages, inadequate for characterizing MRI-driven turbulence. We propose an alternative way to extract spectral information from the turbulent magnetized flow, whose anisotropic character depends strongly on time. This consists of stacking 1D Fourier spectra along three orthogonal directions that exhibit maximum anisotropy in Fourier space. The resulting averaged spectra show that the power along each of the three independent directions differs by several orders of magnitude over most scales, except the largest ones. Our results suggest that a first-principles theory to describe fully developed MRI-driven turbulence will likely have to consider the anisotropic nature of the flow at a fundamental level.

preprint2013arXiv

The Stability of Weakly Collisional Plasmas with Thermal and Composition Gradients

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport phenomena, and long-term evolution of weakly-collisional, magnetized plasmas which are stratified in temperature. These studies have improved our understanding of the physics governing the intra-cluster medium (ICM), but assumed that ICM is a homogeneous. This, however, might not be a good approximation if heavy elements sediment in the inner region of the galaxy cluster. In this paper, we analyze the stability of a weakly-collisional, magnetized plane-parallel atmosphere which is stratified in both temperature and composition. This allows us to discuss for the first time the dynamics of weakly-collisional environments where heat conduction, momentum transport, and ion-diffusion are anisotropic with respect to the direction of the magnetic field. We show that, depending on the relative signs and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous media. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could be alleviated if the mean molecular weight gradient is steep enough, i.e., $(\nabla μ)/μ> (\nabla T)/T$. This study constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of heavy elements in the ICM. (abridged)

preprint2012arXiv

Angular Momentum Transport in Protoplanetary and Black-Hole Accretion Disks: The Role of Parasitic Modes in the Saturation of MHD Turbulence

The magnetorotational instability (MRI) is considered a key process for driving efficient angular momentum transport in astrophysical disks. Understanding its non-linear saturation constitutes a fundamental problem in modern accretion disk theory. The large dynamical range in physical conditions in accretion disks makes it challenging to address this problem only with numerical simulations. We analyze the concept that (secondary) parasitic instabilities are responsible for the saturation of the MRI. Our approach enables us to explore dissipative regimes that are relevant to astrophysical and laboratory conditions that lie beyond the regime accessible to current numerical simulations. We calculate the spectrum and physical structure of parasitic modes that feed off the fastest, exact (primary) MRI mode when its amplitude is such that the fastest parasitic mode grows as fast as the MRI. We argue that this "saturation" amplitude provides an estimate of the magnetic field that can be generated by the MRI before the secondary instabilities suppress its growth significantly. Recent works suggest that the saturation amplitude of the MRI depends mainly on the magnetic Prandtl number. Our results suggest that, as long as viscous effects do not dominate the fluid dynamics, the saturation level of the MRI depends only on the Elsasser number $Λ_η$. We calculate the ratio between the stress and the magnetic energy density, $α_{\rm sat}β_{\rm sat}$, associated with the primary MRI mode. We find that for $Λ_η>1$ Kelvin-Helmholtz modes are responsible for saturation and $α_{\rm sat}β_{\rm sat} = 0.4$, while for $Λ_η< 1$ tearing modes prevail and $α_{\rm sat}β_{\rm sat} \simeq 0.5 \, Λ_η$. Several features of MRI simulations in accretion disks surrounding young stars and compact objects can be interpreted in terms of our findings.

preprint2012arXiv

MHD Turbulence in Accretion Disk Boundary Layers

The physical modeling of the accretion disk boundary layer, the region where the disk meets the surface of the accreting star, usually relies on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity, widely adopted in astrophysics, satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability is inefficient in this inner disk region. I will discuss the results of a recent study on the generation of hydromagnetic stresses and energy density in the boundary layer around a weakly magnetized star. Our findings suggest that although magnetic energy density can be significantly amplified in this region, angular momentum transport is rather inefficient. This seems consistent with the results obtained in numerical simulations and suggests that the detailed structure of turbulent MHD boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

preprint2012arXiv

On Hydromagnetic Stresses in Accretion Disk Boundary Layers

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light into physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

preprint2011arXiv

The Stability of Dilute Plasmas with Thermal and Composition Gradients. I. The Slow Conduction Limit: Overstable Gravity Modes

We analyze the stability of a dilute plasma with thermal and composition gradients in the limit where conduction is slow compared to the dynamical timescale. We find necessary and sufficient conditions for stability when the background magnetic field is either parallel or perpendicular to the thermal and composition gradients that are parallel to the gravitational field. We provide approximate solutions for all the relevant modes involved, which are driven by gravity, conduction, and diffusion. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented.

preprint2009arXiv

Coronae as Consequence of Large Scale Magnetic Fields in Turbulent Accretion Disks

Non-thermal X-ray emission in compact accretion engines can be interpreted to result from magnetic dissipation in an optically thin magnetized corona above an optically thick accretion disk. If coronal magnetic field originates in the disk and the disk is turbulent, then only magnetic structures large enough for their turbulent shredding time to exceed their buoyant rise time survive the journey to the corona. We use this concept and a physical model to constrain the minimum fraction of magnetic energy above the critical scale for buoyancy as a function of the observed coronal to bolometric emission. Our results suggest that a significant fraction of the magnetic energy in accretion disks resides in large scale fields, which in turn provides circumstantial evidence for significant non-local transport phenomena and the need for large scale magnetic field generation. For the example of Seyfert AGN, for which of order 30 per cent of the bolometric flux is in the X-ray band, we find that more than 20 per cent of the magnetic energy must be of large enough scale to rise and dissipate in the corona.

preprint2009arXiv

On the Saturation of the Magnetorotational Instability via Parasitic Modes

We investigate the stability of incompressible, exact, non-ideal magnetorotational (MRI) modes against parasitic instabilities. Both Kelvin-Helmholtz and tearing-mode parasitic instabilities may occur in the dissipative regimes accessible to current numerical simulations. We suppose that a primary MRI mode saturates at an amplitude such that its fastest parasite has a growth rate comparable to its own. The predicted alpha parameter then depends critically on whether the fastest primary and parasitic modes fit within the computational domain and whether non-axisymmetric parasitic modes are allowed. Hence even simulations that resolve viscous and resistive scales may not saturate properly unless the numerical domain is large enough to allow the free evolution of both MRI and parasitic modes. To minimally satisfy these requirements in simulations with vertical background fields, the vertical extent of the domain should accommodate the fastest growing MRI mode while the radial and azimuthal extents must be twice as large. The fastest parasites have horizontal wavelengths roughly twice as long as the vertical wavelength of the primary.