We present evidence of anisotropic galaxy assembly bias in the Baryon Oscillation Spectroscopic Survey Data Release 12 galaxy sample at a level exceeding $5\sigma$. We use measurements of the line-of-sight velocity dispersion $\sigma_\star$ and stellar mass $M_\star$ to perform a simple split into subsamples of galaxies. We show that the amplitude of the monopole and quadrupole moments of the power spectrum depend differently on $\sigma_\star$ and $M_\star$, allowing us to split the galaxy sample into subsets with matching monopoles but significantly different quadrupoles on all scales. Combining data from the LOWZ and CMASS NGC galaxy samples, we find $>5\sigma$ evidence for anisotropic bias on scales $k<0.15\,h\,\rm{Mpc}^{-1}$. We also examine splits using other observed properties. For galaxy samples split using $M_\star$ and projected size $R_0$, we find no significant evidence of anisotropic bias. Galaxy samples selected using additional properties exhibit strongly varying degrees of anisotropic assembly bias, depending on which combination of properties is used to split into subsets. This may explain why previous searches for this effect using the Fundamental Plane found inconsistent results. We conclude that any selection of a galaxy sample that depends on $\sigma_\star$ can give biased and incorrect Redshift Space Distortion measurements.

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3{\sigma} constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6{\sigma} constraint within the same conditions.

We study the effect of large-scale tidal fields on internal halo properties using a set of N-body simulations. We measure significant cross-correlations between large-scale tidal fields and several non-scalar halo properties: shapes, velocity dispersion, and angular momentum. Selection effects that couple to these non-scalar halo properties can produce anisotropic clustering even in real-space. We investigate the size of this effect and show that it can produce a non-zero quadrupole similar in size to the one generated by linear redshift-space distortions (RSD). Finally, we investigate the clustering properties of halos identified in redshift-space and find enormous deviations from the standard linear RSD model, again caused by anisotropic assembly bias. These effects could contaminate the values of cosmological parameters inferred from the observed redshift-space clustering of galaxies, groups, or 21cm emission from atomic hydrogen, if their selection depends on properties affected by halo assembly bias. We briefly discuss ways in which this effect can be measured in existing and future large-scale structure surveys.

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.

We study the kinematics of galaxies within massive clusters, as a probe of the physics of star-formation quenching within clusters. Using N-body simulations, we argue that satellite kinematics provide information about galaxy infall that is complementary to the (instantaneous) spatial distribution of satellites. Comparing the simulation results with measurements of real cluster galaxies, we find evidence that the kinematics of red (quiescent) satellite galaxies are consistent with earlier infall times than that of blue (star-forming) satellites.

Self-Interacting Dark Matter (SIDM) has long been proposed as a solution to small scale problems posed by standard Cold Dark Matter (CDM). We use numerical simulations to study the effect of dark matter interactions on the morphology of disk galaxies falling into galaxy clusters. The effective drag force on dark matter leads to offsets of the stellar disk with respect to the surrounding halo, causing distortions in the disk. For anisotropic scattering cross-sections of 0.5 and 1.0$\,\textrm{cm}^{2}\textrm{g}^{-1}$, we show that potentially observable warps, asymmetries, and thickening of the disk occur in simulations. We discuss observational tests of SIDM with galaxy surveys and more realistic simulations needed to obtain detailed predictions.

The splashback radius is a physical scale in dark matter halos that is set by the gravitational dynamics of recently accreted shells. We use analytical models and N-body simulations to study the dependence of splashback on dark energy and screened modified gravity theories. In modified gravity models, the transition from screened to unscreened regions typically occurs in the cluster outskirts, suggesting potentially observable signatures in the splashback feature. We investigate the location of splashback in both chameleon and Vainshtein screened models and find significant differences compared with $\Lambda$CDM predictions. We also find an interesting interplay between dynamical friction and modified gravity, providing a distinctive signature for modified gravity models in the behavior of the splashback feature as a function of galaxy luminosity.

We study halo assembly bias for cluster-sized halos. Previous work has found little evidence for correlations between large-scale bias and halo mass assembly history for simulated cluster-sized halos, in contrast to the significant correlation found between bias and concentration for halos of this mass. This difference in behavior is surprising, given that both concentration and assembly history are closely related to the same properties of the linear-density peaks that collapse to form halos. Using publicly available simulations, we show that significant assembly bias is indeed found in the most massive halos with $M\sim 10^{15}M_\odot$, using essentially any definition of halo age. For lower halo masses $M\sim 10^{14}M_\odot$, no correlation is found between bias and the commonly used age indicator $a_{0.5}$, the half-mass time. We show that this is a mere accident, and that significant assembly bias exists for other definitions of halo age, including those based on the time when the halo progenitor acquires some fraction $f$ of the ultimate mass at $z=0$. For halos with $M_{\rm vir}\sim 10^{14}M_\odot$, the sense of assembly bias changes sign at $f=0.5$. We explore the origin of this behavior, and argue that it arises because standard definitions of halo mass in halo finders do not correspond to the collapsed, virialized mass that appears in the spherical collapse model used to predict large-scale clustering. Because bias depends strongly on halo mass, these errors in mass definition can masquerade as or even obscure the assembly bias that is physically present. More physically motivated halo definitions using splashback should be free of this particular defect of standard halo finders.

Cosmology is intrinsically intertwined with questions in fundamental physics. The existence of non-baryonic dark matter requires new physics beyond the Standard Model of elemenatary-particle interactions and Einstein's general relativity, as does the accelerating expansion of the universe. Current tensions between various cosmological measurements may be harbingers of yet more new physics. Progress on understanding dark matter and cosmic acceleration requires long term, high-precision measurements and excellent control of systematics, demanding observational programs that are often outside the discovery/characterization mode that drives many areas of astronomy. We outline potential programs through which the Hubble Space Telescope (HST) could have a major impact on issues in fundamental physics in the coming years. To realize this impact, we suggest the introduction of a "HST Fundamental Physics" observational program that would be subject to a modified proposal and review process.

(abridged) We investigate the signatures left by the cosmic neutrino background on the clustering of matter, CDM+baryons and halos in redshift-space using a set of more than 1000 N-body and hydrodynamical simulations with massless and massive neutrinos. We find that the effect neutrinos induce on the clustering of CDM+baryons in redshift-space on small scales is almost entirely due to the change in $\sigma_8$. Neutrinos imprint a characteristic signature in the quadrupole of the matter (CDM+baryons+neutrinos) field on small scales, that can be used to disentangle the effect of $\sigma_8$ and $M_\nu$. We show that the effect of neutrinos on the clustering of halos is very different, on all scales, to the one induced by $\sigma_8$. We find that the effects of neutrinos of the growth rate of CDM+baryons ranges from $\sim0.3\%$ to $2\%$ on scales $k\in[0.01, 0.5]~h{\rm Mpc}^{-1}$ for neutrinos with masses $M_\nu \leqslant 0.15$ eV. We compute the bias between the momentum of halos and the momentum of CDM+baryon and find it to be 1 on large scales for all models with massless and massive neutrinos considered. This point towards a velocity bias between halos and total matter on large scales that it is important to account for in order to extract unbiased neutrino information from velocity/momentum surveys such as kSZ observations. We show that baryonic effects can affect the clustering of matter and CDM+baryons in redshift-space by up to a few percent down to $k=0.5~h{\rm Mpc}^{-1}$. We find that hydrodynamics and astrophysical processes, as implemented in our simulations, only distort the relative effect that neutrinos induce on the anisotropic clustering of matter, CDM+baryons and halos in redshift-space by less than $1\%$. Thus, the effect of neutrinos in the fully non-linear regime can be written as a transfer function with very weak dependence on astrophysics.

Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries.

