Showing 1–50 of 134 results for author: Wood, S

An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models

Authors: Ebraheem Farag, Christopher J. Fontes, F. X. Timmes, Earl P. Bellinger, Joyce A. Guzik, Evan B. Bauer, Suzannah R. Wood, Katie Mussack, Peter Hakel, James Colgan, David P. Kilcrease, Manolo E. Sherrill, Tryston C. Raecke, Morgan T. Chidester

Abstract: We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10… ▽ More We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 28 pages, 12 (13) figures. Accepted for Publication in ApJ

Long-Term Density Trend in the Mesosphere and Lower Thermosphere from Occultations of the Crab Nebula with X-Ray Astronomy Satellites

Authors: Satoru Katsuda, Teruaki Enoto, Andrea N. Lommen, Koji Mori, Yuko Motizuki, Motoki Nakajima, Nathaniel C. Ruhl, Kosuke Sato, Gunter Stober, Makoto S. Tashiro, Yukikatsu Terada, Kent S. Wood

Abstract: We present long-term density trends of the Earth's upper atmosphere at altitudes between 71 and 116 km, based on atmospheric occultations of the Crab Nebula observed with X-ray astronomy satellites, ASCA, RXTE, Suzaku, NuSTAR, and Hitomi. The combination of the five satellites provides a time period of 28 yr from 1994 to 2022. To suppress seasonal and latitudinal variations, we concentrate on the… ▽ More We present long-term density trends of the Earth's upper atmosphere at altitudes between 71 and 116 km, based on atmospheric occultations of the Crab Nebula observed with X-ray astronomy satellites, ASCA, RXTE, Suzaku, NuSTAR, and Hitomi. The combination of the five satellites provides a time period of 28 yr from 1994 to 2022. To suppress seasonal and latitudinal variations, we concentrate on the data taken in autumn (49< doy <111) and spring (235< doy <297) in the northern hemisphere with latitudes of 0--40 degrees. With this constraint, local times are automatically limited either around noon or midnight. We obtain four sets (two seasons times two local times) of density trends at each altitude layer. We take into account variations due to a linear trend and the 11-yr solar cycle using linear regression techniques. Because we do not see significant differences among the four trends, we combine them to provide a single vertical profile of trend slopes. We find a negative density trend of roughly -5 %/decade at every altitude. This is in reasonable agreement with inferences from settling rate of the upper atmosphere. In the 100--110 km altitude, we found an exceptionally high density decline of about -12 %/decade. This peak may be the first observational evidence for strong cooling due to water vapor and ozone near 110 km, which was first identified in a numerical simulation by Akmaev et al. (2006). Further observations and numerical simulations with suitable input parameters are needed to establish this feature. △ Less

Submitted 24 January, 2023; originally announced February 2023.

Comments: 38 pages, 9 figures, accepted for publication in Journal of Geophysical Research - Space Physics

Shear-driven magnetic buoyancy in the solar tachocline: The mean electromotive force due to rotation

Authors: Craig D. Duguid, Paul J. Bushby, Toby S. Wood

Abstract: The leading theoretical paradigm for the Sun's magnetic cycle is an $αω$-dynamo process, in which a combination of differential rotation and turbulent, helical flows produces a large-scale magnetic field that reverses every 11 years. Most $αω$ solar dynamo models rely on differential rotation in the solar tachocline to generate a strong toroidal field. The most problematic part of such models is t… ▽ More The leading theoretical paradigm for the Sun's magnetic cycle is an $αω$-dynamo process, in which a combination of differential rotation and turbulent, helical flows produces a large-scale magnetic field that reverses every 11 years. Most $αω$ solar dynamo models rely on differential rotation in the solar tachocline to generate a strong toroidal field. The most problematic part of such models is then the production of the large-scale poloidal field, via a process known as the $α$-effect. Whilst this is usually attributed to small-scale convective motions under the influence of rotation, the efficiency of this regenerative process has been called into question by some numerical simulations. Motivated by likely conditions within the tachocline, the aim of this paper is to investigate an alternative mechanism for the poloidal field regeneration, namely the magnetic buoyancy instability in a shear-generated, rotating magnetic layer. We use a local, fully compressible model in which an imposed vertical shear winds up an initially vertical magnetic field. The field ultimately becomes buoyantly unstable, and we measure the resulting mean electromotive force (EMF). For sufficiently rapid rotation, we find that a significant component of the mean EMF is aligned with the direction of the mean magnetic field, which is the characteristic feature of the classical $αω$-dynamo model. Our results therefore suggest that magnetic buoyancy could contribute directly to the generation of large-scale poloidal field in the Sun. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS, 16 pages, 12 figures, 3 tables

Validity of sound-proof approximations for magnetic buoyancy

Authors: John B. Moss, Toby S. Wood, Paul J. Bushby

Abstract: The presence of acoustic waves in models of compressible flows can present complications for analytical and numerical analysis. Therefore, several methods have been developed to filter out these waves, leading to various "sound-proof" models, including the Boussinesq, anelastic and pseudo-incompressible models. We assess the validity of each of these approximate models for describing magnetic buoy… ▽ More The presence of acoustic waves in models of compressible flows can present complications for analytical and numerical analysis. Therefore, several methods have been developed to filter out these waves, leading to various "sound-proof" models, including the Boussinesq, anelastic and pseudo-incompressible models. We assess the validity of each of these approximate models for describing magnetic buoyancy in the context of the solar interior. A general sound-proof model is introduced and compared to the fully compressible system in a number of asymptotic regimes, including both non-rotating and rotating cases. We obtain specific constraints that must be satisfied in order that the model captures the leading-order behaviour of the fully compressible system. We then discuss which of the existing sound-proof models satisfy these constraints, and in what parameter regimes. We also present a variational derivation of the pseudo-incompressible MHD model, demonstrating its underlying Hamiltonian structure. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Three-Dimensional Magnetothermal Simulations of Magnetar Outbursts

Authors: Davide De Grandis, Roberto Turolla, Roberto Taverna, Elisa Lucchetta, Toby S. Wood, Silvia Zane

Abstract: The defining trait of magnetars, the most strongly magnetized neutron stars (NSs), is their transient activity in the X/$γ$-bands. In particular, many of them undergo phases of enhanced emission, the so-called outbursts, during which the luminosity rises by a factor $\sim$10$-$1000 in a few hours to then decay over months/years. Outbursts often exhibit a thermal spectrum, associated with the appea… ▽ More The defining trait of magnetars, the most strongly magnetized neutron stars (NSs), is their transient activity in the X/$γ$-bands. In particular, many of them undergo phases of enhanced emission, the so-called outbursts, during which the luminosity rises by a factor $\sim$10$-$1000 in a few hours to then decay over months/years. Outbursts often exhibit a thermal spectrum, associated with the appearance of hotter regions on the surface of the star, which subsequently change in shape and cool down. Here we simulate the unfolding of a sudden, localized heat injection in the external crust of a NS with a 3D magneto-thermal evolution code, finding that this can reproduce the main features of magnetar outbursts. A full 3D treatment allows us to study for the first time the inherently asymmetric hot-spots which appear on the surface of the star as the result of the injection and to follow the evolution of their temperature and shape. We investigate the effects produced by different physical conditions in the heated region, highlighting in particular how the geometry of the magnetic field plays a key role in determining the properties of the event. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 12 pages, 10 figures, accepted for publication in ApJ

