Engineering topics

Miller, M. Coleman

Publications and source records attributed to Miller, M. Coleman.

18 records

DOE OSTI2024

A More Precise Measurement of the Radius of PSR J0740+6620 Using Updated NICER Data

PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be 2.08 ± 0.07 M ⊙ . Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. and Riley et al. reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 2020 April 17, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 2021 December 28 was presented in Salmi et al. Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 April 21, totaling ~1.1 Ms additional exposure compared to the data set analyzed in Miller et al. and Riley et al., and to data from XMM-Newton observations. We find that the equatorial circumferential radius of PSR J0740+6620 is ${12.92}_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty ~83% the width of that reported in Miller et al., in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al., then our 68% credible region would become $R={12.76}_{-1.02}^{+1.49}$ km, which is close to the headline result of Salmi et al. Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

A Search for High-Energy Counterparts to Fast Radio Bursts

We report on a search for high-energy counterparts to fast radio bursts (FRBs) with the Fermi Gamma-ray Burst Monitor, Fermi Large Area Telescope, and the Neil Gehrels Swift Observatory Burst Alert Telescope. We find no significant associations for any of the 23 FRBs in our sample, but report upper limits to the high-energy fluence for each on timescales of 0.1, 1, 10, and 100 s. We report lower limits on the ratio of the radio to high-energy fluence, f(sub r) /f(sub γ), for timescales of 0.1 and 100 s. We discuss the implications of our non-detections on various proposed progenitor models for FRBs, including analogs of giant pulses from the Crab pulsar and hyperflares from magnetars. This work demonstrates the utility of analyses of high-energy data for FRBs in tracking down the nature of these elusive sources..

Cunningham, Virginia↗

NASA NTRS2012

Pulsar Emission Geometry and Accelerating Field Strength

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems, The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models. within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle. maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius. under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry

DeCesar, Megan E.↗

NASA NTRS2012

General Relativistic Simulations of Magnetized Plasmas Around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this paper we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe, total amplification of the magnetic field of approx 2 orders of magnitude which is driven by the accretion onto the binary and that leads to stronger electromagnetic signals than in the force-free regime where such amplifications are not possible.

Giacomazzo, Bruno↗

NASA NTRS2011

PSR J0007+7303 in the CTA1 SNR: New Gamma-ray Results from Two Years of Fermi-LAT Observations

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of gamma-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its gamma-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of I-ray emission in the off-pulse phase interval at the approx. 6sigma level. The flux from this emission in the energy range E > or =::: 100 MeV is F(sub 100) = (1.73+/-0.40) x 10(exp -8) photons/sq cm/s and is best fitted by a power law with a photon index of Gamma = 2.54+/-0.14. The pulsed gamma-ray flux in the same energy range is F(sub 100) = (3.95+/-0.07) x 10(exp -7) photons/sq cm/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of Gamma = 1.41+/-0.23 and a cutoff energy E(sub c) = 4.04+/-0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the gamma-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

Abdo, A. A.↗

NASA NTRS2011

Gamma-Ray Pulsar Light Curves in Vacuum and Force-Free Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

Harding, Alice K.↗

NASA NTRS2011

Gamma-Ray Pulsar Light Curves in Offset Polar Cap Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres, used to model high-energy light curves have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profile~. We find that. corn pared to the profile:-; derived from :-;ymmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines. increases significantly relative to the off-peak emission. formed along leading field lines. The enhanced contrast produces greatly improved slot gap model fits to Fermi pulsar light curves like Vela, which show very little off-peak emIssIon.

Harding, Alice K.↗

NASA NTRS2010

Test of a General Formula for Black Hole Gravitational Wave Kicks

Although the gravitational wave kick velocity in the orbital plane of coalescing black holes has been understood for some time, apparently conflicting formulae have been proposed for the dominant outof- plane kick, each a good fit to different data sets. This is important to resolve because it is only the out-of-plane kicks that can reach more than 500 km s-l and can thus eject merged remnants from galaxies. Using a different ansatz for the out-of-plane kick, we show that we can fit almost all existing data to better than 5%. This is good enough for any astrophysical calculation and shows that the previous apparent conflict was only because the two data sets explored different aspects of the kick parameter space.