We analyze the ability of galaxy and CMB lensing surveys to constrain massive neutrinos and new models of dark radiation. We present a Fisher forecast analysis for neutrino mass constraints with the LSST galaxy survey and the CMB S4 survey. A joint analysis of the three galaxy and shear 2-point functions, along with key systematics parameters and Planck priors, constrains the neutrino masses to $\sum m_\nu = 0.041\,$eV at 1-$\sigma$ level, comparable to constraints expected from Stage 4 CMB lensing. If low redshift information from upcoming spectroscopic surveys like DESI is included, the constraint becomes $\sum m_\nu = 0.032\,$eV. These constraints are derived having marginalized over the number of relativistic species ($N_{\rm eff}$), which is somewhat degenerate with the neutrino mass. We also explore the gain by combining LSST and CMB S4, that is, using the five relevant auto- and cross-correlations of the two datasets. We conclude that advances in modeling the nonlinear regime and the measurements of other parameters are required to ensure a neutrino mass detection. Using the same datasets, we explore the ability of LSST-era surveys to test "nonstandard" models with dark radiation. We find that if evidence for dark radiation is found from $N_{\rm eff}$ measurements, the mass of the dark radiation candidate can be measured at a 1-$\sigma$ level of $0.162\,$eV for fermionic dark radiation, and $0.137\,$eV for bosonic dark radiation, for $\Delta N_{\rm eff} = 0.15$. We also find that the NNaturalness model of Arkani-Hamed et al 2016, with extra light degrees of freedom, has a sub-percent effect on the power spectrum: even more ambitious surveys than the ones considered here will be needed to test such models.

We present a new method for simulating cosmologies that contain massive particles with thermal free streaming motion, such as massive neutrinos or warm/hot dark matter. This method combines particle and fluid descriptions of the thermal species to eliminate the shot noise known to plague conventional N-body simulations. We describe this method in detail, along with results for a number of test cases to validate our method, and check its range of applicability. Using this method, we demonstrate that massive neutrinos can produce a significant scale-dependence in the large-scale biasing of deep voids in the matter field. We show that this scale-dependence may be quantitatively understood using an extremely simple spherical expansion model which reproduces the behavior of the void bias for different neutrino parameters.

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 10<lambda<20, and find that bright (M_i<-21.5) satellites have smaller splashback radii than fainter (M_i>-20) satellites at 99% confidence.

When dark matter halos are accreted by massive host clusters, strong gravitational tidal forces begin stripping mass from the accreted subhalos. This stripping eventually removes all mass beyond a subhalo's tidal radius, but the unbound mass remains in the vicinity of the satellite for at least a dynamical time t_dyn. The N-body subhalo study of Chamberlain et al. verified this picture and pointed out a useful observational consequence: measurements of subhalo correlations beyond the tidal radius are sensitive to the infall time, t_infall, of the subhalo onto its host. We perform this cross-correlation measurement using ~ 160,000 red satellite galaxies in SDSS redMaPPer clusters and find evidence that subhalo correlations do persist well beyond the tidal radius, suggesting that many of the observed satellites fell into their current host less than a dynamical time ago, t_infall < t_dyn. Combined with estimated dynamical times t_dyn ~ 3-5 Gyr and SED fitting results for the time at which satellites stopped forming stars, t_quench ~ 6 Gyr, we infer that for a significant fraction of the satellites, star formation quenched before those satellites entered their current hosts. The result holds for red satellites over a large range of cluster-centric distances 0.1 - 0.6 Mpc/h. We discuss the implications of this result for models of galaxy formation.

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

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\sigma$ 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 $\Lambda$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.

New observations of Jupiter's decametric radio emissions have been made with the Long Wavelength Array Station 1 (LWA1) which is capable of making high quality observations as low as 11 MHz. Full Stokes parameters were determined for bandwidths of 16 MHz. Here we present the first LWA1 results for the study of six Io-related events at temporal resolutions as fine as 0.25 ms. LWA1 data show excellent spectral detail in Jovian DAM such as simultaneous left hand circular (LHC) and right hand circular (RHC) polarized Io-related arcs and source envelopes, modulation lane features, S-bursts structures, narrow band N-events, and interactions between S-bursts and N-events. The sensitivity of the LWA1 combined with the low radio frequency interference environment allow us to trace the start of the LHC Io-C source region to much earlier CMLIII than typically found in the literature. We find the Io-C starts as early as CMLIII = 230 degrees at frequencies near 11 MHz. This early start of the Io-C emission may be valuable for refining models of the emission mechanism. We also detect modulation lane structures that appear continuous across LHC and RHC emissions, suggesting that both polarizations may originate from the same hemisphere of Jupiter. We present a study of rare S-bursts detected during an Io-D event and show drift rates are consistent with those from other Io-related sources. Finally, S-N burst events are seen in high spectral and temporal resolution and our data strongly support the co-spatial origins of these events.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

In recent years, it has become possible to detect individual dark matter subhalos near images of strongly lensed extended background galaxies. Typically, only the most massive subhalos in the strong lensing region may be detected this way. In this work, we show that strong lenses may also be used to constrain the much more numerous population of lower mass subhalos that are too small to be detected individually. In particular, we show that the power spectrum of projected density fluctuations in galaxy halos can be measured using strong gravitational lensing. We develop the mathematical framework of power spectrum estimation, and test our method on mock observations. We use our results to determine the types of observations required to measure the substructure power spectrum with high significance. We predict that deep observations ($\sim10$ hours on a single target) with current facilities can measure this power spectrum at the $3\sigma$ level, with no apparent degeneracy with unknown clumpiness in the background source structure or fluctuations from detector noise. Upcoming ALMA measurements of strong lenses are capable of placing strong constraints on the abundance of dark matter subhalos and the underlying particle nature of dark matter.

Recent work has shown that density profiles in the outskirts of dark matter halos can become extremely steep over a narrow range of radius. This behavior is produced by splashback material on its first apocentric passage after accretion. We show that the location of this splashback feature may be understood quite simply, from first principles. We present a simple model, based on spherical collapse, that accurately predicts the location of splashback without any free parameters. The important quantities that determine the splashback radius are accretion rate and redshift.

