BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Coleman Miller

M. Coleman Miller appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-phastro-ph.COastro-ph.GAnucl-thastro-ph.IMastro-ph.SRhep-thnucl-ex

Catalog footprint

What is connected

70works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

70 published item(s)

preprint2023arXiv

Kilohertz quasiperiodic oscillations in short gamma-ray bursts

Short gamma-ray bursts are associated with binary neutron star mergers, which are multimessenger astronomical events that have been observed both in gravitational waves and in the multiband electromagnetic spectrum. Depending on the masses of the stars in the binary and on details of their largely unknown equation of state, a dynamically evolving and short-lived neutron star may be formed after the merger, existing for approximately 10-300 ms before collapsing to a black hole. Numerical relativity simulations across different groups consistently show broad power spectral features in the 1-5 kHz range in the post-merger gravitational wave signal, which is inaccessible by current gravitational-wave detectors but could be seen by future third generation ground-based detectors in the next decade. This implies the possibility of quasiperiodic modulation of the emitted gamma-rays in a subset of events where a neutron star is formed shortly prior to the final collapse to a black hole. Here we present two such signals identified in the short bursts GRB 910711 and GRB 931101B from archival BATSE data, which are compatible with the predictions from numerical relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The Radius of PSR J0740+6620 from NICER with NICER Background Estimates

We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electromagnetic Counterparts to Massive Black Hole Mergers

The next two decades are expected to open the door to the first coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures associated with massive black hole (MBH) binaries heading for coalescence. These detections will launch a new era of multimessenger astrophysics by expanding this growing field to the low-frequency GW regime and will provide an unprecedented understanding of the evolution of MBHs and galaxies. They will also constitute fundamentally new probes of cosmology and would enable unique tests of gravity. The aim of this Living Review is to provide an introduction to this research topic by presenting a summary of key findings, physical processes and ideas pertaining to EM counterparts to MBH mergers as they are known at the time of this writing. We review current observational evidence for close MBH binaries, discuss relevant physical processes and timescales, and summarize the possible EM counterparts to GWs in the precursor, coalescence, and afterglow stages of a MBH merger. We also describe open questions and discuss future prospects in this dynamic and quick-paced research area.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational-Wave and X-ray Probes of the Neutron Star Equation of State

Neutron stars are a remarkable marriage of Einstein's theory of general relativity with nuclear physics. Their interiors harbor extreme matter that cannot be probed in the laboratory. At such high densities and pressures, their cores may consist predominantly of exotic matter such as free quarks or hyperons. Gravitational wave observations from the Laser Interferometer Gravitational-wave Observatory (LIGO) and from other interferometers, and X-ray observations from the Neutron Star Interior Composition Explorer (NICER), are beginning to pierce through the veil. These observations provide information about neutron star cores, and therefore, about the physics that makes such objects possible. In this review, we discuss what we have learned about the physics of neutron stars from gravitational wave and X-ray observations. We focus on what has been observed with certainty and what should be observable in the near future, with an eye out for the physics that these new observations will teach us.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Observing intermediate-mass black holes and the upper--stellar-mass gap with LIGO and Virgo

Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses $50\mbox{--} 500\, M_\odot$, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at $2σ$, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of $\sim 10-40\%$ at a signal-to-noise ratio of $20$. Focusing on the stellar-mass gap with new tabulations of the $^{12}\text{C}(α, γ)^{16} \text{O}$ reaction rate and its uncertanties, we evolve massive helium core stars using \MESA\, to establish the lower and upper edge of the mass gap as $\simeq$\,59$^{+34}_{-13}$\,$M_{\odot}$ and $\simeq$\,139$^{+30}_{-14}$\,$M_{\odot}$ respectively, where the error bars give the mass range that follows from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16} \text{O}$ nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we re-analyse GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90\% credibility.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Uncertain Future of Massive Binaries Obscures the Origin of LIGO/Virgo Sources

The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial black holes. Here we demonstrate that such conclusions cannot yet be robust, because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss and supernova outbursts. Using standard tools such as the population synthesis codes StarTrack and COMPAS and the detailed stellar evolution code MESA, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system. We show that, despite its well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BH-BH binary (which might be seen in microlensing observations), or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Investigating the I-Love-Q and w-mode Universal Relations Using Piecewise Polytropes

Neutron stars are expected to have a tight relation between their moment of inertia ($I$), tidal deformability ($λ$, which is related to the Love number), and rotational mass quadrupole moment ($Q$) that is nearly independent of the unknown equation of state (EoS) of cold dense matter. These and similar relations are often called "universal", and they have been used for various applications including analysis of gravitational wave data. We extend these studies using piecewise polytropic representations of dense matter, including for so-called twin stars that have a second branch of stability at high central densities. The second-branch relations are less tight, by a factor of $\sim 3$, than the relations found in the first stable branch. We find that the relations on both branches become tighter when we increase the lower limit to the maximum mass for the EoS under consideration. We also propose new empirical relations between $I$, $λ$, $Q$, and the complex frequency $ω=ω_R+iω_I$ of the fundamental axial $w$-mode, and find that they are comparably tight to the I-Love-Q correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Merger rates in primordial black hole clusters without initial binaries

Primordial black holes formed through the collapse of cosmological density fluctuations have been hypothesised as contributors to the dark matter content of the Universe. At the same time, their mergers could contribute to the recently observed population of gravitational-wave sources. We investigate the scenario in which primordial black holes form binaries at late times in the Universe. Specifically, we re-examine the mergers of primordial black holes in small clusters of ~30 objects in the absence of initial binaries. Binaries form dynamically through Newtonian gravitational interactions. These binaries act as heat sources for the cluster, increasing the cluster's velocity dispersion, which inhibits direct mergers through gravitational-wave two-body captures. Meanwhile, three-body encounters of tight binaries are too rare to tighten binaries sufficiently to allow them to merge through gravitational-wave emission. We conclude that in the absence of initial binaries, merger rates of primordial black holes in the Bird et al. (2016) initial cluster configuration are at least an order of magnitude lower than previously suggested, which makes gravitational-wave detections of such sources improbable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Constraining the equation of state of high-density cold matter using nuclear and astronomical measurements

The increasing richness of data related to cold dense matter, from laboratory experiments to neutron-star observations, requires a framework for constraining the properties of such matter that makes use of all relevant information. Here, we present a rigorous but practical Bayesian approach that can include diverse evidence, such as nuclear data and the inferred masses, radii, tidal deformabilities, moments of inertia, and gravitational binding energies of neutron stars. We emphasize that the full posterior probability distributions of measurements should be used rather than, as is common, imposing a cut on the maximum mass or other quantities. Our method can be used with any parameterization of the equation of state (EOS). We use both a spectral parameterization and a piecewise polytropic parameterization with variable transition densities to illustrate the implications of current measurements and show how future measurements in many domains could improve our understanding of cold catalyzed matter. We find that different types of measurements will play distinct roles in constraining the EOS in different density ranges. For example, better symmetry energy measurements will have a major influence on our understanding of matter somewhat below nuclear saturation density but little influence above that density. In contrast, precise radius measurements or multiple tidal deformability measurements of the quality of those from GW170817 or better will improve our knowledge of the EOS over a broader density range.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Star Formation in Accretion Disks and SMBH Growth

