This paper presents data and analysis of SN 2010kd, a low-redshift ($z = 0.101$) H-deficient superluminous supernova (SLSN), based on ultraviolet/optical photometry and optical spectroscopy spanning between $-$28 and +194 days relative to $\mathit{B}$ band maximum light. The $\mathit{B}$ band light curve comparison of SN 2010kd with a subset of well-studied SLSNe I at comparable redshifts indicates that it is a slow-decaying PTF12dam like SLSN. Analytical light-curve modeling using the $\mathtt{Minim}$ code suggests that the bolometric light curve of SN 2010kd favors circumstellar matter interaction for the powering mechanism. $\mathtt{SYNAPPS}$ modeling of the early-phase spectra does not identify broad H or He lines, whereas the photospheric-phase spectra are dominated by O I, O II, C II, C IV and Si II, particularly, presence of both low and high-velocity components of O II and Si II lines. The nebular-phase spectra of SN 2010kd are dominated by O I and Ca II emission lines similar to those seen in other SLSNe I. The line velocities in SN 2010kd exhibit flatter evolution curves similar to SN 2015bn but with comparatively higher values. SN 2010kd shows a higher single-zone local thermodynamic equilibrium temperature in comparison to PTF12dam and SN 2015bn, and it has an upper O I ejected mass limit of $\sim 10~M_\odot$. The host of SN 2010kd is a dwarf galaxy with a high star-formation rate ($\sim 0.18 \pm 0.04~M_\odot$ yr$^{-1}$) and extreme emission lines.

Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le \rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter function, including its characteristic break ($m^*$) to a faint end power-law slope ($\alpha$). Our method enables us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant ($\sim 1.9 \sigma$) evolution in the faint end slope $\alpha$ versus redshift. We also find no dependence in $\alpha$ or $m^*$ with the X-ray inferred cluster masses. However, the amplitude of the LF as a function of cluster mass is constrained to $\sim 20\%$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the model prediction. We estimate the average growth rate since $z = 1.0$ for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the $0.40$ dex growth rate expected from the semi-analytic model. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.

Deep optical images are often crowded with overlapping objects. This is especially true in the cores of galaxy clusters, where images of dozens of galaxies may lie atop one another. Accurate measurements of cluster properties require deblending algorithms designed to automatically extract a list of individual objects and decide what fraction of the light in each pixel comes from each object. We present new software called the Gradient And INterpolation based deblender (GAIN) as a secondary deblender to improve deblending the images of cluster cores. This software relies on using image intensity gradient and using an image interpolation technique usually used to correct flawed terrestrial digital images. We test this software on Dark Energy Survey coadd images. GAIN helps extracting unbiased photometry measurement for blended sources. It also helps improving detection completeness while introducing only a modest amount of spurious detections. For example, when applied to deep images simulated with high level of deblending difficulties, this software improves detection completeness from 91% to 97% for sources above the 10? limiting magnitude at 25.3 mag. We expect this software to be a useful tool for cluster population measurements.

We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samples are derived. The optical LFs at 100s are fitted by assuming that the co-moving GRB rate traces the star-formation rate. The detection function of ROTSE-III is taken into account during the fitting of the optical LFs by using Monte Carlo simulations. We find that the cumulative distribution of optical emission at 100s is well-described with an exponential rise and power-law decay (ERPLD), broken power-law (BPL), and Schechter LFs. A single power-law (SPL) LF, on the other hand, is ruled out with high confidence.

Weak-lensing measurements of the averaged shear profiles of galaxy clusters binned by some proxy for cluster mass are commonly converted to cluster mass estimates under the assumption that these cluster stacks have spherical symmetry. In this paper we test whether this assumption holds for optically selected clusters binned by estimated optical richness. Using mock catalogues created from N-body simulations populated realistically with galaxies, we ran a suite of optical cluster finders and estimated their optical richness. We binned galaxy clusters by true cluster mass and estimated optical richness and measure the ellipticity of these stacks. We find that the processes of optical cluster selection and richness estimation are biased, leading to stacked structures that are elongated along the line-of-sight. We show that weak-lensing alone cannot measure the size of this orientation bias. Weak lensing masses of stacked optically selected clusters are overestimated by up to 3-6 per cent when clusters can be uniquely associated with haloes. This effect is large enough to lead to significant biases in the cosmological parameters derived from large surveys like the Dark Energy Survey, if not calibrated via simulations or fitted simultaneously. This bias probably also contributes to the observed discrepancy between the observed and predicted Sunyaev-Zel'dovich signal of optically-selected clusters.

We present a method to study the photometric properties of galaxies in filaments by stacking the galaxy populations between pairs of galaxy clusters. Using Sloan Digital Sky Survey data, this method can detect the inter-cluster filament galaxy overdensity with a significance of $\sim 5 \sigma$ out to $z=0.40$. Using this approach, we study the $g-r$ color and luminosity distribution of filament galaxies as a function of redshift. Consistent with expectation, filament galaxies are bimodal in their color distribution and contain a larger blue galaxy population than clusters. Filament galaxies are also generally fainter than cluster galaxies. More interestingly, the observed filament population seems to show redshift evolution at $0.12<z<0.40$: the blue galaxy fraction has a trend to increase at higher redshift: a filament "Butcher Oemler Effect". We test the dependence of the observed filament density on the richness of the cluster pair: richer clusters are connected by higher density filaments. We also test the spatial dependence of filament galaxy overdensity: this quantity decreases when moving away from the inter-cluster axis between a cluster pair. This method provides an economical way to probe the photometric properties of filament galaxies and should prove useful for upcoming projects like the Dark Energy Survey.

We present the first results of the SOAR (Southern Astrophysical Research) Gravitational Arc Survey (SOGRAS). The survey imaged 47 clusters in two redshift intervals centered at $z=0.27$ and $z=0.55$, targeting the richest clusters in each interval. Images were obtained in the $g'$, $r'$ and $i'$ bands using the SOAR Optical Imager (SOI), with a median seeing of 0.83, 0.76 and 0.71 arcsec, respectively, in these filters. Most of the survey clusters are located within the Sloan Digital Sky Survey (SDSS) Stripe 82 region and all of them are in the SDSS footprint. Photometric calibration was therefore performed using SDSS stars located in our SOI fields. We reached for galaxies in all fields the detection limits of $g \sim 23.5$, $r \sim 23$ and $i \sim 22.5$ for a signal-to-noise ratio (S/N) = 3. As a by-product of the image processing, we generated a source catalogue with 19760 entries, the vast majority of which are galaxies, where we list their positions, magnitudes and shape parameters. We compared our galaxy shape measurements to those of local galaxies and concluded that they were not strongly affected by seeing. From the catalogue data, we are able to identify a red sequence of galaxies in most clusters in the lower $z$ range. We found 16 gravitational arc candidates around 8 clusters in our sample. They tend to be bluer than the central galaxies in the lensing cluster. A preliminary analysis indicates that $\sim 10%$ of the clusters have arcs around them, with a possible indication of a larger efficiency associated to the high-$z$ systems when compared to the low-$z$ ones. Deeper follow-up images with Gemini strengthen the case for the strong lensing nature of the candidates found in this survey.

Launched on June 11, 2008, the LAT instrument onboard the $Fermi$ Gamma-ray Space Telescope has provided a rare opportunity to study high energy photon emission from gamma-ray bursts. Although the majority of such events (27) have been iden tified by the Fermi LAT Collaboration, four were uncovered by using more sensiti ve statistical techniques (Akerlof et al 2010, Akerlof et al 2011, Zheng et al 2 012). In this paper, we continue our earlier work by finding three more GRBs ass ociated with high energy photon emission, GRB 110709A, 111117A and 120107A. To s ystematize our matched filter approach, a pipeline has been developed to identif y these objects in near real time. GRB 120107A is the first product of this anal ysis procedure. Despite the reduced threshold for identification, the number of GRB events has not increased significantly. This relative dearth of events with low photon number prompted a study of the apparent photon number distribution. W e find an extremely good fit to a simple power-law with an exponent of -1.8 $\pm $ 0.3 for the differential distribution. As might be expected, there is a substa ntial correlation between the number of lower energy photons detected by the GBM and the number observed by the LAT. Thus, high energy photon emission is associ ated with some but not all of the brighter GBM events. Deeper studies of the pro perties of the small population of high energy emitting bursts may eventually yi eld a better understanding of these entire phenomena.

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

Since the launch of the Fermi gamma - ray Space Telescope on June 11, 2008, the LAT instrument has solidly detected more than 20 GRBs with high energy photon emission above 100 MeV. Using the matched filter technique, 3 more GRBs have also shown evidence of correlation with high energy photon emission as demonstrated by Akerlof et al. In this paper, we present another GRB unambiguously detected by the matched filter technique, GRB 081006A. This event is associated with more than 13 high energy photons above 100 MeV. The likelihood analysis code provided by the $Fermi$ Science Support Center (FSSC) generated an independent verification of this detection by comparison of the Test Statistics (TS) value with similar calculations for random LAT data fields. We have performed detailed temporal and spectral analysis of photons from 8 keV up to 0.8 GeV from the GBM and the LAT. The properties of GRB 081006A can be compared to the other two long duration GRBs detected at similar significance, GRB 080825C and GRB 090217A. We find that GRB 081006A is more similar to GRB 080825C with comparable appearances of late high energy photon emission. As demonstrated previously, there appears to be a surprising dearth of faint LAT GRBs, with only one additional GRB identified in a sample of 74 GRBs. In this unique period when both $Swift$ and $Fermi$ are operational, there is some urgency to explore this aspect of GRBs as fully as possible.

We report the discovery of 12 new fossil groups of galaxies, systems dominated by a single giant elliptical galaxy and cluster-scale gravitational potential, but lacking the population of bright galaxies typically seen in galaxy clusters. These fossil groups (FGs), selected from the maxBCG optical cluster catalog, were detected in snapshot observations with the Chandra X-ray Observatory. We detail the highly successful selection method, with an 80% success rate in identifying 12 FGs from our target sample of 15 candidates. For 11 of the systems, we determine the X-ray luminosity, temperature, and hydrostatic mass, which do not deviate significantly from expectations for normal systems, spanning a range typical of rich groups and poor clusters of galaxies. A small number of detected FGs are morphologically irregular, possibly due to past mergers, interaction of the intra-group medium (IGM) with a central AGN, or superposition of multiple massive halos. Two-thirds of the X-ray-detected FGs exhibit X-ray emission associated with the central BCG, although we are unable to distinguish between AGN and extended thermal galaxy emission using the current data. This sample, a large increase in the number of known FGs, will be invaluable for future planned observations to determine FG temperature, gas density, metal abundance, and mass distributions, and to compare to normal (non-fossil) systems. Finally, the presence of a population of galaxy-poor systems may bias mass function determinations that measure richness from galaxy counts. When used to constrain power spectrum normalization and {\Omega}_m, these biased mass functions may in turn bias these results.

