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.

The LIGO/Virgo collaboration has reported the detection of GW190412, a black hole-black hole (BH-BH) merger with the most unequal masses to date: m_1=24.4--34.7 Msun and m_2=7.4--10.1 Msun, corresponding to a mass ratio of q=0.21--0.41 (90 per cent probability range). Additionally, GW190412's effective spin was estimated to be X_eff=0.14--0.34, with the spin of the primary BH in the range a_spin=0.17--0.59. Based on this and prior detections, >10 per cent of BH-BH mergers have q<0.4. Major BH-BH formation channels (i.e., dynamics in dense stellar systems, classical isolated binary evolution, or chemically homogeneous evolution) tend to produce BH-BH mergers with comparable masses (typically with q>0.5). Here we test whether the classical isolated binary evolution channel can produce mergers resembling GW190412. We show that our standard binary evolution scenario, with the typical assumptions on input physics we have used in the past, produces such mergers. We provide an explicit example of an unequal mass BH-BH merger, which forms at low metallicity (Z=0.002) from two massive stars (77 and 35 Msun) and results in a BH-BH merger with m_1=27.0 Msun, m_2=9.9 Msun (q=0.37), with primary black hole spin a_spin=0.19 and an effective spin parameter X_eff=0.335. For this particular model of the input physics the overall BH-BH merger rate density in the local Universe (z=0) is: 73.5 Gpc^-3 yr^-1, while for systems with q<0.21, 0.28, 0.41, and 0.59 the rate density is: 0.01, 0.12, 6.8, and 22.2 Gpc^-3 yr^-1, respectively. The results from our standard model are consistent with the masses and spins of the black holes in GW190412, as well as with the LIGO/Virgo estimate of the fraction of unequal-mass BH-BH mergers.

It is currently not known if repeating fast radio bursts (FRBs) are fundamentally different from those that have not been seen to repeat. One striking difference between repeaters and apparent non-repeaters in the CHIME sample is that the once-off events are typically shorter in duration than sources that have been detected two or more times. We offer a simple explanation for this discrepancy based on a selection effect due to beamed emission, in which highly-beamed FRBs are less easily observed to repeat, but are abundant enough to detect often as once-off events. The explanation predicts that there is a continuous distribution of burst duration---not a bimodal one---with a correlation between repetition rate and width. Pulse width and opening angle may be related by relativistic effects in shocks, where short-duration bursts have small solid angles due to a large common Lorentz factor. Alternatively, the relationship could be a geometric effect where narrow beams sweep past the observer more quickly, as with pulsars. Our model has implications for the FRB emission mechanism and energy scale, volumetric event rates, and the application of FRBs to cosmology. It also provides a way of predicting which sources could be seen to repeat in the future.

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.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

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.

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.

Neutron star mergers can form a hypermassive neutron star (HMNS) remnant, which may be the engine of a short gamma ray burst (SGRB) before it collapses to a black hole, possibly several hundred milliseconds after the merger. During the lifetime of a HMNS, numerical relativity simulations indicate that it will undergo strong oscillations and emit GWs with frequencies of a few kilohertz, which are unfortunately too high for detection to be probable with Advanced LIGO. Here we discuss the current and future prospects for detecting these oscillations as modulation of the SGRB. The understanding of the physical mechanism responsible for the HMNS oscillations will provide information on the equation of state of the hot HMNS, and the observation of these frequencies in the SGRB data would give us insight into the emission mechanism of the SGRB.

The discovery of quasi-periodic brightness oscillations (QPOs) in the X-ray emission accompanying the giant flares of the soft gamma-ray repeaters SGR 1806-20 and SGR 1900+14 has led to intense speculation about their nature and what they might reveal about the interiors of neutron stars. Here we take a fresh look at the giant flare data for SGR 1806-20, and in particular we analyze short segments of the post-peak emission using a Bayesian procedure that has not previously been applied to these data. We find at best weak evidence that any QPO persists for more than $\sim 1$ second; instead, almost all the data are consistent with a picture in which there are numerous independently-excited modes that decay within a few tenths of a second. This has interesting implications for the rapidity of decay of the QPO modes, which could occur by the previously-suggested mechanism of coupling to the MHD continuum. The strongest QPOs favor certain rotational phases, which might suggest special regions of the crust or of the magnetosphere. We also find several previously unreported QPOs in these data, which may help in tracking down their origin.

Pulsating thermal X-ray emission from millisecond pulsars can be used to obtain constraints on the neutron star equation of state, but to date only five such sources have been identified. Of these five millisecond pulsars, only two have well constrained neutron star masses, which improve the determination of the radius via modelling of the X-ray waveform. We aim to find other millisecond pulsars that already have well constrained mass and distance measurements that show pulsed thermal X-ray emission in order to obtain tight constraints on the neutron star equation of state. The millisecond pulsar PSR~J1909--3744 has an accurately determined mass, M = 1.54$\pm$0.03 M$_\odot$ (1 $\sigma$ error) and distance, D = 1.07$\pm$0.04 kpc. We analysed {\em XMM-Newton} data of this 2.95 ms pulsar to identify the nature of the X-ray emission. We show that the X-ray emission from PSR~J1909--3744 appears to be dominated by thermal emission from the polar cap. Only a single component model is required to fit the data. The black-body temperature of this emission is kT=0.26\ud{0.03}{0.02} keV and we find a 0.2--10 keV un-absorbed flux of 1.1 $\times$ 10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ or an un-absorbed luminosity of 1.5 $\times$ 10$^{30}$ erg s$^{-1}$. Thanks to the previously determined mass and distance constraints of the neutron star PSR~J1909--3744, and its predominantly thermal emission, deep observations of this object with future X-ray facilities should provide useful constraints on the neutron star equation of state.

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

The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.

The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs.

Gravitational wave and electromagnetic observations can provide new insights into the nature of matter at supra-nuclear densities inside neutron stars. Improvements in electromagnetic and gravitational wave sensing instruments continue to enhance the accuracy with which they can measure the masses, radii, and tidal deformability of neutron stars. These better measurements place tighter constraints on the equation of state of cold matter above nuclear density. In this article, we discuss a complementary approach to get insights into the structure of neutron stars by providing a model prediction for non-linear fundamental eigenmodes (f-modes) and their decay over time, which are thought to be induced by time-dependent tides in neutron star binaries. Building on pioneering studies that relate the properties of f-modes to the structure of neutron stars, we systematically study this link in the non-perturbative regime using models that utilize numerical relativity. Using a suite of fully relativistic numerical relativity simulations of oscillating TOV stars, we establish blueprints for the numerical accuracy needed to accurately compute the frequency and damping times of f-mode oscillations, which we expect to be a good guide for the requirements in the binary case. We show that the resulting f-mode frequencies match established results from linear perturbation theory, but the damping times within numerical errors depart from linear predictions. This work lays the foundation for upcoming studies aimed at a comparison of theoretical models of f-mode signatures in gravitational waves, and their uncertainties with actual gravitational wave data, searching for neutron star binaries on highly eccentric orbits, and probing neutron star structure at high densities.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

A tentative excess in the electron spectrum at 1.4 TeV was recently reported by the DArk Matter Particle Explorer (DAMPE). A non-astrophysical scenario in which dark matter particles annihilate or decay in a local clump has been invoked to explain the excess. If $e^\pm$ annihilation channels in the final states are mediated by left-handed leptons as a component in the $SU(2)_L$ doublet, neutrinos with similar energies should have been simultaneously produced. We demonstrate that generic dark matter models can be decisively tested by the existing IceCube data. In case of a non-detection, such models would be excluded at the $5\sigma$ level by the five-year data for a point-like source and by the ten-year data for an extended source of dark matter particles with left-handed leptons.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

Graham et al. (2015a) reported a periodically varying quasar and supermassive black hole binary candidate, PG1302-102 (hereafter PG1302), which was discovered in the Catalina Real-Time Transient Survey (CRTS). Its combined Lincoln Near-Earth Asteroid Research (LINEAR) and CRTS optical light curve is well fitted to a sinusoid of an observed period of $\approx 1,884$ days and well modeled by the relativistic Doppler boosting of the secondary mini-disk (D'Orazio et al. 2015). However, the LINEAR+CRTS light curve from MJD $\approx 52700$ to MJD $\approx 56400$ covers only $\sim 2$ cycles of periodic variation, which is a short baseline that can be highly susceptible to normal, stochastic quasar variability (Vaughan et al. 2016). In this Letter, we present a re-analysis of PG1302, using the latest light curve from the All-Sky Automated Survey for Supernovae (ASAS-SN), which extends the observational baseline to the present day (MJD $\approx 58200$), and adopting a maximum likelihood method which searches for a periodic component in addition to stochastic quasar variability. When the ASAS-SN data are combined with the previous LINEAR+CRTS data, the evidence for periodicity decreases. For genuine periodicity one would expect that additional data would strengthen the evidence, so the decrease in significance may be an indication that the binary model is disfavored.

Surface abundance observations of halo stars hint at the occurrence of $r$-process nucleosynthesis at low metallicity ($\rm{[Fe/H]< -3}$), possibly within the first $10^8$ yr after the formation of the first stars. Possible loci of early-Universe $r$-process nucleosynthesis are the ejecta of either black hole--neutron star or neutron star--neutron star binary mergers. Here we study the effect of the inclination--eccentricity oscillations raised by a tertiary (e.g. a star) on the coalescence time scale of the inner compact object binaries. Our results are highly sensitive to the assumed initial distribution of the inner binary semi-major axes. Distributions with mostly wide compact object binaries are most affected by the third object, resulting in a strong increase (by more than a factor of 2) in the fraction of fast coalescences. If instead the distribution preferentially populates very close compact binaries, general relativistic precession prevents the third body from increasing the inner binary eccentricity to very high values. In this last case, the fraction of coalescing binaries is increased much less by tertiaries, but the fraction of binaries that would coalesce within $10^8$ yr even without a third object is already high. Our results provide additional support to the compact object merger scenario for $r$-process nucleosynthesis.

We summarize and reanalyze observations bearing upon missing galactic baryons, where we propose a consistent picture for halo gas in L >~ L* galaxies. The hot X-ray emitting halos are detected to 50-70 kpc, where typically, M_hot(<50 kpc) ~ 5E9 Msun, and with density n \propto r^-3/2. When extrapolated to R200, the gas mass is comparable to the stellar mass, but about half of the baryons are still missing from the hot phase. If extrapolated to 1.9-3 R200, the baryon to dark matter ratio approaches the cosmic value. Significantly flatter density profiles are unlikely for R < 50 kpc and they are disfavored but not ruled out for R > 50 kpc. For the Milky Way, the hot halo metallicity lies in the range 0.3-1 solar for R < 50 kpc. Planck measurements of the thermal Sunyaev-Zeldovich effect toward stacked luminous galaxies (primarily early-type) indicate that most of their baryons are hot, near the virial temperature, and extend beyond R200. This stacked SZ signal is nearly an order of magnitude larger than that inferred from the X-ray observations of individual (mostly spiral) galaxies with M_* > 10^11.3 Msun. This difference suggests that the hot halo properties are distinct for early and late type galaxies, possibly due to different evolutionary histories. For the cooler gas detected in UV absorption line studies, we argue that there are two absorption populations: extended halos; and disks extending to ~50 kpc, containing most of this gas, and with masses a few times lower than the stellar masses. Such extended disks are also seen in 21 cm HI observations and in simulations.

In the past few years wide-field optical and UV transient surveys as well as X-ray telescopes have allowed us to identify a few dozen candidate tidal disruption events (TDEs). While in theory the physical processes in TDEs are expected to be ubiquitous, a few distinct classes of TDEs have been observed. Some TDEs radiate mainly in NUV/optical while others produce prominent X-rays. Moreover, relativistic jets have been observed in only a handful of TDEs. This diversity might be related to the details of the super-Eddington accretion and emission physics relevant to TDE disks. In this Letter, we utilize novel three-dimensional general relativistic radiation magnetohydrodynamics simulations to study the super-Eddington compact disk phase expected in TDEs. Consistent with previous studies, geometrically thick disks, wide-angle optically-thick fast outflows and relativistic jets are produced. The outflow density and velocity depend sensitively on the inclination angle, and hence so does the reprocessing of emission produced from the inner disk. We then use Monte-Carlo radiative transfer to calculate the reprocessed spectra and find that that the observed ratio of optical to X-ray fluxes increases with increasing inclination angle. This naturally leads to a unified model for different classes of TDEs in which the spectral properties of the TDE depend mainly on the viewing-angle of the observer with respect to the orientation of the disk.

This erratum addresses an error in the text of Lo et al. (2013) and more importantly corrects a computational error that was made in all previously published relativistic computations of absolute pulse waveform fluxes produced by hot spots on rotating stars that followed the radiation forward in time from the star to the observer. These computations used an incorrect expression for the differential emitting area on the stellar surface, as measured in the comoving frame. This erratum explains this error, lists all the published papers that contain results affected by this error, and summarizes in a general way why and how the results in Lo et al. (2013) were affected by this error, in a way that makes clear why and how the results in the other papers were affected. The idealized radiation spectra, beaming patterns, and spot shapes assumed in many of the previously reported computations of pulse waveforms very likely introduced systematic errors that are much larger than the relatively small errors introduced by using an incorrect expression for the differential emitting area. This erratum also corrects an error in the caption of Figure 5, an error in the text on page 19, and an incomplete description of our results on page 24.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

Observations of thermonuclear X-ray bursts from accreting neutron stars (NSs) in low-mass X-ray binary systems can be used to constrain NS masses and radii. Most previous work of this type has set these constraints using Planck function fits as a proxy: both the models and the data are fit with diluted blackbody functions to yield normalizations and temperatures which are then compared against each other. Here, for the first time, we fit atmosphere models of X-ray bursting NSs directly to the observed spectra. We present a hierarchical Bayesian fitting framework that uses state-of-the-art X-ray bursting NS atmosphere models with realistic opacities and relativistic exact Compton scattering kernels as a model for the surface emission. We test our approach against synthetic data, and find that for data that are well-described by our model we can obtain robust radius, mass, distance, and composition measurements. We then apply our technique to Rossi X-ray Timing Explorer observations of five hard-state X-ray bursts from 4U 1702-429. Our joint fit to all five bursts shows that the theoretical atmosphere models describe the data well but there are still some unmodeled features in the spectrum corresponding to a relative error of 1-5% of the energy flux. After marginalizing over this intrinsic scatter, we find that at 68% credibility the circumferential radius of the NS in 4U 1702-429 is R = 12.4+-0.4 km, the gravitational mass is M=1.9+-0.3 Msun, the distance is 5.1 < D/kpc < 6.2, and the hydrogen mass fraction is X < 0.09.

The first direct detections of gravitational waves from merging binary black holes open a unique window into the binary black hole formation environment. One promising environmental signature is the angular distribution of the black hole spins; systems formed through dynamical interactions among already-compact objects are expected to have isotropic spin orientations whereas binaries formed from pairs of stars born together are more likely to have spins preferentially aligned with the binary orbital angular momentum. We consider existing gravitational wave measurements of the binary effective spin, the best-measured combination of spin parameters, in the four likely binary black hole detections GW150914, LVT151012, GW151226, and GW170104. If binary black hole spin magnitudes extend to high values we show that the data exhibit a $2.4\sigma$ ($0.015$ odds ratio) preference for an isotropic angular distribution over an aligned one. By considering the effect of 10 additional detections, we show that such an augmented data set would enable in most cases a preference stronger than $5\sigma$ ($2.9 \times 10^{-7}$ odds ratio). The existing preference for either an isotropic spin distribution or low spin magnitudes for the observed systems will be confirmed (or overturned) confidently in the near future.

As gravitational wave instrumentation becomes more sensitive, it is interesting to speculate about subtle effects that could be analyzed using upcoming generations of detectors. One such effect that has great potential for revealing the properties of very dense matter is fluid oscillations of neutron stars. These have been found in numerical simulations of the hypermassive remnants of double neutron star mergers and of highly eccentric neutron star orbits. Here we focus on the latter and sketch out some ideas for the production, gravitational-wave detection, and analysis of neutron star oscillations. These events will be rare (perhaps up to several tens per year could be detected using third-generation detectors such as the Einstein Telescope or the Cosmic Explorer), but they would have unique diagnostic power for the analysis of cold, catalyzed, dense matter. Furthermore, these systems are unusual in that analysis of the tidally excited f-modes of the stars could yield simultaneous measurements of their masses, moments of inertia, and tidal Love numbers, using the frequency, damping time, and amplitude of the modes. They would thus present a nearly unique opportunity to test observationally the I-Love-Q relation. The analysis of such events will require significant further work in nuclear physics and general relativistic nonlinear mode coupling, and thus we discuss further directions that will need to be pursued. For example, we note that for nearly-grazing encounters, numerical simulations show that the energy delivered to the f-modes may be up to two orders of magnitude greater than predicted in the linear theory.

The detection of ultrahigh-energy (UHE) neutrino sources would contribute significantly to solving the decades-old mystery of the origin of the highest-energy cosmic rays. We investigate the ability of a future UHE neutrino detector to identify the brightest neutrino point sources, by exploring the parameter space of the total number of observed events and the angular resolution of the detector. The favored parameter region can be translated to requirements for the effective area, sky coverage and angular resolution of future detectors, for a given source number density and evolution history. Moreover, by studying the typical distance to sources that are expected to emit more than one event for a given diffuse neutrino flux, we find that a significant fraction of the identifiable UHE neutrino sources may be located in the nearby Universe if the source number density is above $\sim10^{-6}\,\rm Mpc^{-3}$. If sources are powerful and rare enough, as predicted in blazar scenarios, they can first be detected at distant locations. Our result also suggests that if UHE cosmic-ray accelerators are neither beamed nor transients, it will be possible to associate the detected UHE neutrino sources with nearby UHE cosmic-ray and gamma-ray sources, and that they may also be observed using other messengers, including ones with limited horizons such as TeV gamma rays, UHE gamma rays and cosmic rays. We find that for a $\gtrsim5\sigma$ detection of UHE neutrino sources with a uniform density, $n_s\sim{10}^{-7}-{10}^{-5}~{\rm Mpc}^{-3}$, at least $\sim100-1000$ events and sub-degree angular resolution are needed, and the results depend on the source evolution model.

The possible Fermi detection of an electromagnetic counterpart to the double black hole merger GW150914 has inspired many theoretical models, some of which propose that the holes spiraled together inside a massive star. However, we show that the heat produced by the dynamical friction on such black hole orbits can exceed the stellar binding energy by a large factor, which means that this heat could destroy the star and thus make it difficult for enough gas to be near the holes at merger to produce detectable photons. These considerations must be taken into account when models are proposed for electromagnetic counterparts to the coalescence of two stellar-mass black holes. We find that only when the two black holes form very close to the center can the star avoid destruction. In that case, dynamical friction can make the black holes coalesce faster than they would in vacuum, which leads to a modification of the gravitational waveform that is potentially observable by advanced LIGO.

Disk and wind signatures are seen in the soft state of Galactic black holes, while the jet is seen in the hard state. Here we study the disk-wind connection in the $\rho$ class of variability in GRS 1915+105 using a joint NuSTAR-Chandra observation. The source shows 50 sec limit cycle oscillations. By including new information provided by the reflection spectrum, and using phase-resolved spectroscopy, we find that the change in the inner disk inferred from the blackbody emission is not matched by reflection measurements. The latter is almost constant, independent of the continuum model. The two radii are comparable only if the disk temperature color correction factor changes, an effect that could be due to the changing opacity of the disk caused by changes in metal abundances. The disk inclination is similar to that inferred from the jet axis, and oscillates by ~10 deg. The simultaneous Chandra data show the presence of two wind components with velocities between 500-5000 km/s, and possibly two more with velocities reaching 20,000 km/s (~0.06 c). The column densities are ~5e22 cm$^{-2}$. An upper limit to the wind response time of 2 sec is measured, implying a launch radius of <6e10 cm. The changes in wind velocity and absorbed flux require the geometry of the wind to change during the oscillations, constraining the wind to be launched from a distance of 290 - 1300 rg from the black hole. Both datasets support fundamental model predictions in which a bulge originates in the inner disk and moves outward as the instability progresses.

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.

The Fermi bubbles are two lobes filled with non-thermal particles that emit gamma rays, extend $\approx$10 kpc vertically from the Galactic center, and formed from either nuclear star formation or accretion activity on Sgr A*. Simulations predict a range of shock strengths as the bubbles expand into the surrounding hot gas halo distribution ($T_{halo} \approx 2 \times 10^6$ K), but with significant uncertainties in the energetics, age, and thermal gas structure. The bubbles should contain thermal gas with temperatures between $10^6$ and $10^8$ K, with potential X-ray signatures. In this work, we constrain the bubbles' thermal gas structure by modeling the OVII and OVIII emission line strengths from archival XMM-Newton and Suzaku data. Our emission model includes a hot thermal volume-filled bubble component cospatial with the gamma-ray region, and a shell of compressed material. We find that a bubble/shell model with $n \approx 1 \times 10^{-3}$ cm$^{-3}$ and with log($T$) $\approx$ 6.60-6.70 is consistent with the observed line intensities. In the framework of a continuous Galactic outflow, we infer a bubble expansion rate, age, and energy injection rate of $490_{-77}^{+230}$ km s$^{-1}$, $4.3_{-1.4}^{+0.8}$ Myr, and $2.3_{-0.9}^{+5.1} \times 10^{42}$ erg s$^{-1}$. These estimates are consistent with the bubbles forming from a Sgr A* accretion event rather than from nuclear star formation.

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.

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.

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.

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.

The hot gaseous halos of galaxies likely contain a large amount of mass and are an integral part of galaxy formation and evolution. The Milky Way has a 2e6 K halo that is detected in emission and by absorption in the OVII resonance line against bright background AGNs, and for which the best current model is an extended spherical distribution. Using XMM-Newton RGS data, we measure the Doppler shifts of the OVII absorption-line centroids toward an ensemble of AGNs. These Doppler shifts constrain the dynamics of the hot halo, ruling out a stationary halo at about 3sigma and a corotating halo at 2sigma, and leading to a best-fit rotational velocity of 183+/-41 km/s for an extended halo model. These results suggest that the hot gas rotates and that it contains an amount of angular momentum comparable to that in the stellar disk. We examined the possibility of a model with a kinematically distinct disk and spherical halo. To be consistent with the emission-line X-ray data the disk must contribute less than 10% of the column density, implying that the Doppler shifts probe motion in the extended hot halo.

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.

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]