Fermi-LAT Gamma-ray Detection of the Recurrent Nova RS Ophiuchi during its 2021 Outburst

Authors: C. C. Cheung, T. J. Johnson, P. Jean, M. Kerr, K. L. Page, J. P. Osborne, A. P. Beardmore, K. V. Sokolovsky, F. Teyssier, S. Ciprini, G. Marti-Devesa, I. Mereu, S. Razzaque, K. S. Wood, S. N. Shore, S. Korotkiy, A. Levina, A. Blumenzweig

Abstract: We report the Fermi-LAT gamma-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ~20-23 GeV, the highest-energy photons detected in… ▽ More We report the Fermi-LAT gamma-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ~20-23 GeV, the highest-energy photons detected in some sub-intervals, are consistent with $π^{\rm 0}$-decay emission from shocks in the ejecta as proposed by Tatischeff & Hernanz (2007) for its previous 2006 outburst. The LAT light-curve displayed a fast rise to its peak >0.1 GeV flux of $\simeq$6x10^-6 ph cm^-2 s^-1 beginning on day 0.745 after its optically-constrained eruption epoch of 2021 August 8.50. The peak lasted for ~1 day, and exhibited a power-law decline up to the final LAT detection on day 45. We analyze the data on shorter timescales at early times and found evidence of an approximate doubling of emission over ~200 minutes at day 2.2, possibly indicating a localized shock-acceleration event. Comparing the data collected by the AAVSO, we measured a constant ratio of ~2.8x10^-3 between the gamma-ray and optical luminosities except for a ~5x smaller ratio within the first day of the eruption likely indicating attenuation of gamma rays by ejecta material and lower high-energy proton fluxes at the earliest stages of the shock development. The hard X-ray emission due to bremsstrahlung from shock-heated gas traced by the Swift-XRT 2-10 keV light-curve peaked at day ~6, later than at GeV and optical energies. Using X-ray derived temperatures to constrain the velocity profile, we find the hadronic model reproduces the observed >0.1 GeV light-curve. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: ApJ, accepted. 21 pages, 10 figures, 4 tables

A Detection of Red Noise in PSR J1824$-$2452A and Projections for PSR B1937+21 using NICER X-ray Timing Data

Authors: Jeffrey S. Hazboun, Jack Crump, Andrea N. Lommen, Sergio Montano, Samantha J. H. Berry, Jesse Zeldes, Elizabeth Teng, Paul S. Ray, Matthew Kerr, Zaven Arzoumanian, Slavko Bogdanov, Julia Deneva, Natalia Lewandowska, Craig B. Markwardt, Scott Ransom, Teruaki Enoto, Kent S. Wood, Keith C. Gendreau, David A. Howe, Aditya Parthasarathy

Abstract: We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar mediu… ▽ More We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long time-scale, correlated variations ("red noise") in the pulse times of arrival from the millisecond pulsars PSR J1824$-$2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency dependent propagation effects of the interstellar medium. Our Bayesian search methodology yields strong evidence (natural log Bayes factor of $9.634 \pm 0.016$) for red noise in PSR J1824$-$2452A, but is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher precision X-ray datasets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5-year mission with a time-of-arrival (TOA) error of 2 microseconds and an observing cadence of 20 observations per month compared to the 5 microsecond TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that an injected stochastic gravitational wave background (GWB) with an amplitude of $A_{\rm GWB}=2\times10^{-15}$ and spectral index of $γ_{\rm GWB}=13/3$ can be detected in a pulsar with similar TOA precision to PSR B1937+21, but with no additional red noise, in a 10-year mission that observes the pulsar 15 times per month and has an average TOA error of 1 microsecond. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 21 pages, 13 figures

NICER Detection of Thermal X-ray Pulsations from the Massive Millisecond Pulsars PSR J0740+6620 and PSR J1614-2230

Authors: M. T. Wolff, S. Guillot, S. Bogdanov, P. S. Ray, M. Kerr, Z. Arzoumanian, K. C. Gendreau, M. C. Miller, A. J. Dittmann, W. C. G. Ho, L. Guillemot, I. Cognard, G. Theureau, K. S. Wood

Abstract: We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X… ▽ More We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X-ray fields. PSR J0740+6620 exhibits a broad double-peaked profile with a separation of ~0.4 in phase. PSR J1614-2230, on the other hand, has a broad single-peak profile. The broad modulations with soft X-ray spectra of both pulsars are indicative of thermal radiation from one or more small regions of the stellar surface. We show the NICER detections of X-ray pulsations for both pulsars and also discuss the phase relationship to their radio pulsations. In the case of PSR J0740+6620, this paper documents the data reduction performed to obtain the pulsation detection and prepare for pulse profile modeling analysis. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: Submitted to The Astrophysical Journal Letters

X-ray Emission from Isolated Neutron Stars revisited: 3D magnetothermal simulations

Authors: Davide De Grandis, Roberto Taverna, Roberto Turolla, Andrea Gnarini, Sergei B. Popov, Silvia Zane, Toby S. Wood

Abstract: X-ray emission from the surface of isolated neutron stars (NSs) has been now observed in a variety of sources. The ubiquitous presence of pulsations clearly indicates that thermal photons either come from a limited area, possibly heated by some external mechanism, or from the entire (cooling) surface but with an inhomogeneous temperature distribution. In a NS the thermal map is shaped by the magne… ▽ More X-ray emission from the surface of isolated neutron stars (NSs) has been now observed in a variety of sources. The ubiquitous presence of pulsations clearly indicates that thermal photons either come from a limited area, possibly heated by some external mechanism, or from the entire (cooling) surface but with an inhomogeneous temperature distribution. In a NS the thermal map is shaped by the magnetic field topology, since heat flows in the crust mostly along the magnetic field lines. Self-consistent surface thermal maps can hence be produced by simulating the coupled magnetic and thermal evolution of the star. We compute the evolution of the neutron star crust in three dimensions for different initial configurations of the magnetic field and use the ensuing thermal surface maps to derive the spectrum and the pulse profile as seen by an observer at infinity, accounting for general-relativistic effects. In particular, we compare cases with a high degree of symmetry with inherently 3D ones, obtained by adding a quadrupole to the initial dipolar field. Axially symmetric fields result in rather small pulsed fractions ($\lesssim 5\%$), while more complex configurations produce higher pulsed fractions, up to $\sim25\%$. We find that the spectral properties of our axisymmetric model are close to those of the bright isolated NS RX~J1856.5-3754 at an evolutionary time comparable with the inferred dynamical age of the source. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 14 pages, 11 figures, accepted for publication in ApJ