We present measurements of two types of cluster galaxy alignments based on a volume limited and highly pure ($\ge$ 90%) sample of clusters from the GMBCG catalog derived from SDSS DR7. We detect a clear BCG alignment (the alignment of major axis of the BCG toward the distribution of cluster satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG absolute magnitude decrease, and becomes weaker as BCG stellar mass decreases. No dependence of the BCG alignment on cluster richness is found. We can detect a statistically significant ($\ge$ 3 sigma) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the BCG) only when we use the isophotal fit position angles (PAs, hereafter), and the satellite alignment depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests the detected satellite alignment based on isophotoal PAs from the SDSS pipeline is possibly due to the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with our current photometric SDSS data. We perform our measurements on both SDSS $r$ band and $i$ band data, but did not observe a passband dependence of the alignments.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

Reducing the scatter between cluster mass and optical richness is a key goal for cluster cosmology from photometric catalogs. We consider various modifications to the red-sequence matched filter richness estimator of Rozo et al. (2009), and evaluate their impact on the scatter in X-ray luminosity at fixed richness. Most significantly, we find that deeper luminosity cuts can reduce the recovered scatter, finding that sigma_lnLX|lambda=0.63+/-0.02 for clusters with M_500c >~ 1.6e14 h_70^-1 M_sun. The corresponding scatter in mass at fixed richness is sigma_lnM|lambda ~ 0.2-0.3 depending on the richness, comparable to that for total X-ray luminosity. We find that including blue galaxies in the richness estimate increases the scatter, as does weighting galaxies by their optical luminosity. We further demonstrate that our richness estimator is very robust. Specifically, the filter employed when estimating richness can be calibrated directly from the data, without requiring a-priori calibrations of the red-sequence. We also demonstrate that the recovered richness is robust to up to 50% uncertainties in the galaxy background, as well as to the choice of photometric filter employed, so long as the filters span the 4000 A break of red-sequence galaxies. Consequently, our richness estimator can be used to compare richness estimates of different clusters, even if they do not share the same photometric data. Appendix 1 includes "easy-bake" instructions for implementing our optimal richness estimator, and we are releasing an implementation of the code that works with SDSS data, as well as an augmented maxBCG catalog with the lambda richness measured for each cluster.

We present the results of a comprehensive study of the gamma-ray burst 080928 and of its afterglow. GRB 080928 was a long burst detected by Swift/BAT and Fermi/GBM. It is one of the exceptional cases where optical emission had already been detected when the GRB itself was still radiating in the gamma-ray band. For nearly 100 seconds simultaneous optical, X-ray and gamma-ray data provide a coverage of the spectral energy distribution of the transient source from about 1 eV to 150 keV. In particular, we show that the SED during the main prompt emission phase agrees with synchrotron radiation. We constructed the optical/near-infrared light curve and the spectral energy distribution based on Swift/UVOT, ROTSE-IIIa (Australia), and GROND (La Silla) data and compared it to the X-ray light curve retrieved from the Swift/XRT repository. We show that its bumpy shape can be modeled by multiple energy-injections into the forward shock.Furthermore, we investigate whether the temporal and spectral evolution of the tail emission of the first strong flare seen in the early X-ray light curve can be explained by large-angle emission (LAE). We find that a nonstandard LAE model is required to explain the observations. Finally, we report on the results of our search for the GRB host galaxy, for which only a deep upper limit can be provided.

We present a large catalog of optically selected galaxy clusters from the application of a new Gaussian Mixture Brightest Cluster Galaxy (GMBCG) algorithm to SDSS Data Release 7 data. The algorithm detects clusters by identifying the red sequence plus Brightest Cluster Galaxy (BCG) feature, which is unique for galaxy clusters and does not exist among field galaxies. Red sequence clustering in color space is detected using an Error Corrected Gaussian Mixture Model. We run GMBCG on 8240 square degrees of photometric data from SDSS DR7 to assemble the largest ever optical galaxy cluster catalog, consisting of over 55,000 rich clusters across the redshift range from 0.1 < z < 0.55. We present Monte Carlo tests of completeness and purity and perform cross-matching with X-ray clusters and with the maxBCG sample at low redshift. These tests indicate high completeness and purity across the full redshift range for clusters with 15 or more members.

The red sequence is an important feature of galaxy clusters and plays a crucial role in optical cluster detection. Measurement of the slope and scatter of the red sequence are affected both by selection of red sequence galaxies and measurement errors. In this paper, we describe a new error corrected Gaussian Mixture Model for red sequence galaxy identification. Using this technique, we can remove the effects of measurement error and extract unbiased information about the intrinsic properties of the red sequence. We use this method to select red sequence galaxies in each of the 13,823 clusters in the maxBCG catalog, and measure the red sequence ridgeline location and scatter of each. These measurements provide precise constraints on the variation of the average red galaxy populations in the observed frame with redshift. We find that the scatter of the red sequence ridgeline increases mildly with redshift, and that the slope decreases with redshift. We also observe that the slope does not strongly depend on cluster richness. Using similar methods, we show that this behavior is mirrored in a spectroscopic sample of field galaxies, further emphasizing that ridgeline properties are independent of environment.

From the launch of the Fermi Gamma-ray Space Telescope to July 9, 2010, the Gamma-ray Burst Monitor (GBM) has detected 497 probable GRB events. Twenty-two of these satisfy the simultaneous requirements of an estimated burst direction within 52^\circ of the Fermi Large Area Telescope (LAT) boresight and a low energy fluence exceeding 5 $\mu$erg/cm^2. Using matched filter techniques, the spatially correlated Fermi/LAT photon data above 100 MeV have been examined for evidence of bursts that have so far evaded detection at these energies. High energy emission is detected with great confidence for one event, GRB 090228A. Since the LAT has significantly better angular resolution than the GBM, real-time application of these methods could open the door to optical identification and richer characterization of a larger fraction of the relatively rare GRBs that include high energy emission.

Since the launch of the Fermi Gamma-ray Space Telescope on June 11, 2008, 55 gamma-ray bursts (GRBs) have been observed at coordinates that fall within 66^\circ of the Fermi Large Area Telescope (LAT) boresight with precise localizations provided by the NASA Swift mission or other satellites. Imposing selection cuts to exclude low Galactic latitudes and high zenith angles reduces the sample size to 41. Using matched filter techniques, the Fermi/LAT photon data for these fields have been examined for evidence of bursts that have so far evaded detection at energies above 100 MeV. Following comparisons with similar random background fields, two events, GRB 080905A and GRB 091208B, stand out as excellent candidates for such an identification. After excluding the six bright bursts previously reported by the LAT team, the remaining 35 events exhibit an excess of LAT "diffuse" photons with a statistical significance greater than 2 sigma, independent of the matched filter analysis. After accounting for the total number of photons in the well-localized fields and including estimates of detection efficiency, one concludes that somewhere in the range of 11% to 19% of all GRBs within the LAT field of view illuminate the detector with two or more energetic photons. These are the most stringent estimates of the high energy photon content of GRBs to date. The two new events associated with high energy photon emission have similar ratios of high to low energy fluences as observed previously. This separates them from bursts with similar low energy fluences by a factor of ten, suggesting a distinct class of events rather than a smooth continuum.

SN 2007if was the third over-luminous SN Ia detected after 2003fg and 2006gz. We present the photometric and spectroscopic observations of the supernova and its host by ROTSE-III, HET and Keck. From the H_alpha line identified in the host spectra, we determine a redshift of 0.0736. At this distance, the supernova reached an absolute magnitude of -20.4, brighter than any other SNe Ia ever observed. If the source of luminosity is radioactive decay, a large amount of radioactive nickel (~1.5 solar masses) is required to power the peak luminosity, more than can be produced realistically in a Chandrasekhar mass progenitor. Low expansion velocity, similar to that of 2003fg, is also measured around the maximum light. The observations may suggest that SN 2007if was from a massive white dwarf progenitor, plausibly exploding with mass well beyond 1.4 solar masses. Alternatively, we investigate circumstellar interaction that may contribute to the excess luminosity.

The optical-infrared afterglow of the LAT-detected long duration burst, GRB 090902B, has been observed by several instruments. The earliest detection by ROTSE-IIIa occurred 80 minutes after detection by the GBM instrument onboard the Fermi Gamma-Ray Space Telescope, revealing a bright afterglow and a decay slope suggestive of a reverse shock origin. Subsequent optical-IR observations followed the light curve for 6.5 days. The temporal and spectral behavior at optical-infrared frequencies is consistent with synchrotron fireball model predictions; the cooling break lies between optical and XRT frequencies ~ 1.9 days after the burst. The inferred electron energy index is $p = 1.8 \pm 0.2$, which would however imply an X-ray decay slope flatter than observed. The XRT and LAT data have similar spectral indices and the observed steeper value of the LAT temporal index is marginally consistent with the predicted temporal decay in the radiative regime of the forward shock model. Absence of a jet break during the first 6 days implies a collimation-corrected $\gamma$-ray energy $E_{\gamma} > 2.2\times10^{52}\rm$ ergs, one of the highest ever seen in a long-duration GRBs. More events combining GeV photon emission with multi-wavelength observations will be required to constrain the nature of the central engine powering these energetic explosions and to explore the correlations between energetic quanta and afterglow emission.

We examine the alignment between Brightest Cluster Galaxies (BCGs) and their host clusters in a sample of 7031 clusters with 0.08<z<0.44 found using a matched-filter algorithm and an independent sample of 5744 clusters with 0.1<z<0.3 selected with the maxBCG algorithm, both extracted from the Sloan Digital Sky Survey Data Release 6 imaging data. We confirm that BCGs are preferentially aligned with the cluster's major axis; clusters with dominant BCGs (>0.65 mag brighter than the mean of the second and third ranked galaxies) show stronger alignment than do clusters with less dominant BCGs at the 4.4 sigma level. Rich clusters show a stronger alignment than do poor clusters at the 2.3 sigma level. Low redshift clusters (z<0.26) show more alignment than do high redshift (z>0.26) clusters, with a difference significant at the 3.0 sigma level. Our results do not depend on the algorithm used to select the cluster sample, suggesting that they are not biased by systematics of either algorithm. The correlation between BCG dominance and cluster alignment may be a consequence of the hierarchical merging process which forms the cluster. The observed redshift evolution may follow from secondary infall at late redshifts.

