Source author record

Abraham Loeb

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

316works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

316 published item(s)

preprint2026arXiv

Increased and Varied Radiation during the Sun's Encounters with Cold Clouds in the last 10 million years

Recent research raises the possibility that 3 and 7 million years ago, the Sun encountered massive clouds that shrank the heliosphere--the solar cocoon protecting our solar system--exposing Earth to its interstellar environment, in agreement with geological evidence from 60Fe and 244Pu isotopes. Here we show that during such encounters Earth was exposed to increased radiation in the form of high-energy particles. During periods of Earth's immersion in the heliosphere, it received particle radiation that we name Heliospheric Energetic Particles (HEPs). The intensity of < 10 MeV protons was at least an order of magnitude more intense than today's most extreme solar energetic particle (SEP) events. SEPs today last minutes to hours, but HEP exposure then lasted for extensive periods of several months, making it a prolonged external driver. During Earth's excursion outside the heliosphere, it was exposed to a galactic cosmic ray radiation with the intensity of < 1 GeV protons at least an order of magnitude more intense than today. Therefore, the space surrounding Earth was permeated by a variable high-energy radiation. We discuss the implications for Earth's climate and biodiversity.

preprint2026arXiv

Rare Near-Opposition Alignment of 3I/ATLAS on 22 January, 2026

We point out that on 22 January 2026, the interstellar object 3I/ATLAS will align to within an exceptionally small angle, alpha= 0.69 degrees, with the Earth-Sun axis. This rare alignment provides unique circumstances for measuring the opposition surge and polarimetric properties of interstellar cometary dust. We characterize the alignment geometry, outline key scientific opportunities, and define the observational requirements for data collection. Observations before and after the alignment time offer an unprecedented opportunity which may not repeat for decades, for characterizing the albedo, structure, and composition of interstellar matter.

preprint2026arXiv

Rotation Period of 3I/ATLAS After Perihelion from Jet Position Angle Wobble and Photometric Variability

We determine the post-perihelion rotation period of 3I/ATLAS using two independent diagnostics: the temporal modulation of the position angle (PA) of a persistent jet-like feature, and a time-series photometric light curve in the Gr (R) band. For the jet morphology, we measure the PA at multiple epochs by applying the Larson-Sekanina Rotational Gradient filter to Hubble Space Telescope images between November 20, 2025 and December 27, 2025, and model the phase-folded PA curve with weighted least-squares Fourier series up to two harmonics while scanning trial periods P to identify minima in \c{hi}2(P). For the photometry, we adopt the best-fit period from an independent 30-minute binned analysis (from a 0.25 meter telescope MPC L92) based on a refined \c{hi}2(P) profile for a sinusoidal model with nightly offsets. We find that the jet-PA modulation gives a period Pjet = 7.20 +/-0.05 h (adopting a conservative uncertainty dominated by sparse sampling and systematics), while the photometry yields Pphot = 7.136+/-0.001 h (formal 1σ uncertainty). Although the periods differ slightly, the offset is plausibly attributable to non-Gaussian systematics and aliasing. The combined data supports a post-perihelion rotation period of 7.1 h associated with precession of the jet structure around the rotation axis by 20°. The rotation axis is aligned with the sunward direction to within 20°

preprint2025arXiv

On the coincidence between the close passage of HD7977 and the Pliocene-Pleistocene transition

The Oort Cloud's dynamical evolution is significantly influenced by both the galactic tide and stellar flybys. This study investigates the particular case of HD7977's close encounter 2.47 Myr ago, which likely repopulated the Inner Oort Cloud and potentially triggered a significant comet shower on Earth. Our results demonstrate that the shower's intensity strongly depends on HD7977's impact parameter ($b$), with possible flyby distances ranging from 2,300 AU to $\sim$ 13,000 AU. For the closest approach ($b \sim 2,300$ AU), the terrestrial impact probability of 1 km comets increases by an order of magnitude compared to the steady state, slightly exceeding the asteroid impact probability at this size scale. We propose an analytical method to compute the probability of comet showers impacting Earth, which saves considerable computation time compared to N-body simulations. We identify a threshold diameter $D_0 = 2.25$ km for which yields $P = 1$ in our model, with $D_0$ following a logarithmic dependence on $b$. These findings suggest that HD7977's flyby may have caused an enhanced comet flux during the Pliocene-Pleistocene transition, which could plausibly be related to the environmental changes at this era.

preprint2025arXiv

Supermassive Stars Match the Spectral Signatures of JWST's Little Red Dots

The James Webb Space Telescope (JWST) has unveiled a population of enigmatic, compact sources at high redshift known as ``Little Red Dots'' (LRDs), whose physical nature remains a subject of intense debate. Concurrently, the rapid assembly of the first supermassive black holes (SMBHs) requires the formation of heavy seeds, for which supermassive stars (SMSs) are leading theoretical progenitors. In this work, we perform the first quantitative test of the hypothesis that LRDs are the direct observational manifestation of these primordial SMSs. We present a novel, first-principles pipeline generating synthetic spectra for a non-rotating, metal-free SMS up to $10^6 \, M_\odot$. We establish that its luminosity ($L_λ\approx 1.7 \times 10^{44} \, \text{erg} \, \text{s}^{-1} \, μ\text{m}^{-1}$ at 4050\,Å) provides a decisive constraint, matching prominent LRDs. Our model self-consistently reproduces their defining spectral features: the V-shaped Balmer break morphology is shown to be an intrinsic photospheric effect, while the complex line phenomenology, strong H$β$ in emission with other Balmer lines in absorption arises from non-LTE effects in a single stellar atmosphere. With wind and macroturbulent broadening, we match LRD spectra at $z=7.76$ and $z=3.55$, including the H$β$ width of MoM-BH*-1 to within 4\%. We predict a luminosity-dependent observability window, $\sim10^{4}$ yr for the most luminous systems and $10^{5}$--$10^{6}$ yr if $L_λ(4050\,\textÅ)$ is lower by 1--2 dex. These results provide a self-consistent alternative to multi-component obscured AGN scenarios and suggest JWST may be witnessing luminous stages of SMBH progenitors before collapse.

preprint2022arXiv

A Hot Subdwarf Model for the 18.18 Minute Pulsar GLEAM-X

We suggest that the recently discovered, enigmatic pulsar with a period of 18.18 minutes, GLEAM-X J162759.5-523504.3, is most likely a hot subdwarf (proto white dwarf). A magnetic dipole model explains the observed period and period-derivative for a highly magnetized ($\sim 10^8$G), hot subdwarf of typical mass $\sim 0.5M_\odot$ and radius $\sim 0.3R_\odot$, and an age of $\sim 3\times 10^4$yr. The subdwarf spin is close to its breakup speed and its spindown luminosity is near its Eddington limit, likely as a result of accretion from a companion.

preprint2022arXiv

A Statistical Detection of Wide Binary Systems in the Ultra-Faint Dwarf Galaxy Reticulum II

Binary stars can inflate the observed velocity dispersion of stars in dark matter dominated systems such as ultra-faint dwarf galaxies (UFDs). However, the population of binaries in UFDs is poorly constrained by observations, with preferred binary fractions for individual galaxies ranging from a few percent to nearly unity. Searching for wide binaries through nearest neighbor (NN) statistics (or the two-point correlation function) has been suggested in the literature, and we apply this method for the first time to detect wide binaries in a UFD. By analyzing the positions of stars in Reticulum II (Ret II) from Hubble Space Telescope images, we search for angularly resolved wide binaries in Ret II. We find that the distribution of their NN distances shows an enhancement at projected separations of $<1$ arcsec relative to a model containing no binaries. We show that such an enhancement can be explained by a wide~binary fraction of $f_b\approx0.007^{+0.008}_{-0.003}$ at separations of more than 3000 AU. Under the assumption that the binary separation distribution is similar to that in the Milky Way, the total binary fraction in Ret II may be on the order of 50%. We also use the observed magnitude distribution of stars in Ret II to constrain the initial mass function over the mass range $0.34-0.78~M_{\odot}$, finding that a shallow power-law slope of $1.01 \le α\le 1.15$ matches the data.

preprint2022arXiv

An Ocean Expedition by the Galileo Project to Retrieve Fragments of the First Large Interstellar Meteor CNEOS 2014-01-08

The earliest confirmed interstellar object, `Oumuamua, was discovered in the Solar System by Pan-STARRS in 2017, allowing for a calibration of the abundance of interstellar objects of its size $\sim 100\;$ m. This was followed by the discovery of Borisov, which allowed for a similar calibration of its size $\sim 0.4 - 1 \mathrm{\; km}$. One would expect a much higher abundance of significantly smaller interstellar objects, with some of them colliding with Earth frequently enough to be noticeable. Based on the CNEOS catalog of bolide events, we identified in 2019 the meteor detected at 2014-01-08 17:05:34 UTC as originating from an unbound hyperbolic orbit with 99.999\% confidence. In 2022, the U.S. Department of Defense has since verified that "the velocity estimate reported to NASA is sufficiently accurate to indicate an interstellar trajectory," making the object the first detected interstellar object and the first detected interstellar meteor. Here, we discuss the dynamical and compositional properties of CNEOS 2014-01-08, and describe our plan for an expedition to retrieve meteoritic fragments from the ocean floor.

preprint2022arXiv

Are the Newly-Discovered $z \sim 13$ Drop-out Sources Starburst Galaxies or Quasars?

The detection of two $z\sim 13$ galaxy candidates has opened a new window on galaxy formation at an era only $330$ Myr after the Big Bang. Here, we investigate the physical nature of these sources: are we witnessing star forming galaxies or quasars at such early epochs? If powered by star formation, the observed ultraviolet (UV) luminosities and number densities can be jointly explained if: (i) these galaxies are extreme star-formers with star formation rates $5-24\times$ higher than those expected from extrapolations of average lower-redshift relations; (ii) the star formation efficiency increases with halo mass and is countered by increasing dust attenuation from $z \sim 10-5$; (iii) they form stars with an extremely top-heavy initial mass function. The quasar hypothesis is also plausible, with the UV luminosity produced by black holes of $\sim 10^8 \, \mathrm{M_{\odot}}$ accreting at or slightly above the Eddington rate ($f_{\rm Edd}\sim 1.0$). This black hole mass at $z\sim13$ would require very challenging, but not implausible, growth parameters. If spectroscopically confirmed, these two sources will represent a remarkable laboratory to study the Universe at previously inaccessible redshifts.

preprint2022arXiv

Eccentricity Evolution in Gaseous Dynamical Friction

We analyse how drag forces modify the orbits of objects moving through extended gaseous distributions. We consider how hydrodynamic (surface area) drag forces and dynamical friction (gravitational) drag forces drive the evolution of orbital eccentricity. While hydrodynamic drag forces cause eccentric orbits to become more circular, dynamical friction drag can cause orbits to become more eccentric. We develop a semi-analytic model that accurately predicts these changes by comparing the total work and torque applied to the orbit at periapse and apoapse. We use a toy model of a radial power-law density profile, $ρ\propto r^{-γ}$, to determine that there is a critical $γ= 3$ power index which separates the eccentricity evolution in dynamical friction: orbits become more eccentric for $γ< 3$ and circularize for $γ> 3$. We apply these findings to the infall of a Jupiter-like planet into the envelope of its host star. The hydrostatic envelopes of stars are defined by steep density gradients near the limb and shallower gradients in the interior. Under the influence of gaseous dynamical friction, an infalling object's orbit will first decrease in eccentricity, then increase. The critical separation that delineates these regimes is predicted by the local density slope and is linearly dependent on polytropic index. More broadly, our findings indicate that binary systems may routinely emerge from common envelope phases with non-zero eccentricities that were excited by the dynamical friction forces that drove their orbital tightening.

preprint2022arXiv

Galactic Kites

I show that interstellar films of material thinner than a micron, drift away from the Galactic plane as a result of stellar radiation pressure. Such films, whether produced naturally by dust coagulation in proto-planetary disks or artificially by technological civilizations, would accumulate over the age of the Milky-Way and hover above the Galactic disk at a scale-height set gravitationally by the dark matter halo. Limits on scattered starlight imply that this population carries a fraction below 2x10^{-3} of the interstellar medium mass.

preprint2022arXiv

Gravitational Redshift for Wide Binaries in Gaia eDR3

The Doppler effect is commonly used to infer the velocity difference between stars based on the relative shifts in the rest-frame wavelengths of their spectral features. In wide binaries, the difference in gravitational redshift from the surfaces of the constituent stars with distinct compactness dominates at separations above 0.01 pc. I suggest that this effect became apparent for wide pairs in the Gaia eDR3 catalogue but incorrectly interpreted as a possible modification of Newtonian gravity in the internal kinematics of very wide binaries.

preprint2022arXiv

Implication of spin constraints by the Event Horizon Telescope on stellar orbits in the Galactic Center

The center of the Milky Way hosts the closest supermassive black hole, SgrA$^*$. Decades of near-infrared observations of our Galactic Center have shown the presence of a small population of stars (the so called S-star cluster) orbiting SgrA$^*$, which were recently reported to be arranged in two orthogonal disks. In this case, the timescale for Lense-Thirring precession of S-stars should be longer than their age, implying a low spin for SgrA$^*$. In contrast, the recent results by the Event Horizon Telescope favor a highly-spinning SgrA$^*$, which seem to suggest that the S-stars could not be arranged in disks. Alternatively, the spin of SgrA$^*$ must be small, suggesting that the models for its observed image are incomplete.

preprint2022arXiv

Implications for the Hubble tension from the ages of the oldest astrophysical objects

We use the ages of old astrophysical objects (OAO) in the redshift range $0 \lesssim z \lesssim 8$ as stringent tests of the late-time cosmic expansion history. Since the age of the Universe at any redshift is inversely proportional to $H_0$, requiring that the Universe be older than the oldest objects it contains at any redshift, provides an upper limit on $H_0$. Using a combination of galaxies imaged from the CANDELS program and various high-$z$ quasars, we construct an age-redshift diagram of $\gtrsim 100$ OAO up to $z \sim 8$. Assuming the $Λ$CDM model at late times, we find the 95\% confidence level upper limit $H_0<73.2\,{\rm km}/{\rm s}/{\rm Mpc}$, in slight disagreement with a host of local $H_0$ measurements. Taken at face value, and assuming that the OAO ages are reliable, this suggests that ultimately a combination of pre- and post-recombination ($z \lesssim 10$) new physics might be required to reconcile cosmic ages with early-time and local $H_0$ measurements. In the context of the Hubble tension, our results motivate the study of either combined global pre- and post-recombination modifications to $Λ$CDM, or local new physics which only affects the local $H_0$ measurements.

preprint2022arXiv

Lensing in the Darkness: A Bayesian Analysis of 22 Chandra Sources at $z \gtrsim 6$ Shows No Evidence of Lensing

More than $200$ quasars have been detected so far at $z > 6$, with only one showing clear signs of strong gravitational lensing. Some studies call for a missing population of lensed high-$z$ quasars, but their existence is still in doubt. A large fraction of high-$z$ quasars being lensed would have a significant effect on the shape of the intrinsic quasar luminosity function (QLF). Here, we perform the first systematic search for lensed X-ray-detected quasars at $z \gtrsim 6$ employing a Bayesian analysis, with the code BAYMAX, to look for morphological evidence of multiple images that may escape a visual inspection. We analyzed a sample of 22 quasars at $z > 5.8$ imaged by the Chandra X-ray observatory and found none with statistically significant multiple images. In the sub-sample of the 8 sources with photon counts $>20$ we exclude multiple images with separations $r>1''$ and count ratios $f>0.4$, or with separations as small as $0.''7$ and $f>0.7$ at $95\%$ confidence level. Comparing this non-detection with predictions from theoretical models suggesting a high and a low lensed fraction, we placed upper limits on the bright-end slope, $β$, of the QLF. Using only the sub-sample with 8 sources, we obtain, in the high-lensing model, a limit $β< 3.38$. Assuming no multiple source is present in the full sample of 22 sources, we obtain $β< 2.89$ and $β< 3.53$ in the high and low lensing models, respectively. These constraints strongly disfavor steep QLF shapes previously proposed in the literature.

preprint2022arXiv

Measuring the Expansion or Contraction of Galaxies

Galaxies lose mass as a result of their luminosity or gaseous outflows. I calculate the resulting radial migration of stars outwards and show that it could potentially be measured with high resolution spectrographs on the next generation of large telescopes. Substantial accretion of matter in dense cosmic environments could trigger inward stellar migration that would be even more easily measurable.

preprint2022arXiv

Merger rates of intermediate-mass black hole binaries in nuclear star clusters

Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravitational waves (GWs), which represent an unprecedented opportunity to reveal IMBHs. In this paper, we analyze the statistical results of a wide range of numerical simulations, which encompass different cluster metallicities, initial BH seed masses, and initial BH spins, and we compute the merger rate of IMBH binaries. We find that merger rates are in the range $0.01$-$10$\,Gpc$^{-3}$\,yr$^{-1}$ depending on IMBH masses. We also compute the number of multi-band detections in ground-based and space-based observatories. Our model predicts that a few merger events per year should be detectable with LISA, DECIGO, ET, and LIGO for IMBHs with masses $\lesssim 1000\msun$, and a few tens of merger events per year with DECIGO, ET, and LIGO only.

preprint2022arXiv

New Constraints on the Composition and Initial Speed of CNEOS 2014-01-08

We study the newly released light curve from the fireball of the first interstellar meteor CNEOS 2014-01-08. The measured velocity and three observed flares down to an altitude of $18.7 \mathrm{\; km}$ imply ambient ram pressure in the range of $113-194$ MPa when the meteor disintegrated. The required yield strength is $\gtrsim 20$ times higher than stony meteorites and $\gtrsim 2$ times larger than iron meteorites. The implied slowdown in the atmosphere suggests an initial speed of about $66.5 \; {\rm km~s^{-1}}$, strengthening the case for an interstellar origin of this meteor and making it an outlier relative to the velocity dispersion of local stars.

preprint2022arXiv

Self-Interacting Dark Matter from Gravitational Scattering

I show that gravitational scattering of dark-matter objects of 10^4 solar masses and speeds of 10 km/s, provides the cross-section per unit mass required in self-interacting dark matter models that alleviate the small-scale structure challenges to the collisionless cold dark matter model. For primordial objects of mass 10^4*(M_4) solar masses, moving at the velocity dispersion characteristic of dwarf galaxies, 10*(v_1) km/s, the cross-section per unit mass for gravitational scattering is 10*[M_4/(v_1)^4] cm^2/g. The steep decline in interaction with increasing velocity explains why self-interaction is not evident in data on massive galaxies and clusters of galaxies.

preprint2022arXiv

Signatures of Population III supernovae at Cosmic Dawn: the case of GN-z11-flash

We illustrate the observability of the end stages of the earliest (Population III) stars at high redshifts $z \gtrsim 10$, using the recently observed transient, GN-z11-flash as an example. We find that the observed spectrum of this transient is consistent with its originating from a shock-breakout in a Population III supernova occurring in the GN-z11 galaxy at $z \sim 11$. The energetics of the explosion indicate a progenitor star of mass $\sim 300 M_{\odot}$ in that galaxy, with of order unity such events expected over an observing timescale of a few years. We forecast the expected number of such transients from $z > 10$ galaxies as a function of their host stellar mass and star formation rate. Our findings are important in the context of future searches to detect and identify the signatures of galaxies at Cosmic Dawn.

preprint2022arXiv

Signatures of supermassive black hole binaries on maser systems

We illustrate a novel signature of black hole binaries from their effect on the kinematics of water maser emission in their environments. With the help of simulations, we establish the condition for clumps to mase based on their coherence lengths calibrated to those of the known maser galaxy NGC 4258. This is then used to identify masing clumps around a binary black hole system, and quantify the kinematic and spectral differences relative to the single black hole case. For some generic circumstances, blue-shifted masers around a binary black hole are found to preferentially follow the Keplerian rotation curve observed in the single black hole case. The redshifted ones, however, are found to visibly deviate from this relation, and also display more scatter with a tendency towards lower absolute values of the velocity along the line-of-sight. The spectrum of the masers as a function of line-of-sight velocity also shows a double peaked structure, reminiscent of recent observations of systems such as Mrk 1. Our results motivate future prospects for identifying binary black hole candidates with the help of water maser emissions.

preprint2022arXiv

Spurious Radial Migration from Relativistic Effects in the Milky-Way Disk

The gradient of the gravitational redshift in the potential of the Milky-Way induces an apparent spurious radial migration. I show that this effect is simply related to the local acceleration, which was measured recently by Gaia eDR3, implying a spectroscopic shift of $[-2.4x10^{-2}/(r/8kpc)] km/s/kpc. The transverse Doppler effect yields a comparable contribution. The spurious radial velocity from both relativistic effects amounts to crossing a major portion of the Milky-Way disk during the age of the universe, and must be corrected for in any future measurement of the actual radial migration of stars.

preprint2022arXiv

Terrestrial Impact from the Passage of the Solar System through a Cold Cloud a Few Million Years Ago

It is expected that as the Sun travels through the interstellar medium (ISM), there will be different filtration of Galactic Cosmic Rays (GCR) that affect Earth. The effect of GCR on Earth's atmosphere and climate is still uncertain. Although the interaction with molecular clouds was previously considered, the terrestrial impact of compact cold clouds was neglected. There is overwhelming geological evidence from 60Fe and 244Pu isotopes that Earth was in direct contact with the ISM 2 million years ago, and the local ISM is home to several nearby cold clouds. Here we show, with a state-of the art simulation that incorporate all the current knowledge about the heliosphere that if the solar system passed through a cloud such as Local Leo Cold Cloud, then the heliosphere which protects the solar system from interstellar particles, must have shrunk to a scale smaller than the Earth's orbit around the Sun (0.22). Using a magnetohydrodynamic simulation that includes charge exchange between neutral atoms and ions, we show that during the heliosphere shrinkage, Earth was exposed to a neutral hydrogen density of up to 3000cm-3. This could have had drastic effects on Earth's climate and potentially on human evolution at that time, as suggested by existing data.

preprint2022arXiv

The 2019 Discovery of a Meteor of Interstellar Origin

The earliest confirmed interstellar object, `Oumuamua, was discovered in the Solar System by Pan-STARRS in 2017, allowing for a calibration of the abundance of interstellar objects of its size $\sim 100\;$ m. This was followed by the discovery of Borisov, which allowed for a similar calibration of its size $\sim 0.4 - 1 \mathrm{\; km}$. One would expect a much higher abundance of significantly smaller interstellar objects, with some of them colliding with Earth frequently enough to be noticeable. Based on the CNEOS catalog of bolide events, we identify the $\sim 0.45$m meteor detected at 2014-01-08 17:05:34 UTC as originating from an unbound hyperbolic orbit with 99.999\% confidence. The U.S. Department of Defense has since verified that "the velocity estimate reported to NASA is sufficiently accurate to indicate an interstellar trajectory." We infer that the meteor had an asymptotic speed of $v_{\infty} \sim 42.1 \pm 5.5\; \mathrm{km \; s^{-1}}$ outside of the solar system. Its origin is approximately towards R.A. $49.4 \pm 4.1^{\circ}$ and declination $11.2 \pm 1.8^{\circ}$, implying that its initial velocity vector was $58\pm6\; \mathrm{km\;s^{-1}}$ away from the velocity of the Local Standard of Rest (LSR). Its high LSR speed implies a possible origin from the deep interior of a planetary system or a star in the thick disk of the Milky Way galaxy. The local number density of its population is $10^{6{^{+0.75}_{-1.5}}} \; \mathrm{AU^{-3}}$ or $9 \times 10^{21{^{+0.75}_{-1.5}}} \; \mathrm{pc^{-3}}$ (necessitating 0.2 -- 20 Earth masses of material to be ejected per local star). We show that the detections of CNEOS 2014-01-08, `Oumuamua, and Borisov collectively imply that the differential size distribution in good agreement with a collisional distribution, with a power-law slope is $q \sim 3.6 \pm 0.5$, where the quoted uncertainty corresponds to $2 σ$.

preprint2022arXiv

Time Evolution of the CMB Quadrupole

I show that the quadrupole of the Cosmic Microwave Background (CMB) evolves more rapidly than previously expected, as a result of the acceleration of the Sun towards the Galactic center. The acceleration, measured most recently by Gaia EDR3, implies a fractional change in the quadrupole of ~10^{-9} per year, an order of magnitude larger than expected from the evolution in the last scattering surface of the CMB.

preprint2022arXiv

Where to find over-massive brown dwarfs: new benchmark systems for binary evolution

Under the right conditions brown dwarfs that gain enough mass late in their lives to cross the hydrogen burning limit will not turn into low-mass stars, but rather remain essentially brown dwarf-like. While these objects, called either beige dwarfs or over-massive brown dwarfs, may exist in principle, it remains unclear exactly how they would form astrophysically. We show that accretion from AGB winds, aided by the wind Roche lobe overflow mechanism, is likely to produce a substantial population of observable overmassive brown dwarfs, though other mechanisms are still plausible. Specifically we predict that sun-like stars born with a massive brown dwarf companion on an orbit with a semi-major axis of order 10 AU will likely produce overmassive brown dwarfs, which may be found today as companions to the donor star's remnant white dwarf. The identification and characterization of such an object would produce unique constraints on binary evolution because there is a solid upper limit on the brown dwarf's initial mass.

preprint2021arXiv

Breakup of a Long-Period Comet as the Origin of the Dinosaur Extinction

The origin of the Chicxulub impactor, which is attributed as the cause of the K/T mass extinction event, is an unsolved puzzle. The background impact rates of main-belt asteroids and long-period comets have been previously dismissed as being too low to explain the Chicxulub impact event. Here, we show that a fraction of long-period comets are tidally disrupted after passing close to the Sun, each producing a collection of smaller fragments that cross the orbit of Earth. This population could increase the impact rate of long-period comets capable of producing Chicxulub impact events by an order of magnitude. This new rate would be consistent with the age of the Chicxulub impact crater, thereby providing a satisfactory explanation for the origin of the impactor. Our hypothesis explains the composition of the largest confirmed impact crater in Earth's history as well as the largest one within the last million years. It predicts a larger proportion of impactors with carbonaceous chondritic compositions than would be expected from meteorite falls of main-belt asteroids.

preprint2021arXiv

Detecting Black Hole Occultations by Stars with Space Interferometric Telescopes

We show that the occultation of Sagittarius A* by stars can be detected with space-based or space-ground very-long-baseline-interferometers (SVLBIs), with an expected event rate that is high due to relativistic precession. We compute the tell-tale signal of an occultation event, and describe methods to flag non-occultation events that can masquerade as the signal.

preprint2021arXiv

Determining the Composition of Relativistic Jets from Polarization Maps

We present a stationary, axisymmetric, self-similar semi-analytic model of magnetically dominated jet plasma based on force-free regions of a relativistic magnetohydrodynamic simulation. We use this model to illustrate how the composition of relativistic jet plasma can be determined, with special attention to the example of M87. In particular, we compute synthetic Stokes maps in e-e+p plasmas with various positron-to-proton ratios using synchrotron emission models scaling the partial pressure of electrons and positrons emitting at the observed frequency to the magnetic pressure, taking into account Faraday rotation and conversion. The lepton-dominated models produce bilaterally asymmetric radio intensity profiles with strong linear polarization and Stokes Q and U maps that are bilaterally asymmetric (but strongly up-down correlated) and antisymmetric (and sometimes up-down anticorrelated), respectively. The hadronic models produce more centrally brightened intensity and polarization maps. Circular polarization provides the cleanest observational tool for distinguishing the plasmas, as it increases outward from the jet core and central axis for highly ionic plasma, and is suppressed for pair dominated plasma. We find a measurable degree of circular polarization V/I of O(10e-3) for sub-equipartition hadronic jet plasmas. Our stationary model predicts that the intensity-normalized autocorrelation functions of Q and U increase and decrease with frequency, respectively. On the other hand, the autocorrelation of V is less sensitive to the frequency. Multi-band polarimetric observations could therefore be used as a novel probe of the composition of jet plasma.

preprint2021arXiv

Distinguishing AGN from starbursts as the origin of double peaked Lyman-Alpha Emitters in the reionization era

We discuss the possible origin of the double-peaked profiles recently observed in Lyman-Alpha Emitters (LAEs) at the epoch of reionization ($z \gtrsim 6.5$) from obscured active galactic nuclei (AGN). Combining the extent of the Lyman-$α$ near-zones estimated from the blue peak velocity offset in these galaxies, with the ionizing emissivity of quasars at $z \gtrsim 6$, we forecast the intrinsic UV and X-ray luminosities of the AGN needed to give rise to their double-peaked profiles. We also estimate the extent of the obscuration of the AGN by comparing their luminosities to those of similar quasar samples at these epochs. Future X-ray and radio observations, as well as those with the ${James \ Webb \ Space \ Telescope}$, will be valuable tools to test the AGN contribution to the intergalactic-scale ionization zones of high-redshift LAEs.

preprint2021arXiv

Implications of recoil kicks for black hole mergers from LIGO/Virgo catalogs

The first and second Gravitational Wave Transient Catalogs by the LIGO/Virgo Collaboration include $50$ confirmed merger events from the first, second, and first half of the third observational runs. We compute the distribution of recoil kicks imparted to the merger remnants and estimate their retention probability within various astrophysical environments as a function of the maximum progenitor spin ($χ_{\rm max}$), assuming that the LIGO/Virgo binary black hole (BBH) mergers were catalyzed by dynamical assembly in a dense star cluster. We find that the distributions of average recoil kicks are peaked at about $150$ km s$^{-1}$, $250$ km s$^{-1}$, $350$ km s$^{-1}$, $600$ km s$^{-1}$, for maximum progenitor spins of $0.1$, $0.3$, $0.5$, $0.8$, respectively. Only environments with escape speed $\gtrsim 100$ km s$^{-1}$, as found in galactic nuclear star clusters as well as in the most massive globular clusters and super star clusters, could efficiently retain the merger remnants of the LIGO/Virgo BBH population even for low progenitor spins ($χ_{\rm max}=0.1$). In the case of high progenitor spins ($χ_{\rm max}\gtrsim 0.5$), only the most massive nuclear star clusters can retain the merger products. We also show that the estimated values of the effective spin and of the remnant spin of GW170729, GW190412, GW190519, and GW190620 can be reproduced if their progenitors were moderately spinning ($χ_{\rm max}\gtrsim 0.3$), while for GW190517 if the progenitors were rapidly spinning ($χ_{\rm max}\gtrsim 0.8$). Alternatively, some of these events could be explained if at least one of the progenitors is already a second-generation BH, originated from a previous merger.

preprint2021arXiv

Positron Effects on Polarized Images and Spectra from Jet and Accretion Flow Models of M87* and Sgr A*

We study the effects of including a nonzero positron-to-electron fraction in emitting plasma on the polarized SEDs and sub-millimeter images of jet and accretion flow models for near-horizon emission from M87* and Sgr A*. For M87*, we consider a semi-analytic fit to the force-free plasma regions of a general relativistic magnetohydrodynamic jet simulation which we populate with power-law leptons with a constant electron-to-magnetic pressure ratio. For Sgr A*, we consider a standard self-similar radiatively inefficient accretion flow where the emission is predominantly from thermal leptons with a small fraction in a power-law tail. In both models, we fix the positron-to-electron ratio throughout the emission region. We generate polarized images and spectra from our models using the general-relativistic ray tracing and radiative transfer from GRTRANS. We find that a substantial positron fraction reduces the circular polarization fraction at infrared and higher frequencies. However, in sub-millimeter images higher positron fractions increase polarization fractions due to strong effects of Faraday conversion. We find a M87* jet model that best matches the available broadband total intensity and 230 GHz polarization data is a sub-equipartition, with positron fraction of $\simeq$ 10%. We show that jet models with significant positron fractions do not satisfy the polarimetric constraints at 230 GHz from the Event Horizon Telescope (EHT). Sgr A* models show similar trends in their polarization fractions with increasing pair fraction. Both models suggest that resolved, polarized EHT images are useful to constrain the presence of pairs at 230 GHz emitting regions of M87* and Sgr A*.

preprint2021arXiv

The Copernican Principle Rules Out BLC1 as a Technological Radio Signal from the Alpha Centauri System

Without evidence for occupying a special time or location, we should not assume that we inhabit privileged circumstances in the Universe. As a result, within the context of all Earth-like planets orbiting Sun-like stars, the origin of a technological civilization on Earth should be considered a single outcome of a random process. We show that in such a Copernican framework, which is inherently optimistic about the prevalence of life in the Universe, the likelihood of the nearest star system, Alpha Centauri, hosting a radio-transmitting civilization is $\sim 10^{-8}$. This rules out, \textit{a priori}, Breakthrough Listen Candidate 1 (BLC1) as a technological radio signal from the Alpha Centauri system, as such a scenario would violate the Copernican principle by about eight orders of magnitude. We also show that the Copernican principle is consistent with the vast majority of Fast Radio Bursts being natural in origin.

preprint2021arXiv

The Extended Habitable Epoch of the Universe for Liquids Other than Water

At high redshifts, the temperature of the cosmic microwave background (CMB) was higher than its value today. We explore the possibility that life may have arisen early because the higher CMB temperature would have supplied the requisite energy for the existence of different solvents on the surfaces of objects. At redshifts of $z \lesssim 70$, after the first stars are predicted to have formed, a number of molecules (but not water) might have existed in liquid form over intervals of $\sim 10$ Myr to $\sim 100$ Myr. We delineate the challenges and prospects for life in the high-redshift Universe, and assess the various candidates for alternative biochemistries in this context -- of the options considered herein, we conclude that ethane is probably the most promising contender.

preprint2021arXiv

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters are the dominant formation channel, we emphasize that many formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior.

preprint2021arXiv

The Search for the Farthest Quasar: Consequences for Black Hole Growth and Seed Models

The quest for high-redshift quasars has led to a series of record-breaking sources, with the current record holder at $z=7.642$. Here, we show how future detections of $z>8$ quasars impact the constraints on the parameters for black hole growth and seed models. Using broad flat priors on the growth parameters (Eddington ratio $\,f_{\rm Edd}$, duty cycle ${\cal D}$, seed mass $M_{\rm \bullet, seed}$ and radiative efficiency $ε$), we show that the large uncertainties in their determination decrease by a factor $\sim 5$ when a quasar's detection redshift goes from $z=9$ to $z=12$. In this high-redshift regime, $ε$ tends to the lowest value allowed, and the distribution for $M_{\rm \bullet, seed}$ peaks well inside the heavy seed domain. Remarkably, two quasars detected at $z > 7$ with low accretion rates (J1243+0100 and J0313-1806) already tighten the available parameter space, requiring $M_{\rm \bullet, seed} > 10^{3.5} \,{\rm M_\odot}$ and $ε< 0.1$. The radiative efficiency is a crucial unknown, with factor $\sim 2$ changes able to modify the predicted mass by $\sim 3$ orders of magnitude already at $z\sim 9$. The competing roles of inefficient accretion (decreasing $ε$) and black hole spin-up (increasing $ε$) significantly impact growth models. Finally, we suggest that yields currently predicted by upcoming quasar surveys (e.g., Euclid) will be instrumental for determining the most-likely seed mass regime. For example, assuming thin-disk accretion, a detection of a quasar with $M_\bullet \sim 10^{10} \,{\rm M_\odot}$ by $z\sim 9-10$ would exclude the entire parameter space available for light seeds and dramatically reduce the one for heavy seeds.

preprint2020arXiv

A Real-Time Search for Interstellar Impacts on the Moon

The discovery of `Oumuamua and CNEOS 2014-01-08 allowed for a calibration of the impact rate of interstellar objects. We propose a new telescope in lunar orbit to study in real-time interstellar meteoroid impacts and to serve as a laboratory for hypervelocity collisions. We show that a telescope with diameter $D_a \gtrsim 2 \mathrm{\; m}$ would be capable of detecting $\gtrsim 1$ interstellar meteoroid impacts (among hundreds of Solar System meteoroid impacts) per year. For each meteoroid, measurements of the reflected sunlight and shadow, as well as the impact's optical flash and crater, would allow for the determination of the 3D velocity, mass, density, and composition, as well as the radiative efficiency.

preprint2020arXiv

An Upper Limit on the Initial Temperature of the Radiation-Dominated Universe

