Search

Observations of the S stars, the cluster of young stars in the inner 0.1 pc of the Galactic center, have been crucial in providing conclusive evidence for a supermassive black hole at the center of our galaxy. Since some of the stars have orbits less than that of a typical human lifetime, it is possible to observe multiple orbits and test the weak-field regime of general relativity. Current calculations of orbits require relatively slow and expensive computations in order to perform numerical integrations for the position and momentum of each star at each observing time. In this paper, we present a computationally efficient, first-order post-Newtonian model for the astrometric and spectroscopic data gathered for the S stars. We find that future, 30 m class telescopes—and potentially even current large telescopes with very high spectroscopic resolution—may be able to detect the Shapiro effect for an S star in the next decade or so. [ABSTRACT FROM AUTHOR]

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which shows a ring-like emission structure with the south side only slightly brighter than the north side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17° for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counterexamples, we show that low inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on. [ABSTRACT FROM AUTHOR]

The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to, but offset from, the black hole over the course of each event. We present a model for plasmoids that form during reconnection events and orbit in the coronal region around a black hole to explain these observations. We utilize general-relativistic radiative-transfer calculations that include effects from finite light travel time, plasmoid motion, particle acceleration, and synchrotron cooling, and obtain a rich structure in the flare light curves. This model can naturally account for the observed motion of the bright regions observed by the GRAVITY experiment and the offset between the center of the centroid motion and the position of the black hole. It also explains why some flares may be double peaked while others have only a single peak and uncovers a correlation between the structure in the light curve and the location of the flare. Finally, we make predictions for future observations of flares from the inner accretion flow of Sgr A* that will provide a test of this model. [ABSTRACT FROM AUTHOR]

The neutron star equation of state is now being constrained from a diverse set of multi-messenger data, including gravitational waves from binary neutron star mergers, X-ray observations of the neutron star radius, and many types of laboratory nuclear experiments. These measurements are often mapped to a common domain for comparison with one another or are used to constrain the predictions of theoretical equations of state. We explore here the statistical biases that can arise when such multi-messenger data are compared or combined across different domains. We find that placing Bayesian priors individually in each domain of measurement can lead to biased constraints. We present a new prescription for defining Bayesian priors consistently across different experiments, which will allow for robust cross-domain comparisons. Using the first two binary neutron star mergers as an example, we show that a uniform prior in the tidal deformability can produce inflated evidence for large radii, while a uniform prior in the radius points toward smaller radii. Finally, using this new prescription, we provide a status update on multi-messenger constraints on the neutron star radius. [ABSTRACT FROM AUTHOR]

The Event Horizon Telescope (EHT) is taking the first images of black holes resolved at horizon scales to measure their shadows and probe accretion physics. A promising avenue for testing the hypothesis that astrophysical black holes are described by the Kerr solution to Einstein's equations is to compare the size and shape of the shadow a black hole casts on the surrounding emission to the predictions of the Kerr metric. We develop here an efficient parametric framework to perform this test. We carry out ray-tracing simulations for several parameterized non-Kerr metrics to create a large data set of non-Kerr shadows that probe the allowed parameter space for the free parameters of each metric. We then perform principal components analysis (PCA) on this set of shadows and show that only a small number of components are needed to accurately reconstruct all shadows within the set. We further show that the amplitude of the PCA components are smoothly related to the free parameters in the metrics and, therefore, that these PCA components can be fit to EHT observations in order to place constraints on the free parameters of these metrics that will help quantify any potential deviations from the Kerr solution. [ABSTRACT FROM AUTHOR]

We use high-resolution cosmological zoom simulations with ∼200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ∼1011-2 × 1012M☼ at z = 2. We present a detailed comparison to spatially and spectrally resolved Hα and other observations of z ≈ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kiloparsec. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics—among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in Hα but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth. [ABSTRACT FROM AUTHOR]

