A joint analysis of data collected by the Planck and BICEP2+Keck teams has previously given $r = 0.09^{+0.06}_{-0.04}$ for BICEP2 and $r = 0.02^{+0.04}_{-0.02}$ for Keck. Analyzing BICEP2 using its published noise estimate, we had earlier (Colley & Gott 2015) found $r = 0.09 \pm 0.04$, agreeing with the final joint results for BICEP2. With the Keck data now available, we have done something the joint analysis did not: a correlation study of the BICEP2 vs. Keck B-mode maps. Knowing the correlation coefficient between the two and their amplitudes allows us to determine the noise in each map (which we check using the E-modes). We find the noise power in the BICEP2 map to be twice the original BICEP2 published estimate, explaining the anomalously high $r$ value obtained by BICEP2. We now find $r = 0.004 \pm 0.04$ for BICEP2 and $r = -0.01 \pm 0.04$ for Keck. Since $r \ge 0$ by definition, this implies a maximum likelihood value of $r = 0$, or no evidence for gravitational waves. Starobinsky Inflation ($r = 0.0036$) is not ruled out, however. Krauss & Wilzcek (2014) have already argued that "measurement of polarization of the CMB due to a long-wavelength stochastic background of gravitational waves from Inflation in the early Universe would firmly establish the quantization of gravity," and, therefore, the existence of gravitons. We argue it would also constitute a detection of gravitational Hawking radiation (explicitly from the causal horizons due to Inflation).

We have analyzed the genus topology of the BICEP2 B-modes and find them to be Gaussian random phase as expected if they have a cosmological origin. These BICEP2 B-modes can be produced by gravity waves in the early universe, but some question has arisen as to whether these B-modes (for 50 < l < 120) may instead be produced by foreground polarized dust emission. The dust emission at 150 GHz observed by BICEP2 should be less in magnitude but have similar structure to that at 353 GHz. We have therefore calculated and mapped the B-modes in the BICEP2 region from the publicly available Q and U 353 GHz preliminary Planck polarization maps. These have a genus curve that is different from that seen in the BICEP2 observations, with features at different locations from those in the BICEP2 map. The two maps show a positive correlation coefficient of 15.2% +/- 3.9% (1-sigma). This requires the amplitude of the Planck (50 < l <120) dust modes to be low in the BICEP2 region, and the majority of the Planck 353 GHz signal in the BICEP2 region in these modes to be noise. We can explain the observed correlation coefficient of 15.2% with a BICEP2 gravity wave signal with an rms amplitude equal to 54% of the total BICEP2 rms amplitude. The gravity wave signal corresponds to a tensor-to-scalar ratio r = 0.11 +/- 0.04 (1-sigma). This is consistent with a gravity wave signal having been detected, at a 2.5-sigma level. The Planck and BICEP2 teams have recently engaged in joint analysis of their combined data|it will be interesting to see if that collaboration reaches similar conclusions.

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen in the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data. A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation. We set 95% confidence limits on non-linearities of -101 < f_{nl} < 107. We also find that the observed low l (l <= 8) modes show a slight anti-correlation with the Galactic foreground, but not exceeding 95% confidence, and that the topology defined by these modes is consistent with that of a Gaussian random-phase distribution (within 95% confidence).

It is useful to have mathematical criteria for evaluating errors in map projections. The Chebyshev criterion for minimizing rms (root mean square) local scale factor errors for conformal maps has been useful in developing conformal map projections of continents. Any local error criterion will be minimized ultimately by map projections with multiple interruptions, on which some pairs of points that are close on the globe are far apart on the map. Since it is as bad to have two points on the map at two times their proper separation as to have them at half their proper separation, it is the rms logarithmic distance, s, between random points in the mapped region that we will minimize. The best previously known projection of the entire sphere for distances is the Lambert equal-area azimuthal with an rms logarithmic distance error of s=0.343. For comparison, the Mercator has s=0.444, and the Mollweide has s=0.390. We present new projections: the "Gott equal-area elliptical" with perfect shapes on the central meridian, the "Gott-Mugnolo equal-area elliptical" and the "Gott-Mugnolo azimuthal" with rms logarithmic distance errors of s=0.365, s=0.348, and s=0.341 respectively, which improve on previous projections of their type. The "Gott-Mugnolo azimuthal" has the lowest distance errors of any map and is produced by a new technique using "forces" between pairs of points on a map which make them move so as to minimize s. The "Gott equal-area elliptical" projection produces a particularly attractive map of Mars, and the "Gott-Mugnolo azimuthal" projection produces an interesting map of the moon.