A long term study of Mars mesospheric clouds seen at twilight based on Mars Express VMC images

Authors: J. Hernandez-Bernal, A. Sanchez-Lavega, T. del Rio-Gaztelurrutia, R. Hueso, E. Ravanis, A. Cardesin-Moinelo, S. Wood, D. Titov

Abstract: We present the first systematic study of clouds observed during twilight on Mars. We analyze images obtained by the Visual Monitoring Camera (VMC) on Mars Express between 2007 and 2020. Using an automated retrieval algorithm we found 407 cases of clouds observed at twilight, in which the geometry of the observations allows to derive the minimum altitude, revealing that many of these clouds are in… ▽ More We present the first systematic study of clouds observed during twilight on Mars. We analyze images obtained by the Visual Monitoring Camera (VMC) on Mars Express between 2007 and 2020. Using an automated retrieval algorithm we found 407 cases of clouds observed at twilight, in which the geometry of the observations allows to derive the minimum altitude, revealing that many of these clouds are in the mesosphere (above 40km and up to 90km). The majority of these mesospheric clouds were detected in mid-latitudes at local autumn and winter, a new trend only hinted at by previous studies. In particular, we find a massive concentration of clouds in the southern mid-latitudes between Terra Cimmeria and Aonia, a region where high altitude events have been previously observed. We propose that there is an unknown mechanism in these regions that enhances the probability to host high altitude clouds around the southern winter solstice. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Journal ref: Volume48, Issue7, April 2021

An Extremely Elongated Cloud over Arsia Mons Volcano on Mars: I. Life Cycle

Authors: J. Hernández-Bernal, A. Sánchez-Lavega, T. del Río-Gaztelurrutia, E. Ravanis, A. Cardesín-Moinelo, K. Connour, D. Tirsch, I. Ordóñez-Etxeberria, B. Gondet, S. Wood, D. Titov, N. M. Schneider, R. Hueso, R. Jaumann, E. Hauber

Abstract: We report a previously unnoticed annually repeating phenomenon consisting of the daily formation of an extremely elongated cloud extending as far as 1800 km westward from Arsia Mons. It takes place in the Solar Longitude (Ls) range of ~220-320, around the Southern solstice. We study this Arsia Mons Elongated Cloud (AMEC) using images from different orbiters, including ESA Mars Express, NASA MAVEN,… ▽ More We report a previously unnoticed annually repeating phenomenon consisting of the daily formation of an extremely elongated cloud extending as far as 1800 km westward from Arsia Mons. It takes place in the Solar Longitude (Ls) range of ~220-320, around the Southern solstice. We study this Arsia Mons Elongated Cloud (AMEC) using images from different orbiters, including ESA Mars Express, NASA MAVEN, Viking 2, MRO, and ISRO Mars Orbiter Mission (MOM). We study the AMEC in detail in Martian Year (MY) 34 in terms of Local Time and Ls and find that it exhibits a very rapid daily cycle: the cloud growth starts before sunrise on the western slope of the volcano, followed by a westward expansion that lasts 2.5 hours with a velocity of around 170 m/s in the mesosphere (~45 km over the areoid). The cloud formation then ceases, it detaches from its formation point, and continues moving westward until it evaporates before the afternoon, when most sun-synchronous orbiters observe. Moreover we comparatively study observations from different years (i.e. MYs 29-34) in search of interannual variations and find that in MY33 the cloud exhibits lower activity, whilst in MY34 the beginning of its formation was delayed compared to other years, most likely due to the Global Dust Storm. This phenomenon takes place in a season known for the general lack of clouds on Mars. In this paper we focus on observations, and a theoretical interpretation will be the subject of a separate paper. △ Less

Submitted 5 March, 2021; originally announced March 2021.

The 2018 Martian Global Dust Storm over the South Polar Region studied with MEx/VMC

Authors: J. Hernández-Bernal, A. Sánchez-Lavega, T. del Río-Gaztelurrutia, R. Hueso, A. Cardesín-Moinelo, E. Ravanis, A. de Burgos-Sierra, D. Titov, S. Wood

Abstract: We study the 2018 Martian Global DustStorm (GDS 2018) over the Southern Polar Region using images obtained by the Visual Monitoring Camera (VMC) on board Mars Express during June and July 2018. Dust penetrated into the polar cap region but never covered the cap completely, and its spatial distribution was nonhomogeneous and rapidly changing. However, we detected long but narrow aerosol curved arcs… ▽ More We study the 2018 Martian Global DustStorm (GDS 2018) over the Southern Polar Region using images obtained by the Visual Monitoring Camera (VMC) on board Mars Express during June and July 2018. Dust penetrated into the polar cap region but never covered the cap completely, and its spatial distribution was nonhomogeneous and rapidly changing. However, we detected long but narrow aerosol curved arcs with a length of 2,000-3,000 km traversing part of the cap and crossing the terminator into the night side. Tracking discrete dust clouds allowed measurements of their motions that were towards the terminator with velocities up to 100 m/s. The images of the dust projected into the Martian limb show maximum altitudes of around 70 km but with large spatial and temporal variations. We discuss these results in the context of the predictions of a numerical model for dust storm scenario. △ Less

Submitted 26 February, 2021; originally announced February 2021.

Journal ref: Geophysical Research Letters, 2019, 46(17-18), 10330-10337

3D Magneto-thermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects

Authors: Andrei P. Igoshev, Konstantinos N. Gourgouliatos, Rainer Hollerbach, Toby S. Wood

Abstract: Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this ca… ▽ More Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magneto-thermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns which consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of $2.5-10\times 10^{47}$ erg have temperatures of hot regions that reach $\approx 0.2$ keV, to be compared with the bulk temperature of $\approx 0.1$ keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of central compact objects. The hot spots produce periodic modulations in light curve with typical amplitudes of $\leq 9-11$ %. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some central compact objects. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 12 pages, 25 figures, accepted for publication ApJ on January 19, 2021

Superconducting phases in a two-component microscale model of neutron star cores

Authors: Toby S. Wood, Vanessa Graber, William G. Newton

Abstract: We identify the possible ground states for a mixture of two superfluid condensates (one neutral, the other electrically charged) using a phenomenological Ginzburg-Landau model. While this framework is applicable to any interacting condensed-matter mixture of a charged and a neutral component, we focus on nuclear matter in neutron star cores, where proton and neutron condensates are coupled via non… ▽ More We identify the possible ground states for a mixture of two superfluid condensates (one neutral, the other electrically charged) using a phenomenological Ginzburg-Landau model. While this framework is applicable to any interacting condensed-matter mixture of a charged and a neutral component, we focus on nuclear matter in neutron star cores, where proton and neutron condensates are coupled via non-dissipative entrainment. We employ the Skyrme interaction to determine the neutron star's equilibrium composition, and hence obtain realistic coefficients for our Ginzburg-Landau model at each depth within the star's core. We then use the Ginzburg-Landau model to determine the ground state in the presence of a magnetic field. In this way, we obtain superconducting phase diagrams for six representative Skyrme models, revealing the microphysical magnetic flux distribution throughout the neutron star core. The phase diagrams are rather complex and the locations of most of the phase transitions can only be determined through numerical calculations. Nonetheless, we find that for all equations of state considered in this work, much of the outer core exhibits type-1.5 superconductivity, rather than type-II superconductivity as is generally assumed. For local magnetic field strengths $\lesssim 10^{14} \, {\rm G}$, the magnetic flux is distributed inhomogeneously, with bundles of magnetic fluxtubes separated by flux-free Meissner regions. We provide an approximate criterion to determine the transition between this type-1.5 phase and the type-I region in the inner core. △ Less