Gravitational waves (GWs) are produced by colliding particles through the gravitational analogue of electromagnetic bremsstrahlung. We calculate the contribution of free-free emission in the radiation-dominated Universe to the stochastic GW background. We find that the energy density of the resulting GW radiation is heavily dependent on the number of elementary particles, $N_{\mathrm{tot}}$, and the maximum initial temperature, $T_{\mathrm{max}}$. We rule out $N_{\mathrm{tot}}\gtrsim N_{\mathrm{SM}}$ for $T_{\mathrm{max}}\sim T_{\mathrm{Planck}}\approx10^{19}$ GeV and $N_{\mathrm{tot}}\gtrsim10^{13}\times N_{\mathrm{SM}}$ for $T_{\mathrm{max}}\sim10^{16}$ GeV, where $N_{\mathrm{SM}}$ is the number of particles in the Standard Model. In the case of inflation, existing cosmological data constrain $T_{\mathrm{max}}\lesssim10^{16}$ GeV. However, alternative models to inflation such as bouncing cosmologies allow for $T_{\mathrm{max}}$ near $T_{\mathrm{Planck}}$. At the energy scales we are considering, the extra number of particles arise naturally in models of extra dimensions.

preprint2020arXiv

Can Planet Nine Be Detected Gravitationally by a Sub-Relativistic Spacecraft?

Planet 9 was proposed as an explanation for the clustering of orbits for some trans-Neuptunian objects. Recently, the use of a sub-relativistic spacecraft was proposed to indirectly probe Planet 9's gravitational influence. Here we study the effects of the drag and electromagnetic forces exerted on a sub-relativistic spacecraft by the interstellar medium (ISM) and compare these forces with the gravitational force induced by Planet 9. We find that the resulting noise due to density and magnetic fluctuations would dominate over Planet 9's gravitational signal at sub-relativistic speeds, $v\gtrsim 0.001~c$. We then identify the parameter space required to overcome the drag and magnetic noise from the ISM turbulence and enable the detection of Planet 9's gravity. Finally, we discuss practical strategies to mitigate the effect of the drag and electromagnetic forces.

preprint2020arXiv

Constraints on Aquatic Photosynthesis for Terrestrial Planets Around Other Stars

Aquatic photosynthesis plays a major role in carbon fixation and O$_2$ production on Earth. In this Letter, we analyze the prospects for oxygenic photosynthesis in aquatic environments on modern Earth-analogs around F-, G-, K- and M-type stars. Our analysis takes into account the spectral type of the host star, attenuation of light by aquatic organisms, and rates of respiration and photosynthesis. We study the compensation depth ($\mathcal{Z}_\mathrm{CO}$) and the critical depth ($\mathcal{Z}_\mathrm{CR}$), defined respectively as the locations where the net growth rates and vertically integrated net growth rates of photoautotrophs become zero. Our analysis suggests that $\mathcal{Z}_\mathrm{CO}$ declines by more than an order of magnitude as one moves from the habitable zones around Sun-like stars to late-type M-dwarfs, but $\mathcal{Z}_\mathrm{CR}$ decreases by only a modest amount ($\sim 40\%$). For M-dwarf exoplanets, we propose that the photosynthetic red edge may constitute a more robust biosignature of aquatic photosynthesis compared to atmospheric O$_2$.

preprint2020arXiv

Constraints on the abundance of $0.01\,c$ stellar engines in the Milky Way

Stellar engines are hypothesized megastructures that extract energy from the host star, typically with the purpose of generating thrust and accelerating the stellar system. We explore the maximum potential speeds that could be realizable by stellar engines, and determine that speeds up to $\sim 0.1\,c$ might perhaps be attainable under optimal conditions. In contrast, natural astrophysical phenomena in the Milky Way are very unlikely to produce such speeds. Hence, astrometric surveys of hypervelocity stars may be utilized to conduct commensal searches for high-speed stellar engines in the Milky Way. It may be possible to derive bounds on their abundance, but this requires certain assumptions regarding the spatiotemporal distribution of such engines, which are not guaranteed to be valid.

preprint2020arXiv

Destruction of molecular hydrogen ice and Implications for 1I/2017 U1 (`Oumuamua)

The first interstellar object observed in our solar system, 1I/2017 U1 (`Oumuamua), exhibited a number of peculiar properties, including extreme elongation and acceleration excess. Recently, \cite{Seligman:2020vb} proposed that the object was made out of molecular hydrogen (H$_{2}$) ice. The question is whether H$_2$ objects could survive their travel from the birth sites to the solar system. Here we study destruction processes of icy H$_2$ objects through their journey from giant molecular clouds (GMCs) to the interstellar medium (ISM) and the solar system, owing to interstellar radiation, gas and dust, and cosmic rays. We find that thermal sublimation due to heating by starlight can destroy `Oumuamua-size objects in less than 10 Myr. Thermal sublimation by collisional heating in GMCs could destroy H$_2$ objects of `Oumuamua-size before their escape into the ISM. Most importantly, the formation of icy grains rich in H$_2$ is unlikely to occur in dense environments because collisional heating raises the temperature of the icy grains, so that thermal sublimation rapidly destroys the H$_2$ mantle before grain growth.

preprint2020arXiv

Detectability of Thermal Emission from Sub-Relativistic Objects

We calculate the surface temperature and the resulting brightness of sub-relativistic objects moving through the Solar system due to collisional heating by gas and radiative heating by solar radiation. The thermal emission from objects of size $\gtrsim 100$ m and speed of $\gtrsim 0.1c$, can be detected by the upcoming {\it James Webb Space Telescope} out to a distance of $\sim 100$ au. Future surveys could therefore set interesting limits on the abundance of fast-moving interstellar objects or spacecraft.

preprint2020arXiv

Detecting Interstellar Objects Through Stellar Occultations

Stellar occultations have been used to search for Kuiper Belt and Oort Cloud objects. We propose a search for interstellar objects based on the characteristic durations ($\sim 0.1 \mathrm{s}$) of their stellar occultation signals and high inclination relative to the ecliptic plane. An all-sky monitoring program of all $\sim 7 \times 10^6$ stars with $R \lesssim 12.5$ using 1-m telescopes with $0.1 \; \mathrm{s}$ cadences is predicted to discover $\sim 1$ interstellar object per year.

preprint2020arXiv

Detecting the Orbital Motion of Nearby Supermassive Black Hole Binaries with Gaia

We show that a 10 year Gaia mission could astrometrically detect the orbital motion of ~1 sub-parsec separation supermassive black hole binary in the heart of nearby, bright active galactic nuclei (AGN). Candidate AGN lie out to a redshift of z=0.02 and in the V-band magnitude range $10 \lesssim m_V \lesssim 13$. The distribution of detectable binary masses peaks at a few times $10^7$ solar masses and is truncated above a few times $10^8$ solar masses.

preprint2020arXiv

Electromagnetic signals from the decay of free neutrons in the first hours of neutron star mergers

The first hours following a neutron star merger are considered to provide several UV/optical/NIR signals: $ β$-decay emission from free neutrons, radioactive decay of shocked heavy elements in the cocoon and cocoon's cooling emission. Here we consider two additional emission sources: $ β$-decay of free neutrons in the cocoon and synchrotron by the $ β$-decay electrons. We present 3D RHD simulations of jets that propagate in a multi-layer ejecta from the merger and calculate semi-analytically the resulting light curves. We find that the free neutrons emission at high latitudes is enhanced by the cocoon by a factor of a few to power a wide ($ \lesssim 60^\circ $) and brief ($ \sim 1 $ hour) UV signal that can reach an absolute magnitude of $\gtrsim$ -15, comparable with the cooling emission. If the ejected neutron matter mass is $ M_n \gtrsim 10^{-4}{\rm M_\odot} $, the synchrotron emission may yield a long ($ \sim 8 $ hours) quasi-isotropic UV/optical signal with an absolute magnitude between -12 and -15, depending on the magnetic field. Such a high mass of a mildly-relativistic component may partly obscure the cocoon's shocked r-process elements, thereby attenuating its radioactive decay emission. Future observations on these timescales, including null detections, may place constraints on the ejected neutron matter mass and shed light on the ejecta and jet-cocoon characteristics.

preprint2020arXiv

Eppur è piatto? The cosmic chronometer take on spatial curvature and cosmic concordance

The question of whether Cosmic Microwave Background (CMB) temperature and polarization data from Planck favor a spatially closed Universe with curvature parameter $Ω_K<0$ has been the subject of recent intense discussions. Attempts to break the geometrical degeneracy combining Planck data with external datasets such as Baryon Acoustic Oscillation (BAO) measurements all point towards a spatially flat Universe, at the cost of significant tensions with Planck, which make the resulting dataset combination problematic. Settling this issue requires identifying a dataset which can break the geometrical degeneracy while not incurring in these tensions. We argue that cosmic chronometers (CC), measurements of the expansion rate $H(z)$ from the relative ages of massive early-type passively evolving galaxies, are the dataset we are after. Furthermore, CC come with the additional advantage of being virtually free of cosmological model assumptions. Combining Planck 2018 CMB temperature and polarization data with the latest CC measurements, we break the geometrical degeneracy and find $Ω_K=-0.0054 \pm 0.0055$, consistent with a spatially flat Universe and competitive with the Planck+BAO constraint. Our results are stable against minimal parameter space extensions and CC systematics, and we find no substantial tension between Planck and CC data within a non-flat Universe, making the resulting combination reliable. Our results allow us to assert with confidence that the Universe is spatially flat to the ${\cal O}(10^{-2})$ level, a finding which might possibly settle the ongoing spatial curvature debate, and lends even more support to the already very successful inflationary paradigm.

preprint2020arXiv

Exploration of the outer solar system with fast and small sailcraft

Two new interplanetary technologies have advanced in the past decade to the point where they may enable exciting, affordable missions that reach further and faster deep into the outer regions of our solar system: (i) small and capable interplanetary spacecraft and (ii) light-driven sails. Combination of these two technologies could drastically reduce travel times within the solar system. We discuss a new paradigm that involves small and fast moving sailcraft that could enable exploration of distant regions of the solar system much sooner and faster than previously considered. We present some of the exciting science objectives for these miniaturized intelligent space systems that could lead to transformational advancements in the space sciences in the coming decade.

preprint2020arXiv

Formation of mass gap objects in highly asymmetric mergers

The LIGO/Virgo Collaboration (LVC) recently reported the detection of GW190814, a merger of a $23^{+1.0}_{-0.9}~M_{\odot}$ primary black hole (BH), and a $2.6^{+0.08}_{-0.08}~M_{\odot}$ secondary. The secondary's mass falls into the mass-gap regime, which refers to the scarcity of compact objects in the mass range of 2-5 $M_{\odot}$. The first clue to the formation of the GW190814 lies in the fact that the primary is a very massive BH. We suggest that the secondary was born as a neutron star (NS) where a significant amount of the supernova ejecta mass from its formation remained bound to the binary due to the presence of the massive BH companion. The bound mass forms a circumbinary accretion disk, and its accretion onto the NS created a mass-gap object. In this scenario, LIGO/Virgo will only detect mass-gap objects in binary mergers with an extreme mass ratio. We also predict a correlation between the mass of the secondary and the mass of the primary in such asymmetric mergers. Our model can be tested with future data from the LVC's third-observing run.

preprint2020arXiv

Gravitational-wave captures by intermediate-mass black holes in galactic nuclei

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.

preprint2020arXiv

Impacts of Dust Grains Accelerated by Supernovae on the Moon

There is evidence that ejecta from nearby supernovae have rained down on Earth in the past. Supernovae can accelerate pre-existing dust grains in the interstellar medium to speeds of $\sim 0.01 \mathrm{\;c}$. We investigate the survival and impact of dust grains from supernovae on the moon, finding that supernova dust grains can form detectable tracks with widths of $\sim 0.01 - 0.07 \mathrm{\; μm}$ and depths of $\sim 0.1 - 0.7 \mathrm{\; mm}$ in lunar rocks. These tracks could potentially shed light on the timings, luminosities, and directions of nearby supernovae.

preprint2020arXiv

Ionization and dissociation induced fragmentation of a tidally disrupted star into planets around a supermassive black hole

We show results from the radiation hydrodynamics (RHD) simulations of tidal disruption of a star on a parabolic orbit by a supermassive black hole (SMBH) based on a three-dimensional smoothed particle hydrodynamics code with radiative transfer. We find that such a tidally disrupted star fragment and form clumps soon after its tidal disruption. The fragmentation results from the endothermic processes of ionization and dissociation that reduce the gas pressure, leading to local gravitational collapse. Radiative cooling is less effective because the stellar debris is still highly optically thick in such an early time. Our simulations reveal that a solar-type star with a stellar density profile of n=3 disrupted by a 10^6 solar mass black hole produces $\sim20$ clumps of masses in the range of 0.1 to 12 Jupiter masses. The mass fallback rate decays with time, with pronounced spikes from early to late time. The spikes provide evidence for the clumps of the returning debris, while the clumps on the unbound debris can be potentially freely-floating planets and brown dwarfs. This ionization and dissociation induced fragmentation on a tidally disrupted star are a promising candidate mechanism to form low-mass stars to planets around an SMBH.

preprint2020arXiv

It is Feasible to Directly Measure Black Hole Masses in the First Galaxies

In the local universe, black hole masses have been inferred from the observed increase in the velocities of stars at the centres of their host galaxies. So far, masses of supermassive black holes in the early universe have only been inferred indirectly, using relationships calibrated to their locally observed counterparts. Here, we use the latest observational constraints on the evolution of stellar masses in galaxies to predict, for the first time, that the region of influence of a central supermassive black hole at the epochs where the first galaxies were formed is $\mathit{directly \ resolvable}$ by current and upcoming telescopes. We show that the existence of the black hole can be inferred from observations of the gas or stellar disc out to $> 0.5$ kpc from the host halo at redshifts $z \gtrsim 6$. Such measurements will usher in a new era of discoveries unraveling the formation of the first supermassive black holes based on subarcsecond-scale spectroscopy with the $JWST$, $ALMA$, and the $SKA$. The measured mass distribution of black holes will allow forecasting of the future detection of gravitational waves from the earliest black hole mergers.

preprint2020arXiv

Merging Black Holes in the Low-mass and High-mass Gaps from 2+2 Quadruple Systems

The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function.

preprint2020arXiv

New empirical constraints on the cosmological evolution of gas and stars in galaxies

We combine the latest observationally motivated constraints on stellar properties in dark matter haloes, along with data-driven predictions for the atomic (HI) and molecular (H$_2$) gas evolution in galaxies, to derive empirical relationships between the build-up of galactic components and their evolution over cosmic time. At high redshift ($z \gtrsim 4$), the frameworks imply that galaxies acquire their cold gas (both atomic and molecular) mostly by accretion, with the fraction of cold gas reaching about 20% of the cosmic baryon fraction. We infer a strong dependence of the star formation rate on the H$_2$ mass, suggesting a near-universal depletion timescale of 0.1-1 Gyr in Milky Way sized haloes (of masses $10^{12} \ M_{\odot}$ at $z = 0$). There is also evidence for a near-universality of the Kennicutt-Schmidt relation across redshifts, with very little dependence on stellar mass, if a constant conversion factor ($α_{\rm CO}$) of CO luminosity to molecular gas mass is assumed. Combining the atomic and molecular gas observations with the stellar build-up illustrates that galactic mass assembly in Milky-Way sized haloes proceeds from smooth accretion at high redshifts, towards becoming merger-dominated at late times ($z \lesssim 0.6$). Our results can be used to constrain numerical simulations of the dominant growth and accretion processes of galaxies over cosmic history.

preprint2020arXiv

Observational Signatures of Sub-Relativistic Meteors

It is currently unknown whether solid particles larger than dust from supernova ejecta rain down on Earth at high speeds. We develop a hydrodynamic and radiative model to explore the detectability of $\gtrsim 1 \mathrm{\; mm}$ sub-relativistic meteors. We find that a large fraction of the meteor energy during its passage through the Earth's upper atmosphere powers the adiabatic expansion of a hot plasma cylinder, giving rise to acoustic shocks detectable by infrasound microphones. Additionally, a global network of several hundred all-sky optical cameras with a time resolution of $\lesssim 10^{-4} \mathrm{\;s}$ would be capable of detecting $\gtrsim 1 \mathrm{\; mm}$ sub-relativistic meteors.

preprint2020arXiv

Observational signatures of the black hole Mass Distribution in the Galactic Center

We simulate the star cluster, made of stars in the main sequence and different black hole (BH) remnants, around SgrA* at the center of the Milky Way galaxy. Tracking stellar evolution, we find the BH remnant masses and construct the BH mass function. We sample 4 BH species and consider the impact of the mass-function in the dynamical evolution of system. Starting from an initial 6 dimensional family of parameters and using an MCMC approach, we find the best fits to various parameters of model by directly comparing the results of the simulations after $t = 10.5$ Gyrs with current observations of the stellar surface density, stellar mass profile and the mass of SgrA*. Using these parameters, we study the dynamical evolution of system in detail. We also explore the mass-growth of SgrA* due to tidally disrupted stars and swallowed BHs. We show that the consumed mass is dominated for the BH component with larger initial normalization as given by the BH mass-function. Assuming that about 10% of the tidally disrupted stars contribute in the growth of SgrA* mass, stars make up the second dominant effect in enhancing the mass of SgrA*. We consider the detectability of the GW signal from inspiralling stellar mass BHs around SgrA* with LISA. Computing the fraction of the lifetime of every BH species in the LISA band, with signal to noise ratio $\gtrsim 8$, to their entire lifetime, and rescaling this number with the total number of BHs in the system, we find that the total expected rate of inspirals per Milky-Way sized galaxy per year is $10^{-5}$. Quite interestingly, the rate is dominated for the BH component with larger initial normalization as dictated by the BH mass-function. We interpret it as the second signature of the BH mass-function.

preprint2020arXiv

On Spin Dependence of the Fundamental Plane of Black Hole Activity

The Fundamental Plane (FP) of Black Hole (BH) Activity in galactic nuclei relates X-ray and radio luminosities to BH mass and accretion rate. However, there is a large scatter exhibited by the data, which motivated us for a new variable. We add BH spin as a new variable and estimate the spin dependence of the jet power and disk luminosity in terms of radio and X-ray luminosities. We assume the Blandford-Znajek process as the main source of the outflow, and find that the jet power depends on BH spin stronger than quadratically at moderate and large spin values. We perform a statistical analysis for 10 AGNs which have sub-Eddington accretion rates and whose spin values are measured independently via the reflection or continuum-fitting methods, and find that the spin-dependent relation describes the data significantly better. This analysis, if supported with more data, could imply not only the spin dependence of the FP relation, but also the Blandford-Znajek process in AGN jets.

preprint2020arXiv

On the Habitable Lifetime of Terrestrial Worlds with High Radionuclide Abundances

The presence of a liquid solvent is widely regarded as an essential prerequisite for habitability. We investigate the conditions under which worlds outside the habitable zones of stars are capable of supporting liquid solvents on their surface over geologically significant timescales via combined radiogenic and primordial heat. Our analysis suggests that super-Earths with radionuclide abundances that are $\gtrsim 10^3$ times higher than Earth can host long-lived water oceans. In contrast, the requirements for long-lived ethane oceans, which have been explored in the context of alternative biochemistries, are less restrictive: relative radionuclide abundances of $\gtrsim 10^2$ could be sufficient. We find that this class of worlds might be detectable ($10σ$ detection over $\sim 10$ days integration time at $12.8$ $μ$m) in principle by the James Webb Space Telescope at distances of $\sim 10$ pc if their ages are $\lesssim 1$ Gyr.

preprint2020arXiv

Possible Transfer of Life by Earth-Grazing Objects to Exoplanetary Systems

Recently, a 30-cm object was discovered to graze the Earth's atmosphere and shift into a Jupiter-crossing orbit. We use the related survey parameters to calibrate the total number of such objects. The number of objects that could have exported terrestrial microbes out of the Solar System is in the range $2 \times 10^{9}$--$3 \times 10^{11}$. We find that $10^{7}$--$10^{9}$ such objects could have been captured by binary star systems over the lifetime of the Solar System. Adopting the fiducial assumption that one polyextremophile colony is picked up by each object, the total number of objects carrying living colonies on them upon capture could be $10$--$10^3$.

preprint2020arXiv

Potential For Liquid Water Biochemistry Deep Under The Surfaces Of The Moon, Mars And Beyond

We investigate the prospects for the past or current existence of habitable conditions deep underneath the surfaces of the Moon and Mars as well as generic bound and free-floating extrasolar rocky objects. We construct a simple model that takes into account the thermal limits of life as well as the size, surface temperature, and relative radionuclide abundance of a given object and yields the spatial extent of the subsurface habitable region. We also investigate the constraint imposed by pressure on habitability, and show that it is unlikely to rule out the prospects for life altogether. We estimate the maximum biomass that might be sustainable in deep subsurface environments as a function of the aforementioned parameters from an energetic perspective. We find that it might be a few percent that of Earth's subsurface biosphere, and three orders of magnitude smaller than Earth's global biomass, under ideal circumstances. We conclude with a brief exposition of the prevalence of rocky objects with deep biospheres and methods for detecting signatures of biological activity through forthcoming missions to visit the Moon and Mars.

preprint2020arXiv

Pre-Common-Envelope Mass Loss from Coalescing Binary Systems

Binary systems undergoing unstable Roche Lobe overflow spill gas into their circumbinary environment as their orbits decay toward coalescence. In this paper, we use a suite of hydrodynamic models of coalescing binaries involving an extended donor and a more compact accretor. We focus on the period of unstable Roche Lobe overflow that ends as the accretor plunges within the envelope of the donor at the onset of a common envelope phase. During this stage, mass is removed from the donor and flung into the circumbinary environment. Across a wide range of binary mass ratios, we find that the mass expelled as the separation decreases from the Roche limit to the donor's original radius is of the order of 25% of the accretor's mass. We study the kinematics of this ejecta and its dependencies on binary properties and find that it assembles into a toroidal circumbinary distribution. These circumbinary tori have approximately constant specific angular momentum due to momentum transport by spiral shocks launched from the orbiting binary. We show that an analytic model with these torus properties captures many of the main features of the azimuthally-averaged profiles of our hydrodynamic simulations. Our results, in particular the simple relationship between accretor mass and expelled mass and its spatial distribution, may be useful in interpreting stellar coalescence transients like luminous red novae, and in initializing hydrodynamic simulations of the subsequent common envelope phase.

preprint2020arXiv

Probing Intermediate Mass Black Holes in M87 through Multi-Wavelength Gravitational Wave Observations

We analyze triple systems composed of the super massive black hole (SMBH) near the center of M87 and a pair of black holes (BHs) with masses in the range $10-10^3$ $M_{\odot}$. We consider the post Newtonian precession as well as the Kozai-Lidov interactions at the quadruple and octupole levels in modeling the evolution of binary black hole (BBH) under the influence of the SMBH. Kozai-Lidov oscillations enhance the gravitational wave (GW) signal in some portions of the parameter space. We identify frequency peaks and examine the detectability of GWs with LISA as well as future observatories such as $μ$Ares and DECIGO. We show examples in which GW signal can be observed with a few or all of these detectors. Multi-Wavelength GW spectroscopy holds the potential to discover stellar to intermediate mass BHs near the center of M87. We estimate the rate, $Γ$, of collisions between the BBHs and flyby stars at the center of M87. Our calculation suggests $Γ< 10$ $\rm{Gyr}^{-1}$ for a wide range of the mass and semi-major axes of the inner binary.

preprint2020arXiv

Propulsion of Spacecrafts to Relativistic Speeds Using Natural Astrophysical Sources

In this paper, we explore from a conceptual standpoint the possibility of using natural astrophysical sources to accelerate spacecrafts to relativistic speeds. We focus on light sails and electric sails, which are reliant on momentum transfer from photons and protons, respectively, because these two classes of spacecrafts are not required to carry fuel on board. The payload is assumed to be stationed near the astrophysical source, and the sail is subsequently unfolded and activated when the source is functional. By considering a number of astrophysical objects such as massive stars, microquasars, supernovae, pulsar wind nebulae, and active galactic nuclei, we show that terminal speeds approaching the speed of light might be realizable under idealized circumstances provided that sufficiently advanced sail materials and control techniques exist. We also investigate the constraints arising from the sail's material properties, the voyage through the ambient source environment, and the passage through the interstellar medium. While all of these considerations pose significant challenges to spacecrafts, our analysis indicates that they are not insurmountable in optimal conditions. Finally, we sketch the implications for carrying out future technosignature searches.

preprint2020arXiv

Prospects for Life on Temperate Planets Around Brown Dwarfs

There is growing evidence that brown dwarfs may be comparable to main-sequence stars in terms of their abundance. In this paper, we explore the prospects for the existence of life on Earth-like planets around brown dwarfs. We consider the following factors: (i) the length of time that planets can exist in the temporally shifting habitable zone, (ii) the minimum photon fluxes necessary for oxygenic photosynthesis, and (iii) the lower limits on the fluxes of ultraviolet radiation to drive prebiotic reactions ostensibly necessary for the origin of life. By taking these effects into consideration, we find that it is unlikely for brown dwarfs with masses $\lesssim 30 M_J$ to host habitable planets over geologically significant timescales. We also briefly discuss some of the major biosignatures that might arise on these planets, assess the likelihood of their detection, and highlight some avenues for further study.

preprint2020arXiv

Reality or Mirage? Observational Test and Implications for the Claimed Extremely Magnified Quasar at $z = 6.3$

In the last two decades, approximately 200 quasars have been discovered at $z>6$, hosting active super-massive black holes with masses $M_{\bullet} \gtrsim 10^9 M_{\odot}$. While these sources reflect only the tip of the iceberg of the black hole mass distribution, their detection challenges standard growth models. The most massive $z>6$ black hole that was inferred thus far (J0100+2802, $M_{\bullet} \approx 1.2\times 10^{10} M_{\odot}$) was recently claimed to be lensed, with a magnification factor $μ=450$. Here, we perform a consistency check of this claim, finding that the detection of such source requires a bright-end slope $β\geq 3.7$ for the intrinsic quasar luminosity function, $Φ(L) \propto L^{-β}$. Commonly used values of $β\sim 2.8$ are rejected at $>3σ$. If the claim is confirmed, it is very unlikely that all the remaining $51$ sources in the SDSS sample are not magnified. Furthermore, it suffices that $\gtrsim 25\%$ of the remaining sources are lensed for the intrinsic luminosity function to differ significantly (i.e., $>3σ$) from the observed one. The presence of additional extremely magnified sources in the sample would lower the requirement to $\sim 4\%$. Our results urge the community to perform more extended multi-wavelength searches targeting $z>6$ lensed quasars, also among known samples. This effort could vitally contribute to solving the open problem of the growth of the brightest $z\sim 7$ quasars.

preprint2020arXiv

Repeated Gravitational Lensing of Gravitational Waves in Hierarchical Black Hole Triples

We consider binary black holes (BBHs) in a hierarchical triple system where a more compact, less-massive binary is emitting detectable gravitational waves (GWs), and the tertiary is a supermassive BH at the center of a nuclear star cluster. As previous works have shown, the orbital motion of the outer binary can generate a detectable relativistic Doppler boost of the GWs emitted by the orbiting inner binary. We show here that for outer-binary orbits with a period of order one year, there can be a non-negligible probability for repeated gravitational lensing of the GWs emitted by the inner binary. Repeating gravitational lensing events could be detected by the LISA observatory as periodic GW amplitude spikes before the BBH enters the LIGO band. Such a detection would confirm the origin of some BBH mergers in nuclear star clusters. GW lensing also offers new testing grounds for strong gravity.

preprint2020arXiv

Repeated Impact-Driven Plume Formation On Enceladus Over Megayear Timescales

Water plumes erupting from the `tiger stripe' features on the south pole of Enceladus are thought to connect to a global subsurface ocean. Proposed origins for the initial stress necessary to form the `tiger stripes' include a giant impact, which would require true polar wander to explain the location of the plumes if the impact did not occur at the South Pole, or tensile stresses, which would require volumetric expansion associated with partial freezing of the subsurface ocean. In addition to the requirement of a fine-tuned impact, true polar wander, or partial freezing of the subsurface ocean, a further issue with these hypotheses is that the `tiger stripes' may be short-lived. We show here that impact resurfacing can seal off plumes and mass loss can lead to their compression and closure over $\sim 1 \mathrm{\;Myr}$. Since plumes are observed at present, a mechanism by which new plumes can be generated every $\sim 1 \mathrm{\;Myr}$ and by which such plumes are most likely to form at the south pole is needed. We propose and investigate the possibility that impacts constitute an adequate repeating source for the continual instigation of fractures and plumes. We find that the rate of impacts on Enceladus suggests the formation of $\sim 10^3$ independent plume systems per Gyr, the vast majority on the south pole, and is consistent with the Cassini-derived age of the south pole for a lunar-like bombardment history, our estimates of fracture lifetimes, and with the needed parameters for parallel fracture propagation. The model favors a bombardment history similar to that of Triton over one more similar to that of the Galilean satellites.

preprint2020arXiv

Runaway Coalescence of Pre-Common-Envelope Stellar Binaries

We study the process of runaway, unstable Roche lobe overflow in coalescing binary systems and its dependence on the properties of the binary involved. We create three-dimensional hydrodynamic models of binary coalescences, and follow them through a phase of increasing Roche lobe overflow until the accretor is engulfed by the donor at the onset of a common envelope phase. In these models, we vary binary properties of mass ratio, donor structure and spin, and equation of state through the gas adiabatic index. We compare the numerical results to semi-analytic models of binary orbit evolution based on mass and angular momentum exchange between two point masses. Using our hydrodynamic simulations, we measure the key parameters: the donor mass loss rate and the angular momentum exchanged per unit mass loss from the donor. Using these calibrations, the semi-analytic model closely reproduces the escalating mass loss and runaway orbital decay observed in the hydrodynamic models. The semi-analytic model accurately reproduces the major differences in orbit evolution that arise with varying mass ratio and donor structure. We encapsulate the semi-analytic model in a publicly-released python package RLOF. We apply this model to the observed period decay and subsequent merger of the binary V1309 Sco, and find that it can simultaneously reproduce the observed orbital decay and time of outburst. We further demonstrate that there is a relationship between period derivative and second derivative that can be a useful metric for evaluating candidate merging binaries.

preprint2020arXiv

Searching for Black Holes in the Outer Solar System with LSST

Planet Nine has been proposed to potentially be a black hole in the outer solar system. We investigate the accretion flares that would result from impacts of small Oort cloud objects, and find that the upcoming LSST observing program will be able to either rule out or confirm Planet Nine as a black hole within a year. We also find that LSST could rule out or confirm the existence of trapped planet-mass black holes out to the edge of the Oort cloud, indirectly probing the dark matter fraction in subsolar mass black holes and potentially improving upon current limits by orders of magnitude.

preprint2020arXiv

Searching for exotic cores with binary neutron star inspirals

We study the feasibility of detecting exotic cores in merging neutron stars with ground-based gravitational-wave detectors. We focus on models with a sharp nuclear/exotic matter interface, and assume a uniform distribution of neutron stars in the mass range $[1,2] M_\odot$. We find that the existence of exotic cores can be confirmed at the 70% confidence level with as few as several tens of detections. Likewise, with such a sample, we find that some models of exotic cores can be excluded {with high confidence}.

preprint2020arXiv

Separating Accretion and Mergers in the Cosmic Growth of Black Holes with X-ray and Gravitational Wave Observations

Black holes across a broad range of masses play a key role in the evolution of galaxies. The initial seeds of black holes formed at $z \sim 30$ and grew over cosmic time by gas accretion and mergers. Using observational data for quasars and theoretical models for the hierarchical assembly of dark matter halos, we study the relative importance of gas accretion and mergers for black hole growth, as a function of redshift ($0<z<10$) and black hole mass ($ 10^3 \, \mathrm{M_{\odot}} < M_{\bullet} < 10^{10} \, \mathrm{M_{\odot}}$). We find that (i) growth by accretion is dominant in a large fraction of the parameter space, especially at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z>6$; and (ii) growth by mergers is dominant at $M_{\bullet} < 10^5 \, \mathrm{M_{\odot}}$ and $z>5.5$, and at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z<2$. As the growth channel has direct implications for the black hole spin (with gas accretion leading to higher spin values), we test our model against $\sim 20$ robust spin measurements available thus far. As expected, the spin tends to decline toward the merger-dominated regime, thereby supporting our model. The next generation of X-ray and gravitational-wave observatories (e.g. Lynx, AXIS, Athena and LISA) will map out populations of black holes up to very high redshift ($z\sim 20)$, covering the parameter space investigated here in almost its entirety. Their data will be instrumental to providing a clear picture of how black holes grew across cosmic time.

preprint2020arXiv

The nearest discovered black hole is likely not in a triple configuration

HR 6819 was recently claimed to be a hierarchical triple system of a Be star in a wide orbit around an inner binary system of a black hole (BH) and a B III type star. We argue that this system is unlikely to be a hierarchical triple due to three reasons: (i) Given that this system is discovered in a magnitude limited Bright Star Catalog, the expected number of such systems in the Milky Way amounts to about $10^4$ while the estimate for the MW budget for such systems is between $10^2-10^3$ systems under generous assumptions. Such a large gap cannot be reconciled as it would otherwise likely overflow the MW budget for BHs; (ii) The dynamical stability of this system sets lower bounds on the orbital separation of the outer Be star, while it not being resolved by Gaia places an upper limit on its projected sky separation. We show that these two constraints would imply a narrow range for the outer orbit without resorting to geometrical fine-tuning; (iii) The triple system should have survived the stellar evolution prior to the formation of the BH in the inner binary. We perform numerical simulations starting with conservative initial conditions of this system and show that a small parameter space for BH progenitor star's mass loss, BH natal kicks, and initial orbital separation can reproduce HR 6819. Therefore, we propose this system is a chance superposition of a Be star with a binary.

preprint2020arXiv

Transfer of Life Between Earth and Venus with Planet-Grazing Asteroids

Recently, phosphine was discovered in the atmosphere of Venus as a potential biosignature. This raises the question: if Venusian life exists, could it be related to terrestrial life? Based on the known rate of meteoroid impacts on Earth, we show that at least $\sim 6 \times 10^5$ asteroids have grazed Earth's atmosphere without being significantly heated and later impacted Venus, and a similar number have grazed Venus's atmosphere and later impacted the Earth, both within a period of $\sim 10^5$ years during which microbes could survive in space. Although the abundance of terrestrial life in the upper atmosphere is unknown, these planet-grazing shepherds could have potentially been capable of transferring microbial life between the atmospheres of Earth and Venus. As a result, the origin of possible Venusian life may be fundamentally indistinguishable from that of terrestrial life.

preprint2020arXiv

Upper Limit on the Dissipation of Gravitational Waves in Gravitationally Bound Systems

It is shown that a gravitationally bound system with a one-dimensional velocity dispersion sigma can at most dissipate a fraction ~36(sigma/c)^3 of the gravitational wave energy propagating through it, even if their dynamical time is shorter than the wave period. The limit is saturated for particles with a mean-free-path equal to the size of the system, such as hot protons in galaxy clusters, strongly-interacting dark matter particles in halos, or massive black holes in clusters. For such systems with random motions and no resonances, the dissipated fraction, <10^{-6}, does not degrade the use of gravitational waves as cosmological probes. At high wave frequencies, the dissipated fraction is additionally suppressed by the square of the ratio between the collision frequency and the wave frequency. The electromagnetic counterparts that result from the dissipation are too faint to be detectable at cosmological distances.

preprint2020arXiv

What's in a name: The etymology of astrobiology

Astrobiology has been gaining increasing scientific prominence and public attention as the search for life beyond Earth continues to make significant headway on multiple fronts. In view of these recent developments, the fascinating and dynamic etymology of astrobiology is elucidated, and thus shown to encompass a plethora of vivid characters drawn from different continents, religions, ideologies and centuries.

preprint2019arXiv

Explosions Driven by the Coalescence of a Compact Object with the Core of a Massive-Star Companion Inside a Common Envelope: Circumstellar Properties, Light Curves, and Population Statistics