We re-analyze brightness data sampled intensively over 5 nights at two epochs separated by the quasar lens time delay, to examine the nature of the observed microlensing. We find strong evidence for a microlensing event with an amplitude of 1% and a time scale of twelve hours. The existence of such rapid microlensing, albeit at low amplitude, imposes constraints on the nature of the quasar and of the baryonic dark matter.

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen by the Wilkinson Microwave Anisotropy Probe (WMAP). A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation.

We report on an observing campaign in March 2001 to monitor the brightness of the later arriving Q0957+561 B image in order to compare with the previously published brightness observations of the (first arriving) A image. The 12 participating observatories provided 3543 image frames which we have analyzed for brightness fluctuations. From our classical methods for time delay determination, we find a 417.09 +/- 0.07 day time delay which should be free of effects due to incomplete sampling. During the campaign period, the quasar brightness was relatively constant and only small fluctuations were found; we compare the structure function for the new data with structure function estimates for the 1995--6 epoch, and show that the structure function is statistically non-stationary. We also examine the data for any evidence of correlated fluctuations at zero lag. We discuss the limits to our ability to measure the cosmological time delay if the quasar's emitting surface is time resolved, as seems likely.

We present a data set of images of the gravitationally lensed quasar Q2237+0305, that was obtained at the Apache Point Observatory (APO) between June 1995 and January 1998. Although the images were taken under variable, often poor seeing conditions and with coarse pixel sampling, photometry is possible for the two brighter quasar images A and B with the help of exact quasar image positions from HST observations. We obtain a light curve with 73 data points for each of the images A and B. There is evidence for a long (>~ 100 day) brightness peak in image A in 1996 with an amplitude of about 0.4 to 0.5 mag (relative to 1995), which indicates that microlensing has been taking place in the lensing galaxy. Image B does not vary much over the course of the observation period. The long, smooth variation of the light curve is similar to the results from the OGLE monitoring of the system (Wozniak et al. 2000a).

An observing campaign with 10 participating observatories has undertaken to monitor the optical brightness of the Q0957 gravitationally lensed quasar for 10 consecutive nights in January 2000. The resulting A image brightness curve has significant brightness fluctuations and makes a photometric prediction for the B image light curve for a second campaign planned for 12-21 March 2001. The ultimate purpose is to determine the gravitational lens time delay to a fraction of an hour, and to seek evidence for rapid microlensing.

If the halo of the lensing galaxy 0957+561 is made of massive compact objects (MACHOs), they must affect the lightcurves of the quasar images Q0957+561 A and B differently. We search for this microlensing effect in the double quasar by comparing monitoring data for the two images A and B - obtained with the 3.5m Apache Point Observatory from 1995 to 1998 - with intensive numerical simulations. This way we test whether the halo of the lensing galaxy can be made of MACHOs of various masses. We can exclude a halo entirely made out of MACHOs with masses between 10e(-6) Msun and 10e(-2) Msun for quasar sizes of less than 3*10e(14)cm (h = 0.6), hereby extending previous limits upwards by one order of magnitude.

