BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

E. Quataert

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.COastro-ph.GAastro-phphysics.plasm-phastro-ph.SRphysics.space-phgr-qcastro-ph.IMhep-phnlin.CDnucl-th

Catalog footprint

What is connected

37works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

37 published item(s)

preprint2022arXiv

Large-Scale Poloidal Magnetic Field Dynamo Leads to Powerful Jets in GRMHD Simulations of Black Hole Accretion with Toroidal Field

Accreting black holes (BHs) launch relativistic collimated jets, across many decades in luminosity and mass, suggesting the jet launching mechanism is universal, robust and scale-free. Theoretical models and general relativistic magnetohydrodynamic (GRMHD) simulations indicate that the key jet-making ingredient is large-scale poloidal magnetic flux. However, its origin is uncertain, and it is unknown if it can be generated in situ or dragged inward from the ambient medium. Here, we use the GPU-accelerated GRMHD code H-AMR to study global 3D BH accretion at unusually high resolutions more typical of local shearing box simulations. We demonstrate that turbulence in a radially-extended accretion disc can generate large-scale poloidal magnetic flux in situ, even when starting from a purely toroidal magnetic field. The flux accumulates around the BH till it becomes dynamically-important, leads to a magnetically arrested disc (MAD), and launches relativistic jets that are more powerful than the accretion flow. The jet power exceeds that of previous GRMHD toroidal field simulations by a factor of 10,000. The jets do not show significant kink or pinch instabilities, accelerate to $γ\sim 10$ over 3 decades in distance, and follow a collimation profile similar to the observed M87 jet.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Dark Energy Camera Search for an Optical Counterpart to the First Advanced LIGO Gravitational Wave Event GW150914

We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5σ$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Dark Energy Camera Search for Missing Supergiants in the LMC After the Advanced LIGO Gravitational Wave Event GW150914

The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A DECam Search for an Optical Counterpart to the LIGO Gravitational Wave Event GW151226

We report the results of a Dark Energy Camera (DECam) optical follow-up of the gravitational wave (GW) event GW151226, discovered by the Advanced LIGO detectors. Our observations cover 28.8 deg$^2$ of the localization region in the $i$ and $z$ bands (containing 3% of the BAYESTAR localization probability), starting 10 hours after the event was announced and spanning four epochs at $2-24$ days after the GW detection. We achieve $5σ$ point-source limiting magnitudes of $i\approx21.7$ and $z\approx21.5$, with a scatter of $0.4$ mag, in our difference images. Given the two day delay, we search this area for a rapidly declining optical counterpart with $\gtrsim 3σ$ significance steady decline between the first and final observations. We recover four sources that pass our selection criteria, of which three are cataloged AGN. The fourth source is offset by $5.8$ arcsec from the center of a galaxy at a distance of 187 Mpc, exhibits a rapid decline by $0.5$ mag over $4$ days, and has a red color of $i-z\approx 0.3$ mag. These properties roughly match the expectations for a kilonova. However, this source was detected several times, starting $94$ days prior to GW151226, in the Pan-STARRS Survey for Transients (dubbed as PS15cdi) and is therefore unrelated to the GW event. Given its long-term behavior, PS15cdi is likely a Type IIP supernova that transitioned out of its plateau phase during our observations, mimicking a kilonova-like behavior. We comment on the implications of this detection for contamination in future optical follow-up observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Stellar Feedback Origin for Neutral Hydrogen in High-Redshift Quasar-Mass Halos

Observations reveal that quasar host halos at z~2 have large covering fractions of cool dense gas (>~60% for Lyman limit systems within a projected virial radius). Most simulations have so far have failed to explain these large observed covering fractions. We analyze a new set of 15 simulated massive halos with explicit stellar feedback from the FIRE project, covering the halo mass range M_h~2x10^12-10^13 Msun at z=2. This extends our previous analysis of the circum-galactic medium of high-redshift galaxies to more massive halos. AGN feedback is not included in these simulations. We find Lyman limit system covering fractions consistent with those observed around quasars. The large HI covering fractions arise from star formation-driven galactic winds, including winds from low-mass satellite galaxies that interact with cosmological filaments. We show that it is necessary to resolve these satellite galaxies and their winds to reproduce the large Lyman limit system covering fractions observed in quasar-mass halos. Our simulations predict that galaxies occupying dark matter halos of mass similar to quasars but without a luminous AGN should have Lyman limit system covering fractions comparable to quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A stringent limit on the amplitude of Alfvénic perturbations in high-beta low-collisionality plasmas