We model explosions driven by the coalescence of a black hole or neutron star with the core of its massive-star companion. Upon entering a common envelope phase, a compact object may spiral all the way to the core. The concurrent release of energy is likely to be deposited into the surrounding common envelope, powering a merger-driven explosion. We use hydrodynamic models of binary coalescence to model the common envelope density distribution at the time of coalescence. We find toroidal profiles of material, concentrated in the binary's equatorial plane and extending to many times the massive star's original radius. We use the spherically-averaged properties of this circumstellar material (CSM) to estimate the emergent light curves that result from the interaction between the blast wave and the CSM. We find that typical merger-driven explosions are brightened by up to three magnitudes by CSM interaction. From population synthesis models we discover that the brightest merger-driven explosions, $M_V \sim -18$ to $-19$, are those involving black holes because they have the most massive and extended CSM. Black hole coalescence events are also common; they represent about 50% of all merger-driven explosions and approximately 0.3% of the core-collapse rate. Merger-driven explosions offer a window into the highly-uncertain physics of common envelope interactions in binary systems by probing the properties of systems that merge rather than eject their envelopes.

preprint2019arXiv

Exporting Terrestrial Life Out of the Solar System with Gravitational Slingshots of Earthgrazing Bodies

Exporting terrestrial life out of the Solar System requires a process that both embeds microbes in boulders and ejects those boulders out of the Solar System. We explore the possibility that Earthgrazing long-period comets and interstellar objects could export life from Earth by collecting microbes from the atmosphere and receiving a gravitational slingshot effect from the Earth. We estimate the total number of exportation events over the lifetime of the Earth to be $\sim 1 - 10$ for long-period comets and $\sim 1 - 50$ for interstellar objects. If life existed above an altitude of 100 km, then the number is dramatically increased up to $\sim 10^5$ exportation events over Earth's lifetime.

preprint2019arXiv

Formation and merging of Mass Gap Black Holes in Gravitational Wave Merger Events from Wide Hierarchical Quadruple Systems

We investigate secular evolution in hierarchical quadruple systems as a formation channel of mass-gap black holes (with masses of about $3-5$ $M_{\odot}$) in systems that will eventually lead to binary black hole mergers detectable by ground-based gravitational wave detectors (LIGO/Virgo). We show that in a 3+1 hierarchical system, two episodes of induced mergers would first cause two neutron stars to merge and form a mass-gap black hole, which will subsequently merge with another (more massive) black hole through a second induced merger. We demonstrate that such systems are stable to flybys, and their formation would predict a high mass ratio and eccentric merger of a mass-gap black hole with a more massive black hole companion. Such a formation channel may explain observed gravitational wave events such as the recently-discovered LIGO/Virgo events S190814bv and S190924h.

preprint2019arXiv

Global Stellar Budget for LIGO Black Holes

The binary black hole mergers observed by LIGO-Virgo gravitational-wave detectors pose two major challenges: (i) how to produce these massive black holes from stellar processes; and (ii) how to bring them close enough to merge within the age of the universe? We derive a fundamental constraint relating the binary separation and the available stellar budget in the universe to produce the observed black hole mergers. We find that $\lesssim 14\%$ of the entire budget contributes to the observed merger rate of $(30 + 30)~\mathrm{M}_\odot$ black holes, if the separation is around the diameter of their progenitor stars. Furthermore, the upgraded LIGO detector and third-generation gravitational-wave detectors are not expected to find stellar-mass black hole mergers at high redshifts. From LIGO's strong constraints on the mergers of black holes in the pair-instability mass-gap ($60-120~\mathrm{M}_\odot$), we find that $\lesssim 0.8\%$ of all massive stars contribute to a remnant black hole population in this gap. Our derived separation$-$budget constraint provides a robust framework for testing the formation scenarios of stellar binary black holes.

preprint2019arXiv

Periodic Fast Radio Bursts from Young Neutron Stars

Fast radio bursts (FRBs) are highly energetic radio pulses from cosmological origins. Despite an abundance of detections, their nature remains elusive. At least a subset of FRBs is expected to repeat, as the daily FRB rate surpasses that of any known cataclysmic event, which has been confirmed by observations. One of the proposed mechanisms to generate repeating FRBs is supergiant pulses from young and highly spinning NSs, in which case FRBs could inherit the periodicity of their parent NS. Here we examine the consequences of such a population of periodic fast radio bursts (PFRBs). We calculate the rate and lifetime of PFRB progenitors, and find that each newly born highly spinning NS has to emit a number $N_{\rm PFRB}\sim 10^2$ of bursts during its active lifetime of $τ\sim 100$ years, after which it becomes too dim and crosses a PFRB "death line" analogous to the pulsar one. We propose several tests of this hypothesis. First, the period of PFRBs would increase over time, and their luminosity would decrease, due to the NS spin-down. Second, PFRBs may show modest amounts of rotation measure, given the lack of expelled matter from the pulsar, as opposed to the magnetar-sourced FRBs proposed to explain the first repeater FRB 121102. As an example, we study whether the second confirmed repeater (FRB 180814) is a PFRB, given the preference for an inter-pulse separation of 13 ms within its sub-bursts. We show that, if confirmed, this period would place FRB 180814 in a different category as FRB 121102. We develop tests that would identify---and characterize---the prospective population of PFRBs.

preprint2019arXiv

Photosynthesis on Exoplanets and Exomoons from Reflected Light

Photosynthesis offers a convenient means of sustaining biospheres. We quantify the constraints for photosynthesis to be functional on the permanent nightside of tidally locked rocky exoplanets via reflected light from their exomoons. We show that the exomoons must be at least half the size of Earth's moon in order for conventional oxygenic photosynthesis to operate. This scenario of photosynthesis is unlikely for exoplanets around late-type M-dwarfs due to the low likelihood of large exomoons and their orbital instability over long timescales. Subsequently, we investigate the prospects for photosynthesis on habitable exomoons via reflected light from the giant planets that they orbit. Our analysis indicates that such photosynthetic biospheres are potentially sustainable on these moons except those around late-type M-dwarfs. We conclude our analysis by delineating certain physiological and biochemical features of photosynthesis and other carbon fixation pathways, and the likelihood of their evolution on habitable planets and moons.

preprint2019arXiv

Probing Extrasolar Planetary Systems with Interstellar Meteors

