Source author record

Roger D. Blandford

Roger D. Blandford 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

24works
10topics
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

24 published item(s)

preprint2022arXiv

Incremental Fermi Large Area Telescope Fourth Source Catalog

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

preprint2020arXiv

Characterizing the Gamma-Ray Variability of the Brightest Flat Spectrum Radio Quasars Observed with the Fermi LAT

Almost 10 yr of $γ$-ray observations with the Fermi Large Area Telescope (LAT) have revealed extreme $γ$-ray outbursts from flat spectrum radio quasars (FSRQs), temporarily making these objects the brightest $γ$-ray emitters in the sky. Yet, the location and mechanisms of the $γ$-ray emission remain elusive. We characterize long-term $γ$-ray variability and the brightest $γ$-ray flares of six FSRQs. Consecutively zooming in on the brightest flares, which we identify in an objective way through Bayesian blocks and a hill-climbing algorithm, we find variability on subhour time scales and as short as minutes for two sources in our sample (3C279, CTA102) and weak evidence for variability at time scales less than the Fermi satellite's orbit of 95 minutes for PKS1510-089 and 3C454.3. This suggests extremely compact emission regions in the jet. We do not find any signs for $γ$-ray absorption in the broad-line region (BLR), which indicates that $γ$-rays are produced at distances greater than hundreds of gravitational radii from the central black hole. This is further supported by a cross-correlation analysis between $γ$-ray and radio/millimeter light curves, which is consistent with $γ$-ray production at the same location as the millimeter core for 3C273, CTA102, and 3C454.3. The inferred locations of the $γ$-ray production zones are still consistent with the observed decay times of the brightest flares if the decay is caused by external Compton scattering with BLR photons. However, the minute-scale variability is challenging to explain in such scenarios.

preprint2020arXiv

The Chiral Puzzle of Life

Biological molecules chose one of two structurally, chiral systems which are related by reflection in a mirror. It is proposed that this choice was made, causally, by magnetically polarized and physically chiral cosmic-rays, which are known to have a large role in mutagenesis. It is shown that the cosmic rays can impose a small, but persistent, chiral bias in the rate at which they induce structural changes in simple, chiral monomers that are the building blocks of biopolymers. A much larger effect should be present with helical biopolymers, in particular, those that may have been the progenitors of RNA and DNA. It is shown that the interaction can be both electrostatic, just involving the molecular electric field, and electromagnetic, also involving a magnetic field. It is argued that this bias can lead to the emergence of a single, chiral life form over an evolutionary timescale. If this mechanism dominates, then the handedness of living systems should be universal. Experiments are proposed to assess the efficacy of this process.

preprint2019arXiv

H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3σ$ tension between early and late-Universe probes

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1σ$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3σ$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Λ$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Λ$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

preprint2016arXiv

Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9σ$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $Λ$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

preprint2016arXiv

Kinetic simulations of the lowest-order unstable mode of relativistic magnetostatic equilibria

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

preprint2016arXiv

Kinetic study of radiation-reaction-limited particle acceleration during the relaxation of unstable force-free equilibria

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

preprint2015arXiv

Probing the inner kpc of massive galaxies with strong gravitational lensing

We examine the prospects of detecting demagnified images of gravitational lenses in observations of strongly lensed mm-wave molecular emission lines with ALMA. We model the lensing galaxies as a superposition of a dark matter component, a stellar component, and a central supermassive black hole and assess the detectability of the central images for a range of relevant parameters (e.g., stellar core, black hole mass, and source size). We find that over a large range of plausible parameters, future deep observations of lensed molecular lines with ALMA should enable detection of the central images at $\gtrsim 3σ$ significance. We use a Fisher analysis to examine the constraints that could be placed on these parameters in various scenarios and find that for large stellar cores, both the core size and the mass of the central SMBHs can be accurately measured. We also study the prospects for detecting binary SMBHs with such observations and find that only under rare conditions and with very long integrations ($\sim$40-hr) the masses of both SMBHs may be measured using the distortions of central images.

preprint2015arXiv

Spontaneous decay of periodic magnetostatic equilibria

In order to understand the conditions which lead a highly magnetized, relativistic plasma to become unstable, and in such cases how the plasma evolves, we study a prototypical class of magnetostatic equilibria where the magnetic field satisfies $\nabla \times\mathbf B = α\mathbf B$, where $α$ is spatially uniform, on a periodic domain. Using numerical solutions we show that generic examples of such equilibria are unstable to ideal modes (including incompressible ones) which are marked by exponential growth in the linear phase. We characterize the unstable mode, showing how it can be understood in terms of merging magnetic and current structures, and explicitly demonstrate its instability using the energy principle. Following the nonlinear evolution of these solutions, we find that they rapidly develop regions with relativistic velocities and electric fields of comparable magnitude to the magnetic field, liberating magnetic energy on dynamical timescales and eventually settling into a configuration with the largest allowable wavelength. These properties make such solutions a promising setting for exploring the mechanisms behind extreme cosmic sources of gamma rays.