It is shown that low-collisionality plasmas cannot support linearly polarized shear-Alfvén fluctuations above a critical amplitude $δB_{\perp}/B_{0} \sim β^{\,-1/2}$, where $β$ is the ratio of thermal to magnetic pressure. Above this cutoff, a developing fluctuation will generate a pressure anisotropy that is sufficient to destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, interrupting the fluctuation before any oscillation. The magnetic field lines rapidly relax into a sequence of angular zig-zag structures. Such a restrictive bound on shear-Alfvén-wave amplitudes has far-reaching implications for the physics of magnetized turbulence in the high-$β$ conditions prevalent in many astrophysical plasmas, as well as for the solar wind at $\sim 1 \mathrm{AU}$ where $β\gtrsim 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strongly Time-Variable Ultra-Violet Metal Line Emission from the Circum-Galactic Medium of High-Redshift Galaxies

We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z=2-4) galaxies. Our simulations cover the halo mass range Mh ~ 2x10^11 - 8.5x10^12 Msun at z=2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures. The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation. Spatially-extended UV metal line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutral hydrogen in galaxy halos at the peak of the cosmic star formation history

We use high-resolution cosmological zoom-in simulations from the FIRE project to make predictions for the covering fractions of neutral hydrogen around galaxies at z=2-4. These simulations resolve the interstellar medium of galaxies and explicitly implement a comprehensive set of stellar feedback mechanisms. Our simulation sample consists of 16 main halos covering the mass range M_h~10^9-6x10^12 Msun at z=2, including 12 halos in the mass range M_h~10^11-10^12 Msun corresponding to Lyman break galaxies (LBGs). We process our simulations with a ray tracing method to compute the ionization state of the gas. Galactic winds increase the HI covering fractions in galaxy halos by direct ejection of cool gas from galaxies and through interactions with gas inflowing from the intergalactic medium. Our simulations predict HI covering fractions for Lyman limit systems (LLSs) consistent with measurements around z~2-2.5 LBGs; these covering fractions are a factor ~2 higher than our previous calculations without galactic winds. The fractions of HI absorbers arising in inflows and in outflows are on average ~50% but exhibit significant time variability, ranging from ~10% to ~90%. For our most massive halos, we find a factor ~3 deficit in the LLS covering fraction relative to what is measured around quasars at z~2, suggesting that the presence of a quasar may affect the properties of halo gas on ~100 kpc scales. The predicted covering fractions, which decrease with time, peak at M_h~10^11-10^12 Msun, near the peak of the star formation efficiency in dark matter halos. In our simulations, star formation and galactic outflows are highly time dependent; HI covering fractions are also time variable but less so because they represent averages over large areas.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Impact of Baryonic Physics on the Structure of Dark Matter Halos: the View from the FIRE Cosmological Simulations

We study the distribution of cold dark matter (CDM) in cosmological simulations from the FIRE (Feedback In Realistic Environments) project, for $M_{\ast}\sim10^{4-11}\,M_{\odot}$ galaxies in $M_{\rm h}\sim10^{9-12}\,M_{\odot}$ halos. FIRE incorporates explicit stellar feedback in the multi-phase ISM, with energetics from stellar population models. We find that stellar feedback, without "fine-tuned" parameters, greatly alleviates small-scale problems in CDM. Feedback causes bursts of star formation and outflows, altering the DM distribution. As a result, the inner slope of the DM halo profile ($α$) shows a strong mass dependence: profiles are shallow at $M_{\rm h}\sim10^{10}-10^{11}\,M_{\odot}$ and steepen at higher/lower masses. The resulting core sizes and slopes are consistent with observations. This is broadly consistent with previous work using simpler feedback schemes, but we find steeper mass dependence of $α$, and relatively late growth of cores. Because the star formation efficiency $M_{\ast}/M_{\rm h}$ is strongly halo mass dependent, a rapid change in $α$ occurs around $M_{\rm h}\sim 10^{10}\,M_{\odot}$ ($M_{\ast}\sim10^{6}-10^{7}\,M_{\odot}$), as sufficient feedback energy becomes available to perturb the DM. Large cores are not established during the period of rapid growth of halos because of ongoing DM mass accumulation. Instead, cores require several bursts of star formation after the rapid buildup has completed. Stellar feedback dramatically reduces circular velocities in the inner kpc of massive dwarfs; this could be sufficient to explain the "Too Big To Fail" problem without invoking non-standard DM. Finally, feedback and baryonic contraction in Milky Way-mass halos produce DM profiles slightly shallower than the Navarro-Frenk-White profile, consistent with the normalization of the observed Tully-Fisher relation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron heat conduction in the solar wind: transition from Spitzer-Härm to the collisionless limit

