Source author record

Christopher Thompson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2023arXiv

Utilization of Impedance Disparity Incurred from Switching Activities to Monitor and Characterize Firmware Activities

The massive trend toward embedded systems introduces new security threats to prevent. Malicious firmware makes it easier to launch cyberattacks against embedded systems. Systems infected with malicious firmware maintain the appearance of normal firmware operation but execute undesirable activities, which is usually a security risk. Traditionally, cybercriminals use malicious firmware to develop possible back-doors for future attacks. Due to the restricted resources of embedded systems, it is difficult to thwart these attacks using the majority of contemporary standard security protocols. In addition, monitoring the firmware operations using existing side channels from outside the processing unit, such as electromagnetic radiation, necessitates a complicated hardware configuration and in-depth technical understanding. In this paper, we propose a physical side channel that is formed by detecting the overall impedance changes induced by the firmware actions of a central processing unit. To demonstrate how this side channel can be exploited for detecting firmware activities, we experimentally validate it using impedance measurements to distinguish between distinct firmware operations with an accuracy of greater than 90%. These findings are the product of classifiers that are trained via machine learning. The implementation of our proposed methodology also leaves room for the use of hardware authentication.

preprint2022arXiv

Radio Emission of Pulsars. I. Slow Tearing of a Quantizing Magnetic Field

The pulsed radio emission of rotating neutron stars is connected to slow resistive instabilities feeding off an inhomogeneous twist profile within the open circuit. This paper considers the stability of a weakly sheared, quantizing magnetic field in which the current is supported by a relativistic particle flow. The electromagnetic field is almost perfectly force-free, and the particles are confined to the lowest Landau state, experiencing no appreciable curvature drift. In a charge-neutral plasma, we find multiple branches of slowly growing tearing modes, relativistic analogs of the double tearing mode, with peak growth rate $s \gtrsim 4π\widetilde k_y J_z/B_z$. Here, $B_z$ is the strong (nearly potential) guide magnetic field, $J_z$ the field-aligned current density, and $\widetilde k_y$ is the mode wavenumber normalized by the current gradient scale. These modes are overstable when the plasma carries net charge, with real frequency $ω\sim s\cdot |n_0^+ - n_0^-|/(n_0^+ + n_0^-)$ proportional to the imbalance in the densities of positive and negative charges. An isolated current sheet thinner than the skin depth supports localized tearing modes with growth rate scaling as (sheet thickness/skin depth)$^{-1/2}$. In a pulsar, the peak growth rate is comparable to the angular frequency of rotation, $s \gtrsim 2\widetilde k_y Ω$, slow compared with the longitudinal oscillations of particles and fields in a polar gap. The resistive modes experience azimuthal drift reminiscent of sub-pulse drift and are a promising driver of pulse-to-pulse flux variations. A companion paper demonstrates a Cerenkov-like instability of current-carrying Alfvén waves in thin current sheets with relativistic particle flow, and proposes coherent curvature emission by these waves as a source of pulsar radio emission.

preprint2020arXiv

Limits on Protoplanet Growth by Accretion of Small Solids

This paper identifies constraints on the growth of a small planetary core (0.3 M$_{\oplus}$) that accretes millimeter-sized pebbles from a gaseous disk. We construct time-dependent spherical envelope models that capture physical processes not included in existing global hydrodynamic simulations, including particle size evolution, dust transport, and realistic gas equations of state. We assume a low enough disk density that pebbles are marginally coupled to the gas and are trapped efficiently near the core Bondi radius. Pebbles then drift rapidly enough to experience erosion by sandblasting, mutual collisions, and sublimation of water ice. We find that pebble fragmentation is more efficient than dust re-sticking. Therefore the large pebble accretion rate $\dot{M}_p$ needed to build a core of mass > M$_{\oplus}$ leads to a high envelope metallicity and grain opacity. Above $\dot{M}_p$ $\sim 10^{-7}$ M$_{\oplus}$/yr, and without other luminosity sources, convective motions expand near the Bondi radius. The warm, dusty, and turbulent envelope buffers the inward drift of pebble debris: given a turbulent concentration factor $f_{turb} \geq$ 1 near the lower convective boundary, the core growth rate is limited to $1 \times 10^{-7} f_{turb}$ M$_{\oplus}$/yr and the e-folding time 3$f_{turb}^{-1}$ Myr. The remainder of the solid debris is expelled as highly processed silicates. Pebble ice never reaches the core, and the envelope contains comparable amounts of H2/He and metals. We interpret our results using simpler steady models and semi-analytical estimates. Future global simulations incorporating the processes modelled here are needed to understand the influence of rotation and vertical disk structure.

