Radiative relativistic magnetohydrodynamic simulations of neutron star column accretion in Cartesian geometry
ABSTRACT

High luminosity accretion on to a strongly magnetized neutron star results in a radiation pressure dominated, magnetically confined accretion column. We investigate the dynamics of these columns using 2D radiative relativistic magnetohydrodynamic simulations, restricting consideration to modest accretion rates where the height of the column is low enough that Cartesian geometry can be employed. The column structure is dynamically maintained through high-frequency oscillations of the accretion shock at ≃ 10–25 kHz. These oscillations arise because it is necessary to redistribute the power released at the accretion shock through bulk vertical motions, both to balance the cooling and to provide vertical pressure support against gravity. Sideways cooling always dominates the loss of internal energy. In addition to the vertical oscillations, photon bubbles form in our simulations and add additional spatial complexity to the column structure. They are not themselves responsible for the oscillations, and they do not appear to affect the oscillation period. However, they enhance the vertical transport of radiation and increase the oscillation amplitude in luminosity. The time-averaged column structure in our simulations resembles the trends in standard 1D stationary models, the main difference being that the time-averaged height of the shock front is lower because of the higher more »

Authors:
; ;
Publication Date:
NSF-PAR ID:
10369592
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
515
Issue:
3
Page Range or eLocation-ID:
p. 4371-4390
ISSN:
0035-8711
Publisher:
Oxford University Press
National Science Foundation
##### More Like this
1. SUMMARY

Protracted episodes of 0.5–7 Hz pre-eruptive volcanic tremor (PVT) are common at active stratovolcanoes. Reliable links to processes related to magma movement consequently enable a potential to use properties of PVT as diagnostic eruptive precursors. A challenging feature of PVT is that generic spectral and amplitude properties of the signal evolve similarly, independent of widely varying volcano structures and conduit geometries on which most physical models rely. The ‘magma wagging’ model introduced in Jellinek & Bercovici (2011) and extended by Bercovici et al. (2013), Liao et al. and Liao & Bercovici (2018) makes progress because it depends on magma dynamics that are only weakly sensitive to volcano architecture: The flow of gas through a permeable foamy annulus of gas bubbles excites, modulates and maintains a wagging oscillation of a central magma column rising in an erupting conduit. ‘Magma wagging’ and resulting PVT are driven through an energy transfer from a ‘Bernoulli mode’ related to azimuthal variations in annular gas flow speeds. Consistent with observations, spectral and amplitude properties of PVT are predicted to evolve before an eruption as the width of the annulus decreases with increased gas fluxes. To confirm this critical Bernoulli-to-wagging energy transfer we use extensive experiments and restricted numericalmore »

2. Abstract

We explore the effect of dust on the growth of seed black holes (BHs) in the early universe. Previous 1D radiation-hydrodynamic (RHD) simulations show that increased radiation pressure on dust further suppresses the accretion rate than the case for the chemically pristine gas. Using the Enzo+Moray code, we perform a suite of 3D RHD simulations of accreting BHs in a dusty interstellar medium (ISM). We use the modified Grackle cooling library to consider dust physics in its nonequilibrium chemistry. The BH goes through an early evolutionary phase, where ionizing BH radiation creates an oscillating Hiiregion as it cycles between accretion and feedback. As the simulations proceed, dense cold gas accumulates outside the ionized region where inflow from the neutral medium meets the outflow driven by radiation pressure. In the late phase, high-density gas streams develop and break the quasi-spherical symmetry of the ionized region, rapidly boosting the accretion rate. The late phase is characterized by the coexistence of strong ionized outflows and fueling high-density gas inflows. The mean accretion rate increases with metallicity reaching a peak atZ∼ 0.01–0.1Z, one order of magnitude higher than the one for pristine gas. However, as the metallicity approaches the solar abundance, the meanmore »

3. ABSTRACT

Using 3D radiation-hydrodynamical simulations, we study the effects of ionizing radiation on the formation of second-generation (SG) stars in globular clusters (GCs) with multiple stellar populations. In particular, we focus on massive ($10^7 \, \mathrm{M}_{\odot }$) and young (40-Myr old) GCs. We consider stellar winds from asymptotic giant branch (AGB) stars, ram pressure, gas accretion on to the cluster, and photo-ionization feedback of binary stars. We find that the stellar luminosity is strong enough to warm and ionize the intracluster medium, but it does not lead to a significant gas expulsion. The cluster can thus retain the ejecta of AGB stars and the accreted pristine gas. In addition, efficient cooling occurs in the central region of the cluster within $50\, \mathrm{Myr}$ from the formation of first generation stars, leading to the formation of SG stars. Our results indicate that the inclusion of photo-ionization does not suppress SG formation, but rather delays it by about $\sim 10\, \mathrm{Myr}$. The time delay depends on the density of the pristine gas, so that a denser medium exhibits a shorter delay in star formation. Moreover, photo-ionization leads to a modest decrease in the total SG mass, compared to a model without it.

4. ABSTRACT Recent ALMA molecular line observations have revealed 3D gas velocity structure in protoplanetary discs, shedding light on mechanisms of disc accretion and structure formation. (1) By carrying out viscous simulations, we confirm that the disc’s velocity structure differs dramatically using vertical stress profiles from different accretion mechanisms. Thus, kinematic observations tracing flows at different disc heights can potentially distinguish different accretion mechanisms. On the other hand, the disc surface density evolution is mostly determined by the vertically integrated stress. The sharp disc outer edge constrained by recent kinematic observations can be caused by a radially varying α in the disc. (2) We also study kinematic signatures of a young planet by carrying out 3D planet–disc simulations. The relationship between the planet mass and the ‘kink’ velocity is derived, showing a linear relationship with little dependence on disc viscosity, but some dependence on disc height when the planet is massive (e.g. 10MJ). We predict the ‘kink’ velocities for the potential planets in DSHARP discs. At the gap edge, the azimuthally averaged velocities at different disc heights deviate from the Keplerian velocity at similar amplitudes, and its relationship with the planet mass is consistent with that in 2D simulations. After removingmore »
5. ABSTRACT

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure.more »