About 30-50% of the baryons in the local Universe are unaccounted for and are likely in a hot phase, 10^5.5-10^8 K. A hot halo (10^6.3 K) is detected around the Milky Way through the O VII and O VIII resonance absorption and emission lines in the soft X-ray band. Current instruments are not sensitive enough to detect this gas in absorption around other galaxies and galaxy groups, the two most likely sites. We show that resonant line absorption by this hot gas can be detected with current technology, with a collecting area exceeding about 300 cm^2 and a resolution R > 2000. For a few notional X-ray telescope configurations that could be constructed as Explorer or Probe missions, we calculate the differential number of O VII and O VIII absorbers as a function of equivalent width through redshift space, dN/dz. The hot halos of individual external galaxies produce absorption that should be detectable out to about their virial radii. For the Milky Way, one can determine the radial distribution of density, temperature, and metallicity, after making optical depth corrections. Spectroscopic observations can determine the rotation of a hot gaseous halo.

Detections of $z \approx$ 0 oxygen absorption and emission lines indicate the Milky Way hosts a hot ($\sim 10^6$ K), low-density plasma extending $\gtrsim$50 kpc into the Mily Way's halo. Current X-ray telescopes cannot resolve the line profiles, but the variation of their strengths on the sky constrains the radial gas distribution. Interpreting the OVII K$\alpha$ absorption line strengths has several complications, including optical depth and line of sight velocity effects. Here, we present model absorption line profiles accounting for both of these effects to show the lines can exhibit asymmetric structures and be broader than the intrinsic Doppler width. The line profiles encode the hot gas rotation curve, the net inflow or outflow of hot gas, and the hot gas angular momentum profile. We show how line of sight velocity effects impact the conversion between equivalent width and the column density, and provide modified curves of growth accounting for these effects. As an example, we analyze the LMC sight line pulsar dispersion measure and OVII equivalent width to show the average gas metallicity is $\gtrsim 0.6 Z_{\odot}$ and $b$ $\gtrsim$ 100 km s$^{-1}$. Determining these properties offers valuable insights into the dynamical state of the Milky Way's hot gas, and improves the line strength interpretation. We discuss future strategies to observe these effects with an instrument that has a spectral resolution of about 3000, a goal that is technically possible today.

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.

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 $\beta$, 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 $\beta$, 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 $\beta> 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.

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.

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.

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.

The Milky Way hosts a hot ($\approx 2 \times 10^6$ K), diffuse, gaseous halo based on detections of z = 0 OVII and OVIII absorption lines in quasar spectra and emission lines in blank-sky spectra. Here we improve constraints on the structure of the hot gas halo by fitting a radial model to a much larger sample of OVII and OVIII emission line measurements from XMM-Newton EPIC-MOS spectra compared to previous studies ($\approx$ 650 sightlines). We assume a modified $\beta$-model for the halo density distribution and a constant-density Local Bubble from which we calculate emission to compare with the observations. We find an acceptable fit to the OVIII emission line observations with $\chi^{2}_{red}$ (dof) = 1.08 (644) for best-fit parameters of $n_o r_c^{3\beta} = 1.35 \pm 0.24$ cm$^{-3}$ kpc$^{3\beta}$ and $\beta = 0.50 \pm 0.03$ for the hot gas halo and negligible Local Bubble contribution. The OVII observations yield an unacceptable $\chi^{2}_{red}$ (dof) = 4.69 (645) for similar best-fit parameters, which is likely due to temperature or density variations in the Local Bubble. The OVIII fitting results imply hot gas masses of $M$(< 50 kpc) = $3.8_{-0.3}^{+0.3} \times 10^{9} M_{\odot}$ and $M$(< 250 kpc) = $4.3_{-0.8}^{+0.9} \times 10^{10} M_{\odot}$, accounting for $\lesssim$ 50% of the Milky Way's missing baryons. We also explore our results in the context of optical depth effects in the halo gas, the halo gas cooling properties, temperature and entropy gradients in the halo gas, and the gas metallicity distribution. The combination of absorption and emission line analyses implies a sub-solar gas metallicity that decreases with radius, but that also must be $\geq 0.3 Z_{\odot}$ to be consistent with the pulsar dispersion measure toward the LMC.

CoGeNT has taken data for over 3 years, with 1136 live days of data accumulated as of April 23, 2013. We report on the results of a maximum likelihood analysis to extract any possible dark matter signal present in the collected data. The maximum likelihood signal extraction uses 2-dimensional probability density functions (PDFs) to characterize the anticipated variations in dark matter interaction rates for given observable nuclear recoil energies during differing periods of the Earth's annual orbit around the Sun. Cosmogenic and primordial radioactivity backgrounds are characterized by their energy signatures and in some cases decay half-lives. A third parameterizing variable -- pulse rise-time -- is added to the likelihood analysis to characterize slow rising pulses described in prior analyses. The contribution to each event category is analyzed for various dark matter signal hypotheses including a dark matter standard halo model and a case with free oscillation parameters (i.e., amplitude, period, and phase). The best-fit dark matter signal is in close proximity to previously reported results. We find that the significance of the extracted dark matter signal remains well below evidentiary at 1.7 $\sigma$.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

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.

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.

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 $\lambda_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 \lambda_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.

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.

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.

Weakly Interacting Massive Particles (WIMPs) are well-established dark matter candidates. WIMP interactions with sensitive detectors are expected to display a characteristic annual modulation in rate. We release a dataset spanning 3.4 years of operation from a low-background germanium detector, designed to search for this signature. A previously reported modulation persists, concentrated in a region of the energy spectrum populated by an exponential excess of unknown origin. Its phase and period agree with phenomenological expectations, but its amplitude is a factor $\sim$4-7 larger than predicted for a standard WIMP galactic halo. We consider the possibility of a non-Maxwellian local halo velocity distribution as a plausible explanation, able to help reconcile recently reported WIMP search anomalies.

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.

Soft, potentially thermal spectral components observed in some ULXs can be fit with models for emission from cool, optically-thick accretion disks. If that description is correct, the low temperatures that are observed imply accretion onto "intermediate-mass" black holes. Subsequent work has found that these components may follow an inverse relationship between luminosity and temperature, implying a non-blackbody origin for this emission. We have re-analyzed numerous XMM-Newton spectra of extreme ULXs. Crucially, observations wherein the source fell on a chip gap were excluded owing to their uncertain flux calibration, and the neutral column density along the line of sight to a given source was jointly determined by multiple spectra. The luminosity of the soft component is found to be positively correlated with temperature, and to be broadly consistent with L ~ T^4 in the measured band pass, as per blackbody emission from a standard thin disk. These results are nominally consistent with accretion onto black holes with masses above the range currently known in Galactic X-ray binaries, though there are important caveats. Emission from inhomogeneous or super-Eddington disks may also be consistent with the data.

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

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.

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.

The Milky Way's million degree gaseous halo contains a considerable amount of mass that, depending on its structural properties, can be a significant mass component. In order to analyze the structure of the Galactic halo, we use XMM-Newton Reflection Grating Spectrometer archival data and measure OVII K alpha absorption-line strengths toward 26 active galactic nuclei, LMC X-3, and two Galactic sources (4U 1820-30 and X1735-444). We assume a beta-model as the underlying gas density profile and find best-fit parameters of n_o = 0.46^{+0.74}_{-0.35} cm^-3, r_c = 0.35^{+0.29}_{-0.27} kpc, and beta = 0.71^{+0.13}_{-0.14}. These parameters result in halo masses ranging between M(18 kpc) = 7.5^{+22.0}_{-4.6} x 10^8 M_sun and M(200 kpc) = 3.8^{+6.0}_{-0.5} x 10^{10} M_sun assuming a gas metallicity of Z = 0.3 Z_sun, which are consistent with current theoretical and observational work. The maximum baryon fraction from our halo model of f_b = 0.07^{+0.03}_{-0.01} is significantly smaller than the universal value of f_b = 0.171, implying the mass contained in the Galactic halo accounts for 10% - 50% of the missing baryons in the Milky Way. We also discuss our model in the context of several Milky Way observables, including ram pressure stripping in dwarf spheroidal galaxies, the observed X-ray emission measure in the 0.5 - 2 keV band, the Milky Way's star formation rate, spatial and thermal properties of cooler gas (~10^5 K) and the observed Fermi bubbles toward the Galactic center. Although the metallicity of the halo gas is a large uncertainty in our analysis, we place a lower limit on the halo gas between the Sun and the Large Magellanic Cloud (LMC). We find that Z >~ 0.2 Z_sun based on the pulsar dispersion measure toward the LMC.

CoGeNT employs p-type point-contact (PPC) germanium detectors to search for Weakly Interacting Massive Particles (WIMPs). By virtue of its low energy threshold and ability to reject surface backgrounds, this type of device allows an emphasis on low-mass dark matter candidates (wimp mass of about 10 GeV/c2). We report on the characteristics of the PPC detector presently taking data at the Soudan Underground Laboratory, elaborating on aspects of shielding, data acquisition, instrumental stability, data analysis, and background estimation. A detailed background model is used to investigate the low energy excess of events previously reported, and to assess the possibility of temporal modulations in the low-energy event rate. Extensive simulations of all presently known backgrounds do not provide a viable background explanation for the excess of low-energy events in the CoGeNT data, or the previously observed temporal variation in the event rate. Also reported on for the first time is a determination of the surface (slow pulse rise time) event contamination in the data as a function of energy. We conclude that the CoGeNT detector technology is well suited to search for the annual modulation signature expected from dark matter particle interactions in the region of WIMP mass and coupling favored by the DAMA/LIBRA results

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.