preprint2020arXiv

QED Phenomena in an Ultrastrong Magnetic Field. I. Electron-Photon Scattering, Pair Creation, and Annihilation

We evaluate several basic electrodynamic processes as modified by the presence of a very strong magnetic field, exceeding $B_{\rm Q} \equiv m^2/e = 4.4\times 10^{13}$ G. These results are needed to build models of dissipative phenomena outside magnetars and some other neutron stars. Differential and total cross sections and rates are presented for electron-photon scattering, the annihilation of an electron-positron pair into two photons, the inverse process of two-photon pair creation, and single-photon pair creation into the lowest Landau state. The relative importance of these interactions changes as the background magnetic field grows in strength. The particle phase space relevant for a given process may be restricted by single-photon pair creation, which also opens up efficient channels for pair multiplication, e.g. in combination with scattering. Our results are presented in the form of compact formulae that allow for relativistic electron (positron) motion, in the regime where Landau excitations can be neglected (corresponding to $10^3B_{\rm Q} \gg B \gg B_{\rm Q}$ for moderately relativistic motion along the magnetic field). Where a direct comparison is possible, our results are tested against earlier calculations, and a brief astrophysical context is provided. A companion paper considers electron-positron scattering, scattering of electrons and positrons by ions, and relativistic electron-ion bremsstrahlung.

preprint2020arXiv

QED Phenomena in an Ultrastrong Magnetic Field. II. Electron-Positron Scattering, $e^\pm$-Ion Scattering, and Relativistic Bremsstrahlung

This paper continues the approach of Kostenko \& Thompson to calculating quantum electrodynamic processes in the ultrastrong magnetic field near some neutron stars, such as magnetars or merging binary neutron stars. Here we consider electron-positron scattering, the Coulomb scattering of electrons and positrons off ions, and relativistic $e^\pm$-ion bremsstrahlung. The evaluation of differential and total cross sections simplifies considerably when the magnetic field lies in the range $10^3B_{\rm Q} \gg B \gg B_{\rm Q}$, where $B_{\rm Q} \equiv m^2/e = 4.4\times 10^{13}$ G. Then, relativistic motion of $e^\pm$ is possible even when restricted to the lowest Landau state. Accurate results for differential and total cross sections are obtained by truncating the sum over intermediate-state Landau levels and otherwise disregarding terms inversely proportional to the magnetic field, which are complicated enough to have inhibited previous attempts to calculate magnetic electron-positron scattering and relativistic bremsstrahlung. A quantitative account is made of the effects of Debye screening.

preprint2016arXiv

Determining Health Utilities through Data Mining of Social Media

'Health utilities' measure patient preferences for perfect health compared to specific unhealthy states, such as asthma, a fractured hip, or colon cancer. When integrated over time, these estimations are called quality adjusted life years (QALYs). Until now, characterizing health utilities (HUs) required detailed patient interviews or written surveys. While reliable and specific, this data remained costly due to efforts to locate, enlist and coordinate participants. Thus the scope, context and temporality of diseases examined has remained limited. Now that more than a billion people use social media, we propose a novel strategy: use natural language processing to analyze public online conversations for signals of the severity of medical conditions and correlate these to known HUs using machine learning. In this work, we filter a dataset that originally contained 2 billion tweets for relevant content on 60 diseases. Using this data, our algorithm successfully distinguished mild from severe diseases, which had previously been categorized only by traditional techniques. This represents progress towards two related applications: first, predicting HUs where such information is nonexistent; and second, (where rich HU data already exists) estimating temporal or geographic patterns of disease severity through data mining.

preprint2015arXiv

Radially Magnetized Protoplanetary Disk: Vertical Profile