Submitted 26 February, 2021; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 33 pages, 15 figures

3D modelling of magneto-thermal evolution of neutron stars: method and test cases

Authors: Davide De Grandis, Roberto Turolla, Toby S. Wood, Silvia Zane, Roberto Taverna, Konstantinos N. Gourgouliatos

Abstract: Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the star's observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emi… ▽ More Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the star's observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emission, obtained with the pseudo-spectral code Parody. We investigate both the secular evolution, especially in connection with the onset of instabilities during the Hall phase, and the short-term evolution following episodes of localised energy injection. Simulations show that a resistive tearing instability develops in about a Hall time if the initial toroidal field exceeds ~$10^{15}$ G. This leads to crustal failures because of the huge magnetic stresses coupled with the local temperature enhancement produced by dissipation. Localised heat deposition in the crust results in the appearance of hot spots on the star surface which can exhibit a variety of patterns. Since the transport properties are strongly influenced by the magnetic field, the hot regions tend to drift away and get deformed following the magnetic field lines while cooling. The shapes obtained with our simulations are reminiscent of those recently derived from NICER X-ray observations of the millisecond pulsar PSR J0030+0451. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 15 pages, 14 figures, accepted for publication in ApJ

The Role of Inhomogeneities in Supernova Shock Breakout Emission

Authors: Chris L. Fryer, Christopher J. Fontes, James S. Warsa, Pete W. A. Roming, Shane X. Coffing, Suzannah R. Wood

Abstract: The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in th… ▽ More The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in the immediate surroundings of a massive Wolf-Rayet star, calculating the effect these interactions have on the shock breakout signal from supernovae. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 14 pages, 12 figures, submitted to ApJ

On the saturation mechanism of the fluctuation dynamo at ${\text{Pr}_\mathrm{M}} \ge 1$

Authors: Amit Seta, Paul J. Bushby, Anvar Shukurov, Toby S. Wood

Abstract: The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynam… ▽ More The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynamo in driven turbulence. We demonstrate that as the magnetic field growth saturates, its amplification and diffusion are both affected by the back-reaction of the Lorentz force upon the flow. The amplification of the magnetic field is reduced due to stronger alignment between the velocity field, magnetic field, and electric current density. Furthermore, we confirm that the amplification decreases due to a weaker stretching of the magnetic field lines. The enhancement in diffusion relative to the field line stretching is quantified by a decrease in the computed local value of the magnetic Reynolds number. Using the Minkowski functionals, we quantify the shape of the magnetic structures produced by the dynamo as magnetic filaments and ribbons in both kinematic and saturated dynamos and derive the scalings of the typical length, width, and thickness of the magnetic structures with the magnetic Reynolds number. We show that all three of these magnetic length scales increase as the dynamo saturates. The magnetic intermittency, strong in the kinematic dynamo (where the magnetic field strength grows exponentially) persists in the statistically steady state, but intense magnetic filaments and ribbons are more volume-filling. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: 17 pages, 19 figures, accepted for publication in Physical Review Fluids

A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars

Authors: Dominick M. Rowan, Zaynab Ghazi, Lauren Lugo, Elizabeth Spano, Andrea Lommen, Alice Harding, Christo Venter, Renee Ludlam, Paul S. Ray, Matthew Kerr, Zaven Arzoumanian, Slavko Bogdanov, Julia Deneva, Sebastien Guillot, Natalia Lewandowska, Craig B. Markwardt, Scott Ransom, Teruaki Enoto, Kent S. Wood, Keith Gendreau

Abstract: We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separatio… ▽ More We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3σ$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

NICER X-ray Observations of Seven Nearby Rotation-Powered Millisecond Pulsars

Authors: S. Guillot, M. Kerr, P. S. Ray, S. Bogdanov, S. Ransom, J. S. Deneva, Z. Arzoumanian, P. Bult, D. Chakrabarty, K. C. Gendreau, W. C. G. Ho, G. K. Jaisawal, C. Malacaria, M. C. Miller, T. E. Strohmayer, M. T. Wolff, K. S. Wood, N. A. Webb, L. Guillemot, I. Cognard, G. Theureau

Abstract: NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7σ$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sam… ▽ More NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7σ$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$σ$ and almost 3$σ$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 25 pages, 11 tables, 4 figures. In press in The Astrophysical Journal Letters

Journal ref: ApJL, 887, L27 (2019)

Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-ray Data Set

Authors: Slavko Bogdanov, Sebastien Guillot, Paul S. Ray, Michael T. Wolff, Deepto Chakrabarty, Wynn C. G. Ho, Matthew Kerr, Frederick K. Lamb, Andrea Lommen, Renee M. Ludlam, Reilly Milburn, Sergio Montano, M. Coleman Miller, Michi Baubock, Feryal Ozel, Dimitrios Psaltis, Ronald A. Remillard, Thomas E. Riley, James F. Steiner, Tod E. Strohmayer, Anna L. Watts, Kent S. Wood, Jesse Zeldes, Teruaki Enoto, Takashi Okajima , et al. (5 additional authors not shown)

Abstract: We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the i… ▽ More We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 28 pages, 9 figures; accepted for publication in the Astrophysical Journal Letters

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Authors: M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, A. K. Harding, W. C. G. Ho, J. M. Lattimer, R. M. Ludlam, S. Mahmoodifar, S. M. Morsink, P. S. Ray, T. E. Strohmayer, K. S. Wood, T. Enoto, R. Foster, T. Okajima, G. Prigozhin, Y. Soong

Abstract: Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mas… ▽ More Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 49 pages, 16 figures, part of The Astrophysical Journal Letters focus issue on the Neutron Star Interior Composition Explorer

Journal ref: ApJL 2019, 887, L24

X-ray and Radio Variabilities of PSR J2032+4127 near Periastron

Authors: C. -Y. Ng, W. C. G. Ho, E. V. Gotthelf, J. P. Halpern, M. J. Coe, B. W. Stappers, A. G. Lyne, K. S. Wood, M. Kerr