We investigate optimization strategies to measure primordial non-Gaussianity with future spectroscopic surveys. We forecast measurements coming from the 3D galaxy power spectrum and compute constraints on primordial non-Gaussianity parameters f_NL and n_NG. After studying the dependence on those parameters upon survey specifications such as redshift range, area, number density, we assume a reference mock survey and investigate the trade-off between number density and area surveyed. We then define the observational requirements to reach the detection of f_NL of order 1. Our results show that while power spectrum constraints on non-Gaussianity from future spectroscopic surveys can be competitive with current CMB limits, measurements from higher-order statistics will be needed to reach a sub unity precision in the measurements of the non-Gaussianity parameter f_NL.

This document presents the results from the Distances subgroup of the Cosmic Frontier Community Planning Study (Snowmass 2013). We summarize the current state of the field as well as future prospects and challenges. In addition to the established probes using Type IA supernovae and baryon acoustic oscillations, we also consider prospective methods based on clusters, active galactic nuclei, gravitational wave sirens and strong lensing time delays.

Satellites within simulated massive clusters are significantly spatially correlated with each other, even when those satellites are not gravitationally bound to each other. This correlation is produced by satellites that entered their hosts relatively recently, and is undetectable for satellites that have resided in their hosts for multiple dynamical timescales. Therefore, a measurement of clustering statistics of cluster satellites may be used to determine the typical accretion redshifts of those satellites into their observed hosts. We argue that such measurements may be used to determine the fraction of satellite galaxies that were quenched by their current hosts, thereby discriminating among models for quenching of star formation in satellite galaxies.

We propose a method to measure the ellipticities of dark matter halos using the lens-shear-shear 3-point correlation function. This method is immune to effects of galaxy-halo misalignments that can potentially limit 2-point galaxy-galaxy lensing measurements of halo anisotropy. Using a simple model for the projected mass distributions of dark matter halos, we construct an ellipticity estimator that sums over all possible triangular configurations of the 3-point function. By applying our estimator to halos from N-body simulations, we find that systematic errors in the recovered ellipticity will be at the <5% fractional level. We estimate that future imaging surveys like LSST will have sufficient statistics to detect halo ellipticities using 3-point lensing.

We perform cosmological N-body simulations with non-Gaussian initial conditions generated from two independent fields. The dominant contribution to the perturbations comes from a purely Gaussian field, but we allow the second field to have local non-Gaussianity that need not be weak. This scenario allows us to adjust the relative importance of non-Gaussian contributions beyond the skewness, producing a scaling of the higher moments different from (and stronger than) the scaling in the usual single field local ansatz. We compare semi-analytic prescriptions for the non-Gaussian mass function, large scale halo bias, and stochastic bias against the simulation results. We discuss applications of this work to large scale structure measurements that can test a wider range of models for the primordial fluctuations than is usually explored.

Recent observations have accumulated compelling evidence that some short gamma-ray bursts (SGRBs) are associated with the mergers of neutron star (NS) binaries. This would indicate that the SGRB event is associated with a gravitational-wave (GW) signal corresponding to the final inspiral of the compact binary. In addition, the radioactive decay of elements produced in NS binary mergers may result in transients visible in the optical and infrared with peak luminosities on hours-days timescales. Simultaneous observations of the inspiral GWs and signatures in the electromagnetic band may allow us to directly and independently determine both the luminosity distance and redshift to a binary. These standard sirens (the GW analog of standard candles) have the potential to provide an accurate measurement of the low-redshift Hubble flow. In addition, these systems are absolutely calibrated by general relativity, and therefore do not experience the same set of astrophysical systematics found in traditional standard candles, nor do the measurements rely on a distance ladder. We show that 15 observable GW and EM events should allow the Hubble constant to be measured with 5% precision using a network of detectors that includes advanced LIGO and Virgo. Measuring 30 beamed GW-SGRB events could constrain H_0 to better than 1%. When comparing to standard Gaussian likelihood analysis, we find that each event's non-Gaussian posterior in H_0 helps reduce the overall measurement errors in H_0 for an ensemble of NS binary mergers.

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

The first station of the Long Wavelength Array (LWA1) was completed in April 2011 and is currently performing observations resulting from its first call for proposals in addition to a continuing program of commissioning and characterization observations. The instrument consists of 258 dual-polarization dipoles, which are digitized and combined into beams. Four independently-steerable dual-polarization beams are available, each with two "tunings" of 16 MHz bandwidth that can be independently tuned to any frequency between 10 MHz and 88 MHz. The system equivalent flux density for zenith pointing is ~3 kJy and is approximately independent of frequency; this corresponds to a sensitivity of ~5 Jy/beam (5sigma, 1 s); making it one of the most sensitive meter-wavelength radio telescopes. LWA1 also has two "transient buffer" modes which allow coherent recording from all dipoles simultaneously, providing instantaneous all-sky field of view. LWA1 provides versatile and unique new capabilities for Galactic science, pulsar science, solar and planetary science, space weather, cosmology, and searches for astrophysical transients. Results from LWA1 will detect or tightly constrain the presence of hot Jupiters within 50 parsecs of Earth. LWA1 will provide excellent resolution in frequency and in time to examine phenomena such as solar bursts, and pulsars over a 4:1 frequency range that includes the poorly understood turnover and steep-spectrum regimes. Observations to date have proven LWA1's potential for pulsar observing, and just a few seconds with the completed 256-dipole LWA1 provide the most sensitive images of the sky at 23 MHz obtained yet. We are operating LWA1 as an open skies radio observatory, offering ~2000 beam-hours per year to the general community.

The inspirals and mergers of compact binaries are among the most promising events for ground-based gravitational-wave (GW) observatories. The detection of electromagnetic (EM) signals from these sources would provide complementary information to the GW signal. It is therefore important to determine the ability of gravitational-wave detectors to localize compact binaries on the sky, so that they can be matched to their EM counterparts. We use Markov Chain Monte Carlo techniques to study sky localization using networks of ground-based interferometers. Using a coherent-network analysis, we find that the Laser Interferometer Gravitational Wave Observatory (LIGO)-Virgo network can localize 50% of their ~8 sigma detected neutron star binaries to better than 50 sq.deg. with 95% confidence region. The addition of the Large Scale Cryogenic Gravitational Wave Telescope (LCGT) and LIGO-Australia improves this to 12 sq.deg.. Using a more conservative coincident detection threshold, we find that 50% of detected neutron star binaries are localized to 13 sq.deg. using the LIGO-Virgo network, and to 3 sq.deg. using the LIGO-Virgo-LCGT-LIGO-Australia network. Our findings suggest that the coordination of GW observatories and EM facilities offers great promise.