We use a statistically significant set of measurements to show that the field-aligned electron heat flux $q_\parallel$ in the solar wind at 1 AU is consistent with the Spitzer-Härm collisional heat flux $q_{sh}$ for temperature gradient scales larger than a few mean free paths $L_T \gtrsim 3.5 ~λ_{fp}$. This represents about 65% of the measured data and corresponds primarily to high $β$, weakly collisional plasma ('slow solar wind'). In the more collisionless regime $λ_{fp}/L_T \gtrsim 0.28$, the electron heat flux is limited to $q_\parallel/q_0 \sim 0.3$, independent of mean free path, where $q_0$ is the 'free-streaming' value; the measured $q_\parallel$ does not achieve the full $q_0$. This constraint $q_\parallel/q_0 \sim 0.3$ might be attributed to wave-particle interactions, an interplanetary electric potential, or inherent flux limitation. We also show a $β_e$ dependence to these results that is consistent with a local radial electron temperature profile $T_e \sim r^{-α}$ that is a function of the thermal electron beta $α= α(β_e)$ and that the $β$ dependence of the collisionless regulation constraint is not obviously consistent with a whistler heat flux instability. We discuss the results in a broader astrophysical context.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Feedback-regulated star formation in molecular clouds and galactic discs

We present a two-zone theory for feedback-regulated star formation in galactic discs, consistently connecting the galaxy-averaged star formation law with star formation proceeding in giant molecular clouds (GMCs). Our focus is on galaxies with gas surface density Sigma_g>~100 Msun pc^-2. In our theory, the galactic disc consists of Toomre-mass GMCs embedded in a volume-filling ISM. Radiation pressure on dust disperses GMCs and most supernovae explode in the volume-filling medium. A galaxy-averaged star formation law is derived by balancing the momentum input from supernova feedback with the gravitational weight of the disc gas. This star formation law is in good agreement with observations for a CO conversion factor depending continuously on Sigma_g. We argue that the galaxy-averaged star formation efficiency per free fall time, epsilon_ff^gal, is only a weak function of the efficiency with which GMCs convert their gas into stars. This is possible because the rate limiting step for star formation is the rate at which GMCs form: for large efficiency of star formation in GMCs, the Toomre Q parameter obtains a value slightly above unity so that the GMC formation rate is consistent with the galaxy-averaged star formation law. We contrast our results with other theories of turbulence-regulated star formation and discuss predictions of our model. Using a compilation of data from the literature, we show that the galaxy-averaged star formation efficiency per free fall time is non-universal and increases with increasing gas fraction, as predicted by our model. We also predict that the fraction of the disc gas mass in bound GMCs decreases for increasing values of the GMC star formation efficiency. This is qualitatively consistent with the smooth molecular gas distribution inferred in local ultra-luminous infrared galaxies and the small mass fraction in giant clumps in high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic Heating, Differential Flow, and the Alpha-to-Proton Temperature Ratio in the Solar Wind

We extend previous theories of stochastic ion heating to account for the motion of ions along the magnetic field. We derive an analytic expression for the ion-to-proton perpendicular temperature ratio in the solar wind for any ion species, assuming that stochastic heating is the dominant ion heating mechanism. This expression describes how this temperature ratio depends upon the average velocity of the ions along the magnetic field direction and the ratio of the parallel proton pressure to the magnetic pressure. We compare our model with previously published measurements of alpha particles and protons from the WIND spacecraft. We find that stochastic heating offers a promising explanation for these measurements when the fractional cross helicity and Alfven ratio at the proton-gyroradius scale have values that are broadly consistent with solar-wind measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The physics of galactic winds driven by active galactic nuclei