van Meter, James R.↗

NASA NTRS2009

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 10(exp 4) mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approximately equal to 5-30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology. This work is supported in part by the Cooperative Education Program at NASA/GSFC.

Alicea-Munoz, E.↗

NASA NTRS2008

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 1 0A4 mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approx. 5 - 30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology.

Alicea-Munoz, E.↗

NASA NTRS2008

Black Hole Mergers as Probes of Structure Formation

Intense structure formation and reionization occur at high redshift, yet there is currently little observational information about this very important epoch. Observations of gravitational waves from massive black hole (MBH) mergers can provide us with important clues about the formation of structures in the early universe. Past efforts have been limited to calculating merger rates using different models in which many assumptions are made about the specific values of physical parameters of the mergers, resulting in merger rate estimates that span a very wide range (0.1 - 104 mergers/year). Here we develop a semi-analytical, phenomenological model of MBH mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z 5 - 30 structure formation era. We do this by generating synthetic LISA observable data (total BH mass, BH mass ratio, redshift, merger rates), which are then analyzed using a Markov Chain Monte Carlo method. This allows us to constrain the physical parameters of the mergers. We find that our methodology works well at estimating merger parameters, consistently giving results within 1- of the input parameter values. We also discover that the number of merger events is a key discriminant among models. This helps our method be robust against observational uncertainties. Our approach, which at this stage constitutes a proof of principle, can be readily extended to physical models and to more general problems in cosmology and gravitational wave astrophysics.

Alicea-Munoz, E.↗

NASA NTRS2008

Modeling Kicks from the Merger of Generic Black-hole Binaries

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum "recoil kick" of up to approx. 4000 km/s. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick parallel to the orbital axis does not scale as approx. eta(sup 2) (where eta is the symmetric mass ratio), as previously proposed, but is more consistent with approx. eta(sup 3). We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters. S

Baker, John G.↗

NASA NTRS2005

QPO Constraints on Neutron Stars

The kilohertz frequencies of QPOs from accreting neutron star systems imply that they are generated in regions of strong gravity, close to the star. This suggests that observations of the QPOs can be used to constrain the properties of neutron stars themselves, and in particular to inform us about the properties of cold matter beyond nuclear densities. Here we discuss some relatively model-insensitive constraints that emerge from the kilohertz QPOs, as well as recent developments that may hint at phenomena related to unstable circular orbits outside neutron stars.

Miller, M. Coleman↗

NASA NTRS2005

Getting a Kick Out of Numerical Relativity

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation; and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. Our estimated kick is 10(exp 5) km/s with an error of less than 10%. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts z greater than or approx. equal to 10, halos with masses less than or approx. equal to 10(exp 9) Solar Mass will have difficulty retaining coalesced black holes after major mergers.

Baker, John G.↗

NASA NTRS1995

Phase lags in Cygnus X-1

Coronae of hot electrons are ubiquitous features in models of the inner regions of accreting black holes and nonmagnetic neutron stars. The scattering optical depth of these coronae inferred from observations is often tau approximately 3, so the energy spectrum of the disk in these regions is likely masked by the spectrum of upscattered photons in the corona and some of the disk properties are thus obscured. Observations of the dependence on not available with the energy spectrum. In our picture, the disk emission is modeled as a photon source injected into a Comptonizing corona; Comptonization may also occur in the disk, but for our purposes the disk emission is simply an input to the corona. We show that, contrary to some claims, the functional dependence of lag on Fourier frequency emerges intact from transit through the corona, modulo a multiplicative factor (which may in principle be negative, so that a phase lag can be changed to a phase lead), even if the properties of the corona vary with time. We also show that any frequency dependence of the lag due to variation in the corona itself is only second order in the amplitude of the variation, and cannot exceed the transit time approximately ms of the corona; thus, the lags of up to 0.1 s seen in several black hole candidates come from lags in the emission from the disk. Finally, we predict that plots of the time lag versus Fourier frequency in black hole candidates should have a 'shelf' of constant lag equal to the coronal lag (approximately 1 ms), with the constant being proportional to In (E(sub 2)/E(sub 1)) for the lag between energies E(sub 2) and E(sub 1). The lack of such a shelf in current observations of several galactic black hole candidates constraints the radii of the coronae to be R less than or approximately = 10(exp 8) tau cm.