We investigate the collapse and internal structure of dark matter halos. We consider halo formation from initially scale-free perturbations, for which gravitational collapse is self-similar. Fillmore and Goldreich (1984) and Bertschinger (1985) solved the one dimensional (i.e. spherically symmetric) case. We generalize their results by formulating the three dimensional self-similar equations. We solve the equations numerically and analyze the similarity solutions in detail, focusing on the internal density profiles of the collapsed halos. By decomposing the total density into subprofiles of particles that collapse coevally, we identify two effects as the main determinants of the internal density structure of halos: adiabatic contraction and the shape of a subprofile shortly after collapse; the latter largely reflects the triaxiality of the subprofile. We develop a simple model that describes the results of our 3D simulations. In a companion paper, we apply this model to more realistic cosmological fluctuations, and thereby explain the origin of the nearly universal (NFW-like) density profiles found in N-body simulations.

We investigate the effect of supersonic relative velocities between baryons and dark matter, recently shown to arise generically at high redshift, on baryonic acoustic oscillation (BAO) measurements at low redshift. The amplitude of the relative velocity effect at low redshift is model-dependent, but can be parameterized by using an unknown bias. We find that if unaccounted, the relative velocity effect can shift the BAO peak position and bias estimates of the dark energy equation-of-state due to its non-smooth, out-of-phase oscillation structure around the BAO scale. Fortunately, the relative velocity effect can be easily modeled in constraining cosmological parameters without substantially inflating the error budget. We also demonstrate that the presence of the relative velocity effect gives rise to a unique signature in the galaxy bispectrum, which can be utilized to isolate this effect. Future dark energy surveys can accurately measure the relative velocity effect and subtract it from the power spectrum analysis to constrain dark energy models with high precision.

The friends-of-friends algorithm (hereafter, FOF) is a percolation algorithm which is routinely used to identify dark matter halos from N-body simulations. We use results from percolation theory to show that the boundary of FOF halos does not correspond to a single density threshold but to a range of densities close to a critical value that depends upon the linking length parameter, b. We show that for the commonly used choice of b = 0.2, this critical density is equal to 81.62 times the mean matter density. Consequently, halos identified by the FOF algorithm enclose an average overdensity which depends on their density profile (concentration) and therefore changes with halo mass contrary to the popular belief that the average overdensity is ~180. We derive an analytical expression for the overdensity as a function of the linking length parameter b and the concentration of the halo. Results of tests carried out using simulated and actual FOF halos identified in cosmological simulations show excellent agreement with our analytical prediction. We also find that the mass of the halo that the FOF algorithm selects crucially depends upon mass resolution. We find a percolation theory motivated formula that is able to accurately correct for the dependence on number of particles for the mock realizations of spherical and triaxial Navarro-Frenk-White halos. However, we show that this correction breaks down when applied to the real cosmological FOF halos due to presence of substructures. Given that abundance of substructure depends on redshift and cosmology, we expect that the resolution effects due to substructure on the FOF mass and halo mass function will also depend on redshift and cosmology and will be difficult to correct for in general. Finally, we discuss the implications of our results for the universality of the mass function.

Motivated by the properties of early universe scenarios that produce observationally large local non-Gaussianity, we perform N-body simulations with non-Gaussian initial conditions from a generalized local ansatz. The bispectra are schematically of the local shape, but with scale-dependent amplitude. We find that in such cases the size of the non-Gaussian correction to the bias of small and large mass objects depends on the amplitude of non-Gaussianity roughly on the scale of the object. In addition, some forms of the generalized bispectrum alter the scale dependence of the non-Gaussian term in the bias by a fractional power of k. These features may allow significant observational constraints on the particle physics origin of any observed local non-Gaussianity, distinguishing between scenarios where a single field or multiple fields contribute to the curvature fluctuations. While analytic predictions for the non-Gaussian bias agree qualitatively with the simulations, we find numerically a stronger observational signal than expected. This suggests that a more precise understanding of halo formation is needed to fully explain the consequences of primordial non-Gaussianity

This paper presents a search for radio transients at a frequency of 73.8 MHz (4 m wavelength) using the all-sky imaging capabilities of the Long Wavelength Demonstrator Array (LWDA). The LWDA was a 16-dipole phased array telescope, located on the site of the Very Large Array in New Mexico. The field of view of the individual dipoles was essentially the entire sky, and the number of dipoles was sufficiently small that a simple software correlator could be used to make all-sky images. From 2006 October to 2007 February, we conducted an all-sky transient search program, acquiring a total of 106 hr of data; the time sampling varied, being 5 minutes at the start of the program and improving to 2 minutes by the end of the program. We were able to detect solar flares, and in a special-purpose mode, radio reflections from ionized meteor trails during the 2006 Leonid meteor shower. We detected no transients originating outside of the solar system above a flux density limit of 500 Jy, equivalent to a limit of no more than about 10^{-2} events/yr/deg^2, having a pulse energy density >~ 1.5 x 10^{-20} J/m^2/Hz at 73.8 MHz for pulse widths of about 300 s. This event rate is comparable to that determined from previous all-sky transient searches, but at a lower frequency than most previous all-sky searches. We believe that the LWDA illustrates how an all-sky imaging mode could be a useful operational model for low-frequency instruments such as the Low Frequency Array, the Long Wavelength Array station, the low-frequency component of the Square Kilometre Array, and potentially the Lunar Radio Array.

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics.

Recent observations support the hypothesis that a large fraction of "short-hard" gamma-ray bursts (SHBs) are associated with compact binary inspiral. Since gravitational-wave (GW) measurements of well-localized inspiraling binaries can measure absolute source distances, simultaneous observation of a binary's GWs and SHB would allow us to independently determine both its luminosity distance and redshift. Such a "standard siren" (the GW analog of a standard candle) would provide an excellent probe of the relatively nearby universe's expansion, complementing other standard candles. In this paper, we examine binary measurement using a Markov Chain Monte Carlo technique to build the probability distributions describing measured parameters. We assume that each SHB observation gives both sky position and the time of coalescence, and we take both binary neutron stars and black hole-neutron star coalescences as plausible SHB progenitors. We examine how well parameters particularly distance) can be measured from GW observations of SHBs by a range of ground-based detector networks. We find that earlier estimates overstate how well distances can be measured, even at fairly large signal-to-noise ratio. The fundamental limitation to determining distance proves to be a degeneracy between distance and source inclination. Overcoming this limitation requires that we either break this degeneracy, or measure enough sources to broadly sample the inclination distribution. (Abridged)

The apparent presence of large core radii in Low Surface Brightness galaxies has been claimed as evidence in favor of warm dark matter. Here we show that WDM halos do not have cores that are large fractions of the halo size: typically, r_core/r_200 < 0.001. This suggests an astrophysical origin for the large cores observed in these galaxies, as has been argued by other authors.