Active galactic nuclei (AGN) drive fast winds in the interstellar medium of their host galaxies. It is commonly assumed that the high ambient densities and intense radiation fields in galactic nuclei imply short cooling times, thus making the outflows momentum-conserving. We show that cooling of high-velocity, shocked winds in AGN is in fact inefficient in a wide range of circumstances, including conditions relevant to ultra-luminous infrared galaxies (ULIRGs), resulting in energy-conserving outflows. We further show that fast energy-conserving outflows can tolerate a large amount of mixing with cooler gas before radiative losses become important. For winds with initial velocity v_in>~10,000 km s^-1, as observed in ultra-violet and X-ray absorption, the shocked wind develops a two-temperature structure. While most of the thermal pressure support is provided by the protons, the cooling processes operate directly only on the electrons. This significantly slows down inverse Compton cooling, while free free cooling is negligible. Slower winds with v_in~1,000 km s^-1, such as may be driven by radiation pressure on dust, can also experience energy-conserving phases but under more restrictive conditions. During the energy-conserving phase, the momentum flux of an outflow is boosted by a factor ~v_in/2v_s by work done by the hot post-shock gas, where v_s is the velocity of the swept-up material. Energy-conserving outflows driven by fast AGN winds (v_in~0.1c) may therefore explain the momentum fluxes Pdot>>L_AGN/c of galaxy-scale outflows recently measured in luminous quasars and ULIRGs. Shocked wind bubbles expanding normal to galactic disks may also explain the large-scale bipolar structures observed in some systems, including around the Galactic Center, and can produce significant radio, X-ray, and gamma-ray emission. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A gas cloud on its way towards the super-massive black hole in the Galactic Centre

Measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a black hole four million times the mass of the Sun. With the exception of modest X-ray and infrared flares, Sgr A* is surprisingly faint, suggesting that the accretion rate and radiation efficiency near the event horizon are currently very low. Here we report the presence of a dense gas cloud approximately three times the mass of Earth that is falling into the accretion zone of Sgr A*. Our observations tightly constrain the cloud's orbit to be highly eccentric, with an innermost radius of approach of only ~3,100 times the event horizon that will be reached in 2013. Over the past three years the cloud has begun to disrupt, probably mainly through tidal shearing arising from the black hole's gravitational force. The cloud's dynamic evolution and radiation in the next few years will probe the properties of the accretion flow and the feeding processes of the super-massive black hole. The kilo-electronvolt X-ray emission of Sgr A* may brighten significantly when the cloud reaches pericentre. There may also be a giant radiation flare several years from now if the cloud breaks up and its fragments feed gas into the central accretion zone.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Physical Model of FeLoBALs: Implications for Quasar Feedback

Photoionization modeling of the low-ionization broad absorption lines of certain quasars, known as FeLoBALs, has recently revealed the number density of the wind absorbers and their distance from the central supermassive black hole. From these, the feedback efficiency of the quasars can in principle be derived. The implied properties of the FeLoBALs are, however, surprising, with the thickness of the absorbers relative to their distance from the black hole, Delta R/R, as small as ~10^-5. Such absorbers are unlikely to survive the journey from the supermassive black hole to their inferred location. We show that the observed FeLoBAL properties are readily explained if they are formed in situ in radiative shocks produced when a quasar blast wave impacts a moderately dense interstellar clump along the line of sight. This physical picture differs significantly from the thin shell approximation often assumed, and implies outflow rates, kinetic luminosities and momentum fluxes that differ correspondingly, in some cases at the order of magnitude level. Using the radiative shock model, we estimate the ratio of the outflow kinetic luminosity to bolometric luminosity for three bright FeLoBAL quasars in the literature. We find Edot/Lbol~2-5% (and corresponding momentum fluxes Pdot~2-15 Lbol/c), similar to what is adopted in models reproducing the M-sigma relation. These outflow properties are also comparable to those recently inferred for molecular outflows in local ultra-luminous infrared galaxies, suggesting that active galactic nuclei are capable of driving such outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nickel-Rich Outflows Produced by the Accretion-Induced Collapse of White Dwarfs: Lightcurves and Spectra

The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can leave behind a rotationally supported disk with mass of up to ~ 0.1 M_sun. The disk is initially composed of free nucleons but as it accretes and spreads to larger radii, the free nucleons recombine to form helium, releasing sufficient energy to unbind the remaining disk. Most of the ejected mass fuses to form Ni56 and other iron group elements. We present spherically symmetric radiative transfer calculations of the transient powered by the radioactive heating of this ejecta. For an ejecta mass of 1e-2 M_sun (3e-3 M_sun), the lightcurve peaks after <~ 1 day with a peak bolometric luminosity ~ 2e41 erg/s (~ 5e40 erg/s), i.e., a "kilonova"; the decay time is ~ 4 (2) days. Overall, the spectra redden with time reaching U-V ~ 4 after ~ 1 day; the optical colors (B-V) are, however, somewhat blue. Near the peak in the lightcurve, the spectra are dominated by Doppler broadened Nickel features, with no distinct spectral lines present. At ~ 3-5 days, strong Calcium lines are present in the infrared, although the Calcium mass fraction is only ~ 1e-4.5. If rotationally supported disks are a common byproduct of AIC, current and upcoming transient surveys such as the Palomar Transient Factory should detect a few AIC per year for an AIC rate of ~ 1e-2 of the Type Ia rate. We discuss ways of distinguishing AIC from other rapid, faint transients, including .Ia's and the ejecta from binary neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Short GRBs with Extended Emission from Magnetar Birth: Jet Formation and Collimation