This paper studies the response of a thin accretion disk to an external radial magnetic field. Our focus is on protoplanetary disks (PPDs), which are exposed during their later evolution to an intense, magnetized wind from the central star. A radial magnetic field is mixed into a thin surface layer, is wound up by the disk shear, and is pushed downward by a combination of turbulent mixing and ambipolar and Ohmic drift. The toroidal field reaches much greater strengths than the seed vertical field that is usually invoked in PPD models, even becoming superthermal. Linear stability analysis indicates that the disk experiences the magnetorotational instability (MRI) at a higher magnetization than a vertically magnetized disk when both the effects of ambipolar and Hall drift are taken into account. Steady vertical profiles of density and magnetic field are obtained at several radii between 0.06 and 1 AU in response to a wind magnetic field $B_r \sim (10^{-4}$-$10^{-2})(r/{\rm AU})^{-2}$ G. Careful attention is given to the radial and vertical ionization structure resulting from irradiation by stellar X-rays. The disk is more strongly magnetized closer to the star, where it can support a higher rate of mass transfer. As a result, the inner $\sim 1$ AU of a PPD is found to evolve toward lower surface density. Mass transfer rates around $10^{-8}\,M_\odot$ yr$^{-1}$ are obtained under conservative assumptions about the MRI-generated stress. The evolution of the disk, and the implications for planet migration, are investigated in the accompanying paper.

preprint2015arXiv

Rotation of Giant Stars

The internal rotation of post-main sequence stars is investigated, in response to the convective pumping of angular momentum toward the stellar core, combined with a tight magnetic coupling between core and envelope. The spin evolution is calculated using model stars of initial mass 1, 1.5 and $5\,M_\odot$, taking into account mass loss on the giant branches. We also include the deposition of orbital angular momentum from a sub-stellar companion, as influenced by tidal drag along with the excitation of orbital eccentricity by a fluctuating gravitational quadrupole moment. A range of angular velocity profiles $Ω(r)$ is considered in the envelope, extending from solid rotation to constant specific angular momentum. We focus on the back reaction of the Coriolis force, and the threshold for dynamo action in the inner envelope. Quantitative agreement with measurements of core rotation in subgiants and post-He core flash stars by Kepler is obtained with a two-layer angular velocity profile: uniform specific angular momentum where the Coriolis parameter ${\rm Co} \equiv Ωτ_{\rm con} \lesssim 1$ (here $τ_{\rm con}$ is the convective time); and $Ω(r) \propto r^{-1}$ where ${\rm Co} \gtrsim 1$. The inner profile is interpreted in terms of a balance between the Coriolis force and angular pressure gradients driven by radially extended convective plumes. Inward angular momentum pumping reduces the surface rotation of subgiants, and the need for a rejuvenated magnetic wind torque. The co-evolution of internal magnetic fields and rotation is considered in Kissin & Thompson, along with the breaking of the rotational coupling between core and envelope due to heavy mass loss.

preprint2015arXiv

Spin and Magnetism of White Dwarfs

The magnetism and rotation of white dwarf (WD) stars are investigated in relation to a hydromagnetic dynamo operating in the progenitor during shell burning phases. The downward pumping of angular momentum in the convective envelope, in combination with the absorption of a planet or tidal spin-up from a binary companion, can trigger strong dynamo action near the core-envelope boundary. Several arguments point to the outer core as the source for a magnetic field in the WD remnant: the outer third of a $\sim 0.55\,M_\odot$ WD is processed during the shell burning phase(s) of the progenitor; the escape of magnetic helicity through the envelope mediates the growth of (compensating) helicity in the core, as is needed to maintain a stable magnetic field in the remnant; and the intense radiation flux at the core boundary facilitates magnetic buoyancy within a relatively thick tachocline layer. The helicity flux into the growing core is driven by a dynamical imbalance with a latitude-dependent rotational stress. The magnetic field deposited in an isolated massive WD is concentrated in an outer shell of mass $\lesssim 0.1\,M_\odot$ and can reach $\sim 10\,$MG. A buried toroidal field experiences moderate ohmic decay above an age $\sim 0.3$ Gyr, which may lead to growth or decay of the external magnetic field. The final WD spin period is related to a critical spin rate below which magnetic activity shuts off, and core and envelope decouple; it generally sits in the range of hours to days. WD periods ranging up to a year are possible if the envelope re-expands following a late thermal pulse.

preprint2014arXiv

Chondrule destruction in nebular shocks

Chondrules are millimeter-sized silicate spherules ubiquitous in primitive meteorites, but whose origin remains mysterious. One of the main proposed mechanisms for producing them is melting of solids in shock waves in the gaseous protoplanetary disk. However, evidence is mounting that chondrule-forming regions were enriched in solids well above solar abundances. Given the high velocities involved in shock models destructive collisions would be expected between differently sized grains after passage of the shock front as a result of differential drag. We investigate the probability and outcome of collisions of particles behind a 1D shock using analytic methods as well as a full integration of the coupled mass, momentum, energy and radiation equations. Destruction of protochondrules seems unavoidable for solid/gas ratios $ε\gtrsim 0.1$, and possibly even for solar abundances because of "sandblasting" by finer dust. A flow with $ε\gtrsim 10$ requires much smaller shock velocities ($\sim 2$ vs 8 km s$^{-1}$) in order to achieve chondrule-melting temperatures, and radiation trapping allows slow cooling of the shocked fragments. Initial destruction would still be extensive; although re-assembly of mm-sized particles would naturally occur by grain sticking afterward, the compositional heterogeneity of chondrules may be difficult to reproduce. We finally note that solids passing through small-scale bow shocks around few-km-sized planetesimals might experience partial melting and yet escape fragmentation.