Abstract: We present X-ray and radio monitoring observations of the gamma-ray binary PSR J2032+4127/MT91 213 during its periastron passage in late 2017. Dedicated Chandra, XMM-Newton,NuSTAR X-ray observations and VLA radio observations of this long orbit (50 years), 143 ms pulsar/Be star system clearly revealed flux and spectral variability during the passage. The X-ray spectrum hardened near periastron, wi… ▽ More We present X-ray and radio monitoring observations of the gamma-ray binary PSR J2032+4127/MT91 213 during its periastron passage in late 2017. Dedicated Chandra, XMM-Newton,NuSTAR X-ray observations and VLA radio observations of this long orbit (50 years), 143 ms pulsar/Be star system clearly revealed flux and spectral variability during the passage. The X-ray spectrum hardened near periastron, with a significant decrease in the power-law photon index from Γ~ 2 to 1.2 and evidence of an increased absorption column density. We identified a possible spectral break at a few keV in the spectrum that suggests synchrotron cooling. A coincident radio and X-ray flare occurred one week after periastron, which is possibly the result of the pulsar wind interacting with the Be stellar disk and generating synchrotron radiation. However, a multi-wavelength comparison indicate that the X-ray and radio spectra cannot be simply connected by a single power-law component. Hence, the emission in these two energy bands must originate from different particle populations. △ Less

Submitted 12 July, 2019; originally announced July 2019.

Comments: Accepted for publication in the ApJ

High-Precision X-ray Timing of Three Millisecond Pulsars with NICER: Stability Estimates and Comparison with Radio

Authors: J. S. Deneva, P. S. Ray, A. Lommen, S. M. Ransom, S. Bogdanov, M. Kerr, K. S. Wood, Z. Arzoumanian, K. Black, J. Doty, K. C. Gendreau, S. Guillot, A. Harding, N. Lewandowska, C. Malacaria, C. B. Markwardt, S. Price, L. Winternitz, M. T. Wolff, L. Guillemot, I. Cognard, P. T. Baker, H. Blumer, P. R. Brook, H. T. Cromartie , et al. (23 additional authors not shown)

Abstract: The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pul… ▽ More The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $σ_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $σ_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $σ_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 30 pages, 15 figures; accepted to ApJ

On energy equipartition between cosmic rays and magnetic fields

Authors: Amit Seta, Anvar Shukurov, Paul J. Bushby, Toby S. Wood

Abstract: Interpretations of synchrotron observations often assume a tight correlation between magnetic and cosmic ray energy densities. We examine this assumption using both test-particle simulations of cosmic rays and MHD simulations which include cosmic rays as a diffusive fluid. We find no spatial correlation between the cosmic rays and magnetic field energy densities at turbulent scales. Moreover, the… ▽ More Interpretations of synchrotron observations often assume a tight correlation between magnetic and cosmic ray energy densities. We examine this assumption using both test-particle simulations of cosmic rays and MHD simulations which include cosmic rays as a diffusive fluid. We find no spatial correlation between the cosmic rays and magnetic field energy densities at turbulent scales. Moreover, the cosmic ray number density and magnetic field energy density are statistically independent. Nevertheless, the cosmic ray spatial distribution is highly inhomogeneous, especially at low energies because the particles are trapped between random magnetic mirrors. These results can significantly change the interpretation of synchrotron observations and thus our understanding of the strength and structure of magnetic fields in the Milky Way and nearby spiral galaxies. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 4 pages, 2 figures, Proceedings of Focus Meeting FM8: New Insights in Extragalactic Magnetic Fields, XXXth IAU General Assembly, Vienna, August 29-31, 2018

Solar Models with dynamic screening and early mass loss tested by helioseismic, astrophysical, and planetary constraints

Authors: Suzannah R. Wood, Katie Mussack, Joyce A. Guzik

Abstract: The faint young Sun paradox remains an open question. Presented here is one possible solution to this apparent inconsistency, a more massive early Sun. Based on conditions on early Earth and Mars, a luminosity constraint is set as a function of mass. We examine helioseismic constraints of these alternative mass-losing models. Additionally, we explore a dynamic electron screening correction in an e… ▽ More The faint young Sun paradox remains an open question. Presented here is one possible solution to this apparent inconsistency, a more massive early Sun. Based on conditions on early Earth and Mars, a luminosity constraint is set as a function of mass. We examine helioseismic constraints of these alternative mass-losing models. Additionally, we explore a dynamic electron screening correction in an effort to improve helioseismic agreement in the core of the an early mass-losing model. △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: 16 pages, 6 figures, accepted Solar Physics

STROBE-X: A probe-class mission for X-ray spectroscopy and timing on timescales from microseconds to years

Authors: Paul S. Ray, Zaven Arzoumanian, Søren Brandt, Eric Burns, Deepto Chakrabarty, Marco Feroci, Keith C. Gendreau, Olivier Gevin, Margarita Hernanz, Peter Jenke, Steven Kenyon, José Luis Gálvez Thomas J. Maccarone, Takashi Okajima, Ronald A. Remillard, Stéphane Schanne, Chris Tenzer, Andrea Vacchi, Colleen A. Wilson-Hodge, Berend Winter, Silvia Zane, David R. Ballantyne, Enrico Bozzo, Laura W. Brenneman, Edward Cackett, Alessandra De Rosa , et al. (8 additional authors not shown)

Abstract: We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The… ▽ More We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration. △ Less

Submitted 3 July, 2018; originally announced July 2018.

Comments: 20 pages, 10 figures, to appear in Proc. SPIE Astronomical Telescopes and Instrumentation 2018

NICER and Fermi GBM Observations of the First Galactic Ultraluminous X-ray Pulsar Swift J0243.6+6124

Authors: C. A. Wilson-Hodge, C. Malacaria, P. A. Jenke, G. K. Jaisawal, M. Kerr, M. T. Wolff, Z. Arzoumanian, D. Chakrabarty, J. P. Doty, K. C. Gendreau, S. Guillot, W. C. G. Ho, B. LaMarr, C. B. Markwardt, F. Ozel, G. Y. Prigozhin, P. S. Ray, M. Ramos-Lerate, R. A. Remillard, T. E. Strohmayer, M. L. Vezie, K. S. Wood

Abstract: Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Bu… ▽ More Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: 35 pages, 7 figures (main text), 5 figures (appendix); Accepted for publication in ApJ on June 19, 2018

A Luminous and Highly-variable Gamma-ray Flare Following the 2017 Periastron of PSR B1259-63/LS 2883

Authors: T. J. Johnson, K. S. Wood, M. Kerr, R. H. D. Corbet, C. C. Cheung, P. S. Ray, N. Omodei