Gas accreting onto a galaxy will be of low metallicity while halo gas due to a galactic fountain will be of near-solar metallicity. We test these predictions by measuring the metal absorption line properties of halo gas 5 kpc above the plane of the edge-on galaxy NGC 891, using observations taken with HST/STIS toward a bright background quasar. Metal absorption lines of Fe II, Mg II, and Mg I in the halo of NGC 891 are clearly seen, and when combined with recent deep H I observations, we are able to place constraints on the metallicity of the halo gas for the first time. The H I line width defines the line broadening, from which we model opacity effects in these metal lines, assuming the absorbing gas is continuously distributed in the halo. The gas-phase metallicities are [Fe/H] = -1.18 +/- 0.07 and [Mg/H] = -0.23 +0.36/-0.27 (statistical errors) and this difference is probably due to differential depletion onto grains. When corrected for such depletion using Galactic gas as a guide, both elements have approximately solar or even supersolar abundances. This suggests that the gas is from the galaxy disk, probably expelled into the halo by a galactic fountain, rather than from accretion of intergalactic gas, which would have a low metallicity. The abundances would be raised by significant amounts if the absorbing gas lies in a few clouds with thermal widths smaller than the rotational velocity of the halo. If this is the case, both the abundances and [Mg/Fe] would be supersolar.

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.

We have carried out a comprehensive study of the molecular conditions and star-forming activities in dark cloud L1174 with multi-wavelength data. Mapping observations of L1174 in $^{13}$CO $J=2-1$ and $^{12}$CO $J=3-2$ were performed using the KOSMA 3-meter telescope. Six molecular cores with masses ranging from 5 to 31 $M_\odot$ and sizes ranging from 0.17 to 0.39 pc are resolved. Large area ahead of a Herbig Be star, HD 200775, is in expanding and core 1 is with collapse signature. Large line widths of $^{13}$CO $J=2-1$ indicate the ubiquity of turbulent motions in this region. Spectra of $^{12}$CO $J=3-2$ prevalently show conspicuously asymmetric double-peaked profiles. In a large area, red-skewed profiles are detected and suggestive of a scenario of global expansion. There is a large cavity around the Herbig Be star HD 200775, the brightest star in L1174. The gas around the cavity has been severely compressed by the stellar winds from HD 200775. Feedbacks from HD 200775 may have helped form the molecular cores around the cavity. Seventeen 2MASS potential young stellar objects were identified according to their 2MASS colour indices. The spatial distribution of the these 2MASS sources indicates that some of them have a triggered origin. All these suggest that feedbacks from a Herbig Ae/Be star may also have the potential to trigger star forming activities.

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.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

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.

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.

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.

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.

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.

Fifteen months of cumulative CoGeNT data are examined for indications of an annual modulation, a predicted signature of Weakly Interacting Massive Particle (WIMP) interactions. Presently available data support the presence of a modulated component of unknown origin, with parameters prima facie compatible with a galactic halo composed of light-mass WIMPs. Unoptimized estimators yield a statistical significance for a modulation of ~2.8 sigma, limited by the short exposure.

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.

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of {\gamma}-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its {\gamma}-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of {\gamma}-ray emission in the off-pulse phase interval at the ~ 6{\sigma} 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 {\Gamma} = 2.54\pm0.14. The pulsed {\gamma}-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 {\Gamma} = 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 {\gamma}-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees.

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.

With the KOSMA 3-m telescope, 54 Herbig Ae/Be stars were surveyed in CO and $^{13}$CO emission lines. The properties of the stars and their circumstellar environments are studied by fitting the SEDs. The mean line width of $^{13}$CO (2-1) lines of this sample is 1.87 km s$^{-1}$. The average column density of H$_{2}$ is found to be $4.9\times10^{21}$ cm$^{-2}$ for the stars younger than $10^{6}$ yr, while drops to $2.5\times10^{21}$ cm$^{-2}$ for those older than $10^{6}$ yr. No significant difference is found among the SEDs of Herbig Ae stars and Herbig Be stars at the same age. The infrared excess decreases with age. The envelope masses and the envelope accretion rates decease with age after $10^{5}$ yr. The average disk mass of the sample is $3.3\times10^{-2} M_{\sun}$. The disk accretion rate decreases more slowly than the envelope accretion rate. A strong correlation between the CO line intensity and the envelope mass is found.

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.

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.

It is now clear that a subset of supernovae display evidence for jets and are observed as gamma-ray bursts. The angular momentum distribution of massive stellar endpoints provides a rare means of constraining the nature of the central engine in core-collapse explosions. Unlike supermassive black holes, the spin of stellar-mass black holes in X-ray binary systems is little affected by accretion, and accurately reflects the spin set at birth. A modest number of stellar-mass black hole angular momenta have now been measured using two independent X-ray spectroscopic techniques. In contrast, rotation-powered pulsars spin-down over time, via magnetic braking, but a modest number of natal spin periods have now been estimated. For both canonical and extreme neutron star parameters, statistical tests strongly suggest that the angular momentum distributions of black holes and neutron stars are markedly different. Within the context of prevalent models for core-collapse supernovae, the angular momentum distributions are consistent with black holes typically being produced in GRB-like supernovae with jets, and with neutron stars typically being produced in supernovae with too little angular momentum to produce jets via magnetohydrodynamic processes. It is possible that neutron stars are imbued with high spin initially, and rapidly spun-down shortly after the supernova event, but the available mechanisms may be inconsistent with some observed pulsar properties.

We have carried out CO J=2-1 and CO J=3-2 observations toward Tycho's supernova remnant (SNR) using the KOSMA 3m-Telescope. From these observations we identified three molecular clouds (MCs) around the SNR. The small cloud in the southwest was discovered for the first time. In the north and east, two MCs (cloud A and cloud B) adjacent in space display a bow-shaped morphology, and have broad emission lines, which provide some direct evidences of the SNR-MCs interaction. The MCs is revealed at -69~-59 km \s, well coincident with the Tycho's SNR. The MCs associated with Tycho's SNR have a mass of ~2130 sun masses. Position-velocity diagrams show the two clouds adjacent in velocity which means possible cloud-cloud collision in this region. The maximum value (0.66) of integrated CO line intensity ratio I(CO J=3-2/CO J=2-1) for the three MCs agrees well with the previous measurement of individual Galactic MCs, implying that the SNR shock just drove into the MCs. The two MCs have a line intensity ratio gradient. The distribution of the ratio appears to indicate that the shock propagates from the southwest to the northeast.

Aims: We aim at deriving the excitation conditions of the interstellar gas as well as the local FUV intensities in the molecular cloud surrounding NGC 3603 to get a coherent picture of how the gas is energized by the central stars. Methods: The NANTEN2-4m submillimeter antenna is used to map the [CI] 1-0, 2-1 and CO 4-3, 7-6 lines in a 2' x 2' region around the young OB cluster NGC 3603 YC. These data are combined with C18O 2-1 data, HIRES-processed IRAS 60 and 100 micron maps of the FIR continuum, and Spitzer/IRAC maps. Results: The NANTEN2 observations show the presence of two molecular clumps located south-east and south-west of the cluster and confirm the overall structure already found by previous CS and C18O observations. We find a slight position offset of the peak intensity of CO and [CI], and the atomic carbon appears to be further extended compared to the molecular material. We used the HIRES far-infrared dust data to derive a map of the FUV field heating the dust. We constrain the FUV field to values of \chi = 3 - 6 \times 10^3 in units of the Draine field across the clouds. Approximately 0.2 to 0.3 % of the total FUV energy is re-emitted in the [CII] 158 {\mu}m cooling line observed by ISO. Applying LTE and escape probability calculations, we derive temperatures (TMM1 = 43 K, TMM2 = 47 K), column densities (N(MM1) = 0.9 \times 10^22 cm^-2, N(MM2) = 2.5 \times 10^22 cm^-2) and densities (n(MM1) = 3 \times 10^3 cm^-3, n(MM2) = 10^3 -10^4 cm^-3) for the two observed molecular clumps MM1 and MM2. Conclusions: The cluster is strongly interacting with the ambient molecular cloud, governing its structure and physical conditions. A stability analysis shows the existence of gravitationally collapsing gas clumps which should lead to star formation. Embedded IR sources have already been observed in the outskirts of the molecular cloud and seem to support our conclusions.

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.

We have performed the submillimeter and millimeter observations in CO lines toward supernova remnant (SNR) IC443. CO molecular shell is well coincident with the partial shell of the SNR detected in radio continuum observations. The broad emission lines and three 1720-MHz OH masers were detected in CO molecular shell. The present observations have provided further evidence in support of the interaction between the SNR and the adjoining molecular clouds (MCs). The total mass of the MCs is 9.26*10^{3} sun masses. The integrated CO line intensity ratio I(CO(3-2))/I(CO(2-1)) for the whole molecular cloud is between 0.79 and 3.40. The average value is 1.58, which is much higher than previous measurements of individual Galactic MCs. Higher line ratios imply that shocks have driven into the MCs. We conclude that high I(CO(3-2))/I(CO(2-1)) is identified as one good signature of SNR-MCs interacting system. Based on $IRAS$ Point Source Catalog and 2MASS near-infrared database, 12 protostellar objects and 1666 young stellar objects (YSOs) candidates (including 154 CTTSs and 419 HAeBe stars) are selected. In the interacting regions, the significant enhancement of the number of protostellar objects and YSOs indicates the presence of some recently formed stars. After comparing the characteristic timescales of star formation with the age of IC443, we conclude that the protostellar objects and YSOs candidates are not triggered by IC443. For the age of stellar winds shell, we have performed calculation on the basis of stellar winds shell expansion model. The results and analysis suggest that the formation of these stars may be triggered by the stellar winds of IC443 progenitor.

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.

A number of neutron star low-mass X-ray binaries have recently been discovered to show broad, asymmetric Fe K emission lines in their X-ray spectra. These lines are generally thought to be the most prominent part of a reflection spectrum, originating in the inner part of the accretion disk where strong relativistic effects can broaden emission lines. We present a comprehensive, systematic analysis of Suzaku and XMM-Newton spectra of 10 neutron star low-mass X-ray binaries, all of which display broad Fe K emission lines. Of the 10 sources, 4 are Z sources, 4 are atolls and 2 are accreting millisecond X-ray pulsars (also atolls). The Fe K lines are well fit by a relativistic line model for a Schwarzschild metric, and imply a narrow range of inner disk radii (6 - 15 GM/c^2) in most cases. This implies that the accretion disk extends close to the neutron star surface over a range of luminosities. Continuum modeling shows that for the majority of observations, a blackbody component (plausibly associated with the boundary layer) dominates the X-ray emission from 8 - 20 keV. Thus it appears likely that this spectral component produces the majority of the ionizing flux that illuminates the accretion disk. Therefore, we also fit the spectra with a blurred reflection model, wherein a blackbody component illuminates the disk. This model fits well in most cases, supporting the idea that the boundary layer is illuminating a geometrically thin disk.

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.

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.

We report on several features present in the energy spectrum from an ultra low-noise germanium detector operated at 2,100 m.w.e. By implementing a new technique able to reject surface events, a number of cosmogenic peaks can be observed for the first time. We discuss several possible causes for an irreducible excess of bulk-like events below 3 keVee, including a dark matter candidate common to the DAMA/LIBRA annual modulation effect, the hint of a signal in CDMS, and phenomenological predictions. Improved constraints are placed on a cosmological origin for the DAMA/LIBRA effect.

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope to detect mergers of intermediate-mass black holes that may have formed through runaway stellar collisions in globular clusters. We find that detection rates of 500 events per year are plausible.

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.

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.

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.

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.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas to move. We show that this "nearly aligned moving spot model" can explain many properties of the accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some are intermittent, and why all detected so far are transients. This model can be tested by comparing with observations the waveform changes it predicts, including the changes with accretion rate.

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.

GRS 1915+105 harbors one of the most massive known stellar black holes in the Galaxy. In May 2007, we observed GRS 1915+105 for 117 ksec in the low/hard state using Suzaku. We collected and analyzed the data with the HXD/PIN and XIS cameras spanning the energy range from 2.3-55 keV. Fits to the spectra with simple models reveal strong disk reflection through an Fe K emission line and a Compton back-scattering hump. We report constraints on the spin parameter of the black hole in GRS 1915+105 using relativistic disk reflection models. The model for the soft X-ray spectrum (i.e. < 10 keV) suggests a/M = 0.56(2) and excludes zero spin at the 4 sigma level of confidence. The model for the full broadband spectrum suggests that the spin may be higher, a/M = 0.98(1) (1 sigma confidence), and again excludes zero spin at the 2 sigma level of confidence. We discuss these results in the context of other spin constraints and inner disk studies in GRS 1915+105.

We search for evidence of a relation between properties of young stellar objects (YSOs) and their parent molecular clouds to understand the initial conditions of high-mass star formation. A sample of 135 sources was selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog, on the basis of their red color to enhance the possibility of discovering young sources. Using the Kolner Observatorium fur SubMillimeter Astronomie (KOSMA) 3-m telescope, a single-point survey in 13CO(2-1) was carried out for the entire sample, and 14 sources were mapped further. Archival mid-infrared (MIR) data were compared with the 13CO emissions to identify evolutionary stages of the sources. A 13CO observed sample was assembled to investigate the correlation between 13CO line width of the clouds and the luminosity of the associated YSOs. We identified 98 sources suitable for star formation analyses for which relevant parameters were calculated. We detected 18 cores from 14 mapped sources, which were identified with eight pre-UC HII regions and one UC HII region, two high-mass cores earlier than pre-UC HII phase, four possible star forming clusters, and three sourceless cores. By compiling a large (360 sources) 13CO observed sample, a good correlation was found between the 13CO line width of the clouds and the bolometric luminosity of the associated YSOs, which can be fitted as a power law: lg(dV13/km/s)=-0.023+0.135lg(Lbol/Lsolar). Results show that luminous (>10^3Lsolar) YSOs tend to be associated with both more massive and more turbulent (dV13>2km/s) molecular cloud structures.

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.

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.

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.

The association of an electromagnetic signal with the merger of a pair of supermassive black holes would have many important implications. For example, it would provide new information about gas and magnetic field interactions in dynamical spacetimes as well as a combination of redshift and luminosity distance that would enable precise cosmological tests. A proposal first made by Bode & Phinney (2007) is that because radiation of gravitational waves during the final inspiral and merger of the holes is abrupt and decreases the mass of the central object by a few percent, there will be waves in the disk that can steepen into shocks and thus increase the disk luminosity in a characteristic way. We evaluate this process analytically and numerically. We find that shocks only occur when the fractional mass loss exceeds the half-thickness (h/r) of the disk, hence significant energy release only occurs for geometrically thin disks which are thus at low Eddington ratios. This strongly limits the effective energy release, and in fact our simulations show that the natural variations in disk luminosity are likely to obscure this effect entirely. However, we demonstrate that the reduction of luminosity caused by the retreat of the inner edge of the disk following mass loss is potentially detectable. This decrease occurs even if the disk is geometrically thick, and lasts for a duration on the order of the viscous time of the modified disk. Observationally, the best prospect for detection would be a sensitive future X-ray instrument with a field of view of on the order of a square degree, or possibly a wide-field radio array such as the Square Kilometer Array, if the disk changes produce or interrupt radio emission from a jet.

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

Electromagnetic observations over the last 15 years have yielded a growing appreciation for the importance of supermassive black holes (SMBH) to the evolution of galaxies, and for the intricacies of dynamical interactions in our own Galactic center. Here we show that future low-frequency gravitational wave observations, alone or in combination with electromagnetic data, will open up unique windows to these processes. In particular, gravitational wave detections in the 10^{-5}-10^{-1} Hz range will yield SMBH masses and spins to unprecedented precision and will provide clues to the properties of the otherwise undetectable stellar remnants expected to populate the centers of galaxies. Such observations are therefore keys to understanding the interplay between SMBHs and their environments.

Intermediate-mass black holes (IMBHs), with masses of hundreds to thousands of solar masses, will be unique sources of gravitational waves for LISA. Here we discuss their context as well as specific characteristics of IMBH-IMBH and IMBH-supermassive black hole mergers and how these would allow sensitive tests of the predictions of general relativity in strong gravity.

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.

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.

Context: The earliest phases of massive star formation are currently much debated. Aims. In an effort to make progress, we took a census of Class0-like protostellar dense cores in the NGC 3576 region, one of the nearest and most luminous embedded sites of high-mass star formation in the Galaxy. Methods: We used the P-ArTeMiS bolometer camera on the APEX telescope to produce the first 450-micron dust continuum map of the filamentary dense clump associated with NGC 3576. Results: Combining our 450-micron observations with existing data at other wavelengths, we have identified seven massive protostellar sources along the NGC 3576 filament and placed them in the M_env - L_bol evolutionary diagram for protostars. Conclusions: Comparison with theoretical evolutionary tracks suggests that these seven protostellar sources will evolve into massive stars with masses M* ~ 15-50 Msun. Four sources are classified as candidate high-mass Class 0 objects, two sources as massive Class I objects, and one source appears to be at an intermediate stage.

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\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^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.