preprint2014arXiv

Non-thermal Gamma-ray Emission from Delayed Pair Breakdown in a Magnetized and Photon-rich Outflow

We consider delayed, volumetric heating in a magnetized outflow that has broken out of a confining medium and expanded to a high Lorentz factor ($Γ\sim 10^2-10^3$) and low optical depth to scattering ($τ_{\rm T} \sim 10^{-3}-10^{-2}$). The energy flux at breakout is dominated by the magnetic field, with a modest contribution from quasi-thermal gamma rays whose spectrum was calculated in Paper I. We focus on the case of extreme baryon depletion in the magnetized material, but allow for a separate baryonic component that is entrained from a confining medium. Dissipation is driven by relativistic motion between these two components, which develops once the photon compactness drops below $ 4\times 10^3(Y_e/0.5)^{-1}$. We first calculate the acceleration of the magnetized component following breakout, showing that embedded MHD turbulence provides significant inertia, the neglect of which leads to unrealistically high estimates of flow Lorentz factor. After re-heating begins, the pair and photon distributions are evolved self-consistently using a one-zone kinetic code that incorporates an exact treatment of Compton scattering, pair production and annihilation, and Coulomb scattering. Heating leads to a surge in pair creation, and the scattering depth saturates at $τ_{\rm T} \sim$ 1-4. The plasma maintains a very low ratio of particle to magnetic pressure, and can support strong anisotropy in the charged particle distribution, with cooling dominated by Compton scattering. High-energy power-law spectra with photon indices in the range observed in GRBs ($-3 < β< -3/2$) are obtained by varying the ratio of heat input to the seed energy in quasi-thermal photons. We contrast our results with those for continuous heating across an expanding photosphere, and show that the latter model produces soft-hard evolution that is inconsistent with observations of GRBs.

preprint2014arXiv

Pulse Structure of Hot Electromagnetic Outflows with Embedded Baryons

Gamma-ray bursts (GRBs) show a dramatic pulse structure that requires bulk relativistic motion, but whose physical origin has remained murky. We focus on a hot, magnetized jet that is emitted by a black hole and interacts with a confining medium. Strongly relativistic expansion of the magnetic field, as limited by a corrugation instability, may commence only after it forms a thin shell. Then the observed $T_{90}$ burst duration is dominated by the curvature delay, and null periods arise from angular inhomogeneities, not the duty cycle of the engine. We associate the $O(1)$ s timescale observed in the pulse width distribution of long GRBs with the collapse of the central 2.5-3$M_\odot$ of a massive stellar core. A fraction of the baryons are shown to be embedded in the magnetized outflow by the hyper-Eddington radiation flux; they strongly disturb the magnetic field after the compactness drops below $\sim 4\times 10^3(Y_e/0.5)^{-1}$. The high-energy photons so created have a compressed pulse structure. Delayed breakout of magnetic field from heavier baryon shells is also a promising approach to X-ray flares. In the second part of the paper, we calculate the imprint of an expanding, scattering photosphere on pulse evolution. Two models for generating the high-energy spectral tail are contrasted: i) pair breakdown due to reheating of an optically thin pair plasma embedded in a thermal radiation field; and ii) continuous heating extending from large to small scattering depth. The second model is strongly inconsistent with the observed hard-to-soft evolution in GRB pulses. The first shows some quantitative differences if the emission is purely spherical, but we show that finite shell width, mild departures from spherical curvature, and latitudinal Lorentz factor gradients have interesting effects.

preprint2013arXiv

Dissipation in Relativistic Outflows: A Multisource Overview