Recent work has shown that at high redshift, the relative velocity between dark matter and baryonic gas is typically supersonic. This relative velocity suppresses the formation of the earliest baryonic structures like minihalos, and the suppression is modulated on large scales. This effect imprints a characteristic shape in the clustering power spectrum of the earliest structures, with significant power on 100 Mpc scales featuring highly pronounced baryon acoustic oscillations. The amplitude of these oscillations is orders of magnitude larger at z=20 than previously expected. This characteristic signature can allow us to distinguish the effects of minihalos on intergalactic gas at times preceding and during reionization. We illustrate this effect with the example of 21 cm emission and absorption from redshifts during and before reionization. This effect can potentially allow us to probe physics on kpc scales using observations on 100 Mpc scales. We present sensitivity forecasts for FAST and Arecibo. Depending on parameters, this enhanced structure may be detectable by Arecibo at redshifts near z=15-20, and with appropriate instrumentation FAST could measure the BAO power spectrum with high precision. In principle, this effect could also pose a serious challenge for efforts to constrain dark energy using observations of the BAO feature at low redshift.

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy.

We present Hubble Space Telescope (HST) WFPC2 photometry of 13 microlensed source stars from the 5.7 year Large Magellanic Cloud (LMC) survey conducted by the MACHO Project. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. None of these sources is coincident with a background galaxy, which rules out the possibility that the MACHO LMC microlensing sample is contaminated with misidentified supernovae or AGN in galaxies behind the LMC. This supports the conclusion that the MACHO LMC microlensing sample has only a small amount of contamination due to non-microlensing forms of variability. We compare the WFPC2 source star magnitudes with the lensed flux predictions derived from microlensing fits to the light curve data. In most cases the source star brightness is accurately predicted. Finally, we develop a statistic which constrains the location of the Large Magellanic Cloud (LMC) microlensing source stars with respect to the distributions of stars and dust in the LMC and compare this to the predictions of various models of LMC microlensing. This test excludes at > 90% confidence level models where more than 80% of the source stars lie behind the LMC. Exotic models that attempt to explain the excess LMC microlensing optical depth seen by MACHO with a population of background sources are disfavored or excluded by this test. Models in which most of the lenses reside in a halo or spheroid distribution associated with either the Milky Way or the LMC are consistent which these data, but LMC halo or spheroid models are favored by the combined MACHO and EROS microlensing results.

We study the effect of primordial nongaussianity on large-scale structure, focusing upon the most massive virialized objects. Using analytic arguments and N-body simulations, we calculate the mass function and clustering of dark matter halos across a range of redshifts and levels of nongaussianity. We propose a simple fitting function for the mass function valid across the entire range of our simulations. We find pronounced effects of nongaussianity on the clustering of dark matter halos, leading to strongly scale-dependent bias. This suggests that the large-scale clustering of rare objects may provide a sensitive probe of primordial nongaussianity. We very roughly estimate that upcoming surveys can constrain nongaussianity at the level |fNL| <~ 10, competitive with forecasted constraints from the microwave background.

We investigate the origin of halo assembly bias, the dependence of halo clustering on assembly history. We relate halo assembly to peak properties measured in the Lagrangian space of the initial linear Gaussian random density field, and show how these same Lagrangian properties determine large-scale bias. We focus on the two regimes where assembly bias has been observed to be significant: at masses very large and very small compared to the nonlinear mass scale. At high masses, we show that assembly bias is expected from the statistics of the peaks of Gaussian random fluctuations, and we show that the extent of assembly bias found in N-body simulations of rare halos is in excellent agreement with our theoretical prediction. At low masses, we argue that assembly bias largely arises from a sub-population of low mass halos whose mass accretion has ceased. Due to their arrested development, these halos naturally become unbiased, in contrast to their anti-biased peers. We show that a simple toy model incorporating these effects can roughly reproduce the bias trends found in N-body simulations.

In recent years, gravitational lensing has been used as a means to detect substructure in galaxy-sized halos, via anomalous flux ratios in quadruply-imaged lenses. In addition to causing anomalous flux ratios, substructure may also perturb the positions of lensed images at observable levels. In this paper, we numerically investigate the scale of such astrometric perturbations using realistic models of substructure distributions. Substructure distributions that project clumps near the Einstein radius of the lens result in perturbations that are the least degenerate with the best-fit smooth macromodel, with residuals at the milliarcsecond scale. Degeneracies between the center of the lens potential and astrometric perturbations suggest that milliarcsecond constraints on the center of the lensing potential boost the observed astrometric perturbations by an order of magnitude compared to leaving the center of the lens as a free parameter. In addition, we discuss methods of substructure detection via astrometric perturbations that avoid full lens modeling in favor of local image observables and also discuss modeling of systems with luminous satellites to constrain the masses of those satellites.

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

Observations of the gravitational radiation from well-localized, inspiraling compact object binaries can measure absolute source distances with high accuracy. When coupled with an independent determination of redshift through an electromagnetic counterpart, these standard sirens can provide an excellent probe of the expansion history of the Universe and the dark energy. Short gamma-ray bursts, if produced by merging neutron star binaries, would be standard sirens with known redshifts detectable by ground-based GW networks such as LIGO-II, Virgo, and AIGO. Depending upon the collimation of these GRBs, a single year of observation of their gravitational waves can measure the Hubble constant to about 2%. When combined with measurement of the absolute distance to the last scattering surface of the cosmic microwave background, this determines the dark energy equation of state parameter w to 9%. Similarly, supermassive binary black hole inspirals will be standard sirens detectable by LISA. Depending upon the precise redshift distribution, 100 sources could measure w at the 4% level.

We study the impact of lens galaxy environments on the image separation distribution of lensed quasars. We account for both environmental convergence and shear, using a joint distribution derived from galaxy formation models calibrated by galaxy-galaxy lensing data and number counts of massive elliptical galaxies. We find that the external field enhances lensing probabilities, particularly at large image separations; the increase is ~30% at \theta=3'' and ~200% at \theta=5'', when we adopt a power-law source luminosity function \Phi(L) \propto L^-2.1. The enhancement is mainly driven by convergence, which boosts both the image separation and magnification bias (for a fixed lens galaxy mass). These effects have been neglected in previous studies of lens statistics. Turning the problem around, we derive the posterior convergence and shear distributions and point out that they are strong functions of image separation; lens systems with larger image separations are more likely to lie in dense environments.

We compute the expected number of quasars multiply imaged by cluster size dark halos for current wide field quasar surveys by carrying out a large ensemble of ray tracing simulations through clusters from a cosmological N-body simulation of the LCDM cosmology. Our calculation predicts ~ 4 quasar lenses with splittings theta > 10" in the SDSS spectroscopic quasar sample, consistent with the recent discovery of the wide separation lens SDSSJ1004+4112 which has theta=14.6". The SDSS faint photometric quasar survey will contain ~12 multiply imaged quasars with splittings theta > 10". Of these, ~ 2 will be lenses with separations theta > 30", and ~ 2 will be at high redshift (z ~ 4).