preprint2013arXiv

Improving the precision of time-delay cosmography with observations of galaxies along the line of sight

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

preprint2013arXiv

Reconstructing the Lensing Mass in the Universe from Photometric Catalogue Data

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2012arXiv

Alignment of Magnetized Accretion Disks and Relativistic Jets with Spinning Black Holes

Accreting black holes (BHs) produce intense radiation and powerful relativistic jets, which are affected by the BH's spin magnitude and direction. While thin disks might align with the BH spin axis via the Bardeen-Petterson effect, this does not apply to jet systems with thick disks. We used fully three-dimensional general relativistic magnetohydrodynamical simulations to study accreting BHs with various BH spin vectors and disk thicknesses with magnetic flux reaching saturation. Our simulations reveal a "magneto-spin alignment" mechanism that causes magnetized disks and jets to align with the BH spin near BHs and further away to reorient with the outer disk. This mechanism has implications for the evolution of BH mass and spin, BH feedback on host galaxies, and resolved BH images for SgrA* and M87.

preprint2012arXiv

General Relativistic Magnetohydrodynamic Simulations of Magnetically Choked Accretion Flows around Black Holes

Black hole (BH) accretion flows and jets are qualitatively affected by the presence of ordered magnetic fields. We study fully three-dimensional global general relativistic magnetohydrodynamic (MHD) simulations of radially extended and thick (height $H$ to cylindrical radius $R$ ratio of $|H/R|\sim 0.2--1$) accretion flows around BHs with various dimensionless spins ($a/M$, with BH mass $M$) and with initially toroidally-dominated ($ϕ$-directed) and poloidally-dominated ($R-z$ directed) magnetic fields. Firstly, for toroidal field models and BHs with high enough $|a/M|$, coherent large-scale (i.e. $\gg H$) dipolar poloidal magnetic flux patches emerge, thread the BH, and generate transient relativistic jets. Secondly, for poloidal field models, poloidal magnetic flux readily accretes through the disk from large radii and builds-up to a natural saturation point near the BH. For sufficiently high $|a/M|$ or low $|H/R|$ the polar magnetic field compresses the inflow into a geometrically thin highly non-axisymmetric "magnetically choked accretion flow" (MCAF) within which the standard linear magneto-rotational instability is suppressed. The condition of a highly-magnetized state over most of the horizon is optimal for the Blandford-Znajek mechanism that generates persistent relativistic jets with $\gtrsim 100$% efficiency for $|a/M|\gtrsim 0.9$. A magnetic Rayleigh-Taylor and Kelvin-Helmholtz unstable magnetospheric interface forms between the compressed inflow and bulging jet magnetosphere, which drives a new jet-disk quasi-periodic oscillation (JD-QPO) mechanism. The high-frequency QPO has spherical harmonic $|m|=1$ mode period of $τ\sim 70GM/c^3$ for $a/M\sim 0.9$ with coherence quality factors $Q\gtrsim 10$. [abridged]

preprint2012arXiv

Nomads of the Galaxy

We estimate that there may be up to ~10^5 compact objects in the mass range 10^{-8} -10^{-2} solar mass per main sequence star that are unbound to a host star in the Galaxy. We refer to these objects as nomads; in the literature a subset of these are sometimes called free-floating or rogue planets. Our estimate for the number of Galactic nomads is consistent with a smooth extrapolation of the mass function of unbound objects above the Jupiter-mass scale, the stellar mass density limit, and the metallicity of the interstellar medium. We analyze the prospects for detecting nomads via Galactic microlensing. The Wide-Field Infrared Survey Telescope (WFIRST) will measure the number of nomads per main sequence star greater than the mass of Jupiter to ~ 13%, and the corresponding number greater than the mass of Mars to ~25%. All-sky surveys such as GAIA and LSST can identify nomads greater than about the mass of Jupiter. We suggest a dedicated drift scanning telescope that covers approximately 100 square degrees in the Southern hemisphere could identify nomads as small as 10^{-8} solar mass via microlensing of bright stars with characteristic lightcurve timescales of a few seconds.

preprint2012arXiv

The Chandra view of the Largest Quasar Lens SDSS J1029+2623

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

preprint2012arXiv

The Hidden Fortress: Structure and substructure of the complex strong lensing cluster SDSS J1029+2623

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be θ_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

preprint2009arXiv