We measure the topology (genus curve) of the galaxy distribution in a mock redshift catalog designed to resemble the upcoming Sloan Digital Sky Survey (SDSS). The catalog, drawn from a large N-body simulation of a Lambda-CDM cos- mological model, mimics the anticipated spectroscopic selection procedures of the SDSS in some detail. Sky maps, redshift slices, and 3-D contour maps of the mock survey reveal a rich and complex structure, including networks of voids and superclusters that resemble the patterns seen in the CfA redshift survey and the Las Campanas Redshift Survey (LCRS). The 3-D genus curve can be measured from the simulated catalog with superb precision; this curve has the general shape predicted for Gaussian, random phase initial conditions, but the error bars are small enough to demonstrate with high significance the subtle departures from this shape caused by non-linear gravitational evolution. These distortions have the form predicted by Matsubara's (1994) perturbative anal- ysis, but they are much smaller in amplitude. We also measure the 3-D genus curve of the radial peculiar velocity field measured by applying distance- indicator relations (with realistic errors) to the mock catalog. This genus curve is consistent with the Gaussian random phase prediction, though it is of relatively low precision because of the large smoothing length required to overcome noise in the measured velocity field. Finally, we measure the 2-D topology in redshift slices, similar to early slices from the SDSS and to slices already observed in the LCRS. The genus curves of these slices are consistent with the observed genus curves of the LCRS, providing further evidence in favor of the inflationary CDM model with Omega_M~0.4. The catalog is publicly available at http://www.astronomy.ohio-state.edu/~dhw/sdss.html.

We have continued our effort to re-reduce archival Q0957+561 brightness monitoring data and present results for 1629 R-band images using the methods for galaxy subtraction and seeing correction reported previously. The new dataset comes from 4 observing runs, several nights apiece, with sampling of typically 5 minutes, which allows the first measurement of the structure function for variations in the R-band from timescales of hours to years. Comparison of our reductions to previous reductions of the same data, and to r-band photometry produced at Apache Point Observatory shows good overall agreement. Two of the data runs, separated by 417 days, permit a sharpened value for the time delay of 417.4 days, valid only if the time delay is close to the now-fashionable 417-day value; our data do not constrain a delay if it is more than three days from this 417-day estimate. Our present results show no unambiguous signature of the daily microlensing, though a suggestive feature is found in the data. Both time delay measurement and microlensing searches suffer from from the lack of sampling at half-day offsets, inevitable at a single observatory, hence the need for round-the-clock monitoring with participation by multiple observatories.

Since the persuasive determination of the time-delay in Q0957+561, much interest has centered around shifting and subtracting the A and B light-curves to look for residuals due to microlensing. Solar mass objects in the lens galaxy produce variations on timescales of decades, with amplitudes of a few tenths of a magnitude, but MACHO's (with masses of order $10^{-3}$ to $10^{-7}M_\odot$) produce variations at only the 5% level. To detect such small variations, highly precise photometry is required. To that end, we have used 200 observations over three nights to examine the effects of seeing on the light-curves. We have determined that seeing itself can be responsible for correlated 5% variations in the light-curves of A and B. We have found, however, that these effects can be accurately removed, by subtracting the light from the lens galaxy, and by correcting for cross contamination of light between the closely juxtaposed A and B images. We find that these corrections improve the variations due to seeing from 5% to a level only marginally detectable over photon shot noise (0.5%).

We have simulated cosmic microwave background (CMB) anisotropy maps for several COBE-DMR-normalized cold dark matter (CDM) cosmogonies, to make predictions for the upcoming MAP experiment. We have studied the sensitivity of the simulated MAP data to cosmology, sky coverage, and instrumental noise. With accurate knowledge of instrumental noise, MAP data will discriminate among the cosmogonies considered, and determine the topology of the initial fluctuations. A correlation function analysis of the simulated MAP data results in a very accurate measurement of the acoustic Hubble radius at decoupling. A low-density open CDM model with Omega_0=0.4 can be distinguished from the Omega_0=1 fiducial CDM model or a Lambda CDM model with > 99% confidence from the location of the acoustic "valley" in the correlation function. A genus analysis of the simulated MAP data indicates that in cosmogonies with Gaussian random-phase initial conditions, a shift of the zero-crossing point of the genus curve near the mean temperature threshold level nu=0 should not exceed |Delta nu| = 0.01 (0.04) when the total effective FWHM smoothing is 0.3 deg (1.0 deg). The asymmetry of the genus curve at the positive and negative threshold levels should not exceed |Delta g/g(nu=1,-1)| = 0.8% (4%) at 0.3 deg (1.0 deg) FWHM smoothing. Deviations of the observed MAP data in excess of these small values will be evidence for non-Gaussian behavior. The amplitude of the genus curve is a measure of the shape of the power spectrum at the smoothing scale. Even with the expected amount of instrumental noise and partial sky coverage (due to the Galaxy), the MAP data should allow discrimination amongst the cosmogonies considered at more than 99% confidence solely from a genus amplitude analysis.