Approximately 1/4-1/2 of short duration Gamma-Ray Bursts are followed by variable X-ray emission lasting ~ 100 s with a fluence comparable or exceeding that of the initial burst itself. The long duration and significant energy of this `extended emission'(EE) poses a major challenge to the standard binary neutron star merger model. metzger08 recently proposed that the EE is powered by the spin-down of a strongly magnetized neutron star. However, the effects of surrounding material on the magnetar outflow have not yet been considered. Here we present time-dependent axisymmetric relativistic magnetohydrodynamic simulations of the interaction of the relativistic proto-magnetar wind with a surrounding 10^-1 10^-3 M_\odot envelope, which represents material ejected during the merger; the supernova following AIC; or via outflows from the accretion disk. The collision between the relativistic magnetar wind and the expanding ejecta produces a magnetized nebula inside the ejecta. A strong toroidal magnetic field builds up in the nebula, which drives a bipolar jet out through the ejecta, similar to the magnetar model developed in the case of long duration GRBs. We quantify the `break-out' time and opening angle of the jet theta_j as a function of the wind energy flux dot{E} and ejecta mass M_ej. We show that dot{E} and theta_j are inversely correlated, such that the beaming-corrected (isotropic) luminosity of the jet is primarily a function of M_ej. Both variability arguments, and the lower limit on the power of magnetar outflows capable of producing bright emission, suggest that the true opening angle of the magnetar jet must be relatively large. The model thus predicts a class of events for which the EE is observable with no associated short GRB. These may appear as long-duration GRBs or X-Ray Flashes, which may be detected by future all-sky X-ray survey missions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Proto-Magnetar Model for Gamma-Ray Bursts

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A faint type of supernova from a white dwarf with a helium-rich companion

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electromagnetic Counterparts of Compact Object Mergers Powered by the Radioactive Decay of R-process Nuclei

The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self-consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ~1e-2 M_sun[1e-3 M_sun] the resulting light curve peaks on a timescale ~ 1 day at a V-band luminosity nu L_nu ~ 3e41[1e41] ergs/s (M_V = -15[-14]); this corresponds to an effective "f" parameter ~3e-6 in the Li-Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early-time dynamics and ejecta composition. Due to the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r-process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Flares from Sgr A* and their emission mechanism

We summarize recent observations and modeling of the brightest Sgr A* flare to be observed simultaneously in (near)-infrared and X-rays to date. Trying to explain the spectral characteristics of this flare through inverse Compton mechanisms implies physical parameters that are unrealistic for Sgr A*. Instead, a "cooling break" synchrotron model provides a more feasible explanation for the X-ray emission. In a magnetic field of about 5-30 Gauss the X-ray emitting electrons cool very quickly on the typical dynamical timescale while the NIR-emitting electrons cool more slowly. This produces a spectral break in the model between NIR and X-ray wavelengths that can explain the differences in the observed spectral indices.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas

We present a theoretical framework for plasma turbulence in astrophysical plasmas (solar wind, interstellar medium, galaxy clusters, accretion disks). The key assumptions are that the turbulence is anisotropic with respect to the mean magnetic field and frequencies are low compared to the ion cyclotron frequency. The energy injected at the outer scale scale has to be converted into heat, which ultimately cannot be done without collisions. A KINETIC CASCADE develops that brings the energy to collisional scales both in space and velocity. Its nature depends on the physics of plasma fluctuations. In each of the physically distinct scale ranges, the kinetic problem is systematically reduced to a more tractable set of equations. In the "inertial range" above the ion gyroscale, the kinetic cascade splits into a cascade of Alfvenic fluctuations, which are governed by the RMHD equations at both the collisional and collisionless scales, and a passive cascade of compressive fluctuations, which obey a linear kinetic equation along the moving field lines associated with the Alfvenic component. In the "dissipation range" between the ion and electron gyroscales, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like Electron RMHD equations and a passive phase-space cascade of ion entropy fluctuations. The latter cascade brings the energy of the inertial-range fluctuations that was damped by collisionless wave-particle interaction at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for these cascades. Astrophysical and space-physical applications are discussed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Buoyancy Instabilities in Galaxy Clusters: Convection Due to Adiabatic Cosmic Rays and Anisotropic Thermal Conduction

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}/ρ^{4/3}$; $ρ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Evidence for X-ray synchrotron emission from simultaneous mid-IR to X-ray observations of a strong Sgr A* flare

