Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pulsar astronomy”

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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Very-high-energy Emission from Pulsars

Abstract Air-Cherenkov telescopes have detected pulsations at energies above 50 GeV from a growing number of Fermi pulsars. These include the Crab, Vela, PSR B1706−44, and Geminga, with the first two having pulsed detections above 1 TeV. In some cases, there appears to be very-high-energy (VHE) emission that is an extension of the Fermi spectra to high energies, while in other cases, additional higher-energy spectral components that require a separate emission mechanism may be present. We present results of broadband spectral modeling using global magnetospheric fields and multiple emission mechanisms that include synchro-curvature (SC) and inverse Compton scattered (ICS) radiation from accelerated particles (primaries) and synchrotron self-Compton (SSC) emission from lower-energy pairs. Our models predict three distinct VHE components: SC from primaries whose high-energy tail can extend to 100 GeV, SSC from pairs that can extend to several TeV, and ICS from primary particles accelerated in the current sheet that scatter pair synchrotron radiation, which appears beyond 10 TeV. Our models suggest that H.E.S.S.-II and MAGIC have detected the high-energy tail of the primary SC component that produces the Fermi spectrum in Vela, Geminga, and PSR B1706−44. We argue that the ICS component peaking above 10 TeV from Vela has been seen by H.E.S.S. Detection of this emission component from the Crab and other pulsars is possible with the High Altitude Water Cherenkov Observatory and Cherenkov Telescope Array, and will directly measure the maximum particle energy in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Eel Pulsar Wind Nebula: A PeVatron-candidate Origin for HAWC J1826-128 and HESS J182-130

HAWC J1826–128 is one of the brightest Galactic TeV γ-ray sources detected by the High Altitude Water Cherenkov (HAWC) observatory, with photon energies extending up to nearly ~100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826–130 and the "Eel" pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826–1256. In the X-ray band, Chandra and XMM-Newton revealed that the Eel PWN is composed of both a compact nebula (~15'') and diffuse X-ray emission (~6' × 2') extending away from the pulsar. Our NuSTAR observation detected hard X-ray emission from the compact PWN up to ~20 keV and evidence of the synchrotron burn-off effect. In addition to the spatial coincidence between HESS J1826–130 and the diffuse X-ray PWN, our multiwavelength spectral energy distribution (SED) analysis using X-ray and γ-ray data establishes a leptonic origin of the TeV emission associated with the Eel PWN. Furthermore, our evolutionary PWN SED model suggests (1) a low PWN B-field of ~1 μG, (2) a significantly younger pulsar age (t ~ 5.7 kyr) than the characteristic age (τ = 14.4 kyr), and (3) a maximum electron energy of ${E}_{\max }=2$ PeV. The low B-field, as well as the putative supersonic motion of the pulsar, may account for the asymmetric morphology of the diffuse X-ray emission. Our results suggest that the Eel PWN may be a leptonic PeVatron particle accelerator powered by the ~6 kyr old pulsar PSR J1826–1256 with a spin-down power of 3.6 × 10 36 erg s –1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Trace Anomaly as Signature of Conformality in Neutron Stars

