Showing 1–50 of 52 results for author: Pessah, M E

Gaseous Dynamical Friction on Elliptical Keplerian Orbits

Authors: David O'Neill, Daniel J. D'Orazio, Johan Samsing, Martin E. Pessah

Abstract: We compute the Gaseous Dynamical Friction (GDF) force experienced by massive perturbers on elliptical Keplerian orbits. In this paper, we investigate the density wake morphology, dynamical friction force, and secular orbital evolution for massive single perturbers as well as equal mass binaries embedded in an homogenous, static background flow. In all cases, the rate-of-change in semi-major axis i… ▽ More We compute the Gaseous Dynamical Friction (GDF) force experienced by massive perturbers on elliptical Keplerian orbits. In this paper, we investigate the density wake morphology, dynamical friction force, and secular orbital evolution for massive single perturbers as well as equal mass binaries embedded in an homogenous, static background flow. In all cases, the rate-of-change in semi-major axis is found to be negative (as expected), whereas the rate-of-change in eccentricity is negative for strictly-subsonic trajectories and positive for strictly-supersonic trajectories. Transonic orbits can experience both positive and negative torques during the course of an orbit, with some growing in eccentricity and others circularising. We observe all initial orbits becoming highly supersonic and eccentric (over sufficiently long timescales) due to a relentless semi-major axis decay increasing the Mach number and subsequent eccentricity driving. We compare our findings to previous studies for rectilinear and circular motion, while also making our data for orbital decay available. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to ApJ. Comments are welcome

MRI turbulence in vertically stratified accretion discs at large magnetic Prandtl numbers

Authors: Loren E. Held, George Mamatsashvili, Martin E. Pessah

Abstract: The discovery of the first binary neutron star merger, GW170817, has spawned a plethora of global numerical relativity simulations. These simulations are often ideal (with dissipation determined by the grid) and/or axisymmetric (invoking ad hoc mean-field dynamos). However, binary neutron star mergers (similar to X-ray binaries and active galactic nuclei inner discs) are characterised by large mag… ▽ More The discovery of the first binary neutron star merger, GW170817, has spawned a plethora of global numerical relativity simulations. These simulations are often ideal (with dissipation determined by the grid) and/or axisymmetric (invoking ad hoc mean-field dynamos). However, binary neutron star mergers (similar to X-ray binaries and active galactic nuclei inner discs) are characterised by large magnetic Prandtl numbers, $\rm Pm$, (the ratio of viscosity to resistivity). $\rm Pm$ is a key parameter determining dynamo action and dissipation but it is ill-defined (and likely of order unity) in ideal simulations. To bridge this gap, we investigate the magnetorotational instability (MRI) and associated dynamo at large magnetic Prandtl numbers using fully compressible, three-dimensional, vertically stratified, isothermal simulations of a local patch of a disc. We find that, within the bulk of the disc ($z\lesssim2H$, where $H$ is the scale-height), the turbulent intensity (parameterized by the stress-to-thermal-pressure ratio $α$), and the saturated magnetic field energy density, $E_\text{mag}$, produced by the MRI dynamo, both scale as a power with Pm at moderate Pm ($4\lesssim \text{Pm} \lesssim 32$): $E_\text{mag} \sim \text{Pm}^{0.74}$ and $α\sim \text{Pm}^{0.71}$, respectively. At larger Pm ($\gtrsim 32$) we find deviations from power-law scaling and the onset of a plateau. Compared to our recent unstratified study, this scaling with Pm becomes weaker further away from the disc mid-plane, where the Parker instability dominates. We perform a thorough spectral analysis to understand the underlying dynamics of small-scale MRI-driven turbulence in the mid-plane and of large-scale Parker-unstable structures in the atmosphere. △ Less

Submitted 9 April, 2024; v1 submitted 30 September, 2023; originally announced October 2023.

Comments: Accepted by MNRAS (17 pages + 2 pages appendices, 14 figures, 3 tables)

Quantifying the Impact of the Dust Torque on the Migration of Low-mass Planets

Authors: Octavio M. Guilera, Pablo Benitez-Llambay, Marcelo M. Miller Bertolami, Martin E. Pessah

Abstract: Disk solids are critical in many planet formation processes, however, their effect on planet migration remains largely unexplored. Here we assess for the first time this important issue by building on the systematic measurements of dust torques on an embedded planet by Benitez-Llambay & Pessah (2018). Adopting standard models for the gaseous disk and its solid content, we quantify the impact of th… ▽ More Disk solids are critical in many planet formation processes, however, their effect on planet migration remains largely unexplored. Here we assess for the first time this important issue by building on the systematic measurements of dust torques on an embedded planet by Benitez-Llambay & Pessah (2018). Adopting standard models for the gaseous disk and its solid content, we quantify the impact of the dust torque for a wide range of conditions describing the disk/planet system. We show that the total torque can be positive and revert inward planet migration for planetary cores with $M_{\rm p} \lesssim 10 M_\oplus$. We compute formation tracks for low-mass embryos for conditions usually invoked when modeling planet formation processes. Our most important conclusion is that dust torques can have a significant impact on the migration and formation history of planetary embryos. The most important implications of our findings are: $\it{i})$ For nominal dust-to-gas mass ratios $ε\simeq 0.01$, low-mass planets migrate outwards beyond the water ice-line if most of the mass in solids is in particles with Stokes numbers St $\simeq 0.1$. $\it{ii})$. For $ε\gtrsim 0.02-0.05$, solids with small Stokes numbers, St $\simeq 0.01$, can play a dominant role if most of the mass is in those particles. $\it{iii})$ Dust torques have the potential to enable low-mass planetary cores formed in the inner disk to migrate outwards and act as the seed for massive planets at distances of tens of au. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: Paper submitted to ApJ after minor corrections required. Feedback from the community is welcome

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

Authors: Oliver Gressel, Martin E. Pessah

Abstract: 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 simulatio… ▽ More 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. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 8 pages, 4 figures, 1 table, accepted for publication in the ApJ

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

Authors: Clément Bonnerot, Martin E. Pessah, Wenbin Lu

Abstract: 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 differ… ▽ More 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. △ Less

Submitted 17 August, 2022; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: 7 pages, 3 figures, accepted by ApJL

AGN as Potential Factories for Eccentric Black Hole Mergers

Authors: J. Samsing, I. Bartos, D. J. D'Orazio, Z. Haiman, B. Kocsis, N. W. C. Leigh, B. Liu, M. E. Pessah, H. Tagawa

Abstract: There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to th… ▽ More There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to the merger, as basic physics would argue for circularization. Here we show that AGN-disk environments can lead to an excess of eccentric mergers, if the interactions between single and binary black holes are frequent, and occur with mutual inclinations of less than a few degrees. We further illustrate that this eccentric population has a different distribution of the inclination between the spin vectors of the black holes and their orbital angular momentum at merger, referred to as the spin-orbit tilt, compared to the remaining circular mergers. △ Less

Submitted 21 April, 2022; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted

Journal ref: Nature, Volume 603, Issue 7900, p.237-240, 2022

Zero net flux MRI-turbulence in disks $-$ sustenance scheme and magnetic Prandtl number dependence