XEUS has been recently selected by ESA for an assessment study. XEUS is a large mission candidate for the Cosmic Vision program, aiming for a launch date as early as 2018. XEUS is a follow-on to ESA's Cornerstone X-Ray Spectroscopy Mission (XMM-Newton). It will be placed in a halo orbit at L2, by a single Ariane 5 ECA, and comprises two spacecrafts. The Silicon pore optics assembly of XEUS is contained in the mirror spacecraft while the focal plane instruments are contained in the detector spacecraft, which is maintained at the focus of the mirror by formation flying. The main requirements for XEUS are to provide a focused beam of X-rays with an effective aperture of 5 m^2 at 1 keV, 2 m^2 at 7 keV, a spatial resolution better than 5 arcsec, a spectral resolution ranging from 2 to 6 eV in the 0.1-8 keV energy band, a total energy bandpass of 0.1-40 keV, ultra-fast timing, and finally polarimetric capabilities. The High Time Resolution Spectrometer (HTRS) is one of the five focal plane instruments, which comprises also a wide field imager, a hard X-ray imager, a cryogenic spectrometer, and a polarimeter. The HTRS is unique in its ability to cope with extremely high count rates (up to 2 Mcts/s), while providing sub-millisecond time resolution and good (CCD like) energy resolution. In this paper, we focus on the specific scientific objectives to be pursued with the HTRS: they are all centered around the key theme "Matter under extreme conditions" of the Cosmic Vision science program. We demonstrate the potential of the HTRS observations to probe strong gravity and matter at supra-nuclear densities. We conclude this paper by describing the current implementation of the HTRS in the XEUS focal plane.

In many theoretical scenarios it is expected that intermediate-mass black holes (IMBHs, with masses M ~ 100-10000 solar masses) reside at the centers of some globular clusters. However, observational evidence for their existence is limited. Several previous numerical investigations have focused on the impact of an IMBH on the cluster dynamics or brightness profile. Here we instead present results from a large set of direct N-body simulations including single and binary stars. These show that there is a potentially more detectable IMBH signature, namely on the variation of the average stellar mass between the center and the half-light radius. We find that the existence of an IMBH quenches mass segregation and causes the average mass to exhibit only modest radial variation in collisionally relaxed star clusters. This differs from when there is no IMBH. To measure this observationally requires high resolution imaging at the level of that already available from the Hubble Space Telescope (HST) for the cores of a large sample of galactic globular clusters. With a modest additional investment of HST time to acquire fields around the half-light radius, it will be possible to identify the best candidate clusters to harbor an IMBH. This test can be applied only to globulars with a half-light relaxation time less than or equal to 1 Gyr, which is required to guarantee efficient energy equipartition due to two-body relaxation.

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.

Mergers between stellar-mass black holes will be key sources of gravitational radiation for ground-based detectors. However, the rates of these events are highly uncertain, given that such systems are invisible. One formation scenario involves mergers in field binaries, where our lack of complete understanding of common envelopes and the distribution of supernova kicks has led to rate estimates that range over a factor of several hundred. A different, and highly promising, channel involves multiple encounters of binaries in globular clusters or young star clusters. However, we currently lack solid evidence for black holes in almost all such clusters, and their low escape speeds raise the possibility that most are ejected because of supernova recoil. Here we propose that a robust environment for mergers could be the nuclear star clusters found in the centers of small galaxies. These clusters have millions of stars, black hole relaxation times well under a Hubble time, and escape speeds that are several times those of globulars, hence they retain most of their black holes. We present simulations of the three-body dynamics of black holes in this environment and estimate that, if most nuclear star clusters do not have supermassive black holes that interfere with the mergers, at least several tens of events per year will be detectable with Advanced LIGO.

Aims: We searched for the molecular outflows from fifteen molecular clouds associated with ultra-compact and compact HII (UCHII/HII) regions and discussed possible gas heating mechanism. Methods: Mapping observations of CO $J=3-2$ and $J=2-1$ lines were carried out with the KOSMA 3m-Telescope towards the 15 HII regions/molecular cloud complexes. Results: Ten molecular outflows were identified out of the fifteen HII region/molecular cloud complexes. The higher outflow detection rate (67%) suggested that such outflows are as common in high mass star forming regions as those in low mass star forming regions, which is consistent with the results of other authors. The observations also showed that the outflow might occur in the HII region. The integrated CO line intensity ratios (${\rm R_{I_{CO(3-2)}/I_{CO(2-1)}}}$) were determined from the core component of the spectra as well as from both the blue and red wings. Maximum line intensity ratios from the wings and core components appeared to be related to the mid-infrared sources imaged by Midcourse Space Experiment (MSX). The relationship between the maximum line intensity ratios and MSX sources indicates that the molecular gas could be heated by the emission of dust associated with massive stars. Based on maser observations reported in the literature, we found that ${\rm H_{2}O}$ masers were only detected in seven regions. The ${\rm H_{2}O}$ masers in these regions are located near the MSX sources and within the maximum line intensity ratio regions, suggesting that ${\rm H_{2}O}$ masers occur in relatively warm environments.

Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud (LMC) is ideally suited to study star formation in such an environment. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR spectral lines of the main carbon carrying species, CO, CI and CII, which originate in the surface layers of molecular clouds illuminated by the UV radiation of the newly formed, young stars. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 telescope of the 12CO J = 4-3, J = 7-6, and 13CO J = 4-3 rotational and [CI] 3P1-3P0 and 3P2-3P1 fine-structure transitions. The 13CO J =4-3 and [CI] 3P2-3P1 transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W region has temperatures of about 80 K and densities of about 10^4 cm^-3. The estimated C to CO abundance ratio close to unity is substantially higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [CII] line intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about \chi ~220 and an average density of the clump ensemble of about 10^5 cm^-3, thus confirming the presence of high density material in the LMC N159W region.

Studying molecular gas in the central regions of the star burst galaxies NGC4945 and Circinus enables us to characterize the physical conditions and compare them to previous local and high-z studies. We estimate temperature, molecular density and column densities of CO and atomic carbon. Using model predictions we give a range of estimated CO/C abundance ratios. Using the new NANTEN2 4m sub-millimeter telescope in Pampa La Bola, Chile, we observed for the first time CO 4-3 and [CI] 3P1-3 P0 at the centers of both galaxies at linear scale of 682 pc and 732 pc respectively. We compute the cooling curves of 12CO and 13CO using radiative transfer models and estimate the physical conditions of CO and [CI]. The centers of NGC4945 and Circinus are very [CI] bright objects, exhibiting [CI] 3P1 - 3 P0 luminosities of 91 and 67Kkms-1kpc2, respectively. The [CI] 3P1-3 P0/CO 4-3 ratio of integrated intensities are large at 1.2 in NGC4945 and 2.8 in Circinus. Combining previous CO J= 1-0, 2-1 and 3-2 and 13CO J= 1-0, 2-1 studies with our new observations, the radiative transfer calculations give a range of densities, n(H2) = 10^3-3*104^cm-3, and a wide range of kinetic temperatures, Tkin = 20 - 100K, depending on the density. Future CO J= 7-6 and [CI] 2-1 observations will be important to resolve the ambiguity in the physical conditions and confirm the model predictions.

We report on a comprehensive analysis of the kilo-Hz (>~600 Hz) quasi-periodic oscillations (kHz QPOs) detected from the neutron star X-ray transient Aquila X-1 (Aql X-1) with the Rossi X-ray Timing Explorer, between 1997 and 2007. Among kHz QPO sources, Aql X-1 is peculiar because so far only one kHz QPO has been reported, whereas in most sources, two kHz QPOs are usually detected (a lower and an upper kHz QPO). The identification of the QPOs reported so far has therefore been ambiguous, although it has been proposed that they were likely to be the lower QPO. Following up on previous work, we confirm the identification of the QPOs previously reported as lower QPOs, because of their high quality factors and the quality factor versus frequency dependency, which are similar to those observed in other sources. Combining all segments of data containing a lower QPO, we detect for the first time an upper kHz QPO. As in other sources for which the neutron star spin frequency is larger than 400 Hz (550.25 Hz in Aql X-1), the frequency difference between the two kHz QPOs is close to half the spin frequency. Based on this result, we re-examine the link between the neutron star spin and the frequency of the kHz QPOs, to show that a model in which the separation of the lower and upper QPOs relates to the neutron star spin frequency is still as good as any comparably simple model.

The Carina region is an excellent astrophysical laboratory for studying the feedback mechanisms of newly born, very massive stars within their natal giant molecular clouds (GMCs) at only 2.35 kpc distance. We use a clumpy PDR model to analyse the observed intensities of atomic carbon and CO and to derive the excitation conditions of the gas. The NANTEN2-4m submillimeter telescope was used to map the [CI] 3P1-3P0, 3P2-3P1 and CO 4-3, 7-6 lines in two 4'x4' regions of Carina where molecular material interfaces with radiation from the massive star clusters. One region is the northern molecular cloud near the compact OB cluster Tr14, and the second region is in the molecular cloud south of etaCar and Tr16. These data were combined with 13CO SEST spectra, HIRES/IRAS 60um and 100um maps of the FIR continuum, and maps of 8um IRAC/Spitzer and MSX emission. We used the HIRES far-infrared dust data to create a map of the FUV field heating the gas. The northern region shows an FUV field of a few 1000 in Draine units while the field of the southern region is about a factor 10 weaker. We constructed models consisting of an ensemble of small spherically symmetric PDR clumps within the 38" beam (0.43pc), which follow canonical power-law mass and mass-size distributions. We find that an average local clump density of 2x10**5 cm-3 is needed to reproduce the observed line emission at two selected interface positions. Stationary, clumpy PDR models reproduce the observed cooling lines of atomic carbon and CO at two positions in the Carina Nebula.

Using Suzaku observations of three neutron star low-mass X-ray binaries (Ser X-1, 4U 1820-30 and GX 349+2) we have found broad, asymmetric, relativistic Fe K emission lines in all three objects. These Fe K lines can be well fit by a model for lines from a relativistic accretion disk ('diskline'), allowing a measurement of the inner radius of the accretion disk, and hence an upper limit on the neutron star radius. These upper limits correspond to 14.5 - 16.5 km for a 1.4 M(solar) neutron star. The inner disk radii we measure with Fe K lines are in good agreement with the inner disk radii implied by kHz QPOs observed in both 4U 1820-30 and GX 349+2, supporting the inner disk nature of kHz QPOs. Additionally, the Fe K lines observed in these neutron stars are narrower than those in the black holes that are thought to be close to maximally spinning, as one would expect if inferences for spin are robust.

Using the IRAM 30 m telescope, a mapping survey in optically thick and thin lines was performed towards 46 high mass star-forming regions. The sample includes UC H{\sc ii} precursors and UC H{\sc ii} regions. Seventeen sources are found to show "blue profiles", the expected signature of collapsing cores. The excess of sources with blue over red profiles ([$N_{\rm blue}$ -- $N_{\rm red}$]/$N_{\rm total}$) is 29% in the HCO$^+$ $J$=1--0 line, with a probability of 0.6% that this is caused by random fluctuations. UC H{\sc ii} regions show a higher excess (58%) than UC H{\sc ii} precursors (17%), indicating that material is still accreted after the onset of the UC H{\sc ii} phase. Similar differences in the excess of blue profiles as a function of evolutionary state are not observed in low mass star-forming regions. Thus, if confirmed for high mass star-forming sites, this would point at a fundamental difference between low- and high-mass star formation. Possible explanations are inadequate thermalization, stronger influence of outflows in massive early cores, larger gas reserves around massive stellar objects or different trigger mechanisms between low- and high- mass star formation.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

We investigate the structure of a field around the position of V838 Mon as seen in the lowest CO rotational transitions. We also measure and analyse emission in the same lines at the position of V838 Mon.Observations have primarily been done in the CO J = 2-1 and J = 3-2 lines using the KOSMA telescope. A field of 3.4 squared degrees has been mapped in the on-the-fly mode in these transitions. Longer integration spectra in the on-off mode have been obtained to study the emission at the position of V838 Mon. Selected positions in the field have also been observed in the CO J = 1-0 transition using the Delingha telescope.In the observed field we have identified many molecular clouds. They can be divided into two groups from the point of view of their observed radial velocities. One, having V(LSR) in the range 18-32 km/s, can be identified with the Perseus Galactic arm. The other one, having V(LSR) between 44-57 km/s, probably belongs to the Norma-Cygnus arm. The radial velocity of V838 Mon is within the second range but the object does not seem to be related to any of the observed clouds. We did not find any molecular buble of a 1 degree dimension around the position of V838 Mon claimed in van Loon et al. An emission has been detected at the position of the object in the CO J = 2-1 and J = 3-2 transitions. The emission is very narrow (FWHM ~ 1.2 km/s) and at V(LSR) = 53.3 km/s. Our analysis of the data suggests that the emission is probably extended.

Black hole binaries with extreme ($\gtrsim 10^4:1$) or intermediate ($\sim 10^2-10^4:1$) mass ratios are among the most interesting gravitational wave sources that are expected to be detected by the proposed Laser Interferometer Space Antenna. These sources have the potential to tell us much about astrophysics, but are also of unique importance for testing aspects of the general theory of relativity in the strong field regime. Here we discuss these sources from the perspectives of astrophysics, data analysis, and applications to testing general relativity, providing both a description of the current state of knowledge and an outline of some of the outstanding questions that still need to be addressed. This review grew out of discussions at a workshop in September 2006 hosted by the Albert Einstein Institute in Golm, Germany.

Recent numerical relativistic simulations of black hole coalescence suggest that in certain alignments the emission of gravitational radiation can produce a kick of several thousand kilometers per second. This exceeds galactic escape speeds, hence unless there a mechanism to prevent this, one would expect many galaxies that had merged to be without a central black hole. Here we show that in most galactic mergers, torques from accreting gas suffice to align the orbit and spins of both black holes with the large-scale gas flow. Such a configuration has a maximum kick speed <200 km/s, safely below galactic escape speeds. We predict, however, that in mergers of galaxies without much gas, the remnant will be kicked out several percent of the time. We also discuss other predictions of our scenario, including implications for jet alignment angles and X-type radio sources.

Analysis of archival RXTE data on neutron stars binaries has shown that for several sources the quality factor (Q) of the lower kilohertz Quasi-Periodic Oscillations (QPO) drops sharply beyond a certain frequency. This is one possible signature of the approach to the general relativistic innermost stable circular orbit (ISCO), but the implications of such an interpretation for strong gravity and dense matter are important enough that it is essential to explore alternate explanations. In this spirit, Mendez has recently proposed that Q depends fundamentally on mass accretion rate (as measured by spectral hardness) rather than the frequency of the QPO. We test this hypothesis for 4U1636-536 by measuring precisely spectral colors simultaneously with the lower QPO frequency and Q after correction for the frequency drift, over a data set spanning eight years of RXTE observations. We find that in this source there is no correlation between Q and spectral hardness. In particular, no apparent changes in hardness are observed when Q reaches its maximum before dropping off. We perform a similar analysis on 4U1608-522; another source showing a sharp drop in the quality factor of its lower kHz QPO. We find that for this source, positive and negative correlations are observed between spectral hardness, frequency and Q. Consequently, if we are to search for a common explanation for the sharp drop in the quality factor seen in both sources, the spectral hardness is not a good candidate for the independent variable whereas the frequency remains. Therefore, we conclude that the ISCO explanation is viable for 4U1636-536, and thus possibly for others.

We investigate the dynamic stability of inspiraling neutron stars by performing multiple-orbit numerical relativity simulations of the binary neutron star inspiral process. By introducing eccentricities in the orbits of the neutron stars, significant changes in orbital separation are obtained within orbital timescales. We find that as the binary system evolves from apastron to periastron (as the binary separation decreases), the central rest mass density of each star decreases, thus stabilizing the stars against individual prompt collapse. As the binary system evolves from periastron to apastron, the central rest mass density increases; the neutron stars re-compress as the binary separation increases.