The first interstellar object, `Oumuamua, was discovered in the Solar System by Pan-STARRS in 2017, allowing for a calibration of the impact rate of interstellar meteors of its size $\sim 100\;$m. The discovery of CNEOS 2014-01-08 allowed for a calibration of the impact rate of interstellar meteors of its size $\sim 1\;$m. Analysis of interstellar dust grains have allowed for calibrations of the impact rate of smaller interstellar meteors down to the size $\sim 10^{-8}\;$m. We analyze the size distribution of interstellar meteors, finding that for smooth power-law fits of the form $N(r)\propto r^{-q}$, the possible values of $q$ are in the range $3.41 \pm 0.17$. We then consider the possibility of analyzing interstellar meteors to learn about their parent planetary systems. We propose a strategy for determining the orbits and chemical compositions of interstellar meteors, using a network of $\sim 600$ all-sky camera systems to track and conduct remote spectroscopy on meteors larger than $\sim 5$cm once every few years. It should also be possible to retrieve meteorites from the impact sites, providing the first samples of materials from other planetary systems.

preprint2016arXiv

3-cm Fine Structure Masers: A Unique Signature of Supermassive Black Hole Formation via Direct Collapse in the Early Universe

The direct collapse black hole (DCBH) scenario describes the isothermal collapse of a pristine gas cloud directly into a massive, M_BH=10^4-10^6 M_sun black hole. In this paper we show that large HI column densities of primordial gas at T~10^4 K with low molecular abundance - which represent key aspects of the DCBH scenario - provide optimal conditions for pumping of the 2p-level of atomic hydrogen by trapped Lyman alpha (Lya) photons. This Lya pumping mechanism gives rise to inverted level population of the 2s_1/2-2p_3/2 transition, and therefore to stimulated fine structure emission at 3.04 cm (rest-frame). We show that simplified models of the DCBH scenario amplify the CMB by up to a factor of 10^5, above which the maser saturates. Hyperfine splitting of the 3-cm transition gives rise to a characteristic broad (FWHM ~ tens of MHz in the observers frame) asymmetric line profile. This signal subtends an angular scale of ~ 1-10 mas, which translates to a flux of ~ 0.3-3 microJy, which is detectable with ultra-deep surveys being planned with SKA1-MID. While challenging, as the signal is visible for a fraction of the collapse time of the cloud, the matching required physical conditions imply that a detection of the redshifted 3-cm emission line would provide direct evidence for the DCBH scenario.

preprint2016arXiv

Blackbody Radiation from Isolated Neptunes

Recent analyses of the orbits of some Kuiper Belt objects hypothesize the presence of an undiscovered Neptune-size planet at a very large separation from the Sun. The energy budget of Neptunes on such distant orbits is dominated by the internal heat released by their cooling rather than solar irradiation (making them effectively "isolated"). The blackbody radiation that these planets emit as they cool may provide the means for their detection. Here we use an analytical toy model to study the cooling and radiation of isolated Neptunes. This model can translate a detection (or a null detection) to a constraint on the size and composition of the hypothesized "Planet Nine". Specifically, the thick gas atmosphere of Neptune-like planets serves as an insulating blanket which slows down their cooling. Therefore, a measurement of the blackbody temperature, $T_{\rm eff}\sim 50\textrm{K}$, at which a Neptune emits can be used to estimate the mass of its atmosphere, $M_{\rm atm}$. Explicitly, we find the relation $T_{\rm eff}\propto M_{\rm atm}^{1/12}$. Despite this weak relation, a measurement of the flux at the Wien tail can constrain the atmospheric mass, at least to within a factor of a few, and provide useful limits to possible formation scenarios of these planets. Finally, we constrain the size and composition of Planet Nine by combining our model with the null results of recent all-sky surveys.

preprint2016arXiv

CEMP stars: possible hosts to carbon planets in the early universe

We explore the possibility of planet formation in the carbon-rich protoplanetary disks of carbon-enhanced metal-poor (CEMP) stars, possible relics of the early Universe. The chemically anomalous abundance patterns ([C/Fe] $\geq$ 0.7) in this subset of low-mass stars suggest pollution by primordial core-collapsing supernovae (SNe) ejecta that are particularly rich in carbon dust grains. By comparing the dust-settling timescale in the protoplanetary disks of CEMP stars to the expected disk lifetime (assuming dissipation via photoevaporation), we determine the maximum distance $r_{max}$ from the host CEMP star at which carbon-rich planetesimal formation is possible, as a function of the host star's [C/H] abundance. We then use our linear relation between $r_{max}$ and [C/H], along with the theoretical mass-radius relation derived for a solid, pure carbon planet, to characterize potential planetary transits across host CEMP stars. Given that the related transits are detectable with current and upcoming space-based transit surveys, we suggest initiating an observational program to search for carbon planets around CEMP stars in hopes of shedding light on the question of how early planetary systems may have formed after the Big Bang.

preprint2016arXiv

Circularization of Tidally Disrupted Stars around Spinning Supermassive Black Holes

We study the circularization of tidally disrupted stars on bound orbits around spinning supermassive black holes by performing three-dimensional smoothed particle hydrodynamic simulations with Post-Newtonian corrections. Our simulations reveal that debris circularization depends sensitively on the efficiency of radiative cooling. There are two stages in debris circularization if radiative cooling is inefficient: first, the stellar debris streams self-intersect due to relativistic apsidal precession; shocks at the intersection points thermalize orbital energy and the debris forms a geometrically thick, ring-like structure around the black hole. The ring rapidly spreads via viscous diffusion, leading to the formation of a geometrically thick accretion disk. In contrast, if radiative cooling is efficient, the stellar debris circularizes due to self-intersection shocks and forms a geometrically thin ring-like structure. In this case, the dissipated energy can be emitted during debris circularization as a precursor to the subsequent tidal disruption flare. The possible radiated energy is up to ~2*10^{52} erg for a 1 Msun star orbiting a 10^6 Msun black hole. We also find that a retrograde (prograde) black hole spin causes the shock-induced circularization timescale to be shorter (longer) than that of a non-spinning black hole in both cooling cases. The circularization timescale is remarkably long in the radiatively efficient cooling case, and is also sensitive to black hole spin. Specifically, Lense-Thirring torques cause dynamically important nodal precession, which significantly delays debris circularization. On the other hand, nodal precession is too slow to produce observable signatures in the radiatively inefficient case. We also discuss the relationship between our simulations and the parabolic TDEs that are characteristic of most stellar tidal disruptions.

preprint2016arXiv

Constraining the CMB Optical Depth Through the Dispersion Measure of Cosmological Radio Transients

The dispersion measure of extragalactic radio transients, such as of recently discovered Fast Radio Burst FRB150418, can be used to measure the column density of free electrons in the intergalactic medium. The same electrons also scatter the Cosmic Microwave Background (CMB) photons, affecting precision measurements of cosmological parameters. We explore the connection between the dispersion measure of radio transients existing during the Epoch of Reionization (EoR) and the total optical depth for the CMB, $τ_{CMB}$, showing that the existence of such transients would provide a new sensitive probe of $τ_{CMB}$. As an example, we consider the population of FRBs. Assuming they exist during the EoR, we show that: (i) such sources can probe the reionization history by measuring $τ_{CMB}$ to sub-percent accuracy, and (ii) they can be detected with high significance by an instrument such as the Square Kilometer Array.

preprint2016arXiv

Cumulative neutrino background from quasar-driven outflows

Quasar-driven outflows naturally account for the missing component of the extragalactic $γ$-ray background through neutral pion production in interactions between protons accelerated by the forward outflow shock and interstellar protons. We study the simultaneous neutrino emission by the same protons. We adopt outflow parameters that best fit the extragalactic $γ$-ray background data and derive a cumulative neutrino background of $\sim10^{-7}\,\rm GeV\,cm^{-2}\,s^{-1}\,sr^{-1}$ at neutrino energies $E_ν\gtrsim 10$ TeV, which naturally explains the most recent IceCube data without tuning any free parameters. The link between the $γ$-ray and neutrino emission from quasar outflows can be used to constrain the high-energy physics of strong shocks at cosmological distances.

preprint2016arXiv

Detection of Gravitational Wave Emission by Supermassive Black Hole Binaries Through Tidal Disruption Flares

Galaxy mergers produce supermassive black hole binaries, which emit gravitational waves prior to their coalescence. We perform three-dimensional hydrodynamic simulations to study the tidal disruption of stars by such a binary in the final centuries of its life. We find that the gas stream of the stellar debris moves chaotically in the binary potential and forms accretion disks around both black holes. The accretion light curve is modulated over the binary orbital period owing to relativistic beaming. This periodic signal allows to detect the decay of the binary orbit due to gravitational wave emission by observing two tidal disruption events that are separated by more than a decade.

preprint2016arXiv

Electromagnetic Counterparts to Black Hole Mergers Detected by LIGO

Mergers of stellar-mass black holes (BHs), such as GW150914 observed by LIGO, are not expected to have electromagnetic counterparts. However, the Fermi GBM detector identified of a gamma-ray transient 0.4 s after the gravitational wave (GW) signal GW150914 with consistent sky localization. I show that the two signals might be related if the BH binary detected by LIGO originated from two clumps in a dumbbell configuration that formed when the core of a rapidly rotating massive star collapsed. In that case, the BH binary merger was followed by a gamma-ray burst (GRB) from a jet that originated in the accretion flow around the remnant BH. A future detection of a GRB afterglow could be used to determine the redshift and precise localization of the source. A population of standard GW sirens with GRB redshifts would provide a new approach for precise measurements of cosmological distances as a function of redshift.

preprint2016arXiv

Evidence for a direct collapse black hole in the Lyman-alpha source CR7

Throughout the epoch of reionization the most luminous Lyα emitters are capable of ionizing their own local bubbles. The CR7 galaxy at $z \approx 6.6$ stands out for its combination of exceptionally bright Lyα and HeII 1640 Angstrom line emission but absence of metal lines. As a result CR7 may be the first viable candidate host of a young primordial starburst or direct collapse black hole. High-resolution spectroscopy reveals a +160 km s$^{-1}$ velocity offset between the Lyα and HeII line peaks while the spatial extent of the Lyα emitting region is $\sim 16$ kpc. The observables are indicative of an outflow signature produced by a strong central source. We present one-dimensional radiation-hydrodynamics simulations incorporating accurate Lyα feedback and ionizing radiation to investigate the nature of the CR7 source. We find that a Population III star cluster with $10^5$ K blackbody emission ionizes its environment too efficiently to generate a significant velocity offset. However, a massive black hole with a nonthermal Compton-thick spectrum is able to reproduce the Lyα signatures as a result of higher photon trapping and longer potential lifetime. For both sources, Lyα radiation pressure turns out to be dynamically important.

preprint2016arXiv

Inferring the Distances of Fast Radio Bursts Through Associated 21-cm Absorption

The distances of Fast Radio Burst (FRB) sources are currently unknown. We show that the 21-cm absorption line of hydrogen can be used to infer the redshifts of FRB sources, and determine whether they are Galactic or extragalactic. We calculate a probability of $\sim 10\%$ for the host galaxy of an FRB to exhibit a 21-cm absorption feature of equivalent width $\gtrsim 10 ~\mathrm{km}~\mathrm{s}^{-1}$. Arecibo, along with several future radio observatories, should be capable of detecting such associated 21-cm absorption signals for strong bursts of $\gtrsim {\rm several}~\mathrm{Jy}$ peak flux densities.

preprint2016arXiv

Lessons from Mayan Astronomy

The Mayan culture collected exquisite astronomical data for over a millennium. However, it failed to come up with the breakthrough ideas of modern astronomy because the data was analyzed within a mythological culture of astrology that rested upon false but mathematically sophisticated theories about the Universe. Have we learned the necessary lessons to prevent our current scientific culture from resembling Mayan Astronomy? Clearly, data collection by itself is not a guarantee for good science as commonly assumed by funding agencies. A vibrant scientific culture should cultivate multiple approaches to analyzing existing data and to collecting new data.

preprint2016arXiv

Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization

The dynamical impact of Lyman-alpha (Lyα) radiation pressure on galaxy formation depends on the rate and duration of momentum transfer between Lyα photons and neutral hydrogen gas. Although photon trapping has the potential to multiply the effective force, ionizing radiation from stellar sources may relieve the Lyα pressure before appreciably affecting the kinematics of the host galaxy or efficiently coupling Lyα photons to the outflow. We present self-consistent Lyα radiation-hydrodynamics simulations of high-$z$ galaxy environments by coupling the Cosmic Lyα Transfer code (COLT) with spherically symmetric Lagrangian frame hydrodynamics. The accurate but computationally expensive Monte-Carlo radiative transfer calculations are feasible under the one-dimensional approximation. The initial starburst drives an expanding shell of gas from the centre and in certain cases Lyα feedback significantly enhances the shell velocity. Radiative feedback alone is capable of ejecting baryons into the intergalactic medium (IGM) for protogalaxies with a virial mass of $M_{\rm vir} \lesssim 10^8~{\rm M}_\odot$. We compare the Lyα signatures of Population III stars with $10^5$ K blackbody emission to that of direct collapse black holes with a nonthermal Compton-thick spectrum and find substantial differences if the Lyα spectra are shaped by gas pushed by Lyα radiation-driven winds. For both sources, the flux emerging from the galaxy is reprocessed by the IGM such that the observed Lyα luminosity is reduced significantly and the time-averaged velocity offset of the Lyα peak is shifted redward.

preprint2016arXiv

On the Habitability of Our Universe

Is life most likely to emerge at the present cosmic time near a star like the Sun? We consider the habitability of the Universe throughout cosmic history, and conservatively restrict our attention to the context of "life as we know it" and the standard cosmological model, LCDM. The habitable cosmic epoch started shortly after the first stars formed, about 30 Myr after the Big Bang, and will end about 10 Tyr from now, when all stars will die. We review the formation history of habitable planets and find that unless habitability around low mass stars is suppressed, life is most likely to exist near 0.1 solar mass stars ten trillion years from now. Spectroscopic searches for biosignatures in the atmospheres of transiting Earth-mass planets around low mass stars will determine whether present-day life is indeed premature or typical from a cosmic perspective.

preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

preprint2016arXiv

Precise Measurement of the Reionization Optical Depth from The Global 21-cm Signal Accounting for Cosmic Heating

As a result of our limited data on reionization, the total optical depth for electron scattering, $τ$, limits precision measurements of cosmological parameters from the Cosmic Microwave Background (CMB). It was recently shown that the predicted 21-cm signal of neutral hydrogen contains enough information to reconstruct $τ$ with sub-percent accuracy, assuming that the neutral gas was much hotter than the CMB throughout the entire epoch of reionization. Here we relax this assumption and use the global 21-cm signal alone to extract $τ$ for realistic X-ray heating scenarios. We test our model-independent approach using mock data for a wide range of ionization and heating histories and show that an accurate measurement of the reionization optical depth at a sub-percent level is possible in most of the considered scenarios even when heating is not saturated during the epoch of reionization, assuming that the foregrounds are mitigated. However, we find that in cases where heating sources had hard X-ray spectra and their luminosity was close to or lower than what is predicted based on low-redshift observations, the global 21-cm signal alone is not a good tracer of the reionization history.

preprint2016arXiv

Predicted Sizes of Pressure-Supported HI Clouds in the Outskirts of the Virgo Cluster

Using data from the ALFALFA Arecibo HI survey of galaxies and the Virgo cluster X-ray pressure profiles from XMM-Newton, we investigate the possibility that starless dark HI clumps, also known as "dark galaxies", are supported by external pressure in the surrounding intercluster medium. We find that the starless HI clump masses, velocity dispersions and positions allow these clumps to be in pressure equilibrium with the X-ray gas near the virial radius of the Virgo cluster. We predict the sizes of these clumps to range from 1-10 kpc, in agreement with the range of sizes found for spatially resolved HI starless clumps outside of Virgo. Based on the predicted HI surface density of the Virgo sources, as well as a sample of other similar resolved ALFALFA HI dark clumps with follow up optical/radio observations, we predict that most of the HI dark clumps are on the cusp of forming stars. These HI sources therefore mark the transition between starless HI clouds and dwarf galaxies with stars.

preprint2016arXiv

Prompt Radiation and Mass Outflows from the Stream-Stream Collisions of Tidal Disruption Events

Stream-stream collisions play an important role for the circularization of highly eccentric streams resulting from tidal disruption events (TDEs). We perform three dimensional radiation hydrodynamic simulations to show that stream collisions can contribute significant optical and ultraviolet light to the flares produced by TDEs, and can sometimes explain the majority of the observed emission. Our simulations focus on the region near the radiation pressure dominated shock produced by a collision and track how the kinetic energy of the stream is dissipated by the associated shock. When the mass flow rate of the stream $\dot{M}$ is a significant fraction of the Eddington accretion rate, $\gtrsim2\%$ of the initial kinetic energy is converted to radiation directly as a result of the collision. In this regime, the collision redistributes the specific kinetic energy into the downstream gas and more than $16\%$ of the mass can become unbound. The fraction of unbound gas decreases rapidly as $\dot{M}$ drops significantly below the Eddington limit, with no unbound gas being produced when $\dot{M}$ drops to $1\%$ of Eddington; we find however that the radiative efficiency increases slightly to $\lesssim 8\%$ in these low $\dot{M}$ cases. The effective radiation temperature and size of the photosphere is determined by the stream velocity and $\dot{M}$, which we find to be a few times $10^4$~K and $10^{14}$~cm in our calculations, comparable to the inferred values of some TDE candidates. The photosphere size is directly proportional to $\dot{M}$, which can explain the rapidly changing photosphere sizes seen in TDE candidates such as PS1-10jh.

preprint2016arXiv

Prospects for Characterizing the Atmosphere of Proxima Centauri b

The newly detected Earth-mass planet in the habitable zone of Proxima Centauri could potentially host life - if it has an atmosphere that supports surface liquid water. We show that thermal phase curve observations with the James Webb Space Telescope (JWST) from 5-12 microns can be used to test the existence of such an atmosphere. We predict the thermal variation for a bare rock versus a planet with 35% heat redistribution to the nightside and show that a JWST phase curve measurement can distinguish between these cases at $4σ$ confidence, assuming photon-limited precision. We also consider the case of an Earth-like atmosphere, and find that the ozone 9.8 micron band could be detected with longer integration times (a few months). We conclude that JWST observations have the potential to put the first constraints on the possibility of life around the nearest star to the Solar System.

preprint2016arXiv

Quasar-driven outflows account for the missing extragalactic gamma-ray background

The origin of the extragalactic $γ$-ray background permeating throughout the Universe remains a mystery forty years after its discovery. The extrapolated population of blazars can account for only half of the background radiation at the energy range of ~ 0.1-10 GeV. Here we show that quasar-driven outflows generate relativistic protons that produce the missing component of the extragalactic $γ$-ray background and naturally match its spectral fingerprint, with a generic break above ~ 1 GeV. The associated $γ$-ray sources are too faint to be detected individually, explaining why they had not been identified so far. However, future radio observations may image their shock fronts directly. Our best fit to the Fermi-LAT observations of extragalactic $γ$-ray background spectrum provides constraints on the outflow parameters that agree with observations of these outflows and theoretical predictions.

preprint2016arXiv

Relative Likelihood for Life as a Function of Cosmic Time

Is life most likely to emerge at the present cosmic time near a star like the Sun? We address this question by calculating the relative formation probability per unit time of habitable Earth-like planets within a fixed comoving volume of the Universe, dP(t)/dt, starting from the first stars and continuing to the distant cosmic future. We conservatively restrict our attention to the context of "life as we know it" and the standard cosmological model, LCDM. We find that unless habitability around low mass stars is suppressed, life is most likely to exist near 0.1 solar-mass stars ten trillion years from now. Spectroscopic searches for biosignatures in the atmospheres of transiting Earth-mass planets around low mass stars will determine whether present-day life is indeed premature or typical from a cosmic perspective.

preprint2016arXiv

Scaling Relations of Halo Cores for Self-Interacting Dark Matter

Using a simple analytic formalism, we demonstrate that significant dark matter self-interactions produce halo cores that obey scaling relations nearly independent of the underlying particle physics parameters such as the annihilation cross section and the mass of the dark matter particle. For dwarf galaxies, we predict that the core density $ρ_c$ and the core radius $r_c$ should obey $ρ_c r_c \approx 41 \,\text{M}_\odot \text{pc}^{-2}$ with a weak mass dependence $\sim M^{0.2}$. Remarkably, such a scaling relation has recently been empirically inferred. Scaling relations involving core mass, core radius, and core velocity dispersion are predicted and agree well with observational data. By calibrating against numerical simulations, we predict the scatter in these relations and find them to be in excellent agreement with existing data. Future observations can test our predictions for different halo masses and redshifts.

preprint2016arXiv

Spectral Distortion of the CMB by the Cumulative CO Emission from Galaxies throughout Cosmic History

We show that the cumulative CO emission from galaxies throughout cosmic history distorts the spectrum of the cosmic microwave background (CMB) at a level that is well above the detection limit of future instruments, such as the Primordial Inflation Explorer (PIXIE). The modeled CO signal has a prominent bump in the frequency interval 100-200 GHz, with a characteristic peak intensity of ~ 2$\times$10$^{-23}$ W m$^{-2}$ Hz$^{-1}$ sr$^{-1}$. Most of the CO foreground originates from modest redshifts, z ~ 2-5, and needs to be efficiently removed for more subtle distortions from the earlier universe to be detected.

preprint2016arXiv

Testing General Relativity with Accretion-Flow Imaging of Sgr A*

The Event Horizon Telescope is a global very-long baseline interferometer capable of probing potential deviations from the Kerr metric, which is believed to provide the unique description of astrophysical black holes. Here we report an updated constraint on the quadrupolar deviation of Sagittarius A* within the context of a radiatively inefficient accretion flow model in a quasi-Kerr background. We also simulate near-future constraints obtainable by the forthcoming eight-station array and show that in this model already a one-day observation can measure the spin magnitude to within $0.005$, the inclination to within $0.09^\circ$, the position angle to within $0.04^\circ$, and the quadrupolar deviation to within $0.005$ at $3σ$ confidence. Thus, we are entering an era of high-precision strong gravity measurements.

preprint2016arXiv

The Distortion of the Cosmic Microwave Background Spectrum Due to Intergalactic Dust

Infrared emission from intergalactic dust might compromise the ability of future experiments to detect subtle spectral distortions in the Cosmic Microwave Background (CMB) from the early Universe. We provide the first estimate of foreground contamination of the CMB signal due to diffuse dust emission in the intergalactic medium. We use models of the extragalactic background light to calculate the intensity of intergalactic dust emission and find that emission by intergalactic dust at redshifts z<0.5 exceeds the sensitivity of the planned Primordial Inflation Explorer (PIXIE) to CMB spectral distortions by 1-3 orders of magnitude. We place an upper limit of 0.23% on the contribution to the far-infrared background from intergalactic dust emission.

preprint2016arXiv

Zipf's Law from Scale-free Geometry

The spatial distribution of people exhibits clustering across a wide range of scales, from household ($\sim 10^{-2}$ km) to continental ($\sim 10^4$ km) scales. Empirical data indicates simple power-law scalings for the size distribution of cities (known as Zipf's law) and the population density fluctuations as a function of scale. Using techniques from random field theory and statistical physics, we show that these power laws are fundamentally a consequence of the scale-free spatial clustering of human populations and the fact that humans inhabit a two-dimensional surface. In this sense, the symmetries of scale invariance in two spatial dimensions are intimately connected to urban sociology. We test our theory by empirically measuring the power spectrum of population density fluctuations and show that the logarithmic slope $α= 2.04 \pm 0.09$, in excellent agreement with our theoretical prediction $α= 2$. The model enables the analytic computation of many new predictions by importing the mathematical formalism of random fields.

preprint2015arXiv

An Empirical Model for the Galaxy Luminosity and Star-Formation Rate Function at High Redshift

Using the most recent measurements of the ultraviolet (UV) luminosity functions (LFs) and dust estimates of early galaxies, we derive updated dust-corrected star-formation rate functions (SFRFs) at z~4-8, which we model to predict the evolution to higher redshifts, z>8. We employ abundance matching techniques to calibrate a relation between galaxy star formation rate (SFR) and host halo mass M$_{h}$ by mapping the shape of the observed SFRFs at z~4-8 to that of the halo mass function. The resulting scaling law remains roughly constant over this redshift range. We apply the average SFR-M$_{h}$ relation to reproduce the observed SFR functions at 4 <= z <= 8 and also derive the expected UV LFs at higher redshifts. At z~9 and z~10 these model LFs are in excellent agreement with current observed estimates. Our predicted number densities and UV LFs at z>10 indicate that JWST will be able to detect galaxies out to z~15 with an extensive treasury sized program. We also derive the redshift evolution of the star formation rate density and associated reionization history by galaxies. Models which integrate down to the current HUDF12/XDF detection limit (M$_{UV}$ ~ -17.7 mag) result in a SFRD that declines as (1+z)$^{-10.4\pm0.3}$ at high redshift and fail to reproduce the observed CMB electron scattering optical depth, $τ$ ~ 0.066, to within 1$σ$. On the other hand, we find that the inclusion of galaxies with SFRs well below the current detection limit (M$_{UV}$ < -5.7 mag) leads to a fully reionized universe by z ~ 6.5 and an optical depth of $τ$ ~ 0.054, consistent with the recently derived Planck value at the 1$σ$ level.

preprint2015arXiv

Density Jumps Near the Virial Radius of Galaxy Clusters

Recent simulations have indicated that the dark matter halos of galaxy clusters should feature steep density jumps near the virial radius. Since the member galaxies are expected to follow similar collisionless dynamics as the dark matter, the galaxy density profile should show such a feature as well. We examine the potential of current datasets to test this prediction by selecting cluster members for a sample of 56 low-redshift (0.1<z<0.3) galaxy clusters, constructing their projected number density profiles, and fitting them with two profiles, one with a steep density jump and one without. Additionally, we investigate the presence of a jump using a non-parametric spline approach. We find that some of these clusters show strong evidence for a model with a density jump. We discuss avenues for further analysis of the density jump with future datasets.

preprint2015arXiv

Detecting Stars at the Galactic Centre via Synchrotron Emission

Stars orbiting within 1$\arcsec$ of the supermassive black hole in the Galactic Centre, Sgr A*, are notoriously difficult to detect due to obscuration by gas and dust. We show that some stars orbiting this region may be detectable via synchrotron emission. In such instances, a bow shock forms around the star and accelerates the electrons. We calculate that around the 10 GHz band (radio) and at 10$^{14}$ Hz (infrared) the luminosity of a star orbiting the black hole is comparable to the luminosity of Sgr A*. The strength of the synchrotron emission depends on a number of factors including the star's orbital velocity. Thus, the ideal time to observe the synchrotron flux is when the star is at pericenter. The star S2 will be $\sim 0.015\arcsec$ from Sgr A* in 2018, and is an excellent target to test our predictions.

preprint2015arXiv

Distortion of the luminosity function of high-redshift galaxies by gravitational lensing

The observed properties of high redshift galaxies depend on the underlying foreground distribution of large scale structure, which distorts their intrinsic properties via gravitational lensing. We focus on the regime where the dominant contribution originates from a single lens and examine the statistics of gravitational lensing by a population of virialized and non-virialized structures using sub-mm galaxies at z ~ 2.6 and Lyman-break galaxies at redshifts z ~ 6 - 15 as the background sources. We quantify the effect of lensing on the luminosity function of the high redshift sources, focusing on the intermediate and small magnifications, mu < 2, which affect the majority of the background galaxies, and comparing to the case of strong lensing. We show that, depending on the intrinsic properties of the background galaxies, gravitational lensing can significantly affect the observed luminosity function even when no obvious strong lenses are present. Finally, we find that in the case of the Lyman-break galaxies it is important to account for the surface brightness profiles of both the foreground and the background galaxies when computing the lensing statistics, which introduces a selection criterion for the background galaxies that can actually be observed. Not taking this criterion into account leads to an overestimation of the number densities of very bright galaxies by nearly two orders of magnitude.

preprint2015arXiv

Evaporation and Accretion of Extrasolar Comets Following White Dwarf Kicks

Several lines of observational evidence suggest that white dwarfs receive small birth kicks due to anisotropic mass loss. If other stars possess extrasolar analogues to the Solar Oort cloud, the orbits of comets in such clouds will be scrambled by white dwarf natal kicks. Although most comets will be unbound, some will be placed on low angular momentum orbits vulnerable to sublimation or tidal disruption. The dusty debris from these comets will manifest itself as an IR excess temporarily visible around newborn white dwarfs; examples of such disks may already have been seen in the Helix Nebula, and around several other young white dwarfs. Future observations with the James Webb Space Telescope may distinguish this hypothesis from alternatives such as a dynamically excited Kuiper Belt analogue. Although competing hypotheses exist, the observation that $\gtrsim 15\%$ of young white dwarfs possess such disks, if interpreted as indeed being cometary in origin, provides indirect evidence that low mass gas giants (thought necessary to produce an Oort cloud) are common in the outer regions of extrasolar planetary systems. Hydrogen abundances in the atmospheres of older white dwarfs can, if sufficiently low, also be used to place constraints on the joint parameter space of natal kicks and exo-Oort cloud models.

preprint2015arXiv

Fast radio bursts: the observational case for a Galactic origin

There are by now ten published detections of fast radio bursts (FRBs), single bright GHz-band millisecond pulses of unknown origin. Proposed explanations cover a broad range from exotic processes at cosmological distances to atmospheric and terrestrial sources. Loeb et al. have previously suggested that FRB sources could be nearby flare stars, and pointed out the presence of a W-UMa-type contact binary within the beam of one out of three FRB fields that they examined. Using time-domain optical photometry and spectroscopy, we now find possible flare stars in additional FRB fields, with one to three such cases among eight FRB fields studied. We evaluate the chance probabilities of these possible associations to be in the range 0.1% to 9%, depending on the input assumptions. Further, we re-analyze the probability that two FRBs recently discovered 3 years apart within the same radio beam are unrelated. Contrary to other claims, we conclude with 99% confidence that the two events are from the same repeating source. The different dispersion measures between the two bursts then rule out a cosmological origin for the dispersion measure, but are consistent with the flare-star scenario with a varying plasma blanket between bursts. Finally, we review some theoretical objections that have been raised against a local flare-star FRB origin, and show that they are incorrect.

preprint2015arXiv

Implications of the Eccentric Kozai-Lidov Mechanism for Stars Surrounding Supermassive Black Hole Binaries

An enhanced rate of stellar tidal disruption events (TDEs) may be an important characteristic of supermassive black hole (SMBH) binaries at close separations. Here we study the evolution of the distribution of stars around a SMBH binary due to the eccentric Kozai-Lidov (EKL) mechanism, including octupole effects and apsidal precession caused by the stellar mass distribution and general relativity. We identify a region around one of the SMBHs in the binary where the EKL mechanism drives stars to high eccentricities, which ultimately causes the stars to either scatter off the second SMBH or get disrupted. For SMBH masses 10^7 Msun and 10^8 Msun, the TDE rate can reach 10^{-2} yr and deplete a region of the stellar cusp around the secondary SMBH in ~0.5 Myr. As a result, the final geometry of the stellar distribution between 0.01 and 0.1 pc around the secondary SMBH is a torus. These effects may be even more prominent in nuclear stellar clusters hosting a supermassive and an intermediate mass black hole.

preprint2015arXiv

Limits on Intergalactic Dust During Reionization

In this Letter, we constrain the dust-to-gas ratio in the intergalactic medium (IGM) at high redshifts. We employ models for dust in the local Universe to contrain the dust-to-gas ratio during the epoch of reionization at redshifts z ~ 6-10. The observed level of reddening of high redshift galaxies implies that the IGM was enriched to an intergalactic dust-to-gas ratio of less than 3% of the Milky Way value by a redshift of z=10.

preprint2015arXiv

Mapping the Dynamics of Cold Gas around Sgr A* through 21 cm Absorption

The presence of a circumnuclear stellar disk around Sgr A* and megamaser systems near other black holes indicates that dense neutral disks can be found in galactic nuclei. We show that depending on their inclination angle, optical depth, and spin temperature, these disks could be observed spectroscopically through 21 cm absorption. Related spectroscopic observations of Sgr A* can determine its HI disk parameters and the possible presence of gaps in the disk. Clumps of dense gas similar to the G2 could could also be detected in 21 cm absorption against Sgr A* radio emission.

preprint2015arXiv

Measuring the Direction and Angular Velocity of a Black Hole Accretion Disk via Lagged Interferometric Covariance

We show that interferometry can be applied to study irregular, rapidly rotating structures, as are expected in the turbulent accretion flow near a black hole. Specifically, we analyze the lagged covariance between interferometric baselines of similar lengths but slightly different orientations. For a flow viewed close to face-on, we demonstrate that the peak in the lagged covariance indicates the direction and angular velocity of the emission pattern from the flow. Even for moderately inclined flows, the covariance robustly estimates the flow direction, although the estimated angular velocity can be significantly biased. Importantly, measuring the direction of the flow as clockwise or counterclockwise on the sky breaks a degeneracy in accretion disk inclinations when analyzing time-averaged images alone. We explore the potential efficacy of our technique using three-dimensional, general relativistic magnetohydrodynamic (GRMHD) simulations, and we highlight several baseline pairs for the Event Horizon Telescope (EHT) that are well-suited to this application. These results indicate that the EHT may be capable of estimating the direction and angular velocity of the emitting material near Sagittarius A*, and they suggest that a rotating flow may even be utilized to improve imaging capabilities.

preprint2015arXiv

Observational Cosmology With Semi-Relativistic Stars

Galaxy mergers lead to the formation of massive black hole binaries which can accelerate background stars close to the speed of light. We estimate the comoving density of ejected stars with a peculiar velocity in excess of $0.1c$ or $0.5c$ to be $\sim 10^{10}$ and $10^5$ Gpc$^{-3}$ respectively, in the present-day Universe. Semi-relativistic giant stars will be detectable with forthcoming telescopes out to a distance of a few Mpc, where their proper motion, radial velocity, and age, can be spectroscopically measured. In difference from traditional cosmological messengers, such as photons, neutrinos, or cosmic-rays, these stars shine and so their trajectories need not be directed at the observer for them to be detected. Tracing the stars to their parent galaxies as a function of speed and age will provide a novel test of the equivalence principle and the standard cosmological parameters. Semi-relativistic stars could also flag black hole binaries as gravitational wave sources for the future eLISA observatory.

preprint2015arXiv

Optical inverse-Compton emission from clusters of galaxies

Shocks around clusters of galaxies accelerate electrons which upscatter the Cosmic Microwave Background photons to higher-energies. We use an analytical model to calculate this inverse Compton (IC) emission, taking into account the effects of additional energy losses via synchrotron and Coulomb scattering. We find that the surface brightness of the optical IC emission increases with redshift and halo mass. The IC emission surface brightness, 32--34~mag~arcsec$^{-2}$, for massive clusters is potentially detectable by the newly developed Dragonfly Telephoto Array.

preprint2015arXiv

Predicting the intensity mapping signal for multi-$J$ CO lines

We present a novel approach to estimating the intensity mapping signal of any CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic temperature, volume density, velocity gradient, molecular abundance, and column density. These parameters, which drive the physics of CO transitions and ultimately dictate the shape and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, mainly the star formation rate (SFR) and the SFR surface density. By further employing an empirically derived SFR-M relation for high redshift galaxies, we can express the LVG parameters, and thus the specific intensity of any CO rotational transition, as functions of the host halo mass M and redshift z. Integrating over the range of halo masses expected to host CO-luminous galaxies, we predict a mean CO(1-0) brightness temperature ranging from ~0.6 μK at z= 6 to ~0.03 μ at z= 10 with brightness temperature fluctuations of $Δ_{CO}^2$ ~ 0.1 and 0.005 μK respectively, at k = 0.1 Mpc$^{-1}$. In this model, the CO emission signal remains strong for higher rotational levels at z = 6, with <T$_{CO}$>~ 0.3 and 0.05 μK for the CO J = 6->5 and CO J = 10->9 transitions respectively. Including the effects of CO photodissociation in these molecular clouds, especially at low metallicities, results in the overall reduction in the amplitude of the CO signal, with the low- and high-J lines weakening by 2-20% and 10-45%, respectively, over the redshift range 4 < z < 10.

preprint2015arXiv

Probing the gaseous halo of galaxies through non-thermal emission from AGN-driven outflows

Feedback from outflows driven by active galactic nuclei (AGN) can affect the distribution and properties of the gaseous halos of galaxies. We study the hydrodynamics and non-thermal emission from the forward outflow shock produced by an AGN-driven outflow. We consider a few possible profiles for the halo gas density, self-consistently constrained by the halo mass, redshift and the disk baryonic concentration of the galaxy. We show that the outflow velocity levels off at $\sim 10^3\,\rm km\, s^{-1}$ within the scale of the galaxy disk. Typically, the outflow can reach the virial radius around the time when the AGN shuts off. We show that the outflows are energy-driven, consistently with observations and recent theoretical findings. The outflow shock lights up the halos of massive galaxies across a broad wavelength range. For Milky Way (MW) mass halos, radio observations by The Jansky Very Large Array (JVLA) and The Square Kilometer Array (SKA) and infrared/optical observations by The James Webb Space Telescope (JWST) and Hubble Space Telescope (HST) can detect the emission signal of angular size $\sim 8"$ from galaxies out to redshift $z\sim5$. Millimeter observations by The Atacama Large Millimeter/submillimeter Array (ALMA) are sensitive to non-thermal emission of angular size $\sim 18"$ from galaxies at redshift $z\lesssim1$, while X-ray observations by Chandra, XMM-Newton and The Advanced Telescope for High Energy Astrophysics (ATHENA) is limited to local galaxies ($z\lesssim 0.1$) with an emission angular size of $\sim2'$. Overall, the extended non-thermal emission provides a new way of probing the gaseous halos of galaxies at high redshifts.

preprint2015arXiv

Probing the Spacetime Around Supermassive Black Holes with Ejected Plasma Blobs

Millimeter-wavelength VLBI observations of the supermassive black holes in Sgr A* and M87 by the Event Horizon Telescope could potentially trace the dynamics of ejected plasma blobs in real time. We demonstrate that the trajectory and tidal stretching of these blobs can be used to test general relativity and set new constraints on the mass and spin of these black holes.

preprint2015arXiv

Quantum coherent oscillations in the early universe

Cosmic inflation is commonly assumed to be driven by quantum fields. Quantum mechanics predicts phenomena such as quantum fluctuations and tunneling of the field. Here we show an example of a quantum interference effect which goes beyond the semi-classical treatment and which may be of relevance in the early universe. We study the quantum coherent dynamics for a tilted, periodic potential, which results in genuine quantum oscillations of the inflaton field, analogous to Bloch oscillations in condensed matter and atomic systems. Our results show that quantum interference phenomena may be of relevance in cosmology.

preprint2015arXiv

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

preprint2015arXiv

SETI via Leakage from Light Sails in Exoplanetary Systems

The primary challenge of rocket propulsion is the burden of needing to accelerate the spacecraft's own fuel, resulting in only a logarithmic gain in maximum speed as propellant is added to the spacecraft. Light sails offer an attractive alternative in which fuel is not carried by the spacecraft, with acceleration being provided by an external source of light. By artificially illuminating the spacecraft with beamed radiation, speeds are only limited by the area of the sail, heat resistance of its material, and power use of the accelerating apparatus. In this paper, we show that leakage from a light sail propulsion apparatus in operation around a solar system analogue would be detectable. To demonstrate this, we model the launch and arrival of a microwave beam-driven light sail constructed for transit between planets in orbit around a single star, and find an optimal beam frequency on the order of tens of GHz. Leakage from these beams yields transients with flux densities of Jy and durations of tens of seconds at 100 pc. Because most travel within a planetary system would be conducted between the habitable worlds within that system, multiply-transiting exoplanetary systems offer the greatest chance of detection, especially when the planets are in projected conjunction as viewed from Earth. If interplanetary travel via beam-driven light sails is commonly employed in our galaxy, this activity could be revealed by radio follow-up of nearby transiting exoplanetary systems. The expected signal properties define a new strategy in the search for extraterrestrial intelligence (SETI).

preprint2015arXiv

Shapiro Delays at the Quadrupole Order for Tests of the No-Hair Theorem Using Pulsars around Spinning Black Holes

One avenue for testing the no-hair theorem is obtained through timing a pulsar orbiting close to a black hole and fitting for quadrupolar effects on the time-of-arrival of pulses. If deviations from the Kerr quadrupole are measured, then the no-hair theorem is invalidated. To this end, we derive an expression for the light travel time delay for a pulsar orbiting in a black-hole spacetime described by the Butterworth-Ipser metric, which has an arbitrary spin and quadrupole moment. We consider terms up to the quadrupole order in the black-hole metric and derive the time-delay expression in a closed analytic form. This allows for fast computations that are useful in fitting time-of-arrival observations of pulsars orbiting close to astrophysical black holes.

preprint2015arXiv

Statistical Signatures of Panspermia in Exoplanet Surveys

A fundamental astrobiological question is whether life can be transported between extrasolar systems. We propose a new strategy to answer this question based on the principle that life which arose via spreading will exhibit more clustering than life which arose spontaneously. We develop simple statistical models of panspermia to illustrate observable consequences of these excess correlations. Future searches for biosignatures in the atmospheres of exoplanets could test these predictions: a smoking gun signature of panspermia would be the detection of large regions in the Milky Way where life saturates its environment interspersed with voids where life is very uncommon. In a favorable scenario, detection of as few as $\sim 25$ biologically active exoplanets could yield a $5σ$ detection of panspermia. Detectability of position-space correlations is possible unless the timescale for life to become observable once seeded is longer than the timescale for stars to redistribute in the Milky Way.

preprint2015arXiv

The Fastest Unbound Stars in the Universe

The discovery of hypervelocity stars (HVS) leaving our galaxy with speeds of nearly $10^{3}$ km s$^{-1}$ has provided strong evidence towards the existence of a massive compact object at the galaxy's center. HVS ejected via the disruption of stellar binaries can occasionally yield a star with $v_{\infty} \lesssim 10^4$ km s$^{-1}$, here we show that this mechanism can be extended to massive black hole (MBH) mergers, where the secondary star is replaced by a MBH with mass $M_2 \gtrsim 10^5 M_{\odot}$. We find that stars that are originally bound to the secondary MBH are frequently ejected with $v_{\infty} > 10^4$ km s$^{-1}$, and occasionally with velocities $\sim 10^5$ km s$^{-1}$ (one third the speed of light), for this reason we refer to stars ejected from these systems as "semi-relativistic" hypervelocity stars (SHS). Bound to no galaxy, the velocities of these stars are so great that they can cross a significant fraction of the observable universe in the time since their ejection (several Gpc). We demonstrate that if a significant fraction of MBH mergers undergo a phase in which their orbital eccentricity is $\gtrsim 0.5$ and their periapse distance is tens of the primary's Schwarzschild radius, the space density of fast-moving ($v_{\infty} > 10^{4}$ km s$^{-1}$) SHS may be as large as $10^{3}$ Mpc$^{-3}$. Hundreds of the SHS will be giant stars that could be detected by future all-sky infrared surveys such as WFIRST or Euclid and proper motion surveys such as LSST, with spectroscopic follow-up being possible with JWST.

preprint2015arXiv

Using Atomic Clocks to Detect Gravitational Waves

Atomic clocks have recently reached a fractional timing precision of $<10^{-18}$. We point out that an array of atomic clocks, distributed along the Earth's orbit around the Sun, will have the sensitivity needed to detect the time dilation effect of mHz gravitational waves (GWs), such as those emitted by supermassive black hole binaries at cosmological distances. Simultaneous measurement of clock-rates at different phases of a passing GW provides an attractive alternative to the interferometric detection of temporal variations in distance between test masses separated by less than a GW wavelength, currently envisioned for the eLISA mission.

preprint2014arXiv

A Simple Physical Model for the Gas Distribution in Galaxy Clusters

The dominant baryonic component of galaxy clusters is hot gas whose distribution is commonly probed through X-ray emission arising from thermal bremsstrahlung. The density profile thus obtained has been traditionally modeled with a beta-profile, a simple function with only three parameters. However, this model is known to be insufficient for characterizing the range of cluster gas distributions, and attempts to rectify this shortcoming typically introduce additional parameters to increase the fitting flexibility. We use cosmological and physical considerations to obtain a family of profiles for the gas with fewer parameters than the beta-model but which better accounts for observed gas profiles over wide radial intervals.

preprint2014arXiv

A size-duration trend for gamma-ray burst progenitors

Gamma-ray bursts (GRBs) show a bimodal distribution of durations, separated at a duration of ~2 s. Observations have confirmed the association of long GRBs with the collapse of massive stars. The origin of short GRBs is still being explored. We examine constraints on the emission region size in short and long GRBs detected by Fermi/GBM. We find that the emission region size during the prompt emission, R, and the burst duration, T$_{90}$, are consistent with the relation R ~ c x T$_{90}$, for both long and short GRBs. We find the characteristic size for the prompt emission region to be ~2 x 10$^{10}$ cm, and ~4 x 10$^{11}$ cm for short and long GRBs, respectively.

preprint2014arXiv

Chaos in the Test Particle Eccentric Kozai-Lidov Mechanism

The Kozai-Lidov mechanism can be applied to a vast variety of astrophysical systems involving hierarchical three-body systems. Here, we study the Kozai-Lidov mechanism systematically in the test particle limit at the octupole level of approximation. We investigate the chaotic and quasiperiodic orbital evolution by studying surfaces of section and the Lyapunov exponents. We find that the resonances introduced by the octupole level of approximation cause orbits to flip from prograde to retrograde and back as well as cause significant eccentricity excitation, and the chaotic behaviors occur when the mutual inclination between the inner and the outer binary is high. We characterize the parameter space that allows large amplitude oscillations in eccentricity and inclination.

preprint2014arXiv

Constraining Sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. II. The Population with Kinematically Offset Broad Balmer Emission Lines

A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]

preprint2014arXiv

Detecting floating black holes as they traverse the gas disk of the Milky Way

A population of intermediate-mass black holes (BHs) is predicted to be freely floating in the Milky Way (MW) halo, due to gravitational wave recoil, ejection from triple BH systems, or tidal stripping in the dwarf galaxies that merged to make the MW. As these BHs traverse the gaseous MW disk, a bow shock forms, producing detectable radio and mm/sub-mm synchrotron emission from accelerated electrons. We calculate the synchrotron flux to be $\sim \rm 0.01-10\, mJy$ at GHz frequency, detectable by Jansky Very Large Array, and $\sim 10-100\,μ\rm Jy$ at $\sim10^{10}-10^{12} \,\rm Hz$ frequencies, detectable by Atacama Large Millimeter/sub-millimter Array. The discovery of the floating BH population will provide insights on the formation and merger history of the MW as well as on the evolution of massive BHs in the early Universe.

preprint2014arXiv

Detecting industrial pollution in the atmospheres of earth-like exoplanets

Detecting biosignatures, such as molecular oxygen in combination with a reducing gas, in the atmospheres of transiting exoplanets has been a major focus in the search for alien life. We point out that in addition to these generic indicators, anthropogenic pollution could be used as a novel biosignature for intelligent life. To this end, we identify pollutants in the Earth's atmosphere that have significant absorption features in the spectral range covered by the James Webb Space Telescope (JWST). We focus on tetrafluoromethane CF4 and trichlorofluoromethane (CCl3F), which are the easiest to detect chlorofluorocarbons (CFCs) produced by anthropogenic activity. We estimate that ~1.2 days (~1.7 days) of total integration time will be sufficient to detect or constrain the concentration of CCl3F (CF4) to ~10 times current terrestrial level.

preprint2014arXiv

Eccentricity growth and orbit flip in coplanar hierarchical three body systems

The dynamical evolution of a hierarchical three body system is well characterized by the eccentric Kozai-Lidov mechanism, where the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40-140 degree) can produce long timescale modulations that drives the eccentricity to extremely large value and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ~180 degree, rolling over its major axis. The ~180 degree flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exo-planetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems.

preprint2014arXiv

Effects of Intermediate Mass Black Holes on Nuclear Star Clusters

Nuclear star clusters (NSCs) are dense stellar clusters observed in galactic nuclei, typically hosting a central massive black hole. Here we study the possible formation and evolution of NSCs through the inspiral of multiple star clusters hosting intermediate mass black holes (IMBHs). Using an N-body code we examine the dynamics of the IMBHs and their effects on the NSC. We find that IMBHs inspiral to the core of the newly formed NSC and segregate there. Although the IMBHs scatter each other and the stars, none of them is ejected from the NSC. The IMBHs are excited to high eccentricities and their radial density profile develops a steep power-law cusp. The stars also develop a power-law cusp (instead of the central core that forms in their absence), but with a shallower slope. The relaxation rate of the NSC is accelerated due to the presence of IMBHs, which act as massive-perturbers. This in turn fills the loss-cone and boosts the tidal disruption rate of stars both by the MBH and the IMBHs to a value excluded by rate estimates based on current observations. Rate estimates of tidal disruptions can therefore provide a cumulative constraint on the existence of IMBHs in NSCs.

preprint2014arXiv

Event-Horizon-Telescope Evidence for Alignment of the Black Hole in the Center of the Milky Way with the Inner Stellar Disk

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

preprint2014arXiv

Finding rocky asteroids around white dwarfs by their periodic thermal emission

Since white dwarfs are small, the contrast between the thermal emission of an orbiting object and a white dwarf is dramatically enhanced compared to a main sequence host. Furthermore, rocky objects much smaller than the moon have no atmospheres and are tidally locked to the white dwarf. We show that this leads to temperature contrasts between their day and night side of order unity that should lead to temporal variations in infrared flux over an orbital period of $\sim 0.2$ to $\sim 2$ days. Ground based telescopes could detect objects with a mass as small as $1\%$ of the lunar mass $M_L$ around Sirius B with a few hours of exposure. The James Webb Space Telescope (JWST) may be able to detect objects as small as $10^{-3} M_L$ around most nearby white dwarfs. The tightest constraints will typically be placed on 12,000 K white dwarfs, whose Roche zone coincides with the dust sublimation zone. Constraining the abundance of minor planets around white dwarfs as a function of their surface temperatures (and therefore age) provides a novel probe for the physics of planetary formation.

preprint2014arXiv

Inflationary schism after Planck2013

Classic inflation, the theory described in textbooks, is based on the idea that, beginning from typical initial conditions and assuming a simple inflaton potential with a minimum of fine-tuning, inflation can create exponentially large volumes of space that are generically homogeneous, isotropic and flat, with nearly scale-invariant spectra of density and gravitational wave fluctuations that are adiabatic, Gaussian and have generic predictable properties. In a recent paper, we showed that, in addition to having certain conceptual problems known for decades, classic inflation is for the first time also disfavored by data, specifically the most recent data from WMAP, ACT and Planck2013. Guth, Kaiser and Nomura and Linde have each recently published critiques of our paper, but, as made clear here, we all agree about one thing: the problematic state of classic inflation. Instead, they describe an alternative inflationary paradigm that revises the assumptions and goals of inflation, and perhaps of science generally.

preprint2014arXiv

Is the Universe Simpler than LCDM?

In the standard cosmological model, the Universe consists mainly of two invisible substances: vacuum energy with constant mass-density rho_v=Λ/(8pi G) (where Lambda is a `cosmological constant' originally proposed by Einstein and G is Newton's gravitational constant) and cold dark matter (CDM) with mass density that is currently rho_{DM,0}\sim 0.3 rho_v. This `LCDM' model has the virtue of simplicity, enabling straightforward calculation of the formation and evolution of cosmic structure against the backdrop of cosmic expansion. Here we review apparent discrepancies with observations on small galactic scales, which LCDM must attribute to complexity in the baryon physics of galaxy formation. Yet galaxies exhibit structural scaling relations that evoke simplicity, presenting a clear challenge for formation models. In particular, tracers of gravitational potentials dominated by dark matter show a correlation between orbital size, R, and velocity, V, that can be expressed most simply as a characteristic acceleration, a_{DM}\sim 1 km^2 s^{-2} pc^{-1} \approx 3 x 10^{-9} cm s^{-2} \approx 0.2c\sqrt{G rho_v}, perhaps motivating efforts to find a link between localized and global manifestations of the Universe's dark components.

preprint2014arXiv

New insights from deep VLA data on the potentially recoiling black hole CID-42 in the COSMOS field

We present deep 3 GHz VLA observations of the potentially recoiling black hole CID-42 in the COSMOS field. This galaxy shows two optical nuclei in the HST/ACS image and a large velocity offset of ~ 1300 km/s between the broad and narrow H beta emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data has a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multi-scale multi-frequency synthesis and spectral windows stacking. The final resolutions and sensitivities of these maps are 0.7" with rms = 4.6 muJy/beam and 0.9" with rms = 4.8 muJy/beam respectively. With a 7 sigma detection we find that the entire observed 3 GHz radio emission can be associated with the South-Eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, VLBA and GMRT data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum AGN in the South-Eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole. Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the North-Western component.

preprint2014arXiv

On the Benefits of Promoting Diversity of Ideas

A common flaw of astronomers is to believe that they know the truth even when data is scarce. I provide ten examples where this approach led to wrong strategic decisions in research plans, causing unnecessary delays in finding the truth. In order to make the discovery process more efficient, I advocate affirmative action for diversity of ideas by telescope-time allocation committees and funding agencies.

preprint2014arXiv

Possible Origin of the G2 Cloud from the Tidal Disruption of a Known Giant Star by Sgr A*

The discovery of the gas cloud G2 on a near-radial orbit about Sgr A* has prompted much speculation on its origin. In this Letter, we propose that G2 formed out of the debris stream produced by the removal of mass from the outer envelope of a nearby giant star. We perform hydrodynamical simulations of the returning tidal debris stream with cooling, and find that the stream condenses into clumps that fall periodically onto Sgr A*. We propose that one of these clumps is the observed G2 cloud, with the rest of the stream being detectable at lower Br-$γ$ emissivity along a trajectory that would trace from G2 to the star that was partially disrupted. By simultaneously fitting the orbits of S2, G2, and $\sim$ 2,000 candidate stars, and by fixing the orbital plane of each candidate star to G2 (as is expected for a tidal disruption), we find that several stars have orbits that are compatible with the notion that one of them was tidally disrupted to produce G2. If one of these stars were indeed disrupted, it last encountered Sgr A* hundreds of years ago, and has likely encountered Sgr A* repeatedly. However, while these stars are compatible with the giant disruption scenario given their measured positions and proper motions, their radial velocities are currently unknown. If one of these stars' radial velocity is measured to be compatible with a disruptive orbit, it would strongly suggest its disruption produced G2.

preprint2014arXiv

Pulsar Timing Constraints on Cumulative and Individual Mass of Stars in the Galactic Center

We consider the time derivatives of the period $P$ of pulsars at the Galactic Center due to variations in their orbital Doppler shifts. We show that in conjunction with a measurement of a pulsar's proper motion and its projected separation from the supermassive black hole, Sgr A*,measuring two of the three derivatives $\dot{P}$, $\ddot{P}$, or $\dddot{P}$ sets a constraint that allows for the recovery of the complete six phase space coordinates of the pulsar's orbit, as well as the enclosed mass within the orbit. Thus, one can use multiple pulsars at different distances from Sgr A* to determine the radial mass distribution of stars and stellar remnants at the Galactic center. Furthermore, we consider the effect of passing stars on the pulsar's period derivatives and show how it can be exploited to measure the characteristic stellar mass in the Galactic Center.

preprint2014arXiv

Pulsations in Short GRBs from Black Hole-Neutron Star Mergers

The precise origin of short gamma ray bursts (SGRBs) remains an important open question in relativistic astrophysics. Increasingly, observational evidence suggests the merger of a binary compact object system as the source for most SGRBs, but it is currently unclear how to distinguish observationally between a binary neutron star progenitor and a black hole-neutron star progenitor. We suggest the quasi-periodic signal of jet precession as an observational signature of SGRBs originating in mixed binary systems, and quantify both the fraction of mixed binaries capable of producing SGRBs, and the distributions of precession amplitudes and periods. The difficulty inherent in disrupting a neutron star outside the horizon of a stellar-mass black hole biases the jet precession signal towards low amplitude and high frequency. Precession periods of ~ 0.01-0.1 s and disk-BH spin misalignments ~10 degrees are generally expected, although sufficiently high viscosity may prevent the accumulation of multiple precession periods during the SGRB. The precessing jet will naturally cover a larger solid angle in the sky than would standard SGRB jets, enhancing observability for both prompt emission and optical afterglows.

preprint2014arXiv

Reionization on Large Scales IV: Predictions for the 21 cm signal incorporating the light cone effect

We present predictions for the 21 cm brightness temperature power spectrum during the Epoch of Reionization (EoR). We discuss the implications of the "light cone" effect, which incorporates evolution of the neutral hydrogen fraction and 21 cm brightness temperature along the line of sight. Using a novel method calibrated against radiation-hydrodynamic simulations, we model the neutral hydrogen density field and 21 cm signal in large volumes ($L = 2$ Gpc/$h$). The inclusion of the light cone effect leads to a relative decrease of about 50% in the 21 cm power spectrum on all scales. We also find that the effect is more prominent at the midpoint of reionization and later. The light cone effect also can introduce an anisotropy along the line of sight. By decomposing the 3D power spectrum into components perpendicular to and along the line of sight, we find that in our fiducial reionization model, there is no significant anisotropy. However, parallel modes can contribute up to 40% more power for shorter reionization scenarios. The scales on which the light cone effect is relevant are comparable to scales where one measures the baryon acoustic oscillation. We argue that due to its large comoving scale and introduction of anisotropy, the light cone effect is important when considering redshift space distortions and future application to the Alcock-Paczynski test for the determination of cosmological parameters.

preprint2014arXiv

Relativistic redshifts in quasar broad lines

The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital Sky Survey DR7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ~30-45 degrees, consistent with simple AGN unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.

preprint2014arXiv

Signatures of the M31-M32 Galactic Collision

The unusual morphologies of the Andromeda spiral galaxy (M31) and its dwarf companion M32 have been characterized observationally in great detail. The two galaxies' apparent proximity suggests that Andromeda's prominent star-forming ring as well as M32's compact elliptical structure may result from a recent collision. Here we present the first self-consistent model of the M31-M32 interaction that simultaneously reproduces observed positions, velocities, and morphologies for both galaxies. Andromeda's spiral structure is resolved in unprecedented detail, showing that a rare head-on orbit is not necessary to match Andromeda's ring-like morphology. The passage of M32 through Andromeda's disk perturbs the disk velocity structure. We find tidal stripping of M32's stars to be inefficient during the interaction, suggesting that some cEs are intrinsically compact. Additionally, the orbital solution implies that M32 is currently closer to the Milky Way than models have typically assumed, a prediction that may be testable with upcoming observations.

preprint2014arXiv

The Habitable Epoch of the Early Universe

In the redshift range 100<(1+z)<137, the cosmic microwave background (CMB) had a temperature of 273-373K (0-100 degrees Celsius), allowing early rocky planets (if any existed) to have liquid water chemistry on their surface and be habitable, irrespective of their distance from a star. In the standard LCDM cosmology, the first star-forming halos within our Hubble volume started collapsing at these redshifts, allowing the chemistry of life to possibly begin when the Universe was merely 10-17 million years old. The possibility of life starting when the average matter density was a million times bigger than it is today argues against the anthropic explanation for the low value of the cosmological constant.

preprint2014arXiv

The kinetic Sunyaev-Zel'dovich effect from the diffuse gas in the Local Group

Since the Local Group (LG) of galaxies moves with a bulk velocity with respect to the cosmic microwave background radiation (CMB), free electrons in its gaseous halo should imprint large-scale non-primordial temperature shifts in the CMB via the kinetic Sunyaev-Zel'dovich (kSZ) effect. By modeling the distribution of gas in the LG halo and using its inferred velocity with respect to the CMB, we calculate the resulting kSZ signal from the diffuse LG medium. We find that it is dominated by a hot spot $\sim 10^\circ$ in size in the direction of M31, where the optical depth of free electrons is the greatest. By performing a correlation analysis, we find no statistical evidence that the kSZ signal from model of the LG halo is embedded in the CMB temperature map measured by the Planck satellite. We constrain the amount of mass in the LG halo by limiting the kSZ temperature shift around the hot spot to be smaller than the observed temperature shift in the Planck map. We find the tightest constraints for models where the halo mass is highly concentrated, with the mass limited to roughly $2.5-5\times 10^{12}$M$_\odot$, but note that halos with such high concentrations are rare.

preprint2014arXiv

The Mutual Interaction Between Population III Stars and Self-Annihilating Dark Matter

We use cosmological simulations of high-redshift minihalos to investigate the effect of dark matter annihilation (DMA) on the collapse of primordial gas. We numerically investigate the evolution of the gas as it assembles in a Population III stellar disk. We find that when DMA effects are neglected, the disk undergoes multiple fragmentation events beginning at ~ 500 yr after the appearance of the first protostar. On the other hand, DMA heating and ionization of the gas speeds the initial collapse of gas to protostellar densities and also affects the stability of the developing disk against fragmentation, depending on the DM distribution. We compare the evolution when we model the DM density with an analytical DM profile which remains centrally peaked, and when we simulate the DM profile using N-body particles (the 'live' DM halo). When utilizing the analytical DM profile, DMA suppresses disk fragmentation for ~ 3500 yr after the first protostar forms, in agreement with earlier work. However, when using a 'live' DM halo, the central DM density peak is gradually flattened due to the mutual interaction between the DM and the rotating gaseous disk, reducing the effects of DMA on the gas, and enabling secondary protostars of mass ~ 1 M_sol to be formed within ~ 900 yr. These simulations demonstrate that DMA is ineffective in suppressing gas collapse and subsequent fragmentation, rendering the formation of long-lived dark stars unlikely. However, DMA effects may still be significant in the early collapse and disk formation phase of primordial gas evolution.

preprint2014arXiv

Tidal capture of a primordial black hole by a neutron star: implications for constraints on dark matter

In a close encounter with a neutron star, a primordial black hole can get gravitationally captured by depositing a considerable amount of energy into nonradial stellar modes of very high angular number $l$. If the neutron-star equation of state is sufficiently stiff, we show that the total energy loss in the point-particle approximation is formally divergent. Various mechanisms -including viscosity, finite-size effects and the elasticity of the crust- can damp high-$l$ modes and regularize the total energy loss. Within a short time, the black hole is trapped inside the star and disrupts it by rapid accretion. Estimating these effects, we predict that the existence of old neutron stars in regions where the dark-matter density rho_{DM}>10^2 sigma/(km/s) GeV/cm^3 (where sigma is the dark-matter velocity dispersion) limits the abundance of primordial black holes in the mass range 10^{17} g < m_{PBH} < 10^{24} g, which was previously unconstrained. In combination with existing limits, our results suggest that primordial black holes cannot be the dominant dark matter constituent.

preprint2013arXiv

A new rare type of supernovae: hypervelocity stellar collisions at galactic centers

When a binary star system is tidally disrupted by a supermassive black hole at a galactic nucleus, one star is ejected at a high speed while the other remains in a tightly bound orbit around the black hole. The cluster of tightly bound stars builds over time, eventually creating a steady state in which the rate of collisions between these stars is similar to the rate of capturing new stars. A large fraction of the collisions occur near the periapsis of the orbits around the black hole, where the kinetic energies are sufficient to generate an explosive disruption of the two stars involved. The typical flare brightens for several days, with a peak luminosity that is comparable to the lower-luminosity end of known supernovae. The explosion lightcurve is followed by a longer flare due to accretion of ejected matter onto the black hole. Dedicated searches in the near universe could observe several such "collisional-supernovae" per year.

preprint2013arXiv

A Predicted New Population of UV-faint Galaxies at z>4

We show that a bursty model of high redshift star formation explains several puzzling observations of the high redshift galaxy population. We begin by pointing out that the observed specific star formation rate requires a duty-cycle of ~10%, which is much lower than found in many hydro-dynamical simulations. This value follows directly from the fact that the observed star formation rate in galaxies integrated over a Hubble time exceeds the observed stellar mass by an order of magnitude. We use the large observed specific star formation rate to calibrate the efficiency of feedback in a model for the high redshift star formation rate which includes merger driven star formation regulated by SNe feedback. This model reproduces the star formation rate density function and the stellar mass function of galaxies at 4<z<7. A prediction of the model is that the specific star formation rate does not evolve very rapidly with either mass or redshift, in agreement with observation. This is in contrast to results from hydrodynamical simulations where the star formation closely follows the accretion rate, and so increases strongly towards high redshift. The bursty star formation model naturally explains the observation that there is not enough stellar mass at z~2-4 to account for all of the star-formation observed, without invoking properties like an evolving initial mass function of stars. The finding of a duty cycle that is ~10% implies that there should be ten times the number of known galaxies at fixed stellar mass that have not yet been detected through standard UV selection at high redshift. We therefore predict the existence of a large undetected population of UV-faint galaxies that accounts for most of the stellar mass density at z~4-8.

preprint2013arXiv

Constraining Primordial Black-Hole Bombs through Spectral Distortions of the Cosmic Microwave Background

We consider the imprint of superradiant instabilities of nonevaporating primordial black holes (PBHs) on the spectrum of the cosmic microwave background (CMB). In the radiation dominated era, PBHs are surrounded by a roughly homogeneous cosmic plasma which endows photons with an effective mass through the plasma frequency. In this setting, spinning PBHs are unstable to a spontaneous spindown through the well-known "black-hole bomb" mechanism. At linear level, the photon density is trapped by the effective photon mass and grows exponentially in time due to superradiance. As the plasma density declines due to cosmic expansion, the associated energy around PBHs is released and dissipated in the CMB. We evaluate the resulting spectral distortions of the CMB in the redshift range 10^3 < z < 2x10^6. Using the existing COBE/FIRAS bounds on CMB spectral distortions, we derive upper limits on the fraction of dark matter that can be associated with spinning PBHs in the mass range 10^{-8}*Msun < M < 0.2*Msin. For maximally-spinning PBHs, our limits are much tighter than those derived from microlensing or other methods. Future data from the proposed PIXIE mission could improve our limits by several orders of magnitude.

preprint2013arXiv

Constraining sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. I. The General Quasar Population

We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]

preprint2013arXiv

Constraining the Minimum Luminosity of High Redshift Galaxies through Gravitational Lensing

We simulate the effects of gravitational lensing on the source count of high redshift galaxies as projected to be observed by the Hubble Frontier Fields program and the James Webb Space Telescope (JWST) in the near future. Taking the mass density profile of the lensing object to be the singular isothermal sphere (SIS) or the Navarro-Frenk-White (NFW) profile, we model a lens residing at a redshift of z_L = 0.5 and explore the radial dependence of the resulting magnification bias and its variability with the velocity dispersion of the lens, the photometric sensitivity of the instrument, the redshift of the background source population, and the intrinsic maximum absolute magnitude (M_{max}) of the sources. We find that gravitational lensing enhances the number of galaxies with redshifts z >= 13 detected in the angular region θ_E/2 <= θ<= 2θ_E (where θ_E is the Einstein angle) by a factor of ~ 3 and 1.5 in the HUDF (df/dν_0 ~ 9 nJy) and medium-deep JWST surveys (df/dν_0 ~ 6 nJy). Furthermore, we find that even in cases where a negative magnification bias reduces the observed number count of background sources, the lensing effect improves the sensitivity of the count to the intrinsic faint-magnitude cut-off of the Schechter luminosity function. In a field centered on a strong lensing cluster, observations of z >= 6 and z >= 13 galaxies with JWST can be used to infer this cut-off magnitude for values as faint as M_{max} ~ -14.4 and -16.1 mag (L_{min} ~ 2.5*10^{26} and 1.2*10^{27} erg s^{-1} Hz^{-1}) respectively, within the range bracketed by existing theoretical models. Gravitational lensing may therefore offer an effective way of constraining the low-luminosity cut-off of high-redshift galaxies.

preprint2013arXiv

Detecting bio-markers in habitable-zone earths transiting white dwarfs

The characterization of the atmospheres of habitable-zone Earth-mass exoplanets that transit across main-sequence stars, let alone the detection of bio-markers in their atmospheres, will be challenging even with future facilities. It has been noted that white dwarfs (WDs) have long-lived habitable zones and that a large fraction of WDs may host planets. We point out that during a transit of an Earth-mass planet across a WD, the planet's atmospheric transmission spectrum obtains a much higher contrast over the stellar background compared to a main-sequence host, because of the small surface area of the WD. The most prominent bio-marker in the present-day terrestrial atmosphere, molecular oxygen, is readily detectable in a WD transit via its A-band absorption at ~0.76 micron. A potentially life-sustaining Earth-like planet transiting a WD can be found by assembling a suitable sample of ~500 WDs and then surveying them for transits using small telescopes. If and when such a transiting case is found, the O_2 absorption in the planetary atmospheric transmission spectrum would be detectable with the James Webb Space Telescope (JWST) in about 5 hours of total exposure time, integrated over 160 2-minute transits. Characterization of the planet atmosphere using other tracers such as water vapour and CO_2 will be considerably easier. We demonstrate this future discovery space by simulating a possible transmission spectrum that would be obtained with JWST.

preprint2013arXiv

Fast radio bursts may originate from nearby flaring stars

Six cases of fast radio bursts (FRBs) have recently been discovered. The FRBs are bright (~0.1 - 1 Jy) and brief (~ 1 ms) pulses of radio emission with dispersion measures (DMs) that exceed Galactic values, and hence FRBs have been interpreted to be at cosmological distances. We propose, instead, that FRBs are rare eruptions of flaring main-sequence stars within ~1 kpc. Rather than associating their excess DM with the intergalactic medium, we relate it to a blanket of coronal plasma around their host stars. We have monitored at optical bands the stars within the radio beams of three of the known FRBs. In one field, we find a bright (V=13.6 mag) variable star (0.2 mag peak-to-trough) with a main-sequence G-type spectrum and a period P = 7.8 hr, likely a W-UMa-type contact binary. Analysis of our data outside of the FRB beams indicates a 5% chance probability of finding a variable star of this brightness and amplitude within the FRB beams. We find no unusual variable stars in the other two FRB fields. Further observations are needed to investigate if similar nearby (<~ 800 pc) stars are the sources of FRBs.

preprint2013arXiv

Finding Core Collapse Supernova from the Epoch of Reionization Behind Cluster Lenses

Current surveys are underway to utilize gravitational lensing by galaxy clusters with Einstein radii >35" in the search for the highest redshift galaxies. Associated supernova from the epoch of reionization would have their fluxes boosted above the detection threshold, extending their duration of visibility. We predict that the James Webb Space Telescope (JWST) will be able to discover lensed core-collapse supernovae at redshifts exceeding z=7-8.

preprint2013arXiv

Global 21cm signal experiments: A designer's guide

[Abridged] The spatially averaged global spectrum of the redshifted 21cm line has generated much experimental interest, for it is potentially a direct probe of the Epoch of Reionization and the Dark Ages. Since the cosmological signal here has a purely spectral signature, most proposed experiments have little angular sensitivity. This is worrisome because with only spectra, the global 21cm signal can be difficult to distinguish from foregrounds such as Galactic synchrotron radiation, as both are spectrally smooth and the latter is orders of magnitude brighter. We establish a mathematical framework for global signal data analysis in a way that removes foregrounds optimally, complementing spectra with angular information. We explore various experimental design trade-offs, and find that 1) with spectral-only methods, it is impossible to mitigate errors that arise from uncertainties in foreground modeling; 2) foreground contamination can be significantly reduced for experiments with fine angular resolution; 3) most of the statistical significance in a positive detection during the Dark Ages comes from a characteristic high-redshift trough in the 21cm brightness temperature; and 4) Measurement errors decrease more rapidly with integration time for instruments with fine angular resolution. We show that if observations and algorithms are optimized based on these findings, an instrument with a 5 degree beam can achieve highly significant detections (greater than 5-sigma) of even extended (high Delta-z) reionization scenarios after integrating for 500 hrs. This is in contrast to instruments without angular resolution, which cannot detect gradual reionization. Abrupt ionization histories can be detected at the level of 10-100's of sigma. The expected errors are also low during the Dark Ages, with a 25-sigma detection of the expected cosmological signal after only 100 hrs of integration.

preprint2013arXiv

Habitable Planets Around White Dwarfs: an Alternate Mission for the Kepler Spacecraft

A large fraction of white dwarfs (WDs) may host planets in their habitable zones. These planets may provide our best chance to detect bio-markers on a transiting exoplanet, thanks to the diminished contrast ratio between the Earth-sized WD and its Earth-sized planets. The JWST is capable of obtaining the first spectroscopic measurements of such planets, yet there are no known planets around WDs. Here we propose to take advantage of the unique capability of the Kepler spacecraft in the 2-Wheels mode to perform a transit survey that is capable of identifying the first planets in the habitable zone of a WD. We propose to obtain Kepler time-series photometry of 10,000 WDs in the SDSS imaging area to search for planets in the habitable zone. Thanks to the large field of view of Kepler, for the first time in history, a large number of WDs can be observed at the same time, which is essential for discovering transits. Our proposed survey requires a total of 200 days of observing time, and will find up to 100 planets in the WD habitable zone. This survey will maintain Kepler's spirit of searching for habitable Earths, but near new hosts. With few-day observations and minute-cadences per field, it will also open up a completely unexplored discovery space. In addition to planets, this survey is sensitive to pulsating WDs, as well as eclipsing short period stellar and substellar companions. These have important implications for constraining the double WD merger rate and their contribution to Type Ia supernovae and the gravitational wave foreground. Given the relatively low number density of our targets, this program can be combined with other projects that would benefit from high cadence and many-fields observations with Kepler, e.g. a transit survey of a magnitude-limited, complete sample of nearby M dwarfs or asteroseismology of variable stars (e.g. RR Lyrae) in the same fields.

preprint2013arXiv

High-Angular-Resolution and High-Sensitivity Science Enabled by Beamformed ALMA

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

preprint2013arXiv

High-n Hydrogen Recombination Lines from the First Galaxies

We investigate the prospects of blind and targeted searches in the radio domain (10 MHz to 1 THz) for high-n hydrogen recombination lines from the first generation of galaxies, at z < 10. The expected optically thin spontaneous alpha-line luminosities are calculated as a function of the absolute AB magnitude of a galaxy at 1500 angstrom. For a blind search, semi-empirical luminosity functions are used to calculate the number of galaxies whose expected flux densities exceed an assumed detectability threshold. Plots of the minimum sky area, within which at least one detectable galaxy is expected at a given observing frequency, in the fiducial instantaneous passband of 10^4 km/s, allow to assess the blind search time necessary for detection by a given facility. We show that the chances for detection are the highest in the mm and submm domains, but finding spontaneous emission in a blind search, especially from redshifts z >> 1, is a challenge even with powerful facilities, such as ALMA and SKA. The probability of success is higher for a targeted search of lines with principal quantum number n ~ 10 in Lyman-break galaxies amplified by gravitational lensing. Detection of more than one hydrogen line in such a galaxy will allow for line identification and a precise determination of the galaxy's redshift.

preprint2013arXiv

Inflationary paradigm in trouble after Planck2013

Recent results from the Planck satellite combined with earlier observations from WMAP, ACT, SPT and other experiments eliminate a wide spectrum of more complex inflationary models and favor models with a single scalar field, as reported by the Planck Collaboration. More important, though, is that all the simplest inflaton models are disfavored statistically relative to those with plateau-like potentials. We discuss how a restriction to plateau-like models has three independent serious drawbacks: it exacerbates both the initial conditions problem and the multiverse-unpredictability problem and it creates a new difficulty that we call the inflationary "unlikeliness problem." Finally, we comment on problems reconciling inflation with a standard model Higgs, as suggested by recent LHC results. In sum, we find that recent experimental data disfavors all the best-motivated inflationary scenarios and introduces new, serious difficulties that cut to the core of the inflationary paradigm. Forthcoming searches for B-modes, non-Gaussianity and new particles should be decisive.

preprint2013arXiv

Measuring the X-ray Background in the Reionization Era with First Generation 21 cm Experiments

The X-ray background during the epoch of reionization is currently poorly constrained. We demonstrate that it is possible to use first generation 21 cm experiments to calibrate it. Using the semi-numerical simulation, 21cmFAST, we calculate the dependence of the 21 cm power spectrum on the X-ray background flux. Comparing the signal to the sensitivity of the Murchison Widefield Array (MWA) we find that in the redshift interval z=8-14 the 21 cm signal is detectable based on the upper limit set by the present-day unresolved soft X-ray background. We show that there is no degeneracy between the X-ray production efficiency and the Lyman-Alpha production efficiency and that the degeneracy with the ionization fraction of the intergalactic medium can be broken.

preprint2013arXiv

Menus for Feeding Black Holes

Black holes are the ultimate prisons of the Universe, regions of spacetime where the enormous gravity prohibits matter or even light to escape to infinity. Yet, matter falling toward the black holes may shine spectacularly, generating the strongest source of radiation. These sources provide us with astrophysical laboratories of extreme physical conditions that cannot be realized on Earth. This chapter offers a review of the basic menus for feeding matter onto black holes and discusses their observational implications.

preprint2013arXiv

Neutrino Signatures on the High Transmission Regions of the Lyman-alpha Forest

We quantify the impact of massive neutrinos on the statistics of low density regions in the intergalactic medium (IGM) as probed by the Lyman-alpha forest at redshifts z=2.2--4. Based on mock but realistic quasar (QSO) spectra extracted from hydrodynamic simulations with cold dark matter, baryons and neutrinos, we find that the probability distribution of weak Lyman-alpha absorption features, as sampled by Lyman-alpha flux regions at high transmissivity, is strongly affected by the presence of massive neutrinos. We show that systematic errors affecting the Lyman-alpha forest reduce but do not erase the neutrino signal. Using the Fisher matrix formalism, we conclude that the sum of the neutrino masses can be measured, using the method proposed in this paper, with a precision smaller than 0.4 eV using a catalog of 200 high resolution (S/N~100) QSO spectra. This number reduces to 0.27 eV by making use of reasonable priors in the other parameters that also affect the statistics of the high transitivity regions of the Lyman-alpha forest. The constraints obtained with this method can be combined with independent bounds from the CMB, large scale structures and measurements of the matter power spectrum from the Lyman-alpha forest to produce tighter upper limits on the sum of the masses of the neutrinos.

preprint2013arXiv

On the Importance of Conceptual Thinking Outside the Simulation Box

Any ambitious construction project requires architects for its design and engineers who apply the design to the real world. As scientific research shifts towards large groups which focus on the engineering aspects of linking data to existing models, architectural skills are becoming rare among young theorists. Senior researchers should mentor qualified students and postdocs to think creatively about the big picture without unwarranted loyalty to ancient blueprints from past generations of architects.

preprint2013arXiv

Reionization on Large Scales II: Detecting Patchy Reionization through Cross Correlation of the Cosmic Microwave Background

We investigate the effect of patchy reionization on the cosmic microwave background temperature. An anisotropic optical depth tau (theta) alters the TT power spectrum on small scales l > 2000. We make use of the correlation between the matter density and the reionization redshift fields to construct full sky maps of tau(theta). Patchy reionization transfers CMB power from large scales to small scales, resulting in a non-zero cross correlation between large and small angular scales. We show that the patchy tau correlator is sensitive to small root mean square values tau_rms ~ 0.003 seen in our maps. We include other secondary anisotropies such as CMB lensing, kinetic and thermal Sunyaev-Zel'dovich terms, as well as the infrared and point source background, and show that patchy reionization may be detected in the low frequency channels ~ 90 GHz, particularly for extended reionization histories. If frequency dependent secondaries can be minimized by a multi-frequency analysis, we show that even small degrees of patchiness can be detected at high significance. In this case, we show that models with different values of tau_rms are distinguishable, even for the same mean value <tau>. The patchy tau correlator thus provides information that is complementary to what may be obtained from the polarization and the kinetic Sunyaev-Zel'dovich power spectra.

preprint2013arXiv

Resonant Post-Newtonian Eccentricity Excitation in Hierarchical Three-body Systems

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

preprint2013arXiv

Scale-free primordial cosmology

The large-scale structure of the universe suggests that the physics underlying its early evolution is scale-free. This was the historic motivation for the Harrison-Zel'dovich-Peebles spectrum and for inflation. Based on a hydrodynamical approach, we identify scale-free forms for the background equation-of-state for both inflationary and cyclic scenarios and use these forms to derive predictions for the spectral tilt and tensor-to-scalar ratio of primordial density perturbations. For the case of inflation, we find three classes of scale-free models with distinct predictions. Including all classes, we show that scale-free inflation predicts tensor-to-scalar ratio $r > 10^{-4}$. We show that the observationally favored class is theoretically disfavored because it suffers from an initial conditions problem and the hydrodynamical form of an unlikeliness problem similar to that identified recently for certain inflaton potentials. We contrast these results with those for scale-free cyclic models.

preprint2013arXiv

Spurious Velocities in Dynamically-Cold Systems Due to the Gravitational Redshifts of their Constituent Stars

The Doppler effect is commonly used to infer the velocity difference between stars based on the relative shifts in the rest-frame wavelengths of their spectral features. In dynamically-cold systems with a low velocity dispersion, such as wide binaries, loose star clusters, cold stellar streams or cosmological mini-halos, the scatter in the gravitational redshift from the surface of the constituent stars needs to be taken into account as well. Gravitational redshifts could be important for wide binaries composed of main sequence stars with separations >0.01 pc or in mini-halos with velocity dispersions <1 km/s. Variable redshift could also lead to a spurious "detection" of low-mass planets around a star with periodic photospheric radius variations.

preprint2013arXiv

Super-luminous X-ray Emission from the Interaction of Supernova Ejecta with Dense Circumstellar Shells

For supernova powered by the conversion of kinetic energy into radiation due to the interactions of the ejecta with a dense circumstellar shell, we show that there could be X-ray analogues of optically super-luminous SNe with comparable luminosities and energetics. We consider X-ray emission from the forward shock of SNe ejecta colliding into an optically-thin CSM shell, derive simple expressions for the X-ray luminosity as a function of the circumstellar shell characteristics, and discuss the different regimes in which the shock will be radiative or adiabatic, and whether the emission will be dominated by free-free radiation or line-cooling. We find that even with normal supernova explosion energies of 10^51 erg, there exists CSM shell configurations that can liberate a large fraction of the explosion energy in X-rays, producing unabsorbed X-ray luminosities approaching 10^44 erg/s events lasting a few months, or even 10^45 erg/s flashes lasting days. Although the large column density of the circumstellar shell can absorb most of the flux from the initial shock, the most luminous events produce hard X-rays that are less susceptible to photoelectric absorption, and can counteract such losses by completely ionizing the intervening material. Regardless, once the shock traverses the entire circumstellar shell, the full luminosity could be available to observers.

preprint2013arXiv

Testing the No-Hair Theorem with Event Horizon Telescope Observations of Sagittarius A*

The advent of the Event Horizon Telescope (EHT), a millimeter-wave very-long baseline interferometric array, has enabled spatially-resolved studies of the sub-horizon-scale structure for a handful of supermassive black holes. Among these, the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), presents the largest angular cross section. Thus far, these studies have focused upon measurements of the black hole spin and the validation of low-luminosity accretion models. However, a critical input into the analysis of EHT data is the structure of the black hole spacetime, and thus these observations provide the novel opportunity to test the applicability of the Kerr metric to astrophysical black holes. Here we present the first simulated images of a radiatively inefficient accretion flow (RIAF) around Sgr A* employing a quasi-Kerr metric that contains an independent quadrupole moment in addition to the mass and spin that fully characterize a black hole in general relativity. We show that these images can be significantly different from the images of a RIAF around a Kerr black hole with the same spin and demonstrate the feasibility of testing the no-hair theorem by constraining the quadrupolar deviation from the Kerr metric with existing EHT data. Equally important, we find that the disk inclination and spin orientation angles are robust to the inclusion of additional parameters, providing confidence in previous estimations assuming the Kerr metric based upon EHT observations. However, at present the limits upon potential modifications of the Kerr metric remain weak.

preprint2013arXiv

The 21-cm Signal from the Cosmological Epoch of Recombination

The redshifted 21-cm emission by neutral hydrogen offers a unique tool for mapping structure formation in the early universe in three dimensions. Here we provide the first detailed calculation of the 21-cm emission signal during and after the epoch of hydrogen recombination in the redshift range of z ~ 500-1,100, corresponding to observed wavelengths of 100-230 meters. The 21-cm line deviates from thermal equilibrium with the cosmic microwave background (CMB) due to the excess Ly-alpha radiation from hydrogen and helium recombinations. The resulting 21-cm signal reaches a brightness temperature of a milli-Kelvin, orders of magnitude larger than previously estimated. Its detection by a future lunar or space-based observatory could improve dramatically the statistical constraints on the cosmological initial conditions compared to existing two-dimensional maps of the CMB anisotropies.

preprint2013arXiv

The energy production rate density of cosmic rays in the local universe is $\sim10^{44-45}\rm erg~Mpc^{-3}~yr^{-1}$ at all particle energies

The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe. It is only $\sim 10$ times larger than the output, $\varepsilon_p^2 Q=0.5\pm0.2\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$, of Ultra High Energy CRs (UHECRs) at energies $10^{10.5}{\rm GeV}<\varepsilon_p<10^{12}\rm GeV$ (obtained assuming they are mostly protons), which in turn is comparable to the lower limit of $\varepsilon_p^2 Q\ge 0.5\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$ of high energy CRs with $10^6{\rm GeV}\lesssim\varepsilon_p\lesssim 10^{8}\rm GeV$ implied by the saturation of the Waxman-Bahcall bound by the neutrino excess recently discovered by IceCube. These similarities are consistent with a flat production spectrum, $\varepsilon_p^2 Q\sim \text{const}$ for CRs at all observed energies. If a flat production spectrum is generated by our galaxy, the observed CR flux in the range $10^{6.5}-10^{9.5}\rm GeV$, above the "knee", is suppressed compared to lower energies due to propagation effects rather than acceleration upper limits. As suggested by Parizot and Aublin, the most exciting possibility is that cosmic rays at all energies are emitted from a single type of (unknown) sources, which can not be supernova remnants.

preprint2013arXiv

The Ratio of CO to Total Gas Mass in High Redshift Galaxies

Walter et al. (20012) have recently identified the J=6-5, 5-4, and 2-1 CO rotational emission lines, and [C_{II}] fine-structure emission line from the star-forming interstellar medium in the high-redshift submillimeter source HDF 850.1, at z = 5.183. We employ large velocity gradient (LVG) modeling to analyze the spectra of this source assuming the [C_{II}] and CO emissions originate from (i) separate unvirialized regions, (ii) separate virialized regions, (iii) uniformly mixed unvirialized region, and (iv) uniformly mixed virialized regions. We present the best fit set of parameters, including for each case the ratio $α$ between the total hydrogen/helium gas mass and the CO(1-0) line luminosity. We also present computations of the ratio of H_{2} mass to [C_{II}] line-luminosity for optically thin conditions, for a range of gas temperatures and densities, for direct conversion of [C_{II}] line-luminosities to "dark-H_{2}" masses. For HDF 850.1 we find that a model in which the CO and C^{+} are uniformly mixed in gas that is shielded from UV radiation, requires a cosmic-ray or X-ray ionization rate of $ζ\approx$ 10^{-13} s^{-1}, plausibly consistent with the large star-formation rate ($\sim$ 10^{3} M$_{\odot}$ yr^{-1}) observed in this source. Enforcing the cosmological constraint posed by the abundance of dark matter halos in the standard $Λ$CDM cosmology and taking into account other possible contributions to the total gas mass, we find that three of these four models are less likely at the 2$σ$ level. We conclude that modeling HDF 850.1's ISM as a collection of unvirialized molecular clouds with distinct CO and C^{+} layers, for which $α$ = 0.6 M$_{\odot}$ (K km s^{-1} pc^{2})^{-1} for the CO to H_{2} mass-to-luminosity ratio, (similar to the standard ULIRG value), is most consistent with the $Λ$CDM cosmology.

preprint2013arXiv

The virialization density of peaks with general density profiles under spherical collapse

We calculate the non-linear virialization density, $Δ_c$, of halos under spherical collapse from peaks with an arbitrary initial and final density profile. This is in contrast to the standard calculation of $Δ_c$ which assumes top-hat profiles. Given our formalism, the non-linear halo density can be calculated once the shape of the initial peak's density profile and the shape of the virialized halo's profile are provided. We solve for $Δ_c$ for halos in an Einstein de-Sitter and $Λ$CDM universe. As examples, we consider power-law initial profiles as well as spherically averaged peak profiles calculated from the statistics of a Gaussian random field. We find that, depending on the profiles used, $Δ_c$ is smaller by a factor of a few to as much as a factor of 10 as compared to the density given by the standard calculation ($\approx 200$). Using our results, we show that, for halo finding algorithms that identify halos through an over-density threshold, the halo mass function measured from cosmological simulations can be enhanced at all halo masses by a factor of a few. This difference could be important when using numerical simulations to assess the validity of analytic models of the halo mass function.

preprint2013arXiv

Unbound Particles in Dark Matter Halos

We investigate unbound dark matter particles in halos by tracing particle trajectories in a simulation run to the far future (a = 100). We find that the traditional sum of kinetic and potential energies is a very poor predictor of which dark matter particles will eventually become unbound from halos. We also study the mass fraction of unbound particles, which increases strongly towards the edges of halos, and decreases significantly at higher redshifts. We discuss implications for dark matter detection experiments, precision calibrations of the halo mass function, the use of baryon fractions to constrain dark energy, and searches for intergalactic supernovae.

preprint2012arXiv

21-cm cosmology

Imaging the Universe during the first hundreds of millions of years remains one of the exciting challenges facing modern cosmology. Observations of the redshifted 21 cm line of atomic hydrogen offer the potential of opening a new window into this epoch. This would transform our understanding of the formation of the first stars and galaxies and of the thermal history of the Universe. A new generation of radio telescopes is being constructed for this purpose with the first results starting to trickle in. In this review, we detail the physics that governs the 21 cm signal and describe what might be learnt from upcoming observations. We also generalize our discussion to intensity mapping of other atomic and molecular lines.

preprint2012arXiv

A Simple Model for the Density Profiles of Isolated Dark Matter Halos

We explore the possibility that the density profiles of elliptical galaxies and cold dark matter (CDM) halos found in cosmological simulations can be understood in terms of the same physical process, collisionless gravitational collapse. To investigate this, we study a simplified model, the collapse of a perfectly cold Plummer sphere. First, we examine an N-body simulation of this model with particles constrained to purely radial orbits. This results in a final state characterized by a profile slightly steeper than ρ\propto r^{-2} at small radii and behaving as ρ\propto r^{-4} at large radii, which can be understood in terms of simple analytic arguments. Next, we repeat our simulation without the restriction of radial orbits. This results in a shallower inner density profile, like those found in elliptical galaxies and CDM halos. We attribute this change to the radial orbit instability (ROI) and propose a form of the distribution function (DF) motivated by a physical picture of collapse. As evidence of the link between our model and CDM halos, we find that our collapse simulation has a final state with pseudo-phase-space density which scales roughly as ρ/σ^3 \propto r^{-1.875}, like that observed in CDM halos from cosmological simulations (Navarro et al. 2010). The velocity anisotropy profile is also qualitatively similar to that found near the centers of these halos. We argue that the discrepancy at large radii (where CDM halos scale as ρ\propto r^{-3}) is due to the presence of the cosmological background or continued infall. This leads us to predict that the outer CDM halo density profile is not "universal," but instead depends on cosmological environment (be it an underdense void or overdense region).

preprint2012arXiv

A suppressed contribution of low mass galaxies to reionization due to supernova feedback

Motivated by recent observations of the star formation rate density function out to z~7, we describe a simple model for the star formation rate density function at high redshift based on the extended Press-Schechter formalism. This model postulates a starburst following each major merger, lasting for a time t_SF and converting at most f_star of galactic gas into stars. We include a simple physical prescription for supernovae feedback that suppresses star formation in low mass galaxies. Constraining t_SF and f_star to describe the observed star formation rate density at high redshifts, we find that individual starbursts were terminated after a time t_SF~10^7 years. This is comparable to the main-sequence lifetimes of supernova progenitors, indicating that high redshift starbursts are quenched once supernovae feedback had time to develop. High redshift galaxies convert ~10% of their mass into stars for galaxies with star formation rates above ~1 solar mass per year, but a smaller fraction for lower luminosity galaxies. Our best fit model successfully predicts the observed relation between star formation rate and stellar mass at z>~4, while our deduced relation between stellar mass and halo mass is also consistent with data on the dwarf satellites of the Milky Way. We find that supernovae feedback lowers the efficiency of star formation in the lowest mass galaxies and makes their contribution to reionization small. As a result, photo-ionization feedback on low mass galaxy formation does not significantly affect the reionization history. Using a semi-analytic model for the reionization history, we infer that approximately half of the ionizing photons needed to complete reionization have already been observed in star-forming galaxies.

preprint2012arXiv

Accumulated Tidal Heating of Stars Over Multiple Pericenter Passages Near SgrA*

We consider the long-term tidal heating of a star by the supermassive black hole at the Galactic center, SgrA*. We show that gravitational interaction with background stars leads to a linear growth of the tidal excitation energy with the number of pericenter passages near SgrA*. The accumulated heat deposited by excitation of modes within the star over many pericenter passages can lead to a runaway disruption of the star at a pericenter distance that is 4-5 times farther than the standard tidal disruption radius. The accumulated heating may explain the lack of massive ($\gtrsim 10M_{\odot}$) S-stars closer than several tens of AU from SgrA*.

preprint2012arXiv

Consequences of Strong Compression in Tidal Disruption Events

The tidal disruption of a star by a supermassive black hole (SMBH) is a highly energetic event with consequences dependent on the degree to which the star plunges inside the SMBH's tidal sphere. We introduce a new analytic model for tidal disruption events (TDEs) to analyze the dependence of these events on beta, the ratio of the tidal radius to the orbital pericenter. We find, contrary to most previous work, that the spread in debris energy for a TDE is largely constant for all beta. This result has important consequences for optical transient searches targeting TDEs, which we discuss. We quantify leading-order general relativistic corrections to this spread in energy and find that they are small. We also examine the role of stellar spin, and find that a combination of spin-orbit misalignment, rapid rotation, and high beta may increase the spread in debris energy. Finally, we quantify for the first time the gravitational wave emission due to the strong compression of a star in a high-beta TDE. Although this signal is unlikely to be detectable for disruptions of main sequence stars, the tidal disruption of a white dwarf by an intermediate mass black hole can produce a strong signal visible to Advanced LIGO at tens of megaparsecs.

preprint2012arXiv

Constraining Sources of Ultra High Energy Cosmic Rays Using High Energy Observations with the Fermi Satellite

We analyze the conditions that enable acceleration of particles to ultra-high energies, ~10^{20} eV (UHECRs). We show that broad band photon data recently provided by WMAP, ISOCAM, Swift and Fermi satellites, yield constraints on the ability of active galactic nuclei (AGN) to produce UHECRs. The high energy (MeV - GeV) photons are produced by Compton scattering of the emitted low energy photons and the cosmic microwave background or extra-galactic background light. The ratio of the luminosities at high and low photon energies can therefore be used as a probe of the physical conditions in the acceleration site. We find that existing data excludes core regions of nearby radio-loud AGN as possible acceleration sites of UHECR protons. However, we show that giant radio lobes are not excluded. We apply our method to Cen A, and show that acceleration of protons to ~10^{20} eV can only occur at distances >~ 100 kpc from the core.

preprint2012arXiv

Constraints on the Nature of CID-42: Recoil Kick or Supermassive Black Hole Pair?

The galaxy CXOC J100043.1+020637, also known as CID-42, is a highly unusual object. An apparent galaxy merger remnant, it displays signatures of both an inspiraling, kiloparsec-scale active galactic nucleus (AGN) pair and of a recoiling AGN with a kick velocity > 1300 km s^-1. Among recoiling AGN candidates, CID-42 alone has both spatial offsets (in optical and X-ray bands) and spectroscopic offsets. In order to constrain the relative likelihood of both scenarios, we develop models using hydrodynamic galaxy merger simulations coupled with radiative transfer calculations. Our gas-rich, major merger models are generally well matched to the galactic morphology and to the inferred stellar mass and star formation rate. We show that a recoiling supermassive black hole (SMBH) in CID-42 should be observable as an AGN at the time of observation. However, in order for the recoiling AGN to produce narrow-line emission, it must be observed shortly after the kick while it still inhabits a dense gaseous region, implying a large total kick velocity (v_k > 2000 km s^-1). For the dual AGN scenario, an unusually large broad-line offset is required, and the best match to the observed morphology requires a galaxy that is less luminous than CID-42. Further, the lack of X-ray emission from one of the two optical nuclei is not easily attributed to an intrinsically quiescent SMBH or to a Compton-thick galactic environment. While the current data do not allow either the recoiling or the dual AGN scenario for CID-42 to be excluded, our models highlight the most relevant parameters for distinguishing these possibilities with future observations. In particular, high-quality, spatially-resolved spectra that can pinpoint the origin of the broad and narrow line features will be critical for determining the nature of this unique source.

preprint2012arXiv

Detection Technique for Artificially-Illuminated Objects in the Outer Solar System and Beyond

Existing and planned optical telescopes and surveys can detect artificially-illuminated objects comparable in total brightness to a major terrestrial city out to the outskirts of the Solar System. Orbital parameters of Kuiper belt objects (KBOs) are routinely measured to exquisite precisions of <10^{-3}. Here we propose to measure the variation of the observed flux F from such objects as a function of their changing orbital distances D. Sunlight-illuminated objects will show a logarithmic slope alpha=(dlogF/dlog D)=-4 whereas artificially-illuminated objects should exhibit alpha=-2. Planned surveys using the proposed LSST will provide superb data that would allow measurement of alpha for thousands of KBOs. If objects with alpha=-2 are found, follow-up observations can measure their spectra to determine if they are illuminated by artificial lighting. The search can be extended beyond the Solar System with future generations of telescopes on the ground and in space, which would be capable of detecting phase modulation due to very strong artificial illumination on the night-side of planets as they orbit their parent stars.

preprint2012arXiv

Deviation of Stellar Orbits from Test Particle Trajectories Around Sgr A* Due to Tides and Winds

Monitoring the orbits of stars around Sgr A* offers the possibility of detecting the precession of their orbital planes due to frame dragging, of measuring the spin and quadrupole moment of the black hole, and of testing the no-hair theorem. Here we investigate whether the deviations of stellar orbits from test-particle trajectories due to wind mass loss and tidal dissipation of the orbital energy compromise such measurements. We find that the effects of stellar winds are, in general, negligible. On the other hand, for the most eccentric orbits (e>0.96) for which an optical interferometer, such as GRAVITY, will detect orbital plane precession due to frame dragging, the tidal dissipation of orbital energy occurs at timescales comparable to the timescale of precession due to the quadrupole moment of the black hole. As a result, this non-conservative effect is a potential source of systematic uncertainty in testing the no-hair theorem with stellar orbits.

preprint2012arXiv

Disruption of a Proto-Planetary Disk by the Black Hole at the Milky Way Centre

Recently, an ionized cloud of gas was discovered plunging toward the supermassive black hole, SgrA*, at the centre of the Milky Way. The cloud is being tidally disrupted along its path to closest approach at ~3100 Schwarzschild radii from the black hole. Here, we show that the observed properties of this cloud of gas can naturally be produced by a proto-planetary disk surrounding a low-mass star, which was scattered from the observed ring of young stars orbiting SgrA*. As the young star approaches the black hole, its disk experiences both photo-evaporation and tidal disruption, producing a cloud. Our model implies that planets form in the Galactic centre, and that tidal debris from proto-planetary disks can flag low mass stars which are otherwise too faint to be detected.

preprint2012arXiv

Double-peaked Narrow-Line Signatures of Dual Supermassive Black Holes in Galaxy Merger Simulations

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN should often be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)

preprint2012arXiv

Early structure formation from cosmic string loops

We examine the effects of cosmic strings on structure formation and on the ionization history of the universe. While Gaussian perturbations from inflation are known to provide the dominant contribution to the large scale structure of the universe, density perturbations due to strings are highly non-Gaussian and can produce nonlinear structures at very early times. This could lead to early star formation and reionization of the universe. We improve on earlier studies of these effects by accounting for high loop velocities and for the filamentary shape of the resulting halos. We find that for string energy scales Gμ> 10^{-7} the effect of strings on the CMB temperature and polarization power spectra can be significant and is likely to be detectable by the Planck satellite. We mention shortcomings of the standard cosmological model of galaxy formation which may be remedied with the addition of cosmic strings, and comment on other possible observational implications of early structure formation by strings.

preprint2012arXiv

Effect of Population III Multiplicity on Dark Star Formation

We numerically study the mutual interaction between dark matter (DM) and Population III (Pop III) stellar systems in order to explore the possibility of Pop III dark stars within this physical scenario. We perform a cosmological simulation, initialized at z ~ 100, which follows the evolution of gas and DM. We analyze the formation of the first minihalo at z ~ 20 and the subsequent collapse of the gas to densities of 10^12 cm^-3. We then use this simulation to initialize a set of smaller-scale `cut-out' simulations in which we further refine the DM to have spatial resolution similar to that of the gas. We test multiple DM density profiles, and we employ the sink particle method to represent the accreting star-forming region. We find that, for a range of DM configurations, the motion of the Pop III star-disk system serves to separate the positions of the protostars with respect to the DM density peak, such that there is insufficient DM to influence the formation and evolution of the protostars for more than ~ 5000 years. In addition, the star-disk system causes gravitational scattering of the central DM to lower densities, further decreasing the influence of DM over time. Any DM-powered phase of Pop III stars will thus be very short-lived for the typical multiple system, and DM will not serve to significantly prolong the life of Pop III stars.

