BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Prateek Sharma

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.HEastro-ph.COastro-phphysics.flu-dynDistributed, Parallel, and Cluster Computingphysics.comp-phphysics.plasm-phastro-ph.IMeess.SYgr-qcmath-phmath.MPmath.NANumerical AnalysisPerformancephysics.ao-phphysics.space-phSystems and Control

Catalog footprint

What is connected

51works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

51 published item(s)

preprint2023arXiv

Dissipation of AGN jets in a clumpy interstellar medium

Accreting supermassive black holes (SMBHs) frequently power jets that interact with the interstellar/circumgalactic medium (ISM/CGM), regulating star-formation in the galaxy. Highly supersonic jets launched by active galactic nuclei (AGN) power a cocoon that confines them and shocks the ambient medium. We build upon the models of narrow conical jets interacting with a smooth ambient medium, to include the effect of dense clouds that are an essential ingredient of a multiphase ISM. The key physical ingredient of this model is that the clouds along the supersonic jet-beam strongly decelerate the jet-head, but the subsonic cocoon easily moves around the clouds without much resistance. We propose scalings for important physical quantities -- cocoon pressure, head & cocoon speed, and jet radius. We obtain, for the first time, the analytic condition on clumpiness of the ambient medium for the jet to dissipate within the cocoon and verify it with numerical simulations of conical jets interacting with a uniform ISM with embedded spherical clouds. A jet is defined to be dissipated when the cocoon speed exceeds the speed of the jet-head. We compare our models to more sophisticated numerical simulations, direct observations of jet-ISM interaction (e.g., quasar J1316+1753), and discuss implications for the Fermi/eROSITA bubbles. Our work also motivates effective subgrid models for AGN jet feedback in a clumpy ISM unresolved by the present generation of cosmological galaxy formation simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmic rays from massive star clusters : A close look at Westerlund 1

We study the effect of cosmic ray (CR) acceleration in the massive compact star cluster Westerlund 1 in light of its recent detection in $γ$-rays. Recent observations reveal a $1/r$ radial distribution of the CR energy density. Here we theoretically investigate whether or not this profile can help to distinguish between (1) continuous CR acceleration in the star cluster stellar wind-driven shocks and (2) discrete CR acceleration in multiple supernovae shocks -- which are often debated in the literature. Using idealized two-fluid simulations and exploring different acceleration sites and diffusion coefficients, we obtain the CR energy density profile and luminosity to find the best match for the $γ$-ray observations. We find that the inferred CR energy density profiles from observations of $γ$-ray luminosity and mass can be much different from the true radial profile. CR acceleration at either the cluster core region or the wind termination shock can explain the observations, if the diffusion coefficient is $κ_{\rm cr}\sim 10^{27}$ cm$^2$ s$^{-1}$ and a fraction of $\approx 10\%-20\%$ of the shock power/post-shock pressure is deposited into the CR component. We also study the possibility of discrete supernovae (SN) explosions being responsible for CR acceleration and find that with an injection rate of 1 SN in every $\sim 0.03$ Myr, one can explain the observed $γ$-ray profile. This multiple SN scenario is consistent with X-ray observations only if the thermal conductivity is close to the Spitzer value.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological evolution of gas and supermassive black holes in idealized isolated halos

We study the evolution of baryonic gas in cosmologically growing dark matter halos. To accurately model both the inner and outer regions of the halos, we use a dark matter density profile that transitions smoothly from the NFW profile within the virial radius to a more realistic flat profile far beyond the halo. We construct a dark matter gravitational potential consistent with this density profile, and we use a "cosmological" potential that accounts for gas evolution consistent with Hubble expansion at large radii. Gas is initialized with a density $\approx$ 0.2 times the dark matter density, consistent with the universal baryon fraction $ρ_{\rm g}/(ρ_{\rm g}+ρ_{\rm DM}) \approx 0.17$. We study the formation of the virial shock and evolution of the baryon fraction, including the effects of radiative cooling and AGN jet feedback. The feedback is powered by the accretion of cold gas onto a central supermassive black hole (SMBH). The cores of the halo exhibit heating and cooling cycles, whose strength and duration depend on the feedback efficiency and the halo mass. The central SMBH initially grows exponentially with time in the early quasar phase, but the growth slows down at later times. The baryon fraction in the core decreases with increasing feedback efficiency and decreasing halo mass. While the halo outskirts evolve self-similarly, the core density is non-evolving, in agreement with cluster observations. We analyze the correlations between the properties of the gas and the central SMBH, and explore the existence of a fundamental plane.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deprojection of X-ray data in galaxy clusters: confronting simulations with observations

Numerical simulations with varying realism indicate an emergent principle -- multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20-30, but the prevalence of cores with $t_{\rm cool}/t_{\rm ff}$ ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3-D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining $t_{\rm cool}/t_{\rm ff}$. We show that while the used methods are able to recover the $t_{\rm cool}/t_{\rm ff}$ ratio for relaxed clusters, they have an uncertainty of a factor of $2-3$ in systems containing large cavities ($\gtrsim 20$ kpc). We also show that the mass estimates from these methods, in the absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass ($M_{200}$) and the concentration parameter ($c$) in the form of $M_{200} c^2 \approx$ constant. Additionally, lack of soft-X-ray ($\lesssim 0.5$ keV) coverage and poor spatial resolution make us overestimate min($t_{\rm cool}/t_{\rm ff}$) by a factor of few in clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$. This bias can largely explain the lack of cool-core clusters with min($t_{\rm cool}/t_{\rm ff}$) $\lesssim 5$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiphase turbulence in galactic halos: effect of the driving

Supernova explosions, active galactic nuclei jets, galaxy--galaxy interactions and cluster mergers can drive turbulence in the circumgalactic medium (CGM) and in the intracluster medium (ICM). However, the exact nature of turbulence forced by these sources and its impact on the different statistical properties of the CGM/ICM and their global thermodynamics is still unclear. To investigate the effects of different types of forcing, we conduct high resolution ($1008^3$ resolution elements) idealised hydrodynamic simulations with purely solenoidal (divergence-free) forcing, purely compressive (curl-free) forcing, and natural mixture forcing (equal fractions of the two components). The simulations also include radiative cooling. We study the impact of the three different forcing modes (sol, comp, mix) on the morphology of the gas, its temperature and density distributions, sources and sinks of enstrophy, i.e., solenoidal motions, as well as the kinematics of hot ($\sim10^7~\mathrm{K}$) X-ray emitting and cold ($\sim10^4~\mathrm{K}$) H$α$ emitting gas. We find that compressive forcing leads to stronger variations in density and temperature of the gas as compared to solenoidal forcing. The cold phase gas forms large-scale filamentary structures for compressive forcing and misty, small-scale clouds for solenoidal forcing. The cold phase gas has stronger large-scale velocities for compressive forcing. The natural mixture forcing shows kinematics and gas distributions intermediate between the two extremes, the cold-phase gas occurs as both large-scale filaments and small-scale misty clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shadow and weak gravitational lensing of a rotating regular black hole in a non-minimally coupled Einstein-Yang-Mills theory in the presence of plasma

