Models of Hydrostatic Magnetar Atmospheres at High Luminosities

Total Page:16

File Type:pdf, Size:1020Kb

Models of Hydrostatic Magnetar Atmospheres at High Luminosities Models of hydrostatic atmospheres of magnetars at high luminosities Thijs van Putten1 Anna Watts1, Caroline D’Angelo1, Matthew Baring2, Chryssa Kouveliotou3 1University of Amsterdam, 2Rice University, 3NASA/Marshall Space Flight Center arXiv:1208.4212 October 30 2012, Fermi Symposium, Monterey image © Don Dixon Magnetars ✤ Neutron stars with inferred dipole magnetic field B ~ 1013-1016 G. ✤ Exhibit pulses (X-ray & radio), soft gamma ray bursts (~1040 erg s-1) and giant flares (~1044 erg s-1). image © NASA Magnetar atmosphere models Thijs van Putten October 30 2012 Magnetar model ✤ What is the equation of state? ✤ How and where is the emission created? ✤ What is the magnetic field configuration? Thompson & Duncan (1995) Magnetar atmosphere models Thijs van Putten October 30 2012 A peculiar magnetar burst Fermi GBM light curve of August 2008 burst from SGR 0501+4516. Light curve and black body fits of X2127 (Smale 2001) Magnetar atmosphere models Thijs van Putten October 30 2012 Photospheric Radius Expansion in magnetars? ✤ PRE in magnetars seems qualitatively possible (Watts et al. 2010) if magnetars have: ✤ Emission from optically thick region ✤ A critical luminosity ✤ Photosphere cooling with expansion ✤ Opacity increasing with radius Fermi GBM light curve of August 2008 burst from SGR 0501+4516. ✤ Observing it would constrain EoS, B and the emission location. Magnetar atmosphere models Thijs van Putten October 30 2012 Nonmagnetic models 0.9999 ✤ PRE requires sequence of extended stable atmospheres. ✤ Nonmagnetic models made by 0.9998 Paczynski & Anderson (1986). Luminosity / Eddington luminosity 0.0002 0.1 1 10 100 Atmosphere height (km) ✤ Stable nonmagnetic atmospheres exist up to r = 200 km. dP GM L r 0.0001 ⇢ 2 1 dr ' r ! Lcr ' Gas pressure / total 0.1 1 10 100 Atmosphere height (km) Magnetar atmosphere models Thijs van Putten October 30 2012 We find: No hydrostatic magnetar atmospheres with photospheric height > 10 m. Nonmagnetic case Our result Magnetar atmosphere models Thijs van Putten October 30 2012 Magnetar atmosphere models Thijs van Putten October 30 2012 Opacity in a super strong magnetic field BRIEF ARTICLE THE AUTHOR 0 Thomson scattering -2 ✤ Two photon polarizations: O-mode ] -1 g -4 (E∥B) and E-mode (E⊥B) 2 -6 10 keV -8 log Opacity [cm σO σTh 1 keV ' -10 2 2 ! T 0.1 keV -12 σE σTh 2 2 ' !C / B 11 12 13 14 15 16 log Magnetic field strength [G] E-mode Rosseland mean electron scattering opacity Magnetar atmosphere models Thijs van Putten October 30 2012 1 Atmosphere models E O 5 / P g 4 T 2 3 / B2 2 T 1 Lcr − ] ; log T [keV] P / -3 1 [g cm 0 L/Lcr =1 6 Pg / P log -1 0.01 0.1 1 10 Atmosphere height [m] Magnetar atmosphere models Thijs van Putten October 30 2012 Conclusions ✤ Hydrostatic extended magnetar atmospheres do not exist in open field line regions. ✤ Photospheric Radius Expansion cannot occur in magnetars. ✤ Spectral models of magnetar bursts considering two distinct components attributed to the E- and O-mode photospheres may have to be reconsidered. ✤ Future work: what does happen when L reaches Lcr in a magnetar atmosphere? Magnetar atmosphere models Thijs van Putten October 30 2012.
Recommended publications
  • Mass – Luminosity Relation for Massive Stars
    MASS – LUMINOSITY RELATION FOR MASSIVE STARS Within the Eddington model β ≡ Pg/P = const, and a star is an n = 3 polytrope. Large mass stars have small β, and hence are dominated by radiation pressure, and the opacity in them is dominated by electron scattering. Let us consider the outer part of such a star assuming it is in a radiative equilibrium. We have the equation of hydrostatic equilibrium: dP GM = − r ρ, (s2.1) dr r2 and the equation of radiative equilibrium dP κρL r = − r . (s2.2) dr 4πcr2 Dividing these equations side by side we obtain dP κL r = r , (s2.3) dP 4πcGMr Near the stellar surface we have Mr ≈ M and Lr ≈ L, and adopting κ ≈ κe = const, we may integrate equation (s2.3) to obtain κeLr Pr − Pr,0 = (P − P0) , (s2.4) 4πcGMr 4 where Pr,0 = P0 = aTeff /6 is the radiation pressure at the surface, i.e. at τ = 0, where the gas pressure Pg = 0. At a modest depth below stellar surface pressure is much larger than P0, we may neglect the integration constants in (s2.4) to obtain P κ L L 1+ X M⊙ L 1 − β = r = e = = , (s2.5) P 4πcGM LEdd 65300 L⊙ M where LEdd is the Eddington luminosity. Equation (s2.5) gives a relation between stellar mass, luminosity and β. The Eddington model gives a relation between stellar mass and β : 1 2 M 18.1 (1 − β) / = 2 2 , (s2.6) M⊙ µ β where µ is a mean molecular weight in units of mass of a hydrogen atom, µ−1 =2X +0.75Y +0.5Z, and X,Y,Z, are the hydrogen, helium, and heavy element abundance by mass fraction.
    [Show full text]
  • Testing Quantum Gravity with LIGO and VIRGO
    Testing Quantum Gravity with LIGO and VIRGO Marek A. Abramowicz1;2;3, Tomasz Bulik4, George F. R. Ellis5 & Maciek Wielgus1 1Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, 00-716, Warszawa, Poland 2Physics Department, Gothenburg University, 412-96 Goteborg, Sweden 3Institute of Physics, Silesian Univ. in Opava, Bezrucovo nam. 13, 746-01 Opava, Czech Republic 4Astronomical Observatory Warsaw University, 00-478 Warszawa, Poland 5Mathematics Department, University of Cape Town, Rondebosch, Cape Town 7701, South Africa We argue that if particularly powerful electromagnetic afterglows of the gravitational waves bursts detected by LIGO-VIRGO will be observed in the future, this could be used as a strong observational support for some suggested quantum alternatives for black holes (e.g., firewalls and gravastars). A universal absence of such powerful afterglows should be taken as a sug- gestive argument against such hypothetical quantum-gravity objects. If there is no matter around, an inspiral-type coalescence (merger) of two uncharged black holes with masses M1, M2 into a single black hole with the mass M3 results in emission of grav- itational waves, but no electromagnetic radiation. The maximal amount of the gravitational wave 2 energy E=c = (M1 + M2) − M3 that may be radiated from a merger was estimated by Hawking [1] from the condition that the total area of all black hole horizons cannot decrease. Because the 2 2 2 horizon area is proportional to the square of the black hole mass, one has M3 > M2 + M1 . In the 1 1 case of equal initial masses M1 = M2 = M, this yields [ ], E 1 1 h p i = [(M + M ) − M ] < 2M − 2M 2 ≈ 0:59: (1) Mc2 M 1 2 3 M From advanced numerical simulations (see, e.g., [2–4]) one gets, in the case of comparable initial masses, a much more stringent energy estimate, ! EGRAV 52 M 2 ≈ 0:03; meaning that EGRAV ≈ 1:8 × 10 [erg]: (2) Mc M The estimate assumes validity of Einstein’s general relativity.
    [Show full text]
  • Astronomy 112: the Physics of Stars Class 10 Notes: Applications and Extensions of Polytropes in the Last Class We Saw That Poly
    Astronomy 112: The Physics of Stars Class 10 Notes: Applications and Extensions of Polytropes In the last class we saw that polytropes are a simple approximation to a full solution to the stellar structure equations, which, despite their simplicity, yield important physical insight. This is particularly true for certain types of stars, such as white dwarfs and very low mass stars. In today’s class we will explore further extensions and applications of simple polytropic models, which we can in turn use to generate our first realistic models of typical main sequence stars. I. The Binding Energy of Polytropes We will begin by showing that any realistic model must have n < 5. Of course we already determined that solutions to the Lane-Emden equation reach Θ = 0 at finite ξ only for n < 5, which already hints that n ≥ 5 is a problem. Nonetheless, we have not shown that, as a matter of physical principle, models of this sort are unacceptable for stars. To demonstrate this, we will prove a generally useful result about the energy content of polytropic stars. As a preliminary to this, we will write down the polytropic relation (n+1)/n P = KP ρ in a slightly different form. Consider a polytropic star, and imagine moving down within it to the point where the pressure is larger by a small amount dP . The corre- sponding change in density dρ obeys n + 1 dP = K ρ1/ndρ. P n Similarly, since P = K ρ1/n, ρ P it follows that P ! K 1 dP ! d = P ρ(1−n)/ndρ = ρ n n + 1 ρ We can regard this relation as telling us how much the ratio of pressure to density changes when we move through a star by an amount such that the pressure alone changes by dP .
    [Show full text]
  • Ultraluminous X-Ray Pulsars
    Ultraluminous X-ray pulsars: the high-B interpretation Juri Poutanen (University of Turku, Finland; IKI, Moscow) Sergey Tsygankov, Alexander Mushtukov, Valery Suleimanov, Alexander Lutovinov, Anna Chashkina, Pavel Abolmasov Plan: • ULX before 2014 • A short history of ULXPs • Supercritical accretion onto a magnetized neutron star • Large or low B? Propeller? • Beamed? • Winds? Models for ULX (before 2014) • Super-Eddington accretion onto a stellar-mass black hole (e.g. King 2001, Begelman et al. 2006, Poutanen et al. 2007) • Sub-Eddington accretion onto intermediate mass black holes (Colbert & Mushotzky 2001) • Young rotation-powered pulsar (Medvedev & Poutanen 2013) Super-Eddington accretion • Slim disk models Accretion rate is large, but most of the released energy is advected towards the BH. ˙ ˙ M (r) =M 0 • Super-disks with winds Accretion rate is large, but most of the mass is blown away by radiation. Only the Eddington rate goes to the BH. ˙ ˙ r M (r) =M 0 Rsph ˙ ˙ M (rin ) =M Edd M˙ R R 0 - spherization radius sph ≈ in ˙ M Edd Shakura & Sunyaev 1973 Super-Eddington accretion • Slim disk models Accretion rate is large, but most of the released energy is advected towards the BH. ˙ ˙ M (r) =M 0 • Super-disks with winds Accretion rate is large, but most of the mass is blown away by radiation. Only the Eddington rate goes to the BH. ! ! r M (r) =M 0 Rsph ! ! M(rin ) =M Edd M! R ≈ R 0 -spherization radius sph S ! M Edd Ohsuga et al. 2005 Super-Eddington accretion Super-disks with winds and advection Accretion rate is large, a lot of mass is blown away by radiation (using fraction εW of available radiative flux), but still a significant fraction goes to the BH.
    [Show full text]
  • Exceeding the Eddington Limit
    The Fate of the Most Massive Stars ASP Conference Series, Vol. 332, 2005 Roberta M. Humphreys and Krzysztof Z. Stanek Exceeding the Eddington Limit Nir J. Shaviv Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel Abstract. We review the current theory describing the existence of steady state super-Eddington atmospheres. The key to the understanding of these at- mospheres is the existence of a porous layer responsible for a reduced e®ective opacity. We show how porosity arises from radiative-hydrodynamic instabili- ties and why the ensuing inhomogeneities reduce the e®ective opacity. We then discuss the appearance of these atmospheres. In particular, one of their funda- mental characteristic is the continuum driven acceleration of a thick wind from regions where the inhomogeneities become transparent. We end by discussing the role that these atmospheres play in the evolution of massive stars. 1. Introduction The notion that astrophysical objects can remain in a super-Eddington steady state for considerable amounts of time, contradicts the common wisdom usually invoked in which the classical Eddington limit, Edd, cannot be exceeded in a steady state because no hydrostatic solution canLthen exist. In other words, if objects do pass Edd, they are expected to be highly dynamic. A huge mass loss should ensueLsince large parts of their atmospheres are then gravitationally unbound and should therefore be expelled. ´ Carinae, and classical novae for that matter, stand in contrast to the above notion. ´ Car was super-Eddington (hereafter SED) during its 20 year long eruption (Davidson and Humphreys 1997). Yet, its observed mass loss and velocity are inconsistent with any standard solution (Shaviv 2000).
    [Show full text]
  • A Supercomputer Models a Blinking, Impossibly Bright 'Monster Pulsar' 8 September 2016
    A supercomputer models a blinking, impossibly bright 'monster pulsar' 8 September 2016 this case X-rays) emitted by this luminous gas are what astronomers actually observe. But as the photons travel away from the center, they push against the incoming gas, slowing the flow of gas towards the center. This force is called the radiation pressure force. As more gas falls onto the object, it becomes hotter and brighter, but if it becomes too bright the radiation pressure slows the infalling gas so much that it creates a "traffic jam." This traffic jam limits the rate at which new gas can add additional energy to the system and prevents it from getting any brighter. This luminosity upper limit, at which the radiation pressure balances the Artist’s impression of the “New Lighthouse Model.”. gravitational force, is called the Eddington Credit: NAOJ luminosity. The central energy source of enigmatic pulsating Ultra Luminous X-ray sources (ULX) could be a neutron star according to numerical simulations performed by a research group led by Tomohisa Kawashima at the National Astronomical Observatory of Japan (NAOJ). ULXs, which are remarkably bright X-ray sources, were thought to be powered by black holes. But in 2014, the X-ray space telescope "NuSTAR" detected unexpected periodic pulsed emissions in a ULX named M82 X-2. The discovery of this object named "ULX-pulsar" has puzzled The new lighthouse model (a snapshot from Movie 1) and simulation results from the present research (inset astrophysicists. Black holes can be massive on the right.) In the simulation results, the red indicates enough to provide the energy needed to create stronger radiation, and the arrows show the directions of ULXs, but black holes shouldn't be able to produce photon flow.
    [Show full text]
  • Luminous Blue Variables & Mass Loss Near the Eddington Limit
    Massive Stars as Cosmic Engines Proceedings IAU Symposium No. 250, 2008 c 2008 International Astronomical Union Fabio Bresolin, Paul Crowther & Joachim Puls, eds. DOI: 00.0000/X000000000000000X Luminous Blue Variables & Mass Loss near the Eddington Limit Stan Owocki and Allard Jan van Marle Bartol Research Institute, Department of Physics & Astronomy University of Delaware, Newark, DE 19350 USA email: [email protected], [email protected] Abstract. During the course of their evolution, massive stars lose a substantial fraction of their initial mass, both through steady winds and through relatively brief eruptions during their Luminous Blue Variable (LBV) phase. This talk reviews the dynamical driving of this mass loss, contrasting the line-driving of steady winds to the potential role of continuum driving for eruptions during LBV episodes when the star exceeds the Eddington limit. A key theme is to emphasize the inherent limits that self-shadowing places on line-driven mass loss rates, whereas continuum driving can in principle drive mass up to the “photon-tiring” limit, for which the energy to lift the wind becomes equal to the stellar luminosity. We review how the “porosity” of a highly clumped atmosphere can regulate continuum-driven mass loss, but also discuss recent time-dependent simulations of how base mass flux that exceeds the tiring limit can lead to flow stagnation and a complex, time-dependent combination of inflow and outflow regions. A general result is thus that porosity-mediated continuum driving in super-Eddington phases can explain the large, near tiring-limit mass loss inferred for LBV giant eruptions. Keywords. stars: early-type, stars: winds, outflows, stars: mass loss, stars: activity 1.
    [Show full text]
  • Supermassive Mass Black Holes
    A142: Black Holes, Winter 2012 10. Supermassive Black Holes Reading: Thorne, Chapter 9 Luminosity of an Accreting Black Hole The luminosity of an accreting black hole is proportional to the rate at which it is gaining mass. A non-spinning black hole accreting gas at a rate M˙ through a thin disk has a luminosity 1 2 L ≈ Mc˙ . 12 −1 −1 Note that M˙ is a rate of change of mass, with units of, for example, M⊙yr (or kg s ). Mc2 is energy, so Mc˙ 2 is energy/time, the units of luminosity (or power). The factor 1/12 comes from calculating the amount of gravitational potential energy gained by the gas as it moves from a large distance to the orbit at 3RSch, after which it spirals in without radiating more energy. −7 −1 For example, a non-spinning black hole accreting at M˙ = 10 M⊙yr through a thin disk has a luminosity 1 2 5 L ≈ Mc˙ = 1.2 × 10 L⊙, 12 where L⊙ is the luminosity of the sun. For a spinning black hole, the gas can get closer before it spirals in, so it can gain more energy and the factor in this equation would be bigger than 1/12. However, if there is too little “friction” in the flow of accreting gas, then it will just heat up instead of radiating the energy that it gains, and it may eventually carry that heat energy with it into the black hole. In this case we would get a “puffy” disk with lower luminosity.
    [Show full text]
  • Modified Theories and Non-Singular Black Holes
    Physics Area { PhD course in Astroparticle Physics Doctoral Thesis Beyond General Relativity: Modified Theories and Non-Singular Black Holes Candidate: Advisor: Francesco Di Filippo Prof. Stefano Liberati Academic Year 2019 { 2020 ii iii Abstract Einstein's general relativity is an extremely elegant and successful theory. Recent obser- vations coming from the LIGO/VIRGO collaboration as well as from the Event Horizon Telescope give us the possibility to perform precision test of general relativity in regimes never tested before. Even though all the current observations are in perfect agreement with the predictions of general relativity, there are several reasons to study extensions of the theory. From the experimental point of view, we are forced to include a dark sector for the matter and energy content of the universe to explain the cosmological data. Whereas from a conceptual point of view, the theory is not perturbatively renormalizable, and it predicts the formation of spacetime singularities. This thesis studies possible modifications of general relativity both considering spe- cific theories of modified gravity and implementing a model independent approach. In the first part of the thesis, we study a specific class of modified theory of gravity which has the peculiarity of propagating the same number of degrees of freedom of general relativity. The existence of these theories apparently challenges the distinctive role of general relativity as the unique non-linear theory of massless spin-2 particles. However, we provide strong evidence that these theory are actually equivalent to general relativity in vacuum. In the second part of the thesis, we focus on the problem of black hole singularities which are unavoidably present in general relativity.
    [Show full text]
  • © in This Web Service Cambridge University
    Cambridge University Press 978-1-107-03731-1 - General Relativity and Gravitation: A Centennial Perspective Edited by Abhay Ashtekar (Editor in Chief), Beverly K. Berger, James Isenberg and Malcolm Maccallum Index More information Index apparent horizons, 352, 358, 364 CDM, 166, 172, 191, 201, 206, 592 area gap, 574, 581, 586, 595 δN formalism, 189–192 area theorem, 28, 142, 614, 618 δN formula, 186 asymptotic flatness, 41, 412, 466 μ-type distortion (of CMB spectrum), 196 asymptotic quantization, 506 σ constant, 486 asymptotic safety, 504–508, 556–567 r-modes, 108 asymptotically flat initial data, 454 2dF survey, 170, 195 asymptotically flat initial data set, 418 asymptotically simple spacetimes, 422 acceleration of the Universe, 35, 166, 659 Atacama Cosmology Telescope, see ACT accretion, 135, 140, 143 averaging problem, 223, 227 disk, 28, 137, 138, 140, 155 AVTD behavior, 450 magnetized, 142 maser radiation, 137 B-mode polarization, 7, 170, 173, 182, 211, 215, 239, acoustic oscillation, 203–205 289, 334 damping, 204 back reaction, 223, 228 ACT, 163, 195, 206 background independence, 555–557, 595, 596 active galactic nuclei, 121, 134, 135, 143, 156, 157 & n-point functions, 597–600 structure resolved by microlensing, 118 & quantum geometry, 571–574 adaptive mesh refinement, 368, 372 BAO, 7–9, 170–171 ADM angular momentum, 622 Bartnik and McKinnon sequence, 377 ADM energy, 418, 419, 423 Bartnik mass, 432 ADM formulation, 18, 297–299, 364, 366, 413, 423, baryon acoustic oscillation, see BAO 433, 443, 569, 570, 593 BATSE, 148 ADM mass, 429,
    [Show full text]
  • Arxiv:1807.06243V1 [Astro-Ph.HE] 17 Jul 2018
    July 18, 2018 Formation of the First Black Holes Chapter 11: Super-Eddington accretion page 195 Chapter 11 Super-Eddington accretion; flow regimes and conditions in high-z galaxies Lucio Mayer Center for Theoretical Astrophysics and Cosmology Institute for Computational Science University of Zurich Winterthurerstrasse 190, CH-8057 Z¨urich,Switzerland [email protected] We review and discuss theoretical studies addressing the possibility of gas accretion onto black holes occurring at rates exceeding the Eddington limit. Our focus is on the applications to the growth of black hole seeds at high redshift. We first present the general notion of Super-Eddington accretion, and then summarize the different models and numerical sim- ulations developed to study such regime. We consider optically thick flows in accretion disks as well as in spherically symmetric envelopes, and devote particular attention to the widely adopted model based on the SLIM disk solution. While attractive for its simplicity, the SLIM disk solution is challenged by the latest generation of three-dimensional radi- ation (magneto)-hydrodynamical simulations, in which radiative losses can be an order of magnitude higher, and the mechanisms of radiation transport is more complex than straight advection as it takes place in a complex turbulent regime. We then discuss the gas supply rate to the sub-pc scale accretion disk or envelope from larger scales, revisiting gas inflow rates in protogalaxies under various conditions. We conclude that in the dense gaseous nuclei of high-z galaxies the conditions nec- essary for the onset of Super Eddington accretion regimes, such as a high optical depth and high gas supply rates from large scales, should be naturally met.
    [Show full text]
  • Luminosity, Selfgravitation and Nonuniqueness of Stationary Accretion
    A&A 479, 161–166 (2008) Astronomy DOI: 10.1051/0004-6361:20078212 & c ESO 2008 Astrophysics Luminosity, selfgravitation and nonuniqueness of stationary accretion J. Karkowski1, E. Malec1,2, and K. Roszkowski1 1 M. Smoluchowski Institute of Physics, Jagiellonian University, Reymonta 4, 30-059 Kraków, Poland e-mail: [email protected] 2 Max Planck Institute for Gravitational Physics (AEI), Am Mühlenberg 1, 14-476 Potsdam, Germany Received 3 July 2007 / Accepted 24 October 2007 ABSTRACT Aims. We show the existence of two branches of solutions bifurcating from a point with maximal luminosity. Methods. We investigate a Newtonian description of accreting compact bodies with hard surfaces, including luminosity and selfgrav- itation of polytropic perfect fluids. This nonlinear integro-differential problem is studied numerically. Its reduced version simplifies (under appropriate boundary conditions) to an algebraic relation between luminosity and the gas abundance in stationary, spherically symmetric flows and it can be dealt with analytically. Results. There exist – for a given luminosity, asymptotic mass and asymptotic temperature – two sub-critical solutions that bifurcate from an extremal point. They differ by the fluid content and the mass of the compact centre. Their relevance to Thorne-Zytkow˙ stars is discussed. Key words. accretion, accretion disks – hydrodynamics – gravitation – instabilities 1. Introduction We choose below the simplest self-contained model that can be interpreted as a radiating system. Our results apply primarily In this paper we study a Newtonian accreting systems consisting to Thorne-Zytkow˙ stars. We show that there is an ambiguity even of a compact body enclosed in a ball of selfgravitating gas.
    [Show full text]