SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
NASA's Laboratory for High Energy Astrophysics and the Dept. of Physics and Astrophysics at the Univ. of Md. collaberated on a graduate level course with this title. This publication is an edited version of notes used as the course text. Topics include stellar evolution, pulsars, binary stars, X-ray signatures, gamma ray sources, and temporal analysis of X-ray data.
The prime objectives of the Rossi X-ray Timing Explorer (RXTE) are the study of astrophysical compact objects: black holes (galactic and extragalactic), many types of neutron stars, and accreting white dwarfs. RXTE's capability for rapid timing is unique among all past and present X-ray observatories. It has enabled qualitatively new discoveries about dynamical timescale phenomena related to neutron stars and black holes. These phenomena probe basic physics in the most extreme environments of gravity, density, and magnetic fields. The observations to date are a public archive for studying questions not addressed in the initial studies. Proposed future observations emphasize the discovery of additional millisecond pulsars and finding the answers to questions about low mass X-ray binaries and black holes.
The prime scientific objectives of the Rossi X-Ray Timing Explorer (RXTE) were the study of astrophysical compact objects: black holes (galactic and extragalactic), many types of neutron stars, and accreting white dwarfs. RXTE was successful in achieving its original observing objectives of large area and high time resolution observations with broadband (2-200 keV) spectra, scheduled flexibly enough to enable observations of targets of opportunity on any timescale greater than a few hours. These capabilities enabled qualitatively new discoveries about dynamical timescale phenomena related to neutron stars and black holes, phenomena which probe basic physics in the most extreme environments of gravity, density, and magnetic fields. RXTE has extended its lifetime by applying the proportional counter area selectively and maintains schedule flexibility by making use of the distribution of targets around the sky. Proposed future observations emphasize opportunity to discover and study additional millisecond pulsars, pursue the high frequency quasi-periodic oscillations in black hole transients, and connect high frequency phenomena with longer term characteristics. RXTE will continue to strongly support, for both galactic and extragalactic targets, combining RXTE observations with other wavelengths (from IR to TeV) or with other capabilities, such as high spectral resolution.
General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.
Nonlinear wave-coupling is studied in a multispecies degenerate astrophysical plasma consisting of two electron species (at different temperatures): a highly degenerate main component plus a smaller classical relativistic flow immersed in a static neutralizing ion background. It is shown that the high frequency electromagnetic waves through their strong nonlinear interactions with the electron-acoustic waves [sustained by a multielectron component (degenerate) plasma surrounding a compact astrophysical object] can scatter to lower frequencies so that the radiation observed faraway will be spectrally shifted downward. It is also shown that, under definite conditions, the electromagnetic waves could settle into stationary solitonic states. It is expected that the effects of such structures may persist as detectable signatures in forms of modulated micropulses in the radiation observed far away from the accreting compact object. Both these effects will advance our abilities to interpret the radiation coming out of the compact objects.
Aims. In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. Methods. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Results. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15–30% discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport.
Progress made in the following areas is discussed: low energy ion and electron experiments; instrument design for current experiments; magnetospheric measurement of particles; ion measurement in the earth plasma sheet; abundance measurement; X-ray data acquisition; high energy physics; extragalactic astronomy; compact object astrophysics; planetology; and high energy photon detector technology.
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope, and high- precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD, and KORAL. Agreement among the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.
We propose to use gravitational waves from extreme mass ratio inspirals (EMRI), composed of a boson star and a supermassive black hole in the center of galaxies, as a new method to search for boson stars. Gravitational waves from EMRI have the advantage of being long-lasting within the frequency band of future space-based interferometer gravitational wave detectors and can accumulate large signal-to-noise ratio (SNR) for very sub-solar mass boson stars. Compared to gravitational waves from boson star binaries, which fall within the LIGO band, we find that much larger ranges of the mass and compactness of boson stars, as well as the underlying particle physics parameter space, can be probed by EMRI. We take tidal disruption of the boson stars into account and distinguish those which dissolve before the inner-most-stable-circular-orbit (ISCO) and those which dissolve after it. Due to the large number of cycles recorded, EMRIs can lead to a very precise mass determination of the boson star and distinguish it from standard astrophysical compact objects in the event of a signal. Tidal effects in inspiralling binary systems, as well as possible correlated electromagnetic signals, can also serve as potential discriminants.
The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves (GWs) and likely production sites for heavy-element nucleosynthesis by means of rapid neutron capture (the r -process). Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of long-duration GRBs associated with compact object mergers and contains a kilonova similar to AT2017gfo, associated with the GW merger GW170817. We obtained James Webb Space Telescope (JWST) mid-infrared imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns, which we interpret as tellurium (atomic mass A = 130) and a very red source, emitting most of its light in the mid-infrared owing to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r -process elements across a broad atomic mass range and play a central role in heavy-element nucleosynthesis across the Universe.
Multimessenger observations of binary neutron star mergers can provide valuable information on the nuclear equation of state (EOS). Here, we investigate the extent to which electromagnetic observations of the associated kilonovae allow us to place constraints on the EOS. For this, we use state-of-the-art three-dimensional general-relativistic magnetohydrodynamics simulations and detailed nucleosynthesis modeling to connect properties of observed light curves to properties of the accretion disk, and hence, the EOS. Using our general approach, we use multimessenger observations of GW170817/AT2017gfo to study the impact of various sources of uncertainty on inferences of the EOS. We constrain the radius of a 1.4M ⊙ neutron star to lie within 10.30 ≤ R 1.4 ≤ 13.0 km and the maximum mass to be M TOV ≤ 3.06M ⊙ .
Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.
We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.
The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.
Magnetars are the most highly magnetized neutron stars in the cosmos (with magnetic field 10 13 –10 15 G). Giant flares from magnetars are rare, short-duration (about 0.1 s) bursts of hard X-rays and soft γ rays. Owing to the limited sensitivity and energy coverage of previous telescopes, no magnetar giant flare has been detected at gigaelectronvolt (GeV) energies. Here, we report the discovery of GeV emission from a magnetar giant flare on 15 April 2020 (refs. 3,4 and A. J. Castro-Tirado et al., manuscript in preparation). The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected GeV γ rays from 19 s until 284 s after the initial detection of a signal in the megaelectronvolt (MeV) band. Our analysis shows that these γ rays are spatially associated with the nearby (3.5 megaparsecs) Sculptor galaxy and are unlikely to originate from a cosmological γ-ray burst. Thus, we infer that the γ rays originated with the magnetar giant flare in Sculptor. Furthermore, we suggest that the GeV signal is generated by an ultra-relativistic outflow that first radiates the prompt MeV-band photons, and then deposits its energy far from the stellar magnetosphere. After a propagation delay, the outflow interacts with environmental gas and produces shock waves that accelerate electrons to very high energies; these electrons then emit GeV γ rays as optically thin synchrotron radiation. This observation implies that a relativistic outflow is associated with the magnetar giant flare, and suggests the possibility that magnetars can power some short γ-ray bursts.
The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star’s size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hypermassive neutron star, we find that the radius of a 1.4 M⊙ neutron star is $R_{1.4M_⊙} = 11.0^{+0.9}_{-0.6} km$ (90% credible interval). Finally, using this constraint, we show that neutron stars are unlikely to be disrupted in neutron star–black hole mergers; subsequently, such events will not produce observable electromagnetic emission.
LIGO and Virgo’s third observing run revealed the first neutron star–black hole (NSBH) merger candidates in gravitational waves. These events are predicted to synthesize r-process elements creating optical/near-infrared ‘kilonova’ emission. The joint gravitational wave and electromagnetic detection of an NSBH merger could be used to constrain the equation of state of dense nuclear matter, and independently measure the local expansion rate of the Universe. Here, we present the optical follow-up and analysis of two of the only three high-significance NSBH merger candidates detected to date, S200105ae and S200115j, with the Zwicky Transient Facility. The Zwicky Transient Facility observed ~48% of S200105ae and ~22% of S200115j’s localization probabilities, with observations sensitive to kilonovae brighter than –17.5 mag fading at 0.5 mag d –1 in the g- and r-bands; extensive searches and systematic follow-up of candidates did not yield a viable counterpart. We present state-of-the-art kilonova models tailored to NSBH systems that place constraints on the ejecta properties of these NSBH mergers. Here, we show that with observed depths of apparent magnitude ~22 mag, attainable in metre-class, wide-field-of-view survey instruments, strong constraints on ejecta mass are possible, with the potential to rule out low mass ratios, high black hole spins and large neutron star radii.
We report a star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events have now been identified in optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three tidal disruption events. The continuum polarization appears independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code, we simulate the polarization degree as a function of disk mass and viewing angle and we show that the predicted levels are compatible with the observations for extended reprocessing envelopes of ~1,000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for tidal disruption event models, and opens an important new line of exploration in the study of tidal disruption events.