This paper reports measurements of Sgr A* made with NACO in L' -band (3.80 um), Ks-band (2.12 um) and H-band (1.66 um) and with VISIR in N-band (11.88 um) at the ESO VLT, as well as with XMM-Newton at X-ray (2-10 keV) wavelengths. On 4 April, 2007, a very bright flare was observed from Sgr A* simultaneously at L'-band and X-ray wavelengths. No emission was detected using VISIR. The resulting SED has a blue slope (beta > 0 for nuL_nu ~ nu^beta, consistent with nuL_nu ~ nu^0.4) between 12 micron and 3.8 micron. For the first time our high quality data allow a detailed comparison of infrared and X-ray light curves with a resolution of a few minutes. The IR and X-ray flares are simultaneous to within 3 minutes. However the IR flare lasts significantly longer than the X-ray flare (both before and after the X-ray peak) and prominent substructures in the 3.8 micron light curve are clearly not seen in the X-ray data. From the shortest timescale variations in the L'-band lightcurve we find that the flaring region must be no more than 1.2 R_S in size. The high X-ray to infrared flux ratio, blue nuL_nu slope MIR to L' -band, and the soft nuL_nu spectral index of the X-ray flare together place strong constraints on possible flare emission mechanisms. We find that it is quantitatively difficult to explain this bright X-ray flare with inverse Compton processes. A synchrotron emission scenario from an electron distribution with a cooling break is a more viable scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fossil Gas and the Electromagnetic Precursor of Supermassive Binary Black Hole Mergers

Using a one-dimensional height integrated model, we calculate the evolution of an unequal mass binary black hole with a coplanar gas disk that contains a gap due to the presence of the secondary black hole. Viscous evolution of the outer circumbinary disk initially hardens the binary, while the inner disk drains onto the primary (central) black hole. As long as the inner disk remains cool and thin at low $\dot{M}_{\rm ext}$ (rather than becoming hot and geometrically thick), the mass of the inner disk reaches an asymptotic mass typically $\sim 10^{-3}-10^{-4}\Msun$. Once the semimajor axis shrinks below a critical value, angular momentum losses from gravitational waves dominate over viscous transport in hardening the binary. The inner disk then no longer responds viscously to the inspiraling black holes. Instead, tidal interactions with the secondary rapidly drive the inner disk into the primary. Tidal and viscous dissipation in the inner disk lead to a late time brightening in luminosity $L\propto t_{\rm minus}^{5/4}$, where $t_{\rm minus}$ is the time prior to the final merger. This late time brightening peaks $\sim 1$ day prior to the final merger at $\sim 0.1 L_{\rm Edd}$. This behavior is relatively robust because of self regulation in the coupled viscous-gravitational evolution of such binary systems. It constitutes a unique electromagnetic signature of a binary supermassive black hole merger and may allow the host galaxy to be identified if used in conjunction with the Laser Interferometric Space Antenna (LISA) localization.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

GRB 080503: Implications of a Naked Short Gamma-Ray Burst Dominated by Extended Emission

We report on observations of GRB 080503, a short gamma-ray burst with very bright extended emission (about 30 times the gamma-ray fluence of the initial spike) in conjunction with a thorough comparison to other short Swift events. In spite of the prompt-emission brightness, however, the optical counterpart is extraordinarily faint, never exceeding 25 mag in deep observations starting at ~1 hr after the BAT trigger. The optical brightness peaks at ~1 day and then falls sharply in a manner similar to the predictions of Li & Paczynski (1998) for supernova-like emission following compact-binary mergers. However, a shallow spectral index and similar evolution in X-rays inferred from Chandra observations are more consistent with an afterglow interpretation. The extreme faintness of this probable afterglow relative to the bright gamma-ray emission argues for a very low-density medium surrounding the burst (a "naked" GRB), consistent with the lack of a coincident host galaxy down to 28.5 mag in deep HST imaging. Our observations reinforce the notion that short GRBs generally occur outside regions of active star formation, but demonstrate that in some cases the luminosity of the extended prompt emission can greatly exceed that of the short spike, which may constrain theoretical interpretation of this class of events. Because most previous BAT short bursts without observed extended emission are too faint for this signature to have been detectable even if it were present at typical level, conclusions based solely on the observed presence or absence of extended emission in the existing Swift sample are premature. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetic fluctuation power near proton temperature anisotropy instability thresholds in the solar wind