With current terrestrial gravitational wave detectors working at initial design sensitivities, and upgrades and space missions planned, it is likely that in the next five to ten years gravitational radiation will be detected directly from a variety of classes of objects. The most confidently expected of these classes is compact binaries, involving neutron stars or black holes. Detection of their coalescence, or their long-term orbits, has the potential to inform us about the evolutionary history of compact binaries and possibly even star formation over the past several billion years. We review what is currently known about compact binaries as sources of gravitational radiation, as well as the current uncertainties and what we expect to learn from future detections of gravitational waves from these systems.

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86--116 km s$^{-1}$, and the most reliable runs give kicks between 86 and 97 km s$^{-1}$. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts $z\gtrsim 10$, halos with masses $\lesssim 10^9 M_\odot$ will have difficulty retaining coalesced black holes after major mergers.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

Motivated by the report by Cottam et al. (2002) of iron resonance scattering lines in the spectra of thermonuclear bursts from EXO 0748-676, we have investigated the information about neutron star structure and the geometry of the emission region that can be obtained by analyzing the profiles of atomic lines formed at the surface of the star. We have calculated the detailed profiles of such lines, taking into account the star's spin and the full effects of special and general relativity, including light-bending and frame-dragging. We discuss the line shapes produced by rotational Doppler broadening and magnetic splitting of atomic lines, for the spin rates and magnetic fields expected in neutron stars in low-mass X-ray binary systems. We show that narrow lines are possible even for rapidly spinning stars, if the emission region or the line of sight are close to the spin axis. For most neutron stars in low-mass systems, magnetic splitting is too small to obscure the effects of special and general relativity. We show that the ratio of the star's mass to its equatorial radius can be determined to within 5% using atomic line profiles, even if the lines are broad and skewed. This is the precision required to constrain strongly the equation of state of neutron star matter. We show further that if the radius and latitude of emission are known to ~ 5-10% accuracy, then frame-dragging has a potentially detectable effect on the profiles of atomic lines formed at the stellar surface.

Recent numerical simulations have suggested the possibility of forming double intermediate mass black holes (IMBHs) via the collisional runaway scenario in young dense star clusters. The two IMBHs formed would exchange into a common binary shortly after their birth, and quickly inspiral and merge. Since space-borne gravitational wave (GW) observatories such as LISA will be able to see the late phases of their inspiral out to several Gpc, and LIGO will be able to see the merger and ringdown out to similar distances, they represent potentially significant GW sources. In this Letter we estimate the rate at which LISA and LIGO will see their inspiral and merger in young star clusters, and discuss the information that can be extracted from the observations. We find that LISA will likely see tens of IMBH--IMBH inspirals per year, while advanced LIGO could see ~10 merger and ringdown events per year, with both rates strongly dependent on the distribution of cluster masses and densities.

We present observations of V838 Mon and its close vicinity in the three lowest rotational transitions of CO. The J = 2-1 and 3-2 data were obtained using the 3 m KOSMA telescope. They include on-the-fly maps covering a large area (~3.4 sq. deg.) around V838 Mon and long integrations on the star position. Complementary observations in the CO J = 1-0 transition were obtained using the 13.7 m Delingha telescope. The star position as well as 25 other points preselected in the near vicinity of the object have been measured in this transition. We report on a detection of two narrow emission components in J = 2-1 and 3-2 transitions at the position of V838 Mon. Lines were found at radial velocities of V_lsr=53.3 km/s and V_lsr=-11.0 km/s. Their origin is unclear. We also shortly discuss results of the observations of the vicinity of V838 Mon.

We study in a systematic way the quality factor of the lower and upper kHz QPOs in a sample of low luminosity neutron star X-ray binaries, showing both QPOs varying over a wide frequency range. The sample includes 4U 1636-536, 4U 1608-522, 4U 1735-44, 4U 1728-34, 4U 1820-303 and 4U 0614+09. We find that all sources except 4U 0614+091 show evidence of a drop in the quality factor of their lower kHz QPOs at high frequency. For 4U 0614+091 only the rising part of the quality factor versus frequency curve has been sampled so far. At the same time, in all sources but 4U 1728-34, the quality factor of the upper kilo-Hz QPO increases all the way to the highest detectable frequencies. We show that the high-frequency behaviours of both the lower and upper kHz QPO quality factors are consistent with what is expected if the drop is produced by the approach of an active oscillating region to the innermost stable circular orbit: the existence of which is a key feature of General Relativity in the strong field regime. Within this interpretation, our results imply gravitational masses around 2 solar masses for the neutron stars in those systems.

A new scheme for testing nuclear matter equations of state (EsoS) at high densities using constraints from neutron star phenomenology and a flow data analysis of heavy-ion collisions is suggested. An acceptable EoS shall not allow the direct Urca process to occur in neutron stars with masses below $1.5~M_{\odot}$, and also shall not contradict flow and kaon production data of heavy-ion collisions. Compact star constraints include the mass measurements of 2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807, of 2.0 +/- 0.1 M_sun from the innermost stable circular orbit for 4U 1636-536, the baryon mass - gravitational mass relationships from Pulsar B in J0737-3039 and the mass-radius relationships from quasiperiodic brightness oscillations in 4U 0614+09 and from the thermal emission of RX J1856-3754. This scheme is applied to a set of relativistic EsoS constrained otherwise from nuclear matter saturation properties with the result that no EoS can satisfy all constraints simultaneously, but those with density-dependent masses and coupling constants appear most promising.

Characterizing the spatial and velocity structure of molecular clouds is a first step towards a better understanding of interstellar turbulence and its link to star formation. We present observations and structure analysis results for a large-scale (~ 7.10 deg^2) 13CO J = 2--1 and 12 CO J = 3--2 survey towards the nearby Perseus molecular cloud observed with the KOSMA 3m telescope. We study the spatial structure of line-integrated and velocity channel maps, measuring the $\Delta$-variance as a function of size scale. We determine the spectral index $\beta$ of the corresponding power spectrum and study its variation across the cloud and across the lines. We find that the spectra of all CO line-integrated maps of the whole complex show the same index, $\beta$ ~ 3.1, for scales between about 0.2 and 3pc, independent of isotopomer and rotational transition. A complementary 2MASS map of optical extinction shows a noticeably smaller index of 2.6. In contrast to the overall region, the CO maps of individual subregions show a significant variation of $\beta$. The 12CO 3--2 data provide e.g. a spread of indices between 2.9 in L1455 and 3.5 in NGC1333. In general, active star forming regions show a larger power-law exponent. We find that the $\Delta$-variance spectra of individual velocity channel maps are very sensitive to optical depth effects clearly indicating self-absorption in the densest regions. When studying the dependence of the channel-map spectra as a function of the velocity channel width, the expected systematic increase of the spectral index with channel width is only detected in the blue line wings.This could be explained by a filamentary, pillar-like structure which is left at low velocities while the overall molecular gas is swept up by a supernova shock wave.

In this letter, we comment on the robustness of putative cool (kT ~ 0.2 keV) accretion disc components in the X-ray spectra of the most luminous (L ~ 10^40 erg/s) ultra-luminous X-ray sources (ULXs) in nearby normal galaxies. When compared to stellar-mass black holes, the low disc temperatures observed in some ULXs may imply intermediate-mass black hole primaries. Recent work has claimed that such soft excesses are unlikely to be actual disc components, based on the lack of variability in these components, and in the overall source flux. Other work has proposed that alternative phenomenological models, and complex Comptonisation models, rule-out cool disc components in ULX spectra. An inspection of the literature on Galactic stellar-mass black holes and black hole candidates demonstrates that the flux behaviours seen in specific ULXs are consistent with phenomena observed in well-known Galactic X-ray binaries. Applying Comptonisation models to simulated disc blackbody plus power-law spectra shows that at the sensitivity achieved in even the best ULX spectra, Comptonisation fits are highly model-dependent, and do not yield meaningful constraints on the accretion flow. In contrast, the need for a soft, thermal component does not appear to be model-dependent. As we have previously noted, soft thermal components in ULX spectra may not represent accretion discs, but present alternatives to this interpretation are not robust.

We discuss millisecond period brightness oscillations and surface atomic spectral lines observed during type I X-ray bursts from a neutron star in a low mass X-ray binary system. We show that modeling of these phenomena can constrain models of the dense cold matter at the cores of neutron stars. We demonstrate that, even for a broad and asymmetric spectral line, the stellar radius-to-mass ratio can be inferred to better than 5%. We also fit our theoretical models to the burst oscillation data of the low mass X-ray binary XTE J1814-338, and find that the 90% confidence lower limit of the neutron star's dimensionless radius-to-mass ratio is 4.2.

Observations of the centers of galaxies continue to evolve, and it is useful to take a fresh look at the constraints that exist on alternatives to supermassive black holes at their centers. We discuss constraints complementary to those of Maoz (1998) and demonstrate that an extremely wide range of other possibilities can be excluded. In particular, we present the new argument that for the velocity dispersions inferred for many galactic nuclei, even binaries made of point masses cannot stave off core collapse because hard binaries are so tight that they merge via emission of gravitational radiation before they can engage in three-body or four-body interactions. We also show that under these conditions core collapse leads inevitably to runaway growth of a central black hole with a significant fraction of the initial mass, regardless of the masses of the individual stars. For clusters of noninteracting low-mass objects (from low-mass stars to elementary particles), relaxation of stars and compact objects that pass inside the dark region will be accelerated by interactions with the dark mass. If the dark region is instead a self-supported object such as a fermion ball, then if stellar-mass black holes exist they will collide with the object, settle, and consume it. The net result is that the keyhole through which alternatives to supermassive black holes must pass is substantially smaller and more contrived than it was even a few years ago.

We present numerical three-body experiments that include the effects of gravitational radiation reaction by using equations of motion that include the 2.5-order post-Newtonian force terms, which are the leading order terms of energy loss from gravitational waves. We simulate binary-single interactions and show that close approach cross sections for three 1 solar mass objects are unchanged from the purely Newtonian dynamics except for close approaches smaller than 1.0e-5 times the initial semimajor axis of the binary. We also present cross sections for mergers resulting from gravitational radiation during three-body encounters for a range of binary semimajor axes and mass ratios including those of interest for intermediate-mass black holes (IMBHs). Building on previous work, we simulate sequences of high-mass-ratio three-body encounters that include the effects of gravitational radiation. The simulations show that the binaries merge with extremely high eccentricity such that when the gravitational waves are detectable by LISA, most of the binaries will have eccentricities e > 0.9 though all will have circularized by the time they are detectable by LIGO. We also investigate the implications for the formation and growth of IMBHs and find that the inclusion of gravitational waves during the encounter results in roughly half as many black holes ejected from the host cluster for each black hole accreted onto the growing IMBH.

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.

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.

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

Mergers of neutron stars with black holes have been suggested as candidates for short gamma-ray bursts. They have also been studied for their potential as gravitational wave sources observable with ground-based detectors. For these purposes, it is important to know under what circumstances such a merger could leave an accretion disk or result in a period of stable mass transfer. We show that, consistent with recent numerical simulations, it is expected that mergers between neutron stars and black holes will be prompt, with no accretion disk and no stable mass transfer, if the black hole has a mass greater than that of the neutron star and is spinning slowly. The reason is that for comparable masses, angular momentum loss to gravitational radiation starts a plunge orbit well outside the innermost stable circular orbit, causing direct merging rather than extended mass transfer. Even when the black hole is spinning rapidly and exactly prograde with respect to the orbit, we show that it is possible within current understanding that no accretion disk will form under any circumstances, but resolution of this will require full general relativistic numerical simulations with no approximations.

I discuss the accuracy requirements on numerical relativity calculations of inspiraling compact object binaries whose extracted gravitational waveforms are to be used as templates for matched filtering signal extraction and physical parameter estimation in modern interferometric gravitational wave detectors. Using a post-Newtonian point particle model for the pre-merger phase of the binary inspiral, I calculate the maximum allowable errors for the mass and relative velocity and positions of the binary during numerical simulations of the binary inspiral. These maximum allowable errors are compared to the errors of state-of-the-art numerical simulations of multiple-orbit binary neutron star calculations in full general relativity, and are found to be smaller by several orders of magnitude. A post-Newtonian model for the error of these numerical simulations suggests that adaptive mesh refinement coupled with second order accurate finite difference codes will {\it not} be able to robustly obtain the accuracy required for reliable gravitational wave extraction on Terabyte-scale computers. I conclude that higher order methods (higher order finite difference methods and/or spectral methods) combined with adaptive mesh refinement and/or multipatch technology will be needed for robustly accurate gravitational wave extraction from numerical relativity calculations of binary coalescence scenarios.

Detailed modeling of the millisecond brightness oscillations during thermonuclear bursts from low mass X-ray binaries can provide important information about neutron star structure. Until now the implementation of this idea has not been entirely successful, largely because of the negligible harmonic content in burst oscillation lightcurves. However, the recent discovery of non-sinusoidal burst oscillation lightcurves from the accreting millisecond pulsar XTE J1814-338 has changed this situation. We, therefore, for the first time, make use of this opportunity to constrain neutron star parameters. In our detailed study of the lightcurves of 22 bursts, we fit the burst oscillation lightcurves with fully general relativistic models that include light-bending and frame-dragging for lightcurve calculation, and compute numerically the structure of neutron stars using realistic equations of state. We find that for our model and parameter grid values, at the 90% confidence level, Rc^2/GM > 4.2 for the neutron star in XTE J1814-338. We also find that the photons from the thermonuclear flash come out through the layers of accreted matter under conditions consistent with Thomson scattering, and show that the secondary companion is a hydrogen burning main sequence star, with possible bloating (probably due to X-ray heating).

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.

Cool thermal emission components have recently been revealed in the X-ray spectra of a small number of ultra-luminous X-ray (ULX) sources with L_X > 1 E+40 erg/s in nearby galaxies. These components can be well fitted with accretion disk models, with temperatures approximately 5-10 times lower than disk temperatures measured in stellar-mass Galactic black holes when observed in their brightest states. Because disk temperature is expected to fall with increasing black hole mass, and because the X-ray luminosity of these sources exceeds the Eddington limit for 10 Msun black holes (L_Edd = 1.3 E+39 erg/s), these sources are extremely promising intermediate-mass black hole candidates (IMBHCs). In this Letter, we directly compare the inferred disk temperatures and luminosities of these ULXs, with the disk temperatures and luminosities of a number of Galactic black holes. The sample of stellar-mass black holes was selected to include different orbital periods, companion types, inclinations, and column densities. These ULXs and stellar-mass black holes occupy distinct regions of a L_X -- kT diagram, suggesting these ULXs may harbor IMBHs. We briefly discuss the important strengths and weaknesses of this interpretation.

We report the results of an analysis of two XMM-Newton/EPIC-pn spectra of the bright ultraluminous X-ray source M81 X-9 (Holmberg IX X-1), obtained in snapshot observations. Soft thermal emission is clearly revealed in spectra dominated by hard power-law components. Depending on the model used, M81 X-9 was observed at a luminosity of L_X = 1.0-1.6 E+40 erg/s (0.3--10.0 keV). The variability previously observed in this source signals that it is an accreting source which likely harbors a black hole. Remarkably, accretion disk models for the soft thermal emission yield very low inner disk temperatures (kT = 0.17-0.29 keV, including 90% confidence errors and variations between observations and disk models), and improve the fit statistic over any single-component continuum model at the 6 sigma level of confidence. Like NGC 1313 X-1, scaling the temperatures measured in M81 X-9 to those commonly seen in stellar-mass Galactic black holes at their highest observed fluxes (kT ~ 1 keV) may imply that M81 X-9 harbors a black hole with a mass on the order of 10^3 Msun; the measured disk component normalization and broad-band luminosity imply black hole masses on the order of 10^2 Msun. It is therefore possible that these sources harbor 10^3 Msun black holes accreting such at L_X ~ 0.1*L_Edd. Radio and optical observations suggest that beaming and anisotropic emission from a stellar-mass black hole are unlikely to account for the implied luminosity. We further argue that the strength of the hard emission in these sources and well-established phenomena frequently observed in stellar-mass black holes near to the Eddington limit suggest that optically-thick photospheres are unlikely to be the origin of the cool thermal emission in bright ULXs.

