Hercules X-1 is an X-ray binary system in which a neutron star is accreting matter from a close companion star through Roche lobe overflow. The observed 35-day periodicity in the X-ray light curve has been interpreted as the result of a precessing, tilted accretion disk. We use 3D hydrodynamic simulations to explore the dynamics of the tilted accretion disk and the interaction of the tidal stream with the disk. We find that, in the absence of a tidal stream, the tilted disk remains flat and precesses at a rate given by tidally driven precession. The disk tilt angle remains roughly constant with a slight orbital dependence. In the presence of the mass transfer tidal stream the disk tilt decreases over the course of several binary orbits, with a strong orbital modulation caused by the stream hitting the disk face or outer edge at different orbital phases relative to the disk tilt.

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

We use a new overset dual spherical grid system to evolve the circumstellar gas flow in the High-Mass X-Ray Binary Vela X-1 to investigate the gravitational capture of the donor star’s wind and subsequent accretion onto the neutron star companion. For a fast, relatively undistorted stellar wind (small filling factor) we find results consistent with classical Hoyle–Lyttleton accretion with the formation of a steady, nearly axisymmetric bow shock. For models with a filling factor near unity, the mass transfer is dominated by a tidal stream feeding into a thin accretion disk. Finally, for slow winds without Roche Lobe Overflow we find the formation of a wind bow shock but with enough post shock angular momentum to form a thin, stable accretion disk.

This repository contains two gzipped tarballs: plottingTools.tar.gz and plottingData.tar.gz. Within plottingData.tar.gz are the data files necessary to reproduce all of the figures in the manuscript: A Parametric Study of the SASI Comparing General Relativistic and Non-Relativistic Treatments. Within plottingTools.tar.gz are the scripts necessary to reproduce the figures. Contents of plottingTools.tar.gz: globalVariables.py: this contains some user-defined variables such as where the data is stored and what extension to use for saved figures (pdf, png, etc.) plotFigure*.py: these scripts generate the various figures in the manuscript publication.sty: this sets the formatting for the figures requirements.txt: this contains the python packages necessary to execute the scripts

We present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25 Solar masses, are all initially non-rotating, and are of two metallicities: zero and Solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains.

M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. We observed a strong Hatchett-McCray effect in M33 X-7 for the full range of wind velocities. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differ significantly from previous estimates. In particular, the masses of the components are considerably reduced to ≈38 M⊙ for the O-star donor and ≈11.4 M⊙ for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. An extended ionization zone is even present during the eclipse due to scattered X-ray photons. The X-ray ionization zone extends close to the photosphere of the donor during inferior conjunction. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region around the system emitting ∼1047 ph s−1. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, possibly resulting in a common envelope of the black hole and the O-star donor. Since the mass ratio is q ≳ 3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge.

The youngest Galactic supernova remnant, G1.9+0.3, probably the result of a Type Ia supernova, shows surprising anomalies in the distribution of its ejecta in space and velocity. In particular, high-velocity shocked iron is seen in several locations far from the remnant center, in some cases beyond prominent silicon and sulfur emission. These asymmetries strongly suggest a highly asymmetric explosion. We present high-resolution hydrodynamic simulations in two and three dimensions of the evolution from ages of 100 s to hundreds of years of two asymmetric Type Ia models, expanding into a uniform medium. At the age of G1.9+0.3 (about 100 yr), our 2D model shows almost no iron shocked to become visible in X-rays. Only in a much higher-density environment could significant iron be shocked, at which time the model's expansion speed is completely inconsistent with the observations of G1.9+0.3. Our 3D model, evolving the most asymmetric of a suite of Type Ia supernova models from Seitenzahl et al. (2013), shows some features resembling G1.9+0.3. We characterize its evolution with images of composition in three classes: C and O, intermediate-mass elements (IMEs), and iron-group elements (IGEs). From ages of 13 to 1800 yr, we follow the evolution of the highly asymmetric initial remnant as the explosion asymmetries decrease in relative strength, to be replaced by asymmetries due to evolutionary hydrodynamic instabilities. At an age of about 100 yr, our 3D model has comparable shocked masses of C+O, IMEs, and IGEs, with about 0.03 M ⊙ each. Evolutionary changes appear to be rapid enough that continued monitoring with the Chandra X-ray Observatory may show significant variations.