Relativistically expanding sources of X-rays and gamma-rays cover an enormous range of (central) compactness and Lorentz factor. The underlying physics is discussed, with an emphasis on how the dominant dissipative mode and the emergent spectrum depend on these parameters. Photons advected outward from high optical depth are a potentially important source of Compton seeds. Their characteristic energy is bounded below by ~1 MeV in pair-loaded outflows of relatively low compactness, and remains near ~1 MeV at very high compactness and low matter loading. This is compared with the characteristic energy of O(1) MeV observed in the rest frame spectra of many sources, including gamma-ray bursts, OSSE jet sources, MeV Blazars, and the intense initial 0.1 s pulse of the March 5 event. Additional topics discussed include the feedback of pair creation on electron heating and the formation of non-thermal spectra, their effectiveness at shielding the dissipative zone from ambient photons, direct Compton damping of irregularities in the outflow, the relative importance of various soft photon sources, and the softening of the emergent spectrum that results from heavy matter loading. The implications of this work for X-ray and optical afterglow from GRB's are briefly considered. Direct synchrotron emission behind the forward shock is inhibited by the extremely low energy density of the ambient magnetic field. Mildly relativistic ejecta off axis from the main gamma-ray emitting cone become optically thin to scattering on a timescale of ~1 day (E/10^{52} erg)^{1/2}, and can be a direct source of afterglow radiation.

preprint2013arXiv

Hot Electromagnetic Outflows I: Acceleration and Spectra

The theory of cold, relativistic, magnetohydrodynamic outflows is generalized by the inclusion of an intense radiation source. In some contexts, such the breakout of a gamma-ray burst jet from a star, the outflow is heated to a high temperature at a large optical depth. Eventually it becomes transparent and is pushed to a higher Lorentz factor by a combination of the Lorentz force and radiation pressure. We obtain its profile, both inside and outside the fast magnetosonic critical point, when the poloidal magnetic field is radial and monopolar. Most of the energy flux is carried by the radiation field and the toroidal magnetic field that is wound up close to the rapidly rotating engine. Although the entrained matter carries little energy, it couples the radiation field to the magnetic field. Then the fast critical point is pushed inward from infinity and, above a critical radiation intensity, the outflow is accelerated mainly by radiation pressure. We identify a distinct observational signature of this hybrid outflow: a hardening of the radiation spectrum above the peak of the seed photon distribution, driven by bulk Compton scattering. The non-thermal spectrum -- obtained by a Monte Carlo method -- is most extended when the Lorentz force dominates the acceleration, and the seed photon beam is wider than the Lorentz cone of the MHD fluid. This effect is a generic feature of hot, magnetized outflows interacting with slower relativistic material. It may explain why some GRB spectra appear to peak at photon energies above the original Amati et al. scaling. A companion paper addresses the case of jet breakout, where diverging magnetic flux surfaces yield strong MHD acceleration over a wider range of Lorentz factor.

preprint2013arXiv

Hot Electromagnetic Outflows II: Jet Breakout

We consider the interaction between radiation, matter and a magnetic field in a compact, relativistic jet. The entrained matter accelerates outward as the jet breaks out of a star or other confining medium. In some circumstances, such as gamma-ray bursts (GRBs), the magnetization of the jet is greatly reduced by an advected radiation field while the jet is optically thick to scattering. Where magnetic flux surfaces diverge rapidly, a strong outward Lorentz force develops and radiation and matter begin to decouple. The increase in magnetization is coupled to a rapid growth in Lorentz factor. We take two approaches to this problem. The first examines the flow outside the fast magnetosonic critical surface, and calculates the flow speed and the angular distribution of the radiation field over a range of scattering depths. The second considers the flow structure on both sides of the critical surface in the optically thin regime, using a relaxation method. In both approaches, we find how the terminal Lorentz factor, and radial profile of the outflow, depend on the radiation intensity and optical depth at breakout. The effect of bulk Compton scattering on the radiation spectrum is calculated by a Monte Carlo method, while neglecting the effects of internal dissipation. The peak of the scattered spectrum sits near the seed peak if radiation pressure dominates the acceleration, but is pushed to a higher frequency if the Lorentz force dominates, and especially if the seed photon cone is broadened by interaction with a slower component of the outflow.

preprint2013arXiv

Hot Electromagnetic Outflows. III. Displaced Fireball in a Strong Magnetic Field