Accretion disks around active galactic nuclei are potentially unstable to star formation at large radii. We note that when the compact objects formed from some of these stars spiral into the central supermassive black hole, there is no radiative feedback and therefore the accretion rate is not limited by radiation forces. Using a set of accretion disk models, we calculate the accretion rate onto the central supermassive black hole in both gas and compact objects. We find that the timescale for a supermassive black hole to double in mass can decrease by factors ranging from $\sim0.7$ to as low as $\sim0.1$ in extreme cases, compared to gas accretion alone. Our results suggest that the formation of extremely massive black holes at high redshift may occur without prolonged super-Eddington gas accretion or very massive seed black holes. We comment on potential observational signatures as well as implications for other observations of active galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A New Method for Finding Point Sources in High-energy Neutrino Data

The IceCube collaboration has reported the first detection of high-energy astrophysical neutrinos including $\sim 50$ high-energy starting events, but no individual sources have been identified. It is therefore important to develop the most sensitive and efficient possible algorithms to identify point sources of these neutrinos. The most popular current method works by exploring a dense grid of possible directions to individual sources, and identifying the single direction with the maximum probability of having produced multiple detected neutrinos. This method has numerous strengths, but it is computationally intensive and, because it focuses on the single best location for a point source, additional point sources are not included in the evidence. We propose a new maximum likelihood method that uses the angular separations between all pairs of neutrinos in the data. Unlike existing autocorrelation methods for this type of analysis, which also use angular separations between neutrino pairs, our method incorporates information about the point spread function and can identify individual point sources. We find that if the angular resolution is a few degrees or better, then this approach reduces both false positive and false negative errors compared to the current method, and is also more computationally efficient up to, potentially, hundreds of thousands of detected neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A test of the nature of the Fe K Line in the neutron star low-mass X-ray binary Serpens X-1

Broad Fe K emission lines have been widely observed in the X-ray spectra of black hole systems, and in neutron star systems as well. The intrinsically narrow Fe K fluorescent line is generally believed to be part of the reflection spectrum originating in an illuminated accretion disk, and broadened by strong relativistic effects. However, the nature of the lines in neutron star LMXBs has been under debate. We therefore obtained the longest, high-resolution X-ray spectrum of a neutron star LMXB to date with a 300 ks Chandra HETGS observation of Serpens X-1. The observation was taken under the "continuous clocking" mode and thus free of photon pile-up effects. We carry out a systematic analysis and find that the blurred reflection model fits the Fe line of Serpens X-1 significantly better than a broad Gaussian component does, implying that the relativistic reflection scenario is much preferred. Chandra HETGS also provides highest spectral resolution view of the Fe K region and we find no strong evidence for additional narrow lines.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Implications of the Gravitational Wave Event GW150914

The era of gravitational-wave astronomy began on 14 September 2015, when the LIGO Scientific Collaboration detected the merger of two $\sim 30 M_\odot$ black holes at a distance of $\sim 400$ Mpc. This event has facilitated qualitatively new tests of gravitational theories, and has also produced exciting information about the astrophysical origin of black hole binaries. In this review we discuss the implications of this event for gravitational physics and astrophysics, as well as the expectations for future detections. In brief: (1) because the spins of the black holes could not be measured accurately and because mergers are not well calculated for modified theories of gravity, the current analysis of GW150914 does not place strong constraints on gravity variants that change only the generation of gravitational waves, but (2) it does strongly constrain alterations of the propagation of gravitational waves and alternatives to black holes. Finally, (3) many astrophysical models for the origin of heavy black hole binaries such as the GW150914 system are in play, but a reasonably robust conclusion that was reached even prior to the detection is that the environment of such systems needs to have a relatively low abundance of elements heavier than helium.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Observational Constraints on Neutron Star Masses and Radii

Precise and reliable measurements of the masses and radii of neutron stars with a variety of masses would provide valuable guidance for improving models of the properties of cold matter with densities above the saturation density of nuclear matter. Several different approaches for measuring the masses and radii of neutron stars have been tried or proposed, including analyzing the X-ray fluxes and spectra of the emission from neutron stars in quiescent low-mass X-ray binary systems and thermonuclear burst sources; fitting the energy-dependent X-ray waveforms of rotation-powered millisecond pulsars, burst oscillations with millisecond periods, and accretion-powered millisecond pulsars; and modeling the gravitational radiation waveforms of coalescing double neutron star and neutron star -- black hole binary systems. We describe the strengths and weaknesses of these approaches, most of which currently have substantial systematic errors, and discuss the prospects for decreasing the systematic errors in each method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Case for PSR J1614-2230 as a NICER Target

The Neutron star Interior Composition ExploreR (NICER) is expected to launch in early 2017 and will gather X-ray data on neutron stars and other high-energy sources from a berth on the International Space Station. Its prime scientific goal is to measure the masses and radii of non-accreting neutron stars via fits to the energy-dependent waveforms produced by the rotation of hot spots on their surfaces. These measurements will provide valuable input to theoretical models of cold matter beyond nuclear density. Here we propose that PSR J1614$-$2230, despite its low count rate, is a promising source to observe with NICER. The reason is that XMM-Newton observations suggest that the fractional oscillation amplitude from PSR J1614$-$2230 could be high enough that this star cannot be very compact. We show that if we analyze 0.5 Ms of NICER data and 0.1 Ms of nearby off-source data and combine that analysis with the known mass of this star, we would find a robust lower limit to the radius with a statistical uncertainty of only $\sim 0.5-0.7$ km. We also show that even if there is an unmodeled nonthermal component modulated at the pulsation frequency, good statistical fits could rule out significant biases. The low count rate will make reliable upper limits on the radius difficult, but the lower limit could rule out some equations of state that are currently being discussed. This analysis would require a good estimate of the non-source background, so Chandra observations of the vicinity of PSR J1614$-$2230 would be helpful.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Role of the Kozai-Lidov Mechanism in Black Hole Binary Mergers in Galactic Centers