We discuss an interpretation that a peak in the sound velocity in neutron star matter, as suggested by the observational data, signifies strongly coupled conformal matter. Here, the normalized trace anomaly is a dimensionless measure of conformality leading to the derivative and the nonderivative contributions to the sound velocity. We find that the peak in the sound velocity is attributed to the derivative contribution from the trace anomaly that steeply approaches the conformal limit. Smooth continuity to the behavior of high-density QCD implies that the matter part of the trace anomaly may be positive definite. We discuss a possible implication of the positivity condition of the trace anomaly on the M–R relation of the neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Multipolar Magnetic Field of the Millisecond Pulsar PSR J0030+0451

Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed at constraining the neutron star mass and radius, has inferred surface hot spots (the magnetic polar caps) that imply significantly nondipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole-plus-quadrupole components using a static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. and Riley et al., we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov Chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions of the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous γ-ray light curves following Kalapotharakos et al.; these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the γ-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and γ-ray light curves. We conclude with a discussion of the broader implications of our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically dominated (${\sigma }_{0}\gg 1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim {\sigma }_{0}{\gamma }_{\mathrm{th}0}$ (here, ${\gamma }_{\mathrm{th}0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are preferentially oriented in the plane perpendicular to the field. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu {F}_{\nu }\propto {\nu }^{s}$ that is much harder than for isotropic particles. Remarkably, for ${\sigma }_{0}\gg 1$ we find that the angle-integrated spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta {B}_{\mathrm{rms}0}^{2}/{B}_{0}^{2}\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Further, our findings may help explain the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of pulsar wind nebulae.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

From Feast to Famine: A Systematic Study of Accretion onto Oblique Pulsars with 3D GRMHD Simulations

Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

NuSTAR Observations of G11.2–0.3

In this work, we present the hard X-ray view of the pulsar wind nebula in G11.2-0.3 and its central pulsar powered pulsar J1811-1925 as seen by NuSTAR. We complement the data with Chandra for a more complete picture and confirm the existence of a hard, power-law component in the shell with photon index ${\rm{\Gamma }}=2.1\pm 0.1$, which we attribute to synchrotron emission. Our imaging observations of the shell show a slightly smaller radius at higher energies, consistent with Chandra results, and we find shrinkage as a function of increased energy along the jet direction, indicating that the electron outflow in the PWN may be simpler than that seen in other young PWNe. Combining NuSTAR with INTEGRAL, we find that the pulsar spectrum can be fit by a power law with ${\rm{\Gamma }}=1.32\pm 0.07$ up to 300 keV without evidence of curvature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Axion-photon conversion in strongly magnetised plasmas

Axion dark matter can resonantly convert to photons in the magnetosphere of neutron stars, possibly giving rise to radio signals observable on Earth. This method for the indirect detection of axion dark matter has recently received significant attention in the literature. The calculation of the radio signal is complicated by a number of effects; most importantly, the gravitational infall of the axions onto the neutron star accelerates them to semi-relativistic speed, and the neutron star magnetosphere is highly anisotropic. Both of these factors complicate the calculation of the conversion of axions to photons. In this work, we present the first fully three-dimensional calculation of the axion-photon conversion in highly magnetised anisotropic media. Depending on the axion trajectory, this calculation leads to orders-of-magnitude differences in the conversion compared to the simplified one-dimensional calculation used so far in the literature, altering the directionality of the produced photons. Our results will have important implications for the radio signal one would observe in a telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Total energy in supernova neutrinos and the tidal deformability and binding energy of neutron stars

The energy radiated in supernova neutrinos is a fundamental quantity that is closely related to the gravitational binding energy of a neutron star. Recently the tidal deformability of neutron stars was constrained by gravitational wave observations. By considering several equations of state, we find a strong correlation between the tidal deformability and neutron star binding energy. We use this correlation to sharpen predictions of the binding energy of neutron stars and the total neutrino energy in supernovae. We find a minimum binding energy for a neutron star formed in a supernova of ~ 1.5 × 10 53 ergs . Should the neutrino energy in a supernova be significantly below this value, it would strongly suggest new unobserved particles are carrying away some of the supernova energy. Alternatively, if the neutrino energy is observed above ~ 6 × 10 53 ergs , it would strongly imply the formation of a (perhaps surprisingly) massive neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evidence for a high-energy tail in the gamma-ray spectra of globular clusters

Millisecond pulsars are very likely the main source of gamma-ray emission from globular clusters. However, the relative contributions of two separate emission processes – curvature radiation from millisecond pulsar magnetospheres versus inverse Compton emission from relativistic pairs launched into the globular cluster environment by millisecond pulsars – have long been unclear. To address this, we search for evidence of inverse Compton emission in 8-yr Fermi–LAT data from the directions of 157 Milky Way globular clusters. We find a mildly statistically significant (3.8σ) correlation between the measured globular cluster gamma-ray luminosities and their photon field energy densities. However, this may also be explained by a hidden correlation between the photon field densities and the stellar encounter rates of globular clusters. Analysed in toto, we demonstrate that the gamma-ray emission of globular clusters can be resolved spectrally into two components: (i) an exponentially cut-off power law and (ii) a pure power law. The latter component – which we uncover at a significance of 8.2σ – has a power index of 2.79 ± 0.25. It is most naturally interpreted as inverse Compton emission by cosmic-ray electrons and positrons injected by millisecond pulsars. Here, we find the luminosity of this power-law component is comparable to, or slightly smaller than, the luminosity of the curved component, suggesting the fraction of millisecond pulsar spin-down luminosity into relativistic leptons is similar to the fraction of the spin-down luminosity into prompt magnetospheric radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints.We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Orbital Modulation of Gamma Rays from PSR J2339-0533

In this work, we report on orbital modulation of the 100–600 MeV gamma-ray emission of the P B = 4.6 hr millisecond pulsar binary PSR J2339-0533 using 11 yr of Fermi Large Area Telescope data. The modulation has high significance (chance probability p ≈ 10 -7 ), is approximately sinusoidal, peaks near pulsar superior conjunction, and is detected only in the low-energy 100–600 MeV band. The modulation is confined to the on-pulse interval, suggesting that the variation is in the 2.9 ms pulsed signal itself. Additionally, this contrasts with the few other known systems exhibiting GeV orbital modulations, as these are unpulsed and generally associated with beamed emission from an intrabinary shock. The origin of the modulated pulsed signal is not yet clear, although we describe several scenarios, including Compton upscattering of photons from the heated companion. This would require high coherence in the striped pulsar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars

Abstract Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, F ν ∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Glassy quantum nuclear pasta in neutron star crusts

We conduct a comprehensive survey of the shape parameter space of the nuclear pasta phases in neutron star crusts by conducting three-dimensional Hartree-Fock + BCS calculations. Spaghetti, waffles, lasagna, bicontinuous phases and cylindrical holes occupy local minima in the resulting constant-pressure Gibbs energy surfaces, implying multiple geometries coexist at a given depth. Notably, the bicontinuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions appears over a large depth range. Our results support the idea that nuclear pasta is a glassy system. At a characteristic temperature, of order 10 8 – 10 9 K, different phases may become frozen into domains whose sizes we estimate to be 1–50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, very little long-range order exists and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature-dependent transport properties in the deep crust layers. We identify four regions distinguished by whether pasta is the true ground state, and whether the pasta structure allows delocalization of protons. The whole pasta region can occupy up to 70% of the crust by mass and 25% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 15% of its thickness.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Isospin equilibration in neutron star mergers

We analyze the isospin equilibration properties of neutrinoless nuclear (npe) matter in the temperature and density range that is relevant to neutron star mergers. Our analysis incorporates neutrino-transparency corrections to the isospin (“beta”) equilibrium condition which become noticeable at T ≳ 1 MeV. We find that the isospin relaxation rate rises rapidly as temperature rises, and at T ≈ 5 MeV it is comparable to the timescale of the density oscillations that occur immediately after the merger. This produces a resonant peak in the bulk viscosity at T ≈ 5 MeV, which causes density oscillations to be damped on the timescale of the merger. Furthermore, our calculations suggest that isospin relaxation dynamics may also be relevant when neutrinos are treated more accurately via neutrino transport schemes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