The null geodesics of the regular and rotating magnetically charged black hole in a non-minimally coupled Einstein-Yang-Mills theory surrounded by a plasma medium is studied. The effect of magnetic charge and Yang-Mills parameter on the effective potential and radius of photon orbits has investigated. We then study the shadow of a regular and rotating magnetically charged black hole along with the observables in presence of the plasma medium. The presence of plasma medium affects the apparent size of the shadow of a regular rotating black hole in comparison to vacuum case. Variation of shadow radius and deformation parameter with Yang-Mills and plasma parameter has examined. Furthermore, the deflection angle of the massless test particles in weak field approximation around this black hole spacetime in presence of homogeneous plasma medium is also investigated. Finally, we have compared the obtained results with Kerr-Newman and Schwarzschild black hole solutions in general relativity (GR).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A method for determining the parameters in a rheological model for viscoelastic materials by minimizing Tikhonov functionals

Mathematical models describing the behavior of viscoelastic materials are often based on evolution equations that measure the change in stress depending on its material parameters such as stiffness, viscosity or relaxation time. In this article, we introduce a Maxwell-based rheological model, define the associated forward operator and the inverse problem in order to determine the number of Maxwell elements and the material parameters of the underlying viscoelastic material. We perform a relaxation experiment by applying a strain to the material and measure the generated stress. Since the measured data varies with the number of Maxwell elements, the forward operator of the underlying inverse problem depends on parts of the solution. By introducing assumptions on the relaxation times, we propose a clustering algorithm to resolve this problem. We provide the calculations that are necessary for the minimization process and conclude with numerical results by investigating unperturbed as well as noisy data. We present different reconstruction approaches based on minimizing a least squares functional. Furthermore, we look at individual stress components to analyze different displacement rates. Finally, we study reconstructions with shortened data sets to obtain assertions on how long experiments have to be performed to identify conclusive material parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks

Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modeling, and future directions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterising the Dynamo in a Radiatively Inefficient Accretion Flow