We have previously found in the Hubble Deep Field a significant angular correlation of faint, high color-redshift objects on scales below one arcsecond, or several kiloparsecs in metric size. We examine the correlation and nearest neighbor statistics to conclude that 38% of these objects in the HDF have a companion within one arcsecond, three times the number expected in a random distribution with the same number of objects. We examine three dynamical scenarios for these object multiplets: 1) the objects are star-forming regions within normal galaxies, whose disks have been relatively dimmed by K-correction and surface brightness dimming; 2) they are fragments merging into large galaxies; 3) they are satellites accreting onto normal L_* galaxies. We find that hypothesis 1 is most tenable. First, large galaxies in the process of a merger formation would have accumulated too much mass in their centers (5e12 M_sun inside 2 kpc) to correspond to present day objects. Second, accretion by dynamical friction occurs with a predictable density vs. radius slope, not seen among the faint HDF objects. Since the dynamical friction time is roughly (1 Gyr), a steady-state should have been reached by redshift z < 5. Star-forming regions within galaxies clearly present no dynamical problems. Since large spirals would still appear as such in the HDF, we favor a scenario in which the faint compact sources in the HDF are giant starforming regions within small normal galaxies, such as Magellanic irregulars. Finally we checked that reduction in mass-to-light from induced star-formation cannot alone explain the luminosity overdensity.

Photometric monitoring of the gravitational lens system 0957+561A,B in the g and r bands with the Apache Point Observatory (APO) 3.5 m telescope during 1996 shows a sharp g band event in the trailing (B) image light curve at the precise time predicted from the observation of an event during 1995 in the leading (A) image with a delay of 415 days. This success confirms the "short delay," and the lack of any feature at a delay near 540 days rejects the "long delay" for this system, resolving a long-standing controversy. A series of statistical analyses of our light curve data yield a best fit delay of 417 +/- 3 days (95% confidence interval). Recent improvements in the modeling of the lens system (consisting of a galaxy and cluster) allow us to derive a value of the global (at z = 0.36) value of Hubble's constant H_0 using Refsdal's method, a simple and direct distance determination based on securely understood physics and geometry. The result is H_0 = 63 +/- 12 km/s/Mpc (for Omega = 1) where this 95% confidence interval is dominated by remaining lens model uncertainties.

We have measured the topology (genus) of the density distribution of large-scale structure observed in the Las Campanas Redshift Survey (LCRS). The LCRS is complete to magnitude 17.5, and contains nearly 24000 galaxies with median redshift of 30000 km/s. The large volume and large number of galaxies allows sampling of nearly 100 independent structures with which to compute the genus topology, a vast improvement over previous studies. We find that the genus is consistent with a random-phase Gaussian distribution of initial density fluctuations, as would be produced naturally in inflationary models. When we combine these results with the genus measurements of the COBE microwave background fluctuations, we find that two orthogonal projections of the three-dimensional distribution of initial density fluctuations are consistent with Gaussian random-phase behavior, in agreement with standard inflationary models. Particular attention is given to statistical significance of the genus test.