Authors: George Mamatsashvili, George Chagelishvili, Martin E. Pessah, Frank Stefani, Gianluigi Bodo

Abstract: We investigate sustenance and dependence on magnetic Prandtl number (${\rm Pm}$) for magnetorotational instability (MRI)-driven turbulence in astrophysical Keplerian disks with zero net magnetic flux using standard shearing box simulations. We focus on the turbulence dynamics in Fourier space, capturing specific/noncanonical anisotropy of nonlinear processes due to disk flow shear. This is a new t… ▽ More We investigate sustenance and dependence on magnetic Prandtl number (${\rm Pm}$) for magnetorotational instability (MRI)-driven turbulence in astrophysical Keplerian disks with zero net magnetic flux using standard shearing box simulations. We focus on the turbulence dynamics in Fourier space, capturing specific/noncanonical anisotropy of nonlinear processes due to disk flow shear. This is a new type of nonlinear redistribution of modes over wavevector orientations in Fourier space -- the nonlinear transverse cascade -- which is generic to shear flows and fundamentally different from usual direct/inverse cascade. The zero flux MRI has no exponentially growing modes, so its growth is transient, or nonmodal. Turbulence self-sustenance is governed by constructive cooperation of the transient growth of MRI and the nonlinear transverse cascade. This cooperation takes place at small wavenumbers (on the flow size scales) referred to as the vital area in Fourier space. The direct cascade transfers mode energy from the vital area to larger wavenumbers. At large ${\rm Pm}$, the transverse cascade prevails over the direct one, keeping most of modes' energy contained in small wavenumbers. With decreasing ${\rm Pm}$, however, the action of the transverse cascade weakens and can no longer oppose the action of direct cascade which more efficiently transfers energy to higher wavenumbers, leading to increased resistive dissipation. This undermines the sustenance scheme, resulting in the turbulence decay. Thus, the decay of zero net flux MRI-turbulence with decreasing ${\rm Pm}$ is attributed to topological rearrangement of the nonlinear processes when the direct cascade begins to prevail over the transverse cascade. △ Less

Submitted 7 October, 2020; v1 submitted 30 September, 2020; originally announced September 2020.

Comments: 26 pages, 16 figures, accepted for publication in the Astrophysical Journal

Suppressed heat conductivity in the intracluster medium: implications for the magneto-thermal instability

Authors: Thomas Berlok, Eliot Quataert, Martin E. Pessah, Christoph Pfrommer

Abstract: In the outskirts of the intracluster medium (ICM) in galaxy clusters, the temperature decreases with radius. Due to the weakly collisional nature of the plasma, these regions are susceptible to the magneto-thermal instability (MTI), which can sustain turbulence and provide turbulent pressure support in the ICM. This instability arises due to heat conduction directed along the magnetic field, with… ▽ More In the outskirts of the intracluster medium (ICM) in galaxy clusters, the temperature decreases with radius. Due to the weakly collisional nature of the plasma, these regions are susceptible to the magneto-thermal instability (MTI), which can sustain turbulence and provide turbulent pressure support in the ICM. This instability arises due to heat conduction directed along the magnetic field, with a heat conductivity which is normally assumed to be given by the Spitzer value. Recent numerical studies of the ion mirror and the electron whistler instability using particle-in-cell codes have shown that microscale instabilities can lead to a reduced value for the heat conductivity in the ICM. This could in turn influence the efficiency with which the MTI drives turbulence. In this paper we investigate the influence of reduced heat transport on the nonlinear evolution of the MTI. We study plane-parallel, initially static atmospheres and employ a subgrid model that mimics the influence of the mirror instability on the heat conductivity. We use this subgrid model to assess the effect of microscales on the large scale dynamics of the ICM. We find that the nonlinear saturation of the MTI is surprisingly robust in our simulations. Over a factor of $\sim 10^3$ in the thermal-to-magnetic pressure ratio and collisionality we find at most modest changes to the saturation of the MTI with respect to reference simulations where heat transport is unsuppressed. △ Less

Submitted 18 March, 2021; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 20 pages, 18 figures, accepted for publication in MNRAS

Transport properties of Keplerian flows in extended local domains with no imposed field

Authors: Farrukh Nauman, Martin E. Pessah

Abstract: We compare transport statistics of elongated incompressible shearing boxes for different Reynolds and magnetic Prandtl numbers, $Re$ and $Pm$, and aspect ratios, $L_z/L_x$. We find that at fixed aspect ratio $L_z/L_x=4$ and $Re = 10,000$, the turbulent stresses for $Pm \lesssim 1$ do not show considerable variation and follow a power law $\sim Pm^{3/2}$ for $Pm > 1$. This is qualitatively consiste… ▽ More We compare transport statistics of elongated incompressible shearing boxes for different Reynolds and magnetic Prandtl numbers, $Re$ and $Pm$, and aspect ratios, $L_z/L_x$. We find that at fixed aspect ratio $L_z/L_x=4$ and $Re = 10,000$, the turbulent stresses for $Pm \lesssim 1$ do not show considerable variation and follow a power law $\sim Pm^{3/2}$ for $Pm > 1$. This is qualitatively consistent with previous results based on net imposed flux and small box $L_z/L_x \sim 1$ simulations but the power law exponent is different. The saturated level of stresses, the ratio of Maxwell stress to the magnetic energy and Reynolds to Maxwell stress ratio are roughly invariant as $L_z/L_x$ is increased. For cases where the boxes are elongated in both the azimuth and vertical direction, the transport coefficient $α\in [0.1,1.0]$ that is $10-100$ times larger than the case with $L_y/L_x = 2$ and large $L_z/L_x$. Overall, our results suggest that the saturated state of turbulence is sensitive to both dissipation coefficients and aspect ratio (both $L_z/L_x$, $L_y/L_x$) motivating further work on this problem. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: 6 pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2018, Volume 480, Issue 1, p.204-209

Streaming Instability for Particle-Size Distributions

Authors: Leonardo Krapp, Pablo Benítez-Llambay, Oliver Gressel, Martin E. Pessah

Abstract: The streaming instability is thought to play a central role in the early stages of planet formation by enabling the efficient bypass of a number of barriers hindering the formation of planetesimals. We present the first study exploring the efficiency of the linear streaming instability when a particle-size distribution is considered. We find that, for a given dust-to-gas mass ratio, the multi-spec… ▽ More The streaming instability is thought to play a central role in the early stages of planet formation by enabling the efficient bypass of a number of barriers hindering the formation of planetesimals. We present the first study exploring the efficiency of the linear streaming instability when a particle-size distribution is considered. We find that, for a given dust-to-gas mass ratio, the multi-species streaming instability grows on timescales much longer than those expected when only one dust species is involved. In particular, distributions that contain close-to-order-unity dust-to-gas mass ratios lead to unstable modes that can grow on timescales comparable, or larger, with those of secular instabilities. We anticipate that processes leading to particle segregation and/or concentration can create favourable conditions for the instability to grow fast. Our findings may have important implications for a large number of processes in protoplanetary disks that rely on the streaming instability as usually envisioned for a unique dust species. Our results suggest that the growth rates of other resonant-drag-instabilities may also decrease considerably when multiple species are considered. △ Less

