skip to main content

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

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 cooling efficiency of the 2D column shape.

more » « less
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 4371-4390
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We present relativistic radiation magnetohydrodynamic simulations of supercritical neutron star accretion columns in Cartesian geometry, including temperature-dependent polarization-averaged Rosseland mean opacities accounting for classical electron scattering in a magnetic field. Just as in our previous pure Thomson scattering simulations, vertical oscillations of the accretion shock and horizontally propagating entropy waves (photon bubbles) are present in all our simulations. However, at high magnetic fields ≳1012 G, the magnetic opacities produce significant differences in the overall structure and dynamics of the column. At fixed accretion rate, increasing the magnetic field strength results in a shorter accretion column, despite the fact that the overall opacity within the column is larger. Moreover, the vertical oscillation amplitude of the column is reduced. Increasing the accretion rate at high magnetic fields restores the height of the column. However, a new, slower instability takes place at these field strengths because they are in a regime where the opacity increases with temperature. This instability causes both the average height of the column and the oscillation amplitude to substantially increase on a time-scale of ∼10 ms. We provide physical explanations for these results, and discuss their implications for the observed properties of these columns, including mixed fan-beam/pencil-beam emission patterns caused by the oscillations.

    more » « less

    We perform 2D axisymmetric radiative relativistic MHD simulations of radiation pressure supported neutron star accretion columns in split-monopole magnetic fields. The accretion columns exhibit quasi-periodic oscillations, which manifest in the luminosity power spectrum as 2–10 kHz peaks, together with broader extensions to somewhat higher frequencies. The peak frequency decreases for wider columns or higher mass accretion rates. In contrast to the case of shorter columns in uniform magnetic fields, pdV work contributes substantially to maintaining the radiation pressure inside the column against sideways radiative cooling. This is in part due to the compression associated with accretion along the converging magnetic field lines towards the stellar surface. Propagating entropy waves which are associated with the slow-diffusion photon bubble instability form in all our simulations. Radial advection of radiation from the oscillation itself as well as the entropy waves is also important in maintaining radiation pressure inside the column. The time-averaged profile of our fiducial simulation accretion is approximately consistent with the classical 1D stationary model provided one incorporates the correct column shape. We also quantify the porosity in all our accretion column simulations so that this may also in principle be used to improve 1D models.

    more » « less

    We use local stratified shearing-box simulations to elucidate the impact of two-temperature thermodynamics on the thermal structure of coronae in radiatively efficient accretion flows. Rather than treating the coronal plasma as an isothermal fluid, we use a simple, parametrized cooling function that models the collisional transfer of energy from the ions to the rapidly cooling leptons. Two-temperature models naturally form temperature inversions, with a hot, magnetically dominated corona surrounding a cold disc. Simulations with net vertical flux (NF) magnetic fields launch powerful magnetocentrifugal winds that would enhance accretion in a global system. The outflow rates are much better converged with increasing box height than analogous isothermal simulations, suggesting that the winds into two-temperature coronae may be sufficiently strong to evaporate a thin disc and form a radiatively inefficient accretion flow under some conditions. We find evidence for multiphase structure in the corona, with broad density and temperature distributions, and we propose criteria for the formation of a multiphase corona. The fraction of cooling in the surface layers of the disc is substantially larger for NF fields compared to zero net-flux configurations, with moderate NF simulations radiating ≳30 per cent of the flow’s total luminosity above two mid-plane scale heights. Our work shows that NF fields may efficiently power the coronae of luminous Seyfert galaxies and quasars, providing compelling motivation for future studies of the heating mechanisms available to NF fields and the interplay of radiation with two-temperature thermodynamics.

    more » « less

    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 numerical simulations on wagging oscillations excited on analogue viscoelastic columns by annular air flows. We also explore sensitivities of the spatial and temporal characters of wagging to asymmetric annular air flows that are intractable in the existing magma wagging model and expected to occur in nature with spatial variations in annulus permeability. From high-resolution time-series of linear and orbital displacements of analogue column tops and time-series of axial deflections and accelerations of the column centre line, we characterize the excitation, evolution, and steady-state oscillations in unprecedented detail over a broad range of conditions. We show that the Bernoulli mode corresponds to the timescale for the buildup of axial elastic bending stresses in response to pressure variations related to air flows over the heights of columns. We identify three distinct wagging modes: (i) rotational (cf. Liao et al. 2018); (ii) mixed-mode and (iii) chaotic. Rotational modes are favoured for symmetric, high intensity forcing and a maximal delivery of mechanical energy to the fundamental magma wagging mode. Mixed-mode oscillations regimes are favoured for a symmetric, intermediate intensity forcing. Chaotic modes, involving the least efficient delivery of energy to the fundamental mode, occur for asymmetric forcing and where the intensity of imposed airflow is low. Numerical simulations also show that where forcing frequencies are comparable to a natural mode of free oscillation, power delivered by peripheral air flows is concentrated at the lowest frequency fundamental mode generally and spread among higher frequency natural modes where air pressure and column elastic forces are comparable. Our combined experimental and numerical results make qualitative predictions for the evolution of the character of volcanic tremor and its expression in seismic or infrasound arrays during natural events that is testable in field-based studies of PVT and syn-eruptive volcanic tremor.

    more » « less
  5. Abstract

    Recent observations have indicated significant modulation of the Madden–Julian oscillation (MJO) by the phase of the stratospheric quasi-biennial oscillation (QBO) during boreal winter. Composites of the MJO show that upper-tropospheric ice cloud fraction and water vapor anomalies are generally collocated, and that an eastward tilt with height in cloud fraction exists. Through radiative transfer calculations, it is shown that ice clouds have a stronger tropospheric radiative forcing than do water vapor anomalies, highlighting the importance of incorporating upper-tropospheric–lower-stratospheric processes into simple models of the MJO. The coupled troposphere–stratosphere linear model previously developed by the authors is extended by including a mean wind in the stratosphere and a prognostic equation for cirrus clouds, which are forced dynamically and allowed to modulate tropospheric radiative cooling, similar to the effect of tropospheric water vapor in previous formulations. Under these modifications, the model still produces a slow, eastward-propagating mode that resembles the MJO. The sign of zonal mean wind in the stratosphere is shown to control both the upward wave propagation and tropospheric vertical structure of the mode. Under varying stratospheric wind and interactive cirrus cloud radiation, the MJO-like mode has weaker growth rates under stratospheric westerlies than easterlies, consistent with the observed MJO–QBO relationship. These results are directly attributable to an enhanced barotropic mode under QBO easterlies. It is also shown that differential zonal advection of cirrus clouds leads to weaker growth rates under stratospheric westerlies than easterlies. Implications and limitations of the linear theory are discussed.

    Significance Statement

    Recent observations have shown that the strength of the Madden–Julian oscillation (MJO), a global-scale envelope of wind and rain that slowly moves eastward in the tropics and dominates global-weather variations on time scales of around a month, is strongly influenced by the direction of the winds in the lower stratosphere, the layer of the atmosphere that lies above where weather occurs. So far, modeling studies have been unable to reproduce this connection in global climate models. The purpose of this study is to investigate the mechanisms through which the stratosphere can modulate the MJO, by using simple theoretical models. In particular, we point to the role that ice clouds high in the atmosphere play in influencing the MJO.

    more » « less