We present a multi-wavelength study of GRB 081008, at redshift 1.967, by Swift, ROTSE-III and GROND. Compared to other Swift GRBs, GRB 081008 has a typical gamma-ray isotropic equivalent energy output (10^53 erg) during the prompt phase, and displayed two temporally separated clusters of pulses. The early X-ray emission seen by the Swift/XRT was dominated by the softening tail of the prompt emission, producing multiple flares during and after the Swift/BAT detections. Optical observations that started shortly after the first active phase of gamma-ray emission showed two consecutive peaks. We interpret the first optical peak as the onset of the afterglow associated with the early burst activities. A second optical peak, coincident with the later gamma-ray pulses, imposes a small modification to the otherwise smooth lightcurve and thus suggests a minimal contribution from a probable internal component. We suggest the early optical variability may be from continuous energy injection into the forward shock front by later shells producing the second epoch of burst activities. These early observations thus provide a potential probe for the transition from prompt to the afterglow phase. The later lightcurve of GRB 081008 displays a smooth steepening in all optical bands and X-ray. The temporal break is consistent with being achromatic at the observed wavelengths. Our broad energy coverage shortly after the break constrains a spectral break within optical. However, the evolution of the break frequency is not observed. We discuss the plausible interpretations of this behavior.

The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sample represents the final set of galaxies observed using the original SDSS target selection criteria. We analyse the clustering of galaxies within this sample, including both the Luminous Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift Survey (2dFGRS) data. Baryon Acoustic Oscillations are observed in power spectra measured for different slices in redshift; this allows us to constrain the distance--redshift relation at multiple epochs. We achieve a distance measure at redshift z=0.275, of r_s(z_d)/D_V(0.275)=0.1390+/-0.0037 (2.7% accuracy), where r_s(z_d) is the comoving sound horizon at the baryon drag epoch, D_V(z)=[(1+z)^2D_A^2cz/H(z)]^(1/3), D_A(z) is the angular diameter distance and H(z) is the Hubble parameter. We find an almost independent constraint on the ratio of distances D_V(0.35)/D_V(0.2)=1.736+/-0.065, which is consistent at the 1.1sigma level with the best fit Lambda-CDM model obtained when combining our z=0.275 distance constraint with the WMAP 5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, but the discrepancy is now of lower significance, a change caused by a revised error analysis and a change in the methodology adopted, as well as the addition of more data. Using WMAP5 constraints on Omega_bh^2 and Omega_ch^2, and combining our BAO distance measurements with those from the Union Supernova sample, places a tight constraint on Omega_m=0.286+/-0.018 and H_0 = 68.2+/-2.2km/s/Mpc that is robust to allowing curvature and non-Lambda dark energy. This result is independent of the behaviour of dark energy at redshifts greater than those probed by the BAO and supernova measurements. (abridged)

We present Chandra snapshot observations of the first large X-ray sample of optically identified fossil groups. For 9 of 14 candidate groups, we are able to determine the X-ray luminosity and temperature, which span a range typical of large ellipticals to rich groups of galaxies. We discuss these initial results in the context of group IGM and central galaxy ISM evolution, and we also describe plans for a deep X-ray follow-up program.

Precision photometric redshifts will be essential for extracting cosmological parameters from the next generation of wide-area imaging surveys. In this paper we introduce a photometric redshift algorithm, ArborZ, based on the machine-learning technique of Boosted Decision Trees. We study the algorithm using galaxies from the Sloan Digital Sky Survey and from mock catalogs intended to simulate both the SDSS and the upcoming Dark Energy Survey. We show that it improves upon the performance of existing algorithms. Moreover, the method naturally leads to the reconstruction of a full probability density function (PDF) for the photometric redshift of each galaxy, not merely a single "best estimate" and error, and also provides a photo-z quality figure-of-merit for each galaxy that can be used to reject outliers. We show that the stacked PDFs yield a more accurate reconstruction of the redshift distribution N(z). We discuss limitations of the current algorithm and ideas for future work.

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find \Omega_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and \Omega_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find \Omega_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_{\nu} < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find \Omega_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