Submitted 10 June, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: 8 pages, 6 figures; accepted for publication in The Astrophysical Journal Letters

Probing the Protosolar Disk Using Dust Filtering at Gaps in the Early Solar System

Authors: Troels Haugbølle, Philipp Weber, Daniel P. Wielandt, Pablo Benítez-Llambay, Martin Bizzarro, Oliver Gressel, Martin E. Pessah

Abstract: Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interd… ▽ More Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interdisciplinary approach, we use a $μ$-X-Ray-fluorescence scanner to search for large CAIs and a scanning electron microscope to search for small CAIs in the ordinary chondrite NWA 5697. We carry out long-term, two-dimensional simulations including gas, dust, and planets to characterize the transport of grains within the viscous $α$-disk framework exploring the scenarios of a stand-alone Jupiter, Jupiter and Saturn \textit{in situ}, or Jupiter and Saturn in a 3:2 resonance. In each case, we find a critical grain size above which drift is halted as a function of the physical conditions in the disk. From the laboratory search we find four CAIs with a largest size of $\approx$200$\,μ$m. \Combining models and data, we provide an estimate for the upper limit of the $α$-viscosity and the surface density at the location of Jupiter, using reasonable assumptions about the stellar accretion rate during inward transport of CAIs, and assuming angular momentum transport to happen exclusively through viscous effects. Moreover, we find that the compound gap structure in the presence of Saturn in a 3:2 resonance favors inward transport of grains larger than CAIs currently detected in ordinary chondrites. △ Less

Submitted 17 July, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: 16 pages, 10 figures, updated to match published version in Astrophysical Journal

Journal ref: The Astronomical Journal, Volume 158, Number 2, 2019

Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D

Authors: Pablo Benítez-Llambay, Leonardo Krapp, Martín E. Pessah

Abstract: We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our… ▽ More We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our implementation, we develop a test suite for an arbitrary number of species, based on analytical or exact solutions of problems related to perfect damping, damped sound waves, shocks, local and global gas-dust radial drift in a disk and, linear streaming instability. In particular, we obtain first-order, steady-state solutions for the radial drift of multiple dust species in protoplanetary disks, in which the pressure gradient is not necessarily small. We additionally present non-linear shearing-box simulations of the streaming instability and compare them with previous results obtained with Lagrangian particles. We successfully validate our implementation by recovering the solutions from the test suite to second- and first-order accuracy in space and time, respectively. From this, we conclude that our scheme is suitable, and very robust, to study the self-consistent dynamics of several fluids. In particular, it can be used for solving the collisions between gas and dust in protoplanetary disks, with any degree of coupling. △ Less

Submitted 10 June, 2019; v1 submitted 19 November, 2018; originally announced November 2018.

Comments: 27 pages, 12 figures; accepted for publication in ApJS

Dust segregation in Hall-dominated turbulent protoplanetary disks

Authors: Leonardo Krapp, Oliver Gressel, Pablo Benítez-Llambay, Turlough P. Downes, Gopakumar Mohandas, Martin E. Pessah

Abstract: Imaging of the dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which c… ▽ More Imaging of the dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which can lead to local dust enhancements. Adopting a radially varying cylindrical disk model, and considering combinations of vertical and azimuthal initial net flux, we perform 3D non-ideal MHD simulations aimed at studying self-organization induced by the Hall effect in turbulent PPDs. To this end, new modules have been incorporated into the NIRVANA-III and FARGO3D MHD codes. We moreover include dust grains, treated in the fluid approximation, in order to study their evolution subject to the emerging zonal flows. In the regime of a dominant Hall effect, we robustly obtain large-scale organized concentrations in the vertical magnetic field that remain stable for hundreds of orbits. For disks with vertical initial net flux alone, we confirm the presence of zonal flows and vortices that introduce regions of super-Keplerian gas flow. Including a moderately strong net-azimuthal magnetic flux can significantly alter the dynamics, partially preventing the self-organization of zonal flows. For plasma beta-parameters smaller than 50, large-scale, near-axisymmetric structures develop in the vertical magnetic flux. In all cases, we demonstrate that the emerging features are capable of accumulating dust grains for a range of Stokes numbers. △ Less

Submitted 30 May, 2019; v1 submitted 23 August, 2018; originally announced August 2018.

Comments: 20 pages, 14 figures. Accepted for publication in The Astrophysical Journal

On the dynamics of pebbles in protoplanetary disks with magnetically-driven winds

Authors: Mohsen Shadmehri, Fazeleh Khajenabi, Martin E. Pessah

Abstract: We present an analytical model to investigate the production of pebbles and their radial transport through a protoplanetary disk (PPD) with magnetically driven winds. While most of the previous analytical studies in this context assume that the radial turbulent coefficient is equal to the vertical dust diffusion coefficient, in the light of the results of recent numerical simulations, we relax thi… ▽ More We present an analytical model to investigate the production of pebbles and their radial transport through a protoplanetary disk (PPD) with magnetically driven winds. While most of the previous analytical studies in this context assume that the radial turbulent coefficient is equal to the vertical dust diffusion coefficient, in the light of the results of recent numerical simulations, we relax this assumption by adopting effective parametrisations of the turbulent coefficients involved in terms of the strength of the magnetic fields driving the wind. Theoretical studies have already pointed out that even in the absence of winds, these coefficients are not necessarily equal, though its consequences regarding pebble production have not been explored. In this paper, we investigate the evolution of the pebble production line, the radial mass flux of the pebbles and their corresponding surface density as a function of the plasma parameter at the disk midplane. Our analysis explicitly demonstrates that the presence of magnetically-driven winds in a PPD leads to considerable reduction of the rate and duration of the pebble delivery. We show that when the wind is strong, the core growth in mass due to the pebble accretion is so slow that it is unlikely that a core could reach a pebble isolation mass during a PPD lifetime. When the mass of a core reaches this critical value, pebble accretion is halted due to core-driven perturbations in the gas. With decreasing wind strength, however, pebble accretion may, in a shorter time, increase the mass of a core to the pebble isolation mass. △ Less

Submitted 27 June, 2018; originally announced June 2018.

Comments: 15 pages, accepted for publication in ApJ

Analytical Models of Exoplanetary Atmospheres. V. Non-gray Thermal Structure with Coherent Scattering

Authors: Gopakumar Mohandas, Martin E. Pessah, Kevin Heng

Abstract: We apply the picket fence treatment to model the effects brought about by spectral lines on the thermal structure of irradiated atmospheres. The lines may be due to purely absorption processes, purely coherent scattering processes or some combination of absorption and scattering. If the lines arise as a pure absorption process, the surface layers of the atmosphere are cooler whereas this surface c… ▽ More We apply the picket fence treatment to model the effects brought about by spectral lines on the thermal structure of irradiated atmospheres. The lines may be due to purely absorption processes, purely coherent scattering processes or some combination of absorption and scattering. If the lines arise as a pure absorption process, the surface layers of the atmosphere are cooler whereas this surface cooling is completely absent if the lines are due to pure coherent isotropic scattering. The lines also lead to a warming of the deeper atmosphere. The warming of the deeper layers is, however, independent of the nature of line formation. Accounting for coherent isotropic scattering in the shortwave and longwave continuum results in anti-greenhouse cooling and greenhouse warming on an atmosphere-wide scale. The effects of coherent isotropic scattering in the line and continuum operate in tandem to determine the resulting thermal structure of the irradiated atmosphere. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: 8 pages, 5 figures, accepted for publication in ApJ