In order to understand the rate of merger of stellar-mass black hole binaries (BHBs) by gravitational wave (GW) emission it is important to determine the major pathways to merger. We use numerical simulations to explore the evolution of BHBs inside the radius of influence of supermassive black holes (SMBHs) in galactic centers. In this region the evolution of binaries is dominated by perturbations from the central SMBH. In particular, as first pointed out by Antonini and Perets, the Kozai-Lidov (KL) mechanism trades relative inclination of the BHB to the SMBH for eccentricity of the BHB, and for some orientations can bring the BHB to an eccentricity near unity. At very high eccentricities, GW emission from the BHB can become efficient, causing the members of the BHB to coalesce. We use a novel combination of two N-body codes to follow this evolution. We are forced to simulate small systems to follow the behavior accurately. We have completed 400 simulations that range from $\sim$ 300 stars around a $10^{3}$ M$_{\odot}$ black hole to $\sim$ 4500 stars around a $10^{4}$ M$_{\odot}$ black hole. These simulations are the first to follow the internal orbit of a binary near a SMBH while also following the changes to its external orbit self-consistently. We find that this mechanism could produce mergers at a maximum rate per volume of $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ or considerably less if the inclination oscillations of the binary remain constant as the BHB inclination to the SMBH changes, or if the binary black hole fraction is small.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Upper Bound on Neutron Star Masses from Models of Short Gamma-ray Bursts

The discovery of two neutron stars with gravitational masses $\approx 2~M_\odot$ has placed a strong lower limit on the maximum mass of nonrotating neutron stars, and with it a strong constraint on the properties of cold matter beyond nuclear density. Current upper mass limits are much looser. Here we note that, if most short gamma-ray bursts are produced by the coalescence of two neutron stars, and if the merger remnant collapses quickly, then the upper mass limit is constrained tightly. If the rotation of the merger remnant is limited only by mass-shedding (which seems probable based on numerical studies), then the maximum gravitational mass of a nonrotating neutron star is $\approx 2-2.2~M_\odot$ if the masses of neutron stars that coalesce to produce gamma-ray bursts are in the range seen in Galactic double neutron star systems. These limits would be increased by $\sim 4$% in the probably unrealistic case that the remnants rotate at $\sim 30$% below mass-shedding, and by $\sim 15$% in the extreme case that the remnants do not rotate at all. Future coincident detection of short gamma-ray bursts with gravitational waves will strengthen these arguments because they will produce tight bounds on the masses of the components for individual events. If these limits are accurate then a reasonable fraction of double neutron star mergers might not produce gamma-ray bursts. In that case, or in the case that many short bursts are produced instead by the mergers of neutron stars with black holes, the implied rate of gravitational wave detections will be increased.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Determining Neutron Star Properties by Fitting Oblate-Star Waveforms To X-ray Burst Oscillations

We describe sophisticated new Bayesian analysis methods that make it possible to estimate quickly the masses and radii of rapidly rotating, oblate neutron stars by fitting oblate-star waveform models to energy-resolved observations of the X-ray oscillations produced by a hot spot on such stars. We conclude that models that take the oblate shape of the star into account should be used for stars with large radii and rotation rates $>300$ Hz. We find that a 25% variation of the temperature of the hot spot with latitude does not significantly bias estimates of the mass $M$ and equatorial radius $R_{\rm eq}$ derived by fitting a model that assumes a uniform-temperature spot. Our results show that fits of oblate-star waveform models to waveform data can simultaneously determine $M$ and $R_{\rm eq}$ with uncertainties $\lesssim\,$7% if (1) the star's rotation rate is $\gtrsim\,$$600$ Hz; (2) the spot center and observer's sightline are both within $30^\circ$ of the star's rotational equator; (3) the oscillations have a fractional rms amplitude $\gtrsim\,$10%; and (4) $\gtrsim$$10^7$ counts are collected from the star. This is a realistic fractional amplitude, and this many counts could be obtained from a single star by the accepted NICER and proposed LOFT and AXTAR space missions by combining data from many X-ray bursts. These uncertainties are small enough to improve substantially our understanding of cold, ultradense matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disk Winds as an Explanation for Slowly Evolving Temperatures in Tidal Disruption Events

Among the many intriguing aspects of optically discovered tidal disruption events (TDEs) is that their temperatures are lower than expected and that the temperature does not evolve as rapidly with decreasing fallback rate as would be expected in standard disk theory. We show that this can be explained qualitatively using an idea proposed by Laor & Davis in the context of normal active galactic nuclei: that larger accretion rates imply stronger winds and thus that the accretion rate through the inner disk only depends weakly on the inflow rate at the outer edge of the disk. We also show that reasonable quantitative agreement with data requires that, as has been suggested in recent papers, the circularization radius of the tidal stream is approximately equal to the semimajor axis of the most bound orbit of the debris rather than twice the pericenter distance as would be expected without rapid angular momentum redistribution. If this explanation is correct, it suggests that the evolution of TDEs may test both non-standard disk theory and the details of the interactions of the tidal stream.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Flows of X-ray gas reveal the disruption of a star by a massive black hole

Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Soft X-ray Temperature Tidal Disruption Events from Stars on Deep Plunging Orbits

One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$~K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$~K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the supermassive black hole, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $β$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to higher temperatures of more compact disks which are super-Eddington and geometrically and optically thick. For a given star and $β$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at the Eddington luminosity. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with $β> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both expectations are consistent with current data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Wind Accretion Model for HLX-1

The brightest ultraluminous X-ray source currently known, HLX-1, has been observed to undergo five outburst cycles. The periodicity of these outbursts, and their high inferred maximum accretion rates of $\sim{\rm few}\times 10^{-4} M_\odot {\rm yr}^{-1}$, naturally suggest Roche lobe overflow at the pericenter of an eccentric orbit. It is, however, difficult for the Roche lobe overflow model to explain the apparent trend of decreasing decay times over the different outbursts while the integrated luminosity also drops. Thus if the trend is real rather than simply being a reflection of the complex physics of accretion disks, a different scenario may be necessary. We present a speculative model in which, within the last decade, a high-mass giant star had most of its envelope tidally stripped by the $\sim 10^{4-5} M_\odot$ black hole in HLX-1, and the remaining core plus low-mass hydrogen envelope now feeds the hole with a strong wind. This model can explain the short decay time of the disk, and could explain the fast decrease in decay time if the wind speed increases with time. A key prediction of this model is that there will be excess line absorption due to the wind; our analysis does in fact find a flux deficit in the $\sim 0.9-1.1$ keV range that is consistent with predictions, albeit at low significance. If this idea is correct, we also expect that within tens of years the bound material from the original disruption will return and will make HLX-1 a persistently bright source.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Observational Signatures of Binary Supermassive Black Holes

Observations indicate that most massive galaxies contain a supermassive black hole, and theoretical studies suggest that when such galaxies have a major merger, the central black holes will form a binary and eventually coalesce. Here we discuss two spectral signatures of such binaries that may help distinguish them from ordinary AGN. These signatures are expected when the mass ratio between the holes is not extreme and the system is fed by a circumbinary disk. One such signature is a notch in the thermal continuum that has been predicted by other authors; we point out that it should be accompanied by a spectral revival at shorter wavelengths and also discuss its dependence on binary properties such as mass, mass ratio, and separation. In particular, we note that the wavelength $λ_n$ at which the notch occurs depends on these three parameters in such a way as to make the number of systems displaying these notches $\propto λ_n^{16/3}$; longer wavelength searches are therefore strongly favored. A second signature, first discussed here, is hard X-ray emission with a Wien-like spectrum at a characteristic temperature $\sim 100$ keV produced by Compton cooling of the shock generated when streams from the circumbinary disk hit the accretion disks around the individual black holes. We investigate the observability of both signatures. The hard X-ray signal may be particularly valuable as it can provide an indicator of black hole merger a few decades in advance of the event.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sowing the seeds of massive black holes in small galaxies: Young clusters as the building blocks of Ultra-Compact-Dwarf Galaxies

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Masses and Spins of Neutron Stars and Stellar-Mass Black Holes