We explore the MRI driven dynamo in a radiatively inefficient accretion flow (RIAF) using the mean field dynamo paradigm. Using singular value decomposition (SVD) we obtain the least squares fitting dynamo coefficients $α$ and $γ$ by comparing the time series of the turbulent electromotive force and the mean magnetic field. Our study is the first one to show the poloidal distribution of these dynamo coefficients in global accretion flow simulations. Surprisingly, we obtain a high value of the turbulent pumping coefficient $γ$ which transports the mean magnetic flux radially outward. This would have implications for the launching of magnetised jets which are produced efficiently in presence a large-scale poloidal magnetic field close to the compact object. We present a scenario of a truncated disc beyond the RIAF where a large scale dynamo-generated poloidal magnetic field can aid jet-launching close to the black hole. Magnitude of all the calculated coefficients decreases with radius. Meridional variations of $α_{ϕϕ}$, responsible for toroidal to poloidal field conversion, is very similar to that found in shearing box simulations using the `test field' (TF) method. By estimating the relative importance of $α$-effect and shear, we conclude that the MRI driven large-scale dynamo, which operates at high latitudes beyond a disc scale height, is essentially of the $α-Ω$ type.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cloud-scale VM Deflation for Running Interactive Applications On Transient Servers

Transient computing has become popular in public cloud environments for running delay-insensitive batch and data processing applications at low cost. Since transient cloud servers can be revoked at any time by the cloud provider, they are considered unsuitable for running interactive application such as web services. In this paper, we present VM deflation as an alternative mechanism to server preemption for reclaiming resources from transient cloud servers under resource pressure. Using real traces from top-tier cloud providers, we show the feasibility of using VM deflation as a resource reclamation mechanism for interactive applications in public clouds. We show how current hypervisor mechanisms can be used to implement VM deflation and present cluster deflation policies for resource management of transient and on-demand cloud VMs. Experimental evaluation of our deflation system on a Linux cluster shows that microservice-based applications can be deflated by up to 50\% with negligible performance overhead. Our cluster-level deflation policies allow overcommitment levels as high as 50\%, with less than a 1\% decrease in application throughput, and can enable cloud platforms to increase revenue by 30\%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cool-Core Cycles and Phoenix

Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling The Temporally Constrained Preemptions of Transient Cloud VMs

Transient cloud servers such as Amazon Spot instances, Google Preemptible VMs, and Azure Low-priority batch VMs, can reduce cloud computing costs by as much as $10\times$, but can be unilaterally preempted by the cloud provider. Understanding preemption characteristics (such as frequency) is a key first step in minimizing the effect of preemptions on application performance, availability, and cost. However, little is understood about temporally constrained preemptions---wherein preemptions must occur in a given time window. We study temporally constrained preemptions by conducting a large scale empirical study of Google's Preemptible VMs (that have a maximum lifetime of 24 hours), develop a new preemption probability model, new model-driven resource management policies, and implement them in a batch computing service for scientific computing workloads. Our statistical and experimental analysis indicates that temporally constrained preemptions are not uniformly distributed, but are time-dependent and have a bathtub shape. We find that existing memoryless models and policies are not suitable for temporally constrained preemptions. We develop a new probability model for bathtub preemptions, and analyze it through the lens of reliability theory. To highlight the effectiveness of our model, we develop optimized policies for job scheduling and checkpointing. Compared to existing techniques, our model-based policies can reduce the probability of job failure by more than $2\times$. We also implement our policies as part of a batch computing service for scientific computing applications, which reduces cost by $5\times$ compared to conventional cloud deployments and keeps performance overheads under $3\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On modelling CC85 wind in an expanding local box

The Chevalier & Clegg (1985) spherical wind is widely used to model galactic outflows. Efforts are underway to understand the production of multiphase gas in such outflows. Since the important scales for the multiphase gas are very small compared to the size of the wind, it is very useful to have a local model of the cold gas interacting with a volume-filling expanding wind. In this note we explicitly derive the governing hydrodynamic equations in such a local box for the first time. We highlight that the expansion in a wind occurs only in the transverse direction, and this results in anisotropic source terms in the momentum equation. Our equations will be very useful in realistic modelling of the small scale structure of multiphase galactic outflows. [original note does not include an abstract]

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Resolving small-scale cold circumgalactic gas in TNG50

We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M* > 10^11 Msun) at intermediate redshift (z~0.5). We discover a significant abundance of small-scale, cold gas structure in the CGM of 'red and dead' elliptical systems, as traced by neutral HI and MgII. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong drho/rho >> 1 gas density perturbations which stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ~10^7 K to radiatively inefficient ~10^4 K clouds, which act as cosmologically long-lived, 'stimulated cooling' seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma beta << 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the ~8x10^4 Msun cell mass of TNG50 enables the ~few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulence in stratified atmospheres: implications for the intracluster medium

The gas motions in the intracluster medium (ICM) are governed by stratified turbulence. Stratified turbulence is fundamentally different from Kolmogorov (isotropic, homogeneous) turbulence; kinetic energy not only cascades from large to small scales, but it is also converted into buoyancy potential energy. To understand the density and velocity fluctuations in the ICM, we conduct high-resolution ($1024^2\times 1536$ grid points) hydrodynamical simulations of subsonic turbulence (with rms Mach number $\mathcal{M}\approx 0.25$) and different levels of stratification, quantified by the Richardson number $\mathrm{Ri}$, from $\mathrm{Ri}=0$ (no stratification) to $\mathrm{Ri}=13$ (strong stratification). We quantify the density, pressure and velocity fields for varying stratification because observational studies often use surface brightness fluctuations to infer the turbulent gas velocities of the ICM. We find that the standard deviation of the logarithmic density fluctuations ($σ_s$), where $s=\ln(ρ/\left<ρ(z)\right>)$, increases with $\mathrm{Ri}$. For weakly stratified subsonic turbulence ($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new $σ_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation, $σ_s^2=\ln(1+b^2\mathcal{M}^4+0.09\mathcal{M}^2\mathrm{Ri}H_P/H_S)$, where $b=1/3$--$1$ is the turbulence driving parameter, and $H_P$ and $H_S$ are the pressure and entropy scale heights respectively. We further find that the power spectrum of density fluctuations, $P(ρ_k/\left<ρ\right>)$, increases in magnitude with increasing $\mathrm{Ri}$, whereas the velocity power spectrum is invariant. Thus, the ratio between density and velocity power spectra strongly depends on $\mathrm{Ri}$. Pressure fluctuations, on the other hand, are independent of stratification and only depend on $\mathcal{M}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Radio halos of star forming galaxies

We study the synchrotron radio emission from extra-planar regions of star forming galaxies. We use ideal magneto-hydrodynamical (MHD) simulations of a rotating Milky Way-type disk galaxy with distributed star formation sites for three star formation rates (SFRs) (0.3, 3, 30 M$_{\odot}$ yr$^{-1}$). From our simulations, we see emergence of galactic-scale magnetised outflows, carrying gas from the disk. We compare the morphology of the outflowing gas with hydrodynamic (HD) simulations. We look at the spatial distribution of magnetic field in the outflows. Assuming that a certain fraction of gas energy density is converted into cosmic ray energy density, and using information about the magnetic field, we obtain synchrotron emissivity throughout the simulation domain. We generate the surface brightness maps at a frequency of 1.4 GHz. The outflows are more extended in the vertical direction than radial and hence have an oblate shape. We further find that the matter right behind the outer shock, shines brighter in these maps than that above or below. To understand whether this feature can be observed, we produce vertical intensity profiles. We convolve the vertical intensity profile with the typical beam sizes of radio telescopes, for a galaxy located at 10 Mpc (similar to NGC 891) in order to estimate the radio scale height to compare with observations. We find that for our SFRs this feature will lie below the RMS noise limit of instruments. The radio scale height is found to be $\sim 300-1200$ pc , depending on the resolution of the telescope. We relate the advection speed of the outer shock with the surface density of star formation as $\rm{v}_{\rm adv} \propto Σ_{\rm SFR}^{0.3}$ which is consistent with earlier observations and analytical estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Realistic modeling of wind and supernovae shocks in star clusters: addressing ${\rm ^{22}Ne/^{20}Ne}$ and other problems in Galactic cosmic rays

Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $γ$-ray bright ($L_γ\gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Clusters of galaxies and the cosmic web with SKA

The intra-cluster and inter-galactic media (ICM, IGM) that pervade the large scale structure of the Universe are known to be magnetised at sub-micro Gauss to micro Gauss levels and to contain cosmic rays (CRs). The acceleration of CRs and their evolution along with that of magnetic fields in these media is still not well understood. Diffuse radio sources of synchrotron origin associated with the ICM such as radio halos, relics and mini-halos are direct probes of the underlying mechanisms of CR acceleration. Observations with radiotelescopes such as the GMRT, the VLA and the WSRT (0.15 - 2 GHz) have revealed scaling relations between the thermal and non-thermal properties of clusters and favour the role of shocks in the formation of radio relics and of turbulent re-acceleration in the formation of radio halos and mini-halos. Due to the limitations of current radio telescopes, wide-band studies and exploration of low mass and supercluster-scale systems is difficult. The Square Kilometer Array (SKA) is a next generation radio telescope that will operate in the frequency range of 0.05 - 20 GHz with unprecedented sensitivities and resolutions. The expected detection limits of SKA will reveal a few hundred to thousand new radio halos, relics and mini-halos providing the first large and comprehensive samples for their study. The wide frequency coverage along with sensitivity to extended structures will be able to constrain the CR acceleration mechanisms. The higher frequency (> 5 GHz) observations will be able to use the Sunyaev-Zel'dovich effect to probe the ICM pressure in addition to the tracers such as lobes of head-tail radio sources. The SKA also opens prospects to detect the "off-state" radio emission from the ICM predicted by the hadronic models and the turbulent re-acceleration models. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cold gas in cluster cores: Global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g., spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in halos critically depends on the ratio of the cooling time to the free-fall time ($t_{cool}/t_{ff}$). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focussed on {\em in situ} condensation of multiphase gas but spherical simulations studied condensation {\em anywhere} in the box. Unlike previous claims, our nonlinear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas condensing depends on the shape of the gravitational potential well; gas has more time to condense if gravitational acceleration decreases toward the center. In our idealized simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive $t_{cool}/t_{ff}$ far beyond the critical value close to 10. Triggered condensation is very prominent in plane-parallel simulations, in which a large amount of cold gas condenses out.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detection of polarized quasi-periodic microstructure emission in millisecond pulsars

Microstructure emission, involving short time scale, often quasi-periodic, intensity fluctuations in subpulse emission, is well known in normal period pulsars. In this letter, we present the first detections of quasi-periodic microstructure emission from millisecond pulsars (MSPs), from Giant Metrewave Radio Telescope (GMRT) observations of two MSPs at 325 and 610 MHz. Similar to the characteristics of microstructure observed in normal period pulsars, we find that these features are often highly polarized, and exhibit quasi-periodic behavior on top of broader subpulse emission, with periods of the order of a few $μ$s. By measuring their widths and periodicities from single pulse intensity profiles and their autocorrelation functions, we extend the microstructure timescale - rotation period relationship by more than an order of magnitude down to rotation periods $\sim$ 5 ms, and find it to be consistent with the relationship derived earlier for normal pulsars. The similarity of behavior is remarkable, given the significantly different physical properties of MSPs and normal period pulsars, and rules out several previous speculations about the possible different characteristics of microstructure in MSP radio emission. We discuss the possible reasons for the non-detection of these features in previous high time resolution MSP studies along with the physical implications of our results, both in terms of a geometric beam sweeping model and temporal modulation model for micropulse production.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Diffuse X-ray emission from star forming galaxies

We study the diffuse X-ray luminosity ($L_X$) of star forming galaxies using 2-D axisymmetric hydrodynamical simulations and analytical considerations of supernovae (SNe) driven galactic outflows. We find that the mass loading of the outflows, a crucial parameter for determining the X-ray luminosity, is constrained by the availability of gas in the central star forming region, and a competition between cooling and expansion. We show that the allowed range of the mass loading factor can explain the observed scaling of $L_X$ with star formation rate (SFR) as $L_X \propto$ SFR$^2$ for SFR $\gtrsim 1$ M$_\odot$yr$^{-1}$, and a flatter relation at low SFRs. We also show that the emission from the hot circumgalactic medium (CGM) in the halo of massive galaxies can explain the sub-linear behaviour of the $L_X-$SFR relation as well as a large scatter in the diffuse X-ray emission for low SFRs ($\lesssim$ few M$_\odot$yr$^{-1}$). Our results point out that galaxies with small SFRs and large diffuse X-ray luminosities are excellent candidates for detection of the elusive CGM.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How radiation affects superbubbles: Through momentum injection in early phase and photo-heating thereafter

Energetic winds and radiation from massive star clusters push the surrounding gas and blow superbubbles in the interstellar medium (ISM). Using 1-D hydrodynamic simulations, we study the role of radiation in the dynamics of superbubbles driven by a young star cluster of mass $10^{6}$ M$_{\odot}$. We have considered a realistic time evolution of the mechanical power as well as radiation power of the star cluster, and detailed heating and cooling processes. We find that the ratio of the radiation pressure on the shell (shocked ISM) to the thermal pressure ($\sim10^{7}$ K) of the shocked wind region is almost independent of the ambient density, and it is greater than unity before $\lesssim 1$ Myr. We explore the parameter space of density and dust opacity of the ambient medium, and find that the size of the hot gas ($\sim$ 10$^{7}$ K) cavity is insensitive to the dust opacity ($σ_{d}\approx(0.1-1.5)\times 10^{-21}$ cm$^{2}$), but the structure of the photoionized ($\sim10^4$ K) gas depends on it. Most of the radiative losses occur at $\sim10^{4}$ K, with sub-dominant losses at $\lesssim 10^3$ K and $\sim10^{6}-10^{8}$ K. The superbubbles can retain as high as $\sim 10\%$ of its input energy, for an ambient density of $10^{3}\,m{\rm_{H}\,cm^{-3}}$. We discuss the role of ionization parameter and recombination-averaged density in understanding the dominant feedback mechanism. Finally, we compare our results with the observations of 30 Doradus.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Molecular outflows in starburst nuclei

Recent observations have detected molecular outflows in a few nearby starburst nuclei. We discuss the physical processes at work in such an environment in order to outline a scenario that can explain the observed parameters of the phenomenon, such as the molecular mass, speed and size of the outflows. We show that outflows triggered by OB associations, with $N_{OB}\ge 10^5$ (corresponding to a star formation rate (SFR)$\ge 1$ M$_{\odot}$ yr$^{-1}$ in the nuclear region), in a stratified disk with mid-plane density $n_0\sim 200\hbox{--}1000$ cm$^{-3}$ and scale height $z_0\ge 200 (n_0/10^2 \, {\rm cm}^{-3})^{-3/5}$ pc, can form molecules in a cool dense and expanding shell. The associated molecular mass is $\ge 10^7$ M$_\odot$ at a distance of a few hundred pc, with a speed of several tens of km s$^{-1}$. We show that a SFR surface density of $10 \le Σ_{SFR} \le 50$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$ favours the production of molecular outflows, consistent with observed values.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spherical accretion: the influence of inner boundary and quasi-periodic oscillations

Bondi accretion assumes that there is a sink of mass at the center -- which in case of a black hole (BH) corresponds to the advection of matter across the event horizon. Other stars, such as a neutron star (NS), have surfaces and hence the infalling matter has to slow down at the surface. We study the initial value problem in which the matter distribution is uniform and at rest at $t=0$. We consider different inner boundary conditions for BHs and NSs: outflow boundary condition (mimicking mass sink at the center) valid for BHs; and {\em reflective} and steady-shock (allowing gas to cross the inner boundary at subsonic speeds) boundary conditions for NSs. We also obtain a similarity solution for cold accretion on to BHs and NSs. 1-D simulations show the formation of an outward propagating and a standing shock in NSs for reflective and steady-shock boundary conditions, respectively. Entropy is the highest at the bottom of the subsonic region for reflective boundary conditions. In 2-D this profile is convectively unstable. Using steady-shock inner boundary conditions, the flow is unstable to the standing accretion shock instability (SASI) in 2-D, which leads to global shock oscillations and may be responsible for quasi-periodic oscillations (QPOs) seen in the lightcurves of accreting systems. For steady accretion in the quiescent state, spherical accretion rate on to a NS can be suppressed by orders of magnitude compared to that on to a BH.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cool core cycles: Cold gas and AGN jet feedback in cluster cores

Using high-resolution 3-D and 2-D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGN) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity $\approx ε\dot{M}_{\rm acc} c^2$; where $\dot{M}_{\rm acc}$ is the mass accretion rate at 1 kpc) as small as $5 \times 10^{-5}$ is sufficient to reduce the cooling/accretion rate by $\sim 10$ compared to a pure cooling flow. This value is smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted on to the supermassive black hole (SMBH). We find hysteresis cycles in all our simulations with cold mode feedback: {\em condensation} of cold gas when the ratio of the cooling-time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) is $\lesssim 10$ leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and {\em overheating} of the hot ICM such that $t_{\rm cool}/t_{\rm ff} > 10$; further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy or $t_{\rm cool}/t_{\rm ff}$. A fewer number of cycles are observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3-D simulations show the formation of a few-kpc scale, rotationally-supported, massive ($\sim 10^{11} M_\odot$) cold gas torus. (abstract abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Long way to go: how outflows from large galaxies propagate through the hot halo gas

Using hydrodynamic simulations, we study the mass loss due to supernova-driven outflows from Milky Way type disk galaxies, paying particular attention to the effect of the extended hot halo gas. We find that the total mass loss at inner radii scales roughly linearly with total mass of stars formed, and that the mass loading factor at the virial radius can be several times its value at inner radii because of the swept up hot halo gas. The temperature distribution of the outflowing material in the inner region ($\sim $10 kpc) is bimodal in nature, peaking at $10^5$ K and $10^{6.5}$ K, responsible for optical and X-ray emission, respectively. The contribution of cold/warm gas with temperature $\le 10^{5.5}$ K to the outflow rate within 10 kpc is $\approx 0.3\hbox{--}0.5$. The warm mass loading factor, $η_{3e5}$ ($T\le 3 \times 10^5$ K) is related to the mass loading factor at the virial radius ($η_{v}$) as $η_{v} \approx 25\, η_{3e5}\, \left(\mbox{SFR}/{\rm M}_\odot{\rm yr}^{-1} \right)^{-0.15}$ for a baryon fraction of 0.1 and a starburst period of 50 Myr. We also discuss the effect of multiple bursts that are separated by both short and long periods. The outflow speed at the virial radius is close to the sound speed in the hot halo, $\lesssim 200$ km s$^{-1}$. We identify two `sequences' of outflowing cold gas at small scales: a fast ($\approx 500$ km~s$^{-1}$) sequence, driven by the unshocked free-wind; and a slow sequence ($\approx \pm 100$ km s$^{-1}$) at the conical interface of the superwind and the hot halo.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-wavelength features of Fermi Bubbles as signatures of a Galactic wind