Characterizing the variable dust permeability of planet-induced gaps

Authors: Philipp Weber, Pablo Benítez-Llambay, Oliver Gressel, Leonardo Krapp, Martin E. Pessah

Abstract: Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size… ▽ More Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size for retention outside a gap as a function of particle size and for various key parameters defining the protoplanetary disk model. To this end, we perform multifluid hydrodynamical simulations in two dimensions, including different dust species, which we treat as pressureless fluids. We initialize the dust outside of the planet's orbit and study under which conditions dust grains are able to cross the gap carved by the planet. In agreement with previous work, we find that the permeability of the gap depends both on dust dynamical properties and the gas disk structure: while small dust follows the viscously accreting gas through the gap, dust grains approaching a critical size are progressively filtered out. Moreover, we introduce and compute a depletion factor that enables us to quantify the way in which higher viscosity, smaller planet mass, or a more massive disk can shift this critical size to larger values. Our results indicate that gap-opening planets may act to deplete the inner reaches of protoplanetary disks of large dust grains -- potentially limiting the accretion of solids onto forming terrestrial planets. △ Less

Submitted 25 January, 2018; v1 submitted 24 January, 2018; originally announced January 2018.

Comments: 18 pages, 14 figures, accepted for publication in ApJ

Torques Induced by Scattered Pebble-flow in Protoplanetary Disks

Authors: Pablo Benítez-Llambay, Martin E. Pessah

Abstract: Fast inward migration of planetary cores is a common problem in the current planet formation paradigm. Even though dust is ubiquitous in protoplanetary disks, its dynamical role in the migration history of planetary embryos has not been assessed. In this Letter, we show that the scattered pebble-flow induced by a low-mass planetary embryo leads to an asymmetric dust-density distribution that is ab… ▽ More Fast inward migration of planetary cores is a common problem in the current planet formation paradigm. Even though dust is ubiquitous in protoplanetary disks, its dynamical role in the migration history of planetary embryos has not been assessed. In this Letter, we show that the scattered pebble-flow induced by a low-mass planetary embryo leads to an asymmetric dust-density distribution that is able to exert a net torque. By analyzing a large suite of multifluid hydrodynamical simulations addressing the interaction between the disk and a low-mass planet on a fixed circular orbit, and neglecting dust feedback onto the gas, we identify two different regimes, gas- and gravity-dominated, where the scattered pebble-flow results in almost all cases in positive torques. We collect our measurements in a first torque map for dusty disks, which will enable the incorporation of the effect of dust dynamics on migration into population synthesis models. Depending on the dust drift speed, the dust-to-gas mass ratio/distribution and the embryo mass, the dust-induced torque has the potential to halt inward migration or even induce fast outward migration of planetary cores. We thus anticipate that dust-driven migration could play a dominant role during the formation history of planets. Because dust torques scale with disk metallicity, we propose that dust-driven outward migration may enhance the occurrence of distant giant planets in higher-metallicity systems. △ Less

Submitted 22 November, 2018; v1 submitted 24 January, 2018; originally announced January 2018.

Comments: 7 pages, 4 figures - Accepted for publication in ApJL

Journal ref: The Astrophysical Journal Letters, 855:L28, 2018

Orbital stability in static axisymmetric fields

Authors: Gopakumar Mohandas, Tobias Heinemann, Martin E. Pessah

Abstract: We investigate the stability of test-particle equilibrium orbits in axisymmetric, but otherwise arbitrary, gravitational and electromagnetic fields. We extend previous studies of this problem to include a toroidal magnetic field. We find that, even though the toroidal magnetic field does not alter the location of the circular orbits, it enters the problem as a gyroscopic force with the potential t… ▽ More We investigate the stability of test-particle equilibrium orbits in axisymmetric, but otherwise arbitrary, gravitational and electromagnetic fields. We extend previous studies of this problem to include a toroidal magnetic field. We find that, even though the toroidal magnetic field does not alter the location of the circular orbits, it enters the problem as a gyroscopic force with the potential to provide gyroscopic stability. This is in essence similar to the situation encountered in the reduced three-body problem where rotation enables stability around the local maxima of the effective potential. Nevertheless, we show that gyroscopic stabilization by a toroidal magnetic field is impossible for axisymmetric force fields in source-free regions because in this case the effective potential does not possess any local maxima. As an example of an axisymmetric force field with sources, we consider the classical problem of a rotating, aligned magnetosphere. By analyzing the dynamics of halo and equatorial particle orbits we conclude that axisymmetric toroidal fields that are antisymmetric about the equator are unable to provide gyroscopic stabilization. On the other hand, a toroidal magnetic field that does not vanish at the equator can provide gyroscopic stabilization for positively charged particles in prograde equatorial orbits. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 11 pages, 3 figures, submitted to Celestial Mechanics and Dynamical Astronomy

Structure of Protoplanetary Discs with Magnetically-driven Winds

Authors: Fazeleh Khajenabi, Mohsen Shadmehri, Martin E. Pessah, Rebecca G. Martin

Abstract: We present a new set of analytical solutions to model the steady state structure of a protoplanetary disc with a magnetically-driven wind. Our model implements a parametrization of the stresses involved and the wind launching mechanism in terms of the plasma parameter at the disc midplane, as suggested by the results of recent, local MHD simulations. When wind mass-loss is accounted for, we find t… ▽ More We present a new set of analytical solutions to model the steady state structure of a protoplanetary disc with a magnetically-driven wind. Our model implements a parametrization of the stresses involved and the wind launching mechanism in terms of the plasma parameter at the disc midplane, as suggested by the results of recent, local MHD simulations. When wind mass-loss is accounted for, we find that its rate significantly reduces the disc surface density, particularly in the inner disc region. We also find that models that include wind mass-loss lead to thinner dust layers. As an astrophysical application of our models, we address the case of HL Tau, whose disc exhibits a high accretion rate and efficient dust settling at its midplane. These two observational features are not easy to reconcile with conventional accretion disc theory, where the level of turbulence needed to explain the high accretion rate would prevent a thin dust layer. Our disc model that incorporates both mass-loss and angular momentum removal by a wind is able to account for HL Tau observational constraints concerning its high accretion rate and dust layer thinness. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS, 13 pages, 8 figures

Spectral Analysis of Non-Ideal MRI Modes: The effect of Hall diffusion

Authors: Gopakumar Mohandas, Martin E. Pessah