Stellar-mass black holes and neutron stars represent extremes in gravity, density, and magnetic fields. They therefore serve as key objects in the study of multiple frontiers of physics. In addition, their origin (mainly in core-collapse supernovae) and evolution (via accretion or, for neutron stars, magnetic spindown and reconfiguration) touch upon multiple open issues in astrophysics. In this review, we discuss current mass and spin measurements and their reliability for neutron stars and stellar-mass black holes, as well as the overall importance of spins and masses for compact object astrophysics. Current masses are obtained primarily through electromagnetic observations of binaries, although future microlensing observations promise to enhance our understanding substantially. The spins of neutron stars are straightforward to measure for pulsars, but the birth spins of neutron stars are more difficult to determine. In contrast, even the current spins of stellar-mass black holes are challenging to measure. As we discuss, major inroads have been made in black hole spin estimates via analysis of iron lines and continuum emission, with reasonable agreement when both types of estimate are possible for individual objects, and future X-ray polarization measurements may provide additional independent information. We conclude by exploring the exciting prospects for mass and spin measurements from future gravitational wave detections, which are expected to revolutionize our understanding of strong gravity and compact objects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Alignment of supermassive black hole binary orbits and spins

Recent studies of accretion onto supermassive black hole binaries suggest that much, perhaps most, of the matter eventually accretes onto one hole or the other. If so, then for binaries whose inspiral from ~1 pc to 0.001 - 0.01 pc is driven by interaction with external gas, both the binary orbital axis and the individual black hole spins can be reoriented by angular momentum exchange with this gas. Here we show that, unless the binary mass ratio is far from unity, the spins of the individual holes align with the binary orbital axis in a time few-100 times shorter than the binary orbital axis aligns with the angular momentum direction of the incoming circumbinary gas; the spin of the secondary aligns more rapidly than that of the primary by a factor ~(m_1/m_2)^{1/2}>1. Thus the binary acts as a stabilizing agent, so that for gas-driven systems, the black hole spins are highly likely to be aligned (or counteraligned if retrograde accretion is common) with each other and with the binary orbital axis. This alignment can significantly reduce the recoil speed resulting from subsequent black hole merger.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Astrophysical Constraints on Dense Matter in Neutron Stars

Ever since the discovery of neutron stars it has been realized that they serve as probes of a physical regime that cannot be accessed in laboratories: strongly degenerate matter at several times nuclear saturation density. Existing nuclear theories diverge widely in their predictions about such matter. It could be that the matter is primarily nucleons, but it is also possible that exotic species such as hyperons, free quarks, condensates, or strange matter may dominate this regime. Astronomical observations of cold high-density matter are necessarily indirect, which means that we must rely on measurements of quantities such as the masses and radii of neutron stars and their surface effective temperatures as a function of age. Here we review the current status of constraints from various methods and the prospects for future improvements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Determining neutron star masses and radii using energy-resolved waveforms of X-ray burst oscillations

Simultaneous, precise measurements of the mass $M$ and radius $R$ of neutron stars can yield uniquely valuable information about the still uncertain properties of cold matter at several times the density of nuclear matter. One method that could be used to measure $M$ and $R$ is to analyze the energy-dependent waveforms of the X-ray flux oscillations seen during some thermonuclear bursts from some neutron stars. These oscillations are thought to be produced by X-ray emission from hotter regions on the surface of the star that are rotating at or near the spin frequency of the star. Here we explore how well $M$ and $R$ could be determined by generating, and analyzing using Bayesian techniques, synthetic energy-resolved X-ray data that we produce assuming a future space mission having 2--30 keV energy coverage and an effective area of 10 m$^2$, such as the proposed \textit{LOFT} or \textit{AXTAR} missions. We find that if the hot spot is within 10$^\circ$ of the rotation equator, both $M$ and $R$ can usually be determined with an uncertainty of about 10% if there are $10^6$ total counts from the spot, whereas waveforms from spots within 20$^\circ$ of the rotation pole provide no useful constraints. These constraints can usually be achieved even if the burst oscillations vary with time and data from multiple bursts must be used to obtain 10$^6$ counts from the hot spot. This is therefore a promising method to constrain $M$ and $R$ tightly enough to discriminate strongly between competing models of cold, high-density matter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low-Frequency Terrestrial Gravitational-Wave Detectors

Direct detection of gravitational radiation in the audio band is being pursued with a network of kilometer-scale interferometers (LIGO, Virgo, KAGRA). Several space missions (LISA, DECIGO, BBO) have been proposed to search for sub-Hz radiation from massive astrophysical sources. Here we examine the potential sensitivity of three ground-based detector concepts aimed at radiation in the 0.1 -- 10\,Hz band. We describe the plethora of potential astrophysical sources in this band and make estimates for their event rates and thereby, the sensitivity requirements for these detectors. The scientific payoff from measuring astrophysical gravitational waves in this frequency band is great. Although we find no fundamental limits to the detector sensitivity in this band, the remaining technical limits will be extremely challenging to overcome.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing the rotating hot spot model using X-ray burst oscillations from 4U 1636-536

Precise and accurate measurements of neutron star masses and radii would provide valuable information about the still uncertain properties of cold matter at supranuclear densities. One promising approach to making such measurements involves analysis of the X-ray flux oscillations often seen during thermonuclear (type 1) X-ray bursts. These oscillations are almost certainly produced by emission from hotter regions on the stellar surface modulated by the rotation of the star. One consequence of the rotation is that the oscillation should appear earlier at higher photon energies than at lower energies. Ford (1999) found compelling evidence for such a hard lead in the tail oscillations of one type 1 burst from Aql X-1. We have therefore analyzed individually the oscillations observed in the tails of the four type 1 bursts from 4U 1636-536 that, when averaged, provided the strongest evidence for a soft lead in the analysis by Muno et al. (2003). We have also analyzed the oscillation observed during the superburst from this star. We find that the data from these five bursts, treated both individually and jointly, are fully consistent with a rotating hot spot model. Unfortunately, the uncertainties in these data are too large to provide interesting constraints on the mass and radius of this star.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An upper limit to the velocity dispersion of relaxed stellar systems without massive black holes