Using hydrodynamical simulations, we show for the first time that an episode of star formation in the center of the Milky Way, with a star-formation-rate (SFR) $\sim 0.5$ M$_\odot$ yr$^{-1}$ for $\sim 30$ Myr, can produce bubbles that resemble the Fermi Bubbles (FBs), when viewed from the solar position. The morphology, extent and multi-wavelength observations of FBs, especially X-rays, constrain various physical parameters such as SFR, age, and the circum-galactic medium (CGM) density. We show that the interaction of the CGM with the Galactic wind driven by a star formation in the central region can explain the observed surface brightness and morphological features of X-rays associated with the Fermi Bubbles. Furthermore, assuming that cosmic ray electrons are accelerated {\it in situ} by shocks and/or turbulence, the brightness and morphology of gamma-ray emission and the microwave haze can be explained. The kinematics of the cold and warm clumps in our model also matches with recent observations of absorption lines through the bubbles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction ($f_{\rm esc}$) of ionizing photons from starburst galaxies. Using 2-D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the center of the disk along different angles through the superbubble and the gas disk. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of $\sim 40 ^\circ$, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scale heights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed $\sim [1- \cos (1 \, {\rm radian})] = 0.5$ from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time- and angle-averaged escape fraction on the mid-plane disk gas density (in the range $n_0=0.15-50$ cm $^{-3}$) and the disk scale height (between $z_0=10-600$ pc). We find that the escape fraction is related to the disk parameters (the mid-plane disk density and scale height) roughly so that $f_{\rm esc}^αn_0^2 z_0^3$ (with $α\approx 2.2$) is a constant. For disks with a given WNM temperature, massive disks have lower escape fraction than low mass galaxies. For Milky Way ISM parameters, we find $f_{\rm esc}\sim 5\%$, and it increases to $\approx 10\%$ for a galaxy ten times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermal Conduction and Multiphase Gas in Cluster Cores