preprint2012arXiv

Finite, Intense Accretion Bursts from Tidal Disruption of Stars on Bound Orbits

We study accretion processes for tidally disrupted stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that there is a critical value of the orbital eccentricity below which all the stellar debris remains bound to the black hole. For high but sub-critical eccentricities, all the stellar mass is accreted onto the black hole in a finite time, causing a significant deviation from the canonical $t^{-5/3}$ mass fallback rate. When a star is on a moderately eccentric orbit and its pericenter distance is deeply inside the tidal disruption radius, there can be several orbit crossings of the debris streams due to relativistic precession. This dissipates orbital energy in shocks, allowing for rapid circularization of the debris streams and formation of an accretion disk. The resultant accretion rate greatly exceeds the Eddington rate and differs strongly from the canonical rate of $t^{-5/3}$. By contrast, there is little dissipation due to orbital crossings for the equivalent simulation with a purely Newtonian potential. This shows that general relativistic precession is crucial for accretion disk formation via circularization of stellar debris from stars on moderately eccentric orbits.

preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: application to supermassive black hole binaries

We study the interaction of a supermassive black hole (SMBH) binary and a standard radiatively efficient thin accretion disk. We examine steady-state configurations of the disk and migrating SMBH system, self-consistently accounting for tidal and viscous torques and heating, radiative diffusion limited cooling, gas and radiation pressure, and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a function of binary parameters, showing regimes in which both the disk structure and migration have a different character. Although massive binaries can create a central gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile-up of gas outside of its orbit, which can lead to the closing of the gap. We find that this spillover occurs at an orbital separation as large as ~200 M_7^{-1/2} gravitational radii, where M = 10^7 M_7 Msun is the total binary mass. If the secondary is less massive than ~10^6 Msun, then the gap is closed before gravitational waves (GWs) start dominating the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. The corresponding migration rate, which we label Type 1.5, is slower than the usual limiting cases known as Type I and II migration. Compared to an unperturbed disk, the steady-state disk in the overflowing regime is up to several hundred times brighter in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic variability of this population of binaries. Our results imply that the circumbinary disks around SMBHs can extend to small radii during the last stages of their merger, when they are detectable by LISA, and may produce coincident electromagnetic (EM) emission similar to active galactic nuclei (AGN).

preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: general theory

Many astrophysical binaries, from planets to black holes, exert strong torques on their circumbinary accretion disks, and are expected to significantly modify the disk structure. Despite the several decade long history of the subject, the joint evolution of the binary + disk system has not been modeled with self-consistent assumptions for arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evolution equations analytically in the strongly perturbed limit, treating the azimuthally-averaged angular momentum exchange between the disk and the binary and the modifications to the density, scale-height, and viscosity self-consistently, including viscous and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay of the binary. We find a solution with a central cavity and a migration rate similar to those previously obtained for Type-II migration, applicable for large masses and binary separations, and near-equal mass ratios. However, we identify a distinct new regime, applicable at smaller separations and masses, and mass ratio in the range 0.001< q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather than creating a cavity, it continuously overflows as in a porous dam. The disk profile is intermediate between a weakly perturbed disk (producing Type-I migration) and a disk with a gap (with Type-II migration). However, the migration rate of the secondary is typically slower than both Type-I and Type-II rates. We term this new regime "Type-1.5" migration.

preprint2012arXiv

Gravitational Wave Heating of Stars and Accretion Disks

We investigate the electromagnetic (EM) counterpart of gravitational waves (GWs) emitted by a supermassive black hole binary (SMBHB) through the viscous dissipation of the GW energy in an accretion disk and stars surrounding the SMBHB. We account for the suppression of the heating rate if the forcing period is shorter than the turnover time of the largest turbulent eddies. We find that the viscous heating luminosity in 0.1 solar mass stars can be significantly higher than their intrinsic luminosity, but still too low to be detected for extragalactic sources. The relative brightening is small for accretion disks.

preprint2012arXiv

How to Nurture Scientific Discoveries Despite Their Unpredictable Nature

The history of science reveals that major discoveries are not predictable. Naively, one might conclude therefore that it is not possible to artificially cultivate an environment that promotes discoveries. I suggest instead that open research without a programmatic agenda establishes a fertile ground for unexpected breakthroughs. Contrary to current practice, funding agencies should allocate a small fraction of their funds to support research in centers of excellence without programmatic reins tied to specific goals.

preprint2012arXiv

Hypervelocity Planets and Transits Around Hypervelocity Stars

The disruption of a binary star system by the massive black hole at the Galactic Centre, SgrA*, can lead to the capture of one star around SgrA* and the ejection of its companion as a hypervelocity star (HVS). We consider the possibility that these stars may have planets and study the dynamics of these planets. Using a direct $N$-body integration code, we simulated a large number of different binary orbits around SgrA*. For some orbital parameters, a planet is ejected at a high speed. In other instances, a HVS is ejected with one or more planets orbiting around it. In these cases, it may be possible to observe the planet as it transits the face of the star. A planet may also collide with its host star. In such cases the atmosphere of the star will be enriched with metals. In other cases, a planet is tidally disrupted by SgrA*, leading to a bright flare.

preprint2012arXiv

Identifying Stars of Mass >150 Msun from Their Eclipse by a Binary Companion

We examine the possibility that very massive stars greatly exceeding the commonly adopted stellar mass limit of 150 Msun may be present in young star clusters in the local universe. We identify ten candidate clusters, some of which may host stars with masses up to 600 Msun formed via runaway collisions. We estimate the probabilities of these very massive stars being in eclipsing binaries to be >30%. Although most of these systems cannot be resolved at present, their transits can be detected at distances of 3 Mpc even under the contamination of the background cluster light, due to the large associated luminosities ~10^7 Lsun and mean transit depths of ~10^6 Lsun. Discovery of very massive eclipsing binaries would flag possible progenitors of pair-instability supernovae and intermediate-mass black holes.

preprint2012arXiv

Jet Launching Structure Resolved Near the Supermassive Black Hole in M87

Approximately 10% of active galactic nuclei exhibit relativistic jets, which are powered by accretion of matter onto super massive black holes. While the measured width profiles of such jets on large scales agree with theories of magnetic collimation, predicted structure on accretion disk scales at the jet launch point has not been detected. We report radio interferometry observations at 1.3mm wavelength of the elliptical galaxy M87 that spatially resolve the base of the jet in this source. The derived size of 5.5 +/- 0.4 Schwarzschild radii is significantly smaller than the innermost edge of a retrograde accretion disk, suggesting that the M87 jet is powered by an accretion disk in a prograde orbit around a spinning black hole.

preprint2012arXiv

Observing Lense-Thirring Precession in Tidal Disruption Flares

When a star is tidally disrupted by a supermassive black hole (SMBH), the streams of liberated gas form an accretion disk after their return to pericenter. We demonstrate that Lense-Thirring precession in the spacetime around a rotating SMBH can produce significant time evolution of the disk angular momentum vector, due to both the periodic precession of the disk and the nonperiodic, differential precession of the bound debris streams. Jet precession and periodic modulation of disk luminosity are possible consequences. The persistence of the jetted X-ray emission in the Swift J164449.3+573451 flare suggests that the jet axis was aligned with the spin axis of the SMBH during this event.

preprint2012arXiv

Pair-Instability Supernovae at the Epoch of Reionization

Pristine stars with masses between ~140 and 260 M_sun are theoretically predicted to die as pair-instability supernovae. These very massive progenitors could come from Pop III stars in the early universe. We model the light curves and spectra of pair-instability supernovae over a range of masses and envelope structures. At redshifts of reionization z >= 6, we calculate the rates and detectability of pair-instability and core collapse supernovae, and show that with the James Webb Space Telescope, it is possible to determine the contribution of Pop III and Pop II stars toward reionization by constraining the stellar initial mass function at that epoch using these supernovae. We also find the rates of Type Ia supernovae, and show that they are not rare during reionization, and can be used to probe the mass function at 4-8 M_sun. If the budget of ionizing photons was dominated by contributions from top-heavy Pop III stars, we predict that the bright end of the galaxy luminosity function will be contaminated by pair-instability supernovae.

preprint2012arXiv

Pair-Instability Supernovae via Collision Runaway in Young Dense Star Clusters

Stars with helium cores between ~64 and 133 M_sun are theoretically predicted to die as pair-instability supernovae. This requires very massive progenitors, which are theoretically prohibited for Pop II/I stars within the Galactic stellar mass limit due to mass loss via line-driven winds. However, the runaway collision of stars in a dense, young star cluster could create a merged star with sufficient mass to end its life as a pair-instability supernova, even with enhanced mass loss at non-zero metallicity. We show that the predicted rate from this mechanism is consistent with the inferred volumetric rate of roughly ~2x10^-9 Mpc^-3 yr^-1 of the two observed pair-instability supernovae, SN 2007bi and PTF 10nmn, neither of which have metal-free host galaxies. Contrary to prior literature, only pair-instability supernovae at low redshifts z<2 will be observable with the Large Synoptic Survey Telescope (LSST). We estimate the telescope will observe ~10^2 such events per year that originate from the collisional runaway mergers in clusters.

preprint2012arXiv

Probing Pre-galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

preprint2012arXiv

Reionization on Large Scales I: A Parametric Model Constructed from Radiation-Hydrodynamic Simulations

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization-redshift field. The parametric model has three free parameters which can be reduced to one free parameter when we fit the two bias parameters to simulations results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization-redshift field in our model is directly compared cell by cell against simulations and preforms well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (>~ 2 Gpc/h) in order to make mock observations and theoretical predictions.

preprint2012arXiv

Reionization on Large Scales III: Predictions for Low-ell Cosmic Microwave Background Polarization and High-ell Kinetic Sunyaev-Zel'dovich Observables

We present new predictions for temperature (on small angular scales) and polarization (on large angular scales) CMB anisotropies induced during the epoch of reionization (EoR). Using a novel method calibrated from Radiation-Hydrodynamic simulations we model the EoR in large volumes (L >~ 2 Gpc/h) in the context of galactic reionization. We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum (patchy kSZ) ranges between ~0.6 - 2.8 muK^2 at ell = 3000, for the parameter space we explored. These patchy kSZ power spectra are calculated from large 15 Deg x 15 Deg maps that are found to be necessary. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. We find that the amplitude of the patchy kSZ power spectrum at ell = 3000 follows simple scalings of D_ell=3000^kSZ propto <z> and D_ell=3000^kSZ propto Delz^0.47 for the mean redshift (<z>) of reionization and the duration (dz). Using the constraints on <z> from WMAP 7-year results and the lower limit on dz from EDGES we find a lower limit of ~ 0.4 muK^2 on the kSZ at ell = 3000. Planck will constrain the mean redshift and the Thomson optical depth from the low-ell polarization power spectrum. Future measurements of the high-ell CMB power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and polarization measurements constrains both <z> and dz. The patchy kSZ maps, power spectra templates and the polarization power spectra will be publicly available.

preprint2012arXiv

Subhaloes in Self-Interacting Galactic Dark Matter Haloes

We present N-body simulations of a new class of self-interacting dark matter models, which do not violate any astrophysical constraints due to a non-power-law velocity dependence of the transfer cross section which is motivated by a Yukawa-like new gauge boson interaction. Specifically, we focus on the formation of a Milky Way-like dark matter halo taken from the Aquarius project and re-simulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (~1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes are altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dSphs in such a velocity-dependent self-interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.

preprint2012arXiv

The Optimal Cosmic Epoch for Precision Cosmology

The statistical uncertainty in measuring the primordial density perturbations on a given comoving scale is dictated by the number of independent regions of that scale that are accessible to an observer. This number varies with cosmic time and diminishes per Hubble volume in the distant past or future of the standard cosmological model. We show that the best constraints on the initial power spectrum of linear density perturbations are accessible (e.g. through 21-cm intensity mapping) at redshifts z~10, and that the ability to constrain the cosmological initial conditions will deteriorate quickly in our cosmic future.

preprint2012arXiv

The Origin of the Microlensing Events Observed Towards the LMC and the Stellar Counterpart of the Magellanic Stream

We introduce a novel theoretical model to explain the long-standing puzzle of the nature of the microlensing events reported towards the Large Magellanic Cloud (LMC) by the MACHO and OGLE collaborations. We propose that a population of tidally stripped stars from the Small Magellanic Cloud (SMC) located ~4-10 kpc behind a lensing population of LMC disk stars can naturally explain the observed event durations (17-71 days), event frequencies and spatial distribution of the reported events. Differences in the event frequencies reported by the OGLE (~0.33 /yr) and MACHO (~1.75 /yr) surveys appear to be naturally accounted for by their different detection efficiencies and sensitivity to faint sources. The presented models of the Magellanic System were constructed without prior consideration of the microlensing implications. These results favor a scenario for the interaction history of the Magellanic Clouds, wherein the Clouds are on their first infall towards the Milky Way and the SMC has recently collided with the LMC 100-300 Myr ago, leading to a large number of faint sources distributed non-uniformly behind the LMC disk. In contrast to self-lensing models, microlensing events are also expected to occur in fields off the LMC's stellar bar since the stellar debris is not expected to be concentrated in the bar region. This scenario leads to a number of observational tests: the sources are low-metallicity SMC stars, they exhibit high velocities relative to LMC disk stars that may be detectable via proper motion studies, and, most notably, there should exist a stellar counterpart to the gaseous Magellanic Stream and Bridge with a V-band surface brightness of > 32 mag/arcsec^2. In particular, the stellar Bridge should contain enough RR Lyrae stars to be detected by the ongoing OGLE survey of this region.

preprint2012arXiv

Tidal disruption flares from stars on eccentric orbits

We study tidal disruption and subsequent mass fallback for stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that the mass fallback rate decays with the expected -5/3 power of time for parabolic orbits, albeit with a slight deviation due to the self-gravity of the stellar debris. For eccentric orbits, however, there is a critical value of the orbital eccentricity, significantly below which all of the stellar debris is bound to the supermassive black hole. All the mass therefore falls back to the supermassive black hole in a much shorter time than in the standard, parabolic case. The resultant mass fallback rate considerably exceeds the Eddington accretion rate and substantially differs from the -5/3 power of time.

preprint2012arXiv

Unusually Large Fluctuations in the Statistics of Galaxy Formation at High Redshift

We show that various milestones of high-redshift galaxy formation, such as the formation of the first stars or the complete reionization of the intergalactic medium, occurred at different times in different regions of the universe. The predicted spread in redshift, caused by large-scale fluctuations in the number density of galaxies, is at least an order of magnitude larger than previous expectations that argued for a sharp end to reionization. This cosmic scatter in the abundance of galaxies introduces new features that affect the nature of reionization and the expectations for future probes of reionization, and may help explain the present properties of dwarf galaxies in different environments. The predictions can be tested by future numerical simulations and may be verified by upcoming observations. Current simulations, limited to relatively small volumes and periodic boundary conditions, largely omit cosmic scatter and its consequences. In particular, they artificially produce a sudden end to reionization, and they underestimate the number of galaxies by up to an order of magnitude at redshift 20.

preprint2011arXiv

Constraining the Stellar Mass Function in the Galactic Center via Mass Loss from Stellar Collisions

The dense concentration of stars and high velocity dispersions in the Galactic centre imply that stellar collisions frequently occur. Stellar collisions could therefore result in significant mass loss rates. We calculate the amount of stellar mass lost due to indirect and direct stellar collisions and find its dependence on the present-day mass function of stars. We find that the total mass loss rate in the Galactic centre due to stellar collisions is sensitive to the present-day mass function adopted. We use the observed x-ray luminosity in the Galactic centre to preclude any present-day mass functions that result in mass loss rates $> 10^{-5} \mathrm{M_{\odot} yr^{-1}}$ in the vicinity of $\sim 1"$. For present-day mass functions of the form, $dN/dM \propto M^{-α}$, we constrain the present-day mass function to have a minimum stellar mass $\lesssim 7\mathrm{M_{\odot}}$ and a power law slope $\gtrsim 1.25$. We also use this result to constrain the initial mass function in the Galactic centre by considering different star formation scenarios.

preprint2011arXiv

Constraining the Structure of Sagittarius A*'s Accretion Flow with Millimeter-VLBI Closure Phases

Millimeter wave Very Long Baseline Interferometry (mm-VLBI) provides access to the emission region surrounding Sagittarius A*, the supermassive black hole at the center of the Milky Way, on sub-horizon scales. Recently, a closure phase of 0+-40 degrees was reported on a triangle of Earth-sized baselines (SMT-CARMA-JCMT) representing a new constraint upon the structure and orientation of the emission region, independent from those provided by the previously measured 1.3mm-VLBI visibility amplitudes alone. Here, we compare this to the closure phases associated with a class of physically motivated, radiatively inefficient accretion flow models, and present predictions for future mm-VLBI experiments with the developing Event Horizon Telescope (EHT). We find that the accretion flow models are capable of producing a wide variety of closure phases on the SMT-CARMA-JCMT triangle, and thus not all models are consistent with the recent observations. However, those models that reproduce the 1.3mm-VLBI visibility amplitudes overwhelmingly have SMT-CARMA-JCMT closure phases between +-30 degrees, and are therefore broadly consistent with all current mm-VLBI observations. Improving station sensitivity by factors of a few, achievable by increases in bandwidth and phasing together multiple antennas at individual sites, should result in physically relevant additional constraints upon the model parameters and eliminate the current 180 degree ambiguity on the source orientation. When additional stations are included, closure phases of order 45--90 degrees are typical. In all cases the EHT will be able to measure these with sufficient precision to produce dramatic improvements in the constraints upon the spin of Sgr A*.

preprint2011arXiv

Constraints on the Redshift Evolution of the L_X-SFR Relation from the Cosmic X-Ray Backgrounds

Observations of local star forming galaxies have revealed a correlation between the rate at which galaxies form stars and their X-Ray luminosity. We combine this correlation with the most recent observational constraints on the integrated star formation rate density, and find that star forming galaxies account for 5-20% of the total soft and hard X-ray backgrounds, where the precise number depends on the energy band and the assumed average X-ray spectral energy distribution of the galaxies below ~20 keV. If we combine the L_X-SFR relation with recently derived star formation rate function, then we find that star forming galaxies whose X-ray flux falls well (more than a factor of 10) below the detection thresholds of the Chandra Deep Fields, can fully account for the unresolved soft X-ray background, which corresponds to ~6% of its total. Motivated by this result, we put limits on the allowed redshift evolution of the parameter c_X \equiv L_X/SFR, and/or its evolution towards lower and higher star formation rates. If we parametrize the redshift evolution of c_X ~ (1+z)^b, then we find that b \leq 1.3 (95% CL). On the other hand, the observed X-ray luminosity functions (XLFs) of star forming galaxies indicate that c_X may be increasing towards higher redshifts and/or higher star formation rates at levels that are consistent with the X-ray background, but possibly at odds with the locally observed L_X-SFR relation.

preprint2011arXiv

Cosmology with Hypervelocity Stars

In the standard cosmological model, the merger remnant of the Milky Way and Andromeda (Milkomeda) will be the only galaxy remaining within our event horizon once the Universe has aged by another factor of ten, ~10^{11} years after the Big Bang. After that time, the only extragalactic sources of light in the observable cosmic volume will be hypervelocity stars being ejected continuously from Milkomeda. Spectroscopic detection of the velocity-distance relation or the evolution in the Doppler shifts of these stars will allow a precise measurement of the vacuum mass density as well as the local matter distribution. Already in the near future, the next generation of large telescopes will allow photometric detection of individual stars out to the edge of the Local Group, and may target the ~10^{5+-1} hypervelocity stars that originated in it as cosmological tracers.

preprint2011arXiv

Demonstrating the Feasibility of Line Intensity Mapping Using Mock Data of Galaxy Clustering from Simulations

Visbal & Loeb (2010) have shown that it is possible to measure the clustering of galaxies by cross correlating the cumulative emission from two different spectral lines which originate at the same redshift. Through this cross correlation, one can study galaxies which are too faint to be individually resolved. This technique, known as intensity mapping, is a promising probe of the global properties of high redshift galaxies. Here, we test the feasibility of such measurements with synthetic data generated from cosmological dark matter simulations. We use a simple prescription for associating galaxies with dark matter halos and create a realization of emitted radiation as a function of angular position and wavelength over a patch of the sky. This is then used to create synthetic data for two different hypothetical instruments, one aboard the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and another consisting of a pair of ground based radio telescopes designed to measure the CO(1-0) and CO(2-1) emission lines. We find that the line cross power spectrum can be measured accurately from the synthetic data with errors consistent with the analytical prediction of Visbal & Loeb (2010). Removal of astronomical backgrounds and masking bright line emission from foreground contaminating galaxies do not prevent accurate cross power spectrum measurements.

preprint2011arXiv

Effect of Streaming Motion of Baryons Relative to Dark Matter on the Formation of the First Stars

We evaluate the effect of a supersonic relative velocity between the baryons and dark matter on the thermal and density evolution of the first gas clouds at z < 50. Through a series of cosmological simulations, initialized at z=100 with a range of relative streaming velocities and minihalo formation redshifts, we find that the typical streaming velocities will have little effect on the gas evolution. Once the collapse begins, the subsequent evolution of the gas will be nearly indistinguishable from the case of no streaming, and star formation will still proceed in the same way, with no change in the characteristic Pop III stellar masses. Reionization is expected to be dominated by halo masses of > 10^8 M_sun, for which the expected effect of streaming is negligible.

preprint2011arXiv

Effects of supermassive binary black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

preprint2011arXiv

Empirical Constraints on the Coevolution of Supermassive Black Holes and their Host Spheroids

We investigate the evolution of the MBH-σ relation by examining the relationship between the intrinsic scatter in the MBH-σ relation and galaxy bolometric nuclear luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our sample is composed of galaxies with classical bulges when possible, of which 38 have dynamically measured BHs masses, and 17 have BHs masses measured by reverberation mapping. In order to obtain the bolometric nuclear luminosity for galaxies with low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra to bolometric luminosity. We find that the scatter in the MBH-σ relation is uncorrelated with nuclear luminosity over seven orders of magnitude in luminosity, with the high luminosity end approaching the Eddington luminosity. This suggests that at the present epoch galaxies evolve along the MBH-σ relation. This conclusion is consistent with the standard paradigm that BHs grow contemporaneously with their host stellar spheroids.

preprint2011arXiv

Evidence Against Dark Matter Halos Surrounding the Globular Clusters MGC1 and NGC 2419

The conjecture that the ancient globular clusters (GCs) formed at the center of their own dark matter halos was first proposed by Peebles (1984), and has recently been revived to explain the puzzling abundance patterns observed within many GCs. In this paper we demonstrate that the outer stellar density profile of isolated GCs is very sensitive to the presence of an extended dark halo. The GCs NGC 2419, located at 90 kpc from the center of our Galaxy, and MGC1, located at ~200 kpc from the center of M31, are ideal laboratories for testing the scenario that GCs formed at the centers of massive dark halos. Comparing analytic models to observations of these GCs, we conclude that these GCs cannot be embedded within dark halos with a virial mass greater than 10^6 Msun, or, equivalently, the dark matter halo mass-to-stellar mass ratio must be Mdm/M_*<1. If these GCs have indeed orbited within weak tidal fields throughout their lifetimes, then these limits imply that these GCs did not form within their own dark halos. Recent observations of an extended stellar halo in the GC NGC 1851 are also interpreted in the context of our analytic models. Implications of these results for the formation of GCs are briefly discussed.

preprint2011arXiv

Formation of galactic nuclei with multiple supermassive black holes at high redshifts

We examine the formation of groups of multiple supermassive black holes (SMBHs) in gas-poor galactic nuclei due to the high merger rate of galaxies at high redshifts. We calculate the relative likelihood of binary, triple, and quadruple SMBH systems, by considering the timescales for relevant processes and combining merger trees with N-body simulations for the dynamics of stars and SMBHs in galactic nuclei. Typical haloes today with mass $M_0\approx 10^{14}$ M$_\odot$ have an average mass $M_{z=6}=5\times 10^{11}$ M$_\odot$ at $z\sim 6$, while rare haloes with current mass $M_0\gtrsim 10^{15}$ M$_\odot$ have an average mass $M_{z=6}=5\times 10^{12}$ M$_\odot$ at that redshift. These cluster-size haloes are expected to host single galaxies at $z\sim 6$. We expect about 30% galaxies within haloes with present-day mass $M_0\approx 10^{14}$ M$_\odot$ to contain more than two SMBHs at redshifts $2\lesssim z\lesssim 6$. For larger present-day haloes, with $M_0\gtrsim 10^{15}$ M$_\odot$, this fraction is almost 60%. The existence of multiple SMBHs at high redshifts can potentially explain the mass deficiencies observed in the cores of massive elliptical galaxies, which are up to 5 times the mass of their central BHs. Multiple SMBHs would also lead to an enhanced rate of tidal disruption of stars, modified gravitational wave signals compared to isolated BH binaries, and slingshot ejection of SMBHs from galaxies at high speeds in excess of 2000 km s$^{-1}$.

preprint2011arXiv

Gauging the Contribution of X-ray Sources to Reionization Through the Kinetic Sunyaev-Zel'dovich Effect

Measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect from instruments such as the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) will soon put improved constraints on reionization. Popular models assume that UV photons alone are responsible for reionization of the intergalactic medium. We explore the effects of a significant contribution of X-rays to reionization on the kSZ signal. Because X-rays have a large mean free path through the neutral intergalactic medium, they introduce partial ionization in between the sharp-edged bubbles created by UV photons. This smooth ionization component changes the power spectrum of the cosmic microwave background (CMB) temperature anisotropies. We quantify this effect by running semi-numerical simulations of reionization. We test a number of different models of reionization without X-rays that have varying physical parameters, but which are constrained to have similar total optical depths to electron scattering. These are then compared to models with varying levels of contribution to reionization from X-rays. We find that models with more than a 10% contribution from X-rays produce a significantly lower power spectrum of temperature anisotropies than all the UV-only models tested. The expected sensitivity of SPT and ACT may be insufficient to distinguish between our models, however, a non-detection of the kSZ signal from the epoch of reionization could result from the contribution of X-rays. It will be important for future missions with improved sensitivity to consider the impact of X-ray sources on reionization.

preprint2011arXiv

Imprint of Accretion Disk-Induced Migration on Gravitational Waves from Extreme Mass Ratio Inspirals

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

preprint2011arXiv

Inflation and the Scale Dependent Spectral Index: Prospects and Strategies

We consider the running of the spectral index as a probe of both inflation itself, and of the overall evolution of the very early universe. Surveying a collection of simple single field inflationary models, we confirm that the magnitude of the running is relatively consistent, unlike the tensor amplitude, which varies by orders of magnitude. Given this target, we confirm that the running is potentially detectable by future large scale structure or 21 cm observations, but that only the most futuristic measurements can distinguish between these models on the basis of their running. For any specified inflationary scenario, the combination of the running index and unknown post-inflationary expansion history induces a theoretical uncertainty in the predicted value of the spectral index. This effect can easily dominate the statistical uncertainty with which Planck and its successors are expected to measure the spectral index. More positively, upcoming cosmological experiments thus provide an intriguing probe of physics between TeV and GUT scales by constraining the reheating history associated with any specified inflationary model, opening a window into the "primordial dark age" that follows the end of inflation.

preprint2011arXiv

Localizing Sagittarius A* and M87 on Microarcsecond Scales with Millimeter VLBI

With the advent of the Event Horizon Telescope (EHT), a millimeter/sub-millimeter very-long baseline interferometer (VLBI), it has become possible to image a handful of black holes with sub-horizon resolutions. However, these images do not translate into microarcsecond absolute positions due to the lack of absolute phase information when an external phase reference is not used. Due to the short atmospheric coherence time at these wavelengths, nodding between the source and phase reference is impractical. However, here we suggest an alternative scheme which makes use of the fact that many of the VLBI stations within the EHT are arrays in their own right. With this we show that it should be possible to absolutely position the supermassive black holes at the centers of the Milky Way (Sgr A*) and M87 relative to nearby objects with precisions of roughly 1 microarcsecond. This is sufficient to detect the perturbations to Sgr A*'s position resulting from interactions with the stars and stellar-mass black holes in the Galactic cusp on year timescales, and severely constrain the astrophysically relevant parameter space for an orbiting intermediate mass black hole, implicated in some mechanisms for producing the young massive stars in the Galactic center. For M87, it allows the registering of millimeter images, in which the black hole may be identified by its silhouette against nearby emission, and existing larger scale radio images, eliminating present ambiguities in the nature of the radio core and inclination, opening angle, and source of the radio jet.

preprint2011arXiv

Measuring the Redshift of Reionization with a Modest Array of Low-Frequency Dipoles

The designs of the first generation of cosmological 21-cm observatories are split between single dipole experiments which integrate over a large patch of sky in order to find the global (spectral) signature of reionization, and interferometers with arcminute-scale angular resolution whose goal is to measure the 3D power spectrum of ionized regions during reionization. We examine whether intermediate scale instruments with complete Fourier (uv) coverage are capable of placing new constraints on reionization. We find that even without using a full power spectrum analysis, the global redshift of reionization, z_reion, can in principle be measured from the variance in the 21-cm signal among multiple beams as a function of frequency at a roughly 1 degree angular scale. At this scale, the beam-to-beam variance in the differential brightness temperature peaks when the average neutral fraction was around 50%, providing a convenient flag of z_reion. We choose a low angular resolution of order 1 degree to exploit the physical size of the ionized regions and maximize the signal-to-noise ratio. Thermal noise, foregrounds, and instrumental effects should also be manageable at this angular scale, as long as the uv coverage is complete within the compact core required for low-resolution imaging. For example, we find that z_reion can potentially be detected to within a redshift uncertainty of less than around 1 in around 500 hours of integration on the existing MWA prototype (with only 32x16 dipoles), operating at an angular resolution of around 1 degree and a spectral resolution of 2.4 MHz.

preprint2011arXiv

Observable Signatures of EMRI Black Hole Binaries Embedded in Thin Accretion Disks

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

preprint2011arXiv

Observational selection effects and the M-sigma relation

We examine the possibility that the observed relation between black-hole mass and host-galaxy stellar velocity dispersion (the M-sigma relation) is biased by an observational selection effect, the difficulty of detecting a black hole whose sphere of influence is smaller than the telescope resolution. In particular, we critically investigate recent claims that the M-sigma relation only represents the upper limit to a broad distribution of black-hole masses in galaxies of a given velocity dispersion. We find that this hypothesis can be rejected at a high confidence level, at least for the early-type galaxies with relatively high velocity dispersions (median 268 km/s) that comprise most of our sample. We also describe a general procedure for incorporating observational selection effects in estimates of the properties of the M-sigma relation. Applying this procedure we find results that are consistent with earlier estimates that did not account for selection effects, although with larger error bars. In particular, (i) the width of the M-sigma relation is not significantly increased; (ii) the slope and normalization of the M-sigma relation are not significantly changed; (iii) most or all luminous early-type galaxies contain central black holes at zero redshift. Our results may not apply to late-type or small galaxies, which are not well-represented in our sample.

preprint2011arXiv

On The GeV & TeV Detections of the Starburst Galaxies M82 & NGC 253

The GeV and TeV emission from M82 and NGC 253 observed by Fermi, HESS, and VERITAS constrains the physics of cosmic rays (CRs) in these dense starbursts. We argue that the gamma rays are predominantly hadronic in origin, as expected by previous studies. The measured fluxes imply that pionic losses are efficient for CR protons in both galaxies: we show that a fraction F_cal ~ 0.2 - 0.4 of the energy injected in high energy primary CR protons is lost to inelastic proton-proton collisions (pion production) before escape, producing gamma rays, neutrinos, and secondary electrons and positrons. We discuss the factor ~2 uncertainties in this estimate, including supernova rate and leptonic contributions to the GeV-TeV emission. We argue that gamma-ray data on ULIRGs like Arp 220 can test whether M82 and NGC 253 are truly calorimetric, and we present upper limits on Arp 220 from the Fermi data. We show that the observed ratio of the GeV to GHz fluxes of the starbursts suggests that non-synchrotron cooling processes are important for cooling the CR electron/positron population. We briefly reconsider previous predictions in light of the gamma-ray detections, including the starburst contribution to the gamma-ray background and CR energy densities. Finally, as a guide for future studies, we list the brightest star-forming galaxies on the sky and present updated predictions for their gamma-ray and neutrino fluxes.

preprint2011arXiv

Photon Trapping Enables Super-Eddington Growth of Black-Hole Seeds in Galaxies at High Redshift

We identify a physical mechanism that would have resulted in rapid, obscured growth of seed super-massive black-holes in galaxies at z>6. Specifically, we find that the density at the centre of typical high redshift galaxies was at a level where the Bondi accretion rate implies a diffusion speed of photons that was slower than the gravitational infall velocity, resulting in photons being trapped within the accretion flow and advected into the black-hole. We show that there is a range of black-hole masses (M_bh ~ 10^3-10^5 solar masses) where the accretion flow traps radiation, corresponding to black-holes that were massive enough to generate a photon trapping accretion flow, but small enough that their Bondi radii did not exceed the isothermal scale height of self-gravitating gas. Under these conditions we find that the accretion reaches levels far in excess of the Eddington rate. A prediction of this scenario is that X-ray number counts of active galactic nuclei at z>6 would exhibit a cutoff at the low luminosities corresponding to black-hole masses below ~10^5 solar masses. At low redshifts we find photon trapping to be unimportant because it could only occur in rare low spin halos, and would require black-hole masses in excess of expectations from the observed black-hole - halo mass relation. The super-Eddington growth of ~10^5 solar mass seed black-holes at high redshift may have provided a natural acceleration towards the growth of super-massive black-holes at z~6-7, less than a billion years after the Big Bang.

preprint2011arXiv

Pulsar-Black Hole Binaries in the Galactic Center

Binaries consisting of a pulsar and a black hole (BH) are a holy grail of astrophysics, both for their significance for stellar evolution and for their potential application as probes of strong gravity. In spite of extensive surveys of our Galaxy and its system of globular clusters, no pulsar-black hole (PSR-BH) binary has been found to date. Clues as to where such systems might exist are therefore highly desirable. We show that if the central parsec around Sgr A* harbors a cluster of ~25,000 stellar BHs (as predicted by mass segregation arguments) and if it is also rich in recycled pulsar binaries (by analogy with globular clusters), then 3-body exchange interactions should produce PSR-BHs in the Galactic center. Simple estimates of the formation rate and survival time of these binaries suggest that a few PSR-BHs should be present in the central parsec today. The proposed formation mechanism makes unique predictions for the PSR-BH properties: 1) the binary would reside within ~1 pc of Sgr A*; 2) the pulsar would be recycled, with a period of ~1 to a few tens of milliseconds, and a low magnetic field B<~10^10 G; 3) the binary would have high eccentricity, e~0.8, but with a large scatter; and 4) the binary would be relatively wide, with semi-major axis a_b~0.1 - >~3 AU. The potential discovery of a PSR-BH binary therefore provides a strong motivation for deep, high-frequency radio searches for recycled pulsars toward the Galactic center.

preprint2011arXiv

Rating Growth of Scientific Knowledge and Risk from Theory Bubbles

In physics the value of a theory is measured by its agreement with experimental data. But how should the physics community gauge the value of an emerging theory that has not been tested experimentally as of yet? With no reality check, a hypothesis like string theory may linger for a while before physicists will know its actual value in describing nature. In this short article, I advocate the need for a website operated by graduate students that will use various measures of publicly available data (such as the growth rate of newly funded experiments, research grants, publications, and faculty jobs) to gauge the future dividends of various research frontiers. The analysis can benefit from past experience (e.g. in research areas that suffered from limited experimental data over long periods of time) and aim to alert the community of the risk from future theory bubbles.

preprint2011arXiv

Recoiled star clusters in the Milky Way halo: N-body simulations and a candidate search through SDSS

During the formation of the Milky Way, > 100 central black holes (BHs) may have been ejected from their small host galaxies as a result of asymmetric gravitational wave emission. We previously showed that many of these BHs are surrounded by a compact cluster of stars that remained bound to the BH during the ejection process. In this paper, we perform long term N-body simulations of these star clusters to determine the distribution of stars in these clusters today. These numerical simulations, reconciled with our Fokker-Planck simulations, show that stellar density profile follows a power-law with slope ~ -2.15, and show that large angle scattering and tidal disruptions remove 20 - 90% of the stars by ~10^10 yr. We then analyze the photometric and spectroscopic properties of recoiled clusters accounting for the small number of stars in the clusters. We use our results to perform a systematic search for candidates in the Sloan Digital Sky Survey. We find no spectroscopic candidates, in agreement with our expectations from the completeness of the survey. Using generic photometric models of present day clusters we identify ~100 recoiling cluster candidates. Follow-up spectroscopy would be able to determine the nature of these candidates.

preprint2011arXiv