Massive black holes have been discovered in all closely examined galaxies with high velocity dispersion. The case is not as clear for lower-dispersion systems such as low-mass galaxies and globular clusters. Here we suggest that above a critical velocity dispersion of roughly 40 km/s, massive central black holes will form in relaxed stellar systems at any cosmic epoch. This is because above this dispersion primordial binaries cannot support the system against deep core collapse. If, as previous simulations show, the black holes formed in the cluster settle to produce a dense subcluster, then given the extremely high densities reached during core collapse the holes will merge with each other. For low velocity dispersions and hence low cluster escape speeds, mergers will typically kick out all or all but one of the holes due to three-body kicks or the asymmetric emission of gravitational radiation. If one hole remains, it will tidally disrupt stars at a high rate. If none remain, one is formed after runaway collisions between stars, then it tidally disrupts stars at a high rate. The accretion rate after disruption is many orders of magnitude above Eddington. If, as several studies suggest, the hole can accept matter at that rate because the generated radiation is trapped and advected, then it will grow quickly and form a massive central black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 Letter, 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 a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Be/X-ray Binary Swift J1626.6-5156 as a Variable Cyclotron Line Source

Swift J1626.6-5156 is a Be/X-ray binary that was in outburst from December 2005 until November 2008. We have examined RXTE/PCA and HEXTE spectra of three long observations of this source taken early in its outburst, when the PCA 2-20 keV count rate was >70 counts/s/PCU, as well as several combined observations from different stages of the outburst. The spectra are best fit with an absorbed cutoff power law with a ~6.4 keV iron emission line and a Gaussian optical depth absorption line at ~10 keV. We present strong evidence that this absorption-like feature is a cyclotron resonance scattering feature, making Swift J1626.6-5156 a new candidate cyclotron line source. The redshifted energy of ~10 keV implies a magnetic field strength of ~8.6(1+z) x 10^11 G in the region of the accretion column close to the magnetic poles where the cyclotron line is produced. Analysis of phase averaged spectra spanning the duration of the outburst suggests a possible positive correlation between the fundamental cyclotron energy and source luminosity. Phase resolved spectroscopy from a long observation reveals a variable cyclotron line energy, with phase dependence similar to a variety of other pulsars, as well as the first harmonic of the fundamental cyclotron line.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal disruptions of separated binaries in galactic nuclei

Several galaxies have exhibited X-ray flares that are consistent with the tidal disruption of a star by a central supermassive black hole. In theoretical treatments of this process it is usually assumed that the star was initially on a nearly parabolic orbit relative to the black hole. Such an assumption leads in the simplest approximation to a $t^{-5/3}$ decay of the bolometric luminosity and this is indeed consistent with the relatively poorly sampled light curves of such flares. We point out that there is another regime in which the decay would be different: if a binary is tidally separated and the star that remains close to the hole is eventually tidally disrupted from a moderate eccentricity orbit, the decay is slower, typically $\sim t^{-1.2}$. As a result, careful sampling of the light curves of such flares could distinguish between these processes and yield insight into the dynamics of binaries as well as single stars in galactic centres. We explore this process using three-body simulations and analytic treatments and discuss the consequences for present-day X-ray detections and future gravitational wave observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implications of high-precision spectra of thermonuclear X-ray bursts for determining neutron star masses and radii

X-ray burst spectra have long been used to estimate neutron star masses and radii. These estimates assumed that burst spectra are accurately described by the model atmosphere spectra developed over the last three decades. We compared RXTE data from a superburst with these spectra and found that the spectra predicted by previously published model atmospheres are strongly inconsistent with these high-precision measurements. In contrast, a simple Bose-Einstein spectrum is fully consistent with the data, as are recently published model atmosphere spectra. We discuss the implications of our results for determinations of neutron star masses and radii via constraints on their surface gravity and redshift, as originally suggested by Majczyna and Madej.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Low-Frequency Oscillations in Global Simulations of Black Hole Accretion

We have identified the presence of large-scale, low-frequency dynamo cycles in a long-duration, global, magnetohydrodynamic (MHD) simulation of black hole accretion. Such cycles had been seen previously in local shearing box simulations, but we discuss their evolution over 1,500 inner disk orbits of a global pi/4 disk wedge spanning two orders of magnitude in radius and seven scale heights in elevation above/below the disk midplane. The observed cycles manifest themselves as oscillations in azimuthal magnetic field occupying a region that extends into a low-density corona several scale heights above the disk. The cycle frequencies are ten to twenty times lower than the local orbital frequency, making them potentially interesting sources of low-frequency variability when scaled to real astrophysical systems. Furthermore, power spectra derived from the full time series reveal that the cycles manifest themselves at discrete, narrow-band frequencies that often share power across broad radial ranges. We explore possible connections between these simulated cycles and observed low-frequency quasi-periodic oscillations (LFQPOs) in galactic black hole binary systems, finding that dynamo cycles have the appropriate frequencies and are located in a spatial region associated with X-ray emission in real systems. Derived observational proxies, however, fail to feature peaks with RMS amplitudes comparable to LFQPO observations, suggesting that further theoretical work and more sophisticated simulations will be required to form a complete theory of dynamo-driven LFQPOs. Nonetheless, this work clearly illustrates that global MHD dynamos exhibit quasi-periodic behavior on timescales much longer than those derived from test particle considerations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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 γ-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its γ-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of γ-ray emission in the off-pulse phase interval at the ~ 6σ level. The flux from this emission in the energy range E \geq 100 MeV is F_100 = (1.73\pm0.40)\times10^(-8) photons/cm^2/s and is best fitted by a power law with a photon index of Γ = 2.54\pm0.14. The pulsed γ-ray flux in the same energy range is F_100 = (3.95\pm0.07)\times10^(-7) photons/cm^2/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of Γ = 1.41 \pm 0.23 and a cutoff energy E_c = 4.04 \pm 0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the γ-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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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 1 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SMBH Formation via Gas Accretion in Nuclear Stellar Clusters

Black holes exceeding a billion solar masses have been detected at redshifts greater than six. The rapid formation of these objects may suggest a massive early seed or a period of growth faster than Eddington. Here we suggest a new mechanism along these lines. We propose that in the process of hierarchical structure assembly, dense star clusters can be contracted on dynamical time scales due to the nearly free-fall inflow of self-gravitating gas with a mass comparable to or larger than that of the clusters. This increases the velocity dispersion to the point that the few remaining hard binaries can no longer effectively heat the cluster, and the cluster goes into a period of homologous core collapse. The cluster core can then reach a central density high enough for fast mergers of stellar-mass black holes and hence the rapid production of a black hole seed that could be $10^5 M_\odot$ or larger.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Effect of Massive Perturbers on Extreme Mass-Ratio Inspiral Waveforms

Extreme mass ratio inspirals, in which a stellar-mass object merges with a supermassive black hole, are prime sources for space-based gravitational wave detectors because they will facilitate tests of strong gravity and probe the spacetime around rotating compact objects. In the last few years of such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive to small perturbations. We show that one potentially detectable perturbation is the presence of a second supermassive black hole within a few tenths of a parsec. The acceleration produced by the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform to deviate systematically from that of an isolated system. If the perturber is a few tenths of a parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher derivatives of motion might also be detectable. In that case, the mass and distance of the perturber can be derived independently, which would allow a new probe of merger dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Deep Chandra Observation of the X-shaped Radio Galaxy 4C +00.58: A Candidate for Merger-induced Reorientation?