The evolution of a dilute electron-positron fireball is calculated in the regime of strong magnetization and very high compactness (l ~10^3-10^8). Heating is applied at a low effective temperature (< 25 keV), and the fireball is allowed to expand, so that the formation of a black-body spectral distribution is inhibited by pair annihilation. The diffusion equation for Compton scattering is coupled to a single-temperature pair gas and an exact (trans-relativistic) cyclo-synchrotron photon source. We find that the photon spectrum develops a quasi-thermal peak, with a power-law slope below it that is characteristic of gamma-ray bursts. The formation of a thermal high-frequency spectrum is checked using the full kinetic equations at l ~ 10^3. These results have several implications for the central engine of GRBs, and the mechanism of energy transport. 1. Baryon rest mass carries less than ~ 10^{-5} of the energy flux at jet breakout inside ~ 10^{12} cm from the engine, with most carried by the magnetic field. 2. This degree of baryon purity points to the presence of an event horizon in the engine, and neutrons play a negligible role in the prompt emission mechanism. 3. X-ray flashes are emitted by outflows carrying enough baryons that the photosphere is pair-depleted, which we show results in faster thermalization. 4. The relation between observed peak frequency and burst luminosity is bounded below by the observed Amati et al. relation if jet Lorentz factor ~ 1/(opening angle) at breakout. 5. Stellar models are used to demonstrate an inconsistency between the highest observed GRB energies, and a hydrodynamic nozzle: magnetic collimation is required. 6. The magnetized pair gas is dilute enough that high-frequency Alfven waves may become charge starved. Finally, we suggest that limitations on magnetic reconnection from plasma collisionality have been overestimated.

preprint2012arXiv

Intense Electromagnetic Outbursts from Collapsing Hypermassive Neutron Stars

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

preprint2009arXiv

Dynamics of a Spherical Accretion Shock with Neutrino Heating and Alpha-Particle Recombination

We investigate the effects of neutrino heating and alpha-particle recombination on the hydrodynamics of core-collapse supernovae. Our focus is on the non-linear dynamics of the shock wave that forms in the collapse, and the assembly of positive energy material below it. To this end, we perform time-dependent hydrodynamic simulations with FLASH2.5 in spherical and axial symmetry. These generalize our previous calculations by allowing for bulk neutrino heating and for nuclear statistical equilibrium between n, p and alpha. The heating rate is freely tunable, as is the starting radius of the shock relative to the recombination radius of alpha-particles. An explosion in spherical symmetry involves the excitation of an overstable mode, which may be viewed as the L=0 version of the `Standing Accretion Shock Instability'. In 2D simulations, non-spherical deformations of the shock are driven by plumes of material with positive Bernoulli parameter, which are concentrated well outside the zone of strong neutrino heating. The non-spherical modes of the shock reach a large amplitude only when the heating rate is also high enough to excite convection below the shock. The critical heating rate that causes an explosion depends sensitively on the initial position of the shock relative to the recombination radius. Weaker heating is required to drive an explosion in 2D than in 1D, but the difference also depends on the size of the shock. Forcing the infalling heavy nuclei to break up into n and p below the shock only causes a slight increase in the critical heating rate, except when the shock starts out at a large radius. This shows that heating by neutrinos (or some other mechanism) must play a significant role in pushing the shock far enough out that recombination heating takes over.

preprint2009arXiv

Stability of a Spherical Accretion Shock with Nuclear Dissociation

We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability' (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.

preprint2007arXiv

Resonant Cyclotron Scattering in Three Dimensions and the Quiescent Non-thermal X-ray Emission of Magnetars

Although the surface of a magnetar is a source of bright thermal X-rays, its spectrum contains substantial non-thermal components. The X-ray emission is pulsed, with pulsed fractions that can be as high as ~ 70%. Several properties of magnetars indicate the presence of persistent, static currents flowing across the stellar surface and closing within the magnetosphere. The charges supporting these currents supply a significant optical depth to resonant cyclotron scattering in the 1-100 keV band. Here we describe a Monte Carlo approach to calculating the redistribution of thermal seed photons in frequency and angle by multiple resonant scattering in the magnetosphere. The calculation includes the full angular dependence of the cyclotron scattering cross section, the relativistic Doppler effect due to the motion of the charges, and allows for an arbitrary particle velocity distribution and magnetic field geometry. We construct synthetic spectra and pulse profiles for arbitrary orientations of the spin axis, magnetic axis, and line of sight, using a self-similar, twisted dipole field geometry, and assuming that the seed photons are supplied by single-temperature black body emission from the stellar surface. Pulse profiles and 1-10 keV spectra typical of AXPs are easily produced by this model, with pulsed fractions of ~ 50%. However, this model cannot reproduce the hard, rising energy spectra that are observed from SGRs during periods of activity, without overproducing the thermal emission peak. This suggests that the 1-100 keV emission of SGRs has a common origin with the hard X-ray emission detected from some AXPs above ~20 keV.