Using stellar structure calculations in the Hartle-Thorne approximation, we derive analytic expressions connecting the ellipticity of the stellar surface to the compactness, the spin angular momentum, and the quadrupole moment of the spacetime. We also obtain empirical relations between the compactness, the spin angular momentum, and the spacetime quadrupole. Our formulae reproduce the results of numerical calculations to within a few percent and help reduce the number of parameters necessary to model the observational appearance of moderately spinning neutron stars. This is sufficient for comparing theoretical spectroscopic and timing models to observations that aim to measure the masses and radii of neutron stars and to determine the equation of state prevailing in their interiors. [ABSTRACT FROM AUTHOR]

We introduce GRay, a massively parallel integrator designed to trace the trajectories of billions of photons in a curved spacetime. This graphics-processing-unit (GPU)-based integrator employs the stream processing paradigm, is implemented in CUDA C/C++, and runs on nVidia graphics cards. The peak performance of GRay using single-precision floating-point arithmetic on a single GPU exceeds 300 GFLOP (or 1 ns per photon per time step). For a realistic problem, where the peak performance cannot be reached, GRay is two orders of magnitude faster than existing central-processing-unit-based ray-tracing codes. This performance enhancement allows more effective searches of large parameter spaces when comparing theoretical predictions of images, spectra, and light curves from the vicinities of compact objects to observations. GRay can also perform on-the-fly ray tracing within general relativistic magnetohydrodynamic algorithms that simulate accretion flows around compact objects. Making use of this algorithm, we calculate the properties of the shadows of Kerr black holes and the photon rings that surround them. We also provide accurate fitting formulae of their dependencies on black hole spin and observer inclination, which can be used to interpret upcoming observations of the black holes at the center of the Milky Way, as well as M87, with the Event Horizon Telescope. [ABSTRACT FROM AUTHOR]

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom simulations of galaxy formation down to z =2 to show that such strong self-regulation is required in the popular spherical Bondi accretion model, but that a plausible alternative model in which black hole growth is limited by galaxy-scale torques does not require self-regulation. Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along the observed scaling relations by relying only on a fixed, 5% mass retention rate onto the black hole from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other implementation details. The torque-limited model allows for high accretion rates at very early epochs (unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while by z ~ 2 it yields Eddington factors of ~1%-10%. This model also yields a less direct correspondence between major merger events and rapid phases of black hole growth. Instead, growth is more closely tied to cosmological disk feeding, which may help explain observational studies showing that, at least at z ≿ 1, active galaxies do not preferentially show merger signatures. [ABSTRACT FROM AUTHOR]

Neutron stars spinning at moderate rates (~300-600 Hz) become oblate in shape and acquire a nonzero quadrupole moment. In this paper, we calculate the profiles of atomic features from such neutron stars using a ray-tracing algorithm in the Hartle-Thorne approximation. We show that line profiles acquire cores that are much narrower than the widths expected from pure Doppler effects for a large range of observer inclinations. As a result, the effects of both the oblateness and the quadrupole moments of neutron stars need to be taken into account when aiming to measure neutron-star radii from rotationally broadened lines. Moreover, the presence of these narrow cores substantially increases the likelihood of detecting atomic lines from rapidly spinning neutron stars. [ABSTRACT FROM AUTHOR]

We develop a new numerical algorithmtomodel collisional cascades in debris disks. Because of the large dynamical range in particle masses, we solve the integro-differential equations describing erosive and catastrophic collisions in a particle-in-a-box approach, while treating the orbital dynamics of the particles in an approximate fashion. We employ a new scheme for describing erosive (cratering) collisions that yields a continuous set of outcomes as a function of colliding masses. We demonstrate the stability and convergence characteristics of our algorithm and compare it with other treatments. We show that incorporating the effects of erosive collisions results in a decay of the particle distribution that is significantly faster than with purely catastrophic collisions. [ABSTRACT FROM AUTHOR]

Time-resolved X-ray spectroscopy of thermonuclear bursts observed from low-mass X-ray binaries offer a unique tool to measure neutron-star masses and radii. In this paper, we continue our systematic analysis of all the X-ray bursts observed with Rossi X-ray Timing Explorer from X-ray binaries. We determine the events that show clear evidence for photospheric radius expansion and measure the Eddington limits for these accreting neutron stars using the bolometric fluxes attained at the touchdown moments of each X-ray burst. We employ a Bayesian technique to investigate the degree to which the Eddington limit for each source remains constant between bursts. We find that for sources with a large number of radius expansion bursts, systematic uncertainties are at a 5%-10% level. Moreover, in six sources with only pairs of Eddington-limited bursts, the distribution of fluxes is consistent with a ~10% fractional dispersion. This indicates that the spectroscopic measurements of neutron-star masses and radii using thermonuclear X-ray bursts can reach the level of accuracy required to distinguish between different neutron-star equations of state, provided that uncertainties related to the overall flux calibration of X-ray detectors are of comparable magnitude. [ABSTRACT FROM AUTHOR]

The masses and radii of low-magnetic field neutron stars can be measured by combining different observable quantities obtained from their X-ray spectra during thermonuclear X-ray bursts. One of these quantities is the apparent radius of each neutron star as inferred from the X-ray flux and spectral temperature measured during the cooling tails of bursts, when the thermonuclear flash is believed to have engulfed the entire star. In this paper, we analyze 13,095 X-ray spectra of 446 X-ray bursts observed from 12 sources in order to assess possible systematic effects in the measurements of the apparent radii of neutron stars. We first show that the vast majority of the observed X-ray spectra are consistent with blackbody functions to within a few percent. We find that most X-ray bursts follow a very well determined correlation between X-ray flux and temperature, which is consistent with the whole neutron-star surface emitting uniformly during the cooling tails. We develop a Bayesian Gaussian mixture algorithm to measure the apparent radii of the neutron stars in these sources, while detecting and excluding a small number of X-ray bursts that show irregular cooling behavior. This algorithm also provides us with a quantitative measure of the systematic uncertainties in the measurements. We find that those errors in the spectroscopic determination of neutron-star radii that are introduced by systematic effects in the cooling tails of X-ray bursts are in the range ⋍3%-8%. Such small errors are adequate to distinguish between different equations of state provided that uncertainties in the distance to each source and the absolute calibration of X-ray detectors do not dominate the error budget. [ABSTRACT FROM AUTHOR]

We carry out a new suite of cosmological radiation hydrodynamic simulations that explores the relative impacts on reionization-epoch star formation of galactic outflows and photoionization heating from a self-consistently grown extragalactic ultraviolet ionizing background (EUVB). We compare the predictions with observational constraints from the cosmic microwave background, the ultraviolet continuum luminosity function of galaxies, and the Lyα forest. By itself, an EUVB suppresses the luminosity function by less than 50% at z = 6 even if it is orders of magnitude stronger than observed. This overproduces the observed galaxy abundance by a factor of 3-5, indicating the need for an additional feedback process. We confirm that outflows readily suppress both the EUVB and the luminosity function into improved agreement with observations. Population I-II star formation can reionize the universe by z = 6 even in the presence of strong feedback from photoheating and outflows. The resulting EUVB suppresses star formation in halos with virial temperatures below 105 K but has a weaker impact on more massive halos. Nonetheless, halos with virial temperatures below 10-5 K contribute up to ~50% of all ionizing photons owing to the EUVB's inhomogeneity. Overall, star formation rate scales with halo mass Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed for halos with Mh = 108.2-1010.2 M⊙. This is a steeper dependence than is often assumed in reionization models, boosting the expected power spectrum of 21 centimeter fluctuations on large scales. The luminosity function rises steeply to at least M1600 = -13 even in models that treat both outflows and an EUVB. indicating that reionization was driven by faint galaxies (M1600 ≥ -15) that have not yet been observed. Outflows and an EUVB interfere with each other's feedback effects in two ways. Outflows weaken the EUVB, limiting Jeans suppression of low-mass halos: this leads to overall de-amplification of suppression at early times (z > 8). Meanwhile, they amplify each other's impact on more massive halos, leading to overall amplification of suppression at later times. Our models cannot simultaneously explain observations of galaxies, the cosmic microwave background, and the intergalactic medium. Correcting for dynamic range limitations and adjusting our physical treatments will alleviate discrepancies, but observations may still require additional physical scalings such as a mass-dependent ionizing escape traction. [ABSTRACT FROM AUTHOR]

The nuclear symmetry energy plays a role in determining both the nuclear properties of terrestrial matter as well as the astrophysical properties of neutron stars. The first measurement of the neutron star tidal deformability, from gravitational-wave event GW170817, provides a new way of probing the symmetry energy. In this work, we report on new constraints on the symmetry energy from GW170817. We focus in particular on the low-order coefficients: namely, the value of the symmetry energy at the nuclear saturation density, S0, and the slope of the symmetry energy, L0. We find that the gravitational-wave data are relatively insensitive to S0, but that they depend strongly on L0 and point to lower values of L0 than have previously been reported, with a peak likelihood near L0 ∼ 23 MeV. Finally, we use the inferred posteriors on L0 to derive new analytic constraints on higher-order nuclear terms. [ABSTRACT FROM AUTHOR]

Electron acceleration during magnetic reconnection is thought to play a key role in time-variable high-energy emission from astrophysical systems. By means of particle-in-cell simulations of transrelativistic reconnection, we investigate electron injection and acceleration mechanisms in low-β electron–proton plasmas. We set up a diversity of density and field structures (e.g., X-points and plasmoids) by varying the guide field strength and choosing whether to trigger reconnection or let it spontaneously evolve. We show that the number of X-points and plasmoids controls the efficiency of electron acceleration, with more X-points leading to a higher efficiency. Using on-the-fly acceleration diagnostics, we also show that the nonideal electric fields associated with X-points play a critical role in the first stages of electron acceleration. As a further diagnostic, we include two populations of test particles that selectively experience only certain components of electric fields. We find that the out-of-plane component of the parallel electric field determines the hardness of the high-energy tail of the electron energy distribution. These results further our understanding of electron acceleration in this regime of magnetic reconnection and have implications for realistic models of black hole accretion flows. [ABSTRACT FROM AUTHOR]

Observations of isolated neutron stars place constraints on the equation of state (EOS) of cold, neutron-rich matter, while nuclear physics experiments probe the EOS of hot, symmetric matter. Many dynamical phenomena, such as core-collapse supernovae, the formation and cooling of proto-neutron stars, and neutron star mergers, lie between these two regimes and depend on the EOS at finite temperatures for matter with varying proton fractions. In this paper, we introduce a new framework to accurately calculate the thermal pressure of neutron–proton–electron matter at arbitrary density, temperature, and proton fraction. This framework can be expressed using a set of five physically motivated parameters that span a narrow range of values for realistic EOS and are able to capture the leading-order effects of degenerate matter on the thermal pressure. We base two of these parameters on a new approximation of the Dirac effective mass, with which we reproduce the thermal pressure to within ≲30% for a variety of realistic EOS at densities of interest. Three additional parameters, which are based on the behavior of the symmetry energy near the nuclear saturation density, allow us to extrapolate any cold EOS in β-equilibrium to arbitrary proton fractions. Our model thus allows a user to extend any cold nucleonic EOS, including piecewise polytropes, to arbitrary temperature and proton fraction for use in calculations and numerical simulations of astrophysical phenomena. We find that our formalism is able to reproduce realistic finite-temperature EOS with errors of ≲20% and offers a 1–3 orders-of-magnitude improvement over existing ideal-fluid models. [ABSTRACT FROM AUTHOR]

The Neutron-star Interior Composition ExploreR (NICER) is collecting data to measure the radii of neutron stars by observing the pulsed emission from their surfaces. The primary targets are isolated, rotation-powered pulsars, in which the surface polar caps are heated by bombardment from magnetospheric currents of electrons and positrons. We investigate various stopping mechanisms for the beams of particles that bombard the atmosphere and calculate the heat deposition, the atmospheric temperature profiles, and the energy spectra and beaming of the emerging radiation. We find that low-energy particles with γ ∼ 2–10 deposit most of their energy in the upper regions of the atmosphere, at low optical depth, resulting in beaming patterns that are substantially different from those of deep-heated, radiative equilibrium models. Only particles with energies γ ≳ 50 penetrate to high optical depths and fulfill the conditions necessary for a deep-heating approximation. We discuss the implications of our work for modeling the pulse profiles from rotation-powered pulsars and for the inference of their radii with NICER observations. [ABSTRACT FROM AUTHOR]

Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer–Lindquist (BL) coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr–Schild (KS) coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian KS coordinates that outperforms calculations that use the seemingly simpler equations in BL coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multicore CPUs as well as on a wide range of graphics processing units hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics. [ABSTRACT FROM AUTHOR]

We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope. We also show that the simulations in the spatial domain can themselves be compactly represented with a PCA-derived basis of eigenimages, which allows for detailed comparisons to be made between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images. [ABSTRACT FROM AUTHOR]

Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late 2017 September in a giant outburst with a peak luminosity of 2 × 1039(d/7 kpc)2 erg s−1 (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 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 ∼ 1038 erg s−1, 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 turnover 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 ∼ 1013 G. [ABSTRACT FROM AUTHOR]

Nonthermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection—where the magnetization σ, defined as the ratio of magnetic energy density to enthalpy density, is ∼1—is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. We focus on the dependence of the electron energy spectrum on σ and the proton β (the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is nonthermal and can be modeled as a power law. At low β, the slope, p, is independent of β and hardens with increasing σ as . Electrons are primarily accelerated by the nonideal electric field at X-points, either in the initial current layer or in current sheets generated between merging magnetic islands. At higher values of β, the electron power law steepens, and the electron spectrum eventually approaches a Maxwellian distribution for all values of σ. At values of β near βmax ≈ 1/4σ, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. These particles are accelerated via a Fermi-like process by bouncing between the reconnection outflow and a stationary magnetic island [ABSTRACT FROM AUTHOR]

The Event Horizon Telescope will generate horizon scale images of the black hole in the center of the Milky Way, Sgr A*. Image reconstruction using interferometric visibilities rests on the assumption of a stationary image. We explore the limitations of this assumption using high-cadence disk- and jet-dominated GRMHD simulations of Sgr A*. We also employ analytic models that capture the basic characteristics of the images to understand the origin of the variability in the simulated visibility amplitudes. We find that, in all simulations, the visibility amplitudes for baselines oriented parallel and perpendicular to the spin axis of the black hole follow general trends that do not depend strongly on accretion-flow properties. This suggests that fitting Event Horizon Telescope observations with simple geometric models may lead to a reasonably accurate determination of the orientation of the black hole on the plane of the sky. However, in the disk-dominated models, the locations and depths of the minima in the visibility amplitudes are highly variable and are not related simply to the size of the black hole shadow. This suggests that using time-independent models to infer additional black hole parameters, such as the shadow size or the spin magnitude, will be severely affected by the variability of the accretion flow. [ABSTRACT FROM AUTHOR]

The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. We find strong agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ∼1.4 M⊙) originates from progenitors with MZAMS ≈ 9–18 M⊙. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ∼1.8 M⊙) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new, promising evidence that these stellar evolution and explosion models capture the majority of relevant stellar, nuclear, and explosion physics involved in the formation of compact objects. [ABSTRACT FROM AUTHOR]

Non-ideal magnetohydrodynamic (MHD) effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal evolution (SANE) disks, we find that in the reconnection sites, the plasma beta ranges from 0.1 to 1000, the magnetization ranges from 10−4 to 1, and the guide fields are weak compared with the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from 10−2 to 103, magnetizations from 10−3 to 10, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. We also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way. [ABSTRACT FROM AUTHOR]

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. Here, we present a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parameterized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, , can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ∼30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass–radius relation can lead to biases of nearly 1 km in radius, toward larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias. [ABSTRACT FROM AUTHOR]

Closure phases along different baseline triangles carry a large amount of information regarding the structures of the images of black holes in interferometric observations with the Event Horizon Telescope. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A* to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of the small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the funnel-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the u–v space that show a high level of variability will need to be identified and excised. [ABSTRACT FROM AUTHOR]

The Event Horizon Telescope (EHT) is a millimeter-wavelength, very-long-baseline interferometry (VLBI) experiment that is capable of observing black holes with horizon-scale resolution. Early observations have revealed variable horizon-scale emission in the Galactic Center black hole, Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools that explicitly consider the variability in both the data and the models. We develop here a Bayesian method to compare time-resolved simulation images to variable VLBI data, in order to infer model parameters and perform model comparisons. We use mock EHT data based on GRMHD simulations to explore the robustness of this Bayesian method and contrast it to approaches that do not consider the effects of variability. We find that time-independent models lead to offset values of the inferred parameters with artificially reduced uncertainties. Moreover, neglecting the variability in the data and the models often leads to erroneous model selections. We finally apply our method to the early EHT data on Sgr A*. [ABSTRACT FROM AUTHOR]

The Neutron-star Interior Composition Explorer is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated. [ABSTRACT FROM AUTHOR]

The merger of a neutron star (NS) binary may result in the formation of a rapidly spinning magnetar. The magnetar can potentially survive for seconds or longer as a supramassive NS before collapsing to a black hole if, indeed, it collapses at all. During this process, a fraction of the magnetar’s rotational energy of ∼1053 erg is transferred via magnetic spin-down to the surrounding ejecta. The resulting interaction between the ejecta and the surrounding circumburst medium powers a year-long or greater synchrotron radio transient. We present a search for radio emission with the Very Large Array following nine short-duration gamma-ray bursts (GRBs) at rest-frame times of ≈1.3–7.6 yr after the bursts, focusing on those events that exhibit early-time excess X-ray emission that may signify the presence of magnetars. We place upper limits of ≲18–32 μJy on the 6.0 GHz radio emission, corresponding to spectral luminosities of ≲(0.05–8.3) × 1039 erg s−1. Comparing these limits to the predicted radio emission from a long-lived remnant and incorporating measurements of the circumburst densities from broadband modeling of short GRB afterglows, we rule out a stable magnetar with an energy of 1053 erg for half of the events in our sample. A supramassive remnant that injects a lower rotational energy of 1052 erg is ruled out for a single event, GRB 050724A. This study represents the deepest and most extensive search for long-term radio emission following short GRBs to date, and thus the most stringent limits placed on the physical properties of magnetars associated with short GRBs from radio observations. [ABSTRACT FROM AUTHOR]

We present Chandra/ACIS-S subarray observations of the quiescent neutron star (NS) low-mass X-ray binaries X7 and X5 in the globular cluster 47 Tuc. The large reduction in photon pile-up compared to previous deep exposures enables a substantial improvement in the spectroscopic determination of the NS radius and mass of these NSs. Modeling the thermal emission from the NS surface with a non-magnetized hydrogen atmosphere and accounting for numerous sources of uncertainties, we obtain for the NS in X7 a radius of km for an assumed stellar mass of M = 1.4 M⊙ (68% confidence level). We argue, based on astrophysical grounds, that the presence of a He atmosphere is unlikely for this source. Due to the excision of data affected by eclipses and variable absorption, the quiescent low-mass X-ray binary X5 provides less stringent constraints, leading to a radius of km, assuming a hydrogen atmosphere and a mass of M = 1.4 M⊙. When combined with all existing spectroscopic radius measurements from other quiescent low-mass X-ray binaries and Type I X-ray bursts, these measurements strongly favor radii in the 9.9–11.2 km range for a ∼1.5 M⊙ NS and point to a dense matter equation of state that is somewhat softer than the nucleonic ones that are consistent with laboratory experiments at low densities. [ABSTRACT FROM AUTHOR]

The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state (EOS). One way to facilitate the mapping of observables to the EOS is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EOS. We use mock EOS that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables of mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EOS with five fiducial densities between ∼1–8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EOS to within less than 0.5 km for the extreme mock EOS and to within less than 0.12 km for 95% of a sample of 42 proposed, physically motivated EOS. Such a parametrization is also able to reproduce the maximum mass to within 0.04 and the moment of inertia of a 1.338 neutron star to within less than 10% for 95% of the proposed sample of EOS. [ABSTRACT FROM AUTHOR]

Many techniques for measuring neutron star radii rely on absolute flux measurements in the X-rays. As a result, one of the fundamental uncertainties in these spectroscopic measurements arises from the absolute flux calibrations of the detectors being used. Using the stable X-ray burster, GS 1826–238, and its simultaneous observations by Chandra HETG/ACIS-S and RXTE/PCA as well as by XMM-Newton EPIC-pn and RXTE/PCA, we quantify the degree of uncertainty in the flux calibration by assessing the differences between the measured fluxes during bursts. We find that the RXTE/PCA and the Chandra gratings measurements agree with each other within their formal uncertainties, increasing our confidence in these flux measurements. In contrast, XMM-Newton EPIC-pn measures 14.0 ± 0.3% less flux than the RXTE/PCA. This is consistent with the previously reported discrepancy with the flux measurements of EPIC-pn, compared with EPIC MOS1, MOS2, and ACIS-S detectors. We also show that any intrinsic time-dependent systematic uncertainty that may exist in the calibration of the satellites has already been implicity taken into account in the neutron star radius measurements. [ABSTRACT FROM AUTHOR]

Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform general relativistic magnetohydrodynamic simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts. [ABSTRACT FROM AUTHOR]

We present a comprehensive study of interstellar X-ray extinction using the extensive Chandra supernova remnant (SNR) archive and use our results to refine the empirical relation between the hydrogen column density and optical extinction. In our analysis, we make use of the large, uniform data sample to assess various systematic uncertainties in the measurement of the interstellar X-ray absorption. Specifically, we address systematic uncertainties that originate from (i) the emission models used to fit SNR spectra; (ii) the spatial variations within individual remnants; (iii) the physical conditions of the remnant such as composition, temperature, and non-equilibrium regions; and (iv) the model used for the absorption of X-rays in the interstellar medium. Using a Bayesian framework to quantify these systematic uncertainties, and combining the resulting hydrogen column density measurements with the measurements of optical extinction toward the same remnants, we find the empirical relation NH = (2.87 ± 0.12) × 1021AV cm−2, which is significantly higher than the previous measurements. [ABSTRACT FROM AUTHOR]

We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of neutron stars, we place strong and quantitative constraints on the properties of the equation of state between times the nuclear saturation density. We find that around the preferred equation of state predicts radii between 10.1 and 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential. [ABSTRACT FROM AUTHOR]

The half opening angle of a Kerr black hole shadow is always equal to (5 ± 0.2)GM/Dc2, where M is the mass of the black hole and D is its distance from the Earth. Therefore, measuring the size of a shadow and verifying whether it is within this 4% range constitutes a null hypothesis test of general relativity. We show that the black hole in the center of the Milky Way, Sgr A*, is the optimal target for performing this test with upcoming observations using the Event Horizon Telescope (EHT). We use the results of optical/IR monitoring of stellar orbits to show that the mass-to-distance ratio for Sgr A* is already known to an accuracy of ∼4%. We investigate our prior knowledge of the properties of the scattering screen between Sgr A* and the Earth, the effects of which will need to be corrected for in order for the black hole shadow to appear sharp against the background emission. Finally, we explore an edge detection scheme for interferometric data and a pattern matching algorithm based on the Hough/Radon transform and demonstrate that the shadow of the black hole at 1.3 mm can be localized, in principle, to within ∼9%. All these results suggest that our prior knowledge of the properties of the black hole, of scattering broadening, and of the accretion flow can only limit this general relativistic null hypothesis test with EHT observations of Sgr A* to ≲10%. [ABSTRACT FROM AUTHOR]

We explore the variability properties of long, high-cadence general relativistic magnetohydrodynamic (GRMHD) simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both standard and normal evolution (SANE) and magnetically arrested disk (MAD) simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3 mm. We find that the SANE models produce short-timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, MAD models generate only slow variability at lower flux levels. Neither set of models shows any X-ray flares, which most likely indicates that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The SANE models show strong, short-lived millimeter/infrared (IR) flares, with short (≲1 hr) time lags between the millimeter and IR wavelengths, that arise from the combination of short-lived magnetic flux tubes and strong-field gravitational lensing near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts. [ABSTRACT FROM AUTHOR]

We calculate the effects of spot size on pulse profiles of moderately rotating neutron stars. Specifically, we quantify the bias introduced in radius measurements from the common assumption that spots are infinitesimally small. We find that this assumption is reasonable for spots smaller than 10°–18° and leads to errors that are ≤10% in the radius measurement, depending on the location of the spot and the inclination of the observer. We consider the implications of our results for neutron star radius measurements with the upcoming and planned X-ray missions NICER and LOFT. We calculate the expected spot size for different classes of sources and investigate the circumstances under which the assumption of a small spot is justified. [ABSTRACT FROM AUTHOR]

Measuring neutron star radii with spectroscopic and timing techniques relies on the combination of multiple observables to break the degeneracies between the mass and radius introduced by general relativistic effects. Here, we explore a previously used frequentist and a newly proposed Bayesian framework to obtain the most likely value and the uncertainty in such a measurement. We find that for the expected range of masses and radii and for realistic measurement errors, the frequentist approach suffers from biases that are larger than the accuracy in the radius measurement required to distinguish between the different equations of state. In contrast, in the Bayesian framework, the inferred uncertainties are larger, but the most likely values do not suffer from such biases. We also investigate ways of quantifying the degree of consistency between different spectroscopic measurements from a single source. We show that a careful assessment of the systematic uncertainties in the measurements eliminates the need for introducing ad hoc biases, which lead to artificially large inferred radii. [ABSTRACT FROM AUTHOR]

A correction to the article "The Power of Imaging: Constraining the Plasma Properties of GRMHD Simulations Using EHT Observations of SGR A*," by Chik-Kwan Chan and colleagues published in a 2015 issue of "The Astrophysical Journal" is presented.

We model the pulse profiles and the phase-resolved spectra of the anomalous X-ray pulsar 1E 1048.1–5937 obtained with XMM-Newton to map its surface temperature distribution during an active and a quiescent epoch. We develop and apply a model that takes into account the relevant physical and geometrical effects on the neutron star surface, magnetosphere, and spacetime. Using this model, we determine the observables at infinity as a function of pulse phase for different numbers and sizes of hot spots on the surface. We show that the pulse profiles extracted from both observations can be modeled with a single hot spot and an antipodal cool component. The size of the hot spot changes from ≈80° in 2007, three months after the onset of a dramatic flux increase, to ≈30° during the quiescent observation in 2011, when the pulsed fraction returned to the pre-outburst ≈65% level. For the 2007 observation, we also find that a model consisting of a single 0.4 keV hot spot with a magnetic field strength of 1.8 × 1014 G accounts for the spectra obtained at three different pulse phases but underpredicts the flux at the pulse minimum, where the contribution to the emission from the cooler component is non-negligible. The inferred temperature of the spot stays approximately constant between different pulse phases, in agreement with a uniform temperature, single hot spot model. These results suggest that the emitting area grows significantly during outbursts but returns to its persistent and significantly smaller size within a timescale of a few years. [ABSTRACT FROM AUTHOR]