Although rapid reorientation of a black hole spin axis (lasting less than a few Myr) has been suggested as a mechanism for the formation of wings in X-shaped radio galaxies (XRGs), to date no convincing case of reorientation has been found in any XRG. Alternative wing formation models such as the hydrodynamic backflow models are supported by observed trends indicating that XRGs form preferentially with jets aligned along the major axis of the surrounding medium and wings along the minor axis. In this Letter, we present a deep Chandra observation of 4C +00.58, an oddball XRG with its jet oriented along the minor axis. By using the X-ray data in tandem with available radio and optical data, we estimate relevant timescales with which to evaluate wing formation models. The hydrodynamic models have difficulty explaining the long wings, whereas the presence of X-ray cavities (suggesting jet activity along a prior axis) and a potential stellar shell (indicating a recent merger) favor a merger-induced reorientation model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 out-of-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 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the ~ 1-10 Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, ~ 100-1000 solar masses. We primarily discuss two different source types -- mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-10 solar masses into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extreme Mass-Ratio Inspirals in the Effective-One-Body Approach: Quasi-Circular, Equatorial Orbits around a Spinning Black Hole

We construct effective-one-body waveform models suitable for data analysis with LISA for extreme-mass ratio inspirals in quasi-circular, equatorial orbits about a spinning supermassive black hole. The accuracy of our model is established through comparisons against frequency-domain, Teukolsky-based waveforms in the radiative approximation. The calibration of eight high-order post-Newtonian parameters in the energy flux suffices to obtain a phase and fractional amplitude agreement of better than 1 radian and 1 % respectively over a period between 2 and 6 months depending on the system considered. This agreement translates into matches higher than 97 % over a period between 4 and 9 months, depending on the system. Better agreements can be obtained if a larger number of calibration parameters are included. Higher-order mass ratio terms in the effective-one-body Hamiltonian and radiation-reaction introduce phase corrections of at most 30 radians in a one year evolution. These corrections are usually one order of magnitude larger than those introduced by the spin of the small object in a one year evolution. These results suggest that the effective-one-body approach for extreme mass ratio inspirals is a good compromise between accuracy and computational price for LISA data analysis purposes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling Extreme Mass Ratio Inspirals within the Effective-One-Body Approach

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between analytical EOB and numerical Teukolsky-based gravitational waveforms can be reduced to less than 10^(-1) rad and less than 2 x 10^(-3), respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of roughly 6-27 rad after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Super-Eddington fluxes during thermonuclear X-ray bursts

It has been known for nearly three decades that the energy spectra of thermonuclear X-ray bursts are often well-fit by Planck functions with temperatures so high that they imply a super-Eddington radiative flux at the emitting surface, even during portions of bursts when there is no evidence of photospheric radius expansion. This apparent inconsistency is usually set aside by assuming that the flux is actually sub-Eddington and that the fitted temperature is so high because the spectrum has been distorted by the energy-dependent opacity of the atmosphere. Here we show that the spectra predicted by currently available conventional atmosphere models appear incompatible with the highest-precision measurements of burst spectra made using the Rossi X-ray Timing Explorer, such as during the 4U 1820-30 superburst and a long burst from GX 17+2. In contrast, these measurements are well-fit by Bose-Einstein spectra with high temperatures and modest chemical potentials. Such spectra are very similar to Planck spectra. They imply surface radiative fluxes more than a factor of three larger than the Eddington flux. We find that segments of many other bursts from many sources are well-fit by similar Bose-Einstein spectra, suggesting that the radiative flux at the emitting surface also exceeds the Eddington flux during these segments. We suggest that burst spectra can closely approximate Bose-Einstein spectra and have fluxes that exceed the Eddington flux because they are formed by Comptonization in an extended, low-density radiating gas supported by the outward radiation force and confined by a tangled magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The drop of the coherence of the lower kHz QPOs is also observed in XTE J1701-462

We investigate the quality factor and RMS amplitude of the lower kHz QPOs from XTE J1701-462, a unique X-ray source which was observed in both the so-called Z and atoll states. Correcting for the frequency drift of the QPO, we show that, as in all sources for which such a correction can be applied, the quality factor and RMS amplitude drops sharply above above a critical frequency. For XTE J1701-462 this frequency is estimated to be ~800 Hz, where the quality factor reaches a maximum of ~200 (e.g. a value consistent with the one observed from more classical systems, such as 4U~1636-536). Such a drop has been interpreted as the signature of the innermost stable circular orbit, and that interpretation is consistent with the observations we report here. The kHz QPOs in the Z state are much less coherent and lower amplitude than they are in the atoll state. We argue that the change of the QPO properties between the two source states is related to the change of the scale height of the accretion disk; a prediction of the toy model proposed by barret et al. (2007). As a by-product of our analysis, we also increased the significance of the upper kHz QPO detected in the atoll phase up to 4.8 sigma (single trial significance), and show that the frequency separation (266.5+/-13.1 Hz) is comparable with the one measured from simultaneous twin QPOs the Z phase.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

kHz Quasi-Periodic Oscillations in the low-mass X-ray binary 4U 0614+09

We report on a comprehensive analysis of the kilohertz (above 300 Hz) quasi-periodic oscillations (kHz QPOs) detected from the neutron star low-mass X-ray binary 4U0614+09 with the Rossi X-ray Timing Explorer (RXTE). With a much larger data set than previously analyzed (all archival data from February 1996 up to October 2007), we first investigate the reality of the 1330 Hz QPO reported by van Straaten et al. (2000). This QPO would be of particular interest since it has the highest frequency reported for any source. A thorough analysis of the same observation fails to confirm the detection. On the other hand, over our extended data set, the highest QPO frequency we measure for the upper kHz QPO is at about 1224 Hz; a value which is fully consistent with the maximum values observed in similar systems. Second, we demonstrate that the frequency dependence of the quality factor and amplitude of the lower and upper kHz QPOs follow the systematic trends seen in similar systems (Barret et al., 2006). In particular, 4U0614+09 shows a drop of the quality factor of the lower kHz QPO above 700 Hz. If this is due to an approach to the innermost stable circular orbit, it implies a neutron star mass of about 1.9 solar masses. Finally, when analyzing the data over fixed durations, we have found a gap in the frequency distribution of the upper QPO, associated with a local minimum of its amplitude. A similar gap is not present in the distribution of the lower QPO frequencies, suggesting some cautions when interpreting frequency ratio distributions, based on the occurrence of the lower QPO only.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mass Segregation in NGC 2298: limits on the presence of an Intermediate Mass Black Hole

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Modeling Flows Around Merging Black Hole Binaries

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Neutron stars in Einstein-aether theory