We report on a complete set of early optical afterglows of gamma-ray bursts (GRBs) obtained with the ROTSE-III telescope network from March 2005 through June 2007. This set is comprised of 12 afterglows with early optical and Swift/XRT observations, with a median ROTSE-III response time of 45 s after the start of gamma-ray emission (8 s after the GCN notice time). These afterglows span four orders of magnitude in optical luminosity, and the contemporaneous X-ray detections allow multi-wavelength spectral analysis. Excluding X-ray flares, the broadband synchrotron spectra show that the optical and X-ray emission originate in a common region, consistent with predictions of the external forward shock in the fireball model. However, the fireball model is inadequate to predict the temporal decay indices of the early afterglows, even after accounting for possible long-duration continuous energy injection. We find that the optical afterglow is a clean tracer of the forward shock, and we use the peak time of the forward shock to estimate the initial bulk Lorentz factor of the GRB outflow, and find 100<Gamma_0<1000, consistent with expectations.

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find \sigma_8(\Omega_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find \sigma_8=0.807\pm 0.020 and \Omega_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

We report the discovery by the Robotic Optical Transient Experiment (ROTSE-IIIb) telescope of SN 2008es, an overluminous supernova (SN) at z=0.205 with a peak visual magnitude of -22.2. We present multiwavelength follow-up observations with the Swift satellite and several ground-based optical telescopes. The ROTSE-IIIb observations constrain the time of explosion to be 23+/-1 rest-frame days before maximum. The linear decay of the optical light curve, and the combination of a symmetric, broad H\alpha emission line profile with broad P Cygni H\beta and Na I \lambda5892 profiles, are properties reminiscent of the bright Type II-L SNe 1979C and 1980K, although SN 2008es is greater than 10 times more luminous. The host galaxy is undetected in pre-supernova Sloan Digital Sky Survey images, and similar to Type II-L SN 2005ap (the most luminous SN ever observed), the host is most likely a dwarf galaxy with M_r > -17. Swift Ultraviolet/Optical Telescope observations in combination with Palomar photometry measure the SED of the SN from 200 to 800 nm to be a blackbody that cools from a temperature of 14,000 K at the time of the optical peak to 6400 K 65 days later. The inferred blackbody radius is in good agreement with the radius expected for the expansion speed measured from the broad lines (10,000 km/s). The bolometric luminosity at the optical peak is 2.8 x 10^44 erg/s, with a total energy radiated over the next 65 days of 5.6 x 10^50 erg. We favor a model in which the exceptional peak luminosity is a consequence of the core-collapse explosion of a progenitor star with a low-mass extended hydrogen envelope and a stellar wind with a density close to the upper limit on the mass-loss rate measured from the lack of an X-ray detection by the Swift X-Ray Telescope. (Abridged).

The optical afterglow of long-duration GRB 071003 is among the brightest yet to be detected from any GRB, with R ~ 12 mag in KAIT observations starting 42 s after the GRB trigger, including filtered detections during prompt emission. However, our high S/N ratio afterglow spectrum displays only extremely weak absorption lines at what we argue is the host redshift of z = 1.60435 - in contrast to the three other, much stronger Mg II absorption systems observed at lower redshifts. Together with Keck adaptive optics observations which fail to reveal a host galaxy coincident with the burst position, our observations suggest a halo progenitor and offer a cautionary tale about the use of Mg II for GRB redshift determination. We present early through late-time observations spanning the electromagnetic spectrum, constrain the connection between the prompt emission and early variations in the light curve (we observe no correlation), and discuss possible origins for an unusual, marked rebrightening that occurs a few hours after the burst: likely either a late-time refreshed shock or a wide-angle secondary jet. Analysis of the late-time afterglow is most consistent with a wind environment, suggesting a massive star progenitor. Together with GRB 070125, this may indicate that a small but significant portion of star formation in the early universe occurred far outside what we consider a normal galactic disk.

Minimizing the scatter between cluster mass and accessible observables is an important goal for cluster cosmology. In this work, we introduce a new matched filter richness estimator, and test its performance using the maxBCG cluster catalog. Our new estimator significantly reduces the variance in the L_X-richness relation, from \sigma_{\ln L_X}^2=(0.86\pm0.02)^2 to \sigma_{\ln L_X}^2=(0.69\pm0.02)^2. Relative to the maxBCG richness estimate, it also removes the strong redshift dependence of the richness scaling relations, and is significantly more robust to photometric and redshift errors. These improvements are largely due to our more sophisticated treatment of galaxy color data. We also demonstrate the scatter in the L_X-richness relation depends on the aperture used to estimate cluster richness, and introduce a novel approach for optimizing said aperture which can be easily generalized to other mass tracers.

We measure the logarithmic scatter in mass at fixed richness for clusters in the maxBCG cluster catalog, an optically selected cluster sample drawn from SDSS imaging data. Our measurement is achieved by demanding consistency between available weak lensing and X-ray measurements of the maxBCG clusters, and the X-ray luminosity--mass relation inferred from the 400d X-ray cluster survey, a flux limited X-ray cluster survey. We find \sigma_{\ln M|N_{200}}=0.45^{+0.20}_{-0.18} (95% CL) at N_{200} ~ 40, where N_{200} is the number of red sequence galaxies in a cluster. As a byproduct of our analysis, we also obtain a constraint on the correlation coefficient between \ln Lx and \ln M at fixed richness, which is best expressed as a lower limit, r_{L,M|N} >= 0.85 (95% CL). This is the first observational constraint placed on a correlation coefficient involving two different cluster mass tracers. We use our results to produce a state of the art estimate of the halo mass function at z=0.23 -- the median redshift of the maxBCG cluster sample -- and find that it is consistent with the WMAP5 cosmology. Both the mass function data and its covariance matrix are presented.

We present a new measurement of the scaling relation between X-ray luminosity and total mass for 17,000 galaxy clusters in the maxBCG cluster sample. Stacking sub-samples within fixed ranges of optical richness, N_200, we measure the mean 0.1-2.4 keV X-ray luminosity, <L_X>, from the ROSAT All-Sky Survey. The mean mass, <M_200>, is measured from weak gravitational lensing of SDSS background galaxies (Johnston et al. 2007). For 9 <= N_200 < 200, the data are well fit by a power-law, <L_X>/10^42 h^-2 erg/s = (12.6+1.4-1.3 (stat) +/- 1.6 (sys)) (<M_200>/10^14 h^-1 M_sun)^1.65+/-0.13. The slope agrees to within 10% with previous estimates based on X-ray selected catalogs, implying that the covariance in L_X and N_200 at fixed halo mass is not large. The luminosity intercent is 30%, or 2\sigma, lower than determined from the X-ray flux-limited sample of Reiprich & Bohringer (2002), assuming hydrostatic equilibrium. This difference could arise from a combination of Malmquist bias and/or systematic error in hydrostatic mass estimates, both of which are expected. The intercept agrees with that derived by Stanek et al. (2006) using a model for the statistical correspondence between clusters and halos in a WMAP3 cosmology with power spectrum normalization sigma_8 = 0.85. Similar exercises applied to future data sets will allow constraints on the covariance among optical and hot gas properties of clusters at fixed mass.

We present measurements of the excess mass-to-light ratio measured aroundMaxBCG galaxy clusters observed in the SDSS. This red sequence cluster sample includes objects from small groups with masses ranging from ~5x10^{12} to ~10^{15} M_{sun}/h. Using cross-correlation weak lensing, we measure the excess mass density profile above the universal mean \Delta \rho(r) = \rho(r) - \bar{\rho} for clusters in bins of richness and optical luminosity. We also measure the excess luminosity density \Delta l(r) = l(r) - \bar{l} measured in the z=0.25 i-band. For both mass and light, we de-project the profiles to produce 3D mass and light profiles over scales from 25 kpc/ to 22 Mpc/h. From these profiles we calculate the cumulative excess mass M(r) and excess light L(r) as a function of separation from the BCG. On small scales, where \rho(r) >> \bar{\rho}, the integrated mass-to-light profile may be interpreted as the cluster mass-to-light ratio. We find the M/L_{200}, the mass-to-light ratio within r_{200}, scales with cluster mass as a power law with index 0.33+/-0.02. On large scales, where \rho(r) ~ \bar{\rho}, the M/L approaches an asymptotic value independent of cluster richness. For small groups, the mean M/L_{200} is much smaller than the asymptotic value, while for large clusters it is consistent with the asymptotic value. This asymptotic value should be proportional to the mean mass-to-light ratio of the universe <M/L>. We find <M/L>/b^2_{ml} = 362+/-54 h (statistical). There is additional uncertainty in the overall calibration at the ~10% level. The parameter b_{ml} is primarily a function of the bias of the L <~ L_* galaxies used as light tracers, and should be of order unity. Multiplying by the luminosity density in the same bandpass we find \Omega_m/b^2_{ml} = 0.02+/-0.03, independent of the Hubble parameter.

This is the first in a series of papers on the weak lensing effect caused by clusters of galaxies in Sloan Digital Sky Survey. The photometrically selected cluster sample, known as MaxBCG, includes ~130,000 objects between redshift 0.1 and 0.3, ranging in size from small groups to massive clusters. We split the clusters into bins of richness and luminosity and stack the surface density contrast to produce mean radial profiles. The mean profiles are detected over a range of scales, from the inner halo (25 kpc/h) well into the surrounding large scale structure (30 Mpc/h), with a significance of 15 to 20 in each bin. The signal over this large range of scales is best interpreted in terms of the cluster-mass cross-correlation function. We pay careful attention to sources of systematic error, correcting for them where possible. The resulting signals are calibrated to the ~10% level, with the dominant remaining uncertainty being the redshift distribution of the background sources. We find that the profiles scale strongly with richness and luminosity. We find the signal within a given richness bin depends upon luminosity, suggesting that luminosity is more closely correlated with mass than galaxy counts. We split the samples by redshift but detect no significant evolution. The profiles are not well described by power laws. In a subsequent series of papers we invert the profiles to three-dimensional mass profiles, show that they are well fit by a halo model description, measure mass-to-light ratios and provide a cosmological interpretation.

Determining the scaling relations between galaxy cluster observables requires large samples of uniformly observed clusters. We measure the mean X-ray luminosity--optical richness (L_X--N_200) relation for an approximately volume-limited sample of more than 17,000 optically-selected clusters from the maxBCG catalog spanning the redshift range 0.1<z<0.3. By stacking the X-ray emission from many clusters using ROSAT All-Sky Survey data, we are able to measure mean X-ray luminosities to ~10% (including systematic errors) for clusters in nine independent optical richness bins. In addition, we are able to crudely measure individual X-ray emission from ~800 of the richest clusters. Assuming a log-normal form for the scatter in the L_X--N_200 relation, we measure \sigma_\ln{L}=0.86+/-0.03 at fixed N_200. This scatter is large enough to significantly bias the mean stacked relation. The corrected median relation can be parameterized by L_X = (e^\alpha)(N_200/40)^\beta 10^42 h^-2 ergs/s, where \alpha = 3.57+/-0.08 and \beta = 1.82+/-0.05. We find that X-ray selected clusters are significantly brighter than optically-selected clusters at a given optical richness. This selection bias explains the apparently X-ray underluminous nature of optically-selected cluster catalogs.

The distribution of galaxies in position and velocity around the centers of galaxy clusters encodes important information about cluster mass and structure. Using the maxBCG galaxy cluster catalog identified from imaging data obtained in the Sloan Digital Sky Survey, we study the BCG-galaxy velocity correlation function. By modeling its non-Gaussianity, we measure the mean and scatter in velocity dispersion at fixed richness. The mean velocity dispersion increases from 202+/-10 km/s for small groups to more than 854+/-102 km/s for large clusters. We show the scatter to be at most 40.5+/-3.5%, declining to 14.9+/-9.4% in the richest bins. We test our methods in the C4 cluster catalog, a spectroscopic cluster catalog produced from the Sloan Digital Sky Survey DR2 spectroscopic sample, and in mock galaxy catalogs constructed from N-body simulations. Our methods are robust, measuring the scatter to well within one-sigma of the true value, and the mean to within 10%, in the mock catalogs. By convolving the scatter in velocity dispersion at fixed richness with the observed richness space density function, we measure the velocity dispersion function of the maxBCG galaxy clusters. Although velocity dispersion and richness do not form a true mass-observable relation, the relationship between velocity dispersion and mass is theoretically well characterized and has low scatter. Thus our results provide a key link between theory and observations up to the velocity bias between dark matter and galaxies.

We interpret and model the statistical weak lensing measurements around 130,000 groups and clusters of galaxies in the Sloan Digital Sky Survey presented by Sheldon et al. 2007 (Paper I). We present non-parametric inversions of the 2D shear profiles to the mean 3D cluster density and mass profiles in bins of both optical richness and cluster i-band luminosity. We correct the inferred 3D profiles for systematic effects, including non-linear shear and the fact that cluster halos are not all precisely centered on their brightest galaxies. We also model the measured cluster shear profile as a sum of contributions from the brightest central galaxy, the cluster dark matter halo, and neighboring halos. We infer the relations between mean cluster virial mass and optical richness and luminosity over two orders of magnitude in cluster mass; the virial mass at fixed richness or luminosity is determined with a precision of 13% including both statistical and systematic errors. We also constrain the halo concentration parameter and halo bias as a function of cluster mass; both are in good agreement with predictions of LCDM models. The methods employed here will be applicable to deeper, wide-area optical surveys that aim to constrain the nature of the dark energy, such as the Dark Energy Survey, the Large Synoptic Survey Telescope and space-based surveys.

We present several cases of optical observations during gamma-ray bursts (GRBs) which resulted in prompt limits but no detection of optical emission. These limits constrain the prompt optical flux densities and the optical brightness relative to the gamma-ray emission. The derived constraints fall within the range of properties observed in GRBs with prompt optical detections, though at the faint end of optical/gamma flux ratios. The presently accessible prompt optical limits do not require a different set of intrinsic or environmental GRB properties, relative to the events with prompt optical detections.

We report on follow-up observations of the GRB 060927 using the ROTSE-IIIa telescope and a suite of larger aperture ground-based telescopes. An optical afterglow was detected 20 s after the burst, the earliest rest-frame detection of optical emission from any GRB. Spectroscopy performed with the VLT about 13 hours after the trigger shows a continuum break at lambda ~ 8070 A produced by neutral hydrogen absorption at z~5.6. We also detect an absorption line at 8158 A which we interpret as SiII at z=5.467. Hence, GRB 060927 is the second most distant GRB with a spectroscopically measured redshift. The shape of the red wing of the spectral break can be fitted by a damped Lyalpha profile with a column density with log(N_HI/cm^-2) ~ 22.5. We discuss the implications of this work for the use of GRBs as probes of the end of the dark ages and draw three main conclusions: i) GRB afterglows originating from z>6 should be relatively easy to detect from the ground, but rapid NIR monitoring is necessary to ensure that they are found; ii) The presence of large HI column densities in some GRBs host galaxies at z>5 makes the use of GRBs to probe the reionization epoch via spectroscopy of the red damping wing challenging; iii) GRBs appear crucial to locate typical star-forming galaxies at z>5 and therefore the type of galaxies responsible for the reionization of the universe.

We introduce a framework for describing the halo selection function of optical cluster finders. We treat the problem as being separable into a term that describes the intrinsic galaxy content of a halo (the Halo Occupation Distribution, or HOD) and a term that captures the effects of projection and selection by the particular cluster finding algorithm. Using mock galaxy catalogs tuned to reproduce the luminosity dependent correlation function and the empirical color-density relation measured in the SDSS, we characterize the maxBCG algorithm applied by Koester et al. to the SDSS galaxy catalog. We define and calibrate measures of completeness and purity for this algorithm, and demonstrate successful recovery of the underlying cosmology and HOD when applied to the mock catalogs. We identify principal components -- combinations of cosmology and HOD parameters -- that are recovered by survey counts as a function of richness, and demonstrate that percent-level accuracies are possible in the first two components, if the selection function can be understood to ~15% accuracy.

We perform a maximum likelihood analysis of the cluster abundance measured in the SDSS using the maxBCG cluster finding algorithm. Our analysis is aimed at constraining the power spectrum normalization $\sigma_8$, and assumes flat cosmologies with a scale invariant spectrum, massless neutrinos, and CMB and supernova priors Omega_m*h^2=0.128+/-0.01 and h=0.72+/-0.05 respectively. Following the method described in the companion paper Rozo et al. 2007, we derive \sigma_8=0.92+/-0.10$ (1-sigma) after marginalizing over all major systematic uncertainties. We place strong lower limits on the normalization, sigma_8>0.76 (95% CL) (>0.68 at 99% CL). We also find that our analysis favors relatively low values for the slope of the Halo Occupation Distribution (HOD), alpha=0.83+/-0.06. The uncertainties of these determinations will substantially improve upon completion of an ongoing campaign to estimate dynamical, weak lensing, and X-ray cluster masses in the SDSS maxBCG cluster sample.

We present a catalog of galaxy clusters selected using the maxBCG redsequence method from Sloan Digital Sky Survey photometric data. This catalog includes 13,823 clusters with velocity dispersions greater than 400 km/s, and is the largest galaxy cluster catalog assembled to date. They are selected in an approximately volume-limited way from a 0.5 Gpc^3 region covering 7500 square degrees of sky between redshifts 0.1 and 0.3. (ABRIGDED)