We present a detailed investigation into which properties of CDM halos make them effective strong gravitational lenses. Strong lensing cross sections of 878 clusters from an N-body simulation are measured by ray tracing through 13,594 unique projections. We measure concentrations, axis ratios, orientations, and the amount of substructure of each cluster, and compare the lensing weighted distribution of each quantity to that of the cluster population as a whole. The concentrations of lensing clusters are on average 34% larger than the typical cluster in the Universe. Despite this bias, the anomalously high concentrations (c >14) recently measured by several groups, appear to be inconsistent with the concentration distribution in our simulations, which predict < 2% of lensing clusters should have concentrations this high. No correlation is found between lensing cross section and the amount of substructure. We introduce several types of simplified dark matter halos, and use them to isolate which properties of CDM clusters make them effective lenses. Projections of halo substructure onto small radii and the large scale mass distribution of clusters do not significantly influence cross sections. The abundance of giant arcs is primarily determined by the mass distribution within an average overdensity of ~ 10,000. A multiple lens plane ray tracing algorithm is used to show that projections of large scale structure increase the giant arc abundance by a modest amount <7%. We revisit the question of whether there is an excess of giant arcs behind high redshift clusters in the RCS survey and find that the number of high redshift (z > 0.6) lenses is in good agreement with LCDM, although our simulations predict more low redshift (z < 0.6) lenses than were observed. (abridged)

Using measured tangential shear profiles and number counts of massive elliptical galaxies, the halo occupation distribution of strong lensing galaxies is constrained. The resulting HOD is then used to populate an N-body simulation with lens galaxies, in order to assess the importance of environment for strong lensing systems. Typical estimated values for the convergence and shear produced by nearby correlated matter are kappa = gamma = 0.03, with much stronger events occurring relatively infrequently. This implies that estimates of quantities like the Hubble constant are not expected to be significantly biased by environmental effects. One puzzle is that predicted values for the external shear at lens galaxies are far below the values obtained by modeling of strong lensing data.

Giant arcs in strong lensing galaxy clusters can provide a purely geometric determination of cosmological parameters, such as the dark energy density and equation of state. We investigate sources of noise in cosmography with giant arcs, focusing in particular on errors induced by density fluctuations along the line-of-sight, and errors caused by modeling uncertainties. We estimate parameter errors in two independent ways, first by developing a Fisher matrix formalism for strong lensing parameters, and next by directly ray-tracing through N-body simulations using a multi-plane lensing code. We show that for reasonable power spectra, density fluctuations from large-scale structure produce > 100% errors in cosmological parameters derived from any single sightline, precluding the use of individual clusters or golden lenses to derive accurate cosmological constraints. Modeling uncertainties similarly can lead to large errors, and we show that the use of parametrized mass models in fitting strong lensing clusters can significantly bias the inferred cosmological parameters. We lastly speculate on means by which these errors may be corrected.

We cross-correlate the sample of Type Ia supernovae from Riess et al. (2004) with the SDSS DR2 photometric galaxy catalogue. In contrast to recent work, we find no detectable correlation between supernova magnitude and galaxy overdensity on scales ranging between 1 and 10 arcminutes. Our results are in accord with theoretical expectations for gravitational lensing of supernovae by large-scale structure. Future supernova surveys like SNAP will be capable of detecting unambiguously the predicted lensing signal.

Using stellar dynamics and strong gravitational lensing as complementary probes, Sand et al. (2002, 2003) have recently claimed strong evidence for shallow dark matter density profiles in several lensing clusters, which may conflict with predictions of the Cold Dark Matter paradigm. However, systematic uncertainties in the analysis weaken the constraints. By re-analyzing their data, we argue that the tight constraints claimed by Sand et al., were driven by prior assumptions. Relaxing the assumptions, we find that no strong constraints may be derived on the dark matter inner profile from the Sand et al. data; we find satisfactory fits (with reasonable parameters) for a wide range of inner slopes 0 < beta < 1.4. Useful constraints on the mass distributions of lensing clusters can still be obtained, but they require moving beyond mere measurements of lensing critical radii into the realm of detailed lens modeling.

We study the expected properties and statistics of giant arcs produced by galaxy clusters in a LambdaCDM universe and investigate how the characteristics of CDM clusters determine the properties of the arcs they generate. Due to the triaxiality and substructure of CDM halos, the giant arc cross section for individual clusters varies by more than an order of magnitude as a function of viewing angle. In addition, the shallow density cusps and triaxiality of CDM clusters cause systematic alignments of giant arcs which should be testable with larger samples from forthcoming lensing surveys. We compute the predicted statistics of giant arcs for the LambdaCDM model and compare to results from previous surveys. The predicted arc statistics are in excellent agreement with the numbers of giant arcs observed around low redshift (0.2 < z < 0.6) clusters from the EMSS sample, however there are hints of a possible excess of arcs observed around high redshift z > 0.6 clusters. This excess, if real, appears to be due to the presence of highly massive or concentrated clusters at high redshifts.

The flux anomalies in four-image gravitational lenses can be interpreted as evidence for the dark matter substructure predicted by cold dark matter (CDM) halo models. In principle, these flux anomalies could arise from alternate sources such as absorption, scattering or scintillation by the interstellar medium (ISM) of the lens galaxy, problems in the ellipsoidal macro models used to fit lens systems, or stellar microlensing. We apply several tests to the data that appear to rule out these alternate explanations. First, the radio flux anomalies show no significant dependence on wavelength, as would be expected for almost any propagation effect in the ISM or microlensing by the stars. Second, the flux anomaly distributions show the characteristic demagnifications of the brightest saddle point relative to the other images expected for low optical depth substructure, which cannot be mimicked by either the ISM or problems in the macro models. Microlensing by stars also cannot reproduce the suppression of the bright saddle points if the radio source sizes are consistent with the Compton limit for their angular sizes. Third, while it is possible to change the smooth lens models to fit the flux anomalies in some systems, we can rule out the necessary changes in all systems where we have additional lens constraints to check the models. Moreover, the parameters of these models are inconsistent with our present observations and expectations for the structure of galaxies. We conclude that low-mass halos remain the best explanation of the phenomenon.

We use a simple statistical test to show that the anomalous flux ratios observed in gravitational lenses are created by gravitational perturbations from substructure rather than propagation effects in the interstellar medium or incomplete models for the gravitational potential of the lens galaxy. We review current estimates that the substructure represents between 0.6% and 7% (90% confidence) of the lens galaxy mass, and outline future observational programs which can improve the results.

Interferometers offer multiple methods for studying microlensing events and determining the properties of the lenses. We investigate the study of microlensing events with optical interferometers, focusing on narrow-angle astrometry, visibility, and closure phase. After introducing the basics of microlensing and interferometry, we derive expressions for the signals in each of these three channels. For various forecasts of the instrumental performance, we discuss which method provides the best means of measuring the lens angular Einstein radius theta_E, a prerequisite for determining the lens mass. If the upcoming generation of large-aperture, AO-corrected long baseline interferometers (e.g. VLTI, Keck, OHANA) perform as well as expected, theta_E may be determined with signal-to-noise greater than 10 for all bright events. We estimate that roughly a dozen events per year will be sufficiciently bright and have long enough durations to allow the measurement of the lens mass and distance from the ground. We also consider the prospects for a VLTI survey of all bright lensing events using a Fisher matrix analysis, and find that even without individual masses, interesting constraints may be placed on the bulge mass function, although large numbers of events would be required.

