Normalized to: Sharma, P.
[2]
oai:arXiv.org:2006.03291 [pdf] - 2115393
Broad-band spectral analysis of LMXB XTE J1710-281 with Suzaku
Submitted: 2020-06-05
This work presents the broad-band time-averaged spectral analysis of neutron
star low-mass X-ray binary, XTE J1710-281 by using the Suzaku archival data.
The source was in a hard or an intermediate spectral state during this
observation. This is the first time that a detailed spectral analysis of the
persistent emission spectra of XTE J1710-281 has been done up to 30 keV with
improved constraints on its spectral parameters. By simultaneously fitting the
XIS (0.6-9.0 keV) and the HXD-PIN (15.0-30.0 keV) data, we have modelled the
persistent spectrum of the source with models comprising a soft component from
accretion disc and/or neutron star surface/boundary layer and a hard
Comptonizing component. The 0.6-30 keV continuum with neutral absorber can be
described by a multi-colour disc blackbody with an inner disc temperature of
$kT_{\rm disc} = 0.28$ keV, which is significantly Comptonized by the hot
electron cloud with electron temperature of $kT_{\rm e} \approx 5$ keV and
described by photon index $\Gamma = 1.86$. A more complex three-component model
comprising a multi-colour disc blackbody $\approx 0.30$ keV, single temperature
blackbody $\approx 0.65$ keV and Comptonization from the disc, partially
absorbed (about 38 per cent) by an ionized absorber (log($\xi$) $\approx$ 4)
describes the broad-band spectrum equally well.
[3]
oai:arXiv.org:2005.09654 [pdf] - 2098838
Resolving small-scale cold circumgalactic gas in TNG50
Submitted: 2020-05-19
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 contrast 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.
[4]
oai:arXiv.org:1909.12888 [pdf] - 2093117
Cool-Core Cycles and Phoenix
Submitted: 2019-09-27, last modified: 2020-04-28
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}.
[6]
oai:arXiv.org:1912.04916 [pdf] - 2084113
Characterising the Dynamo in a Radiatively Inefficient Accretion Flow
Submitted: 2019-12-10, last modified: 2020-04-27
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 $\alpha$ and
$\gamma$ 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
$\gamma$ 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 $\alpha_{\phi \phi}$, 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 $\alpha$-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 $\alpha-\Omega$ type.
[8]
oai:arXiv.org:2001.06494 [pdf] - 2070111
Turbulence in stratified atmospheres: implications for the intracluster
medium
Submitted: 2020-01-17, last modified: 2020-03-25
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
($\sigma_s$), where $s=\ln(\rho/\left<\rho(z)\right>)$, increases with
$\mathrm{Ri}$. For weakly stratified subsonic turbulence
($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new
$\sigma_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation,
$\sigma_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(\rho_k/\left<\rho\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}$.
[17]
oai:arXiv.org:1911.07870 [pdf] - 2026376
Radio halos of star forming galaxies
Submitted: 2019-11-18
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
\Sigma_{\rm SFR}^{0.3}$ which is consistent with earlier observations and
analytical estimates.
[18]
oai:arXiv.org:1910.10168 [pdf] - 2046295
Realistic modeling of wind and supernovae shocks in star clusters:
addressing ${\rm ^{22}Ne/^{20}Ne}$ and other problems in Galactic cosmic rays
Submitted: 2019-10-22
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
$\gamma$-ray bright ($L_\gamma \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.
[19]
oai:arXiv.org:1910.06339 [pdf] - 1979841
On modelling CC85 wind in an expanding local box
Submitted: 2019-10-14
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]
[20]
oai:arXiv.org:1811.12907 [pdf] - 1973647
GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers
Observed by LIGO and Virgo during the First and Second Observing Runs
The LIGO Scientific Collaboration;
the Virgo Collaboration;
Abbott, B. P.;
Abbott, R.;
Abbott, T. D.;
Abraham, S.;
Acernese, F.;
Ackley, K.;
Adams, C.;
Adhikari, R. X.;
Adya, V. B.;
Affeldt, C.;
Agathos, M.;
Agatsuma, K.;
Aggarwal, N.;
Aguiar, O. D.;
Aiello, L.;
Ain, A.;
Ajith, P.;
Allen, G.;
Allocca, A.;
Aloy, M. A.;
Altin, P. A.;
Amato, A.;
Ananyeva, A.;
Anderson, S. B.;
Anderson, W. G.;
Angelova, S. V.;
Antier, S.;
Appert, S.;
Arai, K.;
Araya, M. C.;
Areeda, J. S.;
Arène, M.;
Arnaud, N.;
Arun, K. G.;
Ascenzi, S.;
Ashton, G.;
Aston, S. M.;
Astone, P.;
Aubin, F.;
Aufmuth, P.;
AultONeal, K.;
Austin, C.;
Avendano, V.;
Avila-Alvarez, A.;
Babak, S.;
Bacon, P.;
Badaracco, F.;
Bader, M. K. M.;
Bae, S.;
Baker, P. T.;
Baldaccini, F.;
Ballardin, G.;
Ballmer, S. W.;
Banagiri, S.;
Barayoga, J. C.;
Barclay, S. E.;
Barish, B. C.;
Barker, D.;
Barkett, K.;
Barnum, S.;
Barone, F.;
Barr, B.;
Barsotti, L.;
Barsuglia, M.;
Barta, D.;
Bartlett, J.;
Bartos, I.;
Bassiri, R.;
Basti, A.;
Bawaj, M.;
Bayley, J. C.;
Bazzan, M.;
Bécsy, B.;
Bejger, M.;
Belahcene, I.;
Bell, A. S.;
Beniwal, D.;
Berger, B. K.;
Bergmann, G.;
Bernuzzi, S.;
Bero, J. J.;
Berry, C. P. L.;
Bersanetti, D.;
Bertolini, A.;
Betzwieser, J.;
Bhandare, R.;
Bidler, J.;
Bilenko, I. A.;
Bilgili, S. A.;
Billingsley, G.;
Birch, J.;
Birney, R.;
Birnholtz, O.;
Biscans, S.;
Biscoveanu, S.;
Bisht, A.;
Bitossi, M.;
Bizouard, M. A.;
Blackburn, J. K.;
Blackman, J.;
Blair, C. D.;
Blair, D. G.;
Blair, R. M.;
Bloemen, S.;
Bode, N.;
Boer, M.;
Boetzel, Y.;
Bogaert, G.;
Bondu, F.;
Bonilla, E.;
Bonnand, R.;
Booker, P.;
Boom, B. A.;
Booth, C. D.;
Bork, R.;
Boschi, V.;
Bose, S.;
Bossie, K.;
Bossilkov, V.;
Bosveld, J.;
Bouffanais, Y.;
Bozzi, A.;
Bradaschia, C.;
Brady, P. R.;
Bramley, A.;
Branchesi, M.;
Brau, J. E.;
Briant, T.;
Briggs, J. H.;
Brighenti, F.;
Brillet, A.;
Brinkmann, M.;
Brisson, V.;
Brockill, P.;
Brooks, A. F.;
Brown, D. D.;
Brunett, S.;
Buikema, A.;
Bulik, T.;
Bulten, H. J.;
Buonanno, A.;
Buskulic, D.;
Rosell, M. J. Bustamante;
Buy, C.;
Byer, R. L.;
Cabero, M.;
Cadonati, L.;
Cagnoli, G.;
Cahillane, C.;
Bustillo, J. Calderón;
Callister, T. A.;
Calloni, E.;
Camp, J. B.;
Campbell, W. A.;
Canepa, M.;
Cannon, K. C.;
Cao, H.;
Cao, J.;
Capocasa, E.;
Carbognani, F.;
Caride, S.;
Carney, M. F.;
Carullo, G.;
Diaz, J. Casanueva;
Casentini, C.;
Caudill, S.;
Cavaglià, M.;
Cavalier, F.;
Cavalieri, R.;
Cella, G.;
Cerdá-Durán, P.;
Cerretani, G.;
Cesarini, E.;
Chaibi, O.;
Chakravarti, K.;
Chamberlin, S. J.;
Chan, M.;
Chao, S.;
Charlton, P.;
Chase, E. A.;
Chassande-Mottin, E.;
Chatterjee, D.;
Chaturvedi, M.;
Chatziioannou, K.;
Cheeseboro, B. D.;
Chen, H. Y.;
Chen, X.;
Chen, Y.;
Cheng, H. -P.;
Cheong, C. K.;
Chia, H. Y.;
Chincarini, A.;
Chiummo, A.;
Cho, G.;
Cho, H. S.;
Cho, M.;
Christensen, N.;
Chu, Q.;
Chua, S.;
Chung, K. W.;
Chung, S.;
Ciani, G.;
Ciobanu, A. A.;
Ciolfi, R.;
Cipriano, F.;
Cirone, A.;
Clara, F.;
Clark, J. A.;
Clearwater, P.;
Cleva, F.;
Cocchieri, C.;
Coccia, E.;
Cohadon, P. -F.;
Cohen, D.;
Colgan, R.;
Colleoni, M.;
Collette, C. G.;
Collins, C.;
Cominsky, L. R.;
Constancio, M.;
Conti, L.;
Cooper, S. J.;
Corban, P.;
Corbitt, T. R.;
Cordero-Carrión, I.;
Corley, K. R.;
Cornish, N.;
Corsi, A.;
Cortese, S.;
Costa, C. A.;
Cotesta, R.;
Coughlin, M. W.;
Coughlin, S. B.;
Coulon, J. -P.;
Countryman, S. T.;
Couvares, P.;
Covas, P. B.;
Cowan, E. E.;
Coward, D. M.;
Cowart, M. J.;
Coyne, D. C.;
Coyne, R.;
Creighton, J. D. E.;
Creighton, T. D.;
Cripe, J.;
Croquette, M.;
Crowder, S. G.;
Cullen, T. J.;
Cumming, A.;
Cunningham, L.;
Cuoco, E.;
Canton, T. Dal;
Dálya, G.;
Danilishin, S. L.;
D'Antonio, S.;
Danzmann, K.;
Dasgupta, A.;
Costa, C. F. Da Silva;
Datrier, L. E. H.;
Dattilo, V.;
Dave, I.;
Davier, M.;
Davis, D.;
Daw, E. J.;
DeBra, D.;
Deenadayalan, M.;
Degallaix, J.;
De Laurentis, M.;
Deléglise, S.;
Del Pozzo, W.;
DeMarchi, L. M.;
Demos, N.;
Dent, T.;
De Pietri, R.;
Derby, J.;
De Rosa, R.;
De Rossi, C.;
DeSalvo, R.;
de Varona, O.;
Dhurandhar, S.;
Díaz, M. C.;
Dietrich, T.;
Di Fiore, L.;
Di Giovanni, M.;
Di Girolamo, T.;
Di Lieto, A.;
Ding, B.;
Di Pace, S.;
Di Palma, I.;
Di Renzo, F.;
Dmitriev, A.;
Doctor, Z.;
Donovan, F.;
Dooley, K. L.;
Doravari, S.;
Dorrington, I.;
Downes, T. P.;
Drago, M.;
Driggers, J. C.;
Du, Z.;
Ducoin, J. -G.;
Dupej, P.;
Dwyer, S. E.;
Easter, P. J.;
Edo, T. B.;
Edwards, M. C.;
Effler, A.;
Ehrens, P.;
Eichholz, J.;
Eikenberry, S. S.;
Eisenmann, M.;
Eisenstein, R. A.;
Essick, R. C.;
Estelles, H.;
Estevez, D.;
Etienne, Z. B.;
Etzel, T.;
Evans, M.;
Evans, T. M.;
Fafone, V.;
Fair, H.;
Fairhurst, S.;
Fan, X.;
Farinon, S.;
Farr, B.;
Farr, W. M.;
Fauchon-Jones, E. J.;
Favata, M.;
Fays, M.;
Fazio, M.;
Fee, C.;
Feicht, J.;
Fejer, M. M.;
Feng, F.;
Fernandez-Galiana, A.;
Ferrante, I.;
Ferreira, E. C.;
Ferreira, T. A.;
Ferrini, F.;
Fidecaro, F.;
Fiori, I.;
Fiorucci, D.;
Fishbach, M.;
Fisher, R. P.;
Fishner, J. M.;
Fitz-Axen, M.;
Flaminio, R.;
Fletcher, M.;
Flynn, E.;
Fong, H.;
Font, J. A.;
Forsyth, P. W. F.;
Fournier, J. -D.;
Frasca, S.;
Frasconi, F.;
Frei, Z.;
Freise, A.;
Frey, R.;
Frey, V.;
Fritschel, P.;
Frolov, V. V.;
Fulda, P.;
Fyffe, M.;
Gabbard, H. A.;
Gadre, B. U.;
Gaebel, S. M.;
Gair, J. R.;
Gammaitoni, L.;
Ganija, M. R.;
Gaonkar, S. G.;
Garcia, A.;
García-Quirós, C.;
Garufi, F.;
Gateley, B.;
Gaudio, S.;
Gaur, G.;
Gayathri, V.;
Gemme, G.;
Genin, E.;
Gennai, A.;
George, D.;
George, J.;
Gergely, L.;
Germain, V.;
Ghonge, S.;
Ghosh, Abhirup;
Ghosh, Archisman;
Ghosh, S.;
Giacomazzo, B.;
Giaime, J. A.;
Giardina, K. D.;
Giazotto, A.;
Gill, K.;
Giordano, G.;
Glover, L.;
Godwin, P.;
Goetz, E.;
Goetz, R.;
Goncharov, B.;
González, G.;
Castro, J. M. Gonzalez;
Gopakumar, A.;
Gorodetsky, M. L.;
Gossan, S. E.;
Gosselin, M.;
Gouaty, R.;
Grado, A.;
Graef, C.;
Granata, M.;
Grant, A.;
Gras, S.;
Grassia, P.;
Gray, C.;
Gray, R.;
Greco, G.;
Green, A. C.;
Green, R.;
Gretarsson, E. M.;
Groot, P.;
Grote, H.;
Grunewald, S.;
Gruning, P.;
Guidi, G. M.;
Gulati, H. K.;
Guo, Y.;
Gupta, A.;
Gupta, M. K.;
Gustafson, E. K.;
Gustafson, R.;
Haegel, L.;
Halim, O.;
Hall, B. R.;
Hall, E. D.;
Hamilton, E. Z.;
Hammond, G.;
Haney, M.;
Hanke, M. M.;
Hanks, J.;
Hanna, C.;
Hannam, M. D.;
Hannuksela, O. A.;
Hanson, J.;
Hardwick, T.;
Haris, K.;
Harms, J.;
Harry, G. M.;
Harry, I. W.;
Haster, C. -J.;
Haughian, K.;
Hayes, F. J.;
Healy, J.;
Heidmann, A.;
Heintze, M. C.;
Heitmann, H.;
Hello, P.;
Hemming, G.;
Hendry, M.;
Heng, I. S.;
Hennig, J.;
Heptonstall, A. W.;
Vivanco, Francisco Hernandez;
Heurs, M.;
Hild, S.;
Hinderer, T.;
Hoak, D.;
Hochheim, S.;
Hofman, D.;
Holgado, A. M.;
Holland, N. A.;
Holt, K.;
Holz, D. E.;
Hopkins, P.;
Horst, C.;
Hough, J.;
Howell, E. J.;
Hoy, C. G.;
Hreibi, A.;
Huang, Y.;
Huerta, E. A.;
Huet, D.;
Hughey, B.;
Hulko, M.;
Husa, S.;
Huttner, S. H.;
Huynh-Dinh, T.;
Idzkowski, B.;
Iess, A.;
Ingram, C.;
Inta, R.;
Intini, G.;
Irwin, B.;
Isa, H. N.;
Isac, J. -M.;
Isi, M.;
Iyer, B. R.;
Izumi, K.;
Jacqmin, T.;
Jadhav, S. J.;
Jani, K.;
Janthalur, N. N.;
Jaranowski, P.;
Jenkins, A. C.;
Jiang, J.;
Johnson, D. S.;
Johnson-McDaniel, N. K.;
Jones, A. W.;
Jones, D. I.;
Jones, R.;
Jonker, R. J. G.;
Ju, L.;
Junker, J.;
Kalaghatgi, C. V.;
Kalogera, V.;
Kamai, B.;
Kandhasamy, S.;
Kang, G.;
Kanner, J. B.;
Kapadia, S. J.;
Karki, S.;
Karvinen, K. S.;
Kashyap, R.;
Kasprzack, M.;
Katsanevas, S.;
Katsavounidis, E.;
Katzman, W.;
Kaufer, S.;
Kawabe, K.;
Keerthana, N. V.;
Kéfélian, F.;
Keitel, D.;
Kennedy, R.;
Key, J. S.;
Khalili, F. Y.;
Khan, H.;
Khan, I.;
Khan, S.;
Khan, Z.;
Khazanov, E. A.;
Khursheed, M.;
Kijbunchoo, N.;
Kim, Chunglee;
Kim, J. C.;
Kim, K.;
Kim, W.;
Kim, W. S.;
Kim, Y. -M.;
Kimball, C.;
King, E. J.;
King, P. J.;
Kinley-Hanlon, M.;
Kirchhoff, R.;
Kissel, J. S.;
Kleybolte, L.;
Klika, J. H.;
Klimenko, S.;
Knowles, T. D.;
Koch, P.;
Koehlenbeck, S. M.;
Koekoek, G.;
Koley, S.;
Kondrashov, V.;
Kontos, A.;
Koper, N.;
Korobko, M.;
Korth, W. Z.;
Kowalska, I.;
Kozak, D. B.;
Kringel, V.;
Krishnendu, N.;
Królak, A.;
Kuehn, G.;
Kumar, A.;
Kumar, P.;
Kumar, R.;
Kumar, S.;
Kuo, L.;
Kutynia, A.;
Kwang, S.;
Lackey, B. D.;
Lai, K. H.;
Lam, T. L.;
Landry, M.;
Lane, B. B.;
Lang, R. N.;
Lange, J.;
Lantz, B.;
Lanza, R. K.;
Lartaux-Vollard, A.;
Lasky, P. D.;
Laxen, M.;
Lazzarini, A.;
Lazzaro, C.;
Leaci, P.;
Leavey, S.;
Lecoeuche, Y. K.;
Lee, C. H.;
Lee, H. K.;
Lee, H. M.;
Lee, H. W.;
Lee, J.;
Lee, K.;
Lehmann, J.;
Lenon, A.;
Leroy, N.;
Letendre, N.;
Levin, Y.;
Li, J.;
Li, K. J. L.;
Li, T. G. F.;
Li, X.;
Lin, F.;
Linde, F.;
Linker, S. D.;
Littenberg, T. B.;
Liu, J.;
Liu, X.;
Lo, R. K. L.;
Lockerbie, N. A.;
London, L. T.;
Longo, A.;
Lorenzini, M.;
Loriette, V.;
Lormand, M.;
Losurdo, G.;
Lough, J. D.;
Lousto, C. O.;
Lovelace, G.;
Lower, M. E.;
Lück, H.;
Lumaca, D.;
Lundgren, A. P.;
Lynch, R.;
Ma, Y.;
Macas, R.;
Macfoy, S.;
MacInnis, M.;
Macleod, D. M.;
Macquet, A.;
Magaña-Sandoval, F.;
Zertuche, L. Magaña;
Magee, R. M.;
Majorana, E.;
Maksimovic, I.;
Malik, A.;
Man, N.;
Mandic, V.;
Mangano, V.;
Mansell, G. L.;
Manske, M.;
Mantovani, M.;
Marchesoni, F.;
Marion, F.;
Márka, S.;
Márka, Z.;
Markakis, C.;
Markosyan, A. S.;
Markowitz, A.;
Maros, E.;
Marquina, A.;
Marsat, S.;
Martelli, F.;
Martin, I. W.;
Martin, R. M.;
Martynov, D. V.;
Mason, K.;
Massera, E.;
Masserot, A.;
Massinger, T. J.;
Masso-Reid, M.;
Mastrogiovanni, S.;
Matas, A.;
Matichard, F.;
Matone, L.;
Mavalvala, N.;
Mazumder, N.;
McCann, J. J.;
McCarthy, R.;
McClelland, D. E.;
McCormick, S.;
McCuller, L.;
McGuire, S. C.;
McIver, J.;
McManus, D. J.;
McRae, T.;
McWilliams, S. T.;
Meacher, D.;
Meadors, G. D.;
Mehmet, M.;
Mehta, A. K.;
Meidam, J.;
Melatos, A.;
Mendell, G.;
Mercer, R. A.;
Mereni, L.;
Merilh, E. L.;
Merzougui, M.;
Meshkov, S.;
Messenger, C.;
Messick, C.;
Metzdorff, R.;
Meyers, P. M.;
Miao, H.;
Michel, C.;
Middleton, H.;
Mikhailov, E. E.;
Milano, L.;
Miller, A. L.;
Miller, A.;
Millhouse, M.;
Mills, J. C.;
Milovich-Goff, M. C.;
Minazzoli, O.;
Minenkov, Y.;
Mishkin, A.;
Mishra, C.;
Mistry, T.;
Mitra, S.;
Mitrofanov, V. P.;
Mitselmakher, G.;
Mittleman, R.;
Mo, G.;
Moffa, D.;
Mogushi, K.;
Mohapatra, S. R. P.;
Montani, M.;
Moore, C. J.;
Moraru, D.;
Moreno, G.;
Morisaki, S.;
Mours, B.;
Mow-Lowry, C. M.;
Mukherjee, Arunava;
Mukherjee, D.;
Mukherjee, S.;
Mukund, N.;
Mullavey, A.;
Munch, J.;
Muñiz, E. A.;
Muratore, M.;
Murray, P. G.;
Nagar, A.;
Nardecchia, I.;
Naticchioni, L.;
Nayak, R. K.;
Neilson, J.;
Nelemans, G.;
Nelson, T. J. N.;
Nery, M.;
Neunzert, A.;
Ng, K. Y.;
Ng, S.;
Nguyen, P.;
Nichols, D.;
Nielsen, A. B.;
Nissanke, S.;
Nitz, A.;
Nocera, F.;
North, C.;
Nuttall, L. K.;
Obergaulinger, M.;
Oberling, J.;
O'Brien, B. D.;
O'Dea, G. D.;
Ogin, G. H.;
Oh, J. J.;
Oh, S. H.;
Ohme, F.;
Ohta, H.;
Okada, M. A.;
Oliver, M.;
Oppermann, P.;
Oram, Richard J.;
O'Reilly, B.;
Ormiston, R. G.;
Ortega, L. F.;
O'Shaughnessy, R.;
Ossokine, S.;
Ottaway, D. J.;
Overmier, H.;
Owen, B. J.;
Pace, A. E.;
Pagano, G.;
Page, M. A.;
Pai, A.;
Pai, S. A.;
Palamos, J. R.;
Palashov, O.;
Palomba, C.;
Pal-Singh, A.;
Pan, Huang-Wei;
Pang, B.;
Pang, P. T. H.;
Pankow, C.;
Pannarale, F.;
Pant, B. C.;
Paoletti, F.;
Paoli, A.;
Papa, M. A.;
Parida, A.;
Parker, W.;
Pascucci, D.;
Pasqualetti, A.;
Passaquieti, R.;
Passuello, D.;
Patil, M.;
Patricelli, B.;
Pearlstone, B. L.;
Pedersen, C.;
Pedraza, M.;
Pedurand, R.;
Pele, A.;
Penn, S.;
Perego, A.;
Perez, C. J.;
Perreca, A.;
Pfeiffer, H. P.;
Phelps, M.;
Phukon, K. S.;
Piccinni, O. J.;
Pichot, M.;
Piergiovanni, F.;
Pillant, G.;
Pinard, L.;
Pirello, M.;
Pitkin, M.;
Poggiani, R.;
Pong, D. Y. T.;
Ponrathnam, S.;
Popolizio, P.;
Porter, E. K.;
Powell, J.;
Prajapati, A. K.;
Prasad, J.;
Prasai, K.;
Prasanna, R.;
Pratten, G.;
Prestegard, T.;
Privitera, S.;
Prodi, G. A.;
Prokhorov, L. G.;
Puncken, O.;
Punturo, M.;
Puppo, P.;
Pürrer, M.;
Qi, H.;
Quetschke, V.;
Quinonez, P. J.;
Quintero, E. A.;
Quitzow-James, R.;
Raab, F. J.;
Radkins, H.;
Radulescu, N.;
Raffai, P.;
Raja, S.;
Rajan, C.;
Rajbhandari, B.;
Rakhmanov, M.;
Ramirez, K. E.;
Ramos-Buades, A.;
Rana, Javed;
Rao, K.;
Rapagnani, P.;
Raymond, V.;
Razzano, M.;
Read, J.;
Regimbau, T.;
Rei, L.;
Reid, S.;
Reitze, D. H.;
Ren, W.;
Ricci, F.;
Richardson, C. J.;
Richardson, J. W.;
Ricker, P. M.;
Riemenschneider, G. M.;
Riles, K.;
Rizzo, M.;
Robertson, N. A.;
Robie, R.;
Robinet, F.;
Rocchi, A.;
Rolland, L.;
Rollins, J. G.;
Roma, V. J.;
Romanelli, M.;
Romano, R.;
Romel, C. L.;
Romie, J. H.;
Rose, K.;
Rosińska, D.;
Rosofsky, S. G.;
Ross, M. P.;
Rowan, S.;
Rüdiger, A.;
Ruggi, P.;
Rutins, G.;
Ryan, K.;
Sachdev, S.;
Sadecki, T.;
Sakellariadou, M.;
Salafia, O.;
Salconi, L.;
Saleem, M.;
Salemi, F.;
Samajdar, A.;
Sammut, L.;
Sanchez, E. J.;
Sanchez, L. E.;
Sanchis-Gual, N.;
Sandberg, V.;
Sanders, J. R.;
Santiago, K. A.;
Sarin, N.;
Sassolas, B.;
Sathyaprakash, B. S.;
Saulson, P. R.;
Sauter, O.;
Savage, R. L.;
Schale, P.;
Scheel, M.;
Scheuer, J.;
Schmidt, P.;
Schnabel, R.;
Schofield, R. M. S.;
Schönbeck, A.;
Schreiber, E.;
Schulte, B. W.;
Schutz, B. F.;
Schwalbe, S. G.;
Scott, J.;
Scott, S. M.;
Seidel, E.;
Sellers, D.;
Sengupta, A. S.;
Sennett, N.;
Sentenac, D.;
Sequino, V.;
Sergeev, A.;
Setyawati, Y.;
Shaddock, D. A.;
Shaffer, T.;
Shahriar, M. S.;
Shaner, M. B.;
Shao, L.;
Sharma, P.;
Shawhan, P.;
Shen, H.;
Shink, R.;
Shoemaker, D. H.;
Shoemaker, D. M.;
ShyamSundar, S.;
Siellez, K.;
Sieniawska, M.;
Sigg, D.;
Silva, A. D.;
Singer, L. P.;
Singh, N.;
Singhal, A.;
Sintes, A. M.;
Sitmukhambetov, S.;
Skliris, V.;
Slagmolen, B. J. J.;
Slaven-Blair, T. J.;
Smith, J. R.;
Smith, R. J. E.;
Somala, S.;
Son, E. J.;
Sorazu, B.;
Sorrentino, F.;
Souradeep, T.;
Sowell, E.;
Spencer, A. P.;
Srivastava, A. K.;
Srivastava, V.;
Staats, K.;
Stachie, C.;
Standke, M.;
Steer, D. A.;
Steinke, M.;
Steinlechner, J.;
Steinlechner, S.;
Steinmeyer, D.;
Stevenson, S. P.;
Stocks, D.;
Stone, R.;
Stops, D. J.;
Strain, K. A.;
Stratta, G.;
Strigin, S. E.;
Strunk, A.;
Sturani, R.;
Stuver, A. L.;
Sudhir, V.;
Summerscales, T. Z.;
Sun, L.;
Sunil, S.;
Suresh, J.;
Sutton, P. J.;
Swinkels, B. L.;
Szczepańczyk, M. J.;
Tacca, M.;
Tait, S. C.;
Talbot, C.;
Talukder, D.;
Tanner, D. B.;
Tápai, M.;
Taracchini, A.;
Tasson, J. D.;
Taylor, R.;
Thies, F.;
Thomas, M.;
Thomas, P.;
Thondapu, S. R.;
Thorne, K. A.;
Thrane, E.;
Tiwari, Shubhanshu;
Tiwari, Srishti;
Tiwari, V.;
Toland, K.;
Tonelli, M.;
Tornasi, Z.;
Torres-Forné, A.;
Torrie, C. I.;
Töyrä, D.;
Travasso, F.;
Traylor, G.;
Tringali, M. C.;
Trovato, A.;
Trozzo, L.;
Trudeau, R.;
Tsang, K. W.;
Tse, M.;
Tso, R.;
Tsukada, L.;
Tsuna, D.;
Tuyenbayev, D.;
Ueno, K.;
Ugolini, D.;
Unnikrishnan, C. S.;
Urban, A. L.;
Usman, S. A.;
Vahlbruch, H.;
Vajente, G.;
Valdes, G.;
van Bakel, N.;
van Beuzekom, M.;
Brand, J. F. J. van den;
Broeck, C. Van Den;
Vander-Hyde, D. C.;
van Heijningen, J. V.;
van der Schaaf, L.;
van Veggel, A. A.;
Vardaro, M.;
Varma, V.;
Vass, S.;
Vasúth, M.;
Vecchio, A.;
Vedovato, G.;
Veitch, J.;
Veitch, P. J.;
Venkateswara, K.;
Venugopalan, G.;
Verkindt, D.;
Vetrano, F.;
Viceré, A.;
Viets, A. D.;
Vine, D. J.;
Vinet, J. -Y.;
Vitale, S.;
Vo, T.;
Vocca, H.;
Vorvick, C.;
Vyatchanin, S. P.;
Wade, A. R.;
Wade, L. E.;
Wade, M.;
Walet, R.;
Walker, M.;
Wallace, L.;
Walsh, S.;
Wang, G.;
Wang, H.;
Wang, J. Z.;
Wang, W. H.;
Wang, Y. F.;
Ward, R. L.;
Warden, Z. A.;
Warner, J.;
Was, M.;
Watchi, J.;
Weaver, B.;
Wei, L. -W.;
Weinert, M.;
Weinstein, A. J.;
Weiss, R.;
Wellmann, F.;
Wen, L.;
Wessel, E. K.;
Weßels, P.;
Westhouse, J. W.;
Wette, K.;
Whelan, J. T.;
White, L. V.;
Whiting, B. F.;
Whittle, C.;
Wilken, D. M.;
Williams, D.;
Williamson, A. R.;
Willis, J. L.;
Willke, B.;
Wimmer, M. H.;
Winkler, W.;
Wipf, C. C.;
Wittel, H.;
Woan, G.;
Woehler, J.;
Wofford, J. K.;
Worden, J.;
Wright, J. L.;
Wu, D. S.;
Wysocki, D. M.;
Xiao, L.;
Yamamoto, H.;
Yancey, C. C.;
Yang, L.;
Yap, M. J.;
Yazback, M.;
Yeeles, D. W.;
Yu, Hang;
Yu, Haocun;
Yuen, S. H. R.;
Yvert, M.;
Zadrożny, A. K.;
Zanolin, M.;
Zappa, F.;
Zelenova, T.;
Zendri, J. -P.;
Zevin, M.;
Zhang, J.;
Zhang, L.;
Zhang, T.;
Zhao, C.;
Zhou, M.;
Zhou, Z.;
Zhu, X. J.;
Zimmerman, A. B.;
Zlochower, Y.;
Zucker, M. E.;
Zweizig, J.
Submitted: 2018-11-30, last modified: 2019-10-04
We present the results from three gravitational-wave searches for coalescing
compact binaries with component masses above 1$\mathrm{M}_\odot$ during the
first and second observing runs of the Advanced gravitational-wave detector
network. During the first observing run (O1), from September $12^\mathrm{th}$,
2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary
black hole mergers were detected. The second observing run (O2), which ran from
November $30^\mathrm{th}$, 2016 to August $25^\mathrm{th}$, 2017, saw the first
detection of gravitational waves from a binary neutron star inspiral, in
addition to the observation of gravitational waves from a total of seven binary
black hole mergers, four of which we report here for the first time: GW170729,
GW170809, GW170818 and GW170823. For all significant gravitational-wave events,
we provide estimates of the source properties. The detected binary black holes
have total masses between $18.6_{-0.7}^{+3.2}\mathrm{M}_\odot$, and
$84.4_{-11.1}^{+15.8} \mathrm{M}_\odot$, and range in distance between
$320_{-110}^{+120}$ Mpc and $2840_{-1360}^{+1400}$ Mpc. No neutron star - black
hole mergers were detected. In addition to highly significant
gravitational-wave events, we also provide a list of marginal event candidates
with an estimated false alarm rate less than 1 per 30 days. From these results
over the first two observing runs, which include approximately one
gravitational-wave detection per 15 days of data searched, we infer merger
rates at the 90% confidence intervals of $110\, -\, 3840$
$\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary neutron stars and $9.7\, -\,
101$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary black holes assuming fixed
population distributions, and determine a neutron star - black hole merger rate
90% upper limit of $610$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$.
[23]
oai:arXiv.org:1901.02903 [pdf] - 1916791
Multiphase gas in the circumgalactic medium: relative role of $t_{\rm
cool}/t_{\rm ff}$ and density fluctuations
Submitted: 2019-01-09, last modified: 2019-07-04
We perform a suite of simulations with realistic gravity and thermal balance
in shells to quantify the role of the ratio of cooling time to the free-fall
time ($t_{\rm cool}/t_{\rm ff}$) and the amplitude of density perturbations
($\delta \rho/\rho$) in the production of multiphase gas in the circumgalactic
medium (CGM). Previous idealized simulations, focussing on small amplitude
perturbations in the intracluster medium (ICM), found that cold gas can
condense out of the hot ICM in global thermal balance when the background
$t_{\rm cool}/t_{\rm ff} \lesssim 10$. Recent observations suggest the presence
of cold gas even when the background profiles have somewhat large values of
${t_{\rm cool}/t_{\rm ff}}$. This partly motivates a better understanding of
additional factors such as large density perturbations that can enhance the
propensity for cooling and condensation even when the background ${t_{\rm
cool}/t_{\rm ff}}$ is high. Such large density contrasts can be seeded by
galaxy wakes or dense cosmological filaments. From our simulations, we
introduce a condensation curve in the $(\delta \rho/\rho)$ - min$(t_{\rm
cool}/t_{\rm ff})$ space, that defines the threshold for condensation of
multiphase gas in the CGM. We show that this condensation curve corresponds to
${(t_{\rm cool}/t_{\rm ff})}_{\rm blob} \lesssim 10$ applied to the overdense
blob instead of the background for which $t_{\rm cool}/t_{\rm ff}$ can be
higher. We also study the modification in the condensation curve by varying
entropy stratification. Steeper (positive) entropy gradients shift the
condensation curve to higher amplitudes of perturbations (i.e., make
condensation difficult). A constant entropy core, applicable to the CGM in
smaller halos, shows condensation over a larger range of radii as compared to
the steeper entropy profiles in the ICM.
[25]
oai:arXiv.org:1906.07200 [pdf] - 1902046
Non-uniqueness of cosmic ray two-fluid equations at shocks and possible
remedies
Submitted: 2019-06-17
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 depend on the choice of the numerical methods (spatial reconstruction,
time stepping, and the CFL number) and the adopted discretization.
Nevertheless, all methods converge to a unique result only if the energy
partition between the thermal and non-thermal fluids at the shock is prescribed
a priori. This highlights the closure problem of the two-fluid equations at
shocks. We suggest a robust method where the solutions are insensitive to the
numerical method. Comparison with the currently used methods, critical test
problems, and future directions are discussed.
[30]
oai:arXiv.org:1903.11212 [pdf] - 1859172
Circumgalactic Gas and the Precipitation Limit
Voit, G. M.;
Babul, A.;
Babyk, Iu.;
Bryan, G. L.;
Chen, H. -W.;
Donahue, M.;
Fielding, D.;
Gaspari, M.;
Li, Y.;
McDonald, M.;
O'Shea, B. W.;
Prasad, D.;
Sharma, P.;
Sun, M.;
Tremblay, G.;
Werk, J.;
Werner, N.;
Zahedy, F.
Submitted: 2019-03-26, last modified: 2019-03-29
During the last decade, numerous and varied observations, along with
increasingly sophisticated numerical simulations, have awakened astronomers to
the central role the circumgalactic medium (CGM) plays in regulating galaxy
evolution. It contains the majority of the baryonic matter associated with a
galaxy, along with most of the metals, and must continually replenish the star
forming gas in galaxies that continue to sustain star formation. And while the
CGM is complex, containing gas ranging over orders of magnitude in temperature
and density, a simple emergent property may be governing its structure and
role. Observations increasingly suggest that the ambient CGM pressure cannot
exceed the limit at which cold clouds start to condense out and precipitate
toward the center of the potential well. If feedback fueled by those clouds
then heats the CGM and causes it to expand, the pressure will drop and the
"rain" will diminish. Such a feedback loop tends to suspend the CGM at the
threshold pressure for precipitation. The coming decade will offer many
opportunities to test this potentially fundamental principle of galaxy
evolution.
[33]
oai:arXiv.org:1808.05231 [pdf] - 1842340
A 1-dimensional hydrodynamic model for accretion, cooling and heating of
gas in dark matter halos from $z=6$ to $z=0$
Submitted: 2018-08-15, last modified: 2019-02-25
We study an idealized 1D model for the evolution of hot gas in dark matter
halos for redshifts $z=[0,6]$. We introduce a numerical setup incorporating
cosmological accretion of gas, along with the growth of the halo, based on the
Van den Bosch model for the average growth of halos as a function of cosmic
time. We evolve one-dimensional Lagrangian shells with radiative cooling of the
gas and heating due to feedback from the gas cooling and moving in toward the
center. A simple Bondi accretion model on to a central black hole is used to
include feedback heating. The setup captures some of the key characteristics of
spherically symmetric accretion onto the halos: formation of virial shocks
slightly outside $r_{200}$ and long-term thermal balance in the form of cooling
and heating cycles. The gas density outside our initial halos at $z=6$ is
constrained by requiring that the baryon fraction within the virial radius for
non-radiative evolution be equal to the universal value at almost all times.
The total mass in the cold phase (taken to be $\sim 10^4$ K) within $40$ kpc is
tracked as a function of the halo mass and redshift. We compare the evolution
of the cold gas mass to the observed stellar-mass versus halo mass relations,
following which, we can constrain the feedback energy required for different
halo masses and redshifts. We also compare and match the hot gas density and
temperature profiles for our most massive halo to those of clusters observed
upto redshift $2$. Our model is thus an improvement over the semi-analytic
models in which isothermal condition and $\rho \propto r^{-2}$ are assumed.
[35]
oai:arXiv.org:1810.00018 [pdf] - 1830471
Turbulence in the intracluster medium: simulations, observables &
thermodynamics
Submitted: 2018-09-28, last modified: 2019-02-07
We conduct two kinds of homogeneous isotropic turbulence simulations relevant
for the intracluster medium (ICM): (i) pure turbulence runs without radiative
cooling; (ii) turbulent heating$+$radiative cooling runs with global thermal
balance. For pure turbulence runs in the subsonic regime, the rms density and
surface brightness (SB) fluctuations vary as the square of the rms Mach number
($\mathcal{M}_{\text{rms}}$). However, with thermal balance, the density and SB
fluctuations $(\delta SB/SB)$ are much larger. These scalings have implications
for translating SB fluctuations into a turbulent velocity, particularly for
cool cores. For thermal balance runs with large (cluster core) scale driving,
both the hot and cold phases of the gas are supersonic. For small scale (one
order of magnitude smaller than the cluster core) driving, multiphase gas forms
on a much longer timescale but $\mathcal{M}_{\text{rms}}$ is smaller. Both
small and large scale driving runs have velocities larger than the Hitomi
results from the Perseus cluster. Thus turbulent heating as the dominant
heating source in cool cluster cores is ruled out if multiphase gas is assumed
to condense out from the ICM. Next we perform thermal balance runs in which we
partition the input energy into thermal and turbulent parts and tune their
relative magnitudes. The contribution of turbulent heating has to be $\lesssim
10\%$ in order for turbulence velocities to match Hitomi observations. If the
dominant source of multiphase gas is not cooling from the ICM (but say uplift
from the central galaxy), the importance of turbulent heating cannot be
excluded.
[38]
oai:arXiv.org:1811.12147 [pdf] - 1791780
Plasma processes at the scales of galaxies and clusters of galaxies
Submitted: 2018-11-29
Matter on large scales in the Universe is distributed in form of sheets,
filaments and spherical halos. Most of the known matter in the Universe is in a
diffuse plasma state - in halos around galaxies and clusters of galaxies, and
in the intergalactic medium tracing the large scale cosmological filaments.
Recent multi-wavelength observations have informed us about the state of this
dilute plasma that covers a range of densities and temperatures. Progress in
high resolution numerical simulations allows us to understand this diffuse
plasma, but there are several puzzles that involve a complex interplay of
processes such as cooling, feedback heating, gravity, and turbulence.
[40]
oai:arXiv.org:1806.09629 [pdf] - 1773527
3D global simulations of RIAFs: convergence, effects of azimuthal extent
and dynamo
Submitted: 2018-06-25, last modified: 2018-10-25
We study the long-term evolution of non-radiative geometrically thick
($H/R\approx 0.5$) accretion flows using 3D global ideal MHD simulations and a
pseudo-Newtonian gravity. We find that resolving the scale height with 42 grid
points is adequate to obtain convergence with the product of quality factors
$\langle \langle Q_{\theta} \rangle \rangle \langle \langle Q_{\phi} \rangle
\rangle \geq 300$ and magnetic tilt angle $\theta_B \sim
13^{\circ}-14^{\circ}$. Like previous global isothermal thin disk simulations,
we find stronger mean magnetic fields for the restricted azimuthal domains.
Imposing periodic boundary conditions with the azimuthal extent smaller than
$2\pi$ make the turbulent field at low $m$ appear as a mean field in the runs
with smaller azimuthal extent. But unlike previous works, we do not find a
monotonic trend in turbulence with the azimuthal extent. We conclude that the
minimum azimuthal extent should be $\geq \pi/2$ to capture the flow structure,
but a full $2 \pi$ extent is necessary to study the dynamo. We find an
intermittent dynamo cycle, with $\alpha$-quenching playing an important role in
the nonlinear saturated state. Unlike previous local studies, we find almost
similar values of kinetic and magnetic $\alpha$-s, giving rise to an irregular
distribution of dynamo-$\alpha$. The effects of dynamical quenching are shown
explicitly for the first time in global simulations of accretion flows.
[41]
oai:arXiv.org:1804.05843 [pdf] - 1732639
Astrophysical radio background cannot explain the EDGES 21-cm signal:
constraints from cooling of non-thermal electrons
Submitted: 2018-04-16, last modified: 2018-08-12
Recently the EDGES experiment has claimed the detection of an absorption
feature centered at 78 MHz. When interpreted as a signature of cosmic dawn,
this feature appears at the correct wavelength (corresponding to a redshift
range of $z\approx15-20$) but is larger by at least a factor of two in
amplitude compared to the standard 21-cm models. One way to explain the excess
radio absorption is by the enhancement of the diffuse radio background at $\nu
= 1.42$ GHz ($\lambda=21$ cm) in the rest frame of the absorbing neutral
hydrogen. Astrophysical scenarios, based on the acceleration of relativistic
electrons by accretion on to supermassive black holes (SMBHs) and by supernovae
(SN) from first stars, have been proposed to produce the enhanced radio
background via synchrotron emission. In this Letter we show that either the
synchrotron or the inverse-Compton (IC) cooling time for such electrons is at
least three orders of magnitude shorter than the duration of the EDGES signal
centered at $z \approx 17$, irrespective of the magnetic field strength. The
synchrotron radio emission at 1.42 GHz due to rapidly cooling electrons is
$\sim 10^3$ times smaller than the non-cooling estimate. Thus astrophysical
scenarios for excess radio background proposed to explain the EDGES signal
appear very unlikely.
[42]
oai:arXiv.org:1804.05877 [pdf] - 1728942
Constraining cosmic ray acceleration in young star clusters using
multi-wavelength observations
Submitted: 2018-04-16
We use $1$D and $3$D two-fluid cosmic ray (CR) hydrodynamic simulations to
investigate the role of CRs in the vicinity of a compact young star cluster. We
model a self-gravitating cloud (density profile $\rho \propto r^{-1}$), include
important thermal and non-thermal processes, and explore two different CR
injection scenarios. We show that if internal shocks in the wind-driving region
are the main site for CR acceleration, then the resulting $\gamma$-ray
luminosity ($L_{\rm \gamma}$) can reach $\approx 5\%$ of the mechanical
luminosity ($L_{\rm w}$), independent of the fraction of wind energy ($\sim
1-20\%$) injected into CRs. In contrast, if the forward/reverse shock of a
bubble is the injection site then $L_{\rm \gamma}$ increases linearly with the
CR injection fraction, as expected analytically. We find that the X-ray
luminosity ($L_{\rm x}$) in the forward/reverse shock injection scenario is
$\gtrsim 10^{-3} L_{\rm w}$, which is $\sim 10$ times larger than in the
central wind-driving injection case. We predict the corresponding range of the
synchrotron radio luminosity. We show how multi-wavelength observations can
constrain the CR parameters. Comparing the predicted multi-wavelength
luminosities with those of 30 Doradus we identify the reverse shock as the most
probable CR injection site, and that thermal conduction is important. We do not
find significant dynamical impact of CRs in our models.
[43]
oai:arXiv.org:1712.02367 [pdf] - 1663246
Magnetized SASI: its mechanism and possible connection to some QPOs in
XRBs
Submitted: 2017-12-06, last modified: 2018-04-10
The presence of a surface at the inner boundary, such as in a neutron star or
a white dwarf, allows the existence of a standing shock in steady spherical
accretion.The standing shock can become unstable in 2D or 3D; this is called
the {\em standing accretion shock instability} (SASI).Two mechanisms --
advective-acoustic and purely acoustic -- have been proposed to explain SASI.
Using axisymmetric hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations,
we find that the advective-acoustic mechanism better matches the observed
oscillation timescales in our simulations. The global shock oscillations
present in the accretion flow can explain many observed high frequency
($\gtrsim 100$ Hz) quasi-periodic oscillations (QPOs) in X-ray binaries (XRBs).
The presence of a moderately strong magnetic field adds more features to the
shock oscillation pattern, giving rise to low frequency modulation in the
computed light curve. This low frequency modulation can be responsible for
$\sim 100$ Hz QPOs (known as hHz QPOs). We propose that the appearance of hHz
QPO determines the separation of twin peak QPOs of higher frequencies.
[44]
oai:arXiv.org:1801.04282 [pdf] - 1732566
Cool-Core Clusters : Role of BCG, Star Formation & AGN-Driven Turbulence
Submitted: 2018-01-12, last modified: 2018-03-16
Recent analysis shows that it is important to explicitly include the
gravitational potential of the central brightest central galaxy (BCG) to infer
the acceleration due to gravity ($g$) and the free-fall time ($t_{\rm ff}
\equiv [2r/g]^{1/2}$) in cool cluster cores. Accurately measuring $t_{\rm ff}$
is crucial because according to numerical simulations cold gas condensation and
strong feedback occur in cluster cores with min($t_{\rm cool}/t_{\rm ff}$)
below a threshold value close to 10. Recent observations which include the BCG
gravity show that the observed threshold in min($t_{\rm cool}/t_{\rm ff}$) lies
at a somewhat higher value, close to 10-30; there are only a few clusters in
which this ratio falls much below 10. In this paper we compare numerical
simulations of feedback AGN (Active Galactic Nuclei) jets interacting with the
intracluster medium (ICM), with and without a BCG potential. We find that, for
a fixed feedback efficiency, the presence of a BCG does not significantly
affect the temperature but increases (decreases) the core density (entropy) on
average. Most importantly, min($t_{\rm cool}/t_{\rm ff}$) is only affected
slightly by the inclusion of the BCG gravity. Also notable is that the lowest
value of min($t_{\rm cool}/t_{\rm ff}$) in the NFW+BCG runs are about twice
larger than in the NFW runs. We also look at the role of depletion of cold gas
due to star formation and show that it only affects the rotationally dominant
component (torus), while the radially dominant component (which regulates the
feedback cycle) remains largely unaffected. The distribution of metals due to
AGN jets in our simulations is predominantly along the jet direction and the
radial spread of metals is less. We also show that the turbulence in cool core
clusters is weak, consistent with recent Hitomi results on Perseus cluster.
[45]
oai:arXiv.org:1712.06252 [pdf] - 1634260
Extra-planar X-ray emission from disc-wide outflows in spiral galaxies
Submitted: 2017-12-18, last modified: 2018-01-30
We study the effects of mass and energy injection due to OB associations
spread across the rotating disc of a Milky Way-type galaxy, with the help of 3D
hydrodynamic simulations. We compare the resulting X-ray emission with that
produced from the injection of mass and energy from a central region. We find
that the predicted X-ray image shows a filamentary structure that arises even
in the absence of disc gas inhomogeneity. This structure stems from warm clumps
made of disc material being lifted by the injected gas. We show that as much as
half of the total X-ray emission comes from regions surrounding warm clumps
that are made of a mix of disk and injected gas. This scenario has the
potential to explain the origin of the observed extra-planar X-ray emission
around star forming galaxies and can be used to understand the observed
sublinear relation between the $L_X$ and SFR. We quantify the mass contained in
these `bow-shock' regions. We also show that the top-most region of the outer
shock above the central area emits harder X-rays than the rest. Further, we
find that the mass distribution in different temperature ranges is bimodal,
peaking at $10^4\hbox{-}10^5$ K (in warm clumps) and $10^6\hbox{-}10^7$ K
(X-ray emitting gas). The mass loading factor is found to decrease with
increasing SFR, consistent with previous theoretical estimates and simulations.
[46]
oai:arXiv.org:1705.04444 [pdf] - 1604213
AGN feedback with the Square Kilometer Array (SKA) and implications for
cluster physics and cosmology
Submitted: 2017-05-12, last modified: 2017-10-17
AGN feedback is regarded as an important non-gravitational process in galaxy
clusters, providing useful constraints on large-scale structure formation. It
modifies the structure and energetics of the intra-cluster medium (ICM) and
hence its understanding is crucially needed in order to use clusters as high
precision cosmological probes. In this context, particularly keeping in mind
the upcoming high quality radio data expected from radio surveys like SKA with
its higher sensitivity, high spatial and spectral resolutions, we review our
current understanding of AGN feedback, its cosmological implications and the
impact that SKA can have in revolutionizing our understanding of AGN feedback
in large-scale structures. Recent developments regarding the AGN outbursts and
its possible contribution to excess entropy in the hot atmospheres of groups
and clusters, its correlation with the feedback energy in ICM, quenching of
cooling flows and the possible connection between cool core clusters and radio
mini-halos, are discussed. We describe current major issues regarding modeling
of AGN feedback and its impact on the surrounding medium. With regard to the
future of AGN feedback studies, we examine the possible breakthroughs that can
be expected from SKA observations. In the context of cluster cosmology, for
example, we point out the importance of SKA observations for cluster mass
calibration by noting that most of $z>1$ clusters discovered by eROSITA X-ray
mission can be expected to be followed up through a 1000 hour SKA-1 mid
programme. Moreover, approximately $1000$ radio mini halos and $\sim 2500$
radio halos at $z<0.6$ can be potentially detected by SKA1 and SKA2 and used as
tracers of galaxy clusters and determination of cluster selection function.
[47]
oai:arXiv.org:1702.05487 [pdf] - 1581533
Scalable explicit implementation of anisotropic diffusion with
Runge-Kutta-Legendre super-time-stepping
Submitted: 2017-02-17, last modified: 2017-08-22
An important ingredient in numerical modelling of high temperature magnetised
astrophysical plasmas is the anisotropic transport of heat along magnetic field
lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically
involves solving the hyperbolic set of conservation equations along with the
induction equation. Incorporating anisotropic thermal conduction requires to
also treat parabolic terms arising from the diffusion operator. An explicit
treatment of parabolic terms will considerably reduce the simulation time step
due to its dependence on the square of the grid resolution ($\Delta x$) for
stability. Although an implicit scheme relaxes the constraint on stability, it
is difficult to distribute efficiently on a parallel architecture. Treating
parabolic terms with accelerated super-time stepping (STS) methods has been
discussed in literature but these methods suffer from poor accuracy (first
order in time) and also have difficult-to-choose tuneable stability parameters.
In this work we highlight a second order (in time) Runge Kutta Legendre (RKL)
scheme (first described by Meyer et. al. 2012) that is robust, fast and
accurate in treating parabolic terms alongside the hyperbolic conversation
laws. We demonstrate its superiority over the first order super time stepping
schemes with standard tests and astrophysical applications. We also show that
explicit conduction is particularly robust in handling saturated thermal
conduction. Parallel scaling of explicit conduction using RKL scheme is
demonstrated up to more than $10^4$ processors.
[48]
oai:arXiv.org:1611.02710 [pdf] - 1580561
AGN jet-driven stochastic cold accretion in cluster cores
Submitted: 2016-11-08, last modified: 2017-07-07
Several arguments suggest that stochastic condensation of cold gas and its
accretion onto the central supermassive black hole (SMBH) is essential for
active galactic nuclei (AGN) feedback to work in the most massive galaxies that
lie at the centres of galaxy clusters. Our 3-D hydrodynamic AGN jet-ICM
(intracluster medium) simulations, looking at the detailed angular momentum
distribution of cold gas and its time variability for the first time, show that
the angular momentum of the cold gas crossing $\lesssim 1$ kpc is essentially
isotropic. With almost equal mass in clockwise and counter-clockwise
orientations, we expect a cancellation of angular momentum on roughly the
dynamical time. This means that a compact accretion flow with a short viscous
time ought to form, through which enough accretion power can be channeled into
jet mechanical energy sufficiently quickly to prevent a cooling flow. The
inherent stochasticity, expected in feedback cycles driven by cold gas
condensation, gives rise to a large variation in the cold gas mass at the
centres of galaxy clusters, for similar cluster and SMBH masses, in agreement
with the observations. Such correlations are expected to be much tighter for
the smoother hot/Bondi accretion. The weak correlation between cavity power and
Bondi power obtained from our simulations also match observations.
[49]
oai:arXiv.org:1603.00815 [pdf] - 1579843
How multiple supernovae overlap to form superbubbles
Submitted: 2016-03-02, last modified: 2017-06-20
We explore the formation of superbubbles through energy deposition by
multiple supernovae (SNe) in a uniform medium. We use total energy conserving,
3-D hydrodynamic simulations to study how SNe correlated in space and time
create superbubbles. While isolated SNe fizzle out completely by $\sim 1$ Myr
due to radiative losses, for a realistic cluster size it is likely that
subsequent SNe go off within the hot/dilute bubble and sustain the shock till
the cluster lifetime. For realistic cluster sizes, we find that the bubble
remains overpressured only if, for a given $n_{g0}$, $N_{\rm OB}$ is
sufficiently large. While most of the input energy is still lost radiatively,
superbubbles can retain up to $\sim 5-10\%$ of the input energy in form of
kinetic+thermal energy till 10 Myr for ISM density $n_{g0} \approx 1$
cm$^{-3}$. We find that the mechanical efficiency decreases for higher
densities ($\eta_{\rm mech} \propto n_{g0}^{-2/3}$). We compare the radii and
velocities of simulated supershells with observations and the classical
adiabatic model. Our simulations show that the superbubbles retain only
$\lesssim 10\%$ of the injected energy, thereby explaining the observed smaller
size and slower expansion of supershells. We also confirm that a sufficiently
large ($\gtrsim 10^4$) number of SNe is required to go off in order to create a
steady wind with a stable termination shock within the superbubble. We show
that the mechanical efficiency increases with increasing resolution, and that
explicit diffusion is required to obtain converged results.
[50]
oai:arXiv.org:1705.10448 [pdf] - 1584017
Lack of thermal energy in superbubbles: hint of cosmic rays?
Submitted: 2017-05-29
Using analytic methods and $1$-D two-fluid simulations, we study the effect
of cosmic rays (CRs) on the dynamics of interstellar superbubbles (ISBs) driven
by multiple supernovae (SNe)/stellar winds in OB associations. In addition to
CR advection and diffusion, our models include thermal conduction and radiative
cooling. We find that CR injection at the reverse shock or within a central
wind-driving region can affect the thermal profiles of ISBs and hence their
X-ray properties. Even if a small fraction ($10-20\%$) of the total mechanical
power is injected into CRs, a significant fraction of the ram pressure at the
reverse shock can be transferred to CRs. The energy transfer becomes efficient
if (1) the reverse shock gas Mach number exceeds a critical value ($M_{\rm
th}\gtrsim 12$) and (2) the CR acceleration time scale $\tau_{\rm acc}\sim
\kappa_{\rm cr}/v^2$ is shorter than the dynamical time, where $\kappa_{\rm
cr}$ is CR diffusion constant and $v$ is the upstream velocity. We show that CR
affected bubbles can exhibit a volume averaged hot gas temperature
$1-5\times10^{6}$ K, lower by a factor of $2-10$ than without CRs. Thus CRs can
potentially solve the long-standing problem of the observed low ISB
temperatures.
[51]
oai:arXiv.org:1601.02347 [pdf] - 1579796
Decoding X-ray observations from centres of galaxy clusters using MCMC
Submitted: 2016-01-11, last modified: 2017-05-26
We correct for the use of electron densities instead of total gas density in
the pressure fits which were used to derive the local free-fall times (t$_{\rm
ff}$) in the original paper.
[52]
oai:arXiv.org:1610.00719 [pdf] - 1550336
Clues to the origin of Fermi Bubbles from OVIII/OVII line ratio
Submitted: 2016-10-03, last modified: 2017-03-07
We constrain the origin of Fermi Bubbles using 2D hydrodynamical simulations
of both star formation driven and black hole accretion driven wind models. We
compare our results with recent observations of OVIII to OVII line ratio within
and near Fermi Bubbles. Our results suggest that independent of the driving
mechanisms, a low luminosity ($\mathcal{L} \sim 0.7-1\times 10^{41}$ erg
s$^{-1}$) energy injection best reproduces the observed line ratio for which
the shock temperature is $\approx 3\times 10^6$ K. Assuming the Galactic halo
temperature to be $2\times 10^6$K, we estimate the shock velocity to be $\sim
300$ km s$^{-1}$ for a weak shock. The corresponding estimated age of the Fermi
bubbles is $\sim 15-25$ Myr. Such an event can be produced either by a star
formation rate of $\sim 0.5$ M$_\odot$ yr$^{-1}$ at the Galactic centre or a
very low luminosity jet/accretion wind arising from the central black hole. Our
analysis rules out any activity that generates an average mechanical luminosity
$\gtrsim 10^{41}$ \ergps as a possible origin of the Fermi Bubbles.
[53]
oai:arXiv.org:1611.07330 [pdf] - 1527219
Detection of polarized quasi-periodic microstructure emission in
millisecond pulsars
Submitted: 2016-11-22
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 $\mu$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.
[54]
oai:arXiv.org:1610.08182 [pdf] - 1521422
Clusters of galaxies and the cosmic web with SKA
Kale, Ruta;
Dwarakanath, K. S.;
Lal, Dharam Vir;
Bagchi, Joydeep;
Paul, Surajit;
Malu, Siddharth;
Datta, Abhirup;
Parekh, Viral;
Sharma, Prateek;
Pandey-Pommier, Mamta
Submitted: 2016-10-26
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]
[55]
oai:arXiv.org:1610.03159 [pdf] - 1494915
The Astropy Problem
Muna, Demitri;
Alexander, Michael;
Allen, Alice;
Ashley, Richard;
Asmus, Daniel;
Azzollini, Ruyman;
Bannister, Michele;
Beaton, Rachael;
Benson, Andrew;
Berriman, G. Bruce;
Bilicki, Maciej;
Boyce, Peter;
Bridge, Joanna;
Cami, Jan;
Cangi, Eryn;
Chen, Xian;
Christiny, Nicholas;
Clark, Christopher;
Collins, Michelle;
Comparat, Johan;
Cook, Neil;
Croton, Darren;
Davids, Isak Delberth;
Depagne, Éric;
Donor, John;
Santos, Leonardo A. dos;
Douglas, Stephanie;
Du, Alan;
Durbin, Meredith;
Erb, Dawn;
Faes, Daniel;
Fernández-Trincado, J. G.;
Foley, Anthony;
Fotopoulou, Sotiria;
Frimann, Søren;
Frinchaboy, Peter;
Garcia-Dias, Rafael;
Gawryszczak, Artur;
George, Elizabeth;
Gonzalez, Sebastian;
Gordon, Karl;
Gorgone, Nicholas;
Gosmeyer, Catherine;
Grasha, Katie;
Greenfield, Perry;
Grellmann, Rebekka;
Guillochon, James;
Gurwell, Mark;
Haas, Marcel;
Hagen, Alex;
Haggard, Daryl;
Haines, Tim;
Hall, Patrick;
Hellwing, Wojciech;
Herenz, Edmund Christian;
Hinton, Samuel;
Hlozek, Renee;
Hoffman, John;
Holman, Derek;
Holwerda, Benne Willem;
Horton, Anthony;
Hummels, Cameron;
Jacobs, Daniel;
Jensen, Jens Juel;
Jones, David;
Karick, Arna;
Kelley, Luke;
Kenworthy, Matthew;
Kitchener, Ben;
Klaes, Dominik;
Kohn, Saul;
Konorski, Piotr;
Krawczyk, Coleman;
Kuehn, Kyler;
Kuutma, Teet;
Lam, Michael T.;
Lane, Richard;
Liske, Jochen;
Lopez-Camara, Diego;
Mack, Katherine;
Mangham, Sam;
Mao, Qingqing;
Marsh, David J. E.;
Mateu, Cecilia;
Maurin, Loïc;
McCormac, James;
Momcheva, Ivelina;
Monteiro, Hektor;
Mueller, Michael;
Munoz, Roberto;
Naidu, Rohan;
Nelson, Nicholas;
Nitschelm, Christian;
North, Chris;
Nunez-Iglesias, Juan;
Ogaz, Sara;
Owen, Russell;
Parejko, John;
Patrício, Vera;
Pepper, Joshua;
Perrin, Marshall;
Pickering, Timothy;
Piscionere, Jennifer;
Pogge, Richard;
Poleski, Radek;
Pourtsidou, Alkistis;
Price-Whelan, Adrian M.;
Rawls, Meredith L.;
Read, Shaun;
Rees, Glen;
Rein, Hanno;
Rice, Thomas;
Riemer-Sørensen, Signe;
Rusomarov, Naum;
Sanchez, Sebastian F.;
Santander-García, Miguel;
Sarid, Gal;
Schoenell, William;
Scholz, Aleks;
Schuhmann, Robert L.;
Schuster, William;
Scicluna, Peter;
Seidel, Marja;
Shao, Lijing;
Sharma, Pranav;
Shulevski, Aleksandar;
Shupe, David;
Sifón, Cristóbal;
Simmons, Brooke;
Sinha, Manodeep;
Skillen, Ian;
Soergel, Bjoern;
Spriggs, Thomas;
Srinivasan, Sundar;
Stevens, Abigail;
Streicher, Ole;
Suchyta, Eric;
Tan, Joshua;
Telford, O. Grace;
Thomas, Romain;
Tonini, Chiara;
Tremblay, Grant;
Tuttle, Sarah;
Urrutia, Tanya;
Vaughan, Sam;
Verdugo, Miguel;
Wagner, Alexander;
Walawender, Josh;
Wetzel, Andrew;
Willett, Kyle;
Williams, Peter K. G.;
Yang, Guang;
Zhu, Guangtun;
Zonca, Andrea
Submitted: 2016-10-10
The Astropy Project (http://astropy.org) is, in its own words, "a community
effort to develop a single core package for Astronomy in Python and foster
interoperability between Python astronomy packages." For five years this
project has been managed, written, and operated as a grassroots,
self-organized, almost entirely volunteer effort while the software is used by
the majority of the astronomical community. Despite this, the project has
always been and remains to this day effectively unfunded. Further, contributors
receive little or no formal recognition for creating and supporting what is now
critical software. This paper explores the problem in detail, outlines possible
solutions to correct this, and presents a few suggestions on how to address the
sustainability of general purpose astronomical software.
[56]
oai:arXiv.org:1607.04328 [pdf] - 1470764
Molecular outflows in starburst nuclei
Submitted: 2016-07-14, last modified: 2016-08-20
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 \Sigma_{SFR} \le 50$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$
favours the production of molecular outflows, consistent with observed values.
[57]
oai:arXiv.org:1604.08214 [pdf] - 1451239
Spherical accretion: the influence of inner boundary and quasi-periodic
oscillations
Submitted: 2016-04-27, last modified: 2016-08-06
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.
[58]
oai:arXiv.org:1606.09127 [pdf] - 1457448
How radiation affects superbubbles: Through momentum injection in early
phase and photo-heating thereafter
Submitted: 2016-06-29, last modified: 2016-08-01
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
($\sigma_{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.
[59]
oai:arXiv.org:1601.00394 [pdf] - 1359139
Diffuse X-ray emission from star forming galaxies
Submitted: 2016-01-04, last modified: 2016-01-25
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.
[60]
oai:arXiv.org:1512.01217 [pdf] - 1356850
Cold gas in cluster cores: Global stability analysis and non-linear
simulations of thermal instability
Submitted: 2015-12-03, last modified: 2016-01-20
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.
[61]
oai:arXiv.org:1310.2242 [pdf] - 1283706
Thermal Conduction and Multiphase Gas in Cluster Cores
Submitted: 2013-10-08, last modified: 2015-09-26
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.
[62]
oai:arXiv.org:1504.02215 [pdf] - 1288347
Cool core cycles: Cold gas and AGN jet feedback in cluster cores
Submitted: 2015-04-09, last modified: 2015-09-24
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 \epsilon \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)
[63]
oai:arXiv.org:1505.03634 [pdf] - 1284810
Multi-wavelength features of Fermi Bubbles as signatures of a Galactic
wind
Submitted: 2015-05-14, last modified: 2015-08-05
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.
[64]
oai:arXiv.org:1412.7924 [pdf] - 1019207
Narrow escape: how ionizing photons escape from disc galaxies
Submitted: 2014-12-26, last modified: 2015-05-13
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}^\alpha n_0^2 z_0^3$ (with $\alpha\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.
[65]
oai:arXiv.org:1409.4874 [pdf] - 934811
Long way to go: how outflows from large galaxies propagate through the
hot halo gas
Submitted: 2014-09-17, last modified: 2015-01-27
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, $\eta_{3e5}$
($T\le 3 \times 10^5$ K) is related to the mass loading factor at the virial
radius ($\eta_{v}$) as $\eta_{v} \approx 25\, \eta_{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.
[66]
oai:arXiv.org:1409.1220 [pdf] - 1216677
The cold mode: A phenomenological model for the evolution of density
perturbations in the intracluster medium
Submitted: 2014-09-03, last modified: 2014-12-31
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$.
[67]
oai:arXiv.org:1403.3395 [pdf] - 1208370
Turbulence and cooling in galaxy cluster cores
Submitted: 2014-03-13, last modified: 2014-11-08
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.
[68]
oai:arXiv.org:1402.6695 [pdf] - 1203599
In a hot bubble: why does superbubble feedback work, but isolated
supernovae do not?
Submitted: 2014-02-26, last modified: 2014-10-25
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.
[69]
oai:arXiv.org:1304.1294 [pdf] - 1165732
Radiatively Inefficient Accretion Flow Simulations with Cooling:
Implications for Black Hole Transients
Submitted: 2013-04-04, last modified: 2013-08-03
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.
[70]
oai:arXiv.org:1303.2664 [pdf] - 1165173
Superbubble breakout and galactic winds from disk galaxies
Submitted: 2013-03-11, last modified: 2013-07-13
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.
[71]
oai:arXiv.org:1304.2408 [pdf] - 649422
Astrophysical Coronae: Lessons from Modeling of the Intracluster Medium
Submitted: 2013-04-08
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.
[72]
oai:arXiv.org:1209.5748 [pdf] - 1151644
Isotropic Heating of Galaxy Cluster Cores via Rapidly Reorienting AGN
Jets
Submitted: 2012-09-25, last modified: 2013-03-14
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.
[73]
oai:arXiv.org:1105.2563 [pdf] - 598094
Thermal Instability in Gravitationally-Stratified Plasmas: Implications
for Multi-Phase Structure in Clusters and Galaxy Halos
Submitted: 2011-05-12, last modified: 2012-12-04
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.
[74]
oai:arXiv.org:1206.4314 [pdf] - 1124241
On the Structure of Hot Gas in Halos: Implications for the Lx-Tx
Relation & Missing Baryons
Submitted: 2012-06-19, last modified: 2012-09-16
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.
[75]
oai:arXiv.org:1201.6407 [pdf] - 1116267
Local 2D Particle-in-cell simulations of the collisionless MRI
Submitted: 2012-01-30
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\'en 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*.
[76]
oai:arXiv.org:1201.0754 [pdf] - 483873
The Effects of Anisotropic Viscosity on Turbulence and Heat Transport in
the Intracluster Medium
Submitted: 2012-01-03
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.
[77]
oai:arXiv.org:1110.6063 [pdf] - 1085195
Cause and Effect of Feedback: Multiphase Gas in Cluster Cores Heated by
AGN Jets
Submitted: 2011-10-27
Multiwavelength data indicate that the X-ray emitting plasma in the cores of
galaxy clusters is not cooling catastrophically. To large extent, cooling is
offset by heating due to active galactic nuclei (AGN) via jets. The cool-core
clusters, with cooler/denser plasmas, show multiphase gas and signs of some
cooling in their cores. These observations suggest that the cool core is
locally thermally unstable while maintaining global thermal equilibrium. Using
high-resolution, three-dimensional simulations we study the formation of
multiphase gas in cluster cores heated by highly-collimated bipolar AGN jets.
Our key conclusion is that spatially extended multiphase filaments form only
when the instantaneous ratio of the thermal instability and free-fall
timescales (t_TI/t_ff) falls below a critical threshold of \approx 10. When
this happens, dense cold gas decouples from the hot ICM phase and generates
inhomogeneous and spatially extended Halpha filaments. These cold gas clumps
and filaments `rain' down onto the central regions of the core, forming a cold
rotating torus and in part feeding the supermassive black hole. Consequently,
the self-regulated feedback enhances AGN heating and the core returns to a
higher entropy level with t_TI/t_ff > 10. Eventually the core reaches
quasi-stable global thermal equilibrium, and cold filaments condense out of the
hot ICM whenever t_TI/t_ff \lesssim 10. This occurs despite the fact that the
energy from AGN jets is supplied to the core in a highly anisotropic fashion.
The effective spatial redistribution of heat is enabled in part by the
turbulent motions in the wake of freely-falling cold filaments. Increased AGN
activity can locally reverse the cold gas flow, launching cold filamentary gas
away from the cluster center. Our criterion for the condensation of spatially
extended cold gas is in agreement with observations and previous idealized
simulations.
[78]
oai:arXiv.org:1109.1285 [pdf] - 1083881
Turbulent Pressure Support in the Outer Parts of Galaxy Clusters
Submitted: 2011-09-06, last modified: 2011-10-07
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.
[79]
oai:arXiv.org:1106.4816 [pdf] - 1077505
Thermal Instability & the Feedback Regulation of Hot Halos in Clusters,
Groups, and Galaxies
Submitted: 2011-06-23
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.
[80]
oai:arXiv.org:0909.5426 [pdf] - 951073
Numerical Implementation of Streaming Down the Gradient: Application to
Fluid Modeling of Cosmic Rays and Saturated Conduction
Submitted: 2009-09-29, last modified: 2011-05-03
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 $\Delta t \propto \Delta x^2$ behaves fine. We show
that the implicit methods show stable and converging results with $\Delta t
\propto \Delta 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.
[81]
oai:arXiv.org:1008.0861 [pdf] - 1034107
A Fast Semi-implicit Method for Anisotropic Diffusion
Submitted: 2010-08-04, last modified: 2011-03-07
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.
[82]
oai:arXiv.org:0901.2676 [pdf] - 296765
Spacetime Singularities in String and its Low Dimensional Effective
Theory
Submitted: 2009-01-17, last modified: 2010-12-12
Spacetime singularities are studied in both the $D+d$-dimensional string
theory and its $D$-dimensional effective theory, obtained by the Kaluza-Klein
compactification. It is found that spacetime singularities in the low
dimensional effective theory may or may not remain after lifted to the
$D+d$-dimensional string theory, depending on particular solutions. It is also
found that there exist cases in which spacetime singularities appearing in
high/low dimensional spacetimes do not necessarily happen on the same surfaces.
[83]
oai:arXiv.org:1005.0389 [pdf] - 270052
Time Dependent Models of Flares from Sagittarius A*
Submitted: 2010-05-03, last modified: 2010-10-07
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.
[84]
oai:arXiv.org:1009.2498 [pdf] - 598093
Can conduction induce convection? The non-linear saturation of buoyancy
instabilities in dilute plasmas
Submitted: 2010-09-13
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.
[85]
oai:arXiv.org:0912.1851 [pdf] - 902689
Turbulence in Galaxy Cluster Cores: a Key to Cluster Bimodality?
Submitted: 2009-12-09, last modified: 2010-03-26
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.
[86]
oai:arXiv.org:1003.5546 [pdf] - 1025985
Thermal Instability with Anisotropic Thermal Conduction and Adiabatic
Cosmic Rays: Implications for Cold Filaments in Galaxy Clusters
Submitted: 2010-03-25
Observations of the cores of nearby galaxy clusters show H$\alpha$ 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.
[87]
oai:arXiv.org:0909.0270 [pdf] - 1017333
Turbulence and Mixing in the Intracluster Medium
Submitted: 2009-09-01
The intracluster medium (ICM) is stably stratified in the hydrodynamic sense
with the entropy $s$ increasing outwards. However, thermal conduction along
magnetic field lines fundamentally changes the stability of the ICM, leading to
the "heat-flux buoyancy instability" when $dT/dr>0$ and the "magnetothermal
instability" when $dT/dr<0$. The ICM is thus buoyantly unstable regardless of
the signs of $dT/dr$ and $ds/dr$. On the other hand, these
temperature-gradient-driven instabilities saturate by reorienting the magnetic
field (perpendicular to $\hat{\bf r}$ when $dT/dr>0$ and parallel to $\hat{\bf
r}$ when $dT/dr<0$), without generating sustained convection. We show that
after an anisotropically conducting plasma reaches this nonlinearly stable
magnetic configuration, it experiences a buoyant restoring force that resists
further distortions of the magnetic field. This restoring force is analogous to
the buoyant restoring force experienced by a stably stratified adiabatic
plasma. We argue that in order for a driving mechanism (e.g, galaxy motions or
cosmic-ray buoyancy) to overcome this restoring force and generate turbulence
in the ICM, the strength of the driving must exceed a threshold, corresponding
to turbulent velocities $\gtrsim 10 -100 {km/s}$. For weaker driving, the ICM
remains in its nonlinearly stable magnetic configuration, and turbulent mixing
is effectively absent. We discuss the implications of these findings for the
turbulent diffusion of metals and heat in the ICM.
[88]
oai:arXiv.org:0905.4500 [pdf] - 24694
Anisotropic Thermal Conduction and the Cooling Flow Problem in Galaxy
Clusters
Submitted: 2009-05-27
We examine the long-standing cooling flow problem in galaxy clusters with 3D
MHD simulations of isolated clusters including radiative cooling and
anisotropic thermal conduction along magnetic field lines. The central regions
of the intracluster medium (ICM) can have cooling timescales of ~200 Myr or
shorter--in order to prevent a cooling catastrophe the ICM must be heated by
some mechanism such as AGN feedback or thermal conduction from the thermal
reservoir at large radii. The cores of galaxy clusters are linearly unstable to
the heat-flux-driven buoyancy instability (HBI), which significantly changes
the thermodynamics of the cluster core. The HBI is a convective,
buoyancy-driven instability that rearranges the magnetic field to be
preferentially perpendicular to the temperature gradient. For a wide range of
parameters, our simulations demonstrate that in the presence of the HBI, the
effective radial thermal conductivity is reduced to less than 10% of the full
Spitzer conductivity. With this suppression of conductive heating, the cooling
catastrophe occurs on a timescale comparable to the central cooling time of the
cluster. Thermal conduction alone is thus unlikely to stabilize clusters with
low central entropies and short central cooling timescales. High central
entropy clusters have sufficiently long cooling times that conduction can help
stave off the cooling catastrophe for cosmologically interesting timescales.
[89]
oai:arXiv.org:0901.4786 [pdf] - 900592
Buoyancy Instabilities in Galaxy Clusters: Convection Due to Adiabatic
Cosmic Rays and Anisotropic Thermal Conduction
Submitted: 2009-01-29
Using a linear stability analysis and two and three-dimensional nonlinear
simulations, we study the physics of buoyancy instabilities in a combined
thermal and relativistic (cosmic ray) plasma, motivated by the application to
clusters of galaxies. We argue that cosmic ray diffusion is likely to be slow
compared to the buoyancy time on large length scales, so that cosmic rays are
effectively adiabatic. If the cosmic ray pressure $p_{cr}$ is $\gtrsim 25 %$ of
the thermal pressure, and the cosmic ray entropy ($p_{\rm cr}/\rho^{4/3}$;
$\rho$ is the thermal plasma density) decreases outwards, cosmic rays drive an
adiabatic convective instability analogous to Schwarzschild convection in
stars. Global simulations of galaxy cluster cores show that this instability
saturates by reducing the cosmic ray entropy gradient and driving efficient
convection and turbulent mixing. At larger radii in cluster cores, the thermal
plasma is unstable to the heat flux-driven buoyancy instability (HBI), a
convective instability generated by anisotropic thermal conduction and a
background conductive heat flux. Cosmic-ray driven convection and the HBI may
contribute to redistributing metals produced by Type 1a supernovae in clusters.
Our calculations demonstrate that adiabatic simulations of galaxy clusters can
artificially suppress the mixing of thermal and relativistic plasma;
anisotropic thermal conduction allows more efficient mixing, which may
contribute to cosmic rays being distributed throughout the cluster volume.
[90]
oai:arXiv.org:0804.1353 [pdf] - 314909
Spherical Accretion with Anisotropic Thermal Conduction
Submitted: 2008-04-08, last modified: 2008-07-05
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.
[91]
oai:arXiv.org:0706.3715 [pdf] - 314879
Faraday Rotation in Global Accretion Disk Simulations: Implications for
Sgr A*
Submitted: 2007-06-25, last modified: 2007-09-18
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.
[92]
oai:arXiv.org:astro-ph/0703572 [pdf] - 316975
Electron Heating in Hot Accretion Flows
Submitted: 2007-03-21, last modified: 2007-09-18
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.
[93]
oai:arXiv.org:0707.2616 [pdf] - 3182
Preserving Monotonicity in Anisotropic Diffusion
Submitted: 2007-07-17
We show that standard algorithms for anisotropic diffusion based on centered
differencing (including the recent symmetric algorithm) do not preserve
monotonicity. In the context of anisotropic thermal conduction, this can lead
to the violation of the entropy constraints of the second law of
thermodynamics, causing heat to flow from regions of lower temperature to
higher temperature. In regions of large temperature variations, this can cause
the temperature to become negative. Test cases to illustrate this for centered
asymmetric and symmetric differencing are presented. Algorithms based on slope
limiters, analogous to those used in second order schemes for hyperbolic
equations, are proposed to fix these problems. While centered algorithms may be
good for many cases, the main advantage of limited methods is that they are
guaranteed to avoid negative temperature (which can cause numerical
instabilities) in the presence of large temperature gradients. In particular,
limited methods will be useful to simulate hot, dilute astrophysical plasmas
where conduction is anisotropic and the temperature gradients are enormous,
e.g., collisionless shocks and disk-corona interface.
[94]
oai:arXiv.org:astro-ph/0703542 [pdf] - 90351
Kinetic Effects on Turbulence Driven by the Magnetorotational
Instability in Black Hole Accretion
Submitted: 2007-03-20
Magnetorotational Instability (MRI), the instability causing turbulent
transport in accretion disks, is studied in the kinetic regime. Radiatively
Inefficient Accretion Flows (RIAFs), like the one around the supermassive black
hole in the center of our Galaxy, are believed to be collisionless. Kinetic MHD
formalism, based on the moments of the Vlasov equation, is used for linear
analysis and nonlinear simulations.
ZEUS MHD code is modified to include key kinetic MHD terms: anisotropic
pressure tensor and anisotropic thermal conduction. Simulations use the local
shearing box approximation. Pressure anisotropy is created, because of the
adiabatic invariance ($\mu=p_\perp/B$), as magnetic field is amplified by the
MRI. Larmor radius scale instabilities--mirror, ion-cyclotron, and
firehose--are excited at large pressure anisotropy. Pressure isotropization due
to pitch angle scattering by these instabilities is included as a subgrid
model. A key result of the kinetic MHD simulations is that the anisotropic
(viscous) stress can be as large as the Maxwell stress.
A new numerical method to simulate anisotropic thermal conduction with large
temperature gradients is suggested. Simple tests show that the centered
differencing of anisotropic thermal conduction can result in heat flowing from
lower to higher temperatures, giving rise to negative temperatures. Limiting of
transverse temperature gradients does not accentuate temperature extrema.
[95]
oai:arXiv.org:astro-ph/0508502 [pdf] - 75391
Shearing Box Simulations of the MRI in a Collisionless Plasma
Submitted: 2005-08-23, last modified: 2005-08-25
We describe local shearing box simulations of turbulence driven by the
magnetorotational instability (MRI) in a collisionless plasma. Collisionless
effects may be important in radiatively inefficient accretion flows, such as
near the black hole in the Galactic Center. The MHD version of ZEUS is modified
to evolve an anisotropic pressure tensor. A fluid closure approximation is used
to calculate heat conduction along magnetic field lines. The anisotropic
pressure tensor provides a qualitatively new mechanism for transporting angular
momentum in accretion flows (in addition to the Maxwell and Reynolds stresses).
We estimate limits on the pressure anisotropy due to pitch angle scattering by
kinetic instabilities. Such instabilities provide an effective ``collision''
rate in a collisionless plasma and lead to more MHD-like dynamics. We find that
the MRI leads to efficient growth of the magnetic field in a collisionless
plasma, with saturation amplitudes comparable to those in MHD. In the saturated
state, the anisotropic stress is comparable to the Maxwell stress, implying
that the rate of angular momentum transport may be moderately enhanced in a
collisionless plasma.
[96]
oai:arXiv.org:astro-ph/0501362 [pdf] - 70474
Plasma instabilities and magnetic-field growth in clusters of galaxies
Submitted: 2005-01-17, last modified: 2005-08-03
We show that under very general conditions, cluster plasmas threaded by weak
magnetic fields are subject to very fast growing plasma instabilities driven by
the anisotropy of the plasma pressure (viscous stress) with respect to the
local direction of the magnetic field. Such an anisotropy will naturally arise
in any weakly magnetized plasma that has low collisionality and is subject to
stirring. The magnetic field must be sufficiently weak for the instabilities to
occur, viz., beta>Re^{1/2}. The instabilities are captured by the extended MHD
model with Braginskii viscosity. However, their growth rates are proportional
to the wavenumber down to the ion gyroscale, so MHD equations with Braginskii
viscosity are not well posed and a fully kinetic treatment is necessary. The
instabilities can lead to magnetic fields in clusters being amplified from seed
strength of ~10^{-18} G to dynamically important strengths of ~10 microG on
cosmologically trivial time scales (~10^8 yr). The fields produced during the
amplification stage are at scales much smaller than observed. Predicting the
saturated field scale and structure will require a kinetic theory of magnetized
cluster turbulence.
[97]
oai:arXiv.org:astro-ph/0411781 [pdf] - 69380
Magnetised plasma turbulence in clusters of galaxies
Submitted: 2004-11-29
Cluster plasmas are magnetised already at very low magnetic field strength.
Low collisionality implies that conservation of the first adiabatic invariant
results in an anisotropic viscous stress (Braginskii viscosity) or,
equivalently, anisotropic plasma pressure. This triggers firehose and mirror
instabilities, which have growth rates proportional to the wavenumber down to
scales of the order of ion Larmor radius. This means that MHD equations with
Braginskii viscosity are not well posed and fully kinetic description is
necessary. In this paper, we review the basic picture of small-scale dynamo in
the cluster plasma and attempt to reconcile it with the existence of plasma
instabilities at collisionless scales.
[98]
oai:arXiv.org:astro-ph/0305486 [pdf] - 56959
Transition from collisionless to collisional MRI
Submitted: 2003-05-25, last modified: 2003-06-13
Recent calculations by Quataert et al. (2002) found that the growth rates of
the magnetorotational instability (MRI) in a collisionless plasma can differ
significantly from those calculated using MHD. This can be important in hot
accretion flows around compact objects. In this paper we study the transition
from the collisionless kinetic regime to the collisional MHD regime, mapping
out the dependence of the MRI growth rate on collisionality. A kinetic closure
scheme for a magnetized plasma is used that includes the effect of collisions
via a BGK operator. The transition to MHD occurs as the mean free path becomes
short compared to the parallel wavelength $2\pi/k_{\Par}$. In the weak magnetic
field regime where the Alfv\'en and MRI frequencies $\omega$ are small compared
to the sound wave frequency $k_{\Par} c_0$, the dynamics are still effectively
collisionless even if $\omega \ll \nu$, so long as the collision frequency $\nu
\ll k_{\Par} c_{0}$; for an accretion flow this requires $\nu \lsim \Omega
\sqrt{\beta}$. The low collisionality regime not only modifies the MRI growth
rate, but also introduces collisionless Landau or Barnes damping of long
wavelength modes, which may be important for the nonlinear saturation of the
MRI.