Measurements of galaxy cluster abundances, clustering properties, and mass to- light ratios in current and future surveys can provide important cosmological constraints. Digital wide-field imaging surveys, the recently-demonstrated fidelity of red-sequence cluster detection techniques, and a new generation of realistic mock galaxy surveys provide the means for construction of large, cosmologicallyinteresting cluster samples, whose selection and properties can be understood in unprecedented depth. We present the details of the "maxBCG" algorithm, a cluster-detection technique tailored to multi-band CCD-imaging data. MaxBCG primarily relies on an observational cornerstone of massive galaxy clusters: they are marked by an overdensity of bright, uniformly red galaxies. This detection scheme also exploits classical brightest cluster galaxies (BCGs), which are often found at the center of these same massive clusters. (ABRIDGED)

The robotic ROTSE-III telescope network detected prompt optical emission contemporaneous with the gamma-ray emission of Swift events GRB051109A and GRB051111. Both datasets have continuous coverage at high signal-to-noise levels from the prompt phase onwards, thus the early observations are readily compared to the Swift XRT and BAT high energy detections. In both cases, the optical afterglow is established, declining steadily during the prompt emission. For GRB051111, there is evidence of an excess optical component during the prompt emission. The component is consistent with the flux spectrally extrapolated from the gamma-rays, using the gamma-ray spectral index. A compilation of spectral information from previous prompt detections shows that such a component is unusual. The existence of two prompt optical components - one connected to the high-energy emission, the other to separate afterglow flux, as indicated in GRB051111 - is not compatible with a simple ``external-external'' shock model for the GRB and its afterglow.

The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi-step process which involves data from three different telescopes: the 1.0-m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which was used to establish the SDSS standard star network); the SDSS 0.5-m Photometric Telescope (PT) at the Apache Point Observatory (APO), New Mexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and the SDSS 2.5-m telescope at APO (which obtains the imaging data for the SDSS proper). In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data from the single-CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0-m and the PT). We also describe transformation equations that convert photometry on the USNO-1.0m u'g'r'i'z' system to photometry the SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the semi-automated PT factory, which runs MTPIPE in the day-to-day standard SDSS operations at Fermilab. Finally, we discuss the use of MTPIPE in current SDSS-related projects, including the Southern u'g'r'i'z' Standard Star project, the u'g'r'i'z' Open Star Clusters project, and the SDSS extension (SDSS-II).

We use data from the Northern Sky Variability Survey (NSVS), obtained from the first generation Robotic Optical Transient Search Experiment (ROTSE-I), to identify and study RR Lyrae variable stars in the solar neighborhood. We initially identified 1197 RRab (RR0) candidate stars brighter than the ROTSE median magnitude V = 14. Periods, amplitudes, and mean V magnitudes are determined for a subset of 1188 RRab stars with well defined light curves. Metallicities are determined for 589 stars by the Fourier parameter method and by the relationship between period, amplitude, and [Fe/H]. We comment upon the difficulties of clearly classifying RRc (RR1) variables in the NSVS dataset. Distances to the RRab stars are calculated using an adopted luminosity-metallicity relation with corrections for interstellar extinction. The 589 RRab stars in our final sample are used to study the properties of the RRab population within 5 kpc of the Sun. The Bailey diagram of period versus amplitude shows that the largest component of this sample belongs to Oosterhoff type I. Metal-rich ([Fe/H] > -1) RRab stars appear to be associated with the Galactic disk. Our metal-rich RRab sample may include a thin disk as well as a thick disk population, although the uncertainties are too large to establish this. There is some evidence among the metal-rich RRab stars for a decline in scale height with increasing [Fe/H], as was found by Layden (1995). The distribution of RRab stars with -1 < [Fe/H] < -1.25 indicates that within this metallicity range the RRab stars are a mixture of stars belonging to halo and disk populations.

The ROTSE-IIIc telescope at the H.E.S.S. site, Namibia, obtained the earliest detection of optical emission from a Gamma-Ray Burst (GRB), beginning only 21.8 s from the onset of Swift GRB 050801. The optical lightcurve does not fade or brighten significantly over the first ~250 s, after which there is an achromatic break and the lightcurve declines in typical power-law fashion. The Swift/XRT also obtained early observations starting at 69 s after the burst onset. The X-ray lightcurve shows the same features as the optical lightcurve. These correlated variations in the early optical and X-ray emission imply a common origin in space and time. This behavior is difficult to reconcile with the standard models of early afterglow emission.

We present the unfiltered ROTSE-III light curve of the optical transient associated with GRB 050319 beginning 4 s after the cessation of gamma-ray activity. We fit a power-law function to the data using the revised trigger time given by Chincarini et al. (2005), and a smoothly broken power-law to the data using the original trigger disseminated through the GCN notices. Including the RAPTOR data from Wozniak et al. (2005), the best fit power-law indices are alpha=-0.854 (+/- 0.014) for the single power-law and alpha_1=-0.364 (+/- 0.020), alpha_2= -0.881 (+/- 0.030), with a break at t_b = 418 (+/- 30) s for the smoothly broken fit. We discuss the fit results with emphasis placed on the importance of knowing the true start time of the optical transient for this multi-peaked burst. As Swift continues to provide prompt GRB locations, it becomes more important to answer the question, "when does the afterglow begin" to correctly interpret the light curves.

In this work we describe a large new sample of contact binary stars extracted in a uniform manner from sky patrol data taken by the ROTSE-I telescope. Extensive ROTSE-I light curve data is combined with J, H, and K band near-infrared data taken from the Two Micron All-Sky Survey (2MASS) to add color information. Contact binaries candidates are selected using the observed period-color relation. Candidates are confirmed by visual examination of the light curves. To enhance the utility of this catalog, we derive a new J-H period-color-luminosity relation and use this to estimate distances for the entire catalog. From these distance estimates we derive an estimated contact binary space density of (1.7 +/- 0.6) x 10^-5 pcs^-3.

Contact binary stars are common variable stars which are all believed to emit relatively large fluxes of x-rays. In this work we combine a large new sample of contact binary stars derived from the ROTSE-I telescope with x-ray data from the ROSAT All-Sky Survey (RASS) to estimate the x-ray volume emissivity of contact binary stars in the galaxy. We obtained x-ray fluxes for 140 contact binaries from the RASS, as well as 2 additional stars observed by the XMM-Newton observatory. From these data we confirm the emission of x-rays from all contact binary systems, with typical luminosities of approximately 1.0 x 10^30 erg s^-1. Combining calculated luminosities with an estimated contact binary space density, we find that contact binaries do not have strong enough x-ray emission to account for a significant portion of the galactic x-ray background.

We present lightcurves of the afterglow of GRB050502A, including very early data at t-t_{GRB} < 60s. The lightcurve is composed of unfiltered ROTSE-IIIb optical observations from 44s to 6h post-burst, R-band MDM observations from 1.6 to 8.4h post-burst, and PAIRITEL J H K_s observations from 0.6 to 2.6h post-burst. The optical lightcurve is fit by a broken power law, where t^{alpha} steepens from alpha = -1.13 +- 0.02 to alpha = -1.44 +- 0.02 at \~5700s. This steepening is consistent with the evolution expected for the passage of the cooling frequency nu_c through the optical band. Even in our earliest observation at 44s post-burst, there is no evidence that the optical flux is brighter than a backward extrapolation of the later power law would suggest. The observed decay indices and spectral index are consistent with either an ISM or a Wind fireball model, but slightly favor the ISM interpretation. The expected spectral index in the ISM interpretation is consistent within 1 sigma with the observed spectral index beta = -0.8 +- 0.1; the Wind interpretation would imply a slightly (~2 sigma) shallower spectral index than observed. A small amount of dust extinction at the source redshift could steepen an intrinsic spectrum sufficiently to account for the observed value of beta. In this picture, the early optical decay, with the peak at or below 4.7e14 Hz at 44s, requires very small electron and magnetic energy partitions from the fireball.

The ROTSE-IIIa telescope at Siding Spring Observatory, Australia, detected prompt optical emission from Swift GRB 050401. In this letter, we present observations of the early optical afterglow, first detected by the ROTSE-IIIa telescope 33 s after the start of gamma-ray emission, contemporaneous with the brightest peak of this emission. This GRB was neither exceptionally long nor bright. This is the first prompt optical detection of a GRB of typical duration and luminosity. We find that the early afterglow decay does not deviate significantly from the power-law decay observable at later times, and is uncorrelated with the prompt gamma-ray emission. We compare this detection with the other two GRBs with prompt observations, GRB 990123 and GRB 041219a. All three bursts exhibit quite different behavior at early times.

Imaging data from the Sloan Digital Sky Survey is used to measure the empirical size-richness relation for a large sample of galaxy clusters. Using population subtraction methods, we determine the radius at which the cluster galaxy number density is 200/Omega_m times the mean galaxy density, without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are unbiased on Mpc scales, this galaxy-density-based R_200 reflects the characteristic radii of clusters. We measure the scaling of this characteristic radius with richness over an order of magnitude in cluster richness, from rich clusters to poor groups. We use this information to examine the radial profiles of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy distribution decreases with richness and is systematically lower than the concentrations measured for dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of the cluster luminosity function, and find that it is well matched by a Schechter function for galaxies brighter than M_r = -18 only after the central galaxy has been removed. We find that the luminosity function varies with richness and with distance from the cluster center, underscoring the importance of using an aperture that scales with cluster mass to compare physically equivalent regions of these different systems. We note that the lowest richness systems in our catalog have properties consistent with those expected of the earliest-forming halos; our cluster-finding algorithm, in addition to reliably finding clusters, may be efficient at finding fossil groups.

We present a new non-parametric method for determining mean 3D density and mass profiles from weak lensing measurements around stacked samples of galaxies or clusters, that is, from measurement of the galaxy-shear or cluster-shear correlation functions. Since the correlation function is statistically isotropic, this method evades problems, such as projection of large-scale structure along the line of sight or halo asphericity, that complicate attempts to infer masses from weak lensing measurements of individual objects. We demonstrate the utility of this method in measuring halo profiles, galaxy-mass and cluster-mass cross-correlation functions, and cluster virial masses. We test this method on an N-body simulation and show that it correctly and accurately recovers the 3D density and mass profiles of halos. We find no evidence of problems due to a mass sheet degeneracy in the simulation. Cross-correlation lensing provides a powerful method for calibrating the mass-observable relation for use in measurement of the cluster mass function in large surveys. It can also be used on large scales to measure and remove the halo bias and thereby provide a direct measurement of Omega_m * sigma_8.