Abstract: The effect of magnetic field diffusion on the stability of accretion disks is a problem that has attracted considerable interest of late. In particular, the Hall effect has the potential to bring about remarkable changes in the dynamical behavior of disks that are without parallel. In this paper, we conduct a systematic examination of the linear eigenmodes in a weakly magnetized differentially rot… ▽ More The effect of magnetic field diffusion on the stability of accretion disks is a problem that has attracted considerable interest of late. In particular, the Hall effect has the potential to bring about remarkable changes in the dynamical behavior of disks that are without parallel. In this paper, we conduct a systematic examination of the linear eigenmodes in a weakly magnetized differentially rotating gas with special focus on Hall diffusion. We first develop a geometrical representation of the eigenmodes and provide a detailed quantitative description of the polarization properties of the oscillatory modes under the combined influence of the Coriolis and Hall effects. We also analyze the effects of magnetic diffusion on the structure of the unstable modes and derive analytical expressions for the kinetic and magnetic stresses and energy densities associated with the non-ideal MRI. Our analysis explicitly demonstrates that, if the dissipative effects are relatively weak, the kinetic stresses and energies make up the dominant contribution to the total stress and energy density when the equilibrium angular momentum and magnetic field vectors are anti-parallel. This is in sharp contrast to what is observed in the case of the ideal or dissipative MRI. We conduct shearing box simulations and find very good agreement with the results derived from linear analysis. As the modes in consideration are also exact solutions of the non-linear equations, the unconventional nature of the kinetic and magnetic stresses may have significant implications for the non-linear evolution in some regions of protoplanetary disks. △ Less

Submitted 16 February, 2017; originally announced February 2017.

Comments: 16 pages, 11 figures, accepted by ApJ

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

Authors: Farrukh Nauman, Martin E. Pessah

Abstract: 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 syste… ▽ More 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. △ Less

Submitted 27 September, 2016; originally announced September 2016.

Comments: 5 pages, 6 figures, submitted to ApJ, changes made in response to referee report

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

Authors: Thomas Berlok, Martin E. Pessah

Abstract: 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… ▽ More 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. △ Less

Submitted 3 June, 2020; v1 submitted 12 August, 2016; originally announced August 2016.

Comments: 15 pages, 9 figures, published in ApJ, animation of Figure 5 at http://tberlok.dk/movies/icm_quasi_global.html

Stability Analysis of Convection in the Intracluster Medium

Authors: Himanshu Gupta, Shailendra K. Rathor, Martin E. Pessah, Sagar Chakraborty

Abstract: 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… ▽ More 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. △ Less

Submitted 31 May, 2016; originally announced May 2016.

Comments: 10 pages, 1 figure

Local Simulations of Instabilities Driven by Composition Gradients in the ICM

Authors: Thomas Berlok, Martin E. Pessah

Abstract: 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 th… ▽ More 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. △ Less

Submitted 3 June, 2020; v1 submitted 9 March, 2016; originally announced March 2016.

Comments: 13 pages, 13 figures, published in ApJ, animation of Figure 7 at http://tberlok.dk/movies/icm_local.html

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

Authors: Thomas Berlok, Martin E. Pessah

Abstract: 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 collisio… ▽ More 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. △ Less

Submitted 12 October, 2015; originally announced October 2015.

Comments: Accepted by ApJ. 15 pages, 17 figures

On the axisymmetric stability of stratified and magnetized accretion disks

Authors: Gopakumar Mohandas, Martin E. Pessah

Abstract: 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… ▽ More 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. △ Less

Submitted 9 October, 2015; originally announced October 2015.

Comments: submitted to ApJ

Characterizing the mean-field dynamo in turbulent accretion disks

Authors: Oliver Gressel, Martin E. Pessah

Abstract: 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 th… ▽ More 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. △ Less

Submitted 30 July, 2015; originally announced July 2015.

Comments: 15 pages, 13 figures, accepted for publication in ApJ

Particle Control in Phase Space by Global K-Means Clustering

Authors: J. Trier Frederiksen, G. Lapenta, M. E. Pessah

Abstract: We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decreasing or increasing the entire particle population, based on k-means clustering of the data. In essence… ▽ More We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decreasing or increasing the entire particle population, based on k-means clustering of the data. In essence the procedure amounts to merging or splitting particles by statistical means, throughout the entire simulation volume in question, while minimizing a 6-dimensional total distance measure to preserve the physics. Particle merging is by far the most demanding procedure when considering conservation laws of physics; it amounts to lossy compression of particle phase space data. We demonstrate that our k-means approach conserves energy and momentum to high accuracy, even for high compression ratios, $\mathcal{R} \approx 3$ --- \emph{i.e.}, $N_{f} \lesssim 0.33N_{i}$. Interestingly, we find that an accurate particle splitting step can be performed using k-means as well; this from an argument of symmetry. The split solution, using k-means, places splitted particles optimally, to obtain maximal spanning on the phase space manifold. Implementation and testing is done using an electromagnetic PIC code, the \ppcode. Nonetheless, the k-means framework is general; it is not limited to Vlasov-Maxwell type PIC codes. We discuss advantages and drawbacks of this optimal phase space reconstruction. △ Less

Submitted 30 November, 2015; v1 submitted 15 April, 2015; originally announced April 2015.

Comments: Revision 1. Major revisions. Added discussion. 18 pages, 22 figures, submitted to Journal of Computational Physics

An assessment of the pulsar outer gap model. II: Implications for the predicted $γ$-ray spectra

Authors: Daniele Viganò, Diego F. Torres, Kouichi Hirotani, Martín E. Pessah

Abstract: One of the most important predictions of any gap model for pulsar magnetospheres is the predicted $γ$-ray spectra. In the outer gap model, the properties of the synchro-curvature radiation are sensitive to many parameters, whose realistic ranges have been studied in detail in an accompanying paper. There we demonstrated that the uncertainty in the radius of curvature, the magnetic field geometry,… ▽ More One of the most important predictions of any gap model for pulsar magnetospheres is the predicted $γ$-ray spectra. In the outer gap model, the properties of the synchro-curvature radiation are sensitive to many parameters, whose realistic ranges have been studied in detail in an accompanying paper. There we demonstrated that the uncertainty in the radius of curvature, the magnetic field geometry, and the X-ray surface flux may affect by orders of magnitude the predicted flux and spectral peak in the $γ$-ray regime. Here, we present a systematic, numerical study of the impact of the different parameters on the particle dynamics along the gap and calculate the emitted synchro-curvature radiation along the trajectory. By integrating the emitted radiation along the gap and convolving it with a parametrized particle distribution, we discuss how the comparison with the wealth of {\em Fermi}-LAT data can be used to constrain the applicability of the model. The resulting spectra show very different energy peaks, fluxes and shapes, qualitatively matching the great variety of the observed {\em Fermi}-LAT pulsars. In particular, if we see a large fraction of photons emitted from the initial part of the trajectory, we show that the spectra will be flatter at the low-energy {\it Fermi}-LAT regime (100 MeV -- 1 GeV). This provides a solution for such observed flat spectra, while still maintain synchro-curvature radiation as the origin of these photons. △ Less

Submitted 1 April, 2015; v1 submitted 3 December, 2014; originally announced December 2014.

Comments: 9 pages, 6 figures, 1 table. Accepted for pubblication in MNRAS main journal. References updated

An assessment of the pulsar outer gap model. I: Assumptions, uncertainties, and implications for the gap size and the accelerating field