We report on the gravitational wave signal computed in the context of a three-dimensional simulation of a core collapse supernova explosion of a 15 Solar mass star. The simulation was performed with our neutrino hydrodynamics code Chimera. We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (< 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the SASI and high-frequency emission (> 600 Hz) emanates from the proto-neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto-neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical-acoustic origin, to it. Our dominant proto-neutron star gravitational wave emission is not well characterized by emission from surface g-modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto-neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors.

We provide a detailed description of the Chimera code, a code developed to model core collapse supernovae in multiple spatial dimensions. The core collapse supernova explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the foundational physics of the neutrino-stellar core weak interactions and the equations of state of all stellar core constituents -particularly, the nuclear equation of state associated with nucleons, both free and bound in nuclei. Chimera, by incorporating detailed neutrino transport, realistic neutrino-matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with core collapse supernovae. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. Chimera is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

We have created a new and unique database, the APOGEE–Galaxy Evolution Explorer (GALEX)–Gaia catalog, to study white dwarfs in close binaries, and, in particular, to: (1) create a minimally biased sample of white dwarf-main sequence (WDMS) systems derived from the combination of GALEX, Gaia, and the dual hemisphere, multi-epoch, high-resolution, infrared APOGEE-1 and -2 spectroscopic surveys, and (2) exploit this catalog to vastly increase the number of well-characterized WDMS systems having derived dynamical parameters (e.g., masses, separation, orbital periods, etc.) as well as precision multi-element chemical compositions by exploiting the high resolution, multi-epoch data from APOGEE. This catalog will significantly further our understanding of compact binary evolution.

Vela X is a prototypical example of a pulsar wind nebula whose morphology and detailed structure have been affected by interaction with the reverse shock of its host supernova remnant. The resulting complex of filamentary structure and mixed-in ejecta embedded in a nebula that is offset from the pulsar provides the best example we have of this middle-age state that characterizes a significant fraction of composite supernova remnants (SNRs), and perhaps all of the large-diameter pulsar wind nebulae (PWNe) seen as TeV sources. Here, we report on an XMM-Newton (hereafter XMM ) Large Project study of Vela X, supplemented by additional Chandra observations. Through broad spectral modeling, as well as detailed spectral investigations of discrete emission regions, we confirm previous studies that report evidence for ejecta material within Vela X, and show that equivalent-width variations of O vii and O viii are consistent with temperature maps within the PWN that show low-temperature regions where the projected SNR emission appears to dominate emission from the ejecta. We identify spectral variations in the nonthermal emission, with hard emission being concentrated near the pulsar. We carry out investigations of the Vela X “cocoon” structure, and with hydrodynamical simulations, show that its overall properties are consistent with structures formed in the late-phase evolution of a composite SNR expanding into a surrounding medium with a density gradient, with ejecta material being swept beyond the pulsar and compressed into an elongated structure in the direction opposite the high external density.

We investigate the effect of progenitor rotation on the standing accretion shock instability (SASI) using two-and three-dimensional hydrodynamic simulations. We find that the growth rate of the SASI is a near-linearly increasing function of the specific angular momentum in the accreting gas. Both the growth rate and the angular frequency in the two-dimensional model with cylindrical geometry agree well with previous linear stability analyses. When excited by very small random perturbations, a one-armed spiral mode dominates the small rotation rates predicted by current stellar evolution models, while progressively higher-order modes are seen as the specific angular momentum increases.

For pulsars born in supernovae, the expansion of the shocked pulsar wind nebula is initially in the freely expanding ejecta of the supernova. While the nebula is in the inner flat part of the ejecta density profile, the swept-up, accelerating shell is subject to the Rayleigh–Taylor instability. We carried out two- and three-dimensional simulations showing that the instability gives rise to filamentary structure during this initial phase but does not greatly change the dynamics of the expanding shell. The flow is effectively self-similar. If the shell is powered into the outer steep part of the density profile, the shell is subject to a robust Rayleigh–Taylor instability in which the shell is fragmented and the shocked pulsar wind breaks out through the shell. The flow is not self-similar in this phase. For a wind nebula to reach this phase requires that the deposited pulsar energy be greater than the supernova energy, or that the initial pulsar period be in the ms range for a typical 10 51 erg supernova. These conditions are satisfied by some magnetar models for Type I superluminous supernovae. We also consider the Crab Nebula, which may be associated with a low energy supernova for which this scenario applies.

In this Letter, we report on the gravitational wave signal computed in the context of an $ab$ $initio$, three-dimensional simulation of a core collapse supernova explosion, beginning with a 15M$_\odot$ star and using state-of-the-art weak interactions. The simulation was performed with our neutrino hydrodynamics code Chimera. We discuss the potential for detection of our predicted gravitational signal by the current generation of gravitational wave detectors.

We present the first three-dimensional measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12-year baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 "tufts" of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line of sight velocity, we use two different methods: a non-equilibrium ionization model fit to the strong Si and S lines in the 1.2-2.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the red or blue shift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km s^-1, with a mean of 4430 km s^-1. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km s^-1. Some Type Ia supernova explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and discuss our findings in light of various explosion models, favoring those delayed detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

The interaction between a pulsar wind nebula (PWN) and its host supernova remnant (SNR) can produce a vast array of observable structures. Asymmetry present within these structures derives from the complexity of the composite system, where many factors take turns playing a dominating hand throughout the stages of composite SNR evolution. Of particular interest are systems characterized by blastwave expansion within a nonuniform interstellar medium (ISM), which contain an active pulsar having a substantial "kick"velocity (upward of 300 km s-1), because these systems tend to produce complex morphologies. We present a numerical model that employs these and several other factors in an effort to generate asymmetry similar to that seen in various X-ray and radio observations. We find that the main parameters driving structure are ISM uniformity and total pulsar spindown energy, with secondary contributions from factors such as pulsar trajectory and initial spin-down luminosity. We also investigate the dynamics behind PWN "tails,"which may form to link active pulsars to a crushed, relic nebula as the reverse shock passes. We find that the directions of such tails are not good indicators of pulsar motion, but direction does reveal the flow of ejecta created by the passage of a reverse shock.

In this work we report briefly on the gravitational wave (GW) signal computed in the context of a self-consistent, 3D simulation of a core-collapse supernova (CCSN) explosion of a 15M_⊙ progenitor star. We present a short overview of the GW signal, including signal amplitude, frequency distribution, and the energy emitted in the form of GWs for each phase of explosion, along with neutrino luminosities, and discuss correlations between them.