Weak lensing of high-redshift Type Ia supernovae induces an external dispersion in their observed standard candle brightnesses, comparable in magnitude to the intrinsic dispersion for redshifts z>1. The same matter fluctuations responsible for the magnification of distant supernovae also generate shear in the images of background galaxies. We investigate the possibility of using lensing shear maps constructed from galaxies surrounding the supernovae as a means of correcting the lensing-induced magnification dispersion. We find that a considerable fraction of the lensing dispersion derives from sub-arcminute scales, which are not probed by shear maps smoothed on arcminute scales. We thus find that weak lensing shear maps will be of only limited value in reducing the weak lensing magnification fluctuations of supernovae.

We devise a method to measure the abundance of satellite halos in gravitational lens galaxies, and apply our method to a sample of 7 lens systems. After using Monte Carlo simulations to verify the method, we find that substructure comprises fraction f=0.02 (median, 0.006<f<0.07 at 90% confidence) of the mass of typical lens galaxies, in excellent agreement with predictions of CDM simulations. We estimate a characteristic critical radius for the satellites of 0.0001<b/arcsec<0.006 (90% confidence). For a satellite mass function of dn/dM M^x with x=-1.8 and M_l<M<M_h, the critical radius provides an estimate that the upper mass limit is 10^6Msun < M_h < 10^9Msun. Our measurement confirms a generic prediction of CDM models, and may obviate the need to invoke alternatives to CDM like warm dark matter or self-interacting dark matter.

We place limits on the linear power spectrum on small scales (k > 50 h/Mpc) using measurements of substructure in gravitational lens galaxies. We find excellent agreement with the simplest LambdaCDM models, and in conjunction with other cosmological probes, place constraints on the neutrino mass m_nu, tilt of the primordial power spectrum n, and mass of the dark matter particle m. We find n>0.94, and for a Harrison-Zeldovich spectrum, find m_nu<0.74 eV and m>5.2 keV, at 95% confidence.

The properties of multiple image gravitational lenses require a fractional surface mass density in satellites of f=0.02 (0.006 < f < 0.07 at 90% confidence) that is consistent with the expectations for CDM. The characteristic satellite mass scale, 10^6-10^9 Msun, is also consistent with the expectations for CDM. The agreement between the observed and expected density of CDM substructure shows that most low mass galactic satellites fail to form stars, and this absence of star formation explains the discrepancy between the number of observed Galactic satellites and CDM predictions rather than any modification to the CDM theory such as self-interacting dark matter or a warm dark matter component.

Roughly half of the observed gamma-ray bursts (GRBs) may arise from the shocking of an ultra-relativistic shell of protons with the interstellar medium (ISM). Any neutrons originally present in the GRB fireball may, depending on the characteristics of the central engine, dynamically decouple as the fireball accelerates. This leads to outflow consisting of separate fast proton and slow neutron components. We derive detailed implications of neutron decoupling for the observed lightcurves of short bursts. We show that the collision of a neutron decayed shell with a decelerating outer shell is expected to result in an observable second peak in the GRB lightcurve. There may be substantial optical emission associated with such an event, so the upcoming Swift satellite may be able to place constraints on models for short bursts. We also discuss interesting inferences about central engine characteristics allowed by existing BATSE data and a consideration of neutron decoupling.

Dark energy models which alter the relative scaling behavior of dark energy and matter could provide a natural solution to the cosmic coincidence problem - why the densities of dark energy and dark matter are comparable today. A generalized class of dark energy models is introduced which allows non-canonical scaling of the ratio of dark matter and dark energy with the Robertson-Walker scale factor a(t). Upcoming observations, such as a high redshift supernova survey, application of the Alcock-Paczynski test to quasar pairs, and cluster evolution, will strongly constrain the relative scaling of dark matter and dark energy as well as the equation of state of the dark energy. Thus, whether there actually is a coincidence problem, and the extent of cosmic coincidence in the universe's recent past can be answered observationally in the near future. Determining whether today is a special time in the history of the universe will be a SNAP.

If gamma-ray burst (GRB) emission is strongly collimated then GRBs occur throughout the Universe at a rate much higher than is detected. Since the emission from the optical afterglow is thought to be more isotropic than the gamma-ray emission, it has been hypothesized that a search for orphan afterglows (those without the triggering GRB) would allow strong constraints to be placed on the degree of GRB collimation. We show here that, within the context of leading models of GRB jet evolution, measurement of the GRB beaming angle using optical orphan searches is extremely difficult, perhaps impossible in practice. This is because in the leading model of GRB jets, the effective afterglow beaming angle scales with the jet angle for small angles, and so the ratio of detected orphan afterglows to GRBs is independent of the jet opening angle. Thus, the number of expected afterglow detections is the same for moderate jet angles (e.g. 20 deg) as for arbitrarily small jet angles (<< 0.1 deg). For nearly isotropic GRB geometry, or for radio afterglow searches in which the jet has become non-relativistic, the ratio of afterglows to GRBs may give information on collimation. However, using a simple model we estimate the expected number of orphan detections in current supernova surveys, and find this number to be less than one, for all jet opening angles. Even for future supernova surveys, the small detection rate and lack of dependence on collimation angle appear to ruin the prospects of determining GRB beaming by this method. Radio searches may provide the best hope to find the missing orphans.

We report the discovery of eight new R Coronae Borealis (RCB) stars in the Large Magellanic Cloud (LMC) using the MACHO project photometry database. The discovery of these new stars increases the number of known RCB stars in the LMC to thirteen. We have also discovered four stars similar to the Galactic variable DY Per. These stars decline much more slowly and are cooler than the RCB stars. The absolute luminosities of the Galactic RCB stars are unknown since there is no direct measurement of the distance to any Galactic RCB star. Hence, the importance of the LMC RCB stars. We find a much larger range of absolute magnitudes (M(V) = -2.5 to -5 mag) than inferred from the small pre-MACHO sample of LMC RCB stars. It is likely that there is a temperature - M(V)relationship with the cooler stars being intrinsically fainter. Cool (~5000 K) RCB stars are much more common than previously thought based on the Galactic RCB star sample. Using the fairly complete sample of RCB stars discovered in the MACHO fields, we have estimated the likely number of RCB stars in the Galaxy to be ~3,200. The SMC MACHO fields were also searched for RCB stars but none were found.

The Space Interferometry Mission (SIM) will observe sources in crowded fields. Recent work has shown that source crowding can induce significant positional errors in SIM's astrometric measurements, even for targets many magnitudes brighter than all other crowding sources. Here we investigate whether the spectral decomposition of the fringe pattern may be used to disentangle the overlapping fringes from multiple blended sources, effectively by performing synthesis imaging with two baselines. We find that spectrally dispersed fringes enable SIM to identify and localize a limited number of field sources quite robustly, thereby removing their effect from SIM astrometry and reducing astrometry errors to near photon noise levels. We simulate SIM measurements of the LMC, and show that (a) SIM astrometry will not be corrupted by blending and (b) extremely precise imaging of mildly crowded fields may be performed using only two orthogonal baseline orientations, allowing microarcsecond positional measurements. We lastly illustrate the method's potential with the example of astrometric microlensing, showing that SIM's mass and distance measurements of lenses will be untainted by crowding.