Authors: Daniele Viganò, Diego F. Torres, Kouichi Hirotani, Martín E. Pessah

Abstract: The popular outer gap model of magnetospheric emission from pulsars has been widely applied to explain the properties observed in $γ$-rays. However, its quantitative predictions rely on a number of approximations and assumptions that are usually overlooked. Here we examine them, reviewing the main ingredients entering in the model, evaluating their range of uncertainties. Usually, in the quantitat… ▽ More The popular outer gap model of magnetospheric emission from pulsars has been widely applied to explain the properties observed in $γ$-rays. However, its quantitative predictions rely on a number of approximations and assumptions that are usually overlooked. Here we examine them, reviewing the main ingredients entering in the model, evaluating their range of uncertainties. Usually, in the quantitative applications of the model, key parameters like the radius of curvature and the energies of the interacting photons are taken to be a fixed, single value. Instead, here we explore their realistic ranges, and the impact of these on the consistency of the model itself. We conclude that the popular evaluation of the trans-field size of the gap as a function of period and period derivative, is unreliable and affected by a huge dispersion. Last, the exploration of the possible values for the radius of curvature, the local magnetic field and other quantities deserve more attention for quantitative applications of the outer gap model, like the calculation of $γ$-ray spectra, which is the subject of an accompanying paper. △ Less

Submitted 4 December, 2014; v1 submitted 3 December, 2014; originally announced December 2014.

Comments: 19 pages, 5 figures, 3 tables. Accepted for pubblication in MNRAS main journal. References updated

Compact formulae, dynamics and radiation of charged particles under synchro-curvature losses

Authors: Daniele Viganò, Diego F. Torres, Kouichi Hirotani, Martín E. Pessah

Abstract: We consider the fundamental problem of charged particles moving along and around a curved magnetic field line, revising the synchro-curvature radiation formulae introduced by Cheng and Zhang (1996). We provide more compact expressions to evaluate the spectrum emitted by a single particle, identifying the key parameter that controls the transition between the curvature-dominated and the synchrotron… ▽ More We consider the fundamental problem of charged particles moving along and around a curved magnetic field line, revising the synchro-curvature radiation formulae introduced by Cheng and Zhang (1996). We provide more compact expressions to evaluate the spectrum emitted by a single particle, identifying the key parameter that controls the transition between the curvature-dominated and the synchrotron-dominated regime. This parameter depends on the local radius of curvature of the magnetic field line, the gyration radius, and the pitch angle. We numerically solve the equations of motion for the emitting particle by considering self-consistently the radiative losses, and provide the radiated spectrum produced by a particle when an electric acceleration is balanced by its radiative losses, as it is assumed to happen in the outer gaps of pulsar's magnetospheres. We compute the average spectrum radiated throughout the particle trajectory finding that the slope of the spectrum before the peak depends on the location and size of the emission region. We show how this effect could then lead to a variety of synchro-curvature spectra. Our results reinforce the idea that the purely synchrotron or curvature losses are, in general, inadequate to describe the radiative reaction on the particle motion, and the spectrum of emitted photons. Finally, we discuss the applicability of these calculations to different astrophysical scenarios. △ Less

Submitted 9 December, 2014; v1 submitted 21 November, 2014; originally announced November 2014.

Comments: 9 pages, 5 figures, 2 tables. Accepted for publication in MNRAS main journal. References updated

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

Authors: Gareth C. Murphy, Martin E. Pessah

Abstract: 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 distr… ▽ More 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. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 13 pages, 13 figures, submitted to ApJ

Journal ref: ApJ, 802, 139 (2015)

On Vertically Global, Horizontally Local Models for Astrophysical Disks

Authors: Colin P. McNally, Martin E. Pessah

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2015; v1 submitted 18 June, 2014; originally announced June 2014.

Comments: 24 pages, 8 figures, updated to match published version

Journal ref: 2015 ApJ 811 121

The Stability of Weakly Collisional Plasmas with Thermal and Composition Gradients

Authors: Martin E. Pessah, Sagar Chakraborty

Abstract: 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 g… ▽ More 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) △ Less

Submitted 14 February, 2013; originally announced February 2013.

Comments: 16 pages, 7 figures. This article supersedes arXiv:1111.3372 (5 pages, 3 figures). The present version of this article is published in ApJ

Journal ref: ApJ, 764 (2013) 13-29

MHD Turbulence in Accretion Disk Boundary Layers

Authors: Chi-kwan Chan, Martin E. Pessah

Abstract: 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 beh… ▽ More 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. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 4 pages, 2 figures, proceedings IAU symposium 294: Solar and Astrophysical Dynamos and Magnetic Activity

Report number: NORDITA-2012-95

On Hydromagnetic Stresses in Accretion Disk Boundary Layers

Authors: Martin E. Pessah, Chi-kwan Chan

Abstract: 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 s… ▽ More 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. △ Less

Submitted 23 May, 2012; v1 submitted 17 November, 2011; originally announced November 2011.

Comments: 6 pages, 4 figures, uses emulateapj, updated to match published version in ApJ

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

Authors: Martin E. Pessah, Sagar Chakraborty

Abstract: 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 sol… ▽ More 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. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 5 pages, 3 figures. Submitted to ApJ

Sustained Magnetorotational Turbulence in Local Simulations of Stratified Disks with Zero Net Magnetic Flux

Authors: S. W. Davis, J. M. Stone, M. E. Pessah

Abstract: We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, cont… ▽ More We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles. △ Less

Submitted 8 September, 2009; originally announced September 2009.

Comments: 13 pages, 18 figures, submitted to ApJ

Journal ref: Astrophys.J.713:52-65,2010

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

Authors: Martin E. Pessah

Abstract: 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 ana… ▽ More 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. △ Less

Submitted 23 May, 2012; v1 submitted 13 August, 2009; originally announced August 2009.

Comments: 11 figures, 16 pages, uses emulateapj, updated to match published version in ApJ

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

Authors: Eric G. Blackman, Martin E. Pessah

Abstract: 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… ▽ More 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. △ Less

Submitted 13 July, 2009; originally announced July 2009.

Comments: submitted to ApJL, 2 figs

Journal ref: Astrophys.J.704:L113-L117,2009

On the Saturation of the Magnetorotational Instability via Parasitic Modes

Authors: Martin E. Pessah, Jeremy Goodman

Abstract: 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 i… ▽ More 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. △ Less

Submitted 19 May, 2009; v1 submitted 4 February, 2009; originally announced February 2009.

Comments: 5 pages, 2 figures, uses emulateapj. Modified text in order to address referee's comments and suggestions. References added. ApJ Letters accepted

Journal ref: Astrophys.J.698:L72-L76,2009

Viscous, Resistive Magnetorotational Modes

Authors: Martin E. Pessah, Chi-kwan Chan

