Showing 1–50 of 75 results for author: Thompson, C

Quantum Plasma Creation near a Magnetar

Authors: Jonathan Zhang, Christopher Thompson

Abstract: Magnetars in quiescent states continue to emit hard X-rays with a power far exceeding the loss of rotational energy. It has recently been noted that this hard X-ray continuum may bear a direct signature of quantum electrodynamic (QED) effects in magnetic fields stronger than the Schwinger field ($B_{\rm Q} = 4.4\times 10^{13}$ G). When the current flowing into the magnetosphere is driven by narrow… ▽ More Magnetars in quiescent states continue to emit hard X-rays with a power far exceeding the loss of rotational energy. It has recently been noted that this hard X-ray continuum may bear a direct signature of quantum electrodynamic (QED) effects in magnetic fields stronger than the Schwinger field ($B_{\rm Q} = 4.4\times 10^{13}$ G). When the current flowing into the magnetosphere is driven by narrow structures in the solid crust, the $e^\pm$ pair plasma supporting the current relaxes to a collisional and trans-relativistic state. The decay of a pair into two photons produces a broad, bremsstrahlung-like spectrum of hard X-rays, similar to that observed and extending up to $0.5-1$ MeV. The conversion of two gamma rays to a pair is further enhanced by a factor $\sim B/B_{\rm Q}$. Monte Carlo calculations of pair creation in a dipole magnetic field are presented. Non-local particle injection is found to be strong enough to suppress the high voltage that otherwise would accompany polar magnetic twist; the hard X-rays are mostly emitted away from the magnetic poles. Some of the pairs annihilate in an optically thin surface layer. The prototypical anomalous X-ray pulsar 1E 2259$+$586, which shows a hard X-ray continuum but relatively weak torque noise, slow spindown, and no radio emission, is a Rosetta Stone for understanding the magnetar circuit, consistent with the picture advanced here. For a $15-60$ keV luminosity as low as $10^{34}$ erg s$^{-1}$, the polar flux of sub-relativistic pairs produces an optical depth $3-30$ to electron cyclotron scattering in the $1-10$ keV band, reducing the net X-ray polarization. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 23 pages, 16 figures, submitted to the Astrophysical Journal

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Direct Emission of Strong Radio Pulses During Magnetar Flares

Authors: Christopher Thompson

Abstract: The emission of intense radio pulses by flaring magnetars is investigated. Small-scale current gradients can be imprinted into a strongly magnetized outflow by the same processes that source fireball radiation in the closed magnetosphere. This structure arises from a combination of crustal yielding, internal tearing, and turbulent cascade. We consider the quasi-linear development of weak, small-sc… ▽ More The emission of intense radio pulses by flaring magnetars is investigated. Small-scale current gradients can be imprinted into a strongly magnetized outflow by the same processes that source fireball radiation in the closed magnetosphere. This structure arises from a combination of crustal yielding, internal tearing, and turbulent cascade. We consider the quasi-linear development of weak, small-scale currents as (i) they are stretched out and frozen by relativistic expansion and then (ii) pass through a shock. In particular, we derive the amplitudes of the ordinary and fast waves that emerge downstream of a relativistically magnetized shock in response to a mode that is frozen into the upstream flow (a frozen Alfvén wave or entropy wave). An upstream mode with comoving wavelength exceeding the skin depth can linearly convert to a secondary mode propagating above the plasma frequency. A simple and accurate treatment of shocks with extreme magnetization is developed, and the formation of internal shocks in the outflow from a bursting, rotating magnetar is outlined. The emission process described here does not require a strong shock or cool $e^\pm$ pairs (in contrast with the electromagnetic maser shock instability). In some cases, a high-frequency wave is reflected back to the observer, but with a minuscule amplitude that makes it subdominant to other emission channels. The dominant secondary electromagnetic mode is superluminal at emission, is subject to weak induced scattering within the outflow, and can reach the observer in the radio band. △ Less

Submitted 2 December, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: MNRAS in press, 24 pages, 4 figures. Expanded comparison with maser shock instability, treatment of MHD cascade in relativistic outflow

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

Measuring the stability of fundamental constants with a network of clocks

Authors: G. Barontini, L. Blackburn, V. Boyer, F. Butuc-Mayer, X. Calmet, J. R. Crespo Lopez-Urrutia, E. A. Curtis, B. Darquie, J. Dunningham, N. J. Fitch, E. M. Forgan, K. Georgiou, P. Gill, R. M. Godun, J. Goldwin, V. Guarrera, A. C. Harwood, I. R. Hill, R. J. Hendricks, M. Jeong, M. Y. H. Johnson, M. Keller, L. P. Kozhiparambil Sajith, F. Kuipers, H. S. Margolis , et al. (19 additional authors not shown)

Abstract: The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the… ▽ More The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the goal of the recently launched QSNET project: A network of clocks for measuring the stability of fundamental constants. QSNET will include state-of-the-art atomic clocks, but will also develop next-generation molecular and highly charged ion clocks with enhanced sensitivity to variations of fundamental constants. We describe the technological and scientific aims of QSNET and evaluate its expected performance. We show that in the range of parameters probed by QSNET, either we will discover new physics, or we will impose new constraints on violations of fundamental symmetries and a range of theories beyond the Standard Model, including dark matter and dark energy models. △ Less

Submitted 11 May, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Journal ref: EPJ Quantum Technology volume 9, Article number: 12 (2022)

Radio Emission of Pulsars. II. Coherence Catalyzed by Cerenkov-Unstable Shear Alfvén Waves

Authors: Christopher Thompson