Cosmic Evolution of Black Holes and Spheroids. IV. The BH Mass - Spheroid Luminosity Relation

From high-resolution images of 23 Seyfert-1 galaxies at z=0.36 and z=0.57 obtained with the Near Infrared Camera and Multi-Object Spectrometer on board the Hubble Space Telescope (HST), we determine host-galaxy morphology, nuclear luminosity, total host-galaxy luminosity and spheroid luminosity. Keck spectroscopy is used to estimate black hole mass (M_BH). We study the cosmic evolution of the M_BH-spheroid luminosity (L_sph) relation. In combination with our previous work, totaling 40 Seyfert-1 galaxies, the covered range in BH mass is substantially increased, allowing us to determine for the first time intrinsic scatter and correct evolutionary trends for selection effects. We re-analyze archival HST images of 19 local reverberation-mapped active galaxies to match the procedure adopted at intermediate redshift. Correcting spheroid luminosity for passive luminosity evolution and taking into account selection effects, we determine that at fixed present-day V-band spheroid luminosity, M_BH/L_sph \propto (1+z)^(2.8+/-1.2). When including a sample of 44 quasars out to z=4.5 taken from the literature, with luminosity and BH mass corrected to a self-consistent calibration, we extend the BH mass range to over two orders of magnitude, resulting in M_BH/L_sph \propto (1+z)^(1.4+/-0.2). The intrinsic scatter of the relation, assumed constant with redshift, is 0.3+/-0.1 dex (<0.6 dex at 95% CL). The evolutionary trend suggests that BH growth precedes spheroid assembly. Interestingly, the M_BH-total host-galaxy luminosity relation is apparently non-evolving. It hints at either a more fundamental relation or that the spheroid grows by a redistribution of stars. However, the high-z sample does not follow this relation, indicating that major mergers may play the dominant role in growing spheroids above z~1.

preprint2009arXiv

Direct Observation of Cosmic Strings via their Strong Gravitational Lensing Effect: II. Results from the HST/ACS Image Archive

We have searched 4.5 square degrees of archival HST/ACS images for cosmic strings, identifying close pairs of similar, faint galaxies and selecting groups whose alignment is consistent with gravitational lensing by a long, straight string. We find no evidence for cosmic strings in five large-area HST treasury surveys (covering a total of 2.22 square degrees), or in any of 346 multi-filter guest observer images (1.18 square degrees). Assuming that simulations ccurately predict the number of cosmic strings in the universe, this non-detection allows us to place upper limits on the unitless Universal cosmic string tension of G mu/c^2 < 2.3 x 10^-6, and cosmic string density of Omega_s < 2.1 x 10^-5 at the 95% confidence level (marginalising over the other parameter in each case). We find four dubious cosmic string candidates in 318 single filter guest observer images (1.08 square degrees), which we are unable to conclusively eliminate with existing data. The confirmation of any one of these candidates as cosmic strings would imply G mu/c^2 ~ 10^-6 and Omega_s ~ 10^-5. However, we estimate that there is at least a 92% chance that these string candidates are random alignments of galaxies. If we assume that these candidates are indeed false detections, our final limits on G mu/c^2 and Omega_s fall to 6.5 x 10^-7 and 7.3 x 10^-6. Due to the extensive sky coverage of the HST/ACS image archive, the above limits are universal. They are quite sensitive to the number of fields being searched, and could be further reduced by more than a factor of two using forthcoming HST data.

preprint2009arXiv

On the Energy Spectra of GeV/TeV Cosmic Ray Leptons

Recent observations of cosmic ray electrons from several instruments have revealed various degrees of deviation in the measured electron energy distribution from a simple power-law, in a form of an excess around TeV energies. An even more prominent deviation has been observed in the fraction of cosmic ray positrons around 100 GeV energies. In this paper we show that the observed excesses in the electron spectrum may be easily re-produced without invoking any unusual sources other than the general diffuse Galactic components of cosmic rays. The primary physical effect involved is the Klein-Nishina suppression of the electron cooling rate around TeV energies. With a very reasonable choice of the model parameters characterizing the local interstellar medium, we can reproduce the most recent observations by Fermi and HESS experiments. We also find that high positron fraction increasing with energy, as claimed by the PAMELA experiment, cannot be explained in our model with the conservative set of the model parameters. We are able, however, to reproduce the PAMELA results assuming high values of the starlight and interstellar gas densities, which would be more appropriate for vicinities of supernova remnants. A possible solution to this problem may be that cosmic rays undergo most of their interactions near their sources due to the efficient trapping in the far upstream of supernova shocks by self-generated, cosmic ray-driven turbulence.

preprint2009arXiv

Understanding the Geometry of Astrophysical Magnetic Fields