We present the results of a search for untriggered gamma-ray burst (GRB) afterglows with the Robotic Optical Transient Search Experiment-III (ROTSE-III) telescope array. This search covers observations from September 2003 to March 2005. We have an effective coverage of 1.74 deg^2 yr for rapidly fading transients that remain brighter than ~ 17.5 magnitude for more than 30 minutes. This search is the first large area survey to be able to detect typical untriggered GRB afterglows. Our background rate is very low and purely astrophysical. We have found 4 previously unknown cataclysmic variables (CVs) and 1 new flare star. We have not detected any candidate afterglow events or other unidentified transients. We can place an upper limit on the rate of fading optical transients with quiescent counterparts dimmer than ~ 20th magnitude at a rate of less than 1.9 deg^-2 yr-1 with 95% confidence. This places limits on the optical characteristics of off-axis (orphan) GRB afterglows. As a byproduct of this search, we have an effective ~ 52 deg^2 yr of coverage for very slowly decaying transients, such as CVs. This implies an overall rate of outbursts from high galactic latitude CVs of 0.1 deg^2 yr^-1.

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We describe the application of the `shapelet' linear decomposition of galaxy images to multi-wavelength morphological classification using the $u,g,r,i,$ and $z$-band images of 1519 galaxies from the Sloan Digital Sky Survey. We utilize elliptical shapelets to remove to first-order the effect of inclination on morphology. After decomposing the galaxies we perform a principal component analysis on the shapelet coefficients to reduce the dimensionality of the spectral morphological parameter space. We give a description of each of the first ten principal component's contribution to a galaxy's spectral morphology. We find that galaxies of different broad Hubble type separate cleanly in the principal component space. We apply a mixture of Gaussians model to the 2-dimensional space spanned by the first two principal components and use the results as a basis for classification. Using the mixture model, we separate galaxies into three classes and give a description of each class's physical and morphological properties. We find that the two dominant mixture model classes correspond to early and late type galaxies, respectively. The third class has, on average, a blue, extended core surrounded by a faint red halo, and typically exhibits some asymmetry. We compare our method to a simple cut on $u-r$ color and find the shapelet method to be superior in separating galaxies. Furthermore, we find evidence that the $u-r=2.22$ decision boundary may not be optimal for separation between early and late type galaxies, and suggest that the optimal cut may be $u-r \sim 2.4$.

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

Weak gravitational lensing provides a unique method to directly map the dark matter in the universe and measure cosmological parameters. Current weak lensing surveys are limited by the atmospheric seeing from the ground and by the small field of view of existing space telescopes. We study how a future wide-field space telescope can measure the lensing power spectrum and skewness, and set constraints on cosmological parameters. The lensing sensitivity was calculated using detailed image simulations and instrumental specifications studied in earlier papers in this series. For instance, the planned SuperNova/Acceleration Probe (SNAP) mission will be able to measure the matter density parameter Omega_m and the dark energy equation of state parameter w with precisions comparable and nearly orthogonal to those derived with SNAP from supernovae. The constraints degrade by a factor of about 2 if redshift tomography is not used, but are little affected if the skewness only is dropped. We also study how the constraints on these parameters depend upon the survey geometry and define an optimal observing strategy.

We study the accuracy with which weak lensing measurements could be made from a future space-based survey, predicting the subsequent precisions of 3-dimensional dark matter maps, projected 2-dimensional dark matter maps, and mass-selected cluster catalogues. As a baseline, we use the instrumental specifications of the Supernova/Acceleration Probe (SNAP) satellite. We first compute its sensitivity to weak lensing shear as a function of survey depth. Our predictions are based on detailed image simulations created using `shapelets', a complete and orthogonal parameterization of galaxy morphologies. We incorporate a realistic redshift distribution of source galaxies, and calculate the average precision of photometric redshift recovery using the SNAP filter set to be Delta z=0.034. The high density of background galaxies resolved in a wide space-based survey allows projected dark matter maps with a rms sensitivity of 3% shear in 1 square arcminute cells. This will be further improved using a proposed deep space-based survey, which will be able to detect isolated clusters using a 3D lensing inversion techniques with a 1 sigma mass sensitivity of approximately 10^13 solar masses at z~0.25. Weak lensing measurements from space will thus be able to capture non-Gaussian features arising from gravitational instability and map out dark matter in the universe with unprecedented resolution.

We present galaxy-galaxy lensing measurements over scales 0.025 to 10 Mpc/h in the Sloan Digital Sky Survey. Using a flux-limited sample of 127,001 lens galaxies with spectroscopic redshifts and mean luminosity <L> = L_* and 9,020,388 source galaxies with photometric redshifts, we invert the lensing signal to obtain the galaxy-mass correlation function xi_{gm}. We find xi_{gm} is consistent with a power-law, xi_{gm} = (r/r_0)^{-gamma}, with best-fit parameters gamma = 1.79 +/- 0.06 and r_0 = (5.4+/-0.7)(0.27/Omega_m)^{1/gamma} Mpc/h. At fixed separation, the ratio xi_{gg}/xi_{gm} = b/r where b is the bias and r is the correlation coefficient. Comparing to the galaxy auto-correlation function for a similarly selected sample of SDSS galaxies, we find that b/r is approximately scale independent over scales 0.2-6.7 Mpc/h, with mean <b/r> = (1.3+/-0.2)(Omega_m/0.27). We also find no scale dependence in b/r for a volume limited sample of luminous galaxies (-23.0 < M_r < -21.5). The mean b/r for this sample is <b/r>_{Vlim} = (2.0+/-0.7)(Omega_m/0.27). We split the lens galaxy sample into subsets based on luminosity, color, spectral type, and velocity dispersion, and see clear trends of the lensing signal with each of these parameters. The amplitude and logarithmic slope of xi_{gm} increases with galaxy luminosity. For high luminosities (L ~5 L_*), xi_{gm} deviates significantly from a power law. These trends with luminosity also appear in the subsample of red galaxies, which are more strongly clustered than blue galaxies.