Recoiling Black Holes in Merging Galaxies

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed (vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to about 0.7 vesc may be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will depend on the timing of the BH merger relative to the change in vesc. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (v > 500 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of up to about 10 - 100 Myr. (5) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher central star formation rates in the merger remnant, thereby extending the starburst phase of the merger and creating a denser, more massive stellar cusp.

preprint2011arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes and Merger Remnant Properties

Central supermassive black holes (SMBHs) are a ubiquitous feature of locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and their host galaxies is closely linked. However, in the event of a merger, gravitational-wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. To explore the consequences of this phenomenon, we use hydrodynamic simulations of gaseous galaxy mergers that include a range of BH recoil velocities. We have generated a suite of over 200 simulations with more than 60 merger models, enabling us to identify systematic trends in the behavior of recoiling BHs -- specifically (i) their dynamics, (ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution. (i) Recoiling BH trajectories depend heavily on the gas content of the host galaxy; maximal BH displacements from the center may vary by up to an order of magnitude between gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the timing of the BH merger relative to the formation of the galaxy merger remnant. (ii) Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kinematic offsets (v > 800 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of about 1 - 100 Myr. In addition, recoil events affect the total AGN lifetime. GW recoil generally reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. (iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. Furthermore, recoil events displace AGN feedback from the galactic center, which enhances central star formation rates. This results in longer starburst phases and higher central stellar densities in merger remnants.

preprint2011arXiv

Supernova Shock Breakout Through a Wind

The breakout of a supernova shock wave through the progenitor star's outer envelope is expected to appear as an X-ray flash. However, if the supernova explodes inside an optically-thick wind, the breakout flash is delayed. We present a simple model for estimating the conditions at shock breakout in a wind based on the general observable quantities in the X-ray flash lightcurve: the total energy E_X, and the diffusion time after the peak, t_diff. We base the derivation on the self-similar solution for the forward-reverse shock structure expected for an ejecta plowing through a pre-existing wind at large distances from the progenitor's surface. We find simple quantitative relations for the shock radius and velocity at breakout. By relating the ejecta density profile to the pre-explosion structure of the progenitor, the model can also be extended to constrain the combination of explosion energy and ejecta mass. For the observed case of XRO08109/SN2008D, our model provides reasonable constraints on the breakout radius, explosion energy, and ejecta mass, and predicts a high shock velocity which naturally accounts for the observed non-thermal spectrum.

preprint2011arXiv

The cosmic X-ray and gamma-ray background from dark matter annihilation

(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-rays. We here create full-sky maps of the radiation from both of these contributions using the high-resolution Millennium-II simulation. We use upper limits on the contributions of unknown sources to the EBL to constrain the intrinsic properties of dark matter using a model-independent approach that can be employed as a template to test different particle physics models (including those with a Sommerfeld enhancement). These upper limits are based on observations spanning eight orders of magnitude in energy (from soft X-rays measured by CHANDRA to gamma-rays measured by Fermi), and on expectations for the contributions from blazars and star forming galaxies. To exemplify this approach, we analyze a set of benchmark Sommerfeld-enhanced models that give the correct dark matter abundance, satisfy CMB constraints, and fit the cosmic ray spectra measured by PAMELA and Fermi without any contribution from local subhalos. We find that these models are in conflict with the EBL constraints unless the contribution of unresolved subhalos is small and the annihilation signal dominates the EBL. We conclude that provided the collisionless cold dark matter paradigm is accurate, even for conservative estimates of the contribution from unresolved subhalos and astrophysical backgrounds, the EBL is at least as sensitive a probe of these types of scenarios as the CMB. Our results disfavor an explanation of the positron excess measured by PAMELA based only on dark matter annihilation in the smooth Galactic halo.

preprint2011arXiv

The imprint of the relative velocity between baryons and dark matter on the 21-cm signal from reionization

The post-recombination streaming of baryons through dark matter keeps baryons out of low mass (<10^6 solar masses) halos coherently on scales of a few comoving Mpc. It has been argued that this will have a large impact on the 21-cm signal before and after reionization, as it raises the minimum mass required to form ionizing sources. Using a semi-numerical code, we show that the impact of the baryon streaming effect on the 21-cm signal during reionization (redshifts z approximately 7-20) depends strongly on the cooling scenario assumed for star formation, and the corresponding virial temperature or mass at which stars form. For the canonical case of atomic hydrogen cooling at 10^4 Kelvin, the minimum mass for star formation is well above the mass of halos that are affected by the baryon streaming and there are no major changes to existing predictions. For the case of molecular hydrogen cooling, we find that reionization is delayed by a change in redshift of approximately 2 and that more relative power is found in large modes at a given ionization fraction. However, the delay in reionization is degenerate with astrophysical assumptions, such as the production rate of UV photons by stars.

preprint2011arXiv

The Shocking Truth: The small contribution to hydrogen reionization from gravitational infall

It is commonly thought that stars are responsible for reionizing the Universe. However, deep searches for star-forming galaxies during the epoch of reionization have not yet found sufficient galaxies to provide the necessary ionising flux. Fast accretion shocks associated with gravitational infall of baryons during the formation of galaxies have recently been proposed as an alternative method of generating the required ionising photons. In this Letter we assess the contribution to hydrogen reionization from shocked gas associated with gravitational infall. We find that shocks can ionize no more than a few percent of the cosmic hydrogen by z~6. However, the small fraction of ionizing radiation produced by fast accretion shocks would be significantly more biased than that associated with stars, leading to a modification of the luminosity weighted source clustering by ~10%. This modification of the bias may be measurable with future precision experiments utilising the redshifted 21cm line to study the distribution of hydrogen during the reionization era.

preprint2011arXiv

Tidal Disruption Flares of Stars From Moderately Recoiled Black Holes

We analyze stellar tidal disruption events as a possible observational signature of gravitational wave induced recoil of supermassive black holes. As a black hole wanders through its galaxy, it will tidally disrupt bound and unbound stars at rates potentially observable by upcoming optical transient surveys. To quantify these rates, we explore a broad range of host galaxy and black hole kick parameters. We find that emission from a transient accretion disk can produce ~1 event per year which LSST would identify as spatially offset, while super-Eddington tidal flares, if they exist, are likely to produce ~10 spatially offset events per year. A majority of tidal disruption flares, and a large majority of flares with an observable spatial offset, are due to bound rather than unbound stars. The total number of disruption events due to recoiled black holes could be almost 1% of the total stellar tidal disruption rate.

preprint2010arXiv

Constraining Parity Violation in Gravity with Measurements of Neutron-Star Moments of Inertia

Neutron stars are sensitive laboratories for testing general relativity, especially when considering deviations where velocities are relativistic and gravitational fields are strong. One such deviation is described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modified through the addition of the gravitational parity-violating Pontryagin density coupled to a field. This four-dimensional effective theory arises naturally both in perturbative and non-perturbative string theory, loop quantum gravity, and generic effective field theory expansions. We calculate here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neutron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of the metric to leading order, thus introducing modifications to the moment of inertia but not to the mass-radius relation. We show that an observational determination of the moment of inertia to an accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a constraint on the Chern-Simons coupling constant of ξ^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

preprint2010arXiv

Constraining the Minimum Mass of High-Redshift Galaxies and Their Contribution to the Ionization State of the IGM

We model the latest HST WFPC3/IR observations of > 100 galaxies at redshifts z=7-8 in terms of a hierarchical galaxy formation model with starburst activity. Our model provides a distribution of UV luminosities per dark matter halo of a given mass and a natural explanation for the fraction of halos hosting galaxies. The observed luminosity function is best fit with a minimum halo mass per galaxy of 10^{9.4+0.3-0.9} Msun, corresponding to a virial temperature of 10^{4.9+0.2-0.7} K. Extrapolating to faint, undetected galaxies, the total production rate of ionizing radiation depends critically on this minimum mass. Future measurements with JWST should determine whether the entire galaxy population can comfortably account for the UV background required to keep the intergalactic medium ionized.

preprint2010arXiv

Constraining the unexplored period between reionization and the dark ages with observations of the global 21 cm signal

Observations of the frequency dependence of the global brightness temperature of the redshifted 21 cm line of neutral hydrogen may be possible with single dipole experiments. In this paper, we develop a Fisher matrix formalism for calculating the sensitivity of such instruments to the 21 cm signal from reionization and the dark ages. We show that rapid reionization histories with duration delta z< 2 can be constrained, provided that local foregrounds can be well modelled by low order polynomials. It is then shown that observations in the range nu = 50 - 100 MHz can feasibly constrain the Lyman alpha and X-ray emissivity of the first stars forming at z = 15 - 25, provided that systematic temperature residuals can be controlled to less than 1 mK. Finally, we demonstrate the difficulty of detecting the 21 cm signal from the dark ages before star formation.

preprint2010arXiv

Cores in Dwarf Galaxies from Dark Matter with a Yukawa Potential

We show that cold dark matter particles interacting through a Yukawa potential could naturally explain the recently observed cores in dwarf galaxies without affecting the dynamics of objects with a much larger velocity dispersion, such as clusters of galaxies. The velocity dependence of the associated cross-section as well as the possible exothermic nature of the interaction alleviates earlier concerns about strongly interacting dark matter. Dark matter evaporation in low-mass objects might explain the observed deficit of satellite galaxies in the Milky Way halo and have important implications for the first galaxies and reionization.

preprint2010arXiv

Dynamics and Magnetization in Galaxy Cluster Cores Traced by X-ray Cold Fronts

Cold fronts (CFs) - density and temperature plasma discontinuities - are ubiquitous in cool cores of galaxy clusters, where they appear as X-ray brightness edges in the intracluster medium, nearly concentric with the cluster center. We analyze the thermodynamic profiles deprojected across core CFs found in the literature. While the pressure appears continuous across these CFs, we find that all of them require significant centripetal acceleration beneath the front. This is naturally explained by a tangential, nearly sonic bulk flow just below the CF, and a tangential shear flow involving a fair fraction of the plasma beneath the front. Such shear should generate near-equipartition magnetic fields on scales ~<50 pc from the front, and could magnetize the entire core. Such fields would explain the apparent stability of cool-core CFs and the recently reported CF-radio minihalo association.

preprint2010arXiv

Evidence for Low Black Hole Spin and Physically Motivated Accretion Models from Millimeter VLBI Observations of Sagittarius A*

Millimeter very-long baseline interferometry (mm-VLBI) provides the novel capacity to probe the emission region of a handful of supermassive black holes on sub-horizon scales. For Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, this provides access to the region in the immediate vicinity of the horizon. Broderick et al. (2009) have already shown that by leveraging spectral and polarization information as well as accretion theory, it is possible to extract accretion-model parameters (including black hole spin) from mm-VLBI experiments containing only a handful of telescopes. Here we repeat this analysis with the most recent mm-VLBI data, considering a class of aligned, radiatively inefficient accretion flow (RIAF) models. We find that the combined data set rules out symmetric models for Sgr A*'s flux distribution at the 3.9-sigma level, strongly favoring length-to-width ratios of roughly 2.4:1. More importantly, we find that physically motivated accretion flow models provide a significantly better fit to the mm-VLBI observations than phenomenological models, at the 2.9-sigma level. This implies that not only is mm-VLBI presently capable of distinguishing between potential physical models for Sgr A*'s emission, but further that it is sensitive to the strong gravitational lensing associated with the propagation of photons near the black hole. Based upon this analysis we find that the most probable magnitude, viewing angle, and position angle for the black hole spin are a=0.0(+0.64+0.86), theta=68(+5+9)(-20-28) degrees, and xi=-52(+17+33)(-15-24) east of north, where the errors quoted are the 1-sigma and 2-sigma uncertainties.

preprint2010arXiv

Extrapolating the Evolution of Galaxy Sizes to the Epoch of Reionization

We use data on the high-redshift evolution of the size distribution and luminosity function of galaxies to constrain the relationship between their star formation efficiency and starburst lifetime. Based on the derived scaling relations, we predict the angular sizes and average surface brightnesses of faint galaxies that will be discovered with JWST. We find that JWST will be able to resolve galaxies at the magnitude limit m<31 out to a redshift of z~14. The next generation of large ground-based telescopes will resolve all galaxies discovered with JWST, provided they are sufficiently clumpy to enable detection above the bright thermal sky. We combine our constraints with simple models for self regulation of star formation, and show that feedback from supernovae at redshifts z>3 is likely mediated through momentum transfer, with the starburst timescale set by the lifetime of the massive stars rather than the dynamical time in the host galactic disk.

preprint2010arXiv

Galaxy Statistics in Pencil-beam Surveys at High Redshifts

Surveys of faint galaxies at high redshifts often result in a "pencil-beam" geometry that is much longer along the line-of-sight than across the sky. We explore the effects of this geometry on the abundance and clustering of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) in current and future surveys based on cosmological N-body simulations which adequately describe the nonlinear growth of structure on small scales and compare to linear theory. We find that the probability distribution of the LBG abundance is skewed toward low values since the narrow transverse dimension of the survey is more likely to probe underdense regions. Over a range that spans 1--2 orders of magnitude in galaxy luminosities, the variance in the number of objects differs from the commonly used analytic prediction and is not dominated by Poisson noise. Additionally, nonlinear bias on small scales results in a one-dimensional power spectrum of LAEs using a James Webb Space Telescope field-of-view that is relatively flat, markedly different from the expectation of linear perturbation theory. We discuss how these results may affect attempts to measure the UV background at high redshifts, estimate the relationship between halo mass and galaxy luminosity, and probe reionization by measuring the power-modulating effect of ionized regions.

preprint2010arXiv

Identifying Supermassive Black Hole Binaries with Broad Emission Line Diagnosis

Double-peaked broad emission lines in Active Galactic Nuclei (AGNs) may indicate the existence of a bound supermassive black hole (SMBH) binary where two distinct broad line regions (BLRs) contribute together to the line profile. An alternative interpretation is a disk emitter origin for the double-peaked line profile. Using simple BLR models, we calculate the expected broad line profile for a SMBH binary at different separations. Under reasonable assumptions that both BLRs are illuminated by the two active SMBHs and that the ionizing flux at the BLR location is roughly constant, we confirm the emergence of double-peaked features and radial velocity drifts of the two peaks due to the binary orbital motion. However, such a clear double-peaked feature only arises in a particular stage of the binary evolution when the two BHs are close enough such that the line-of-sight orbital velocity difference is larger than the FWHM of the individual broad components, while the two BLRs are still mostly distinct. Prior to this stage, the velocity splitting due to the orbit motion of the binary is too small to separate the emission from the two BLRs, leading to asymmetric broad line profiles in general. When the two BHs are even closer such that the two BLRs can no longer be distinct, the line profile becomes more complex and the splitting of the peaks does not correspond to the orbital motion of the binary. In this regime there are no coherent radial velocity drifts in the peaks with time. Asymmetric line profiles are probably a far more common signature of binary SMBHs than are double-peaked profiles. We discuss the temporal variations of the broad line profile for binary SMBHs and highlight the different behaviors of reverberation mapping in the binary and disk emitter cases, which may serve as a feasible tool to disentangle these two scenarios.

preprint2010arXiv

Measuring the 3D Clustering of Undetected Galaxies Through Cross Correlation of their Cumulative Flux Fluctuations from Multiple Spectral Lines

We discuss a method for detecting the emission from high redshift galaxies by cross correlating flux fluctuations from multiple spectral lines. If one can fit and subtract away the continuum emission with a smooth function of frequency, the remaining signal contains fluctuations of flux with frequency and angle from line emitting galaxies. Over a particular small range of observed frequencies, these fluctuations will originate from sources corresponding to a series of different redshifts, one for each emission line. It is possible to statistically isolate the fluctuations at a particular redshift by cross correlating emission originating from the same redshift, but in different emission lines. This technique will allow detection of clustering fluctuations from the faintest galaxies which individually cannot be detected, but which contribute substantially to the total signal due to their large numbers. We describe these fluctuations quantitatively through the line cross power spectrum. As an example of a particular application of this technique, we calculate the signal-to-noise ratio for a measurement of the cross power spectrum of the OI(63 micron) and OIII(52 micron) fine structure lines with the proposed Space Infrared Telescope for Cosmology and Astrophysics. We find that the cross power spectrum can be measured beyond a redshift of z=8. Such observations could constrain the evolution of the metallicity, bias, and duty cycle of faint galaxies at high redshifts and may also be sensitive to the reionization history through its effect on the minimum mass of galaxies. As another example, we consider the cross power spectrum of CO line emission measured with a large ground based telescope like CCAT and 21-cm radiation originating from hydrogen in galaxies after reionization with an interferometer similar in scale to MWA, but optimized for post-reionization redshifts.

preprint2010arXiv

Prompt Tidal Disruption of Stars as an Electromagnetic Signature of Supermassive Black Hole Coalescence

A precise electromagnetic measurement of the sky coordinates and redshift of a coalescing black hole binary holds the key for using its gravitational wave (GW) signal to constrain cosmological parameters and to test general relativity. Here we show that the merger of ~10^{6-7}M_sun black holes is generically followed over a period of years by multiple electromagnetic flares from tidally disrupted stars. The sudden recoil imparted to the merged black hole by GW emission promptly fills its loss cone and results in a tidal disruption rate of stars as high as ~0.1 per year. The prompt disruption of a star within a single galaxy over a short period provides a unique electromagnetic flag of a recent black hole coalescence event, and sequential disruptions could be used on their own to calibrate the expected rate of GW sources for pulsar timing arrays or the proposed Laser Interferometer Space Antenna (LISA).

preprint2010arXiv

Recoiling Black Holes in Merging Galaxies: Relationship to AGN Lifetimes, Starbursts, and the M-sigma Relation

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. This enables us to follow recoiling BHs in self-consistent, evolving merger remnants. In contrast to recent studies on similar topics, we conduct a large parameter study, generating a suite of over 200 simulations with more than 60 merger models and a range of recoil velocities (vk). Our main results are as follows. (1) BHs kicked at nearly the central escape speed (vesc) may oscillate on large orbits for up to a Hubble time, but in gas-rich mergers, BHs kicked with up to ~ 0.7 vesc may be confined to the central few kpc of the galaxy, owing to gas drag and steep central potentials. (2) vesc in gas-rich mergers may increase rapidly during final coalescence, in which case trajectories may depend on the timing of the BH merger relative to the formation of the potential well. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN), but may actually extend AGN lifetimes at lower luminosities. (4) Kinematically-offset AGN (v > 800 km s^-1) may be observable for up to ~ 10 Myr either immediately after the recoil or during pericentric passages through a gas-rich remnant. (5) Spatially-offset AGN (R > 1 kpc) generally have low luminosities and lifetimes of ~ 1 - 100 Myr. (6) Rapidly-recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. This lowers the normalization of the M-sigma relation and contributes to both intrinsic and overall scatter. (7) Finally, the displacement of AGN feedback after a recoil event enhances central star formation rates, thereby extending the starburst phase of the merger and creating a denser stellar cusp. [Abridged.]

preprint2010arXiv

Rotation Speed of the First Stars

We estimate the rotation speed of Population III (Pop III) stars within a minihalo at z ~ 20 using a smoothed particle hydrodynamics (SPH) simulation, beginning from cosmological initial conditions. We follow the evolution of the primordial gas up to densities of 10^12 cm^-3. Representing the growing hydrostatic cores with accreting sink particles, we measure the velocities and angular momenta of all particles that fall onto these protostellar regions. This allows us to record the angular momentum of the sinks and estimate the rotational velocity of the Pop III stars expected to form within them. The rotation rate has important implications for the evolution of the star, the fate encountered at the end of its life, and the potential for triggering a gamma-ray burst (GRB). We find that there is sufficient angular momentum to yield rapidly rotating stars (> 1000 km s^-1, or near break-up speeds). This indicates that Pop III stars likely experienced strong rotational mixing, impacting their structure and nucleosynthetic yields. A subset of them was also likely to result in hypernova explosions, and possibly GRBs.

preprint2010arXiv

Scale-Dependent Bias of Galaxies from Baryonic Acoustic Oscillations

Baryonic acoustic oscillations (BAOs) modulate the density ratio of baryons to dark matter across large regions of the Universe. We show that the associated variation in the mass-to-light ratio of galaxies should generate an oscillatory, scale-dependent bias of galaxies relative to the underlying distribution of dark matter. A measurement of this effect would calibrate the dependence of the characteristic mass-to-light ratio of galaxies on the baryon mass fraction in their large scale environment. This bias, though, is unlikely to significantly affect measurements of BAO peak positions.

preprint2010arXiv

Taking "The Road Not Taken'': On the Benefits of Diversifying Your Academic Portfolio

It is common practice among young astrophysicists these days to invest research time conservatively in mainstream ideas that have already been explored extensively in the literature. This tendency is driven by peer pressure and job market prospects, and is occasionally encouraged by senior researchers. Although the same phenomenon existed in past decades, it is alarmingly more prevalent today because a growing fraction of observational and theoretical projects are pursued in large groups with rigid research agendas. In addition, the emergence of a ``standard model'' in cosmology (albeit with unknown dark components) offers secure ``bonds'' for a safe investment of research time. In this short essay, which summarizes a banquet lecture at a recent conference, I give examples for both safe and risky topics in astrophysics (which I split into categories of ``bonds,'' ``stocks,'' and ``venture capital''), and argue that young researchers should always allocate a small fraction of their academic portfolio to innovative projects with risky but potentially highly profitable returns. In parallel, selection and promotion committees must find new strategies for rewarding candidates with creative thinking.

preprint2010arXiv

Ultra High Energy Cosmic Ray Acceleration in Engine-driven Relativistic Supernovae

The origin of the highest energy cosmic rays remains an enigma. They offer a window to new physics, including tests of physical laws relevant to their propagation and interactions, at energies unattainable by terrestrial accelerators. They must be accelerated locally, as otherwise background radiations would severely suppress the flux of protons and nuclei, at energies above the Greisen-Zatsepin-Kuzmin (GZK) limit. Nearby Gamma Ray Bursts (GRBs), Hypernovae, Active Galactic Nuclei (AGNs) and their flares, have all been suggested and debated as possible sources. A local sub-population of type Ibc supernovae (SNe) with mildly relativistic outflows have been detected as sub-energetic GRBs or X-Ray Flashes (XRFs) and recently as radio afterglows without detected GRB counterparts. We measure the size-magnetic field evolution, baryon loading and energetics, using the observed radio spectra of SN 2009bb. We place such engine-driven SNe above the Hillas line and establish that they can readily explain the post-GZK UHECRs.

preprint2010arXiv

Using Radio Halos and Minihalos to Measure the Distributions of Magnetic Fields and Cosmic-Rays in Galaxy Clusters

Some galaxy clusters show diffuse radio emission in the form of giant halos (GHs) on Mpc scales or minihalos (MHs) on smaller scales. Comparing VLA and XMM radial profiles of several such clusters, we find a universal linear correlation between radio and X-ray surface brightness, valid in both types of halos. It implies a halo central emissivity nu*j_nu=10^{-31.4+-0.2}*(n/10^-2cm^-3)^2*(T/T_0)^{0.2+-0.5} erg/s/cm^3, where T and T_0 are the local and central temperatures, and n is the electron number density. We argue that the tight correlation and the scaling of j_nu, combined with morphological and spectral evidence, indicate that both GHs and MHs arise from secondary electrons and positrons, produced in cosmic-ray ion (CRI) collisions with a strongly magnetized, B >3 muG intracluster gas. When the magnetic energy density drops below that of the microwave background, the radio emission weakens considerably, producing halos with a clumpy morphology (e.g. RXCJ2003.5-2323 and A2255) or a distinct radial break. We thus measure a magnetic field B=3 muG at a radius r~110 kpc in A2029 and r~50 kpc in Perseus. The spectrum of secondaries, produced from hadronic collisions of ~20 GeV CRIs, reflects the energy dependence of the collision cross section. We use the observed spectra of halos, in particular where they steepen with increasing radius or frequency, to (i) measure B~10(nu/700 MHz) muG, with nu the spectral break frequency; (ii) identify a correlation between the average spectrum and the central magnetic field; and (iii) infer a CRI spectral index s <-2.7 and energy fraction xi_p~10^{-3.6+-0.2} at particle energies above 10 GeV. Our results favor a model where CRIs diffuse away from their sources (which are probably supernovae, according to a preliminary correlation with star formation), whereas the magnetic fields are generated by mergers in GHs and by core sloshing in MHs.

preprint2009arXiv

Acceleration of Galactic Supershells by Lyman Alpha Radiation

Scattering of Lyman Alpha (hereafter Lya) photons by neutral hydrogen gas in a single outflowing 'supershell' around star forming regions often explains the shape and offset of the observed Lya emission line from galaxies. We compute the radiation pressure that is exerted by this scattered Lya radiation on the outflowing material. We show that for reasonable physical parameters, Lya radiation pressure alone can accelerate supershells to velocities in the range v_sh=200-400 km/s. These supershells possibly escape from the gravitational potential well of their host galaxies and contribute to the enrichment of the intergalactic medium. We compute the physical properties of expanding supershells that are likely to be present in a sample of known high-redshift (z=2.7-5.0) galaxies, under the assumption that they are driven predominantly by Lya radiation pressure. We predict ranges of radii r_sh=0.1-10 kpc, ages t_sh=1-100 Myr, and energies E_sh=1e53-1e55 ergs, which are in reasonable agreement with the properties of local galactic supershells. Furthermore, we find that the radius, r_sh, of a Lya-driven supershell of constant mass depends uniquely on the intrinsic Lya luminosity of the galaxy, L_alpha, the HI column density of the supershell, N_HI, and the shell speed, v_sh, through the scaling relation r_sh ~ L_alpha/(N_HI v_sh^2). We derive mass outflow rates in supershells that reach ~10-100% of the star formation rates of their host galaxies.

preprint2009arXiv

An Electromagnetic Signature of Galactic Black Hole Binaries That Enter Their Gravitational-Wave Induced Inspiral

Mergers of gas-rich galaxies lead to black hole binaries that coalesce as a result of dynamical friction on the ambient gas. Once the binary tightens to <10^3 Schwarzschild radii, its merger is driven by the emission of gravitational waves (GWs). We show that this transition occurs generically at orbital periods of ~1-10 years and an orbital velocity V of a few thousand km/s, with a very weak dependence on the supply rate of gas (V proportional to Mdot^{1/8}). Therefore, as binaries enter their GW-dominated inspiral, they inevitably induce large periodic shifts in the broad emission lines of any associated quasar(s). The probability of finding a binary in tighter configurations scales as V^{-8} owing to their much shorter lifetimes. Systematic monitoring of the broad emission lines of quasars on timescales of months to decades can set a lower limit on the expected rate of GW sources for LISA.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

Concentrating the Dark Matter in Galaxy Clusters through Tidal Stripping of Baryonically-Compressed Galactic Halos

Gravitational lensing observations of massive X-ray clusters imply a steep characteristic density profile marked by a central concentration of dark matter. The observed mass fraction within a projected radius of 150 kpc is twice that found in state-of-the-art dark matter simulations of the standard Lambda-CDM cosmology. A central baryon enhancement that could explain this discrepancy is not observed, leaving a major puzzle. We propose a solution based on the merger histories of clusters. A significant fraction of the final dark matter content of a cluster halo originates within galaxy-sized halos, in which gas can cool and compress the dark matter core to high densities. The subsequent tidal stripping of this compressed dark matter occurs in denser regions that are closer to the center of the cluster halo. Eventually, the originally cooled gas must be dispersed into the intracluster medium through feedback, for consistency with observations that do not find central baryon enhancements in clusters. Still, the early adiabatic compression of the galactic dark matter leaves a net effect on the cluster. Using a simple model for this process, we show that the central cluster profile is substantially modified, potentially explaining the observed discrepancy.

preprint2009arXiv

Constraining reionization using 21 cm observations in combination with CMB and Lyman-alpha forest data

In this paper, we explore the constraints on the reionization history that are provided by current observations of the Lyman-alpha forest and the CMB. Rather than using a particular semi-analytic model, we take the novel approach of parametrizing the ionizing sources with arbitrary functions, and perform likelihood analyses to constrain possible reionization histories. We find model independent conclusions that reionization is likely to be mostly complete by z=8 and that the IGM was 50% ionized at z=9-10. Upcoming low-frequency observations of the redshifted 21 cm line of neutral hydrogen are expected to place significantly better constraints on the hydrogen neutral fraction at 6<z<12. We use our constraints on the reionization history to predict the likely amplitude of the 21 cm power spectrum and show that observations with the highest signal-to-noise ratio will most likely be made at frequencies corresponding to z=9-10. This result provides an important guide to the upcoming 21 cm observations. Finally, we assess the impact that measurement of the neutral fraction will have on our knowledge of reionization and the early source population. Our results show that a single measurement of the neutral fraction mid-way through the reionization era will significantly enhance our knowledge of the entire reionization history.

preprint2009arXiv

Detecting Changing Polarization Structures in Sagittarius A* with High Frequency VLBI

Sagittarius A* is the source of near infrared, X-ray, radio, and (sub)millimeter emission associated with the supermassive black hole at the Galactic Center. In the submillimeter regime, Sgr A* exhibits time-variable linear polarization on timescales corresponding to <10 Schwarzschild radii of the presumed 4 million solar mass black hole. In previous work, we demonstrated the potential for total-intensity (sub)millimeter-wavelength VLBI to detect time-variable -- and periodic -- source structure changes in the Sgr A* black hole system using nonimaging analyses. Here we extend this work to include full polarimetric VLBI observations. We simulate full-polarization (sub)millimeter VLBI data of Sgr A* using a hot-spot model that is embedded within an accretion disk, with emphasis on nonimaging polarimetric data products that are robust against calibration errors. Although the source-integrated linear polarization fraction in the models is typically only a few percent, the linear polarization fraction on small angular scales can be much higher, enabling the detection of changes in the polarimetric structure of Sgr A* on a wide variety of baselines. The shortest baselines track the source-integrated linear polarization fraction, while longer baselines are sensitive to polarization substructures that are beam-diluted by connected-element interferometry. The detection of periodic variability in source polarization should not be significantly affected even if instrumental polarization terms cannot be calibrated out. As more antennas are included in the (sub)mm-VLBI array, observations with full polarization will provide important new diagnostics to help disentangle intrinsic source polarization from Faraday rotation effects in the accretion and outflow region close to the black hole event horizon.

preprint2009arXiv

Estimating the Parameters of Sgr A*'s Accretion Flow Via Millimeter VLBI

Recent millimeter-VLBI observations of Sagittarius A* (Sgr A*) have, for the first time, directly probed distances comparable to the horizon scale of a black hole. This provides unprecedented access to the environment immediately around the horizon of an accreting black hole. We leverage both existing spectral and polarization measurements and our present understanding of accretion theory to produce a suite of generic radiatively inefficient accretion flow (RIAF) models of Sgr A*, which we then fit to these recent millimeter-VLBI observations. We find that if the accretion flow onto Sgr A* is well described by a RIAF model, the orientation and magnitude of the black hole's spin is constrained to a two-dimensional surface in the spin, inclination, position angle parameter space. For each of these we find the likeliest values and their 1-sigma & 2-sigma errors to be a=0(+0.4+0.7), inclination=50(+10+30)(-10-10) degrees, and position angle=-20(+31+107)(-16-29) degrees, when the resulting probability distribution is marginalized over the others. The most probable combination is a=0(+0.2+0.4), inclination=90(-40-50) degrees and position angle=-14(+7+11)(-7-11) degrees, though the uncertainties on these are very strongly correlated, and high probability configurations exist for a variety of inclination angles above 30 degrees and spins below 0.99. Nevertheless, this demonstrates the ability millimeter-VLBI observations, even with only a few stations, to significantly constrain the properties of Sgr A*.

preprint2009arXiv

Imaging the Black Hole Silhouette of M87: Implications for Jet Formation and Black Hole Spin

The silhouette cast by the horizon of the supermassive black hole in M87 can now be resolved with the emerging millimeter very-long baseline interferometry (VLBI) capability. Despite being ~2000 times farther away than SgrA* (the supermassive black hole at the center of the Milky-Way and the primary target for horizon-scale imaging), M87's much larger black hole mass results in a horizon angular scale roughly half that of SgrA*'s, providing another practical target for direct imaging. However, unlike SgrA*, M87 exhibits a powerful radio jet, providing an opportunity to study jet formation physics on horizon scales. We employ a simple, qualitatively correct force-free jet model to explore the expected high-resolution images of M87 at wavelengths of 1.3mm and 0.87mm (230GHz and 345GHz), for a variety of jet parameters. We show that future VLBI data will be able to constrain the size of the jet footprint, the jet collimation rate, and the black hole spin. Polarization will further probe the structure of the jet's magnetic field and its effect on the emitting gas. Horizon-scale imaging of M87 and SgrA* will enable for the first time the empirical exploration of the relationship between the mass and spin of a black hole and the characteristics of the gas inflow/outflow around it.

preprint2009arXiv

Lyman Alpha Blobs as an Observational Signature of Cold Accretion Streams into Galaxies

Recent hydrodynamic simulations of galaxy formation reveal streams of cold (T ~ 1e4 K) gas flowing into the centers of dark matter halos as massive as 1e12-1e13.5 M_sun at redshifts z~1-3. In this paper we show that if > 20% of the gravitational binding energy of the gas is radiated away, then the simulated cold flows are spatially extended Lyman Alpha (Lya) sources with luminosities, Lya line widths, and number densities that are comparable to those of observed Lya blobs. Furthermore, the filamentary structure of the cold flows can explain the wide range of observed Lya blob morphologies. Since the most massive halos form in dense environments, the association of Lya blobs with overdense regions arise naturally. We argue that Lya blobs - even those which are clearly associated with starburst galaxies or quasars - provide direct observational support for the cold accretion mode of galaxies. We discuss various testable predictions of this association.

preprint2009arXiv

Probing the Epoch of Reionization with Milky Way Satellites

While the connection between high-redshift star formation and the local universe has recently been used to understand the observed population of faint dwarf galaxies in the Milky Way (MW) halo, we explore how well these nearby objects can probe the epoch of first light. We construct a detailed, physically motivated model for the MW satellites based on the state-of-the-art Via Lactea II dark-matter simulations. Our model incorporates molecular hydrogen (H_2) cooling in low-mass systems and inhomogeneous photo-heating feedback during the internal reionization of our own galaxy. We find that the existence of MW satellites fainter than M_V ~ -5 is strong evidence for H_2 cooling in low-mass halos, while satellites with -5 > M_V > -9 were affected by hydrogen cooling and photoheating feedback. The age of stars in very low-luminosity systems and the minimum luminosity of these satellites are key predictions of our model. Most of the stars populating the brightest MW satellites could have formed after the epoch of reionization. Our models also predict a significantly larger dispersion in M_300 values than observed and a number of luminous satellites with M_300 as low as 10^6 solar masses.

preprint2009arXiv

Signatures of Relativistic Helical Motion in the Rotation Measures of AGN Jets

Polarization has proved an invaluable tool for probing magnetic fields in relativistic jets. Maps of the intrinsic polarization vectors have provided the best evidence to date for uniform, toroidally dominated magnetic fields within jets. More recently, maps of the rotation measure (RM) in jets have for the first time probed the field geometry of the cool, moderately relativistic surrounding material. In most cases, clear signatures of toroidal magnetic field are detected, corresponding to gradients in RM profiles transverse to the jet. However, in many objects these profiles also display marked asymmetries which are difficult to explain in simple helical jet models. Furthermore, in some cases the RM profiles are strongly frequency and/or time dependent. Here we show that these features may be naturally accounted for by including relativistic helical motion in the jet model. In particular, we are able to reproduce bent RM profiles observed in a variety of jets, frequency dependent RM profile morphologies and even the time dependence of the RM profiles of knots in 3C 273. Finally, we predict that some sources may show reversals in their RM profiles at sufficiently high frequencies, depending upon the the ratio of the components of jet sheath velocity transverse and parallel to the jet. Thus, multi-frequency RM maps promise a novel way in which to probe the velocity structure of relativistic outflows.

preprint2009arXiv

The 21cm Power Spectrum After Reionization

We discuss the 21cm power spectrum (PS) following the completion of reionization. In contrast to the reionization era, this PS is proportional to the PS of mass density fluctuations, with only a small modulation due to fluctuations in the ionization field on scales larger than the mean-free-path of ionizing photons. We derive the form of this modulation, and demonstrate that its effect on the 21cm PS will be smaller than 1% for physically plausible models of damped Ly-alpha systems. In contrast to the 21cm PS observed prior to reionization, in which HII regions dominate the ionization structure, the simplicity of the 21cm PS after reionization will enhance its utility as a cosmological probe by removing the need to separate the PS into physical and astrophysical components. As a demonstration, we consider the Alcock-Paczynski test and show that the next generation of low-frequency arrays could measure the angular distortion of the PS at the percent level for z~3-5.

preprint2009arXiv

The Event Horizon of Sagittarius A*

Black hole event horizons, causally separating the external universe from compact regions of spacetime, are one of the most exotic predictions of General Relativity (GR). Until recently, their compact size has prevented efforts to study them directly. Here we show that recent millimeter and infrared observations of Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, all but requires the existence of a horizon. Specifically, we show that these observations limit the luminosity of any putative visible compact emitting region to below 0.4% of Sgr A*'s accretion luminosity. Equivalently, this requires the efficiency of converting the gravitational binding energy liberated during accretion into radiation and kinetic outflows to be greater than 99.6%, considerably larger than those implicated in Sgr A*, and therefore inconsistent with the existence of such a visible region. Finally, since we are able to frame this argument entirely in terms of observable quantities, our results apply to all geometric theories of gravity that admit stationary solutions, including the commonly discussed f(R) class of theories.

preprint2009arXiv

The Kinetic Sunyaev-Zel'dovich effect of the Milky Way Halo

We calculate the expected imprint of the ionized gas in the Milky-Way halo on the Cosmic Microwave Background (CMB) through the kinetic Sunyaev-Zel'dovich (kSZ) effect. Unlike other Galactic foregrounds, the halo kSZ signature covers the full sky, generates anisotropies on large angular scales, is not accompanied by spectral distortions, and could therefore be confused with primordial CMB anisotropies. We construct theoretical models for various halo components, including smooth diffuse gas, filaments of cold inflowing gas and high velocity clouds. We find that the kSZ effect for all components is above the sensitivity of the Planck satellite, over a range of angular scales. However, the typical halo contribution is well below the cosmic variance noise in the primordial CMB power spectrum. High velocity clouds could dominate the halo contribution and better observational data is required to mask them out. We derive expected kSZ maps based on existing data from tracers of the halo gas distribution, such as 21cm maps of neutral hydrogen and H_alpha maps of recombining gas. The cross-correlation of these maps with the WMAP5 data does not yield any statistically significant signal.

preprint2009arXiv

The Pulsar Contribution to the Gamma-Ray Background

We estimate the contribution of Galactic pulsars, both ordinary and millisecond pulsars (MSPs), to the high-energy (>100 MeV) gamma-ray background. We pay particular attention to the high-latitude part of the background that could be confused with an extragalactic component in existing analyses that subtract a Galactic cosmic-ray model. Our pulsar population models are calibrated to the results of large-scale radio surveys and we employ a simple empirical gamma-ray luminosity calibration to the spin-down rate that provides a good fit to existing data. We find that while ordinary pulsars are expected to contribute only a fraction ~10^-3 of the high-latitude gamma-ray intensity (I_X~1x10^-5 ph s^-1 cm^-2 sr^-1), MSPs could provide a much larger contribution and even potentially overproduce it, depending on the model parameters. We explore these dependences using a range of MSP models as a guide to how gamma-ray measurements can usefully constrain the MSP population. Existing gamma-ray background measurements and source counts already rule out several models. Finally, we show how fluctuations in the gamma-ray sky can be used to distinguish between different sources of the background.

preprint2009arXiv

Upper Limit on Dimming of Cosmological Sources by Intergalactic Grey Dust from the Soft X-ray Background

Active Galactic Nuclei (AGN) produce a dominant fraction (~80%) of the Soft X-ray background (SXB) at photon energies 0.5<E<2 keV. If dust pervaded throughout the intergalactic medium, its scattering opacity would have produced diffuse X-ray halos around AGN. Taking account of known galaxies and galaxy clusters, only a fraction F_halo <10% of the SXB can be in the form of diffuse X-ray halos around AGN. We therefore limit the intergalactic opacity to optical/infrared photons from large dust grains (with radii in the range a=0.2-2.0 mum) to a level tau_GD<0.15(F_halo/10%) to a redshift z~1. Our results are only weakly dependent on the grain size distribution or the redshift evolution of the intergalactic dust. Stacking X-ray images of AGN can be used to improve our constraints and diminish the importance of dust as a source of systematic uncertainty for future supernova surveys which aim to improve the precision on measuring the redshift evolution of the dark energy equation-of-state.

preprint2009arXiv

Will the Large Synoptic Survey Telescope detect extra-solar planetesimals entering the solar system?

Planetesimal formation is a common by-product of the star formation process. Taking the dynamical history of the Solar system as a guideline -- in which the planetesimal belts were heavily depleted due to graviational perburbation with the giant planets -- and assuming similar processes have take place in other planetary systems, one would expect the interestellar space to be filled with extra-solar planetesimals. However, not a single one of these objects has been detected so far entering the Solar system, even though it would clearly be distinguishable from a Solar system comet due to its highly hyperbolic orbit. The Large Synoptic Survey Telescope (LSST) will provide wide coverage maps of the sky to a very high sensitivity, ideal to detect moving objects like comets, both active and inactive. In anticipation of these observations, we estimate how many inactive "interstellar comets" might be detected during the duration of the survey. The calculation takes into account estimates (from observations and models) of the number density of stars, the amount of solids available to form planetesimals, the frequency of planet and planetesimal formation, the efficiency of planetesimal ejection, and the possible size distribution of these small bodies.

preprint2008arXiv

Brightening of an Accretion Disk Due to Viscous Dissipation of Gravitational Waves During the Coalescence of Supermassive Black Holes

Mergers of supermassive black hole binaries release peak power of up to ~10^57 erg/s in gravitational waves (GWs). As the GWs propagate through ambient gas, they induce shear and a small fraction of their power is dissipated through viscosity. The dissipated heat appears as electromagnetic (EM) radiation, providing a prompt EM counterpart to the GW signal. For thin accretion disks, the GW heating rate exceeds the accretion power at distances farther than 10^3 Schwarzschild radii, independently of the accretion rate and viscosity coefficient.

preprint2008arXiv

Evidence for Merger-Driven Activity in the Clustering of High Redshift Quasars

Recently, a very large clustering length has been measured for quasars at a redshift of z~4. In combination with the observed quasar luminosity function we assess the implications of this clustering for the relationship between quasar luminosity and dark matter halo mass. Our analysis allows for non-linearity and finite scatter in the relation between quasar luminosity and halo mass, as well as a luminosity dependent quasar lifetime. The additional novel ingredient in our modelling is the allowance for an excess in the observed bias over the underlying halo bias owing to the merger driven nature of quasar activity. We find that the observations of clustering and luminosity function can be explained only if both of the following conditions hold: (i) The luminosity to halo mass ratio increases with halo mass; (ii) The observed clustering amplitude is in excess of that expected solely from halo bias. The latter result is statistically significant at the 99% level. Taken together, the observations provide compelling evidence for merger driven quasar activity, with a black-hole growth that is limited by feedback. In difference from previous analyses, we show that there could be scatter in the luminosity halo mass relation of up to 1 dex, and that quasar clustering can not be used to estimate the quasar lifetime.

preprint2008arXiv

Methods for detecting flaring structures in Sagittarius A* with high frequency VLBI

The super massive black hole candidate, Sagittarius A*, exhibits variability from radio to X-ray wavelengths on time scales that correspond to < 10 Schwarzschild radii. We survey the potential of millimeter-wavelength VLBI to detect and constrain time variable structures that could give rise to such variations, focusing on a model in which an orbiting hot spot is embedded in an accretion disk. Non-imaging algorithms are developed that use interferometric closure quantities to test for periodicity, and applied to an ensemble of hot-spot models that sample a range of parameter space. We find that structural periodicity in a wide range of cases can be detected on most potential VLBI arrays using modern VLBI instrumentation. Future enhancements of mm/sub-mm VLBI arrays including phased array processors to aggregate VLBI station collecting area, increased bandwidth recording, and addition of new VLBI sites all significantly aid periodicity detection. The methods described herein can be applied to other models of Sagittarius A*, including jet outflows and Magneto-Hydrodynamic accretion simulations.

preprint2007arXiv

Baryonic Acoustic Oscillations in 21cm Emission: A Probe of Dark Energy out to High Redshifts

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21cm emission after reionization, from hydrogen in unresolved damped Ly-alpha absorbers (such as gas rich galaxies) down to a redshift z~3.5. The inferred power spectrum of 21cm fluctuations at all redshifts will show acoustic oscillations, whose co-moving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few percent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency band-passes of >8MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below one percent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5<z<12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (<1%) of the acoustic peak to be made at more moderate redshifts (1.5<z<3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 square degrees. Together with other data sets, observations of 21cm fluctuations will allow full coverage of the acoustic scale from the present time out to z~12.

preprint2007arXiv

Fluctuations in 21cm Emission After Reionization

The fluctuations in the emission of redshifted 21cm photons from neutral inter-galactic hydrogen will provide an unprecedented probe of the reionization era. Conventional wisdom assumes that this 21cm signal disappears as soon as reionization is complete, when little atomic hydrogen is left through most of the volume of the inter-galactic medium (IGM). However observations of damped Ly-alpha absorbers indicate that the fraction of hydrogen in its neutral form is significant by mass at all redshifts. Here we use a physically-motivated model to show that residual neutral gas, confined to dense regions in the IGM with a high recombination rate, will generate a significant post-reionization 21cm signal. We show that the power-spectrum of fluctuations in this signal will be detectable by the first generation of low-frequency observatories at a signal-to-noise that is comparable to that achievable in observations of the reionization era. The statistics of 21cm fluctuations will therefore probe not only the pre-reionization IGM, but rather the entire process of HII region overlap, as well as the appearance of the diffuse ionized IGM.

preprint2007arXiv

GRB Cosmology

Current observations are about to open up a direct window into the final frontier of cosmology: the first billion years in cosmic history when the first stars and galaxies formed. Even before the launch of the James Webb Space Telescope, it might be possible to utilize Gamma-ray Bursts (GRBs) as unique probes of cosmic star formation and the state of the intergalactic medium (IGM) up to redshifts of several tens, when the first (Population III) stars had formed. The Swift mission, or future satellites such as EXIST, might be the first observatories to detect individual Population III stars, provided that massive metal-free stars were able to trigger GRBs. Spectroscopic follow-up observations of the GRB afterglow emission would allow to probe the ionization state and metal enrichment of the IGM as a function of redshift.

preprint2007arXiv

Properties of the Radio-Emitting Gas Around SgrA*

We show that the radial profiles of the temperature and density of the electrons as well as the magnetic field strength around the massive black hole at the Galactic center, SgrA*, may be constrained directly from existing radio data without any need to make prior assumptions about the dynamics of the emitting gas. The observed spectrum and wavelength-dependent angular size of SgrA* indicate that the synchrotron emission originates from an optically-thick plasma of quasi-thermal electrons. We find that the electron temperature rises above the virial temperature within tens of Schwarzschild radii from the black hole, suggesting that the emitting plasma may be outflowing. Constraints on the electron density profile are derived from polarization measurements. Our best-fit results differ from expectations based on existing theoretical models. However, these models cannot be ruled out as of yet due to uncertainties in the source size measurements. Our constraints could tighten considerably with future improvements in the size determination and simultaneous polarization measurements at multiple wavelengths.

preprint2006arXiv

Dwarf Galaxy Formation Was Suppressed By Cosmic Reionization

A large number of faint galaxies, born less than a billion years after the big bang, have recently been discovered. The fluctuations in the distribution of these galaxies contributed to a scatter in the ionization fraction of cosmic hydrogen on scales of tens of Mpc, as observed along the lines of sight to the earliest known quasars. Theoretical simulations predict that the formation of dwarf galaxies should have been suppressed after cosmic hydrogen was reionized, leading to a drop in the cosmic star formation rate. Here we present evidence for this suppression. We show that the post-reionization galaxies which produced most of the ionizing radiation at a redshift z~5.5, must have had a mass in excess of ~10^{10.6+/-0.4} solar masses or else the aforementioned scatter would have been smaller than observed. This limiting mass is two orders of magnitude larger than the galaxy mass that is thought to have dominated the reionization of cosmic hydrogen (~10^8 solar masses). We predict that future surveys with space-based infrared telescopes will detect a population of smaller galaxies that reionized the Universe at an earlier time, prior to the epoch of dwarf galaxy suppression.

preprint2006arXiv

Eavesdropping on Radio Broadcasts from Galactic Civilizations with Upcoming Observatories for Redshifted 21cm Radiation

The question of whether intelligent life exists elsewhere is one of the fundamental unknowns about our Universe. Over the past decade more than 200 extra-solar planets have been discovered, providing new urgency for addressing this question in these or other planetary systems. Independently of this perspective, new radio observatories for cosmology are currently being constructed with the goal of detecting 21cm emission from cosmic hydrogen in the redshift range 6<z<15. For example, the Low-Frequency Demonstrator (LFD) of the Mileura Wide-Field Array (MWA), will cover in 8 kHz bins the entire frequency range of 80-300 MHz, which is perfectly matched to the band over which our civilization emits most of its radio power. We show that this and other low-frequency observatories (culminating with the Square Kilometer Array [SKA]) will be able to detect radio broadcast leakage from an Earth-like civilization out to a distance of ~10^{1-2.7}pc, within a spherical volume containing 10^{3-8}x(Omega_b/4*pi)^alpha stars, where alpha=1 (or 1.5) for a radar beam of solid angle Omega_b that remains steady (or sweeps) across the sky. Such a radio signal will show-up as a series of narrow spectral lines that do not coincide with known atomic or molecular lines. The high spectral resolution attainable with the upcoming observatories will allow to monitor the periodic Doppler shift of the broadcasted lines over the planet's orbital period around the parent star. Determination of the parent star mass through observations of its spectrum could then be used to infer the inclination, semi-major axis and eccentricity of the planet's orbit. This, in turn, will allow to estimate the temperature on the planet's surface and to assess whether it can support liquid water or life as we know it.

preprint2006arXiv

GRB Cosmology and the First Stars

Gamma-Ray Bursts (GRBs) are unique probes of the cosmic star formation history and the state of the intergalactic medium up to the redshifts of the first stars. In particular, the ongoing {\it Swift} mission might be the first observatory to detect individual Population~III stars, provided that the massive, metal-free stars were able to trigger GRBs. {\it Swift} will empirically constrain the redshift at which Population~III star formation was terminated, thus providing crucial input to models of cosmic reionization and metal enrichment.

preprint2006arXiv

Smooth Boundaries to Cosmological HII Regions from Galaxy Clustering

The HII regions around quasars and galaxies at redshifts beyond the epoch of reionisation will provide prime targets for upcoming 21cm campaigns using a new generation of low-frequency radio observatories. Here we show that the boundaries of these HII regions will not be sharp. Rather, the clustering of sources near massive galaxies results in a neutral fraction that rises gradually towards large radii from an interior value near zero. A neutral fraction corresponding to the global background value is typically reached at a distance of 2-5 times the radius of the HII region around the central massive galaxy.

preprint2006arXiv

Thermal Evaporation of Gas from X-ray Clusters

A fraction of the thermal protons in the outer envelope of an X-ray cluster have velocities that exceed the local escape speed from the cluster gravitational potential. The Coulomb mean-free-path of these protons is larger than the virial radius of the cluster at temperatures >2.5 keV. The resulting leakage of suprathermal particles generates a collisionless shock in neighboring voids and fills them with heat and magnetic fields. The momentum flux of suprathermal particles cannot be confined by magnetic tension at the typical field strength in the periphery of cluster halos (<<1 micro-gauss). Over a Hubble time, thermal evaporation could drain up to a tenth of the cluster gas at its virial temperature. The evaporated fraction could increase dramatically if additional heat is deposited into the gas by cluster mergers, active galactic nuclei or supernovae. Thermal evaporation is not included in existing cosmological simulations since they are based on the fluid approximation. Measurements of the baryon mass fraction in the outer envelopes of hot clusters (through their Sunyaev-Zel'dovich effect or X-ray emission) can be used to empirically constrain their evaporation rate.

preprint2005arXiv

Constraints on the Process that Regulates the Growth of Supermassive Black Holes Based on the Intrinsic Scatter in the M_bh-sigma Relation

We show that the observed scatter in the relations between the mass of supermassive black holes (SMBHs), M_bh, and the velocity dispersion sigma or mass M_sph of their host spheroid, place interesting constraints on the process that regulates SMBH growth in galaxies. When combined with the observed properties of early-type SDSS galaxies, the observed intrinsic scatters imply that SMBH growth is regulated by the spheroid velocity dispersion rather than its mass. The M_bh-M_sph relation is therefore a by-product of a more fundamental M_bh-sigma relation. We construct a theoretical model for the scatter among baryon modified dark matter halo profiles, out of which we generate a population of spheroid hosts and show that these naturally lead to a relation between effective radius and velocity dispersion of the form R_sph ~ sigma^1.5 with a scatter of ~0.2dex, in agreement with the corresponding projection of the fundamental plane for early type galaxies in SDSS. At the redshift of formation, our model predicts the minimum scatter that SMBHs can have at fixed velocity dispersion or spheroid mass under different formation scenarios. We also estimate the additional scatter that is introduced into these relations through collisionless mergers of purely stellar spheroids at z<1. We find that the observed scatter in the M_bh-sigma and M_bh-M_sph relations preclude the properties of dark matter halos from being the governing factor in SMBH growth. Finally, we show that SMBH growth governed by the properties of the host spheroid can lead to the observed values of scatter in the M_bh-sigma and M_bh-M_sph relations, only if the SMBH growth is limited by momentum or energy feedback over the dynamical time of the host spheroid.

preprint2004arXiv

A Size of ~10 Mpc for the Ionized Bubbles at the End of Cosmic Reionization

The first galaxies to appear in the universe at redshifts z>20 created ionized bubbles in the intergalactic medium of neutral hydrogen left over from the Big-Bang. It is thought that the ionized bubbles grew with time, surrounded clusters of dwarf galaxies and eventually overlapped quickly throughout the universe over a narrow redshift interval near z~6. This event signaled the end of the reionization epoch when the universe was a billion years old. Measuring the hitherto unknown size distribution of the bubbles at their final overlap phase is a focus of forthcoming observational programs aimed at highly redshifted 21cm emission from atomic hydrogen. Here we show that the combined constraints of cosmic variance and causality imply an observed bubble size at the end of the overlap epoch of ~10 physical Mpc, and a scatter in the observed redshift of overlap along different lines-of-sight of ~0.15. This scatter is consistent with observational constraints from recent spectroscopic data on the farthest known quasars. Our novel result implies that future radio experiments should be tuned to a characteristic angular scale of ~0.5 degrees and have a minimum frequency band-width of ~8 MHz for an optimal detection of 21cm flux fluctuations near the end of reionization.

preprint2004arXiv

Cosmic Hydrogen Was Significantly Neutral a Billion Years After the Big Bang

The ionization fraction of cosmic hydrogen, left over from the big bang, provides crucial fossil evidence for when the first stars and quasar black holes formed in the infant universe. Spectra of the two most distant quasars known show nearly complete absorption of photons with wavelengths shorter than the Ly-alpha transition of neutral hydrogen, indicating that hydrogen in the intergalactic medium (IGM) had not been completely ionized at a redshift z~6.3, about a billion years after the big bang. Here we show that the radii of influence of ionizing radiation from these quasars imply that the surrounding IGM had a neutral hydrogen fraction of tens of percent prior to the quasar activity, much higher than previous lower limits of ~0.1%. When combined with the recent inference of a large cumulative optical depth to electron scattering after cosmological recombination from the WMAP data, our result suggests the existence of a second peak in the mean ionization history, potentially due to an early formation episode of the first stars.

preprint2004arXiv

Element segregation in giant galaxies and X-ray clusters

We examine the process of element segregation by gravity in giant elliptical galaxies and X-ray clusters. We solve the full set of flow equations developed by Burgers for a multi-component fluid, under the assumption that magnetic fields slow the diffusion by a constant factor F_B. Compared to the previous calculations that neglected the residual heat flow terms, we find that diffusion is faster by ~20%. In clusters we find that the diffusion changes the local abundance typically by factors of 1+0.3(T/10^8 K)^1.5/F_B and 1+0.15(T/10^8 K)^1.5/F_B for helium and heavy elements, respectively, where T is the gas temperature. In elliptical galaxies, the corresponding factors are 1+0.2(T/10^7 K)^1.5/F_B and 1+0.1(T/10^7 K)^1.5/F_B, respectively. If the suppression factor F_B is modest, diffusion could significantly affect observational properties of hot X-ray clusters and cD galaxies. In particular, diffusion steepens the baryon distribution, increases the total X-ray luminosity, and changes the spectrum and evolution of stars that form out of the helium-rich gas. Detection of these diffusion signatures would allow to gauge the significance of the magnetic fields, which also inhibit thermal heat conduction as a mechanism for eliminating cooling flows in the same environments.

preprint2004arXiv

Implications of the Measured Image Size for the Radio Afterglow of GRB 030329

We use data on the image size of the radio afterglow of GRB 030329 (Taylor et al. 2004) to constrain the physical parameters of this explosion. Together with the observed broad band spectrum, this data over-constrains the physical parameters, thus enabling to test different GRB jet models for consistency. We consider two extreme models for the lateral spreading of the jet: model 1 with relativistic expansion in the local rest frame, and model 2 with little lateral expansion as long as the jet is highly relativistic. We find that both models are consistent with the data for a uniform external medium, while for a stellar wind environment model 1 is consistent with the data but model 2 is disfavored by the data. Our derivations can be used to place tighter constraints on the dynamics and structure of GRB jets in future afterglows, following a denser monitoring campaign for the temporal evolution of their image size.

preprint2004arXiv

Measuring the Small-Scale Power Spectrum of Cosmic Density Fluctuations Through 21 cm Tomography Prior to the Epoch of Structure Formation

The thermal evolution of the cosmic gas decoupled from that of the cosmic microwave background (CMB) at a redshift z~200. Afterwards and before the first stars had formed, the cosmic neutral hydrogen absorbed the CMB flux at its resonant 21cm spin-flip transition. We calculate the evolution of the spin temperature for this transition and the resulting anisotropies that are imprinted on the CMB sky due to linear density fluctuations during this epoch. These anisotropies at an observed wavelength of 10.56[(1+z)/50] meters, contain an amount of information that is orders of magnitude larger than any other cosmological probe. Their detection, although challenging, could tightly constrain any possible running of the spectral index from inflation (as suggested by WMAP), small deviations from Gaussianity, or any significant contribution from neutrinos or warm dark matter to the cosmic mass budget.

preprint2004arXiv

Probing the Magnetic Field Structure in Gamma-Ray Bursts through Dispersive Plasma Effects on the Afterglow Polarization

(Abr) The origin and structure of magnetic fields in Gamma-Ray Burst (GRB) fireball plasmas are two of the most important open questions in all GRB models. We show that the structure and strength of the magnetic field may be constrained by radio and IR observations of the early afterglow, where plasma effects on the polarization of propagating radiation are significant. We calculate these propagation effects for cold and relativistic plasmas, and find that in the presence of a uniform equipartition field the degree of linear polarization is suppressed, and circular polarization prevails at low frequencies, nu < 1-3 GHz, (2x10^11 Hz < nu < few x 10^14 Hz) in the forward (reverse) shock. At higher frequencies linear polarization dominates. At the frequency of the transition between circular and linear polarization, the net level of polarization is minimal, ~10-20%. These features are nearly independent of the circumburst density. The transition frequency is smaller by a factor of ~10 when the uniform field is much weaker than equipartition. The dependence of these results on viewing geometry, outflow collimation and magnetic field orientation is discussed. When the configuration of the field is entangled over length scales much smaller than the extent of the emitting plasma, the aforementioned effects should not be observed and a linear polarization at the few % level is expected. Polarimetric observations during the early afterglow, and particularly of the reverse shock emission, may therefore place strong constraints on the structure and strength of the magnetic field within the fireball plasma.

preprint2003arXiv

Detection of Gravitational Waves from the Coalescence of Population-III Remnants with Advanced LIGO

The comoving mass density of massive black hole (MBH) remnants from pre-galactic star formation could have been similar in magnitude to the mass-density of supermassive black holes (SMBHs) in the present-day universe. We show that the fraction of MBHs that coalesce during the assembly of SMBHs can be extracted from the rate of ring-down gravitational waves that are detectable by Advanced LIGO. Based on the SMBH formation history inferred from the evolution of the quasar luminosity function, we show that an observed event rate of 1 per year will constrain the SMBH mass fraction that was contributed by MBHs coalescence down to a level of ~10^-6 for 20 solar mass MBH remnants (or ~10^-4 for 260 solar mass remnants).

preprint2003arXiv

Spectral Signature of Cosmological Infall of Gas Around the First Quasars

Recent observations have shown that, only a billion years after the Big Bang, the Universe was already lit up by bright quasars fuelled by the infall of gas onto supermassive black holes. The masses of these early black holes are inferred from their luminosities to be > 10^9 solar masses, which is a difficult theoretical challenge to explain. Like nearby quasars, the early objects could have formed in the central cores of massive host galaxies. The formation of these hosts could be explained if, like local large galaxies, they were assembled gravitationally inside massive (> 10^12 solar mass) halos of dark matter. There has hitherto been no observational evidence for the presence of these massive hosts or their surrounding halos. Here we show that the cosmic gas surrounding each halo must respond to its strong gravitational pull, where absorption by the infalling hydrogen produces a distinct spectral signature. That signature can be seen in recent data.

preprint2003arXiv

Was the Universe Reionized by Massive Population-III Stars?

The WMAP satellite has measured a large optical depth to electron scattering after cosmological recombination of 0.17+-0.04, implying significant reionization of the primordial gas only ~200 million years after the big bang. However, the most recent overlap of intergalactic HII regions must have occured at z<9 based on the Lyman-alpha forest constraint on the thermal history of the intergalactic medium. Here we argue that a first generation of metal-free stars with a heavy (rather than Salpeter) mass function is therefore required to account for much of the inferred optical depth. This conclusion holds if feedback regulates star formation in early dwarf galaxies as observed in present-day dwarfs.

preprint2002arXiv

Are X-ray Clusters Cooled by Heat Conduction to the Surrounding Intergalactic Medium?

We show that X-ray clusters would have cooled substantially over a Hubble time by transport of heat from their hot interior to the their envelope, if the heat conductivity had not been heavily suppressed relative to the Spitzer value due to magnetic fields. The suppression is required in order for the observed abundance of hot X-ray clusters to be consistent with predictions from popular cosmological models. If a similar or stronger suppression factor applies to cluster cores, then thermal conduction can not be the mechanism that prevents cooling flows there.

preprint2002arXiv

Gravitational lenses magnify up to one third of the most distant quasars

Exceptionally bright quasars with redshifts up to z=6.28 have recently been discovered. Quasars are thought to be powered by the accretion of gas onto supermassive black holes at the centres of galaxies. Their maximum (Eddington) luminosity is proportional to the mass of the black hole, and so these bright quasars are inferred to have black holes with masses of more than a few billion solar masses. The existence of such massive black holes poses a challenge to models for the formation of structures in the early Universe, as it requires that the black holes would grow so massive in less than a billion years after the Big Bang. Here we show that up to a third of known quasars with z~6 will have their observed flux magnified by a factor of 10 or more through gravitational lensing by galaxies along the line of sight. The inferred abundance of quasar host galaxies, as well as the luminosity density provided by the quasars, are therefore substantially overestimated.

preprint2002arXiv

Measuring the Size of Quasar Broad-Line Clouds Through Time Delay Light-Curve Anomalies of Gravitational Lenses

Intensive monitoring campaigns have recently attempted to measure the time delays between multiple images of gravitational lenses. Some of the resulting light-curves show puzzling low-level, rapid variability which is unique to individual images, superimposed on top of (and concurrent with) longer time-scale intrinsic quasar variations which repeat in all images. We demonstrate that both the amplitude and variability time-scale of the rapid light-curve anomalies, as well as the correlation observed between intrinsic and microlensed variability, are naturally explained by stellar microlensing of a smooth accretion disk which is occulted by optically-thick broad-line clouds. The rapid time-scale is caused by the high velocities of the clouds (~5x10^3 km/s), and the low amplitude results from the large number of clouds covering the magnified or demagnified parts of the disk. The observed amplitudes of variations in specific lenses implies that the number of broad-line clouds that cover ~10% of the quasar sky is ~10^5 per 4 pi steradian. This is comparable to the expected number of broad line clouds in models where the clouds originate from bloated stars.

preprint2002arXiv

The Expected Redshift Distribution of Gamma-Ray Bursts

We predict the redshift distribution of Gamma-Ray Bursts (GRBs) assuming that they trace the cosmic star formation history. We find that a fraction >~50% of all GRBs on the sky originate at a redshift z >~5, even though the fraction of the total stellar mass formed by z~5 is only ~15%. These two fractions are significantly different because they involve different cosmological factors when integrating the star formation rate over redshift. Hence, deep observations of transient events, such as GRB afterglows or supernovae, provide an ideal strategy for probing the high-redshift universe. We also show that the intrinsic distribution of GRB durations is bimodal but significantly narrower and shifted towards shorter durations than the observed distribution.

preprint2001arXiv

In the Beginning: The First Sources of Light and the Reionization of the Universe

The formation of the first stars and quasars marks the transformation of the universe from its smooth initial state to its clumpy current state. In popular cosmological models, the first sources of light began to form at redshift 30 and reionized most of the hydrogen in the universe by redshift 7. Current observations are at the threshold of probing the hydrogen reionization epoch. The study of high-redshift sources is likely to attract major attention in observational and theoretical cosmology over the next decade.

preprint2000arXiv

Gamma-Ray Background from Structure Formation in the Intergalactic Medium

The universe is filled with a diffuse and isotropic extragalactic background of gamma-ray radiation, containing roughly equal energy flux per decade in photon energy between 3 MeV-100 GeV. The origin of this background is one of the unsolved puzzles in cosmology. Less than a quarter of the gamma-ray flux can be attributed to unresolved discrete sources, but the remainder appears to constitute a truly diffuse background whose origin has hitherto been mysterious. Here we show that the shock waves induced by gravity during the formation of large-scale structure in the intergalactic medium, produce a population of highly-relativistic electrons with a maximum Lorentz factor above 10^7. These electrons scatter a small fraction of the microwave background photons in the present-day universe up to gamma-ray energies, thereby providing the gamma-ray background. The predicted diffuse flux agrees with the observed background over more than four decades in photon energy, and implies a mean cosmological density of baryons which is consistent with Big-Bang nucleosynthesis.

preprint1999arXiv

Radio Detection of Old GRB Remnants in the Local Universe

We investigate the hydrodynamic evolution of Gamma--Ray Burst (GRB) remnants at late times (>> 1 year), by extending the standard fireball model into the nonrelativistic Sedov--Taylor regime. We calculate the associated synchrotron luminosity as a function of remnant age. The radio luminosity of old remnants depends strongly on the power-law index p= (dlogN/dlogE) of the kinetic energy distribution of shock-accelerated electrons. In addition, the detection probability of remnants in nearby galaxies depends on the GRB rate in the local universe. We predict that if p<2.6 and the BATSE-calibrated GRB rate per comoving volume is independent of redshift (despite the strong evolution in the cosmic star formation rate) then the VLA can discover a 10^3 years-old GRB remnant in the Virgo cluster. Nearby GRB remnants could be searched for and resolved at radio frequencies, and then identified as distinct from supernova remnants based on the detection of the extended emission of high-ionization lines in follow--up optical observations.

preprint1997arXiv

Getting Around Cosmic Variance

Cosmic microwave background (CMB) anisotropies probe the primordial density field at the edge of the observable Universe. There is a limiting precision (``cosmic variance'') with which anisotropies can determine the amplitude of primordial mass fluctuations. This arises because the surface of last scatter (SLS) probes only a finite two-dimensional slice of the Universe. Probing other SLSs observed from different locations in the Universe would reduce the cosmic variance. In particular, the polarization of CMB photons scattered by the electron gas in a cluster of galaxies provides a measurement of the CMB quadrupole moment seen by the cluster. Therefore, CMB polarization measurements toward many clusters would probe the anisotropy on a variety of SLSs within the observable Universe, and hence reduce the cosmic-variance uncertainty.

preprint1997arXiv

Gravitational Lensing of Quasars by Spiral Galaxies

Gravitational lensing by a spiral galaxy occurs when the line-of-sight to a background quasar passes within a few kpc from the center of the galactic disk. Since galactic disks are rich in neutral hydrogen, the quasar spectrum will likely be marked by a damped Lyman-alpha absorption trough at the lens redshift. Therefore, the efficiency of searches for gravitational lensing with sub-arcsecond splitting can be enhanced by 1-2 orders of magnitude by focusing on a subset of all bright quasars which show low-redshift (z<1) strong Lyman-alpha absorption (N>10^{21} cm^{-2}}). The double-image signature of lensing could, in principle, be identified spectroscopically and without the need for high-resolution imaging. The absorption spectrum of a spiral lens would show a generic double-step profile due to the superposition of the two absorption troughs of the different images. Finally, we note that searches for microlensing signatures of quasars with damped Lyman-alpha absorption could calibrate the fraction of MACHOs in galactic halos at high redshift.

preprint1995arXiv

COSMOLOGICAL ORIGIN OF QUASARS

Quasars are the earliest objects known to exist. We examine their origin in the context of popular models for structure formation in the universe. We show that seeds for quasar black holes could have originated from the initial cosmological collapse of overdense regions with unusually small rotation. Most of these seeds have a mass of order 10^6 solar masses, just above the cosmological Jeans mass. For Cold Dark Matter cosmologies, we find of order one seed black hole per bright galaxy. After the galaxy forms, the seeds inside its bulge sink to the center by dynamical friction. We also describe a few empirical methods to study the properties of quasars and their environments. These include: probing quasar hosts through [C II] emission, finding quasar lifetimes from the ``proximity effect'' along two lines of sight, and measuring the amplitude of clustering at high redshifts through the detection of Lyman-alpha clouds beyond the quasar redshift.

preprint1994arXiv

Empirical Constraints on Cosmological Gamma-Ray Bursts

We place empirical constraints on the physical properties of $γ-$ray burst events at cosmological distances. In particular we derive probability distributions for the radiation energy $E_γ$, the minimum Lorentz factor $γ_{\rm min}$, the maximum baryonic mass $M_{\rm max}$ and the upper bound on the surrounding gas density $n_{\rm max}$ in the events, based on 169 bursts from the first BATSE catalog. Using peak flux as a distance indicator we probe bursts where the constraints are stronger than average. The resulting variance and skewness of the cosmological probability distributions are calculated in addition to their mean values: $\langle E_γ\rangle=4\times 10^{51} h^{-2} {\rm erg}$, $\langle γ_{\rm min}\rangle= 5\times 10^2$, $\langle M_{\rm max}\rangle=10^{-5}ξ^{-1} M_\odot$, and $\langle (γ/γ_{\rm min})^{-5} n_{\rm max}\rangle\lsim 10^{2} ξ^{-1}{\rm cm^{-3}}$, where $ξ$ is the fraction of the total energy which is converted to $γ-$rays. The distribution of burst energies ends at about $10^{53} {\rm erg}$, close to the binding energy of a neutron star.

preprint1993arXiv

Are Gamma-Ray Bursts at Cosmological Distances Optically-Thin?

The observed spatial distribution of $γ-$ray bursts indicates that they probably originate at cosmological distances. At this distance scale their variability timescale and flux above MeV imply an initial optical-depth to pair production $> 10^{10}$. This appears to be in conflict with their highly non-thermal spectra. We show that this difficulty can be removed if axion bursts from supernovae are converted to $γ-$rays over cosmological distances. Nonthermal bursts with the relevant flux, duration, variability and spectra are obtained just for the range of axion masses of $10^{-5}-10^{-4}{\rm eV}$ that accounts for the cold dark matter in the universe. The observed rate of $γ-$ray bursts implies that axions should be converted efficiently to photons in only one out of $\sim 10^{4}$ supernovae.