Abstract: Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $\sim30$ M$_{\odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $\sim3$-week long int… ▽ More Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $\sim30$ M$_{\odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $\sim3$-week long interval immediately after periastron, which was followed by an interval of no emission. Sporadic flares were detected starting 40 days post-periastron and lasted approximately 50 days, during which the emission displayed significant spectral curvature, variability on timescales as short as 1.5 minutes, and peak flux levels well in excess of the pulsar spin-down power. By contrast, during the 2010 and 2014 periastron passages, significant gamma-ray emission was not observed with the LAT until 30 and 32 days post-periastron, respectively. The previous flares did not exhibit spectral curvature, showed no short term variability, and did not exceed the pulsar spin-down power. The high flux and short timescales observed in 2017 suggest significant beaming of the emission is required and constrain the size of the emission region. The flares occur long enough after periastron that the neutron star should already have passed through the extended disk-like outflow, thus constraining options for target material and seed photon sources for inverse Compton models. △ Less

Submitted 12 July, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: Accepted for publication in the Astrophysical Journal (2 July 2018), v2 matches accepted version

Journal ref: Astrophysical Journal, 2018, vol. 863, p. 27

A self-consistent model of the solar tachocline

Authors: Toby S. Wood, Nicholas H. Brummell

Abstract: We present a local but fully nonlinear model of the solar tachocline, using three-dimensional direct numerical simulations. The tachocline forms naturally as a statistically steady balance between Coriolis, pressure, buoyancy and Lorentz forces beneath a turbulent convection zone. Uniform rotation is maintained in the radiation zone by a primordial magnetic field, which is confined by meridional f… ▽ More We present a local but fully nonlinear model of the solar tachocline, using three-dimensional direct numerical simulations. The tachocline forms naturally as a statistically steady balance between Coriolis, pressure, buoyancy and Lorentz forces beneath a turbulent convection zone. Uniform rotation is maintained in the radiation zone by a primordial magnetic field, which is confined by meridional flows in the tachocline and convection zone. Such balanced dynamics has previously been found in idealised laminar models, but never in fully self-consistent numerical simulations. △ Less

Submitted 8 January, 2018; originally announced January 2018.

Comments: Accepted for publication in the Astrophysical Journal

The NRL Program in X-ray Navigation

Authors: Kent S. Wood, Paul S. Ray

Abstract: This chapter describes the development of X-ray Navigation at the Naval Research Laboratory (NRL) within its astrophysics research programs. The prospects for applications emerged from early discoveries of X-ray source classes and their properties. Starting around 1988 some NRL X-ray astronomy programs included navigation as one of the motivations. The USA experiment (1999) was the first flight pa… ▽ More This chapter describes the development of X-ray Navigation at the Naval Research Laboratory (NRL) within its astrophysics research programs. The prospects for applications emerged from early discoveries of X-ray source classes and their properties. Starting around 1988 some NRL X-ray astronomy programs included navigation as one of the motivations. The USA experiment (1999) was the first flight payload with an explicit X-ray navigation theme. Subsequently, NRL has continued to work in this area through participation in DARPA and NASA programs. Throughout, the general concept of X-ray navigation (XRNAV) has been broad enough to encompass many different uses of X-ray source observations for attitude determination, position determination, and timekeeping. Pulsar-based X-ray navigation (XNAV) is a special case. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 30 pages, 15 figures, to appear in Proceedings of the 593. WE-Heraeus Seminar on Autonomous Spacecraft Navigation, ed. W. Becker

Characterization of Pulsar Sources for X-ray Navigation

Authors: Paul S. Ray, Kent S. Wood, Michael T. Wolff

Abstract: The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron Star Interior Composition Explorer (NICER) mission, which is scheduled to launch in 2017 and will be hosted as an externally attached payload on the International Space Station (ISS). During NICER's 18-month baseline science mission to understand ultra-dense matter thr… ▽ More The Station Explorer for X-ray Timing and Navigation Technology (SEXTANT) is a technology demonstration enhancement to the Neutron Star Interior Composition Explorer (NICER) mission, which is scheduled to launch in 2017 and will be hosted as an externally attached payload on the International Space Station (ISS). During NICER's 18-month baseline science mission to understand ultra-dense matter through observations of neutron stars in the soft X-ray band, SEXTANT will, for the first-time, demonstrate real-time, on-board X-ray pulsar navigation. Using NICER/SEXTANT as an example, we describe the factors that determine the measurement errors on pulse times of arrival, including source and background count rates, and pulse profile shapes. We then describe properties of the SEXTANT navigation pulsar catalog and prospects for growing it once NICER launches. Finally, we describe the factors affecting the prediction of pulse arrival times in advance, including variable interstellar propagation effect and red timing noise. Together, all of these factors determine how well a particular realization of an X-ray pulsar-based navigation system will perform. △ Less

Submitted 22 November, 2017; originally announced November 2017.

Comments: 16 pages, 7 figures, to appear in Proceedings of the 593. WE-Heraeus Seminar on Autonomous Spacecraft Navigation, ed. W. Becker

Relative distribution of cosmic rays and magnetic fields

Authors: Amit Seta, Anvar Shukurov, Toby S. Wood, Paul J. Bushby, Andrew P. Snodin

Abstract: Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight… ▽ More Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight correlation to hold across all scales. We use test particle simulations, tracing the propagation of charged particles (protons) through a random magnetic field, to study the cosmic ray distribution at scales comparable to the correlation scale of the turbulent flow in the interstellar medium ($\simeq 100\,{\rm pc}$ in spiral galaxies). In these simulations, we find that there is no spatial correlation between the cosmic ray number density and the magnetic field energy density. In fact, their distributions are approximately statistically independent. We find that low-energy cosmic rays can become trapped between magnetic mirrors, whose location depends more on the structure of the field lines than on the field strength. △ Less

Submitted 10 October, 2017; v1 submitted 24 August, 2017; originally announced August 2017.

Comments: 14 pages, 15 figures, 1 table, Accepted for publication in MNRAS, Video showing trajectory of a cosmic ray particle in a magnetic trap- https://youtu.be/EzabLlio3FE

Compton Edge probing basic physics at Jefferson Laboratory: light speed isotropy and Lorentz invariance

Authors: Vahe Gurzadyan, David Gaskell, Vanik Kakoyan, Cynthia Keppel, Amur Margaryan, Harutyun Khachatryan, Sergey Mirzoyan, Dipangkar Dutta, Branislav Vlahovic, Steve Wood

Abstract: We propose to study of the light speed isotropy and Lorentz invariance at Jefferson Laboratory by means of the measurements of the Compton Edge using of the Hall A/C existing experimental setup. Methodologically the same experiment has already been successfully elaborated at GRAAL experiment at the European Synchrotron Radiation Facility in Grenoble with 6 GeV electron beam. This Proposal states t… ▽ More We propose to study of the light speed isotropy and Lorentz invariance at Jefferson Laboratory by means of the measurements of the Compton Edge using of the Hall A/C existing experimental setup. Methodologically the same experiment has already been successfully elaborated at GRAAL experiment at the European Synchrotron Radiation Facility in Grenoble with 6 GeV electron beam. This Proposal states two goals expected to be reached at Jefferson Laboratory, both on Lorentz invariance: (a) the one-way light speed isotropy testing accuracy, following from conservative evaluations at numerical simulations, to about an order of magnitude better than was GRAAL's; (b) the dependence of the light speed on the velocity of the apparatus (Kennedy-Thorndike measurement) will be traced to an accuracy about 3 orders of magnitudes better than the available limits. △ Less

Submitted 27 June, 2017; originally announced June 2017.

Comments: Submitted to PAC45, Jefferson Laboratory, 38 pages, 14 figures

ALMA observations of the interaction of a radio jet with molecular gas in Minkowski's Object

Authors: Mark Lacy, Steve Croft, Chris Fragile, Sarah Wood, Kristina Nyland

Abstract: We use ALMA to detect and image CO (1-0) emission from Minkowski's Object, a dwarf galaxy that is interacting with a radio jet from a nearby elliptical galaxy. These observations are the first to detect molecular gas in Minkowski's Object. We estimate the range in the mass of molecular gas in Minkowski's Object assuming two different values of the ratio of the molecular gas mass to the CO luminosi… ▽ More We use ALMA to detect and image CO (1-0) emission from Minkowski's Object, a dwarf galaxy that is interacting with a radio jet from a nearby elliptical galaxy. These observations are the first to detect molecular gas in Minkowski's Object. We estimate the range in the mass of molecular gas in Minkowski's Object assuming two different values of the ratio of the molecular gas mass to the CO luminosity, $α_{\rm CO}$. For the Milky Way value of $α_{\rm CO}=4.6~M_{\odot}{\rm (K~km~s^{-1}~pc^2)^{-1}}$ we obtain a molecular gas mass of $M_{\rm H_2} =3.0 \times 10^7~M_{\odot}$, 6% of the HI gas mass. We also use the prescription of Narayanan et al. (2012) to estimate an $α_{\rm CO}=27~M_{\odot}{\rm (K~km~s^{-1}~pc^2)^{-1}}$, in which case we obtain $M_{\rm H_2} =1.8 \times 10^8~M_{\odot}$, 36% of the HI mass. The observations are consistent with previous claims of star formation being induced in Minkowski's Object via the passage of the radio jet, and it therefore being a rare local example of positive feedback from an AGN. In particular, we find highly efficient star formation, with gas depletion timescales $\sim 5\times 10^7 - 3\times 10^8$yr (for assumed values of $α_{\rm CO}=4.6$ and $27~M_{\odot}{\rm (K~km~s^{-1}~pc^2)^{-1}}$, respectively) in the upstream regions of Minkowski's Object that were struck first by the jet, and less efficient star formation downstream. We discuss the implications of this observation for models of jet induced star formation and radio mode feedback in massive galaxies. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Comments: ApJ, in press

Search for extended sources in the Galactic Plane using 6 years of Fermi-Large Area Telescope Pass 8 data above 10 GeV

Authors: The Fermi LAT Collaboration, M. Ackermann, M. Ajello, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, M. Caragiulo, P. A. Caraveo, D. Castro, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, C. C. Cheung , et al. (95 additional authors not shown)

Abstract: The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two… ▽ More The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain. △ Less

Submitted 11 April, 2018; v1 submitted 1 February, 2017; originally announced February 2017.

Comments: 33 pages, 22 figures & 3 tables. Published by The Astrophysical Journal. Available on the Fermi Science Support Center (FSSC) together with the 3FHL catalog

The second catalog of flaring gamma-ray sources from the Fermi All-sky Variability Analysis

Authors: S. Abdollahi, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, S. Buson, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman , et al. (102 additional authors not shown)

Abstract: We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data relea… ▽ More We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6σ$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2. △ Less

Submitted 12 September, 2017; v1 submitted 9 December, 2016; originally announced December 2016.

Comments: Accepted for publication by The Astrophysical Journal Supplement. Online analysis results available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/FAVA/index.php . Corresponding authors: M. Giomi, R. Buehler, D. Kocevski, and M. Ajello

Gamma-ray blazar spectra with H.E.S.S. II mono analysis: the case of PKS 2155-304 and PG 1553+113

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, F. Brun , et al. (311 additional authors not shown)

Abstract: The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft sp… ▽ More The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 600, A89 (2017)

The NuSTAR X-ray Spectrum of Hercules X-1: A Radiation-Dominated Radiative Shock

Authors: Michael T. Wolff, Peter A. Becker, Amy M. Gottlieb, Felix Fürst, Paul B. Hemphill, Diana M. Marcu-Cheatham, Katja Pottschmidt, Fritz-Walter Schwarm, Jörn Wilms, Kent S. Wood

Abstract: We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of t… ▽ More We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters. △ Less

Submitted 31 August, 2016; originally announced August 2016.

Comments: Accepted for publication in the Astrophysical Journal

Oscillatory convection and limitations of the Boussinesq approximation

Authors: Toby S Wood, Paul J Bushby

Abstract: We determine the asymptotic conditions under which the Boussinesq approximation is valid for oscillatory convection in a rapidly rotating fluid. In the astrophysically relevant parameter regime of small Prandtl number, we show that the Boussinesq prediction for the onset of convection is valid only under much more restrictive conditions than those that are usually assumed. In the case of an ideal… ▽ More We determine the asymptotic conditions under which the Boussinesq approximation is valid for oscillatory convection in a rapidly rotating fluid. In the astrophysically relevant parameter regime of small Prandtl number, we show that the Boussinesq prediction for the onset of convection is valid only under much more restrictive conditions than those that are usually assumed. In the case of an ideal gas, we recover the Boussinesq results only if the ratio of the domain height to a typical scale height is much smaller than the Prandtl number. This requires an extremely shallow domain in the astrophysical parameter regime. Other commonly-used "sound-proof" approximations generally perform no better than the Boussinesq approximation. The exception is a particular implementation of the pseudo-incompressible approximation, which predicts the correct instability threshold beyond the range of validity of the Boussinesq approximation. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Comments: 13 pages, 1 figure, accepted for publication in JFM

Searching the Gamma-ray Sky for Counterparts to Gravitational Wave Sources: Fermi GBM and LAT Observations of LVT151012 and GW151226

Authors: J. L. Racusin, E. Burns, A. Goldstein, V. Connaughton, C. A. Wilson-Hodge, P. Jenke, L. Blackburn, M. S. Briggs, J. Broida, J. Camp, N. Christensen, C. M. Hui, T. Littenberg, P. Shawhan, L. Singer, J. Veitch, P. N. Bhat, W. Cleveland, G. Fitzpatrick, M. H. Gibby, A. von Kienlin, S. McBreen, B. Mailyan, C. A. Meegan, W. S. Paciesas , et al. (116 additional authors not shown)

Abstract: We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techn… ▽ More We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: 27 pages, 1 table, 11 figures, Submitted to ApJ

Sound speed and oscillation frequencies for solar models evolved with Los Alamos ATOMIC opacities

Authors: Joyce A. Guzik, C. J. Fontes, P. Walczak, S. R. Wood, K. Mussack, E. Farag

Abstract: Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LEDCOP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compar… ▽ More Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LEDCOP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compare with models evolved using the Lawrence Livermore National Laboratory OPAL (Iglesias and Rogers 1996) opacities. We use the solar abundance mixture of Asplund et al. (2009). The new Los Alamos ATOMIC opacities have steeper opacity derivatives than those of OPAL for temperatures and densities of the solar interior radiative zone. We compare the calculated nonadiabatic solar oscillation frequencies and solar interior sound speed to observed frequencies and helioseismic inferences. The calculated sound-speed profiles are similar for models evolved using either the updated Iben evolution code (see \cite{Guzik2010}), or the MESA evolution code (Paxton et al., 2015). The LANL ATOMIC opacities partially mitigate the "solar abundance problem". △ Less

Submitted 14 May, 2016; originally announced May 2016.

Comments: Proceedings paper for IAU Focus Meeting 17, Advances in Stellar Physics from Asteroseismology, XXIXth IAU General Assembly, Honolulu, Hawaii, August 2015. Focus meeting proceedings edited by C.S. Jeffery, J.A. Guzik, and K. Pollard

Fermi LAT Gamma-ray Detections of Classical Novae V1369 Centauri 2013 and V5668 Sagittarii 2015

Authors: C. C. Cheung, P. Jean, S. N. Shore, L. Stawarz, R. H. D. Corbet, J. Knodlseder, S. Starrfield, D. L. Wood, R. Desiante, F. Longo, G. Pivato, K. S. Wood

Abstract: We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after th… ▽ More We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) gamma-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii 2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of gamma-ray onsets beginning ~2 days after their first optical peaks. Significant gamma-ray emission was found extending to 39-55 days after their initial LAT detections, with systematically fainter and longer duration emission compared to previous gamma-ray detected classical novae. These novae were distinguished by multiple bright optical peaks that encompassed the timespans of the observed gamma rays. The gamma-ray light curves and spectra of the two novae are presented along with representative hadronic and leptonic models, and comparisons to other novae detected by the LAT are discussed. △ Less

Submitted 13 May, 2016; originally announced May 2016.

Comments: 13 pages, 6 figures, 4 tables, ApJ accepted

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

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

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

Submitted 21 July, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: For the main Letter, see arXiv:1602.08492

Report number: LIGO-P1600137-v2

Journal ref: The Astrophysical Journal Supplement Series, 225:8 (15pp), 2016 July

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

Authors: M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprini , et al. (90 additional authors not shown)

Abstract: We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interac… ▽ More We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data. △ Less

Submitted 13 April, 2016; v1 submitted 12 April, 2016; originally announced April 2016.

Comments: Published by PRD, 16 pages, 11 figures, corresponding authors: F. Loparco and M. N. Mazziotta

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

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

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

Submitted 21 July, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: For Supplement, see https://arxiv.org/abs/1604.07864

Report number: LIGO-P1500227-v12

Journal ref: The Astrophysical Journal Letters, 826:L13 (8pp), 2016 July 20

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Authors: F. Acero, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi , et al. (109 additional authors not shown)

Abstract: Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Ga… ▽ More Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: Accepted for publication in ApJS

Journal ref: ApJS 223 26 (2016)

Six New Millisecond Pulsars from Arecibo Searches of Fermi Gamma-Ray Sources

Authors: H. T. Cromartie, F. Camilo, M. Kerr, J. S. Deneva, S. M. Ransom, P. S. Ray, E. C. Ferrara, P. F. Michelson, K. S. Wood

Abstract: We have discovered six radio millisecond pulsars (MSPs) in a search with the Arecibo telescope of 34 unidentified gamma-ray sources from the Fermi Large Area Telescope (LAT) 4-year point source catalog. Among the 34 sources, we also detected two MSPs previously discovered elsewhere. Each source was observed at a center frequency of 327 MHz, typically at three epochs with individual integration tim… ▽ More We have discovered six radio millisecond pulsars (MSPs) in a search with the Arecibo telescope of 34 unidentified gamma-ray sources from the Fermi Large Area Telescope (LAT) 4-year point source catalog. Among the 34 sources, we also detected two MSPs previously discovered elsewhere. Each source was observed at a center frequency of 327 MHz, typically at three epochs with individual integration times of 15 minutes. The new MSP spin periods range from 1.99 to 4.66 ms. Five of the six pulsars are in interacting compact binaries (period < 8.1 hr), while the sixth is a more typical neutron star-white dwarf binary with an 83-day orbital period. This is a higher proportion of interacting binaries than for equivalent Fermi-LAT searches elsewhere. The reason is that Arecibo's large gain afforded us the opportunity to limit integration times to 15 minutes, which significantly increased our sensitivity to these highly accelerated systems. Seventeen of the remaining 26 gamma-ray sources are still categorized as strong MSP candidates, and will be re-searched. △ Less

Submitted 20 January, 2016; originally announced January 2016.

Comments: Accepted for publication in ApJ; 8 pages, 3 figures, 5 tables

Very-high-energy gamma-rays from the Universe's middle age: detection of the z=0.940 blazar PKS 1441+25 with MAGIC

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, B. Biassuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, A. Chatterjee, R. Clavero, P. Colin, E. Colombo , et al. (229 additional authors not shown)

Abstract: The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies… ▽ More The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy. △ Less

Submitted 12 January, 2018; v1 submitted 14 December, 2015; originally announced December 2015.

Comments: Corresponding Authors: J. Becerra (josefa.becerra@nasa.gov), M. Nievas Rosillo (mnievas@ucm.es), M. Manganaro (manganaro@iac.es), F. Tavecchio (fabrizio.tavecchio@brera.inaf.it) Published in The Astrophysical Journal Letters (ApJL), 2015ApJ...815L..23A, DOI: 10.1088/2041-8205/815/2/L23

Journal ref: ApJL, 2015, 815, L23

Standing Shock Instability in Advection-Dominated Accretion Flows

Authors: Truong Le, Kent S. Wood, Michael T. Wolff, Peter A. Becker, Joy Putney

Abstract: Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura li… ▽ More Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks. △ Less

Submitted 25 March, 2016; v1 submitted 17 November, 2015; originally announced November 2015.

Comments: accepted for publication in ApJ. The paper contains 42 pages and 54 figures (v1); last revised 25 Mar. 2016 (this version, v2))

Journal ref: ApJ (2016), 819, 112

The Transient Accereting X-Ray Pulsar XTE J1946+274: Stability of the X-Ray Properties at Low Flux and Updated Orbital Solution

Authors: Diana M. Marcu-Cheatham, Katja Pottschmidt, Matthias Kühnel, Sebastian Müller, Sebastian Falkner, Isabel Caballero, Mark H. Finger, Peter J. Jenke, Colleen A. Wilson-Hodge, Felix Fürst, Victoria Grinberg, Paul B. Hemphill, Ingo Kreykenbohm, Dmitry Klochkov, Richard E. Rothschild, Yukikatsu Terada, Teruaki Enoto, Wataru Iwakiri, Michael T. Wolff, Peter A. Becker, Kent S. Wood, Jöern Wilms

Abstract: We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution us… ▽ More We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities. △ Less

Submitted 16 October, 2015; originally announced October 2015.

Comments: Accepted, 13 pages, 8 figures