We examine the role of thermal conduction and magnetic fields in cores of galaxy clusters through global simulations of the intracluster medium (ICM). In particular, we study the influence of thermal conduction, both isotropic and anisotropic, on the condensation of multiphase gas in cluster cores. Previous hydrodynamic simulations have shown that cold gas condenses out of the hot ICM in thermal balance only when the ratio of the cooling time ($t_{\rm cool}$) and the free-fall time ($t_{\rm ff}$) is less than $\approx 10$. Since thermal conduction is significant in the ICM and it suppresses local cooling at small scales, it is imperative to include thermal conduction in such studies. We find that anisotropic (along local magnetic field lines) thermal conduction does not influence the condensation criterion for a general magnetic geometry, even if thermal conductivity is large. However, with isotropic thermal conduction cold gas condenses only if conduction is suppressed (by a factor $\lesssim 0.3$) with respect to the Spitzer value.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

In a hot bubble: why does superbubble feedback work, but isolated supernovae do not?

Using idealized 1-D Eulerian hydrodynamic simulations, we contrast the behavior of isolated supernovae with the superbubbles driven by multiple, collocated supernovae. Continuous energy injection via successive supernovae going off within the hot/dilute bubble maintains a strong termination shock. This strong shock keeps the superbubble over-pressured and drives the outer shock well after it becomes radiative. Isolated supernovae, in contrast, with no further energy injection, become radiative quite early ($\lesssim 0.1$ Myr, 10s of pc), and stall at scales $\lesssim 100$ pc. We show that isolated supernovae lose almost all of their mechanical energy by a Myr, but superbubbles can retain up to $\sim 40\%$ of the input energy in form of mechanical energy over the lifetime of the star cluster (few 10s of Myr). These conclusions hold even in the presence of realistic magnetic fields and thermal conduction. We also compare various recipes for implementing supernova feedback in numerical simulations. For various feedback prescriptions we derive the spatial scale below which the energy needs to be deposited for it to couple to the interstellar medium (ISM). We show that a steady thermal wind within the superbubble appears only for a large number ($\gtrsim 10^4$) of supernovae. For smaller clusters we expect multiple internal shocks instead of a smooth, dense thermalized wind.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The cold mode: A phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a nonlinear phenomenological model for the evolution of density perturbations in the ICM due to local thermal instability and gravity. We have analyzed and extended a model for the evolution of an over-dense blob in the ICM. We find two regimes in which the over-dense blobs can cool to thermally stable low temperatures. One for large $t_{\rm {cool}} / t_{\rm {ff}}$ ($t_{\rm {cool}}$ is the cooling time and $t_{\rm {ff}}$ is the free fall time), where a large initial over-density is required for thermal runaway to occur; this is the regime which was previously analyzed in detail. We discover a second regime for $t_{\rm {cool}} / t_{\rm {ff}} \lesssim 1$ (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if $t_{\rm {cool}} / t_{\rm {ff}} \lesssim 10$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be $\gtrsim$ the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to {\em mixing} with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Astrophysical Coronae: Lessons from Modeling of the Intracluster Medium

Coronae exist in most astrophysical objects: stars, accretion disks, and individual galaxies and clusters of galaxies. Coronae in these varied systems have some common properties: 1) hydrostatic equilibrium in background gravity is a good assumption; 2) they are optically thin, i.e., photons escape as soon as they are born; 3) they have cooling times shorter than their ages, and thus require heating for sustenance. Generally the coronal heating mechanisms are quite complex but the structure of the corona is tightly constrained by the interplay of cooling, {\em global} heating, and background gravity. We briefly summarize the results from our studies of the intracluster medium (the cluster corona) and draw inferences which should apply to most astrophysical coronae.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Isotropic Heating of Galaxy Cluster Cores via Rapidly Reorienting AGN Jets

AGN jets carry more than sufficient energy to stave off catastrophic cooling of the intracluster medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, the ICM must be heated in a near-isotropic fashion and narrow bipolar jets with $P_{\rm jet}=10^{44-45}$ ergs/s, typical of radio AGNs at cluster centres, are inefficient at heating the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster cores, the SMBHs will, in addition to accreting gas via radiatively inefficient flows, experience short stochastic episodes of enhanced accretion via thin discs. In general, the orientation of these accretion discs will be misaligned with the spin axis of the black holes and the ensuing torques will cause the black hole's spin axis (and therefore, the jet axis) to slew and rapidly change direction. This model not only explains recent observations showing successive generations of jet-lobes-bubbles in individual cool-core clusters that are offset from each other in the angular direction with respect to the cluster center, but also shows that AGN jets {\it can} heat the cluster core nearly isotropically on the gas cooling timescale. Our model {\it does} require that the SMBHs at the centers of cool-core clusters be spinning relatively slowly. Torques from individual misaligned discs are ineffective at tilting rapidly spinning black holes by more than a few degrees. Additionally, since SMBHs that host thin accretion discs will manifest as quasars, we predict that roughly 1--2 rich clusters within $z<0.5$ should have quasars at their centers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiatively Inefficient Accretion Flow Simulations with Cooling: Implications for Black Hole Transients

We study the effects of optically thin radiative cooling on the structure of radiatively inefficient accretion flows (RIAFs). The flow structure is geometrically thick, and independent of the gas density and cooling, if the cooling time is longer than the viscous timescale (i.e., $t_{\rm cool} \gtrsim t_{\rm visc}$). For higher densities, the gas can cool before it can accrete and forms the standard geometrically thin, optically thick Shakura-Sunyaev disk. For usual cooling processes (such as bremsstrahlung), we expect an inner hot flow and an outer thin disk. For a short cooling time the accretion flow separates into two phases: a radiatively inefficient hot coronal phase and a cold thin disk. We argue that there is an upper limit on the density of the hot corona corresponding to a critical value of $t_{\rm cool}/t_{\rm ff} (\sim 10-100)$, the ratio of the cooling time and the free-fall time. Based on our simulations, we have developed a model for transients observed in black hole X-ray binaries (XRBs). An XRB in a quiescent hot RIAF state can transition to a cold black-body dominated state because of an increase in the mass accretion rate. The transition from a thin disk to a RIAF happens because of mass exhaustion due to accretion; the transition happens when the cooling time becomes longer than the viscous time at inner radii. Since the viscous timescale for a geometrically thin disk is quite long, the high-soft state is expected to be long-lived. The different timescales in black hole transients correspond to different physical processes such as viscous evolution, cooling, and free-fall. Our model captures the overall features of observed state transitions in XRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Superbubble breakout and galactic winds from disk galaxies