Faraday rotation measurements have provided an invaluable technique with which to measure the properties of astrophysical magnetized plasmas. Unfortunately, typical observations provide information only about the density-weighted average of the magnetic field component parallel to the line of sight. As a result, the magnetic field geometry along the line of sight, and in many cases even the location of the rotating material, is poorly constrained. Frequently, interpretations of Faraday rotation observations are dependent upon underlying models of the magnetic field being probed (e.g., uniform, turbulent, equipartition). However, we show that at sufficiently low frequencies, specifically below roughly 13 (RM/rad m^-2)^(1/4) (B/G)^(1/2) MHz, the character of Faraday rotation changes, entering what we term the ``super-adiabatic regime'' in which the rotation measure is proportional to the integrated absolute value of the line-of-sight component of the field. As a consequence, comparing rotation measures at high frequencies with those in this new regime provides direct information about the geometry of the magnetic field along the line of sight. Furthermore, the frequency defining the transition to this new regime, nu_SA, depends directly upon the local electron density and magnetic field strength where the magnetic field is perpendicular to the line of sight, allowing the unambiguous distinction between Faraday rotation within and in front of the emission region. Typical values of nu_SA range from 10 kHz to 10 GHz, depending upon the details of the Faraday rotating environment. In particular, for resolved AGN, including the black holes at the center of the Milky Way (Sgr A*) and M81, nu_SA ranges from roughly 10 MHz to 10 GHz, and thus can be probed via existing and up-coming ground-based radio observatories.

preprint2008arXiv

Automated detection of galaxy-scale gravitational lenses in high resolution imaging data

Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]

preprint2008arXiv

Hot self-similar relativistic MHD flows

We consider axisymmetric relativistic jets with a toroidal magnetic field and an ultrarelativistic equation of state, with the goal of studying the lateral structure of jets whose pressure is matched to the pressure of the medium through which they propagate. We find all self-similar steady-state solutions of the relativistic MHD equations for this setup. One of the solutions is the case of a parabolic jet being accelerated by the pressure gradient as it propagates through a medium with pressure declining as p(z)\propto z^{-2}. As the jet material expands due to internal pressure gradients, it runs into the ambient medium resulting in a pile-up of material along the jet boundary, while the magnetic field acts to produce a magnetic pinch along the axis of the jet. Such jets can be in a lateral pressure equilibrium only if their opening angle θ_j at distance z is smaller than about 1/γ, where γis the characteristic bulk Lorentz-factor at this distance; otherwise, different parts of the jet cannot maintain causal contact. We construct maps of optically thin synchrotron emission from our self-similar models. We suggest that the boundary pile-up may be the reason for the limb-brightening of the sub-parsec jet of M87. We find that if the synchrotron emissivity falls with the distance from the jet axis, the polarization fraction rises toward the edge, as seen in 3C273 and Mkn501. Projection effects and the emissivity pattern of the jet have a strong effect on the observed polarization signal, so the interpretation of the polarization data in terms of the geometry of magnetic fields is rather uncertain. For example, jets with toroidal magnetic fields display the `spine-sheath' polarization angle pattern seen in some BL Lac objects.

preprint2008arXiv

What is the largest Einstein radius in the universe?

The Einstein radius plays a central role in lens studies as it characterises the strength of gravitational lensing. The distribution of Einstein radii near the upper cutoff should probe the largest mass concentrations in the universe. Adopting a triaxial halo model, we compute expected distributions of large Einstein radii. To assess the cosmic variance, we generate a number of all-sky Monte-Carlo realisations. We find that the expected largest Einstein radius in the universe is sensitive to the cosmological model: for a source redshift z=1, they are 42^{+9}_{-7}, 35^{+8}_{-6}, and 54^{+12}_{-7} arcseconds, assuming best-fit parameters of the WMAP5, WMAP3 and WMAP1 data, respectively. These values are broadly consistent with current observations given their incompleteness. For the same source redshift, we expect in all-sky 35 (WMAP5), 15 (WMAP3), and 150 (WMAP1) clusters that have Einstein radii larger than 20". Whilst the values of the largest Einstein radii are almost unaffected by the primordial non-Gaussianity currently of interest, the abundance of large lens clusters should probe non-Gaussianity competitively with CMB, but only if other cosmological parameters are well-measured. We also find that these "superlens" clusters constitute a highly biased population. For instance, a substantial fraction of these superlens clusters have major axes preferentially aligned with the line-of-sight. As a consequence, the projected mass distributions of the clusters are rounder by an ellipticity of 0.2 and have 40%-60% larger concentrations compared with typical clusters with similar redshifts and masses. We argue that the large concentration measured in A1689 is consistent with our model prediction at the 1.2σlevel. (Abridged)