Evidence is mounting that some Ultra-luminous X-ray sources (ULXs) may contain accreting intermediate-mass black holes (IMBHs). We review the current observational evidence for IMBH-ULXs. While low-luminosity ULXs with L_X <~ 10^39.5 erg/s (assuming isotropic emission) are consistent with mildly X-ray beamed high-mass X-ray binaries, there are a considerable number of ULXs with larger X-ray luminosities that are not easily explained by these models. Recent high-S/N XMM X-ray spectra are showing an increasing number of ULXs with ``cool disks'' -- accretion disks with multi-color blackbody inner disk temperatures kT_in ~ 0.1-0.2 keV, consistent with accreting IMBHs. Optical emission-line studies of ULX nebulae provide useful measurements of X-ray energetics, and can thus determine if the X-rays are emitted isotropically. Analysis of an optical spectrum of the Ho II ULX nebulae implies an X-ray energy source with ~10^40 erg/s is present, suggesting an isotropically-emitting IMBH. The spatial coincidence of ULXs with dense star clusters (young clusters and globular clusters) suggests that IMBHs formed in these clusters could be the compact objects in the associated ULXs. Quasi-periodic oscillations and frequency breaks in XMM power-density spectra of ULXs also suggest that the black hole masses are more consistent with IMBHs than stellar-mass black holes. Since _all_ of these ULXs with evidence for IMBHs are high-luminosity ULXs, i.e., L_X >~ 10^40 ergs, we suggest that this class of ULXs is generally powered by accreting IMBHs.

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.

We present results of numerical simulations of sequences of binary-single scattering events of black holes in dense stellar environments. The simulations cover a wide range of mass ratios from equal mass objects to 1000:10:10 solar masses and compare purely Newtonian simulations to simulations in which Newtonian encounters are interspersed with gravitational wave emission from the binary. In both cases, the sequence is terminated when the binary's merger time due to gravitational radiation is less than the arrival time of the next interloper. We find that black hole binaries typically merge with a very high eccentricity (0.93 < e < 0.95 pure Newtonian; 0.85 < e < 0.90 with gravitational wave emission) and that adding gravitational wave emission decreases the time to harden a binary until merger by ~ 30% to 40%. We discuss the implications of this work for the formation of intermediate-mass black holes and gravitational wave detection.

The high time resolution and large area of the Rossi X-ray Timing Explorer have been essential in the detection and characterization of high-frequency quasi-periodic variability in the flux from neutron stars in low-mass X-ray binaries. An unknown phenomenon prior to RXTE, kilohertz quasi-periodic oscillations (QPOs) have now been detected from more than twenty systems. Their high frequencies (up to 1330 Hz) imply that they are generated close to the neutron star, where general relativistic effects are expected to play an important role. I summarize current models for the kilohertz QPO phenomenon. In particular, I show that there is a significant domain of agreement among the models that can be used to constrain neutron star structure and look for signatures of highly curved spacetime in the properties of the QPOs.

This paper reports on our effort in modeling realistic astrophysical neutron star binaries in general relativity. We analyze under what conditions the conformally flat quasiequilibrium (CFQE) approach can generate ``astrophysically relevant'' initial data, by developing an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. We show that the CFQE assumptions significantly violate the Einstein field equations for corotating neutron stars at orbital separations nearly double that of the innermost stable circular orbit (ISCO) separation, thus calling into question the astrophysical relevance of the ISCO determined in the CFQE approach. With the need to start numerical simulations at large orbital separation in mind, we push for stable and long term integrations of the full Einstein equations for the binary neutron star system. We demonstrate the stability of our numerical treatment and analyze the stringent requirements on resolution and size of the computational domain for an accurate simulation of the system.

KHz QPOs have now been detected in more than twenty accreting neutron stars in low-mass binary systems. Two kHz QPOs are usually detected in each star. Burst oscillations and two kHz QPOs have recently been detected in the 401 Hz accretion-powered X-ray pulsar SAX J1808.4-3658. In this star the frequency of the burst oscillation is approximately equal to the star's spin frequency nu_spin whereas the frequency separation of the two kHz QPOs is approximately nu_spin/2. If as expected the frequency of the burst oscillations in other stars is also approximately nu_spin, the frequency separation is approximately nu_spin in some stars but approximately nu_spin/2 in others. A frequency separation approximately equal to nu_spin/2 is unexplained in all existing models of the kHz QPOs. Here we propose a modified version of the sonic-point beat-frequency model that can explain within a single framework why the frequency separation is close to nu_spin in some stars but close to nu_spin/2 in others. As in the original sonic-point model, the frequency nu_QPO2 of the upper kHz QPO is close to the orbital frequency nu_orb at the radius r_{sp} of the sonic point in the disk flow. We show that magnetic and radiation fields rotating with the star will preferentially excite vertical motions in the disk at the "spin-resonance'' radius r_{sr} where nu_orb - nu_spin is equal to the vertical epicyclic frequency, producing vertical motions in the disk that modulate the X-ray flux at approximately nu_QPO2 - nu_spin or approximately nu_QPO2 - nu_spin/2, depending on whether the disk flow at r_{sr} is smooth or clumped. This sonic-point and spin-resonance model can also explain quantitatively the decrease of the kHz QPO frequency separation with increasing accretion rate that is observed in many sources.

Evidence has been mounting for the existence of black holes with masses from 10^2 to 10^4 M_Solar associated with stellar clusters. Such intermediate-mass black holes (IMBHs) will encounter other black holes in the dense cores of these clusters. The binaries produced in these interactions will be perturbed by other objects as well thus changing the orbital characteristics of the binaries. These binaries and their subsequent mergers due to gravitational radiation are important sources of gravitational waves. We present the results of numerical simulations of high mass ratio encounters, which help clarify the interactions of intermediate-mass black holes in globular clusters and help determine what types of detectable gravitational wave signatures are likely.

Black holes with hundreds to thousands of solar masses are more massive than can be formed from a single star in the current universe, yet the best candidates for these objects are not located in gas-rich environments where gradual accretion could build up the mass. Three main formation scenarios have been suggested in the literature: that intermediate-mass black holes are the remnants of the first, metal-poor, stars; that they result from direct collisions in young stellar clusters; or that they are produced by gradual interactions and mergers of compact objects in old dense clusters. We discuss each of these in turn and speculate on future observations that may help sharpen our understanding of the formation of intermediate-mass black holes.

One often-used approximation in the study of binary compact objects (i.e., black holes and neutron stars) in general relativity is the instantaneously circular orbit assumption. This approximation has been used extensively, from the calculation of innermost circular orbits to the construction of initial data for numerical relativity calculations. While this assumption is inconsistent with generic general relativistic astrophysical inspiral phenomena where the dissipative effects of gravitational radiation cause the separation of the compact objects to decrease in time, it is usually argued that the timescale of this dissipation is much longer than the orbital timescale so that the approximation of circular orbits is valid. Here, we quantitatively analyze this approximation using a post-Newtonian approach that includes terms up to order ({Gm/(rc^2)})^{9/2} for non-spinning particles. By calculating the evolution of equal mass black hole / black hole binary systems starting with circular orbit configurations and comparing them to the more astrophysically relevant quasicircular solutions, we show that a minimum initial separation corresponding to at least 6 (3.5) orbits before plunge is required in order to bound the detection event loss rate in gravitational wave detectors to < 5% (20%). In addition, we show that the detection event loss rate is > 95% for a range of initial separations that include all modern calculations of the innermost circular orbit (ICO).

An approach to general relativity based on conformal flatness and quasiequilibrium (CFQE) assumptions has played an important role in the study of the inspiral dynamics and in providing initial data for fully general relativistic numerical simulations of coalescing compact binaries. However, the regime of validity of the approach has never been established. To this end, we develop an analysis that determines the violation of the CFQE approximation in the evolution of the binary described by the full Einstein theory. With this analysis, we show that the CFQE assumption is significantly violated even at relatively large orbital separations in the case of corotational neutron star binaries. We also demonstrate that the innermost stable circular orbit (ISCO) determined in the CFQE approach for corotating neutron star binaries may have no astrophysical significance.

We have analyzed an XMM-Newton observation of the nearby spiral galaxy NGC 1313, which contains two "ultra-luminous" X-ray (ULX) sources. We measure isotropic luminosities of L_X = 2.0 * 10^(40) erg/s and L_X = 6.6 * 10^(39) erg/s for NGC 1313 X-1 and X-2 (0.2-10.0 keV, assuming a distance of 3.7 Mpc). The spectra statistically require soft and hard spectral components to describe the continuum emission; some prior studies of ULXs have claimed cool soft components with lower statistics. The improvement over several single-component models exceeds the 8 sigma level of confidence for X-1; the improvement for X-2 is significant at the 3 sigma level. The soft components in these ULX spectra are well-fit by multi-color disk blackbody models with color temperatures of kT = 150 eV. This temperature differs markedly from those commonly measured in the spectra of stellar-mass (10 M_sun) black holes in their brightest states (kT ~ 1 keV). It is expected that the temperature of an accretion disk orbiting a black hole should decrease with increasing black hole mass. If the soft components we measure are due to emission from the inner region of an accretion disk, and disks extend close to the innermost stable circular orbit at the accretion rates being probed, the low color temperatures may be interpreted as spectroscopic evidence of black holes with intermediate masses: M_BH ~ 10^(3) M_sun. Simple Eddington scaling arguments suggest a mimum mass of M_BH ~ 10^(2) M_sun. NGC 1313 X-1 and X-2 are found in optical nebulae (Pakull and Mirioni 2002), which may indicate that anisotropic emission geometries are unlikely to account for the fluxes observed.

Recent X-ray observations of galaxies with ROSAT, ASCA, and Chandra have revealed numerous bright off-center point sources which, if isotropic emitters, are likely to be intermediate-mass black holes, with hundreds to thousands of solar masses. The origin of these objects is under debate, but observations suggest that a significant number of them currently reside in young high-density stellar clusters. There is also growing evidence that some Galactic globular clusters harbor black holes of similar mass, from observations of stellar kinematics. In such high-density stellar environments, the interactions of intermediate-mass black holes are promising sources of gravitational waves for ground-based and space-based detectors. Here we explore the signal strengths of binaries containing intermediate-mass black holes or stellar-mass black holes in dense stellar clusters. We estimate that a few to tens per year of these objects will be detectable during the last phase of their inspiral with the advanced LIGO detector, and up to tens per year will be seen during merger, depending on the spins of the black holes. We also find that if these objects reside in globular clusters then tens of sources will be detectable with LISA from the Galactic globular system in a five year integration, and similar numbers will be detectable from more distant galaxies. The signal strength depends on the eccentricity distribution, but we show that there is promise for strong detection of pericenter precession and Lense-Thirring precession of the orbital plane. We conclude by discussing what could be learned about binaries, dense stellar systems, and strong gravity if such signals are detected.

We describe a new numerical 3D relativistic hydrodynamical code, the results of validation tests, and a comparison with earlier, 2D studies. The 3D code has been used to study the deflection and precession of relativistic flows. We find that even quite fast jets (gamma~10) can be significantly influenced by impinging on an oblique density gradient, exhibiting a rotation of the Mach disk in the jet's head. The flow is bent via a potentially strong, oblique internal shock that arises due to asymmetric perturbation of the flow by its cocoon. In extreme cases this cocoon can form a marginally relativistic flow orthogonal to the jet, leading to large scale dynamics quite unlike that normally associated with astrophysical jets. Exploration of a gamma=5 flow subject to a large amplitude precession (semi-angle 11.25dg) shows that it retains its integrity, with modest reduction in Lorentz factor and momentum flux, for almost 50 jet-radii, but thereafter, the collimated flow is disrupted. The flow is approximately ballistic, with velocity vectors not aligned with the local jet `wall'. We consider simple estimators of the flow emissivity in each case and conclude that a) while the oblique internal shocks which mediate a small change in the direction of the deflected flows have little impact on the global dynamics, significantly enhanced flow emission (by a factor of 2-3) may be associated with such regions; and b) the convolution of rest frame emissivity and Doppler boost in the case of the precessed jet invariably leads to a core-jet-like structure, but that intensity fluctuations in the jet cannot be uniquely associated with either change in internal conditions or Doppler boost alone, but in general are a combination of both factors.

The spin periods of the ten observed Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters lie in the very narrow range 6-12 s. We use a point likelihood technique to assess the constraints this clustering imposes on the birth period and on the final period of such systems. We consider a general law for their spin evolution described by a constant braking index. We find that, for positive values of the braking index, the observed clustering requires an upper cut-off period that is very close to the maximum observed period of about 12 s. We also show that the constraint on the birth period depends very strongly on the assumed value of the braking index n, ranging from a few milliseconds for n>2 to a few seconds for n<2. We discuss the possible ways of tightening these constraints based on similarities with the population of radio pulsars and with future observations of such sources with current X-ray observatories.

In the high density cores of globular clusters, multibody interactions are expected to be common, with the result that black holes in binaries are hardened by interactions. It was shown by Sigurdsson & Hernquist (1993) and others that 10 solar mass black holes interacting exclusively by three-body encounters do not merge in the clusters themselves, because recoil kicks the binaries out of the clusters before the binaries are tight enough to merge. Here we consider a new mechanism, involving four-body encounters. Numerical simulations by a number of authors suggest that roughly 20-50% of binary-binary encounters will eject one star but leave behind a stable hierarchical triple. If the orbital plane of the inner binary is strongly tilted with respect to the orbital plane of the outer object, a secular Kozai resonance, first investigated in the context of asteroids in the Solar System, can increase the eccentricity of the inner body significantly. We show that in a substantial fraction of cases the eccentricity is driven to a high enough value that the inner binary will merge by gravitational radiation, without a strong accompanying kick. Thus the merged object remains in the cluster; depending on the binary fraction of black holes and the inclination distribution of newly-formed hierarchical triples, this mechanism may allow massive black holes to accumulate through successive mergers in the cores of globular clusters. It may also increase the likelihood that stellar-mass black holes in globular clusters will be detectable by their gravitational radiation.

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.

As increasingly precise information about the spectrum of the cosmic microwave background fluctuations is gathered with balloon and satellite experiments, interest has grown in foreground sources of opacity affecting these observations. One potentially important source is electron scattering produced by a post-recombination luminosity source, which would significantly attenuate the higher harmonics in the spectrum. If such an ionization source exists, then it would also heat the universe, hence increasing the Jeans mass and suppressing early gravitational structure formation. Here we consider the effects of such heating. We concentrate on one type of ionization source: luminosity generated by accretion onto primordial compact objects. We show that if such objects generate enough luminosity to affect the CMB power spectrum, then they would produce enough heat to prevent the formation of 1 sigma collapsed objects until redshifts of about 5, significantly less than the redshift at which baryonic collapse could otherwise occur. Such processes would leave signatures detectable by upcoming instruments such as NGST, SIRTF, and SWIFT.

The discovery of numerous non-nuclear X-ray point sources with luminosities L>10^39 erg/s in several starburst galaxies has stimulated speculation about their nature and origin. The strong variability seen in several sources points to massive black holes as the central engines. If the flux is isotropic, the luminosities range up to roughly 10^41 erg/s, implying masses of M>10^3 Msun if the luminosity is sub-Eddington. Here we explore a model for these sources. We suggest that in some tens of percent of globular clusters a very massive black hole, M>50 Msun, is formed. This black hole sinks in <10^6 yr to the center of the cluster, where in the 10^10 yr lifetime of the cluster it accretes roughly 10^3 Msun, primarily in the form of lighter black holes. Unlike less massive black holes in binaries, which are flung from clusters by recoil before they can merge gravitationally, a >50 Msun black hole has enough inertia that it remains bound to the cluster. We suggest that 10^3 Msun black holes may be common in the centers of dense globular clusters, and may therefore exist in some tens of percent of current globulars. If the cluster later merges with its host galaxy, accretion from young star clusters in molecular clouds by the black hole can generate luminosity consistent with that observed. We also consider the detectability of massive black holes in globular clusters with gravitational wave detectors such as LISA and LIGO, and speculate on future observations that may test our predictions.

We present an algorithm for solving the general relativistic initial value equations for a corotating polytropic star in quasicircular orbit with a nonspinning black hole. The algorithm is used to obtain initial data for cases where the black hole mass is 1, 3, and 10 times larger than the mass of the star. By analyzing sequences of constant baryon mass, constant black hole mass initial data sets and carefully monitoring the numerical error, we find innermost stable circular orbit (ISCO) configuration for these cases. While these quasiequilibrium, conformally flat sequences of initial data sets are not true solutions of the Einstein equations (each set, however, solves the full initial value problem), and thus, we do not expect the ISCO configurations found here to be completely consistent with the Einstein equations, they will be used as convenient starting points for future numerical evolutions of the full 3+1 Einstein equations.