Abstract: We carry out a comprehensive analysis of the behavior of the magnetorotational instability (MRI) in viscous, resistive plasmas. We find exact, non-linear solutions of the non-ideal magnetohydrodynamic (MHD) equations describing the local dynamics of an incompressible, differentially rotating background threaded by a vertical magnetic field when disturbances with wavenumbers perpendicular to the… ▽ More We carry out a comprehensive analysis of the behavior of the magnetorotational instability (MRI) in viscous, resistive plasmas. We find exact, non-linear solutions of the non-ideal magnetohydrodynamic (MHD) equations describing the local dynamics of an incompressible, differentially rotating background threaded by a vertical magnetic field when disturbances with wavenumbers perpendicular to the shear are considered. We provide a geometrical description of these viscous, resistive MRI modes and show how their physical structure is modified as a function of the Reynolds and magnetic Reynolds numbers. We demonstrate that when finite dissipative effects are considered, velocity and magnetic field disturbances are no longer orthogonal (as it is the case in the ideal MHD limit) unless the magnetic Prandtl number is unity. We generalize previous results found in the ideal limit and show that a series of key properties of the mean Reynolds and Maxwell stresses also hold for the viscous, resistive MRI. In particular, we show that the Reynolds stress is always positive and the Maxwell stress is always negative. Therefore, even in the presence of viscosity and resistivity, the total mean angular momentum transport is always directed outwards. We also find that, for any combination of the Reynolds and magnetic Reynolds numbers, magnetic disturbances dominate both the energetics and the transport of angular momentum and that the total mean energy density is an upper bound for the total mean stress responsible for angular momentum transport. The ratios between the Maxwell and Reynolds stresses and between magnetic and kinetic energy densities increase with decreasing Reynolds numbers for any magnetic Reynolds number; the lowest limit of both ratios is reached in the ideal MHD regime. △ Less

Submitted 29 January, 2008; originally announced January 2008.

Comments: ApJ, submitted, 17 pages (emulateapj style), 12 figures. Comments welcomed

Angular Momentum Transport in Accretion Disks: Scaling Laws in MRI-driven Turbulence

Authors: Martin E. Pessah, Chi-kwan Chan, Dimitrios Psaltis

Abstract: We present a scaling law that predicts the values of the stresses obtained in numerical simulations of saturated MRI-driven turbulence in non-stratified shearing boxes. It relates the turbulent stresses to the strength of the vertical magnetic field, the sound speed, the vertical size of the box, and the numerical resolution and predicts accurately the results of 35 numerical simulations perform… ▽ More We present a scaling law that predicts the values of the stresses obtained in numerical simulations of saturated MRI-driven turbulence in non-stratified shearing boxes. It relates the turbulent stresses to the strength of the vertical magnetic field, the sound speed, the vertical size of the box, and the numerical resolution and predicts accurately the results of 35 numerical simulations performed for a wide variety of physical conditions. We use our result to show that the saturated stresses in simulations with zero net magnetic flux depend linearly on the numerical resolution and would become negligible if the resolution were set equal to the natural dissipation scale in astrophysical disks. We conclude that, in order for MRI-driven turbulent angular momentum transport to be able to account for the large value of the effective alpha viscosity inferred observationally, the disk must be threaded by a significant vertical magnetic field and the turbulent magnetic energy must be in near equipartition with the thermal energy. This result has important implications for the spectra of accretion disks and their stability. △ Less

Submitted 20 August, 2007; v1 submitted 2 May, 2007; originally announced May 2007.

Comments: Accepted for publication in ApJ Letters. Two color figures. References added

The fundamental difference between shear alpha viscosity and turbulent magnetorotational stresses

Authors: Martin E. Pessah, Chi-kwan Chan, Dimitrios Psaltis

Abstract: Numerical simulations of turbulent, magnetized, differentially rotating flows driven by the magnetorotational instability are often used to calculate the effective values of alpha viscosity that is invoked in analytical models of accretion discs. In this paper we use various dynamical models of turbulent magnetohydrodynamic stresses, as well as numerical simulations of shearing boxes, to show th… ▽ More Numerical simulations of turbulent, magnetized, differentially rotating flows driven by the magnetorotational instability are often used to calculate the effective values of alpha viscosity that is invoked in analytical models of accretion discs. In this paper we use various dynamical models of turbulent magnetohydrodynamic stresses, as well as numerical simulations of shearing boxes, to show that angular momentum transport in MRI-driven accretion discs cannot be described by the standard model for shear viscosity. In particular, we demonstrate that turbulent magnetorotational stresses are not linearly proportional to the local shear and vanish identically for angular velocity profiles that increase outwards. △ Less

Submitted 18 February, 2008; v1 submitted 15 December, 2006; originally announced December 2006.

Comments: Updated to match published version in MNRAS, 8 pages, 1 figure

A Local Model for Angular Momentum Transport in Accretion Disks Driven by the Magnetorotational Instability

Authors: Martin E. Pessah, Chi-kwan Chan, Dimitrios Psaltis

Abstract: We develop a local model for the exponential growth and saturation of the Reynolds and Maxwell stresses in turbulent flows driven by the magnetorotational instability. We first derive equations that describe the effects of the instability on the growth and pumping of the stresses. We highlight the relevance of a new type of correlations that couples the dynamical evolution of the Reynolds and Ma… ▽ More We develop a local model for the exponential growth and saturation of the Reynolds and Maxwell stresses in turbulent flows driven by the magnetorotational instability. We first derive equations that describe the effects of the instability on the growth and pumping of the stresses. We highlight the relevance of a new type of correlations that couples the dynamical evolution of the Reynolds and Maxwell stresses and plays a key role in developing and sustaining the magnetorotational turbulence. We then supplement these equations with a phenomenological description of the triple correlations that lead to a saturated turbulent state. We show that the steady-state limit of the model describes successfully the correlations among stresses found in numerical simulations of shearing boxes. △ Less

Submitted 18 October, 2006; originally announced October 2006.

Comments: 4 pages, 2 figures. Physical Review Letters, in press

Journal ref: Phys.Rev.Lett.97:221103,2006

The signature of the magnetorotational instability in the Reynolds and Maxwell stress tensors in accretion discs

Authors: Martin E. Pessah, Chi-kwan Chan, Dimitrios Psaltis

Abstract: The magnetorotational instability is thought to be responsible for the generation of magnetohydrodynamic turbulence that leads to enhanced outward angular momentum transport in accretion discs. Here, we present the first formal analytical proof showing that, during the exponential growth of the instability, the mean (averaged over the disc scale-height) Reynolds stress is always positive, the me… ▽ More The magnetorotational instability is thought to be responsible for the generation of magnetohydrodynamic turbulence that leads to enhanced outward angular momentum transport in accretion discs. Here, we present the first formal analytical proof showing that, during the exponential growth of the instability, the mean (averaged over the disc scale-height) Reynolds stress is always positive, the mean Maxwell stress is always negative, and hence the mean total stress is positive and leads to a net outward flux of angular momentum. More importantly, we show that the ratio of the Maxwell to the Reynolds stresses during the late times of the exponential growth of the instability is determined only by the local shear and does not depend on the initial spectrum of perturbations or the strength of the seed magnetic. Even though we derived these properties of the stress tensors for the exponential growth of the instability in incompressible flows, numerical simulations of shearing boxes show that this characteristic is qualitatively preserved under more general conditions, even during the saturated turbulent state generated by the instability. △ Less