We attempt to determine whether the MACHO microlensing source stars are drawn from the average population of the LMC or from a population behind the LMC by examining the HST color-magnitude diagram (CMD) of microlensing source stars. We present WFPC2 HST photometry of eight MACHO microlensing source stars and the surrounding fields in the LMC. The microlensing source stars are identified by deriving accurate centroids in the ground-based MACHO images using difference image analysis (DIA) and then transforming the DIA coordinates to the HST frame. We consider in detail a model for the background population of source stars based on that presented by Zhao, Graff & Guhathakurta. In this model, the source stars have an additional reddening <E(B-V)> = 0.13 mag and a slightly larger distance modulus <Delta mu> ~ 0.3 mag than the average LMC population. We also investigate a series of source star models, varying the relative fraction of source stars drawn from the average and background populations and the displacement of the background population from the LMC. Due to the small number of analyzed events the distribution of probabilities of different models is rather flat. A shallow maximum occurs at a fraction s_LMC ~ 0.8 of the source stars in the LMC. This is consistent with the interpretation that a significant fraction of observed microlensing events are due to lenses in the Milky Way halo, but does not definitively exclude other models.

We report on a search for long duration microlensing events towards the Large Magellanic Cloud. We find none, and therefore put limits on the contribution of high mass objects to the Galactic dark matter. At 95% confidence level we exclude objects in the mass range 0.3 solar masses to 30.0 solar masses from contributing more than 4 times 10^11 solar masses to the Galactic halo. Combined with earlier results, this means that objects with masses under 30 solar masses cannot make up the entire dark matter halo if the halo is of typical size. For a typical dark halo, objects with masses under 10 solar masses contribute less than 40% of the dark matter.

We investigate the so-called magnification relations of gravitational lensing models. We show that multidimensional residue integrals provide a simple explanation for the existence of these relations, and an effective method of computation. We illustrate the method with several examples, thereby deriving new magnification relations for galaxy lens models and microlensing (point mass lensing).

In light of recent evidence suggesting a nonzero present-day cosmological constant, Adams, Mbonye, & Laughlin (1999) have considered the evolution of black holes in the presence of vacuum energy. Using the assumption that Lambda remains constant with time, and a conjecture based on a paper by Mallett (1986), they reach the remarkable conclusion that black holes with current mass greater than about 2*10^(-9) Msun will not Hawking evaporate in the distant future, but will instead absorb vacuum energy and grow to roughly the de Sitter horizon size. In this letter we reexamine black hole evaporation in the presence of vacuum energy, and find instead that all known black holes will eventually evaporate.

We report on our search for microlensing towards the Large Magellanic Cloud (LMC). Analysis of 5.7 years of photometry on 11.9 million stars in the LMC reveals 13 - 17 microlensing events. This is significantly more than the $\sim$ 2 to 4 events expected from lensing by known stellar populations. The timescales ($\that$) of the events range from 34 to 230 days. We estimate the microlensing optical depth towards the LMC from events with $2 < \that < 400$ days to be 1.2 ^{+0.4}_ {-0.3} \ten{-7}$, with an additional 20% to 30% of systematic error. The spatial distribution of events is mildly inconsistent with LMC/LMC disk self-lensing, but is consistent with an extended lens distribution such as a Milky Way or LMC halo. Interpreted in the context of a Galactic dark matter halo, consisting partially of compact objects, a maximum likelihood analysis gives a MACHO halo fraction of 20% for a typical halo model with a 95% confidence interval of 8% to 50%. A 100% MACHO halo is ruled out at the 95% C.L. for all except our most extreme halo model. Interpreted as a Galactic halo population, the most likely MACHO mass is between $ 0.15 \msun$ and $ 0.9 \msun$, depending on the halo model, and the total mass in MACHOs out to 50 kpc is found to be 9+4-3 10^{10} msun, independent of the halo model. These results are marginally consistent with our previous results, but are lower by about a factor of two. Besides a larger data set, this work also includes an improved efficiency determination, improved likelihood analysis, and more thorough testing of systematic errors, especially with respect to the treatment of potential backgrounds to microlensing, such as supernovae in galaxies behind the LMC. [Abridged]

All of the proposed explanations for the microlensing events observed towards the LMC have difficulties. One of these proposed explanations, LMC self-lensing, which invokes ordinary LMC stars as the long sought-after lenses, has recently gained considerable popularity as a possible solution to the microlensing conundrum. In this paper, we carefully examine the set of LMC self-lensing models. In particular, we review the pertinent observations made of the LMC, and show how these observations place limits on such self-lensing models. We find that, given current observational constraints, no purely LMC disk models are capable of producing optical depths as large as that reported in the MACHO collaboration 2-year analysis. Besides pure disk, we also consider alternate geometries, and present a framework which encompasses the previous studies of LMC self-lensing. We discuss which model parameters need to be pushed in order for such models to succeed. For example, like previous workers, we find that an LMC halo geometry may be able to explain the observed events. However, since all known LMC tracer stellar populations exhibit disk-like kinematics, such models will have difficulty being reconciled with observations. For SMC self-lensing, we find predicted optical depths differing from previous results, but more than sufficient to explain all observed SMC microlensing. In contrast, for the LMC we find a self-lensing optical depth contribution between 0.47e-8 and 7.84e-8, with 2.44e-8 being the value for the set of LMC parameters most consistent with current observations.

We introduce a new method of searching for and characterizing extra-solar planets. We show that by monitoring the center-of-light motion of microlensing alerts using the next generation of high precision astrometric instruments the probability of detecting a planet orbiting the lens is high. We show that adding astrometric information to the photometric microlensing lightcurve greatly helps in determining the planetary mass and semi-major axis. We introduce astrometric maps as a new way for calculating astrometric motion and planet detection probabilities. Finite source effects are important for low mass planets, but even Earth mass planets can give detectable signals.

We demonstrate that for several of the gravitational lens models used to describe galaxies, there exists a quantity we dub the magnification invariant, equaling the sum of the signed magnifications of the images, that is a constant when the image multiplicity is a maximum. This invariant is independent of most of the model parameters, and is independent of the source position as long as the source lies inside the caustic. It is not necessary to solve the lens equations to compute this invariant. For quad lenses, it is usually easy to assign image parities and thus one can compute the sum of the signed fluxes, compare with the magnification invariant of different models, and infer the model magnification factors of the images without fitting. We evaluate the invariant for simple galaxy models, apply it to known cases of quadruple lenses, and discuss implications for the ability of these lens models to reproduce observed image brightnesses.