We study the conditions for disk galaxies to produce superbubbles that can break out of the disk and produce a galactic wind. We argue that the threshold surface density of supernovae rate for seeding a wind depends on the ability of superbubble energetics to compensate for radiative cooling. We first adapt Kompaneets formalism for expanding bubbles in a stratified medium to the case of continuous energy injection and include the effects of radiative cooling in the shell. With the help of hydrodynamic simulations, we then study the evolution of superbubbles evolving in stratified disks with typical disk parameters. We identify two crucial energy injection rates that differ in their effects, the corresponding breakout ranging from being gentle to a vigorous one. (a) Superbubbles that break out of the disk with a Mach number of order 2-3 correspond to an energy injection rate of order 10^{-4} erg cm^{-2} s^{-1}, which is relevant for disk galaxies with synchrotron emitting gas in the extra-planar regions. (b) A larger energy injection threshold, of order 10^{-3} erg cm^{-2} s^{-1}, or equivalently, a star formation surface density of \sim 0.1 solar mass yr^{-1} kpc^{-2}, corresponds to superbubbles with a Mach number \sim 5-10. While the milder superbubbles can be produced by large OB associations, the latter kind requires super-starclusters. These derived conditions compare well with observations of disk galaxies with winds and the existence of multiphase halo gas. Furthermore, we find that contrary to the general belief that superbubbles fragment through Rayleigh-Taylor (RT) instability when they reach a vertical height of order the scale height, the superbubbles are first affected by thermal instability for typical disk parameters and that RT instability takes over when the shells reach a distance of approximately twice the scale height.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local 2D Particle-in-cell simulations of the collisionless MRI

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the rate of Coulomb collisions between particles is very small, making the disk evolve essentially as a collisionless plasma. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell (PIC) plasma simulations. In this initial study we focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field until the Alfvén velocity, v_A, is comparable to the speed of light, c (independent of the initial value of v_A/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the perpendicular pressure being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission -- from the radio to the gamma-rays -- of systems such as Sgr A*.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Structure of Hot Gas in Halos: Implications for the Lx-Tx Relation & Missing Baryons

We present one-dimensional models of the hot gas in dark-matter halos, which both predict the existence of cool cores and explain their structure. Our models are directly applicable to semi-analytic models (SAMs) of galaxy formation. We have previously argued that filaments of cold (~10^4 K) gas condense out of the intracluster medium (ICM) in hydrostatic and thermal equilibrium when the ratio of the thermal instability timescale to the free-fall time $t_{TI}/t_{ff}$ falls below 5-10. This criterion corresponds to an upper limit on the density of the ICM and motivates a model in which a density core forms wherever $t_{TI}/t_{ff} \lesssim 10$. Consistent with observations and numerical simulations, this model predicts larger and more tenuous cores for lower-mass halos---while the core density in a cluster may be as large as ~ 0.1 cm^{-3}, the core density in the Galactic halo should not exceed ~ 10^{-4} cm^{-3}. We can also explain the large densities in smaller mass halos (galactic `coronae') if we include the contribution of the central galaxy to the gravitational potential. Our models produce a favorable match to the observational X-ray luminosity-temperature (Lx-Tx) relation. For halo masses $\lesssim 10^{13}$ Msun the core size approaches the virial radius. Thus, most of the baryons in such halos cannot be in the hot ICM, but either in the form of stars or in the form of hot gas beyond the virial radius. Because of the smaller mass in the ICM and much larger mass available for star formation, the majority of the baryons in low mass halos ($\lesssim 10^{13}$ Msun) can be expelled beyond the virial radius due to supernova feedback. This can account for the baryons `missing' from low mass halos, such as the Galactic halo.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Effects of Anisotropic Viscosity on Turbulence and Heat Transport in the Intracluster Medium

In the intracluster medium (ICM) of galaxy clusters, heat and momentum are transported almost entirely along (but not across) magnetic field lines. We perform the first fully self-consistent Braginskii-MHD simulations of galaxy clusters including both of these effects. Specifically, we perform local and global simulations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability (HBI) and assess the effects of viscosity on their saturation and astrophysical implications. We find that viscosity has only a modest effect on the saturation of the MTI. As in previous calculations, we find that the MTI can generate nearly sonic turbulent velocities in the outer parts of galaxy clusters, although viscosity somewhat suppresses the magnetic field amplification. At smaller radii in cool-core clusters, viscosity can decrease the linear growth rates of the HBI. However, it has less of an effect on the HBI's nonlinear saturation, in part because three-dimensional interchange motions (magnetic flux tubes slipping past each other) are not damped by anisotropic viscosity. In global simulations of cool core clusters, we show that the HBI robustly inhibits radial thermal conduction and thus precipitates a cooling catastrophe. The effects of viscosity are, however, more important for higher entropy clusters. We argue that viscosity can contribute to the global transition of cluster cores from cool-core to non cool-core states: additional sources of intracluster turbulence, such as can be produced by AGN feedback or galactic wakes, suppress the HBI, heating the cluster core by thermal conduction; this makes the ICM more viscous, which slows the growth of the HBI, allowing further conductive heating of the cluster core and a transition to a non cool-core state.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal Instability in Gravitationally-Stratified Plasmas: Implications for Multi-Phase Structure in Clusters and Galaxy Halos

We study the interplay among cooling, heating, conduction, and magnetic fields in gravitationally stratified plasmas using simplified, plane-parallel numerical simulations. Since the physical heating mechanism remains uncertain in massive halos such as groups or clusters, we adopt a simple, observationally-motivated prescription which enforces global thermal equilibrium when averaged over large scales. The plasma remains susceptible to local thermal instability, however, and cooling drives an inward flow of material. In contrast to previous results, we argue that the thermal stability of the plasma is independent of its convective stability. We find that the ratio of the cooling timescale to the dynamical timescale t_cool/t_ff controls the saturation of the thermal instability: when t_cool/t_ff < 1, the plasma develops extended multi-phase structure, whereas when t_cool / t_ff > 1 it does not. (In a companion paper, we show that the criterion for thermal instability in a spherical potential is somewhat less stringent, t_cool / t_ff < 10.) When thermal conduction is anisotropic with respect to the magnetic field, the criterion for multi-phase structure is essentially independent of the thermal conductivity of the plasma. Our criterion for local thermal instability to produce multi-phase structure is an extension of the cold vs. hot accretion modes in galaxy formation that applies at all radii in hot halos, not just to the virial shock. We show that this criterion is consistent with data on multi-phase gas in the ACCEPT sample of clusters; in addition, when t_cool / t_ff > 1, the net cooling rate to low temperatures and the mass flux to small radii are suppressed enough relative to models without heating to be qualitatively consistent with star formation rates and x-ray line emission in groups and clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Fast Semi-implicit Method for Anisotropic Diffusion

Simple finite differencing of the anisotropic diffusion equation, where diffusion is only along a given direction, does not ensure that the numerically calculated heat fluxes are in the correct direction. This can lead to negative temperatures for the anisotropic thermal diffusion equation. In a previous paper we proposed a monotonicity-preserving explicit method which uses limiters (analogous to those used in the solution of hyperbolic equations) to interpolate the temperature gradients at cell faces. However, being explicit, this method was limited by a restrictive Courant-Friedrichs-Lewy (CFL) stability timestep. Here we propose a fast, conservative, directionally-split, semi-implicit method which is second order accurate in space, is stable for large timesteps, and is easy to implement in parallel. Although not strictly monotonicity-preserving, our method gives only small amplitude temperature oscillations at large temperature gradients, and the oscillations are damped in time. With numerical experiments we show that our semi-implicit method can achieve large speed-ups compared to the explicit method, without seriously violating the monotonicity constraint. This method can also be applied to isotropic diffusion, both on regular and distorted meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Implementation of Streaming Down the Gradient: Application to Fluid Modeling of Cosmic Rays and Saturated Conduction

The equation governing the streaming of a quantity down its gradient superficially looks similar to the simple constant velocity advection equation. In fact, it is the same as an advection equation if there are no local extrema in the computational domain or at the boundary. However, in general when there are local extrema in the computational domain it is a non-trivial nonlinear equation. The standard upwind time evolution with a CFL-limited time step results in spurious oscillations at the grid scale. These oscillations, which originate at the extrema, propagate throughout the computational domain and are undamped even at late times. These oscillations arise because of unphysically large fluxes leaving (entering) the maxima (minima) with the standard CFL-limited explicit methods. Regularization of the equation shows that it is diffusive at the extrema; because of this, an explicit method for the regularized equation with $Δt \propto Δx^2$ behaves fine. We show that the implicit methods show stable and converging results with $Δt \propto Δx$; however, surprisingly, even implicit methods are not stable with large enough timesteps. In addition to these subtleties in the numerical implementation, the solutions to the streaming equation are quite novel: non-differentiable solutions emerge from initially smooth profiles; the solutions show transport over large length scales, e.g., in form of tails. The fluid model for cosmic rays interacting with a thermal plasma (valid at space scales much larger than the cosmic ray Larmor radius) is similar to the equation for streaming of a quantity down its gradient, so our method will find applications in fluid modeling of cosmic rays.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal Instability & the Feedback Regulation of Hot Halos in Clusters, Groups, and Galaxies

Observations of clusters and groups imply that such halos are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally-motivated heating prescriptions, we show that local thermal instability (TI) can produce a multi-phase medium---with ~ 10000 K cold filaments condensing out of the hot ICM---only when the ratio of the TI timescale in the hot plasma (t_{TI}) to the free-fall timescale (t_{ff}) satisfies t_{TI}/t_{ff} <~ 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling, and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of ~ 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling-flow is ~ 0.001 for clusters and decreases for lower mass halos; supernova heating may be energetically sufficient to balance cooling in galactic halos. We further argue that the ICM self-adjusts so that t_{TI}/t_{ff} >~ 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives t_{TI}/t_{ff} >~ 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminosity-temperature (L_X-T_X) relation. The high-velocity clouds observed in the Galactic halo can also be due to local TI producing multi-phase gas close to the virial radius.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Turbulent Pressure Support in the Outer Parts of Galaxy Clusters

We use three-dimensional MHD simulations with anisotropic thermal conduction to study turbulence due to the magnetothermal instability (MTI) in the intracluster medium (ICM) of galaxy clusters. The MTI grows on timescales of ~1 Gyr and is capable of driving vigorous, sustained turbulence in the outer parts of galaxy clusters if the temperature gradient is maintained in spite of the rapid thermal conduction. If this is the case, turbulence due to the MTI can provide up to 5-30% of the pressure support beyond r_500 in galaxy clusters, an effect that is strongest for hot, massive clusters. The turbulence driven by the MTI is generally additive to other sources of turbulence in the ICM, such as that produced by structure formation. This new source of non-thermal pressure support reduces the observed Sunyaev-Zel'dovich (SZ) signal and X-ray pressure gradient for a given cluster mass and introduces a cluster mass and temperature gradient-dependent bias in SZ and X-ray mass estimates of clusters. This additional physics may also need to be taken into account when estimating the matter power spectrum normalization, sigma-8, through simulation templates from the observed amplitude of the SZ power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Can conduction induce convection? The non-linear saturation of buoyancy instabilities in dilute plasmas

We study the effects of anisotropic thermal conduction on low-collisionality, astrophysical plasmas using two and three-dimensional magnetohydrodynamic simulations. For weak magnetic fields, dilute plasmas are buoyantly unstable for either sign of the temperature gradient: the heat-flux-driven buoyancy instability (HBI) operates when the temperature increases with radius while the magnetothermal instability (MTI) operates in the opposite limit. In contrast to previous results, we show that, in the presence of a sustained temperature gradient, the MTI drives strong turbulence and operates as an efficient magnetic dynamo (akin to standard, adiabatic convection). Together, the turbulent and magnetic energies contribute up to ~10% of the pressure support in the plasma. In addition, the MTI drives a large convective heat flux, ~1.5% of rho c_s^3. These findings are robust even in the presence of an external source of strong turbulence. Our results on the nonlinear saturation of the HBI are consistent with previous studies but we explain physically why the HBI saturates quiescently by re-orienting the magnetic field (suppressing the conductive heat flux through the plasma), while the MTI saturates by generating sustained turbulence. We also systematically study how an external source of turbulence affects the saturation of the HBI: such turbulence can disrupt the HBI only on scales where the shearing rate of the turbulence is faster than the growth rate of the HBI. In particular, our results provide a simple mapping between the level of turbulence in a plasma and the effective isotropic thermal conductivity. We discuss the astrophysical implications of these findings, with a particular focus on the intracluster medium of galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermal Instability with Anisotropic Thermal Conduction and Adiabatic Cosmic Rays: Implications for Cold Filaments in Galaxy Clusters

Observations of the cores of nearby galaxy clusters show H$α$ and molecular emission line filaments. We argue that these are the result of {\em local} thermal instability in a {\em globally} stable galaxy cluster core. We present local, high resolution, two-dimensional magnetohydrodynamic simulations of thermal instability for conditions appropriate to the intracluster medium (ICM); the simulations include thermal conduction along magnetic field lines and adiabatic cosmic rays. Thermal conduction suppresses thermal instability along magnetic field lines on scales smaller than the Field length ($\gtrsim$10 kpc for the hot, diffuse ICM). We show that the Field length in the cold medium must be resolved both along and perpendicular to the magnetic field in order to obtain numerically converged results. Because of negligible conduction perpendicular to the magnetic field, thermal instability leads to fine scale structure in the perpendicular direction. Filaments of cold gas along magnetic field lines are thus a natural consequence of thermal instability with anisotropic thermal conduction. Nonlinearly, filaments of cold ($\sim 10^4$ K) gas should have lengths (along the magnetic field) comparable to the Field length in the cold medium $\sim 10^{-4}$ pc! Observations show, however, that the atomic filaments in clusters are far more extended, $\sim 10$ kpc. Cosmic ray pressure support (or a small scale turbulent magnetic pressure) may resolve this discrepancy: even a small cosmic ray pressure in the diffuse ICM, $\sim 10^{-4}$ of the thermal pressure, can be adiabatically compressed to provide significant pressure support in cold filaments. This is qualitatively consistent with the large population of cosmic rays invoked to explain the atomic and molecular line ratios observed in filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Time Dependent Models of Flares from Sagittarius A*

The emission from Sgr A*, the supermassive black hole in the Galactic Center, shows order of magnitude variability ("flares") a few times a day that is particularly prominent in the near-infrared (NIR) and X-rays. We present a time-dependent model for these flares motivated by the hypothesis that dissipation of magnetic energy powers the flares. We show that episodic magnetic reconnection can occur near the last stable circular orbit in time-dependent magnetohydrodynamic simulations of black hole accretion - the timescales and energetics of these events are broadly consistent with the flares from Sgr A*. Motivated by these results, we present a spatially one-zone time-dependent model for the electron distribution function in flares, including energy loss due to synchrotron cooling and adiabatic expansion. Synchrotron emission from transiently accelerated particles can explain the NIR/X-ray lightcurves and spectra of a luminous flare observed 4 April 2007. A significant decrease in the magnetic field strength during the flare (coincident with the electron acceleration) is required to explain the simultaneity and symmetry of the simultaneous lightcurves. Our models predict that the NIR and X-ray spectral indices differ by 0.5 and that there is only modest variation in the spectral index during flares. We also explore implications of this model for longer wavelength (radio-submm) emission seemingly associated with X-ray and NIR flares; we argue that a few hour decrease in the submm emission is a more generic consequence of large-scale magnetic reconnection than delayed radio emission from adiabatic expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulence in Galaxy Cluster Cores: a Key to Cluster Bimodality?

We study the effects of externally imposed turbulence on the thermal properties of galaxy cluster cores, using three-dimensional numerical simulations including magnetic fields, anisotropic thermal conduction, and radiative cooling. The imposed "stirring" crudely approximates the effects of galactic wakes, waves generated by galaxies moving through the intracluster medium (ICM), and/or turbulence produced by a central active galactic nucleus. The simulated clusters exhibit a strong bimodality. Modest levels of turbulence, ~100 km/s (~10% of the sound speed), suppress the heat-flux-driven buoyancy instability (HBI), resulting in an isotropically tangled magnetic field and a quasi-stable, high entropy, thermal equilibrium with no cooling catastrophe. Thermal conduction dominates the heating of the cluster core, but turbulent mixing is critical because it suppresses the HBI and (to a lesser extent) the thermal instability. Lower levels of turbulent mixing (approximately less than 100 km/s) are insufficient to suppress the HBI, rapidly leading to a thermal runaway and a cool-core cluster. Remarkably, then, small fluctuations in the level of turbulence in galaxy cluster cores can initiate transitions between cool-core (low entropy) and non cool-core (high entropy) states.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Spherical Accretion with Anisotropic Thermal Conduction

We study the effects of anisotropic thermal conduction on magnetized spherical accretion flows using global axisymmetric MHD simulations. In low collisionality plasmas, the Bondi spherical accretion solution is unstable to the magnetothermal instability (MTI). The MTI grows rapidly at large radii where the inflow is subsonic. For a weak initial field, the MTI saturates by creating a primarily radial magnetic field, i.e., by aligning the field lines with the background temperature gradient. The saturation is quasilinear in the sense that the magnetic field is amplified by a factor of $\sim 10-30$ independent of the initial field strength (for weak fields). In the saturated state, the conductive heat flux is much larger than the convective heat flux, and is comparable to the field-free (Spitzer) value (since the field lines are largely radial). The MTI by itself does not appreciably change the accretion rate $\dot M$ relative to the Bondi rate $\dot M_B$. However, the radial field lines created by the MTI are amplified by flux freezing as the plasma flows in to small radii. Oppositely directed field lines are brought together by the converging inflow, leading to significant resistive heating. When the magnetic energy density is comparable to the gravitational potential energy density, the plasma is heated to roughly the virial temperature; the mean inflow is highly subsonic; most of the energy released by accretion is transported to large radii by thermal conduction; and the accretion rate $\dot M \ll \dot M_B$. The predominantly radial magnetic field created by the MTI at large radii in spherical accretion flows may account for the stable Faraday rotation measure towards Sgr A* in the Galactic Center.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Electron Heating in Hot Accretion Flows

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational instability in a collisionless plasma show that angular momentum transport by pressure anisotropy ($p_\perp \ne p_\parallel$, where the directions are defined with respect to the local magnetic field) is comparable to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a large-scale viscosity, arises because of adiabatic invariants related to $p_\perp$ and $p_\parallel$ in a fluctuating magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show that $\sim 50$ % of the gravitational potential energy is directly converted into heat at large scales by the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and magnetic energy). We show that electrons receive a significant fraction ($\sim [T_e/T_i]^{1/2}$) of this dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow is greater than 0.5% for $\dot{M} \gtrsim 10^{-4} \dot{M}_{Edd}$. Thus a low accretion rate, rather than just a low radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of $\approx 3 \times 10^{-8} M_\odot {\rm yr^{-1}}$ and an electron temperature of $\approx 3 \times 10^{10}$ K at $\approx 10$ Schwarzschild radii; the latter is consistent with the brightness temperature inferred from VLBI observations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Faraday Rotation in Global Accretion Disk Simulations: Implications for Sgr A*

These Faraday rotation calculations of hot, thick accretion flows are motivated by the measured steady rotation measure (RM) of $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A*. In our numerical simulations, the quasi-steady state structure of the accretion flow, and the RM it produces, depends on the initial magnetic field. In spite of this dependence, we can draw several robust conclusions about Faraday rotation produced by geometrically thick accretion disks: i) the time averaged RM does not depend that sensitively on the viewing angle, but the stability of the RM can. Equatorial viewing angles show significant variability in RM (including sign reversals), while polar viewing angles are relatively stable if there is a large scale magnetic field threading the disk at large radii. ii) Most of the RM is produced at small radii for polar viewing angles while all radii contribute significantly near the midplane of the disk. Our simulations confirm previous analytic arguments that the accretion rate onto Sgr A* must satisfy $\dot M_{\rm in} \ll \dot M_{\rm Bondi} \sim 10^{-5} \mpy$ in order to not over-produce the measured RM. We argue that the steady RM $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A* has two plausible explanations: 1) it is produced at $\sim 100$ Schwarzschild radii, requires $\dot{M}_{\rm in} \approx 3 \times 10^{-8} M_\odot$ yr$^{-1}$, and we view the flow at an angle of $\sim 30^\circ$ relative to the rotation axis of the disk; in our simulations, the variation in RM across a finite-sized source is sufficient to depolarize the emission below $\approx$ 100 GHz, consistent with observations. 2) Alternatively, the RM may be produced in the relatively spherical inflowing plasma near the circularization radius at $\sim 10^3-10^4$ Schwarzschild radii.

0 citations0 reviewsNo code linkedOpen access