The Northern Sky Variability Survey (NSVS) is a temporal record of the sky over the optical magnitude range from 8 to 15.5. It was conducted in the course of the first generation Robotic Optical Transient Search Experiment (ROTSE-I) using a robotic system of four co-mounted unfiltered telephoto lenses equipped with CCD cameras. The survey was conducted from Los Alamos, NM, and primarily covers the entire northern sky. Some data in southern fields between declinations 0 and -38 deg is also available, although with fewer epochs and noticeably lesser quality. The NSVS contains light curves for approximately 14 million objects. With a one year baseline and typically 100-500 measurements per object, the NSVS is the most extensive record of stellar variability across the bright sky available today. In a median field, bright unsaturated stars attain a point to point photometric scatter of ~0.02 mag and position errors within 2 arcsec. At Galactic latitudes |b| < 20 deg the data quality is limited by severe blending due to ~14 arcsec pixel size. We present basic characteristics of the data set and describe data collection, analysis, and distribution. All NSVS photometric measurements are available for on-line public access from the Sky Database for Objects in Time-Domain (SkyDOT; http://skydot.lanl.gov) at LANL. Copies of the full survey photometry may also be requested on tape.

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

The ROTSE-IIIa telescope and the SSO-40 inch telescope, both located at Siding Spring Observatory, imaged the early time afterglow of GRB 030418. In this report we present observations of the early afterglow, first detected by the ROTSE-IIIa telescope 211 s after the start of the burst, and only 76 s after the end of the gamma-ray activity. We detect optical emission that rises for ~600 s, slowly varies around R=17.3 mag for ~1400 s, and then fades as a power law of index alpha=-1.36. Additionally, the ROTSE-IIIb telescope, located at McDonald Observatory, imaged the early time afterglow of GRB 030723. The behavior of this light curve was qualitatively similar to that of GRB 030418, but two magnitudes dimmer. These two afterglows are dissimilar to other afterglows such as GRB 990123 and GRB 021211. We investigate whether the early afterglow can be attributed to a synchrotron break in a cooling synchrotron spectrum as it passes through the optical band, but find this model is unable to accurately describe the early light curve. We present a simple model for gamma-ray burst emission emerging from a wind medium surrounding a massive progenitor star. This model provides an effective description of the data, and suggests that the rise of the afterglow can be ascribed to extinction in the local circumburst environment. In this interpretation, these events provide further evidence for the connection between gamma-ray bursts and the collapse of massive stars.

We describe application of the `shapelet' linear decomposition of galaxy images to morphological classification using images of $\sim$ 3000 galaxies from the Sloan Digital Sky Survey. After decomposing the galaxies we perform a principal component analysis to reduce the number of dimensions of the shapelet space to nine. We find that each of these nine principal components contains unique morphological information, and give a description of each principal component's contribution to a galaxy's morphology. We find that galaxies of differing Hubble type separate cleanly in the shapelet space. We apply a Gaussian mixture model to the 9-dimensional space spanned by the principal components and use the results as a basis for classification. Using the mixture model, we separate galaxies into seven classes and give a description of each class's physical and morphological properties. We find that several of the mixture model classes correlate well with the traditional Hubble types both in their morphology and their physical parameters (e.g., color, velocity dispersions, etc.). In addition, we find an additional class of late-type morphology but with high velocity dispersions and very blue color; most of these galaxies exhibit post-starburst activity. This method provides an objective and quantitative alternative to traditional and subjective visual classification.

Using two identical telescopes at widely separated longitudes, the ROTSE-III network observed decaying emission from the remarkably bright afterglow of GRB 030329. In this report we present observations covering 56% of the period from 1.5-47 hours after the burst. We find that the light curve is piecewise consistent with a powerlaw decay. When the ROTSE-III data are combined with data reported by other groups, there is evidence for five breaks within the first 20 hours after the burst. Between two of those breaks, observations from 15.9-17.1 h after the burst at 1-s time resolution with McDonald Observatory's 2.1-m telescope reveal no evidence for fluctuations or deviations from a simple power law. Multiple breaks may indicate complex structure in the jet. There are also two unambiguous episodes at 23 and 45 hours after the burst where the intensity becomes consistent with a constant for several hours, perhaps indicating multiple injections of energy into the GRB/afterglow system.

The distribution in color and absolute magnitude is presented for 55,158 galaxies taken from the Sloan Digital Sky Survey in the redshift range $0.08<z<0.12$, as a function of galaxy number overdensity in a line-of-sight cylinder of transverse radius $1 h^{-1} \mathrm{Mpc}$. In all environments, bulge-dominated galaxies (defined to be those with radial profiles best fit with large S\'ersic indices) form a narrow, well defined color--magnitude relation. Although the most luminous galaxies reside preferentially in the highest density regions, there is only a barely detectable variation in the color (zero-point) or slope of the color--magnitude relation ($<0.02 \mathrm{mag}$ in $^{0.1}[g-r]$ or $[B-V]$). These results constrain variations with environmental density in the ages or metallicities of typical bulge-dominated galaxies to be under 20 percent.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present a catalog of 799 clusters of galaxies in the redshift range z_est = 0.05 - 0.3 selected from ~400 deg^2 of early SDSS commissioning data along the celestial equator. The catalog is based on merging two independent selection methods -- a color-magnitude red-sequence maxBCG technique (B), and a Hybrid Matched-Filter method (H). The BH catalog includes clusters with richness \Lambda >= 40 (Matched-Filter) and N_gal >= 13 (maxBCG), corresponding to typical velocity dispersion of \sigma_v >~ 400 km s^{-1} and mass (within 0.6 h^{-1) Mpc radius) >~ 5*10^{13} h^{-1} M_sun. This threshold is below Abell richness class 0 clusters. The average space density of these clusters is 2*10^{-5} h^3 Mpc^{-3}. All NORAS X-ray clusters and 53 of the 58 Abell clusters in the survey region are detected in the catalog; the 5 additional Abell clusters are detected below the BH catalog cuts. The cluster richness function is determined and found to exhibit a steeply decreasing cluster abundance with increasing richness. We derive observational scaling relations between cluster richness and observed cluster luminosity and cluster velocity dispersion; these scaling relations provide important physical calibrations for the clusters. The catalog can be used for studies of individual clusters, for comparisons with other sources such as X-ray clusters and AGNs, and, with proper correction for the relevant selection functions, also for statistical analyses of clusters.

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on the sky, we detect about 3000 satellites with absolute blue magnitudes going down to M_B = -14; most of the satellites have M_B=-16 to -18, comparable to the magnitudes of M32 and the Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the line-of-sight velocity dispersion of satellites declines with distance to the primary. For an L* galaxy the r.m.s. line-of-sight velocity changes from ~120 km/s at 20 kpc to ~60 km/s at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern cosmological models predict that the density of dark matter in the peripheral parts of galaxies declines as rho_DM propto r^{-3}. Thus, for the first time we find direct observational evidence of the density decline predicted by cosmological models; we also note that this result contradicts alternative theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 100 kpc scales with the absolute magnitude of the central galaxy as sigma propto L^{0.3}; this is very close to the Tully-Fisher relation for normal spiral galaxies.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

The observation of a prompt optical flash from GRB990123 convincingly demonstrated the value of autonomous robotic telescope systems. Pursuing a program of rapid follow-up observations of gamma-ray bursts, the Robotic Optical Transient Search Experiment (ROTSE) has developed a next-generation instrument, ROTSE-III, that will continue the search for fast optical transients. The entire system was designed as an economical robotic facility to be installed at remote sites throughout the world. There are seven major system components: optics, optical tube assembly, CCD camera, telescope mount, enclosure, environmental sensing & protection and data acquisition. Each is described in turn in the hope that the techniques developed here will be useful in similar contexts elsewhere.

We report the discovery of a highly luminous optical transient (OT), SDSS J124602.54+011318.8, associated with a galaxy at a redshift of 0.385. In this paper we consider the possibility that the OT may be a GRB afterglow. Three sets of images and two sets of spectra were obtained as part of the normal operations of the Sloan Digital Sky Survey (SDSS). In the first two image sets, observed two nights apart, the object appears as a point source at $r^{*}\approx 17$. The third image set, observed about 410 days later, shows an extended source which is more than 2.5 magnitudes fainter. The spectra were observed about 400 and 670 days after the first two image sets, and both show an apparently normal galaxy at a redshift of 0.385. Associating the OT with the galaxy, the absolute magnitude was $M_{r^*}=-24.8$, which is over 4 magnitudes brighter than the most luminous supernova ever measured. The spectral energy distributions of the galaxy-subtracted OT derived from the first two image sets are well-fit by single power-laws with indices of $\beta_{\nu}=-0.92$ and -1.29 respectively, similar to most GRB afterglows. Based upon the luminosity of the OT, non-detections in contemporaneous ROTSE-I images, and the change in spectral slope, the OT, if an afterglow, was likely discovered early during a ``plateau'' or slowly-fading phase. The discovery of a GRB afterglow at this stage of the SDSS is consistent with expectations, but only if the optical emission is much less strongly beamed than the gamma-rays. We emphasize that other explanations for the OT cannot be ruled out; a recent follow-up study by [galyam02] provides strong evidence that this source is in fact an unusual AGN.

We present an untriggered search for optical bursts with the ROTSE-I telephoto array. Observations were taken which monitor an effective 256 square degree field continuously over 125 hours to m_{ROTSE}=15.7. The uniquely large field, moderate limiting magnitude and fast cadence of $\sim$10 minutes permits transient searches in a new region of sensitivity. Our search reveals no candidate events. To quantify this result, we simulate potential optical bursts with peak magnitude, m_{p}, at t=10 s, which fade as f=(\frac{t}{t_{0}}) ^{\alpha_{t}}, where \alpha_t < 0. Simple estimates based on observational evidence indicate that a search of this sensitivity begins to probe the possible region occupied by GRB orphan afterglows. Our observing protocol and image sensitivity result in a broad region of high detection efficiency for light curves to the bright and slowly varying side of a boundary running from [\alpha_{t},m_{p}]=[-2.0,6.0] to [-0.3,13.2]. Within this region, the integrated rate of brief optical bursts is less than 1.1\times 10^{-8} {\rm s}^{-1} {\rm deg}^{-2}. At $\sim$22 times the observed GRB rate from BATSE, this suggests a limit on \frac{\theta_{opt}}{\theta_{\gamma}}\lesssim 5 where \theta_{opt} and \theta_{\gamma} are the optical and gamma-ray collimation angles, respectively. Several effects might explain the absence of optical bursts, and a search of the kind described here but more sensitive by about 4 magnitudes should offer a more definitive probe.

We report on the current operating status of the ROTSE-IIIa telescope, currently undergoing testing at Los Alamos National Laboratories in New Mexico. It will be shipped to Siding Spring Observatory, Australia, in first quarter 2002. ROTSE-IIIa has been in automated observing mode since early October, 2001, after completing several weeks of calibration and check-out observations. Calibrated lists of objects in ROTSE-IIIa sky patrol data are produced routinely in an automated pipeline, and we are currently automating analysis procedures to compile these lists, eliminate false detections, and automatically identify transient and variable objects. The manual application of these procedures has already led to the detection of a nova that rose over six magnitudes in two days to a maximum detected brightness of m_R~13.9 and then faded two magnitudes in two weeks. We also readily identify variable stars, includings those suspected to be variables from the Sloan Digital Sky Survey. We report on our system to allow public monitoring of the telescope operational status in real time over the WWW.

Galaxy masses can be estimated by a variety of methods; each applicable in different circumstances, and each suffering from different systematic uncertainties. Confirmation of results obtained by one technique with analysis by another is particularly important. Recent SDSS weak lensing measurements of the projected-mass correlation function reveal a linear relation between galaxy luminosities and the depth of their dark matter halos (measured on 260 \hinv kpc scales). In this work we use an entirely independent dynamical method to confirm these results. We begin by assembling a sample of 618 relatively isolated host galaxies, surrounded by a total of 1225 substantially fainter satellites. We observe the mean dynamical effect of these hosts on the motions of their satellites by assembling velocity difference histograms. Dividing the sample by host properties, we find significant variations in satellite velocity dispersion with host luminosity. We quantify these variations using a simple dynamical model, measuring \mtsd a dynamical mass within 260 \hinv kpc. The appropriateness of this mass reconstruction is checked by conducting a similar analysis within an N-body simulation. Comparison between the dynamical and lensing mass-to-light scalings shows reasonable agreement, providing some quantitative confirmation for the lensing results.

We present the 158 standard stars that define the u'g'r'i'z' photometric system. These stars form the basis for the photometric calibration of the Sloan Digital Sky Survey (SDSS). The defining instrument system and filters, the observing process, the reduction techniques, and the software used to create the stellar network are all described. We briefly discuss the history of the star selection process, the derivation of a set of transformation equations for the UBVRcIc system, and plans for future work.

We report on the discovery of an outburst from a previously unknown high- latitude CV we designate ROTSE3 J015118.59-022300.1. The object was first detected on 2001 Oct 13.291 at m-ROTSE3=14.71+-0.06, after which it faded by more than 2 magnitudes over 13 d. Spectral analysis confirms a Galactic location, despite unusual properties that render classification difficult.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

We present moderate resolution Keck spectroscopy of quasars at z=5.82, 5.99 and 6.28, discovered by the Sloan Digital Sky Survey (SDSS). We find that the Ly Alpha absorption in the spectra of these quasars evolves strongly with redshift. To z~5.7, the Ly Alpha absorption evolves as expected from an extrapolation from lower redshifts. However, in the highest redshift object, SDSSp J103027.10+052455.0 (z=6.28), the average transmitted flux is 0.0038+-0.0026 times that of the continuum level over 8450 A < lambda < 8710 A (5.95<z(abs)<6.16), consistent with zero flux. Thus the flux level drops by a factor of >150, and is consistent with zero flux in the Ly Alpha forest region immediately blueward of the Ly Alpha emission line, compared with a drop by a factor of ~10 at z(abs)~5.3. A similar break is seen at Ly Beta; because of the decreased oscillator strength of this transition, this allows us to put a considerably stronger limit, tau(eff) > 20, on the optical depth to Ly Alpha absorption at z=6. This is a clear detection of a complete Gunn-Peterson trough, caused by neutral hydrogen in the intergalactic medium. Even a small neutral hydrogen fraction in the intergalactic medium would result in an undetectable flux in the Ly Alpha forest region. Therefore, the existence of the Gunn-Peterson trough by itself does not indicate that the quasar is observed prior to the reionization epoch. However, the fast evolution of the mean absorption in these high-redshift quasars suggests that the mean ionizing background along the line of sight to this quasar has declined significantly from z~5 to 6, and the universe is approaching the reionization epoch at z~6.

We present new measurements of scaling laws relating the luminosity of galaxies to the amplitude and shape of their dark matter halos. Early imaging and spectroscopic data from the Sloan Digital Sky Survey are used to make weak lensing measurements of the surface mass density contrast Delta Sigma_+ around classes of lens objects. This surface mass density contrast as a function of radius is a measure of the galaxy-mass correlation function (GMCF). Because spectroscopic redshifts are available for all lens objects, the mass and distance scales are well constrained. The GMCF measured around ~31,000 lenses is well fit by a power law of the form Delta Sigma_+ = (2.5+0.7-0.6) (R/1 Mpc)^{-0.8+-0.2} h M_sun pc^-2. We compare this GMCF to galaxy luminosity, type, and environment, and find that it varies strongly with all three. We quantify these variations by comparing the normalization of a fit to the inner 260 h^-1 kpc, M_260, to the galaxy luminosity. While M_260 is not strongly related to luminosity in bluest band u', there is a simple, linear relation between M_260 and luminosity in redder bands (g', r', i', and z'). We test the universality of these mass-to-light scalings by independently measuring them for spiral and elliptical galaxies,and for galaxies in a variety of environments. We find remarkable consistency in these determinations in the red bands, especially i' and z'. This consistency across a wide range of systems suggests that the measured scaling represents an excellent cosmic average, and that the integrated star formation history of galaxies is strongly related to the dark matter environments in which they form.

The X-ray nova XTE J1118+480 exhibited two outbursts in the early part of 2000. As detected by the Rossi X-ray Timing Explorer (RXTE), the first outburst began in early January and the second began in early March. Routine imaging of the northern sky by the Robotic Optical Transient Search Experiment (ROTSE) shows the optical counterpart to XTE J1118+480 during both outbursts. These data include over 60 epochs from January to June 2000. A search of the ROTSE data archives reveal no previous optical outbursts of this source in selected data between April 1998 and January 2000. While the X-ray to optical flux ratio of XTE J1118+480 was low during both outbursts, we suggest that they were full X-ray novae and not mini-outbursts based on comparison with similar sources. The ROTSE measurements taken during the March 2000 outburst also indicate a rapid rise in the optical flux that preceded the X-ray emission measured by the RXTE by approximately 10 days. Using these results, we estimate a pre-outburst accretion disk inner truncation radius of 1.2 x 10^4 Schwarzschild radii.

Gamma-ray bursts of short duration may harbor vital clues to the range of phenomena producing bursts. However, recent progress from the observation of optical counterparts has not benefitted the study of short bursts. We have searched for early optical emission from six gamma-ray bursts using the ROTSE-I telephoto array. Three of these events were of short duration, including GRB 980527 which is among the brightest short bursts yet observed. The data consist of unfiltered CCD optical images taken in response to BATSE triggers delivered via the GCN. For the first time, we have analyzed the entire 16 degree by 16 degree field covered for five of these bursts. In addition, we discuss a search for the optical counterpart to GRB 000201, a well-localized long burst. Single image sensitivities range from 13th to 14th magnitude around 10 s after the initial burst detection, and 14 - 15.8 one hour later. No new optical counterparts were discovered in this analysis suggesting short burst optical and gamma-ray fluxes are uncorrelated.

We present a lensing study of 42 galaxy clusters imaged in Sloan Digital Sky Survey (SDSS) commissioning data. Cluster candidates are selected optically from SDSS imaging data and confirmed for this study by matching to X-ray sources found independently in the ROSAT all sky survey (RASS). Five color SDSS photometry is used to make accurate photometric redshift estimates that are used to rescale and combine the lensing measurements. The mean shear from these clusters is detected to 2 h-1 Mpc at the 7-sigma level, corresponding to a mass within that radius of 4.2 +/- 0.6 x 10^14 h-1 M_sun. The shear profile is well fit by a power law with index -0.9 +/- 0.3, consistent with that of an isothermal density profile. This paper demonstrates our ability to measure ensemble cluster masses from SDSS imaging data.

We outline here the next generation of cluster-finding algorithms. We show how advances in Computer Science and Statistics have helped develop robust, fast algorithms for finding clusters of galaxies in large multi-dimensional astronomical databases like the Sloan Digital Sky Survey (SDSS). Specifically, this paper presents four new advances: (1) A new semi-parametric algorithm - nicknamed ``C4'' - for jointly finding clusters of galaxies in the SDSS and ROSAT All-Sky Survey databases; (2) The introduction of the False Discovery Rate into Astronomy; (3) The role of kernel shape in optimizing cluster detection; (4) A new determination of the X-ray Cluster Luminosity Function which has bearing on the existence of a ``deficit'' of high redshift, high luminosity clusters. This research is part of our ``Computational AstroStatistics'' collaboration (see Nichol et al. 2000) and the algorithms and techniques discussed herein will form part of the ``Virtual Observatory'' analysis toolkit.

We present two weak lensing reconstructions of the nearby ($z_{cl}=0.055$) merging cluster Abell 3667, based on observations taken $\sim 1$ year apart under different seeing conditions. This is the lowest redshift cluster with a weak lensing mass reconstruction to date. The reproducibility of features in the two mass maps demonstrate that weak lensing studies of low-redshift clusters are feasible. These data constitute the first results from an X-ray luminosity-selected weak lensing survey of 19 low-redshift ($z< 0.1$) southern clusters.

In order to observe nearly simultaneous emission from Gamma-ray Bursts (GRBs), the Robotic Optical Transient Search Experiment (ROTSE) receives triggers via the GRB Coordinates Network (GCN). Since beginning operations in March, 1998, ROTSE has also taken useful data for 10 SGR events: 8 from SGR 1900+14 and 2 from SGR 1806-20. We have searched for new or variable sources in the error regions of these SGRs and no optical counterparts were observed. Limits are in the range m_ROTSE ~ 12.5 - 15.5 during the period 20 seconds to 1 hour after the observed SGR events.

The Robotic Optical Transient Search Experiment (ROTSE) seeks to measure simultaneous and early afterglow optical emission from gamma-ray bursts (GRBs). A search for optical counterparts to six GRBs with localization errors of 1 square degree or better produced no detections. The earliest limiting sensitivity is m(ROTSE) > 13.1 at 10.85 seconds (5 second exposure) after the gamma-ray rise, and the best limit is m(ROTSE) > 16.0 at 62 minutes (897 second exposure). These are the most stringent limits obtained for GRB optical counterpart brightness in the first hour after the burst. Consideration of the gamma-ray fluence and peak flux for these bursts and for GRB990123 indicates that there is not a strong positive correlation between optical flux and gamma-ray emission.

The ROTSE-I experiment has generated CCD photometry for the entire Northern sky in two epochs nightly since March 1998. These sky patrol data are a powerful resource for studies of astrophysical transients. As a demonstration project, we present first results of a search for periodic variable stars derived from ROTSE-I observations. Variable identification, period determination, and type classification are conducted via automatic algorithms. In a set of nine ROTSE-I sky patrol fields covering about 2000 square degrees we identify 1781 periodic variable stars with mean magnitudes between m_v=10.0 and m_v=15.5. About 90% of these objects are newly identified as variable. Examples of many familiar types are presented. All classifications for this study have been manually confirmed. The selection criteria for this analysis have been conservatively defined, and are known to be biased against some variable classes. This preliminary study includes only 5.6% of the total ROTSE-I sky coverage, suggesting that the full ROTSE-I variable catalog will include more than 32,000 periodic variable stars.

The Robotic Optical Transient Search Experiment (ROTSE) seeks to measure contemporaneous and early afterglow optical emission from gamma-ray bursts (GRBs). The ROTSE-I telescope array has been fully automated and responding to burst alerts from the GRB Coordinates Network since March 1998, taking prompt optical data for 30 bursts in its first year. We will briefly review observations of GRB990123 which revealed the first detection of an optical burst occurring during the gamma-ray emission, reaching 9th magnitude at its peak. In addition, we present here preliminary optical results for seven other gamma-ray bursts. No other optical counterparts were seen in this analysis, and the best limiting sensitivities are m(V) > 13.0 at 14.7 seconds after the gamma-ray rise, and m(V) > 16.4 at 62 minutes. These are the most stringent limits obtained for GRB optical counterpart brightness in the first hour after the burst. This analysis suggests that there is not a strong correlation between optical flux and gamma-ray emission.

We describe first results of a project to create weak lensing mass maps for a complete, X-ray luminosity-limited sample of 19 nearby (z < 0.1) southern galaxy clusters scheduled for Sunyaev-Zel'dovich observations by the Viper Telescope at the South Pole. We have collected data on 1/3 of the sample and present motivation for the project as well as projected mass maps of two clusters.

The origin of gamma-ray bursts (GRBs) has been enigmatic since their discovery. The situation improved dramatically in 1997, when the rapid availability of precise coordinates for the bursts allowed the detection of faint optical and radio afterglows - optical spectra thus obtained have demonstrated conclusively that the bursts occur at cosmological distances. But, despite efforts by several groups, optical detection has not hitherto been achieved during the brief duration of a burst. Here we report the detection of bright optical emission from GRB990123 while the burst was still in progress. Our observations begin 22 seconds after the onset of the burst and show an increase in brightness by a factor of 14 during the first 25 seconds; the brightness then declines by a factor of 100, at which point (700 seconds after the burst onset) it falls below our detection threshold. The redshift of this burst, approximately 1.6, implies a peak optical luminosity of 5 times 10^{49} erg per second. Optical emission from gamma-ray bursts has been generally thought to take place at the shock fronts generated by interaction of the primary energy source with the surrounding medium, where the gamma-rays might also be produced. The lack of a significant change in the gamma-ray light curve when the optical emission develops suggests that the gamma-rays are not produced at the shock front, but closer to the site of the original explosion.

We present deep images of the field of gamma-ray burst (GRB) 990123 obtained in a broad-band UV/visible bandpass with the Hubble Space Telescope, and deep near-infrared images obtained with the Keck-I 10-m telescope. Both the HST and Keck images show that the optical transient (OT) is clearly offset by 0.6 arcsec from an extended object, presumably the host galaxy. This galaxy is the most likely source of the metallic-line absorption at z = 1.6004 seen in the spectrum of the OT. With magnitudes V_{C} ~ 24.6 +/- 0.2 and K = 21.65 +/- 0.30 mag this corresponds to an L ~ 0.7 L_* galaxy, assuming that it is located at z = 1.6. The estimated unobscured star formation rate is SFR ~ 6 M_sun/yr, which is not unusually high for normal galaxies at comparable redshifts. The strength of the observed metallic absorption lines is suggestive of a relatively high metallicity of the gas, and thus of a chemically evolved system which may be associated with a massive galaxy. It is also indicative of a high column density of the gas, typical of damped Ly-alpha systems at high redshifts. We conclude that this is the host galaxy of GRB 990123. No other obvious galaxies are detected within the same projected radius from the OT. There is thus no evidence for strong gravitational lensing magnification of this burst, and some alternative explanation for its remarkable energetics may be required. The observed offset of the OT from the center of its apparent host galaxy, 5.5 +/- 0.9 proper kpc (projected) in the galaxy's rest-frame, both refutes the possibility that GRBs are related to galactic nuclear activity and supports models of GRBs which involve the death and/or merger of massive stars. Further, the HST image suggests an intimate connection of GRB 990123 and a star-forming region.

A next generation atmospheric Cherenkov observatory is described based on the Whipple Observatory $\gamma$-ray telescope. A total of nine such imaging telescopes will be deployed in an array that will permit the maximum versatility and give high sensitivity in the 50 GeV - 50 TeV band (with maximum sensitivity from 100 GeV to 10 TeV).

While the strategy for the first applications of weak lensing has been to ``go deep'' it is equally interesting to use one's telescope time to instead ``go wide''. The Sloan Survey (SDSS) provides a natural framework for a very wide area weak lensing survey.