As current and future experiments probe strong gravitational regimes around neutron stars and black holes, it is desirable to have theoretically sound alternatives to general relativity against which to test observations. Here we study the consequences of one such generalization, Einstein-aether theory, for the properties of non-rotating neutron stars. This theory has a parameter range that satisfies all current weak-field tests. We find that within this range it leads to lower maximum neutron star masses, as well as larger surface redshifts at a particular mass, for a given nuclear equation of state. For non-rotating black holes and neutron stars, the innermost stable circular orbit is only slightly modified in this theory.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Origin of intermittent accretion-powered X-ray oscillations in neutron stars with millisecond spin periods

We have shown previously that many of the properties of persistent accretion-powered millisecond pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving spot model" may also explain the intermittent accretion-powered pulsations that have been detected in three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very close to the spin axis to about 10 degrees away can increase the fractional rms amplitude from less than about 0.5 percent, which is usually undetectable with current instruments, to a few percent, which is easily detectable. The second harmonic of the spin frequency usually would not be detected, in agreement with observations. The model produces intermittently detectable oscillations for a range of emitting area sizes and beaming patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in stars that are more compact. In addition to explaining the sudden appearance of accretion-powered millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion-powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars but not others and suggests mechanisms that could explain the occasional temporal association of intermittent accretion-powered oscillations with thermonuclear X-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Systematic variation in the apparent burning area of thermonuclear bursts and its implication for neutron star radius measurement

Precision measurements of neutron star radii can provide a powerful probe of the properties of cold matter beyond nuclear density. Beginning in the late 1970s it was proposed that the radius could be obtained from the apparent or inferred emitting area during the decay portions of thermonuclear (type I) X-ray bursts. However, this apparent area is generally not constant, preventing reliable measurement of the source radius. Here we report for the first time a correlation between the variation of the inferred area and the burst properties, measured in a sample of almost 900 bursts from 43 sources. We found that the rate of change of the inferred area during decay is anticorrelated with the burst decay duration. A Spearman rank correlation test shows that this relation is significant at the <10^{-45} level for our entire sample, and at the 7x10^{-37} level for the 625 bursts without photospheric radius expansion. This anticorrelation is also highly significant for individual sources exhibiting a wide range of burst durations, such as 4U 1636-536 and Aql X-1. We suggest that variations in the colour factor, which relates the colour temperature resulted from the scattering in the neutron star atmosphere to the effective temperature of the burning layer, may explain the correlation. This in turn implies significant variations in the composition of the atmosphere between bursts with long and short durations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Chandra View of Nearby X-shaped Radio Galaxies