We compute the two-point angular correlation function and number-magnitude relation of Hubble Deep Field sources in order to assess their nature. We find that the correlation peaks between 0.25 arcsec and 0.4 arcsec with amplitude of 2 or greater, and much more for the smallest objects. This angular scale corresponds to physical scales of order 1 kpc for redshifts z > 1. The correlation must therefore derive from objects with subgalaxian separations. At faint magnitudes, the counts satisfy the relation Number prop. to 1/flux, expected for images which are subdivisions of larger ones. A conservative conjecture may explain these results. Since high redshift space (z > 0.5) dominates the volume of the sample, observational redshift effects are important. The K-correction and surface brightness dimming of diffuse sources enhances the prominence of compact, unresolved, UV-bright objects, such as star-forming regions within normal gas-rich galaxies. If such regions appear as individual sources, they could explain the subgalaxian correlation.

We have measured the topology (genus) of the fluctuations in the cosmic microwave background seen in the recently completed (four-year) data set produced by the COBE satellite. We find that the genus is consistent with that expected from a random-phase Gaussian distribution, as might be produced naturally in inflationary models.

We have used a set of equations developed by Pringle (1992) to follow the evolution of a viscous twisted disk in a galaxy-like potential which is stationary or tumbling relative to inertial space. In an axisymmetric potential, the disk settles to the equatorial plane at a rate determined largely by the coefficient $\nu_2$ associated with shear perpendicular to the local disk plane. If the disk is initially close to the galaxy equator, the rate at which the inclination decays is well described by the analytic formula of Steiman-Cameron and Durisen (1988); in a highly inclined disk, `breaking waves' of curvature steepen as they propagate through the disk, rendering the numerical treatment untrustworthy. In a triaxial potential which is stationary in inertial space, settling is faster than in an oblate or prolate galaxy, since the disk twists simultaneously about two perpendicular axes. If the figure of the potential tumbles about one of its principal axes, the viscous disk can settle into a warped state in which gas at each radius follows a stable tilted orbit which precesses so as to remain stationary relative to the underlying galaxy.

A unique reconstruction of the image of a high redshift source galaxy responsible for multiple long arcs in the z = 0.4 cluster 0024+1654 is obtained by inverse lensing. Deep B and I imaging with the Hubble Space Telescope (Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by AURA, under NASA contract NAS 5-26555) enables high resolution of the arcs due to strong gravitational lensing of the background source. The gravitational lens in the foreground cluster is thus used to obtain a magnified view of the distant source. Four strongly lensed images of the source lead to a unique reconstruction. Each of the long arcs, when unlensed, leads to the same reconstructed source image exhibiting a beaded ring-like morphology. The U luminosity of the ring alone is equivalent to a normal galaxy. This is likely a galaxy in formation.

CCD photometry of the gravitational lens system 0957+561A,B in the g and r bands was obtained on alternate nights, weather permitting, from December 1994 through May 1995 using the Double Imaging Spectrograph (DIS) on the Apache Point Observatory (APO) 3.5-meter telescope. The remote observing and fast instrument change capabilities of this facility allowed accumulation of light curves sampled frequently and consistently. The Honeycutt ensemble photometry algorithm was applied to the data set and yielded typical relative photometric errors of approximately 0.01 magnitudes. Image A exhibited a sharp drop of about 0.1 magnitudes in late December 1994; no other strong features were recorded in either image. This event displays none of the expected generic features of a microlensing-induced flux variation and is likely to be intrinsic to the quasar; if so, it should also be seen in the B image with the lensing differential time delay. We give the expected 1996 image B light curves based on two values of the time delay and brightness ratio which have been proposed and debated in the literature. Continued monitoring of the system in the first half of 1996 should easily detect the image B event and thus resolve the time-delay controversy.

We consider the feasibility of directly observing gravitational microlensing in extra-galactic sources, whose stars are not generally resolved. This precludes use of the simple optical depth to microlensing formulation, which is applicable only to resolved stars. We, instead, extend this method to consider observational constraints, such as seeing, sky background and minimum detectable change in magnitude. Our analysis provides quantitative relations between these constraints and the expected observational results, event duration and number of detections per year. We find that an ambitious ground-based observer should detect several events per year in M31. We also consider detection of microlensing in visual binary galaxies. We find that, although these objects may present hundreds of events per year, extremely short durations would yield poor prospects for observation.