Miller, M. Coleman↗

NASA NTRS1995

Critical radiation fluxes and luminosities of black holes and relativistic stars

The critial luminosity at which the outward force of radiation balances the inward force of gravity plays an important role in many astrophysical systems. We present expressions for the radiation force on particles with arbitrary cross sections and analyze the radiation field produced by radiating matter, such as a disk, ring, boundary layer, or stellar surface, that rotates slowly around a slowly rotating gravitating mass. We then use these results to investigate the critical radiation flux and, where possible, the critical luminosity of such a system in genral relativity. We demonstrate that if the radiation source is axisymmetric and emission is back-front symmetric with repect to the local direction of motion of the radiating matter, as seen in the comoving frame, then the radial component of the radiation flux and the diagonal components of the radiation stress-energy tensor outside the source are the same, to first order in the rotation rates, as they would be if the radiation source and gravitating mass were not rotating. We argue that the critical radiation flux for matter at rest in the locally nonrotating frame is often satisfactory as an astrophysical benchmark flux and show that if this benchmark is adopted, many of the complications potentially introduced by rotation of the radiation source and the gravitating mass are avoided. We show that if the radiation field in the absence of rotation would be spherically symmetric and the opacity is independent of frequency and direction, one can define a critical luminosity for the system that is independent of frequency and direction, one can define a critical luminosity for the system that is independent of the spectrum and angular size of the radiation source and is unaffected by rotation of the source and mass and orbital motion of the matter, to first order. Finally, we analyze the conditions under which the maximum possible luminosity of a star or black hole powered by steady spherically symmetric radial accretion is the same in general relativity as in the Newtonian limit.

Lamb, Frederick K.↗

NASA NTRS1994

Electrodynamics of disk-accreting magnetic neutron stars

We have investigated the electrodynamics of magnetic neutron stars accreting from Keplerian disks and the implications for particle acceleration and gamma-ray emission by such systems. We argue that the particle density in the magnetospheres of such stars is larger by orders of magnitude than the Goldreich-Julian density, so that the formation of vacuum gaps is unlikely. We show that even if the star rotates slowly, electromotive forces (EMFs) of order 10(exp 15) V are produced by the interaction of plasma in the accretion disk with the magnetic field of the neutron star. The resistance of the disk-magnetosphere-star circuit is small, and hence these EMFs drive very large conduction currents. Such large currents are likely to produce magnetospheric instabilities, such as relativistic double layers and reconnection events, that can accelerate electrons or ions to very high energies.

Miller, M. Coleman↗

NASA NTRS1994

Disk-accreting magnetic neutron stars as high-energy particle accelerators

Interaction of an accretion disk with the magnetic field of a neutron star produces large electromotive forces, which drive large conduction currents in the disk-magnetosphere-star circuit. Here we argue that such large conduction currents will cause microscopic and macroscopic instabilities in the magnetosphere. If the minimum plasma density in the magnetosphere is relatively low is less than or aproximately 10(exp 9)/cu cm, current-driven micro-instabilities may cause relativistic double layers to form, producing voltage differences in excess of 10(exp 12) V and accelerating charged particles to very high energies. If instead the plasma density is higher (is greater than or approximately = 10(exp 9)/cu cm, twisting of the stellar magnetic field is likely to cause magnetic field reconnection. This reconnection will be relativistic, accelerating plasma in the magnetosphere to relativistic speeds and a small fraction of particles to very high energies. Interaction of these high-energy particles with X-rays, gamma-rays, and accreting plasma may produce detectable high-energy radiation.

Hamilton, Russell J.↗