The proton temperature anisotropy in the solar wind is known to be constrained by the theoretical thresholds for pressure anisotropy-driven instabilities. Here we use approximately 1 million independent measurements of gyroscale magnetic fluctuations in the solar wind to show for the first time that these fluctuations are enhanced along the temperature anisotropy thresholds of the mirror, proton oblique firehose, and ion cyclotron instabilities. In addition, the measured magnetic compressibility is enhanced at high plasma beta ($β_\parallel \gtrsim 1$) along the mirror instability threshold but small elsewhere, consistent with expectations of the mirror mode. The power in this frequency (the 'dissipation') range is often considered to be driven by the solar wind turbulent cascade, an interpretation which should be qualified in light of the present results. In particular, we show that the short wavelength magnetic fluctuation power is a strong function of collisionality, which relaxes the temperature anisotropy away from the instability conditions and reduces correspondingly the fluctuation power.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetized Relativistic Jets and Long-Duration GRBs from Magnetar Spindown during Core-Collapse Supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spindown of a newly formed millisecond, B ~ 10^{15} G, magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum, and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spindown powers (~ 10^{51}-10^{52} ergs/s), the magnetar wind is super-fast at almost all latitudes, while for lower spindown powers (~ 10^{50} erg/s), the wind is sub-fast but still super-Alfvenic. In all cases, the rates at which the neutron star loses mass, angular momentum, and energy are very similar to the corresponding free wind values (<~ 30% differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated (~5-10 deg.) relativistic jet out along the rotation axis of the star. Nearly all of the spindown power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nickel-Rich Outflows from Accretion Disks Formed by the Accretion-Induced Collapse of White Dwarfs

A white dwarf (WD) approaching the Chandrasekhar mass may in several cases undergo accretion-induced collapse (AIC) to a neutron star (NS) before a thermonuclear explosion ensues. It has generally been assumed that AIC does not produce a detectable supernova (SN). If, however, the progenitor WD is rapidly rotating (as may be expected due to its prior accretion), a centrifugally supported disk forms around the NS upon collapse. We calculate the subsequent evolution of this accretion disk using time-dependent height-integrated simulations with initial conditions taken from the AIC calculations of Dessart et al. (2006). Initially, the disk is cooled by neutrinos and its composition is driven neutron-rich (electron fraction Ye ~ 0.1) by electron captures. However, as the disk viscously spreads, it is irradiated by neutrinos from the central proto-NS, which dramatically alters its neutron-to-proton ratio. We find that electron neutrino captures increase Ye to ~ 0.5 by the time that weak interactions in the disk freeze out. Because the disk becomes radiatively inefficient and begins forming alpha-particles soon after freeze out, powerful winds blow away most of the disk's remaining mass. These Ye ~ 0.5 outflows synthesize up to a few times 1e-2 Msun in 56Ni. As a result, AIC may be accompanied by a radioactively powered SN-like transient that peaks on a timescale of ~ 1 day. Since few intermediate mass elements are likely synthesized, these Ni-rich explosions should be spectroscopically distinct from other SNe. PanSTARRs and the Palomar Transient Factory should detect a few AIC transients per year if their true rate is ~1/100 of the Type Ia rate, and LSST should detect hundreds per year. High cadence observations (< 1 day) are optimal for the detection and follow-up of AIC (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The effects of r-process heating on fall-back accretion in compact object mergers

We explore the effects of r-process nucleosynthesis on fall-back accretion in neutron star(NS)-NS and black hole-NS mergers, and the resulting implications for short-duration gamma-ray bursts (GRBs). Though dynamically important, the energy released during the r-process is not yet taken into account in merger simulations. We use a nuclear reaction network to calculate the heating (due to beta-decays and nuclear fission) experienced by material on the marginally-bound orbits nominally responsible for late-time fall-back. Since matter with longer orbital periods t_orb experiences lower densities, for longer periods of time, the total r-process heating rises rapidly with t_orb, such that material with t_orb > 1 seconds can become completely unbound. Thus, r-process heating fundamentally changes the canonical prediction of an uninterrupted power-law decline in the fall-back rate dM/dt at late times. When the timescale for r-process to complete is > 1 second, the heating produces a complete cut-off in fall-back accretion after ~ 1 second; if robust, this would imply that fall-back accretion cannot explain the late-time X-ray flaring observed following some short GRBs. However, for a narrow, but physically plausible, range of parameters, fall-back accretion can resume after ~ 10 s, despite having been strongly suppressed for ~ 1-10 s after the merger. This suggests the intriguing possibility that the gap observed between the prompt and extended emission in short GRBs is a manifestation of r-process heating.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Turbulence and Mixing in the Intracluster Medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Neutron-Rich Freeze-Out in Viscously Spreading Accretion Disks Formed from Compact Object Mergers

Accretion disks with masses ~0.001-0.1 Msun form during the merger of neutron star (NS)-NS and black hole-NS binaries. Initially, such hyper-accreting disks cool efficiently by neutrino emission and their composition is driven neutron-rich by pair captures under degenerate conditions. However, as the disk viscously spreads and its temperature drops, cooling becomes inefficient and the disk becomes advective. Analytic arguments and numerical simulations suggest that once this occurs, powerful winds likely drive away most of the disk's remaining mass. We calculate the thermal evolution and nuclear composition of viscously spreading accretion disks formed from compact object mergers using one-dimensional height-integrated simulations. We show that freeze-out from weak equilibrium necessarily accompanies the disk's late-time transition to an advective state. As a result, hyper-accreting disks generically freeze out neutron-rich (with electron fraction Ye ~ 0.2-0.4), and their late-time outflows robustly synthesize rare neutron-rich isotopes. Using the measured abundances of these isotopes in our solar system, we constrain the compact object merger rate in the Milky Way to be < 1e-5 (M_d,0/0.1 Msun)^(-1) per year, where M_d,0 is the average initial mass of the accretion disk. Thus, either the NS-NS merger rate is at the low end of current estimates or the average disk mass produced during a typical merger is << 0.1 Msun. We also show that if most short duration gamma-ray bursts (GRBs) are produced by compact object mergers, their beaming fraction must exceed f_b ~ 0.13(M_d,0/0.1 Msun), corresponding to a jet half-opening angle > 30(M_d,0/0.1 Msun)^(1/2) degrees. This is consistent with other evidence that short duration GRB outflows are less collimated than those produced in long duration GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

From Thin to Thick: The Impact of X-ray Irradiation on Accretion Disks in AGN

We argue that the X-ray and UV flux illuminating the parsec-scale accretion disk around luminous active galactic nuclei (AGN) is super-Eddington with respect to the local far-infrared dust opacity. The far infrared opacity may be larger than in the interstellar medium of the Milky Way due to a combination of supersolar metallicity and the growth of dust grains in the dense accretion disk. Because of the irradiating flux, the outer accretion disk puffs up with a vertical thickness $h\sim R$. This provides a mechanism for generating a geometrically thick obscuring region from an intrinsically thin disk. We find obscuring columns $\sim 10^{22} - 10^{23} {\rm cm}^{-2}$, in reasonable agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Photon Bubbles in the Circumstellar Envelopes of Young Massive Stars

We show that the optically-thick dusty envelopes surrounding young high-mass stars are subject to the photon bubble instability. The infrared radiation passing through the envelope amplifies magnetosonic disturbances, with growth rates in our local numerical radiation MHD calculations that are consistent with a linear analysis. Modes with wavelengths comparable to the gas pressure scale height grow by more than two orders of magnitude in a thousand years, reaching non-linear amplitudes within the envelope lifetime. If the magnetic pressure in the envelope exceeds the gas pressure, the instability develops into trains of propagating shocks. Radiation escapes readily through the low-density material between the shocks, enabling accretion to continue despite the Eddington limit imposed by the dust opacity. The supersonic motions arising from the photon bubble instability can help explain the large velocity dispersions of hot molecular cores, while conditions in the shocked gas are suitable for maser emission. We conclude that the photon bubble instability may play a key role in the formation of massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Relativistic MHD Winds from Rotating Neutron Stars

We solve the time-dependent dynamics of axisymmetric, general relativistic MHD winds from rotating neutron stars. The mass loss rate is obtained self consistently as a solution of the MHD equations, subject to a finite thermal pressure at the stellar surface. Conditions are chosen to be representative of the neutrino driven phase in newly born magnetars, which have been considered as a possible engine for GRBs. We compute the angular momentum and energy losses as a function of $σ$ and compare them with the analytic expectation from the classical theory of pulsar winds. We observe the convergence to the force-free limit in the energy loss and we study the evolution of the closed zone for increasing magnetization. Results also show that the dipolar magnetic field and the presence of a closed zone do not modify significantly the acceleration and collimation properties of the wind.

0 citations0 reviewsNo code linkedOpen access