Submitted 20 July, 2006; v1 submitted 8 March, 2006; originally announced March 2006.

Comments: 9 pages, 4 figures. Minor revisions. Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:183-190,2006

Mass Measurements of AGN from Multi-Lorentzian Models of X-ray Variability. I. Sampling Effects in Theoretical Models of the rms^2-M_BH Correlation

Authors: Martin E. Pessah

Abstract: Recent X-ray variability studies suggest that the log of the square of the fractional rms variability amplitude, rms^2, seems to correlate with the log of the AGN black-hole mass, M_BH, with larger black holes being less variable for a fixed time interval. This has motivated the theoretical modeling of the rms^2-M_BH correlation with the aim of constraining AGN masses based on X-ray variability.… ▽ More Recent X-ray variability studies suggest that the log of the square of the fractional rms variability amplitude, rms^2, seems to correlate with the log of the AGN black-hole mass, M_BH, with larger black holes being less variable for a fixed time interval. This has motivated the theoretical modeling of the rms^2-M_BH correlation with the aim of constraining AGN masses based on X-ray variability. A viable approach to addressing this problem is to assume an underlying power spectral density with a suitable mass dependence, derive the functional form of the rms^2-M_BH correlation for a given sampling pattern, and investigate whether the result is consistent with the observations. For simplicity, previous studies, inspired by the similarities shared by the timing properties of AGN and X-ray binaries, have explored model power spectral densities characterized by broken power laws. and ignored, in general, the distorting effects that the particular sampling pattern imprints in the observed power spectral density. Motivated by the latest timing results from X-ray binaries, obtained with RXTE, we propose that AGN broad-band noise spectra consist of a small number of Lorentzian components. This assumption allows, for the first time, to fully account for sampling effects in theoretical models of X-ray variability in an analytic manner. We show that, neglecting sampling effects when deriving the fractional rms from the model power spectral density can lead to underestimating it by a factor of up to 80% with respect to its true value for the typical sampling patterns used to monitor AGN. We discuss the implications of our results for the derivation of AGN masses using theoretical models of the rms^2-M_BH correlation. (Abridged) △ Less

Submitted 11 October, 2006; v1 submitted 13 September, 2005; originally announced September 2005.

Comments: The Astrophysical Journal, in press, 11 pages, 6 figures

Journal ref: Astrophys.J.655:66-76,2007

The Stability of Magnetized Rotating Plasmas with Superthermal Fields

Authors: Martin E. Pessah, Dimitrios Psaltis

Abstract: During the last decade it has become evident that the magnetorotational instability is at the heart of the enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the evolution of the magnetorotational instability beyond the weak-field l… ▽ More During the last decade it has become evident that the magnetorotational instability is at the heart of the enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the evolution of the magnetorotational instability beyond the weak-field limit. We show that, when superthermal toroidal fields are considered, the effects of both compressibility and magnetic tension forces, which are related to the curvature of toroidal field lines, should be taken fully into account. We demonstrate that the presence of a strong toroidal component in the magnetic field plays a non-trivial role. When strong fields are considered, the strength of the toroidal magnetic field not only modifies the growth rates of the unstable modes but also determines which modes are subject to instabilities. We find that, for rotating configurations with Keplerian laws, the magnetorotational instability is stabilized at low wavenumbers for toroidal Alfven speeds exceeding the geometric mean of the sound speed and the rotational speed. We discuss the significance of our findings for the stability of cold, magnetically dominated, rotating fluids and argue that, for these systems, the curvature of toroidal field lines cannot be neglected even when short wavelength perturbations are considered. We also comment on the implications of our results for the validity of shearing box simulations in which superthermal toroidal fields are generated. △ Less

Submitted 13 April, 2005; v1 submitted 2 June, 2004; originally announced June 2004.

Comments: 24 pages, 12 figures. Accepted for publication in ApJ. Sections 2 and 5 substantially expanded, added Appendix A and 3 figures with respect to previous version. Animations are available at http://www.physics.arizona.edu/~mpessah/research/

Journal ref: Astrophys.J. 628 (2005) 879-901

Gravitational microlensing of gamma-ray blazars

Authors: Diego F. Torres, Gustavo E. Romero, Ernesto F. Eiroa, Joachim Wambsganss, M. E. Pessah

Abstract: We present a detailed study of the effects of gravitational microlensing on compact and distant $γ$-ray blazars. These objects have $γ$-ray emitting regions which are small enough as to be affected by microlensing effects produced by stars lying in intermediate galaxies. We analyze the temporal evolution of the gamma-ray magnification for sources moving in a caustic pattern field, where the comb… ▽ More We present a detailed study of the effects of gravitational microlensing on compact and distant $γ$-ray blazars. These objects have $γ$-ray emitting regions which are small enough as to be affected by microlensing effects produced by stars lying in intermediate galaxies. We analyze the temporal evolution of the gamma-ray magnification for sources moving in a caustic pattern field, where the combined effects of thousands of stars are taken into account using a numerical technique. We propose that some of the unidentified $γ$-ray sources (particularly some of those lying at high galactic latitude whose gamma-ray statistical properties are very similar to detected $γ$-ray blazars) are indeed the result of gravitational lensing magnification of background undetected Active Galactic Nuclei (AGNs). △ Less

Submitted 24 May, 2002; originally announced May 2002.

Comments: 30 pages, 27 figures. Four figures are being submitted only as .gif files, and should be printed separately. The abstract below has been shortened from the actual version appearing in the paper

Journal ref: Mon.Not.Roy.Astron.Soc.339:335-352,2003

Statistical mechanics and the description of the early universe II. Principle of detailed balance and primordial 4He formation

Authors: M. E. Pessah, Diego F. Torres

Abstract: If the universe is slightly non-extensive, and the distribution functions are not exactly given by those of Boltzmann-Gibbs, the primordial production of light elements will be non-trivially modified. In particular, the principle of detailed balance (PDB), of fundamental importance in the standard analytical analysis, is no longer valid, and a non-extensive correction appears. This correction is… ▽ More If the universe is slightly non-extensive, and the distribution functions are not exactly given by those of Boltzmann-Gibbs, the primordial production of light elements will be non-trivially modified. In particular, the principle of detailed balance (PDB), of fundamental importance in the standard analytical analysis, is no longer valid, and a non-extensive correction appears. This correction is computed and its influence is studied and compared with previous works, where, even when the universe was considered as an slightly non-extensive system, the PDB was assumed valid. We analytically track the formation of Helium and Deuterium, and study the kind of deviation one could expect from the standard regime. The correction to the capture time, the moment in which Deuterium can no longer be substantially photo-disintegrated, is also presented. This allows us to take into account the process of the free decay of neutrons, which was absent in all previous treatments of the topic. We show that even when considering a first (linear) order correction in the quantum distribution functions, the final output on the primordial nucleosynthesis yields cannot be reduced to a linear correction in the abundances. We finally obtain new bounds upon the non-extensive parameter, both comparing the range of physical viability of the theory, and using the latest observational data. △ Less

Submitted 4 May, 2001; originally announced May 2001.

Comments: 24 pages, to appear in Physica A (2001)

Journal ref: Physica A297 (2001) 201-228