We perform fully relativistic hydrodynamic simulations of the deceleration and lateral expansion of a relativistic jet as it expands into an ambient medium. The hydrodynamic calculations use a 2D adaptive mesh refinement (AMR) code, which provides adequate resolution of the thin shell of matter behind the shock. We find that the sideways propagation is different than predicted by simple analytic models. The physical conditions at the sides of the jet are found to be significantly different than at the front of the jet, and most of the emission occurs within the initial opening angle of the jet. The light curves, as seen by observers at different viewing angles with respect to the jet axis, are then calculated assuming synchrotron emission. For an observer along the jet axis, we find a sharp achromatic `jet break' in the light curve at frequencies above the typical synchrotron frequency, at $t_{jet}\approx 5.8(E_{52}/n_1)^{1/3}(\theta_0/0.2)^{8/3}$ days, while the temporal decay index $\alpha$ ($F_{\nu}\propto t^{\alpha}$) after the break is steeper than $-p$ ($\alpha=-2.85$ for $p=2.5$). At larger viewing angles $t_{jet}$ increases and the jet break becomes smoother.

X-ray timing observations of neutron stars and black holes are among the few available probes of ultrastrong magnetic fields, strong gravity, high densities, and the propagation of thermonuclear burning. Here we review the evidence for these effects revealed with data from the Rossi Explorer in the last five years. We also discuss the exciting prospects for making the first quantitative tests of strong-gravity general relativistic predictions with a large-area X-ray timing mission.

The first extrasolar planets were discovered in 1992 around the millisecond pulsar PSR 1257+12. We show that recent developments in the study of accretion onto magnetized stars, plus the existence of the innermost, moon-sized planet in the PSR 1257+12 system, suggest that the pulsar was born with approximately its current rotation frequency and magnetic moment. If so, this has important implications for the formation and evolution of neutron star magnetic fields as well as for the formation of planets around pulsars. In particular, it suggests that some and perhaps all isolated millisecond pulsars may have been born with high spin rates and low magnetic fields instead of having been recycled by accretion.

Previous work on the sonic-point beat-frequency (SPBF) model of the kilohertz quasi-periodic oscillations (QPOs) observed in the X-ray flux from neutron stars in low-mass binary systems has shown that it naturally explains many properties of these QPOs. These include the existence of just two principal QPOs in a given source, the commensurability of the frequency separation \Dnu\ of the two kilohertz QPOs and the spin frequency \nus inferred from burst oscillations, and the high frequencies, coherence, and amplitudes of these QPOs. Here we show that the SPBF model predicts that \Dnu is less than but close to \nus, consistent with the observed differences between \Dnu and \nus. It also explains naturally the decrease in \Dnu with increasing QPO frequency seen in some sources and the plateau in the QPO frequency--X-ray flux observed in 4U 1820-30. The model fits well the QPO frequency behavior observed in Sco X-1, 4U 1608-52, 4U 1728-34, and 4U 1820-30 (chisqdof = 0.4-2.1), giving masses ranging from 1.59 to 2.0 Msun and spin rates ranging from 279 to 364 Hz. In the SPBF model, the kilohertz QPOs are effects of strong-field gravity. Thus, if the model is validated, the kilohertz QPOs can be used not only to determine the properties of neutron stars but also to explore quantitatively general relativistic effects in the strong-field regime.

Kilohertz QPO sources are reasonably well-characterized observationally, but many questions remain about the theoretical framework for these sources and the consequent implications of the observations for disk physics, strong gravity, and dense matter. We contrast the predictions and implications of the most extensively studied class of kilohertz QPO models, the beat-frequency models, with those of alternative classes of models. We also discuss the expected impact of new observations of these sources with satellites such as Chandra, XMM, and Constellation-X.

Observations with RXTE have revealed kilohertz quasi-periodic brightness oscillations (QPOs) from nearly twenty different neutron-star low-mass X-ray binaries (LMXBs). These frequencies often appear as a pair of kilohertz QPOs in a given power density spectrum. In many models the higher-frequency of these QPOs is a beaming oscillation at the frequency of a nearly circular orbit at some radius near the neutron star. In such models it is expected that there will also be beaming oscillations at the stellar spin frequency and at overtones of the orbital frequency, but no strong QPOs have been detected at these frequencies. We therefore examine the processes that can attenuate beaming oscillations near neutron stars, and in doing so extend the work on this subject that was initiated by the discovery of lower-frequency QPOs from LMXBs. Among our main results are (1)in a spherical scattering cloud, all overtones of rotationally modulated beaming oscillations are attenuated strongly, not just the even harmonics, and (2)it is possible to have a relatively high-amplitude modulation near the star at, e.g., the stellar spin frequency, even if no peak at that frequency is detectable in a power density spectrum taken at infinity. We discuss the application of these results to modeling of kilohertz QPOs.

Many lines of evidence suggest that nonbaryonic dark matter constitutes roughly 30% of the critical closure density, but the composition of this dark matter is unknown. One class of candidates for the dark matter is compact objects formed in the early universe, with typical masses M between 0.1 and 1 solar masses to correspond to the mass scale of objects found with microlensing observing projects. Specific candidates of this type include black holes formed at the epoch of the QCD phase transition, quark stars, and boson stars. Here we show that accretion onto these objects produces substantial ionization in the early universe, with an optical depth to Thomson scattering out to z=1100 of approximately tau=2-4 [f_CO\epsilon_{-1}(M/Msun)]^{1/2} (H_0/65)^{-1}, where \epsilon_{-1} is the accretion efficiency \epsilon\equiv L/{\dot M}c^2 divided by 0.1 and f_CO is the fraction of matter in the compact objects. The current upper limit to the scattering optical depth, based on the anisotropy of the microwave background, is approximately 0.4. Therefore, if accretion onto these objects is relatively efficient, they cannot be the main component of nonbaryonic dark matter.

The discovery of high-amplitude brightness oscillations during type I X-ray bursts from six low-mass X-ray binaries has provided a powerful new tool to study the properties of matter at supranuclear densities, the effects of strong gravity, and the propagation of thermonuclear burning. There is substantial evidence that these brightness oscillations are produced by spin modulation of one or two localized hot spots confined to the stellar surface. It is therefore important to calculate the expected light curves produced by such hot spots under various physical assumptions, so that comparison with the observed light curves may most sensitively yield information about the underlying physical quantities. In this paper we make general relativistic calculations of the light curves and oscillation amplitudes produced by a rotating neutron star with one or two hot spots as a function of spot size, stellar compactness, rotational velocity at the stellar surface, spot location, orientation of the line of sight of the observer, and the angular dependence of the surface specific intensity. We find that stellar rotation and beaming of the emission tend to increase the observed oscillation amplitudes whereas greater compactness and larger spot size tend to decrease them. By applying these results to 4U 1636--536, we show that this source must have two emitting spots and place strong constraints on the neutron star's magnetic field geometry. We also show that the data on the phase lags between photons of different energies in the persistent pulsations in SAX J1808--58 can be fit well with a model in which the observed hard leads are due to Doppler beaming.

We analyze 426 observations of the bursting pulsar GRO J1744-28 by Ulysses and BATSE. Triangulating each burst, and statistically combining the triangulation annuli, we obtain a 3 sigma error ellipse whose area is 532 sq. arcsec. The accuracy of this statistical method has been independently verified with observations of the soft gamma repeater SGR1900+14. The ellipse is fully contained within the 1 ' radius ASCA error circle of the soft X-ray counterpart, and partially overlaps the 10 " radius ROSAT error circle of a source which may also be the soft X-ray counterpart. A variable source which has been proposed as a possible IR counterpart lies at the edge of the 3 sigma error ellipse, making it unlikely from a purely statistical point of view to be associated with the bursting pulsar.

We describe fully relativistic three dimensional calculations of the slowing down and spreading of a relativistic jet by an external medium like the ISM. We calculate the synchrotron spectra and light curves using the conditions determined by the hydrodynamic calculations. Preliminary results with a moderate resolution are presented here. Higher resolution calculations are in progress.

Many thermonuclear X-ray bursts exhibit brightness oscillations. The brightness oscillations are thought to be due to the combined effects of non-uniform nuclear burning and rotation of the neutron star. The waveforms of the oscillations contain information about the size and number of burning regions. They also contain substantial information about the mass and radius of the star, and hence about strong gravity and the equation of state of matter at supranuclear densities. We have written general relativistic ray-tracing codes that compute the waveforms and spectra of rotating hot spots as a function of photon energy. Using these codes, we survey the effect on the oscillation waveform and amplitude of parameters such as the compactness of the star, the spot size, the surface rotation velocity, and whether there are one or two spots. We also fit phase lag versus photon energy curves to data from the millisecond X-ray pulsar, SAX J1808--3658.

Shapiro put forth a conjecture stating that neutron stars in head-on collisions (infalling from infinity) will not collapse to black holes before neutrino cooling, independent of the mass of the neutron stars. In a previous paper we carried out a numerical simulation showing a counter example based on 1.4 $M_{\odot}$ neutron stars, and provided an analysis explaining why Shapiro's argument was not applicable for this case. A recent paper by Shapiro put forth an argument suggesting that numerical simulations of the 1.4 $M_{\odot}$ collisions could not disprove the conjecture with the accuracy that is presently attainable. We show in this paper that this argument is not applicable for the same reason that the Shapiro conjecture is not applicable to the 1.4 $M_{\odot}$ neutron star collision, namely, the collision is too dynamical to be treated by quasi-equilibrium arguments.

We study the question of prompt vs. delayed collapse in the head-on collision of two neutron stars. We show that the prompt formation of a black hole is possible, contrary to a conjecture of Shapiro which claims that collapse is delayed until after neutrino cooling. We discuss the insight provided by Shapiro's conjecture and its limitation. An understanding of the limitation of the conjecture is provided in terms of the many time scales involved in the problem. General relativistic simulations in the Einstein theory with the full set of Einstein equations coupled to the general relativistic hydrodynamic equations are carried out in our study.

The discovery of nearly coherent brightness oscillations during thermonuclear X-ray bursts from six neutron-star low-mass X-ray binaries has opened up a new way to study the propagation of thermonuclear burning, and may ultimately lead to greater understanding of thermonuclear propagation in other astrophysical contexts, such as in Type Ia supernovae. Here we report detailed analyses of the ~580 Hz brightness oscillations during bursts from 4U 1636-536. We investigate the bursts as a whole and, in more detail, the initial portions of the bursts. We analyze the ~580 Hz oscillations in the initial 0.75 seconds of the five bursts that were used in a previous search for a brightness oscillation at the expected ~290 Hz spin frequency, and find that if the same frequency model describes all five bursts there is insufficient data to require more than a constant frequency or, possibly, a frequency plus a frequency derivative. Therefore, although it is appropriate to use an arbitrarily complicated model of the ~580 Hz oscillations to generate a candidate waveform for the ~290 Hz oscillations, models with more than two parameters are not required by the data. For the bursts as a whole we show that the characteristics of the brightness oscillations vary greatly from burst to burst. We find, however, that in at least one of the bursts, and possibly in three of the four that have strong brightness oscillations throughout the burst, the oscillation frequency reaches a maximum several seconds into the burst and then decreases. This behavior has not been reported previously for burst brightness oscillations, and it poses a challenge to the standard burning layer expansion explanation for the frequency changes.

We study the horizontal branch oscillations (HBO) and the two simultaneous kilohertz QPOs comparing their properties in five Z sources with the predictions of the magnetospheric beat-frequency and Lense-Thirring precession models. We find that the variation of the HBO frequency with accretion rate predicted by the magnetospheric beat-frequency model for a purely dipolar stellar magnetic field and a radiation-pressure-dominated inner accretion disk agrees well with the observed variation. This model implies that the neutron stars in the Z sources are near magnetic spin equilibrium, that their magnetic field strengths are ~10^9-10^10 G, and that the critical fastness parameter for these sources is >0.8. If the frequency of the upper kilohertz QPO is an orbital frequency in the accretion disk, the magnetospheric beat-frequency model requires that a small fraction of the gas in the disk does not couple strongly to the stellar magnetic field at 3-4 stellar radii but instead drifts slowly inward in nearly circular orbits until it is within a few kilometers of the neutron star surface. The Lense-Thirring precession model is consistent with the observed magnitudes of the HBO frequencies only if the moments of inertia of the neutron stars in the Z sources are 4-5 times larger than the largest values predicted by realistic neutron-star equations of state. The trend of the correlation between the HBO and upper kHz QPO frequencies is not consistent with this model. We argue that the change in the slope of the correlation between the frequency of the horizontal branch oscillation and the frequency of the upper kilohertz QPO, when the latter is greater than 850 Hz, is directly related to the varying frequency separation of the kilohertz QPOs.

We calculate the expected number of multiply-imaged galaxies in the Hubble Deep Field (HDF), using photometric redshift information for galaxies with $m_I < 27$ that were detected in all four HDF passbands. A comparison of these expectations with the observed number of strongly lensed galaxies places a lower limit on the current value of $\Omega_m-\Omega_{\Lambda}$, where $\Omega_m$ is the cosmological mass density of the universe and $\Omega_\Lambda$ is the normalized cosmological constant. Based on current estimates of the HDF luminosity function and associated uncertainties in individual parameters, our 95% confidence lower limit on $\Omega_m-\Omega_{\Lambda}$ is between -0.44, if there are no strongly lensed galaxies in the HDF, and -0.73, if there are two strongly lensed galaxies in the HDF. If the only lensed galaxy in the HDF is the one presently viable candidate, then, in a flat universe ($\Omega_m+\Omega_\Lambda=1$), $\Omega_{\Lambda} < 0.79$ (95% C.L.). These lower limits are compatible with estimates based on high-redshift supernovae and with previous limits based on gravitational lensing.

We calculate the expected number of multiply-imaged galaxies in the Hubble Deep Field (HDF), using photometric redshift information for galaxies with m_I < 27 that were detected in all four HDF passbands. A comparison of these expectations with the observed number of strongly lensed galaxies constrains the current value of Omega_m-Omega_Lambda, where Omega_m is the mean mass density of the universe and Omega_Lambda is the normalized cosmological constant. Based on current estimates of the HDF luminosity function and associated uncertainties in individual parameters, our 95% confidence lower limit on Omega_m-Omega_Lambda ranges between -0.44, if there are no strongly lensed galaxies in the HDF, and -0.73, if there are two strongly lensed galaxies in the HDF. If the only lensed galaxy in the HDF is the one presently viable candidate, then, in a flat universe (Omega_m+Omega_Lambda=1), Omega_Lambda < 0.79 (95% C.L.). These limits are compatible with estimates based on high-redshift supernovae and with previous limits based on gravitational lensing.

This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and convergence of our general relativistic hydrodynamic treatment and its coupling to the spacetime evolutions described by the full set of Einstein equations with a perfect fluid source. The numerical treatment of the general relativistic hydrodynamic equations is based on high resolution shock capturing schemes. These schemes rely on the characteristic information of the system. A spectral decomposition for general relativistic hydrodynamics suitable for a general spacetime metric is presented. Evolutions based on three different approximate Riemann solvers coupled to four different discretizations of the Einstein equations are studied and compared. The coupling between the hydrodynamics and the spacetime (the right and left hand side of the Einstein equations) is carried out in a treatment which is second order accurate in {\it both} space and time. Convergence tests for all twelve combinations with a variety of test beds are studied, showing consistency with the differential equations and correct convergence properties. The test-beds examined include shocktubes, Friedmann-Robertson-Walker cosmology tests, evolutions of self-gravitating compact (TOV) stars, and evolutions of relativistically boosted TOV stars. Special attention is paid to the numerical evolution of strongly gravitating objects, e.g., neutron stars, in the full theory of general relativity, including a simple, yet effective treatment for the surface region of the star (where the rest mass density is abruptly dropping to zero).

Strong, coherent, quasi-periodic brightness oscillations (QPOs) with frequencies ranging from about 300 Hz to 1200 Hz have been discovered with the Rossi X-ray Timing Explorer in the X-ray emission from some fifteen neutron stars in low-mass binary systems. Two simultaneous kilohertz QPOs differing in frequency by 250 to 350 Hertz have been detected in twelve of the fifteen sources. Here we propose a model for these QPOs. In this model the X-ray source is a neutron star with a surface magnetic field of 10^7 to 10^10 G and a spin frequency of a few hundred Hertz, accreting gas via a Keplerian disk. The frequency of the higher-frequency QPO in a kilohertz QPO pair is the Keplerian frequency at a radius near the sonic point at the inner edge of the Keplerian flow whereas the frequency of the lower-frequency QPO is approximately the difference between the Keplerian frequency at a radius near the sonic point and the stellar spin frequency. This model explains naturally many properties of the kilohertz QPOs, including their frequencies, amplitudes, and coherence. We show that if the frequency of the higher-frequency QPO in a pair is an orbital frequency, as in the sonic-point model, the frequencies of these QPOs place interesting upper bounds on the masses and radii of the neutron stars in the kilohertz QPO sources and provide new constraints on the equation of state of matter at high densities. Further observations of these QPOs may provide compelling evidence for the existence of a marginally stable orbit, confirming a key prediction of general relativity in the strong-field regime.

The discovery of high-frequency brightness oscillations in thermonuclear X-ray bursts from several neutron-star low-mass X-ray binaries has important implications for the beat frequency model of kilohertz quasi-periodic brightness oscillations, the propagation of nuclear burning, the structure of the subsurface magnetic fields in neutron stars, and the equation of state of high-density matter. These implications depend crucially on whether the observed frequency is the stellar spin frequency or its first overtone. Here we report an analysis of five bursts from 4U 1636-536 which exhibit strong oscillations at approximately 580 Hz. We show that combining the data from the first 0.75 seconds of each of the five bursts yields a signal at 290 Hz that is significant at the $4\times 10^{-5}$ level when the number of trials is taken into account. This strongly indicates that 290 Hz is the spin frequency of this neutron star and that 580 Hz is its first overtone, in agreement with other arguments about this source but in contrast to suggestions in the literature that 580 Hz is the true spin frequency. The method used here, which is an algorithm for combining time series data from the five bursts so that the phases of the 580 Hz oscillations are aligned, may be used in any source to search for weak oscillations that have frequencies related in a definite way to the frequency of a strong oscillation.

The lack of bright host galaxies in several recently examined gamma-ray burst (GRB) error boxes suggests that the redshifts of cosmological GRBs may be significantly higher than previously believed. On the other hand, the non-detection of multiple images in the BATSE 4B catalog implies an upper limit on the average redshift <z> of GRBs. Here, we calculate an upper limit to <z>, independent of the physical model for GRBs, using a new statistical lensing method that removes distance ambiguities, and thus permits accurate computation of the lensing rate at high z. The upper limit on <z> depends directly on the cosmological parameters Omega and Lambda. If there are no multiple images among the brightest 80% of the first 1802 bursts in the BATSE 4B catalog, then, at the 95% confidence level, <z> is less than 2.2, 2.8, 4.3, or 5.3 for (Omega,Lambda) values of (0.3, 0.7), (0.5, 0.5), (0.5, 0.0), or (1.0, 0.0), respectively. The 68% upper limit to the average redshift is comparable to or less than the median redshift of GRBs in scenarios in which the GRB rate is proportional to the rate of star formation, for any cosmology. The uncertainty in the lensing rate --- arising from uncertainties in the cosmological parameters and in the number density and average velocity dispersion of galaxies --- will be reduced significantly in the next few years by a new generation of experiments and surveys. Moreover, the continued increase in the number of GRBs observed by BATSE will greatly constrain their redshift distribution.

Quasi-periodic X-ray brightness oscillations (QPOs) with frequencies around a kilohertz have now been discovered in more than a dozen neutron stars in low-mass X-ray binary systems using the Rossi X-ray Timing Explorer. There is strong evidence that the frequencies of the kilohertz oscillations are the orbital frequencies of accreting gas in nearly circular orbits around these stars. Some stars that produce kilohertz QPOs may have spin frequencies greater than 400 Hertz. For spin rates this high, first-order analytic treatments of the effects of the star's rotation on its structure and the spacetime are inaccurate. Here we use the results of a large number of fully relativistic, self-consistent numerical calculations of the stellar structure of rapidly rotating neutron stars and the interior and exterior spacetime to investigate the constraints on the properties of such stars that can be derived if stable circular orbits of various frequencies are observed. We have computed the equatorial radius of the star, the radius of the innermost stable circular orbit, and the frequency of the highest-frequency stable circular orbit as functions of the stellar spin rate, for spin rates up to the maximum possible and for several illustrative equations of state. Our calculations show that the upper bounds on the stiffness of neutron star matter implied by a given orbital frequency are typically significantly stricter for stars with spin frequencies more than 400 Hertz than for slowly rotating stars.

One of the most dramatic discoveries made with the Rossi X-Ray Timing Explorer is that many accreting neutron stars in low-mass binary systems produce strong, remarkably coherent, high-frequency X-ray brightness oscillations. The 325-1200 Hz quasi-periodic oscillations (QPOs) observed in the accretion-powered emission are thought to be produced by gas orbiting very close to the neutron star, whereas the 360-600 Hz brightness oscillations seen during thermonuclear X-ray bursts are produced by one or two hot spots rotating with the star and have frequencies equal to the stellar spin frequency or its first overtone. The oscillations constrain the masses and radii of these neutron stars, which are thought to be the progenitors of the millisecond pulsars. Modeling indicates that the stars have spin frequencies 250-350 Hz and magnetic fields 1e7-5e9 G.

One of the most dramatic discoveries made so far with the Rossi X-Ray Timing Explorer is that many accreting neutron stars with weak magnetic fields generate strong, remarkably coherent, high-frequency X-ray brightness oscillations. The 325-1200 Hz quasi-periodic oscillations (QPOs) observed in the accretion-powered emission are almost certainly produced by gas orbiting very close to the stellar surface and have frequencies related to the orbital frequencies of the gas. The 360-600 Hz brightness oscillations seen during thermonuclear X-ray bursts are produced by one or two hotter regions on the stellar surface and have frequencies equal to the stellar spin frequency or its first overtone. Measurements of these oscillations are providing tight upper bounds on the masses and radii of neutron stars, and important new constraints on the equation of state of neutron star matter.

The dramatic discovery with the Rossi X-Ray Timing Explorer satellite of remarkably coherent $\sim$300--1200 Hz oscillations in the X-ray brightness of some sixteen neutron stars in low-mass binary systems has spurred theoretical modeling of these oscillations and investigation of their implications for the neutron stars and accretion flows in these systems. High-frequency oscillations are observed both during thermonuclear X-ray bursts and during intervals of accretion-powered emission and appear to be a characteristic feature of disk-accreting neutron stars with weak magnetic fields. In this review we focus on the high-frequency quasi-periodic oscillations (QPOs) seen in the accretion-powered emission. We first summarize the key properties of these kilohertz QPOs and then describe briefly the models that have been proposed to explain them. The existing evidence strongly favors beat-frequency models. We mention several of the difficulties encountered in applying the magnetospheric beat-frequency model to the kilohertz QPOs. The most fully developed and successful model is the sonic-point beat-frequency model. We describe the work on this model in some detail. We then discuss observations that could help to distinguish between models. We conclude by noting some of the ways in which study of the kilohertz QPOs may advance our understanding of dense matter and strong gravitational fields.

Gravitomagnetic precession near neutron stars and black holes has received much recent attention, particularly as a possible explanation of 15--60 Hz quasi-periodic brightness oscillations (QPOs) from accreting neutron stars in low-mass X-ray binaries, and of somewhat higher-frequency QPOs from accreting stellar-mass black holes. Previous analyses of this phenomenon have either ignored radiation forces or assumed for simplicity that the radiation field is isotropic, and in particular that there is no variation of the radiation field with angular distance from the rotational equatorial plane of the compact object. However, in most realistic accretion geometries (e.g., those in which the accretion proceeds via a geometrically thin disk) the radiation field depends on latitude. Here we show that in this case radiation forces typically have an important, even dominant, effect on the precession frequency of test particles in orbits that are tilted with respect to the star's rotational equator. Indeed, we find that even for accretion luminosities only a few percent of the Eddington critical luminosity, the precession frequency near a neutron star can be changed by factors of up to $\sim 10$. Radiation forces must therefore be included in analyses of precession frequencies near compact objects, in such varied contexts as low-frequency QPOs, warp modes of disks, and trapped oscillation modes. We discuss specifically the impact of radiation forces on models of low-frequency QPOs involving gravitomagnetic precession, and show that such models are rendered much less plausible by the effects of radiation forces.

The frequencies of the highest-frequency kilohertz QPOs recently discovered in some sixteen neutron stars in low-mass X-ray binary systems are most likely the orbital frequencies of gas in Keplerian orbit around these neutron stars. If so, these QPOs provide tight upper bounds on the masses and radii of these neutron stars and important new constraints on the equation of state of neutron star matter. If the frequency of a kilohertz QPO can be established as the orbital frequency of gas at the innermost stable circular orbit, this would confirm one of the key predictions of general relativity in the strong-gravity regime. If the spin frequency of the neutron star can also be determined, the frequency of the QPO would fix the mass of the neutron star for each assumed equation of state. Here we show how to derive mass and radius bounds, using the kilohertz QPOs, for nonrotating and slowly rotating stars, and discuss how these bounds are affected by rapid stellar rotation and radial radiation forces. We also describe observational results that would be strong evidence for the presence of an innermost stable circular orbit. No such strong evidence is present in current data, but future prospects are excellent.

The frequencies of the highest-frequency kilohertz QPOs recently discovered in some sixteen neutron stars in low-mass X-ray binary systems are most likely orbital frequencies. If so, these QPOs provide tight upper bounds on the masses and radii of these neutron stars and interesting new constraints on the equation of state of neutron star matter. If the frequency of a kilohertz QPO can be established as the orbital frequency of the innermost stable circular orbit, this would confirm one of the key predictions of general relativity in the strong-field regime. If the spin frequency of the neutron star can also be determined, the frequency of the QPO would fix the mass of the neutron star for each assumed equation of state. Here we describe how bounds on the stellar mass and radius can be derived and how these bounds are affected by the stellar spin. We also discuss detection of the innermost stable circular orbit.

The remarkable discovery of kilohertz quasi-periodic brightness oscillations (QPOs) in the accretion-powered emission from some sixteen neutron-star low-mass X-ray binary systems has led to much speculation about and theoretical modeling of the origin of these oscillations. It has also led to intense study of the implications that these QPOs have for the properties of neutron stars and of the accretion flow onto them. In this review we describe the strengths and weaknesses of the models that have been proposed for the kilohertz QPOs observed in the accretion-powered emission. We conclude that beat-frequency models, and in particular the sonic-point model, are currently the most promising. If these models are correct, the kilohertz QPOs provide the strongest constraints to date on the masses and radii of neutron stars in low-mass X-ray binaries, and on the equation of state of the dense matter in all neutron stars.

The discovery of high-amplitude brightness oscillations at the spin frequency or its first overtone in six neutron stars in low-mass X-ray binaries during type~1 X-ray bursts provides a powerful new way to constrain the compactness of these stars, and hence to constrain the equation of state of the dense matter in all neutron stars. Here we present the results of general relativistic calculations of the maximum fractional rms amplitudes that can be observed during bursts. In particular, we determine the dependence of the amplitude on the compactness of the star, the angular dependence of the emission from the surface, the rotational velocity at the stellar surface, and whether there are one or two emitting poles. We show that if two poles are emitting, as is strongly indicated by independent evidence in 4U 1636-536 and KS 1731-26, the resulting limits on the compactness of the star can be extremely restrictive. We also discuss the expected amplitudes of X-ray color oscillations and the observational signatures necessary to derive convincing constraints on neutron star compactness from the amplitudes of burst oscillations.

The energetics, spectrum, and composition of cosmic rays with energies below about $10^{15}$ eV are fairly well explained by models involving supernova shocks. In contrast, no widely accepted theory exists for the origin of ultra-high energy cosmic rays (UHECRs), which have energies above $10^{15}$ eV. Instead of proposing a specific model, here we place strong constraints on any model of UHECRs involving isolated neutron stars (no companions). We consider the total power requirements and show that the only viable power source associated with isolated neutron stars is rotation. Mechanisms based on accretion from the interstellar medium fall short of the necessary power despite the most optimistic assumptions. Power considerations also demonstrate that not enough rotational energy is tapped by a "propeller"-like acceleration of interstellar matter. The most promising source of energy is rotational spindown via magnetic braking. We examine microphysical energy loss processes near magnetized neutron stars and conclude that the most likely site for yielding UHECRs from isolated neutron stars is near or beyond the light cylinder.

Quasi-periodic brightness oscillations (QPOs) with frequencies in the range 300-1200 Hz have been detected from at least nine neutron star low-mass X-ray binaries. Here we summarize the sonic-point model for these brightness oscillations, which we present in detail elsewhere. If the sonic-point interpretation of kilohertz QPOs is confirmed, measurements of kilohertz QPO frequencies in low-mass X-ray binaries will provide new bounds on the masses and radii of neutron stars in these systems and new constraints on the equation of state of matter at high densities.

We present results from a series of optical (g and r-band) and near-infrared (K'-band) observations of the region of the sky including the entire XTE and ROSAT error circles for the ``Bursting Pulsar'' GRO J1744-28. These data were taken with the Astrophysical Research Consortium's 3.5-m telescope at Apache Point Observatory and with the 2.2-m telescope at the European Southern Observatory. We see no new object, nor any significant brightening of any known object, in these error circles, with the exception of an object detected in our 8 February 1996 image. This object has already been proposed as a near-infrared counterpart to GRO J1744-28. While it is seen in only two of our ten 8 February frames, there is no evidence that this is an instrumental artifact, suggesting the possibility of near-infrared flares from GRO J1744-28, similar to those that have been reported from the Rapid Burster. The distance to the ``Bursting Pulsar'' must be more than 2 kpc, and we suggest that it is more than 7 kpc.

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.

Using a Hartree-Fock formalism, we estimate energy levels and photon cross sections for atomic iron in magnetic fields B ~ 10^13 G. Computing ionization equilibrium and normal mode opacities with these data, we construct LTE neutron star model atmospheres at 5.5 < Log(T_eff) < 6.5 and compute emergent spectra. We examine the dependence of the emergent spectra on T_eff and B. We also show the spectral variation with the angle between the magnetic field and the atmosphere normal and describe the significant limb darkening in the X-ray band. These results are compared with recent detailed computations of neutron star H model atmospheres in high fields and with low field Fe and H model atmospheres constructed from detailed opacities. The large spectral differences for different surface compositions may be discernible with present X-ray data; we also note improvements needed to allow comparison of Fe models with high quality spectra.

We argue that the bursting, transient X-ray source GRO J1744-28 is a binary consisting of a very low-mass ($M \approx 0.2 M_\sun$), highly evolved giant star that is transferring mass by Roche lobe overflow onto a high-mass ($M \approx 1.8 M_\sun$) neutron star. We explore a picture in which the bursts are due to thermonuclear flashes in matter that has accreted onto the neutron star. We attribute the unusually hard spectra of the bursts and the high burst rate compared to normal Type I X-ray burst sources -- indeed the existence of unstable nuclear burning itself -- to the unusually strong surface magnetic field, $ B_s \approx 1 \times 10^{13}$ G, of the neutron star.

The association of all three soft gamma-ray repeaters (SGRs) with supernova remnants has made possible estimates of the distance to and luminosity of the sources of SGRs, which have provided a starting point for detailed modeling. One of the most popular classes of models involves strongly magnetized neutron stars, with surface dipole fields $B\sim 10^{14}-10^{15}$ G. In these ``magnetar" models, many otherwise negligible processes can play an important role. Here we consider the spectral effects of strong-field modifications to Compton scattering, in particular those related to the contribution of vacuum polarization to the dielectric tensor. Vacuum polarization introduces a density-dependent photon frequency, called the second vacuum frequency, at which the normal modes of polarization become nonorthogonal and the mean free path of photons decreases sharply. Monte Carlo simulations of photon propagation through a magnetized plasma show that this effect leads, under a wide range of physical conditions, to a broad absorption-like feature in the energy range $\sim$5 keV---40 keV.