We present new and archival Chandra X-ray Observatory observations of X-shaped radio galaxies within z < 0.1 alongside a comparison sample of normal double-lobed FR I and II radio galaxies. By fitting elliptical distributions to the observed diffuse hot X-ray emitting atmospheres, we find that the ellipticity and the position angle of the hot gas follows that of the stellar light distribution for radio galaxy hosts in general. Moreover, compared to the control sample, we find a strong tendency for X-shaped morphology to be associated with wings directed along the minor axis of the hot gas distribution. Taken at face value, this result favors the hydrodynamic backflow models for the formation of X-shaped radio galaxies which naturally explain the geometry; the merger-induced rapid reorientation models make no obvious prediction about orientation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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 $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\simη^2$ (where $η$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\simη^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On the time variability of geometrically-thin black hole accretion disks I : the search for modes in simulated disks

We present a detailed temporal analysis of a set of hydrodynamic and magnetohydrodynamic (MHD) simulations of geometrically-thin (h/r~0.05) black hole accretion disks. The black hole potential is approximated by the Paczynski-Wiita pseudo-Newtonian potential. In particular, we use our simulations to critically assess two widely discussed models for high-frequency quasi-periodic oscillations, global oscillation modes (diskoseismology) and parametric resonance instabilities. We find that initially disturbed hydrodynamic disks clearly display the trapped global g-mode oscillation predicted by linear perturbation theory. In contrast, the sustained turbulence produced in the simulated MHD disks by the magneto-rotational instability does not excite these trapped g-modes. We cannot say at present whether the MHD turbulence actively damps the hydrodynamic g-mode. Our simulated MHD disks also fail to display any indications of a parametric resonance instability between the vertical and radial epicyclic frequencies. On the other hand, we do see characteristic frequencies at any given radius in the disk corresponding to local acoustic waves. We also conduct a blind search for any quasi-periodic oscillation in a proxy lightcurve based on the instantaneous mass accretion rate of the black hole, and place an upper limit of 2% on the total power in any such feature. We highlight the importance of correcting for secular changes in the simulated accretion disk when performing temporal analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

On the Time Variability of Geometrically-Thin Black Hole Accretion Disks II: Viscosity-Induced Global Oscillation Modes in Simulated Disks

We examine the evolution and influence of viscosity-induced diskoseismic modes in simulated black hole accretion disks. Understanding the origin and behavior of such oscillations will help us to evaluate their potential role in producing astronomically observed high-frequency quasi-periodic oscillations in accreting black hole binary systems. Our simulated disks are geometrically-thin with a constant half-thickness of five percent the radius of the innermost stable circular orbit. A pseudo-Newtonian potential reproduces the relevant effects of general relativity, and an alpha-model viscosity achieves angular momentum transport and the coupling of orthogonal velocity components in an otherwise ideal hydrodynamic numerical treatment. We find that our simulated viscous disks characteristically develop and maintain trapped global mode oscillations with properties similar to those expected of trapped g-modes and inner p-modes in a narrow range of frequencies just below the maximum radial epicyclic frequency. Although the modes are driven in the inner portion of the disk, they generate waves that propagate at the trapped mode frequencies out to larger disk radii. This finding is contrasted with the results of global magnetohydrodynamic disk simulations, in which such oscillations are not easily identified. Such examples underscore fundamental physical differences between accretion systems driven by the magneto-rotational instability and those for which alpha viscosity serves as a proxy for the physical processes that drive accretion, and we explore potential approaches to the search for diskoseismic modes in full magnetohydrodynamic disks.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar-mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass between ~50 and ~350 solar masses may be detectable by the planned advanced generation of ground-based GW interferometers. Such intermediate mass ratio inspirals (IMRIs) are most likely to be found in globular clusters. We analyze four possible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH binary via three-body interactions, (2) hardening via Kozai resonance in a hierarchical triple system, (3) direct capture, and (4) inspiral of a compact object from a tidally captured main-sequence star; we also discuss tidal effects when the inspiraling object is an NS. For each mechanism we predict the typical eccentricities of the resulting IMRIs. We find that IMRIs will have largely circularized by the time they enter the sensitivity band of ground-based detectors. Hardening of a binary via three-body interactions, which is likely to be the dominant mechanism for IMRI formation, yields eccentricities under 10^-4 when the GW frequency reaches 10 Hz. Even among IMRIs formed via direct captures, which can have the highest eccentricities, around 90% will circularize to eccentricities under 0.1 before the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences in globular clusters and the sensitivity of a network of three Advanced LIGO detectors to the resulting GWs. We show that this detector network may see up to tens of IMRIs per year, although rates of one to a few per year may be more plausible. We also estimate the loss in signal-to-noise ratio that will result from using circular IMRI templates for data analysis and find that, for the eccentricities we expect, this loss is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

An abrupt drop in the coherence of the lower kilohertz QPO in 4U 1636-536

Using archival data from the Rossi X-ray Timing Explorer (RXTE), we study in a systematic way the variation of the quality factor and amplitude of the lower and upper kHz quasi-periodic oscillations (QPO) in the low-mass X-ray binary 4U1636-536, over a frequency range from ~550 Hz to ~1200 Hz. When represented in a quality factor versus frequency diagram, the upper and lower QPOs follow two different tracks, suggesting that they are distinct phenomena, although not completely independent because the frequency difference of the two QPOs, when detected simultaneously, remains within ~60 Hz of half the neutron star spin frequency (at 581 Hz). The quality factor of the lower kHz QPO increases with frequency up to a maximum of ~200 at ~850 Hz, then drops precipitously to Q~50 at the highest detected frequencies ~920 Hz. A ceiling of the lower QPO frequencies at 920 Hz is also clearly seen in a frequency versus count rate diagram. The quality factor provides a measure of the coherence of the underlying oscillator. For exponentially damped sinusoidal shots, the highest Q observed corresponds to an oscillator coherence time of ~ 0.1 seconds. All existing QPO models face challenges in explaining such a long coherence time and the significantly different behaviours of the quality factors of the upper and lower QPOs reported here. It is therefore difficult to be certain of the implications of the abrupt change in the lower QPO at ~850 Hz. We discuss various possible causes, including that the drop in coherence is ultimately caused by effects related to the innermost stable circular orbit (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Binary Encounters With Supermassive Black Holes: Zero-Eccentricity LISA Events

Current simulations of the rate at which stellar-mass compact objects merge with supermassive black holes (called extreme mass ratio inspirals, or EMRIs) focus on two-body capture by emission of gravitational radiation. The gravitational wave signal of such events will likely involve a significant eccentricity in the sensitivity range of the Laser Interferometer Space Antenna (LISA). We show that tidal separation of stellar-mass compact object binaries by supermassive black holes will instead produce events whose eccentricity is nearly zero in the LISA band. Compared to two-body capture events, tidal separations have a high cross section and result in orbits that have a large pericenter and small apocenter. Therefore, the rate of interactions per binary is high and the resulting systems are very unlikely to be perturbed by other stars into nearly radial plunges. Depending on the fraction of compact objects that are in binaries within a few parsecs of the center, the rate of low-eccentricity LISA events could be comparable to or larger than the rate of high-eccentricity events.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Drop of coherence of the lower kilo-Hz QPO in neutron stars: is there a link with the innermost stable circular orbit?

Using all available archival data from the Rossi X-ray Timing Explorer (RXTE), we follow the frequency of the kilo-Hz QPOs in three low luminosity neutron star low mass X-ray binaries; namely 4U 1636-536, 4U 1608-522, and 4U1735-44. Following earlier work, we focus our analysis on the lower kilo-Hz QPO, for which we study the dependency of its quality factor (Q) amplitude as a function of frequency over a range covering from 500 Hz to 1000 Hz. As previously found for 4U 1636-536, we show that the quality factor of the lower kilo-Hz increases with frequency up to a maximum frequency around 800 Hz, beyond which an abrupt drop of its coherence is observed down to a limiting frequency where the QPO disappears completely. Simultaneously the amplitude of the QPOs is almost constant below the peak frequency and starts to decrease smoothly afterwards. The peak frequency is 850 Hz, 820 Hz, 740 Hz whereas the limiting frequency is 920 Hz, 900 Hz and 830 Hz for 4U 1636-536, 4U 1608-522 and 4U 1735-44 respectively. A ceiling of the lower QPO frequencies is also seen clearly in a frequency versus count rate diagram for all sources. This behavior is reproducible within an object and between objects. We suggest here that the drop of coherence of the lower QPO may be a geometry-related effect, which could be related to the last stable circular orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Intermediate-Mass Black Holes

The mathematical simplicity of black holes, combined with their links to some of the most energetic events in the universe, means that black holes are key objects for fundamental physics and astrophysics. Until recently, it was generally believed that black holes in nature appear in two broad mass ranges: stellar-mass (roughly 3-20 solar masses), which are produced by the core collapse of massive stars, and supermassive (millions to billions of solar masses), which are found in the centers of galaxies and are produced by a still uncertain combination of processes. In the last few years, however, evidence has accumulated for an intermediate-mass class of black holes, with hundreds to thousands of solar masses. If such objects exist they have important implications for the dynamics of stellar clusters, the formation of supermassive black holes, and the production and detection of gravitational waves. We review the evidence for intermediate-mass black holes and discuss future observational and theoretical work that will help clarify numerous outstanding questions about these objects.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Probing General Relativity With Mergers of Supermassive and Intermediate-Mass Black Holes

Recent observations and stellar dynamics simulations suggest that thousand solar mass black holes can form in compact massive young star clusters. Any such clusters in the bulge of their host galaxy will spiral to the center within a few hundred million years, where their intermediate-mass black holes are likely to merge eventually with the galaxy's supermassive black hole. If such mergers are common, then future space-based gravitational wave detectors such as the Laser Interferometer Space Antenna will detect them with such a high signal to noise ratio that towards the end of the inspiral the orbits will be visible in a simple power density spectrum, without the need for matched filtering. We discuss the astrophysics of the inspiral of clusters in the nuclear region of a galaxy and the subsequent merger of intermediate-mass with supermassive black holes. We also examine the prospects for understanding the spacetime geometry of rotating black holes, based on phase connection of the strong signals visible near the end of these extreme mass ratio inspirals.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Spectral Effects of the Vacuum Resonance in Soft Gamma-Ray Repeaters

The association of all three soft gamma-ray repeaters (SGRs) with supernova remnants has established that SGRs are young neutron stars, and has given us a starting point for detailed modeling. One of the most popular classes of models involves strongly magnetised neutron stars, with surface dipole fields B~ 10^{14}-10^{15} Gauss. In such strong magnetic fields, many otherwise negligible processes can play an important role. Here we consider the effects of vacuum polarisation on Compton scattering. Vacuum polarisation introduces a characteristic density-dependent photon frequency at which the normal modes of polarisation become nonorthogonal and the mean free path of photons decreases sharply. Our analytic results and Monte Carlo simulations of photon propagation through a magnetised plasma show that this effect leads, under a wide range of physical conditions, to a broad absorption-like feature in the energy range ~5 keV---40 keV. We discuss this effect in light of the spectra from SGR 1806-20.

0 citations0 reviewsNo code linkedOpen access