Abstract: This paper explores small-scale departures from force-free electrodynamics around a rotating neutron star, extending our treatment of resistive instability in a quantizing magnetic field. A secondary, Cerenkov instability is identified: relativistic particles flowing through thin current sheets excite propagating charge perturbations that are localized near the sheets. Growth is rapid at wavenumbe… ▽ More This paper explores small-scale departures from force-free electrodynamics around a rotating neutron star, extending our treatment of resistive instability in a quantizing magnetic field. A secondary, Cerenkov instability is identified: relativistic particles flowing through thin current sheets excite propagating charge perturbations that are localized near the sheets. Growth is rapid at wavenumbers below the inverse ambient skin depth $k_{p,\rm ex}$. Small-scale Alfvénic wavepackets are promising sources of coherent curvature radiation. When the group Lorentz factor $γ_{\rm gr} \lesssim (k_{p,\rm ex}R_c)^{1/3} \sim 100$, where $R_c$ is the magnetic curvature radius, a fraction $\sim 10^{-3}$-$10^{-2}$ of the particle kinetic energy is radiated into the extraordinary mode at a peak frequency $\sim 10^{-2}ck_{p,\rm ex}$. Consistency with observations requires a high pair multiplicity ($\sim 10^{3-5}$) in the pulsar magnetosphere. Neither the primary, slow resistive instability nor the secondary, Alfvénic instability depend directly on the presence of magnetospheric `gaps', and may activate where the mean current is fully supplied by outward drift of the corotation charge. The resistive mode is overstable and grows at a rate comparable to the stellar spin frequency; the model directly accommodates strong pulse-to-pulse radio flux variations and coordinated subpulse drift. Alfvén mode growth can track the local plasma conditions, allowing for lower-frequency emission from the outer magnetosphere. Beamed radio emission from charged packets with $γ_{\rm gr} \sim 50-100$ also varies on sub-millisecond timescales. The modes identified here will be excited inside the magnetosphere of a magnetar, and may mediate Taylor relaxation of the magnetic twist. △ Less

Submitted 6 February, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 32 pages, 8 figures, Astrophysical Journal in press. Minor revisions to description of method

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

Authors: Christopher Thompson

Abstract: 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 confin… ▽ More 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. △ Less

Submitted 6 February, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 25 pages, 9 figures, Astrophysical Journal in press. Additional comparison to previous literature added; additional panel added to Fig. 1; formatting error corrected in Fig. 5

Compact Starburst Galaxies with Fast Outflows: Central Escape Velocities and Stellar Mass Surface Densities from Multi-band Hubble Space Telescope Imaging

Authors: Aleksandar M. Diamond-Stanic, John Moustakas, Paul H. Sell, Christy A. Tremonti, Alison L. Coil, Julie D. Davis, James E. Geach, Sophia C. W. Gottlieb, Ryan C. Hickox, Amanda Kepley, Charles Lipscomb, Joshua Rines, Gregory H. Rudnick, Cristopher Thompson, Kingdell Valdez, Christian Bradna, Jordan Camarillo, Eve Cinquino, Senyo Ohene Serena Perrotta, Grayson C. Petter, David S. N. Rupke, Chidubem Umeh, Kelly E. Whalen

Abstract: We present multi-band Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (0.3-1.1 $μ$m) for 12 compact starburst galaxies at z=0.4-0.8. These massive galaxies (M_stellar ~ 10^11 M_Sun) are driving very fast outflows ($v_{max}$=1000-3000 km/s), and their light profiles are dominated by an extremely compact starburst component (half-light radius ~… ▽ More We present multi-band Hubble Space Telescope imaging that spans rest-frame near-ultraviolet through near-infrared wavelengths (0.3-1.1 $μ$m) for 12 compact starburst galaxies at z=0.4-0.8. These massive galaxies (M_stellar ~ 10^11 M_Sun) are driving very fast outflows ($v_{max}$=1000-3000 km/s), and their light profiles are dominated by an extremely compact starburst component (half-light radius ~ 100 pc). Our goal is to constrain the physical mechanisms responsible for launching these fast outflows by measuring the physical conditions within the central kiloparsec. Based on our stellar population analysis, the central component typically contributes $\approx$25% of the total stellar mass and the central escape velocities $v_{esc,central}\approx900$ km/s are a factor of two smaller than the observed outflow velocities. This requires physical mechanisms that can accelerate gas to speeds significantly beyond the central escape velocities, and it makes clear that these fast outflows are capable of traveling into the circumgalactic medium, and potentially beyond. We find central stellar densities comparable to theoretical estimates of the Eddington limit, and we estimate $Σ_1$ surface densities within the central kpc comparable to those of compact massive galaxies at $0.5<z<3.0$. Relative to "red nuggets" and "blue nuggets" at $z\sim2$, we find significantly smaller $r_e$ values at a given stellar mass, which we attribute to the dominance of a young stellar component in our sample and the better physical resolution for rest-frame optical observations at $z\sim0.6$ versus $z\sim2$. We compare to theoretical scenarios involving major mergers and violent disc instability, and we speculate that our galaxies are progenitors of power-law ellipticals in the local universe with prominent stellar cusps. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Comments: accepted for publication in ApJ

Pair Plasma in Super-QED Magnetic Fields and the Hard X-ray/Optical Emission of Magnetars

Authors: Christopher Thompson, Alexander Kostenko

Abstract: The photon spectrum emitted by a transrelativistic pair plasma is calculated in the presence of an ultrastrong magnetic field, and is shown to bear a remarkable resemblance to the rising hard X-ray spectra of quiescent magnetars. This emission is powered by pair annihilation which, in contrast with a weakly magnetized pair plasma, shows an extended low-frequency tail similar to bremsstrahlung. Cro… ▽ More The photon spectrum emitted by a transrelativistic pair plasma is calculated in the presence of an ultrastrong magnetic field, and is shown to bear a remarkable resemblance to the rising hard X-ray spectra of quiescent magnetars. This emission is powered by pair annihilation which, in contrast with a weakly magnetized pair plasma, shows an extended low-frequency tail similar to bremsstrahlung. Cross sections for electron-positron annihilation/scattering, two-photon pair creation, and photon-$e^\pm$ scattering are adopted from our earlier ab initio QED calculations in the regime $10α_{\rm em}^{-1}B_{\rm Q} \gg B \gg B_{\rm Q}$. Careful attention is given to the $u$-channel scattering resonance. Magnetospheric arcades anchored in zones of intense crustal shear and extending to about twice the magnetar radius are identified as the sites of the persistent hard X-ray emission. We deduce a novel and stable configuration for the magnetospheric circuit, with a high plasma density sustained by ohmic heating and in situ pair creation. Pairs are sourced non-locally by photon collisions in zones with weak currents, such as the polar cap. Annihilation bremsstrahlung extends to the optical-IR band, where the plasma cutoff is located. The upper magnetar atmosphere experiences strong current-driven growth of ion-acoustic turbulence, which may limit positron diffusion. Coherent optical-IR emission is bounded near the observed flux by induced scattering. This model accommodates the rapid X-ray brightening of an activating magnetar, concentrated thermal hotspots, and the subdominant thermal X-ray emission of some active magnetars. Current injection is ascribed to continuous magnetic braiding, as seen in the global yielding calculations of Thompson, Yang & Ortiz. △ Less

Submitted 5 October, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 28 pages, 18 figures, accepted for publication in the Astrophysical Journal after minor revisions

Structural study of analogues of Titan s haze by trapped ion mobility coupled with a Fourier transform ion cyclotron mass spectrometer

Authors: Christopher P. Rueger, Julien Maillard, Johann Le Maitre, Mark Ridgeway, Christopher J. Thompson, Isabelle Schmitz-Afonso, Thomas Gautier, Nathalie Carrasco, Melvin A. Park, Pierre Giusti, Carlos Afonso

Abstract: The aerosols present in the atmosphere of the Saturn s moon Titan are of particular planetary science interest and several spacecraft missions already allowed to gather spectroscopic data. Titan haze s analogs, so-called tholins, were produced on earth to push forward the comprehension of their formation and properties. In this study, this highly complex mixture was analyzed here for the first tim… ▽ More The aerosols present in the atmosphere of the Saturn s moon Titan are of particular planetary science interest and several spacecraft missions already allowed to gather spectroscopic data. Titan haze s analogs, so-called tholins, were produced on earth to push forward the comprehension of their formation and properties. In this study, this highly complex mixture was analyzed here for the first time by trapped ion mobility spectrometry coupled to ultra-high resolution mass spectrometry (FTICR MS). Electrospray ionization revealed the characteristic CHNx-class components, with CHN5-6 and DBE 6-7 most abundant. Deploying specialized visualization, enabled by accurate mass measurements and elemental composition assignments, the adapted Kendrick mass defect analysis highlights the C2H3N homolog series, whereas the nitrogen-modified van Krevelen diagram exhibits a clear trend towards H/C 1.5 and N/C 0.5. More interestingly, the representation of m/z versus collision cross section (CCS) allowed hypothesizing a ramified N-PAH structural motif. State-of-the-art IMS is currently not able to resolve the isomeric continuum of ultra-complex mixtures; thus peak parameter other than the CCS-value are explored. As such, analyzing the mobility peak width versus m/z shows a linear increase in isomeric diversity between m/z 170 and 350 and a near plateau in diversity at higher m/z for the N-PAH-like structure. Due to the high complexity of the sample, these structural insights are only to be revealed by TIMS-FTICR MS. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: Journal of The American Society for Mass Spectrometry, 2019

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

Authors: Alexander Kostenko, Christopher Thompson

Abstract: 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 a… ▽ More 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. △ Less

Submitted 28 March, 2020; v1 submitted 5 April, 2019; originally announced April 2019.

Comments: 26 pages, 7 figures. Astrophysical Journal 2019, 875. Erratum added

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

Authors: Alexander Kostenko, Christopher Thompson

Abstract: 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… ▽ More 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. △ Less

Submitted 28 March, 2020; v1 submitted 5 April, 2019; originally announced April 2019.

Comments: 28 pages and 12 figures, Astrophysical Journal 2018, 869, 44. Erratum added

Limits on Protoplanet Growth by Accretion of Small Solids

Authors: Mohamad Ali-Dib, Christopher Thompson

Abstract: 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 lo… ▽ More 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. △ Less

Submitted 9 July, 2020; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 41 pages, 17 figures. Accepted for publication in the Astrophysical Journal. Expanded analytic treatment of threshold for refluxing. Corrected treatment of ice sublimation in the envelope. Effect added of heating due to the accretion of a fraction of pebble debris onto the core

Journal ref: The Astrophysical Journal, Volume 900, Number 2 (2020)

Impulsive Electromagnetic Emission near a Black Hole

Authors: Christopher Thompson

Abstract: The electromagnetic signature of a point explosion near a Kerr black hole (BH) is evaluated. The first repetitions produced by gravitational lensing are not periodic in time; periodicity emerges only as the result of multiple circuits of the prograde and retrograde light rings and is accompanied by exponential dimming. Gravitational focusing creates a sequence of concentrated caustic features and… ▽ More The electromagnetic signature of a point explosion near a Kerr black hole (BH) is evaluated. The first repetitions produced by gravitational lensing are not periodic in time; periodicity emerges only as the result of multiple circuits of the prograde and retrograde light rings and is accompanied by exponential dimming. Gravitational focusing creates a sequence of concentrated caustic features and biases the detection of a repeating source toward alignment of the BH spin with the plane of the sky. We consider the polarization pattern in the case of emission by the Lorentz upboosting and reflection of a magnetic field near the explosion site. Then the polarized fraction of the detected pulse approaches unity, and rays propagating near the equatorial plane maintain a consistent polarization direction. Near a slowly accreting supermassive BH (SMBH), additional repetitions are caused by reflection off annular fragments of an orbiting disk that has passed through an ionization instability. These results are applied to the repeating fast radio burst (FRB) source 121102, giving a concrete and predictive example of how FRB detectability may be biased by lensing. A gravitational lensing delay of 10-30 s, and reflection delay up to $\sim 10^4$ s, are found for emission near the innermost stable circular orbit of a $3\times 10^5\,M_\odot$ SMBH; these effects combine to produce interesting correlations between delay time and burst fluence. A similar repetitive pulse envelope could be seen in the gravitational wave signal produced by a collision between compact stars near a SMBH. △ Less

Submitted 7 April, 2019; v1 submitted 30 October, 2018; originally announced October 2018.

Comments: Published in the Astrophysical Journal, 32 pages, 17 figures

Journal ref: ApJ 2019, 874, 48

Rotation and Magnetism of Massive Stellar Cores

Authors: Yevgeni Kissin, Christopher Thompson

Abstract: The internal rotation and magnetism of evolved massive stars are considered in response to i) the inward pumping of angular momentum through deep and slowly rotating convective layers; and ii) the winding up of a helical magnetic field in radiative layers. Field winding can transport angular momentum effectively even when the toroidal field is limited by kinking. Magnetic helicity is pumped into a… ▽ More The internal rotation and magnetism of evolved massive stars are considered in response to i) the inward pumping of angular momentum through deep and slowly rotating convective layers; and ii) the winding up of a helical magnetic field in radiative layers. Field winding can transport angular momentum effectively even when the toroidal field is limited by kinking. Magnetic helicity is pumped into a growing radiative layer from an adjacent convective envelope (or core). The receding convective envelope that forms during the early accretion phase of a massive star is the dominant source of helicity in its core, yielding a $\sim 10^{13}$ G polar magnetic field in a collapsed neutron star (NS) remnant. Using MESA models of various masses, we find that the NS rotation varies significantly, from $P_{\rm NS} \sim 0.1-1$ s in a 13$\,M_\odot$ model to $P_{\rm NS} \sim 2$ ms in a $25\,M_\odot$ model with an extended core. Stronger inward pumping of angular momentum is found in more massive stars, due to the growing thickness of the convective shells that form during the later stages of thermonuclear burning. On the other hand, stars that lose enough mass to form blue supergiants in isolation end up as very slow rotators. The tidal spin-up of a 40$\,M_\odot$ star by a massive binary companion is found to dramatically increase the spin of the remnant black hole, allowing a rotationally supported torus to form during the collapse. The implications for post-collapse decay or amplification of the magnetic field are also considered. △ Less

Submitted 26 February, 2018; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 21 pages, 23 figures, minor revisions including expanded comparison with previous work, Astrophysical Journal, in press

Giant Primeval Magnetic Dipoles

Authors: Christopher Thompson

Abstract: Macroscopic magnetic dipoles are considered as cosmic dark matter. Permanent magnetism in relativistic field structures can involve some form of superconductivity, one example being current-carrying string loops (`springs') with vanishing net tension. We derive the cross section for free classical dipoles to collide, finding it depends weakly on orientation when mutual precession is rapid. The col… ▽ More Macroscopic magnetic dipoles are considered as cosmic dark matter. Permanent magnetism in relativistic field structures can involve some form of superconductivity, one example being current-carrying string loops (`springs') with vanishing net tension. We derive the cross section for free classical dipoles to collide, finding it depends weakly on orientation when mutual precession is rapid. The collision rate of `spring' loops with tension ${\cal T} \sim 10^{-8}c^4/G$ in galactic halos approaches the measured rate of fast radio bursts (FRBs) if the loops comprise most of the dark matter. A large superconducting dipole (LSD) with mass $\sim 10^{20}$ g and size $\sim 1$ mm will form a $\sim 100$ km magnetosphere moving through interstellar plasma. Although hydromagnetic drag is generally weak, it is strong enough to capture some LSDs into long-lived rings orbiting supermassive black holes (SMBHs) that form by the direct collapse of massive gas clouds. Repeated collisions near young SMBHs could dominate the global collision rate, thereby broadening the dipole mass spectrum. Colliding LSDs produce tiny, hot electromagnetic explosions. The accompanying paper shows that these explosions couple effectively to propagating low-frequency electromagnetic modes, with output peaking at 0.01-1 THz. We describe several constraints on, and predictions of, LSDs as cosmic dark matter. The shock formed by an infalling LSD triggers self-sustained thermonuclear burning in a C/O (ONeMg) white dwarf (WD) of mass $\gtrsim 1\,M_\odot$ ($1.3\,M_\odot$). The spark is generally located well off the center of the WD. The rate of LSD-induced explosions matches the observed rate of Type Ia supernovae. △ Less

Submitted 29 May, 2017; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: 23 pages, 19 figures, Astrophysical Journal in press

Tiny Electromagnetic Explosions

Authors: Christopher Thompson

Abstract: This paper considers electromagnetic transients of a modest total energy (${\cal E} \gtrsim 10^{40}$ erg) and small initial size (${\cal R} \gtrsim 10^{-1}$ cm). They could be produced during collisions between relativistic field structures (e.g. macroscopic magnetic dipoles) that formed around, or before, cosmic electroweak symmetry breaking. The outflowing energy has a dominant electromagnetic c… ▽ More This paper considers electromagnetic transients of a modest total energy (${\cal E} \gtrsim 10^{40}$ erg) and small initial size (${\cal R} \gtrsim 10^{-1}$ cm). They could be produced during collisions between relativistic field structures (e.g. macroscopic magnetic dipoles) that formed around, or before, cosmic electroweak symmetry breaking. The outflowing energy has a dominant electromagnetic component; a subdominant thermal component (temperature $> 1$ GeV) supplies inertia in the form of residual $e^\pm$. A thin shell forms that expands subluminally, attaining a Lorentz factor $\sim 10^{6-7}$ before decelerating. Drag is supplied by the reflection of an ambient magnetic field, and by deflection of ambient free electrons. Emission of low-frequency (GHz-THz) superluminal waves takes place through three channels: i) reflection of the ambient magnetic field; ii) direct linear conversion of the embedded magnetic field into a superluminal mode; and iii) excitation outside the shell by corrugation of its surface. The escaping electromagnetic pulse is very narrow (a few wavelengths) and so the width of the detected transient is dominated by propagation effects. GHz radio transients are emitted from i) the dark matter halos of galaxies and ii) the near-horizon regions of supermassive black holes that formed by direct gas collapse and now accrete slowly. Brighter and much narrower 0.01-1 THz pulses are predicted at a rate at least comparable to fast radio bursts, experiencing weaker scattering and absorption. The same explosions also accelerate protons up to $\sim 10^{19}$ eV and heavier nuclei up to $10^{20-21}$ eV. △ Less

Submitted 6 June, 2017; v1 submitted 1 March, 2017; originally announced March 2017.

Comments: 25 pages, 16 figures, Astrophysical Journal, in press

Global Crustal Dynamics of Magnetars in Relation to their Bright X-ray Outbursts

Authors: Christopher Thompson, Huan Yang, Néstor Ortiz

Abstract: This paper considers the yielding response of a neutron star crust to smooth, unbalanced Maxwell stresses imposed at the core-crust boundary, and the coupling of the dynamic crust to the external magnetic field. Stress buildup and yielding in a magnetar crust is a global phenomenon: an elastic distortion radiating from one plastically deforming zone is shown to dramatically increase the creep rate… ▽ More This paper considers the yielding response of a neutron star crust to smooth, unbalanced Maxwell stresses imposed at the core-crust boundary, and the coupling of the dynamic crust to the external magnetic field. Stress buildup and yielding in a magnetar crust is a global phenomenon: an elastic distortion radiating from one plastically deforming zone is shown to dramatically increase the creep rate in distant zones. Runaway creep to dynamical rates is shown to be possible, being enhanced by in situ heating and suppressed by thermal conduction and shearing of an embedded magnetic field. A global and time-dependent model of elastic, plastic, magnetic, and thermal evolution is developed. Fault-like structures develop naturally, and a range of outburst timescales is observed. Transient events with time profiles similar to giant magnetar flares (millisecond rise, $\sim$ 0.1 s duration, and decaying power-law tails) result from runaway creep that starts in localized sub-km-sized patches and spreads across the crust. A one-dimensional model of stress relaxation in the vertically stratified crust shows that a modest increase in applied stress allows embedded magnetic shear to escape the star over $\sim$ 3-10 ms, dissipating greater energy if the exterior field is already sheared. Several such zones coupled to each other naturally yield a burst of duration $\sim$ 0.1 s, as is observed over a wide range of burst energies. The collective interaction of many plastic zones forces an overstability of global elastic modes of the crust, consistent with QPO activity extending over $\gtrsim$ 100 s. Giant flares probably involve sudden meltdown in localized zones, with high-frequency ($\gg$ 100 Hz) QPOs corresponding to standing Alfvén waves within these zones. △ Less

Submitted 8 August, 2016; originally announced August 2016.

Comments: 25 pages, 20 figures, submitted to the Astrophysical Journal

Journal ref: ApJ 841, 54, 2017

Radially Magnetized Protoplanetary Disk: Vertical Profile

Authors: Matthew Russo, Christopher Thompson

Abstract: 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 ambi… ▽ More 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. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: 25 pp, 20 figures, Astrophysical Journal, in press

Constrained Evolution of a Radially Magnetized Protoplanetary Disk: Implications for Planetary Migration

Authors: Matthew Russo, Christopher Thompson

Abstract: We consider the inner $\sim$ AU of a protoplanetary disk (PPD), at a stage where angular momentum transport is driven by the mixing of a radial magnetic field into the disk from a T-Tauri wind. Because the radial profile of the imposed magnetic field is well constrained, a deterministic calculation of the disk mass flow becomes possible. The vertical disk profiles obtained in Paper I imply a stron… ▽ More We consider the inner $\sim$ AU of a protoplanetary disk (PPD), at a stage where angular momentum transport is driven by the mixing of a radial magnetic field into the disk from a T-Tauri wind. Because the radial profile of the imposed magnetic field is well constrained, a deterministic calculation of the disk mass flow becomes possible. The vertical disk profiles obtained in Paper I imply a stronger magnetization in the inner disk, faster accretion, and a secular depletion of the disk material. Inward transport of solids allows the disk to maintain a broad optical absorption layer even when the grain abundance becomes too small to suppress its ionization. Thus a PPD may show a strong middle-to-near infrared spectral excess even while its mass profile departs radically from the minimum-mass solar nebula. The disk surface density is buffered at $\sim 30$ g cm$^{-2}$: below this, X-rays trigger strong enough magnetorotational turbulence at the midplane to loft mm-cm sized particles high in the disk, followed by catastrophic fragmentation. A sharp density gradient bounds the inner depleted disk, and propagates outward to $\sim 1$-2 AU over a few Myr. Earth-mass planets migrate through the inner disk over a similar timescale, whereas the migration of Jupiters is limited by the supply of gas. Gas-mediated migration must stall outside 0.04 AU, where silicates are sublimated and the disk shifts to a much lower column. A transition disk emerges when the dust/gas ratio in the MRI-active layer falls below $X_d \sim 10^{-6}(a_d/μ{\rm m})$, where $a_d$ is the grain size. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: 22 pp, 18 figures, Astrophysical Journal, in press

Rotation of Giant Stars

Authors: Yevgeni Kissin, Christopher Thompson

Abstract: 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 orbi… ▽ More 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. △ Less

Submitted 9 June, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 15 pages, 21 figures, ApJ in press, figure colours improved, minor corrections

Spin and Magnetism of White Dwarfs

Authors: Yevgeni Kissin, Christopher Thompson

Abstract: 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 arg… ▽ More 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. △ Less

Submitted 31 July, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 23 pages, 27 figures, accepted to the ApJ, typos and figures corrected

Chondrule destruction in nebular shocks

Authors: Emmanuel Jacquet, Christopher Thompson

Abstract: 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… ▽ More 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. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 18 pages, 12 figures. Accepted to the Astrophysical Journal

Pulse Structure of Hot Electromagnetic Outflows with Embedded Baryons

Authors: Christopher Thompson, Ramandeep Gill

Abstract: 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 o… ▽ More 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. △ Less

Submitted 19 June, 2014; originally announced June 2014.

Comments: Submitted to the ApJ, 21 pages, 22 figures

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

Authors: Ramandeep Gill, Christopher Thompson

Abstract: 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. W… ▽ More 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. △ Less

Submitted 18 June, 2014; originally announced June 2014.

Comments: Submitted to the ApJ, 27 pages, 19 figures

Dissipation in Relativistic Outflows: A Multisource Overview

Authors: Christopher Thompson

Abstract: 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 ene… ▽ More 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. △ Less

Submitted 9 October, 2013; originally announced October 2013.

Comments: 23 pages, August 1997, invited review, Cracow Conference on Relativistic Jets in AGN, eds. M. Ostrowski, M. Sikora, G. Madejski, and M. Begelman, Jagellonian University Press, p. 63 [cited in accompanying astro-ph preprint, but not available on NASA ADS]

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

Authors: Christopher Thompson, Ramandeep Gill

Abstract: 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… ▽ More 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. △ Less

Submitted 9 October, 2013; originally announced October 2013.

Comments: 29 pages, 34 figures, submitted to the ApJ

Hot Electromagnetic Outflows II: Jet Breakout

Authors: Matthew Russo, Christopher Thompson

Abstract: 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 magn… ▽ More 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. △ Less

Submitted 7 March, 2013; v1 submitted 6 March, 2013; originally announced March 2013.

Comments: Submitted to the Astrophysical Journal; revised in response to the referee's comments [v1 incorrectly duplicated Figure 8a]

Hot Electromagnetic Outflows I: Acceleration and Spectra

Authors: Matthew Russo, Christopher Thompson

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 6 March, 2013; originally announced March 2013.

Comments: To be published in the Astrophysical Journal

The Pulsar Search Collaboratory: Discovery and Timing of Five New Pulsars

Authors: R. Rosen, J. Swiggum, M. A. McLaughlin, D. R. Lorimer, M. Yun, S. Heatherly, J. Boyles, R. Lynch, V. I. Kondratiev, S. Scoles, S. M. Ransom, M. L. Moniot, A. Cottrill, M. Weaver, A. Snider, C. Thompson, M. Raycraft, J. Dudenhoefer, L. Allphin, J. Thorley, B. Meadows, G. Marchiny, A. Liska, A. M. O'Dwyer, B. Butler , et al. (46 additional authors not shown)

Abstract: We present the discovery and timing solutions of five new pulsars by students involved in the Pulsar Search Collaboratory (PSC), a NSF-funded joint program between the National Radio Astronomy Observatory and West Virginia University designed to excite and engage high-school students in Science, Technology, Engineering, and Mathematics (STEM) and related fields. We encourage students to pursue STE… ▽ More We present the discovery and timing solutions of five new pulsars by students involved in the Pulsar Search Collaboratory (PSC), a NSF-funded joint program between the National Radio Astronomy Observatory and West Virginia University designed to excite and engage high-school students in Science, Technology, Engineering, and Mathematics (STEM) and related fields. We encourage students to pursue STEM fields by apprenticing them within a professional scientific community doing cutting edge research, specifically by teaching them to search for pulsars. The students are analyzing 300 hours of drift-scan survey data taken with the Green Bank Telescope at 350 MHz. These data cover 2876 square degrees of the sky. Over the course of five years, more than 700 students have inspected diagnostic plots through a web-based graphical interface designed for this project. The five pulsars discovered in the data have spin periods ranging from 3.1 ms to 4.8 s. Among the new discoveries are - PSR J1926-1314, a long period, nulling pulsar; PSR J1821+0155, an isolated, partially recycled 33-ms pulsar; and PSR J1400-1438, a millisecond pulsar in a 9.5-day orbit whose companion is likely a white dwarf star. △ Less

Submitted 14 February, 2013; v1 submitted 18 September, 2012; originally announced September 2012.

Journal ref: 2013 ApJ, 768, 85

Intense Electromagnetic Outbursts from Collapsing Hypermassive Neutron Stars

Authors: Luis Lehner, Carlos Palenzuela, Steven L. Liebling, Christopher Thompson, Chad Hanna

Abstract: 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 ex… ▽ More 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. △ Less

Submitted 5 September, 2012; v1 submitted 12 December, 2011; originally announced December 2011.

Comments: 28 pages, 20 figures (quality lowered to reduce sizes). Improved initial data and matching condition results in a lower, but still important, energy emission. Added appendix with a discussion on effects of transition layer

Journal ref: Phys. Rev. D86,104035 (2012)

Teraflop per second gravitational lensing ray-shooting using graphics processing units

Authors: Alexander C. Thompson, Christopher J. Fluke, David G. Barnes, Benjamin R. Barsdell

Abstract: Gravitational lensing calculation using a direct inverse ray-shooting approach is a computationally expensive way to determine magnification maps, caustic patterns, and light-curves (e.g. as a function of source profile and size). However, as an easily parallelisable calculation, gravitational ray-shooting can be accelerated using programmable graphics processing units (GPUs). We present our imp… ▽ More Gravitational lensing calculation using a direct inverse ray-shooting approach is a computationally expensive way to determine magnification maps, caustic patterns, and light-curves (e.g. as a function of source profile and size). However, as an easily parallelisable calculation, gravitational ray-shooting can be accelerated using programmable graphics processing units (GPUs). We present our implementation of inverse ray-shooting for the NVIDIA G80 generation of graphics processors using the NVIDIA Compute Unified Device Architecture (CUDA) software development kit. We also extend our code to multiple-GPU systems, including a 4-GPU NVIDIA S1070 Tesla unit. We achieve sustained processing performance of 182 Gflop/s on a single GPU, and 1.28 Tflop/s using the Tesla unit. We demonstrate that billion-lens microlensing simulations can be run on a single computer with a Tesla unit in timescales of order a day without the use of a hierarchical tree code. △ Less

Submitted 14 May, 2009; originally announced May 2009.

Comments: 21 pages, 4 figures, submitted to New Astronomy

Journal ref: New Astron.15:16-23,2010

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

Authors: Rodrigo Fernández, Christopher Thompson

Abstract: 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… ▽ More 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. △ Less

Submitted 20 September, 2009; v1 submitted 30 December, 2008; originally announced December 2008.

Comments: Published in ApJ. This version includes some proof corrections. Minor changes relative to submitted version (rephrasing of some statements, and updated figures 5, 6a, 9, 10, 12, 13, 14, 17). Version with high resolution figures at http://sns.ias.edu/~rafernan/alpha_movies.html

Journal ref: Astrophys.J.703:1464-1485,2009

Stability of a Spherical Accretion Shock with Nuclear Dissociation

Authors: Rodrigo Fernández, Christopher Thompson

Abstract: 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 us… ▽ More 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. △ Less

Submitted 1 March, 2009; v1 submitted 17 November, 2008; originally announced November 2008.

Comments: Version accepted by ApJ. Consistent parameter range, and additional plots with eigenfrequencies in units of inverse postshock advection timescale. Minor changes otherwise

Journal ref: Astrophys.J.697:1827-1841,2009

Electrodynamics of Magnetars IV: Self-Consistent Model of the Inner Accelerator, with Implications for Pulsed Radio Emission

Authors: Christopher Thompson

Abstract: We consider the voltage structure in the open-field circuit and outer magnetosphere of a magnetar. The standard polar-cap model for radio pulsars is modified significantly when the polar magnetic field exceeds 1.8x10^{14} G. Pairs are created by accelerated particles via resonant scattering of thermal X-rays, followed by the nearly instantaneous conversion of the scattered photon to a pair. A su… ▽ More We consider the voltage structure in the open-field circuit and outer magnetosphere of a magnetar. The standard polar-cap model for radio pulsars is modified significantly when the polar magnetic field exceeds 1.8x10^{14} G. Pairs are created by accelerated particles via resonant scattering of thermal X-rays, followed by the nearly instantaneous conversion of the scattered photon to a pair. A surface gap is then efficiently screened by e+- creation, which regulates the voltage in the inner part of the circuit to ~10^9 V. We also examine the electrostatic gap structure that can form when the magnetic field is somewhat weaker, and deduce a voltage 10-30 times larger over a range of surface temperatures. We examine carefully how the flow of charge back to the star above the gap depends on the magnitude of the current that is extracted from the surface of the star, on the curvature of the magnetic field lines, and on resonant drag. The rates of different channels of pair creation are determined self-consistently, including the non-resonant scattering of X-rays, and collisions between gamma rays and X-rays. We find that the electrostatic gap solution has too small a voltage to sustain the observed pulsed radio output of magnetars unless i) the magnetic axis is nearly aligned with the rotation axis and the light of sight; or ii) the gap is present on the closed as well as the open magnetic field lines. Several properties of the radio magnetars -- their rapid variability, broad pulses, and unusually hard radio spectra -- are consistent with a third possibility, that the current in the outer magnetosphere is strongly variable, and a very high rate of pair creation is sustained by a turbulent cascade. △ Less

Submitted 19 February, 2008; originally announced February 2008.

Comments: 32 pages, submitted to the Astrophysical Journal

Electrodynamics of Magnetars III: Pair Creation Processes in an Ultrastrong Magnetic Field and Particle Heating in a Dynamic Magnetosphere

Authors: Christopher Thompson

Abstract: We consider the details of the QED processes that create electron-positron pairs in magnetic fields approaching and exceeding 10^{14} G. The formation of free and bound pairs is addressed, and the importance of positronium dissociation by thermal X-rays is noted. We calculate the collision cross section between an X-ray and a gamma ray, and point out a resonance in the cross section when the gam… ▽ More We consider the details of the QED processes that create electron-positron pairs in magnetic fields approaching and exceeding 10^{14} G. The formation of free and bound pairs is addressed, and the importance of positronium dissociation by thermal X-rays is noted. We calculate the collision cross section between an X-ray and a gamma ray, and point out a resonance in the cross section when the gamma ray is close to the threshold for pair conversion. We also discuss how the pair creation rate in the open-field circuit and the outer magnetosphere can be strongly enhanced by instabilities near the light cylinder. When the current has a strong fluctuating component, a cascade develops. We examine the details of particle heating, and show that a high rate of pair creation can be sustained close to the star, but only if the spin period is shorter than several seconds. The dissipation rate in this turbulent state can easily accommodate the observed radio output of the transient radio-emitting magnetars, and even their infrared emission. Finally, we outline how a very high rate of pair creation on the open magnetic field lines can help to stabilize a static twist in the closed magnetosphere and to regulate the loss of magnetic helicity by reconnection at the light cylinder. △ Less

Submitted 19 February, 2008; originally announced February 2008.

Comments: 25 pages, submitted to the Astrophysical Journal

Efficient approximations of neutrino physics for three-dimensional simulations of stellar core collapse

Authors: M. Liebendoerfer, U. -L. Pen, C. Thompson

Abstract: Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the fluid-instabilities in the post-bounce phase is indispensable to model supernova explosions. While much effort is required to develop a reliable neutrino transport technique i… ▽ More Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the fluid-instabilities in the post-bounce phase is indispensable to model supernova explosions. While much effort is required to develop a reliable neutrino transport technique in axisymmetry, we explore neutrino physics approximations and parameterizations for an efficient three-dimensional simulation of the fluid-instabilities in the shock-heated matter that accumulates between the accretion shock and the protoneutron star. We demonstrate the reliability of a simple parameterization scheme in the collapse phase and extend our 3D magneto-hydrodynamical collapse simulations to a preliminary postbounce evolution. The growth of magnetic fields is investigated. △ Less

Submitted 24 September, 2006; originally announced September 2006.

Comments: 5 pages, 4 figures, in Proceedings of "Nuclei in the Cosmos IX, Geneva, Jun 25-30", associated movies are displayed at http://www.physik.unibas.ch/~liebend/display

Journal ref: PoSNIC-IX:132,2006

Magnetar Corona

Authors: A. M. Beloborodov, C. Thompson

Abstract: Persistent high-energy emission of magnetars is produced by a plasma corona around the neutron star, with total energy output of ~10^{36}erg/s. The corona forms as a result of occasional starquakes that twist the external magnetic field of the star and induce electric currents in the closed magnetosphere. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. T… ▽ More Persistent high-energy emission of magnetars is produced by a plasma corona around the neutron star, with total energy output of ~10^{36}erg/s. The corona forms as a result of occasional starquakes that twist the external magnetic field of the star and induce electric currents in the closed magnetosphere. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. The voltage along the magnetic field lines is maintained near threshold for ignition of pair production, in the regime of self-organized criticality. The voltage is found to be about 1 GeV which is in agreement with the observed dissipation rate \~10^{36}erg/s. The coronal particles impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona also emits curvature radiation up to 10^{14}Hz and can supply the observed IR-optical luminosity. △ Less

Submitted 17 August, 2006; originally announced August 2006.

Comments: 9 pages, 6 figures; talk at "Isolated Neutron Stars: from the Interior to the Surface," London, April 2006

Journal ref: Astrophys.SpaceSci.308:631-639,2007

Thermalization in Relativistic Outflows and the Correlation between Spectral Hardness and Apparent Luminosity in Gamma-ray Bursts

Authors: C. Thompson, P. Meszaros, M. J. Rees

Abstract: We present an interpretation of the phenomenological relations between the spectral peak, isotropic luminosity and duration of long gamma ray bursts that have been discovered by Amati et al., Ghirlanda et al., Firmani et al., and Liang & Zhang. In our proposed model, a jet undergoes internal dissipation which prevents its bulk Lorentz factor from exceeding 1/theta (theta being the jet opening an… ▽ More We present an interpretation of the phenomenological relations between the spectral peak, isotropic luminosity and duration of long gamma ray bursts that have been discovered by Amati et al., Ghirlanda et al., Firmani et al., and Liang & Zhang. In our proposed model, a jet undergoes internal dissipation which prevents its bulk Lorentz factor from exceeding 1/theta (theta being the jet opening angle) until it escapes from the core of its progenitor star, at a radius of order 1e10 cm; dissipation may continue at larger radii. The dissipated radiation will be partially thermalized, and we identify its thermal peak (Doppler boosted by the outflow) with E_pk. The radiation comes, in effect, from within the jet photosphere. The non-thermal, high energy part of the GRB emission arises from Comptonization of this radiation by relativistic electrons and positrons outside the effective photosphere. This model can account naturally not only for the surprisingly small scatter in the various claimed correlations, but also for the normalization, as well as the slopes. It then has further implications for the jet energy, the limiting jet Lorentz factor, and the relation of the energy, opening angle and burst duration to the mass and radius of the stellar stellar progenitor. The observed relation between pulse width and photon frequency can be reproduced by inverse-Compton emission at ~ 1e14 cm from the engine, but there are significant constraints on the energy distribution and isotropy of the radiating particles. △ Less

Submitted 25 April, 2007; v1 submitted 14 August, 2006; originally announced August 2006.

Comments: 29 pages, Astrophysical Journal, vol. 664, in press

Journal ref: Astrophys.J.666:1012-1023,2007

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

Authors: Rodrigo Fernandez, Christopher Thompson

Abstract: 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 thes… ▽ More 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. △ Less

Submitted 8 April, 2007; v1 submitted 14 August, 2006; originally announced August 2006.

Comments: 29 pages, 23 figures, published in the Astrophysical Journal

Journal ref: Astrophys.J.660:615-640,2007

Corona of Magnetars

Authors: Andrei M. Beloborodov, Christopher Thompson

Abstract: We develop a theoretical model that explains the formation of hot coronae around strongly magnetized neutron stars -- magnetars. The starquakes of a magnetar shear its external magnetic field, which becomes non-potential and is threaded by an electric current. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particl… ▽ More We develop a theoretical model that explains the formation of hot coronae around strongly magnetized neutron stars -- magnetars. The starquakes of a magnetar shear its external magnetic field, which becomes non-potential and is threaded by an electric current. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time of 10^{-4} s and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. Using a direct numerical experiment, we show that the corona self-organizes quickly (on a millisecond timescale) into a quasi-steady state, with voltage ~1 GeV along the magnetic lines. The heating rate of the corona is ~10^{36} erg/s, in agreement with the observed persistent, high-energy output of magnetars. We deduce that a static twist that is suddenly implanted into the magnetosphere will decay on a timescale of 1-10 yrs. The particles accelerated in the corona impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona emits curvature radiation and can supply the observed IR-optical luminosity. (Abridged) △ Less

Submitted 15 August, 2006; v1 submitted 19 February, 2006; originally announced February 2006.

Comments: 70 pages, 14 figures, accepted to ApJ

Journal ref: Astrophys.J.657:967-993,2007

Deceleration of a Relativistic, Photon-Rich Shell: End of Preacceleration, Damping of MHD Turbulence, and the Emission Mechanism of Gamma-Ray Bursts

Authors: Christopher Thompson

Abstract: (Abridged) We consider the interaction of a relativistically-moving shell, composed of thermal photons, a reversing magnetic field and a small admixture of charged particles, with a dense Wolf-Rayet wind. A thin outer layer of Wolf-Rayet material is entrained by the jet head; it cools and becomes Rayleigh-Taylor unstable, thereby providing an additional source of inertia and variability. Pair cr… ▽ More (Abridged) We consider the interaction of a relativistically-moving shell, composed of thermal photons, a reversing magnetic field and a small admixture of charged particles, with a dense Wolf-Rayet wind. A thin outer layer of Wolf-Rayet material is entrained by the jet head; it cools and becomes Rayleigh-Taylor unstable, thereby providing an additional source of inertia and variability. Pair creation in the wind material, and the associated pre-acceleration, defines a characteristic radiative compactness at the point where the reverse shock has completed its passage back through the shell. We argue that the prompt gamma-ray emission is triggered by this external braking, at an optical depth ~1 to electron scattering. Torsional waves, excited by the forced reconnection of the reversing magnetic field, carry a fluctuating current, and are damped at high frequencies by the electrostatic acceleration of electrons and positrons. We show that inverse Compton radiation by the accelerated charges is stronger than their synchrotron emission, and is beamed along the magnetic field. Thermal radiation that is advected out from the base of the jet cools the particles. The observed relation between peak energy and isotropic luminosity is reproduced if the blackbody seeds are generated in a relativistic jet core that is subject to Kelvin-Helmholtz instabilities with the Wolf-Rayet envelope. This relation is predicted to soften to E_peak ~ L_iso^{1/4} below an isotropic luminosity L_iso ~ 3x10^{50} ergs/s. The duration of spikes in the inverse-Compton emission is narrower at higher frequencies, in agreement with the observed relation. The transition from prompt gamma-ray emission to afterglow can be explained by the termination of the thermal X-ray seed and the onset of synchrotron-self-Compton emission. △ Less

Submitted 25 April, 2006; v1 submitted 18 July, 2005; originally announced July 2005.

Comments: 32 pages, 5 figures, in press in the Astrophysical Journal, vol. 647

Journal ref: Astrophys.J.651:333-365,2006

Force-Free Magnetohydrodynamic Waves: Non-Linear Interactions and Effects of Strong Gravity

Authors: Parker Troischt, Christopher Thompson

Abstract: The propagation and non-linear interactions of magnetohydrodynamic waves are considered in the force-free limit, where the inertia of the conducting matter which enforces the MHD condition E.B = 0 can be neglected in comparison with the inertia of the electromagnetic field. By extending the analysis beyond the WKB approximation, we are able to study the non-linearities induced by a gravitational… ▽ More The propagation and non-linear interactions of magnetohydrodynamic waves are considered in the force-free limit, where the inertia of the conducting matter which enforces the MHD condition E.B = 0 can be neglected in comparison with the inertia of the electromagnetic field. By extending the analysis beyond the WKB approximation, we are able to study the non-linearities induced by a gravitational field. We treat the perturbed electromagnetic field as a fluid of infinite conductivity. We calculate the scattering of a torsional (Alfven) wave by a gravitational potential, and demonstrate a nonlinear coupling with a compressive (fast) wave which is second order in the amplitude of the Alfven wave. In a cylindrically symmetric spacetime with slow rotation, the coupling is second order in g_{tφ} and first order in the amplitude of the wave. We also give a fresh analysis of the non-linear interactions between compressive and torsional waves in Minkowski space, with a focus on the relative strengths of their three- and four-mode interactions. In contrast with non-relativistic magnetofluids, the effects of compression are always present. In the case of colliding fast waves, a net displacement of the field lines across (at least) one of the colliding wavepackets is shown to have a strong effect on the outgoing waveform, and to have a qualitatively different interpretation than was previously suggested for colliding Alfven waves. Finally, we show how spacetime curvature modifies the collision between two torsional waves, in both the weak- and strong-field regimes. △ Less

Submitted 16 July, 2005; originally announced July 2005.

Comments: 29 pages

Report number: CITA-2004-19

Journal ref: Phys.Rev.D70:124030,2004

Magnetospheric Eclipses in the Binary Pulsar J0737-3039

Authors: Maxim Lyutikov, Christopher Thompson

Abstract: (Abridged) In the binary radio pulsar system J0737-3039, the faster pulsar A is eclipsed once per orbit. We construct a simple geometric model which successfully reproduces the eclipse light curves, based on the idea that the radio pulses are attenuated by synchrotron absorption on the closed magnetic field lines of pulsar B. The model explains most of the properties of the eclipse: its asymmetr… ▽ More (Abridged) In the binary radio pulsar system J0737-3039, the faster pulsar A is eclipsed once per orbit. We construct a simple geometric model which successfully reproduces the eclipse light curves, based on the idea that the radio pulses are attenuated by synchrotron absorption on the closed magnetic field lines of pulsar B. The model explains most of the properties of the eclipse: its asymmetric form, the nearly frequency-independent duration, and the modulation of the brightness of pulsar A at both once and twice the rotation frequency of pulsar B in different parts of the eclipse. This detailed agreement confirms the dipolar structure of the star's poloidal magnetic field. The model makes clear predictions for the degree of linear polarization of the transmitted radiation. The weak frequency dependence of the eclipse duration implies that the absorbing plasma is relativistic, with a density much larger than the corotation charge density. Such hot, dense plasma can be effectively stored in the outer magnetosphere, where cyclotron cooling is slow. The gradual loss of particles inward through the cooling radius is compensated by an upward flux driven by a fluctuating component of the current, and by the pumping of magnetic helicity on the closed field lines. The trapped particles are heated to relativistic energies by the damping of magnetospheric turbulence and, at a slower rate, by the absorption of the radio emission of the companion pulsar. △ Less

Submitted 17 February, 2005; v1 submitted 17 February, 2005; originally announced February 2005.

Comments: 45 pages, 13 figures, submitted to ApJ; minor typo corrected

Journal ref: Astrophys.J. 634 (2005) 1223-1241

An exceptionally bright flare from SGR1806-20 and the origins of short-duration gamma-ray bursts

Authors: K. Hurley, S. E. Boggs, D. M. Smith, R. C. Duncan, R. Lin, A. Zoglauer, S. Krucker, G. Hurford, H. Hudson, C. Wigger, W. Hajdas, C. Thompson, I. Mitrofanov, A. Sanin, W. Boynton, C. Fellows, A. von Kienlin, G. Lichti, A. Rau, T. Cline

Abstract: Soft-gamma-ray repeaters (SGRs) are galactic X-ray stars that emit numerous short-duration (about 0.1 s) bursts of hard X-rays during sporadic active periods. They are thought to be magnetars: strongly magnetized neutron stars with emissions powered by the dissipation of magnetic energy. Here we report the detection of a long (380 s) giant flare from SGR 1806-20, which was much more luminous tha… ▽ More Soft-gamma-ray repeaters (SGRs) are galactic X-ray stars that emit numerous short-duration (about 0.1 s) bursts of hard X-rays during sporadic active periods. They are thought to be magnetars: strongly magnetized neutron stars with emissions powered by the dissipation of magnetic energy. Here we report the detection of a long (380 s) giant flare from SGR 1806-20, which was much more luminous than any previous transient event observed in our Galaxy. (In the first 0.2 s, the flare released as much energy as the Sun radiates in a quarter of a million years.) Its power can be explained by a catastrophic instability involving global crust failure and magnetic reconnection on a magnetar, with possible large-scale untwisting of magnetic field lines outside the star. From a great distance this event would appear to be a short-duration, hard-spectrum cosmic gamma-ray burst. At least a significant fraction of the mysterious short-duration gamma-ray bursts therefore may come from extragalactic magnetars. △ Less

Submitted 11 May, 2005; v1 submitted 16 February, 2005; originally announced February 2005.

Comments: 21 pages, 5 figures. Published in Nature

Journal ref: Nature 434, 1098, 2005

High Energy Emission from Magnetars

Authors: C. Thompson, A. M. Beloborodov

Abstract: The recently discovered soft gamma-ray emission from the anomalous X-ray pulsar 1E 1841-045 has a luminosity L_g ~ 10^{36} ergs/s. This luminosity exceeds the spindown power by three orders of magnitude and must be fed by an alternative source of energy such as an ultrastrong magnetic field. A gradual release of energy in the stellar magnetosphere is expected if it is twisted and a strong electr… ▽ More The recently discovered soft gamma-ray emission from the anomalous X-ray pulsar 1E 1841-045 has a luminosity L_g ~ 10^{36} ergs/s. This luminosity exceeds the spindown power by three orders of magnitude and must be fed by an alternative source of energy such as an ultrastrong magnetic field. A gradual release of energy in the stellar magnetosphere is expected if it is twisted and a strong electric current is induced on the closed field lines. We examine two mechanisms of gamma-ray emission associated with the gradual dissipation of this current. (1) A thin surface layer of the star is heated by the downward beam of current-carrying charges, which excite Langmuir turbulence in the layer. As a result, it can reach a temperature kT ~ 100 keV and emit bremsstrahlung photons up to this characteristic energy. (2) The magnetosphere is also a source of soft gamma rays at a distance of ~100 km from the star, where the electron cyclotron energy is in the keV range. A large electric field develops in this region in response to the outward drag force felt by the current-carrying electrons from the flux of keV photons leaving the star. A seed positron injected in this region undergoes a runaway acceleration and upscatters keV photons above the threshold for pair creation. The created pairs emit a synchrotron spectrum consistent with the observed 20-100 keV emission. This spectrum is predicted to extend to higher energies and reach a peak at ~1 MeV. △ Less

Submitted 30 November, 2005; v1 submitted 30 August, 2004; originally announced August 2004.

Comments: Astrophysical Journal, 634, 565 (2005)

Journal ref: Astrophys.J. 634 (2005) 565-569

Approaching the dynamics of hot nucleons in supernovae

Authors: M. Liebendoerfer, U. -L. Pen, C. Thompson

Abstract: All recent numerical simulations agree that stars in the main sequence mass range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the outer layers after core collapse and bounce. Rather they suggest that stellar core collapse and supernova explosion are dynamically distinct astrophysical events, separated by an unspectacular accretion phase of at least ~40 ms duration. As l… ▽ More All recent numerical simulations agree that stars in the main sequence mass range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the outer layers after core collapse and bounce. Rather they suggest that stellar core collapse and supernova explosion are dynamically distinct astrophysical events, separated by an unspectacular accretion phase of at least ~40 ms duration. As long as the neutrinospheres remain convectively stable, the explosion dynamics is determined by the neutrons, protons, electrons and neutrinos in the layer of impact-heated matter piling up on the protoneutron star. The crucial role of neutrino transport in this regime has been emphasized in many previous investigations. Here, we search for efficient means to address the role of magnetic fields and fluid instabilities in stellar core collapse and the postbounce phase. △ Less

Submitted 9 August, 2004; originally announced August 2004.

Comments: 4 pages, contribution to Nuclei in the Cosmos VIII, Jul. 19-23, submitted to Nucl. Phys. A

Soft Gamma Repeaters and Anomalous X-ray Pulsars: Magnetar Candidates

Authors: P. M. Woods, C. Thompson

Abstract: This article is a review of Soft Gamma Repeaters and Anomalous X-ray Pulsars. It contains a brief historical record of the emergence of these classes of neutron stars, a thorough overview of the observational data, a succinct summary of the magnetar model, and suggested directions for future research in this field. This article is a review of Soft Gamma Repeaters and Anomalous X-ray Pulsars. It contains a brief historical record of the emergence of these classes of neutron stars, a thorough overview of the observational data, a succinct summary of the magnetar model, and suggested directions for future research in this field. △ Less

Submitted 4 November, 2004; v1 submitted 4 June, 2004; originally announced June 2004.

Comments: 51 pages, 15 figures, and 5 tables. To appear in "Compact Stellar X-ray Sources", eds. W.H.G. Lewin and M. van der Klis

The INTEGRAL view of the Soft Gamma-ray Repeater SGR 1806-20

Authors: D. Gotz, S. Mereghetti, I. F. Mirabel, K. Hurley, S. Brandt, N. Lund, P. Ubertini, M. Del Santo, E. Costa, M. Feroci, A. Castro-Tirado, A. Gimenez, J. -L. Atteia, M. Boer, T. Cline, F. Frontera, G. Pizzichini, A. von Kienlin, E. Gougus, C. Kouveliotou, M. Finger, C. Thompson, H. Pedersen, M. van der Klis

Abstract: We present the results obtained by INTEGRAL on the Soft-Gamma Ray Repeater SGR 1806-20. In particular we report on the temporal and spectral properties of the bursts detected during a moderately active period of the source in September and October 2003 and on the search for quiescent emission. We present the results obtained by INTEGRAL on the Soft-Gamma Ray Repeater SGR 1806-20. In particular we report on the temporal and spectral properties of the bursts detected during a moderately active period of the source in September and October 2003 and on the search for quiescent emission. △ Less

Submitted 7 May, 2004; originally announced May 2004.

Comments: To appear in the proceedings (ESA-SP) of the 5th INTEGRAL Workshop, "The INTEGRAL UNIVERSE", Munich, 16-20 February 2004

Magnetars: Time Evolution, Superfluid Properties, and Mechanism of Magnetic Field Decay

Authors: P. Arras, A. Cumming, C. Thompson

Abstract: We calculate the coupled thermal evolution and magnetic field decay in relativistic model neutron stars threaded by superstrong magnetic fields (B > 10^{15} G). Our main goal is to evaluate how such ``magnetars'' evolve with time and how field decay modifies the transitions to core superfluidity and cooling dominated by surface X-ray emission. Observations of a thermal X-ray spectral component a… ▽ More We calculate the coupled thermal evolution and magnetic field decay in relativistic model neutron stars threaded by superstrong magnetic fields (B > 10^{15} G). Our main goal is to evaluate how such ``magnetars'' evolve with time and how field decay modifies the transitions to core superfluidity and cooling dominated by surface X-ray emission. Observations of a thermal X-ray spectral component and fast timing noise place strong constraints on the presence of a superfluid core. We find that the transition to core superfluidity can be significantly delayed by field decay in the age range ~ 10^3-10^5 yrs. The mechanism of Hall drift is related to the stability of the core magnetic field, and to currents flowing outward through the crust. The heating effect is enhanced if it is continuous rather than spasmodic. Condensation of a heavy element layer at the surface is shown to cause only modest changes in the outward conduction of heat. △ Less

Submitted 26 January, 2004; originally announced January 2004.

Comments: Submitted to Ap.J. Letters 1 December 2003; revised 16 January 2004

Journal ref: Astrophys.J. 608 (2004) L49-L52