We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

Quasars have long been known as intrinsically variable sources, but the physical mechanism underlying the temporal optical/UV variability is still not well understood. We propose a novel nonparametric method for modeling and forecasting the optical variability of quasars utilizing an autoencoder neural network to gain insight into the underlying processes. The autoencoder is trained with ~15,000 decade-long quasar light curves obtained by the Catalina Real-time Transient Survey selected with negligible flux contamination from the host galaxy. The autoencoder's performance in forecasting the temporal flux variation of quasars is superior to that of the damped random walk process. We find a temporal asymmetry in the optical variability and a novel relation - the amplitude of the variability asymmetry decreases as luminosity and/or black hole mass increases - is suggested with the help of autoencoded features. The characteristics of the variability asymmetry are in agreement with those from the self-organized disk instability model, which predicts that the magnitude of the variability asymmetry decreases as the ratio of the diffusion mass to inflow mass in the accretion disk increases.

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

We present extensive datasets for a class of intermediate-luminosity optical transients known as "luminous red novae" (LRNe). They show double-peaked light curves, with an initial rapid luminosity rise to a blue peak (at -13 to -15 mag), which is followed by a longer-duration red peak that sometimes is attenuated, resembling a plateau. The progenitors of three of them (NGC4490-2011OT1, M101-2015OT1, and SNhunt248), likely relatively massive blue to yellow stars, were also observed in a pre-eruptive stage when their luminosity was slowly increasing. Early spectra obtained during the first peak show a blue continuum with superposed prominent narrow Balmer lines, with P Cygni profiles. Lines of Fe II are also clearly observed, mostly in emission. During the second peak, the spectral continuum becomes much redder, Halpha is barely detected, and a forest of narrow metal lines is observed in absorption. Very late-time spectra (~6 months after blue peak) show an extremely red spectral continuum, peaking in the infrared (IR) domain. Halpha is detected in pure emission at such late phases, along with broad absorption bands due to molecular overtones (such as TiO, VO). We discuss a few alternative scenarios for LRNe. Although major instabilities of single massive stars cannot be definitely ruled out, we favour a common envelope ejection in a close binary system, with possibly a final coalescence of the two stars. The similarity between LRNe and the outburst observed a few months before the explosion of the Type IIn SN 2011ht is also discussed.

Merger simulations predict that tidally induced gas inflows can trigger kpc-scale dual active galactic nuclei (dAGN) in heavily obscured environments. Previously with the Very Large Array, we have confirmed four dAGN with redshifts between $0.04 < z < 0.22$ and projected separations between 4.3 and 9.2 kpc in the SDSS Stripe 82 field. Here, we present $Chandra$ X-ray observations that spatially resolve these dAGN and compare their multi-wavelength properties to those of single AGN from the literature. We detect X-ray emission from six of the individual merger components and obtain upper limits for the remaining two. Combined with previous radio and optical observations, we find that our dAGN have properties similar to nearby low-luminosity AGN, and they agree well with the black hole fundamental plane relation. There are three AGN-dominated X-ray sources, whose X-ray hardness-ratio derived column densities show that two are unobscured and one is obscured. The low obscured fraction suggests these dAGN are no more obscured than single AGN, in contrast to the predictions from simulations. These three sources show an apparent X-ray deficit compared to their mid-infrared continuum and optical [OIII] line luminosities, suggesting higher levels of obscuration, in tension with the hardness-ratio derived column densities. Enhanced mid-infrared and [OIII] luminosities from star formation may explain this deficit. There is ambiguity in the level of obscuration for the remaining five components since their hardness ratios may be affected by non-nuclear X-ray emissions, or are undetected altogether. They require further observations to be fully characterized.

Intense, millisecond-duration bursts of radio waves have been detected from beyond the Milky Way [1]. Their extragalactic origins are evidenced by their large dispersion measures, which are greater than expected for propagation through the Milky Way interstellar medium alone, and imply contributions from the intergalactic medium and potentially host galaxies [2]. Although several theories exist for the sources of these fast radio bursts, their intensities, durations and temporal structures suggest coherent emission from highly magnetised plasma [3,4]. Two sources have been observed to repeat [5,6], and one repeater (FRB 121102) has been localised to the largest star-forming region of a dwarf galaxy at a cosmological redshift of 0.19 [7, 8]. However, the host galaxies and distances of the so far non-repeating fast radio bursts are yet to be identified. Unlike repeating sources, these events must be observed with an interferometer with sufficient spatial resolution for arcsecond localisation at the time of discovery. Here we report the localisation of a fast radio burst (FRB 190523) to a few-arcsecond region containing a single massive galaxy at a redshift of 0.66. This galaxy is in stark contrast to the host of FRB 121102, being a thousand times more massive, with a greater than hundred times lower specific star-formation rate. The properties of this galaxy highlight the possibility of a channel for FRB production associated with older stellar populations.

We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$\beta$/[OIII] ratio. We find 36 sources with declining H$\beta$ flux, 37 sources with increasing H$\beta$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$\beta$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

In this work we use structural properties of Milky Way's outer halo ($R_G > 25\,\mathrm{kpc}$) satellites (dwarf spheroidal galaxies, ultra-faint dwarf galaxies and globular clusters) derived from deep, wide-field and homogeneous data, to present evidence of a correlation in the S\'ersic index v/s effective radius plane followed by a large fraction of outer halo globular clusters and satellite dwarf galaxies. We show that this correlation can be entirely reproduced by fitting empirical relations in the central surface brightness v/s absolute magnitude and S\'ersic index v/s absolute magnitude parameter spaces, and by assuming the existence of two types of outer halo globular clusters: one of high surface brightness (HSB group), with properties similar to inner halo clusters; and another of low surface brightness (LSB group), which share characteristics with dwarf spheroidal and ultra-faint dwarf galaxies. Given the similarities of LSB clusters with dwarf spheroidal and ultra-faint dwarf galaxies, we discuss the possibility that outer halo clusters also originated inside dark matter halos and that tidal forces from different galaxy host's potentials are responsible for the different properties between HSB and LSB clusters.

Aims: In this work, we aim to provide a reliable list of gravitational lens (GL) candidates based on a search performed over the entire Gaia Data Release 2 (Gaia DR2). We also show that the sole astrometric and photometric informations coming from the Gaia satellite yield sufficient insights for supervised learning methods to automatically identify GL candidates with an efficiency that is comparable to methods based on image processing. Methods: We simulated 106,623,188 lens systems composed of more than two images, based on a regular grid of parameters characterizing a non-singular isothermal ellipsoid lens model in the presence of an external shear. These simulations are used as an input for training and testing our supervised learning models consisting of Extremely Randomized Trees. The latter are finally used to assign to each of the 2,129,659 clusters of celestial objects a discriminant value that reflects the ability of our simulations to match the observed relative positions and fluxes from each cluster. Once complemented with additional constraints, these discriminant values allowed us to identify GL candidates out of the list of clusters. Results: We report the discovery of 15 new quadruply-imaged lens candidates with angular separations less than 6" and assess the performance of our approach by recovering 12 out of the 13 known quadruply-imaged systems with all their components detected in Gaia DR2 with a misclassification rate of fortuitous clusters of stars as lens systems that is below one percent. Similarly, the identification capability of our method regarding quadruply-imaged systems where three images are detected in Gaia DR2 is assessed by recovering 10 out of the 13 known quadruply-imaged systems having one of their constituting images discarded. The associated misclassification rate varying then between 5.8% and 20%, depending on the image we decided to remove.

We present an analysis of nine years of Catalina Surveys optical photometry for 1.4 million spectroscopically confirmed SDSS galaxies. We find 717 outburst events that were not reported by ongoing transient surveys. These events have timescales ranging from weeks to years. More than two thirds of these new events are found in starforming galaxies, while such galaxies only constitute ~20% of our sample. Based on the properties of the hosts and events, we find that almost all of the new events are likely to be associated with regular supernovae. However, a small number of long-timescale events are found among the galaxies containing AGN. These events have similar properties to those recently found in the analyses of light curves of large samples of AGN. Given the lack of such events among the more than a million passive galaxies in the sample, we suggest that the long outbursts are associated with super-massive black holes or their environments.

Galaxies are believed to experience star formation and black hole driven nuclear activity symbiotically. The symbiosis may be more extreme in the distant universe, as far-infrared photometry with the Herschel Space Observatory has found many cases of ultra-luminous cool dust emission in z>1 radio galaxies and quasars, which could have its origin in the central black hole activity, or in extreme starbursts. We here present strong evidence for an extreme circumnuclear starburst in the z=1.439 quasar 3C298. Our unparalleled 0.18 arcsecond resolution ALMA image at rest-frame 410micrometer wavelength shows that the ~40K dust in its host galaxy resides in an asymmetric circumnuclear structure. The morphology of this structure implies a starburst origin and a symbiotic physical relation with the AGN driven radio source. The symbiosis is likely to be a general property of distant massive galaxies.

We present an updated catalog of 4680 northern eclipsing binaries (EBs) with Algol-type light curve morphology (i.e., with well-defined beginning and end of primary and secondary eclipses), using data from the Catalina Sky Surveys. Our work includes revised period determinations, phenomenological parameters of the light curves, and system morphology classification based on machine learning techniques. While most of the new periods are in excellent agreement with those provided in the original Catalina catalogs, improved values are now available for ~10% of the stars. A total of 3456 EBs were classified as detached and 449 as semi-detached, while 145 cannot be classified unambiguously into either subtype. The majority of the SD systems seems to be comprised of short-period Algols. By applying color criteria, we searched for K- and M-type dwarfs in these data, and present a subsample of 609 EB candidates for further investigation. We report 119 EBs (2.5% of the total sample) that show maximum quadrature light variations over long timescales, with periods bracketing the range 4.5-18 yrs and fractional luminosity variance of 0.04-0.13. We discuss possible causes for this, making use of models of variable starspot activity in our interpretation of the results.

We present the discovery of a gravitationally lensed dust-reddened QSO at $z=2.517$ discovered in a survey for red QSOs by infrared selection. $Hubble~Space~Telescope$ imaging in the WFC3/IR F160W and F125W filters reveals a quadruply lensed system in a cusp configuration. We find that compared to the central image of the cusp, the nearby, brightest image is anomalous by a factor of $\sim7-11$. Although the source is extremely bright in the mid-infrared, a magnification by a factor of $\sim50-120$ places an upper limit of 1.35 mJy on the intrinsic mid-infrared brightness, well below the $WISE~W4$ detection limit of 6 mJy. We find that this QSO is moderately reddened, with $E(B-V)=0.7$ and that $\sim1\%$ of the intrinsic spectrum is leaked back into the line of sight resulting in an upturn in its UV spectrum. We conclude that the QSO's reddening is intrinsic and not due to the lens. Consistent with previous red quasar samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. Depending on how $L_{\rm bol}$ is computed, the quasar's accretion rate may be as high as $0.26~L_{\rm Edd}$. We detect two Lyman limit systems, at $z=2.102$ and $z=2.431$, with absorption by metal lines likely at small impact parameter to the QSO, and a putative lens redshift of $z=0.599$. Given the rarity of quad lenses, the discovery of this source allows detailed study of a less luminous, more typical infrared-selected quasar at high redshift.

We describe a deep, systematic imaging study of satellites in the outer halo of the Milky Way. Our sample consists of 58 stellar overdensities --- i.e., substructures classified as either globular clusters, classical dwarf galaxies, or ultra-faint dwarf galaxies --- that are located at Galactocentric distances of R$_{\rm GC}$ > 25 kpc (outer halo) and out to ~400 kpc. This includes 44 objects for which we have acquired deep, wide-field, $g-$ and $r-$band imaging with the MegaCam mosaic cameras on the 3.6m Canada-France-Hawaii Telescope and the 6.5m Magellan-Clay telescope. These data are supplemented by archival imaging, or published $gr$ photometry, for an additional 14 objects, most of which were discovered recently in the Dark Energy Survey (DES). We describe the scientific motivation for our survey, including sample selection, observing strategy, data reduction pipeline, calibration procedures, and the depth and precision of the photometry. The typical 5$\sigma$ point-source limiting magnitudes for our MegaCam imaging --- which collectively covers an area of ~52 deg$^{2}$ --- are $g_{\rm lim}$ ~25.6 and $r_{\rm lim}$ ~25.3 AB mag. These limits are comparable to those from the coadded DES images and are roughly a half-magnitude deeper than will be reached in a single visit with LSST. Our photometric catalog thus provides the deepest and most uniform photometric database of Milky Way satellites available for the foreseeable future. In other papers in this series, we have used these data to explore the blue straggler populations in these objects, their density distributions, star formation histories, scaling relations and possible foreground structures.

We present structural parameters from a wide-field homogeneous imaging survey of Milky Way satellites carried out with the MegaCam imagers on the 3.6m Canada-France-Hawaii Telescope (CFHT) and 6.5m Magellan-Clay telescope. Our survey targets an unbiased sample of "outer halo" satellites (i.e., substructures having Galactocentric distances greater than 25 kpc) and includes classical dSph galaxies, ultra-faint dwarfs, and remote globular clusters. We combine deep, panoramic $gr$ imaging for 44 satellites and archival $gr$ imaging for 14 additional objects (primarily obtained with the DECam instrument as part of the Dark Energy Survey), to measure photometric and structural parameters for 58 outer halo satellites. This is the largest and most uniform analysis of Milky Way satellites undertaken to date and represents roughly three quarters ($58/81 \simeq$72\%) of all known outer halo satellites. We use a maximum-likelihood method to fit four density laws to each object in our survey: exponential, Plummer, King and Sersic models. We examine systematically the isodensity contour maps and color magnitude diagrams for each of our program objects, present a comparison with previous results, and tabulate our best-fit photometric and structural parameters, including ellipticities, position angles, effective radii, Sersic indices, absolute magnitudes, and surface brightness measurements. We investigate the distribution of outer halo satellites in the size-magnitude diagram, and show that the current sample of outer halo substructures spans a wide range in effective radius, luminosity and surface brightness, with little evidence for a clean separation into star cluster and galaxy populations at the faintest luminosities and surface brightnesses.

Thanks to its spatial resolution the ESA/Gaia space mission offers a unique opportunity to discover new multiply-imaged quasars and to study the already known lensed systems at sub-milliarcsecond astrometric precisions. In this paper, we address the detection of the known multiply-imaged quasars from the Gaia Data Release 2 and determine the astrometric and photometric properties of the individually detected images found in the Gaia DR2 catalogue. We have compiled an exhaustive list of quasar gravitational lenses from the literature to search for counterparts in the Gaia Data Release 2. We then analyze the astrometric and photometric properties of these Gaia's detections. To highlight the tremendous potential of Gaia at the sub-milliarcsecond level we finally perform a simple Bayesian modeling of the well-known gravitational lens system HE0435-1223, using Gaia Data Release 2 and HST astrometry. From 478 known multiply imaged quasars, 200 have at least one image found in the Gaia Data Release 2. Among the 41 known quadruply-imaged quasars of the list, 26 have at least one image in the Gaia Data Release 2, 12 of which are fully detected, 6 have three counterparts, 7 have two and 1 has only one. As expected, the modeling of HE0435-1223 shows that the model parameters are significantly better constrained when using Gaia astrometry compared to HST astrometry, in particular the relative positions of the background quasar source and the centroid of the deflector. The Gaia sub-milliarcsecond astrometry also significantly reduces the parameter correlations. Besides providing an up-to-date list of multiply imaged quasars and their detection in the Gaia DR2, this paper shows that more complex modeling scenarios will certainly benefit from Gaia sub-milliarcsecond astrometry.

We present a spectroscopically complete sample of 147 infrared-color-selected AGN down to a 22 $\mu$m flux limit of 20 mJy over the $\sim$270 deg$^2$ of the SDSS Stripe 82 region. Most of these sources are in the QSO luminosity regime ($L_{\rm bol} \gtrsim 10^{12} L_\odot$) and are found out to $z\simeq3$. We classify the AGN into three types, finding: 57 blue, unobscured Type-1 (broad-lined) sources; 69 obscured, Type-2 (narrow-lined) sources; and 21 moderately-reddened Type-1 sources (broad-lined and $E(B-V) > 0.25$). We study a subset of this sample in X-rays and analyze their obscuration to find that our spectroscopic classifications are in broad agreement with low, moderate, and large amounts of absorption for Type-1, red Type-1 and Type-2 AGN, respectively. We also investigate how their X-ray luminosities correlate with other known bolometric luminosity indicators such as [O III] line luminosity ($L_{\rm [OIII]}$) and infrared luminosity ($L_{6 \mu{\rm m}}$). While the X-ray correlation with $L_{\rm [OIII]}$ is consistent with previous findings, the most infrared-luminous sources appear to deviate from established relations such that they are either under-luminous in X-rays or over-luminous in the infrared. Finally, we examine the luminosity function (LF) evolution of our sample, and by AGN type, in combination with the complementary, infrared-selected, AGN sample of Lacy et al. (2013), spanning over two orders of magnitude in luminosity. We find that the two obscured populations evolve differently, with reddened Type-1 AGN dominating the obscured AGN fraction ($\sim$30%) for $L_{5 \mu{\rm m}} > 10^{45}$ erg s$^{-1}$, while the fraction of Type-2 AGN with $L_{5 \mu{\rm m}} < 10^{45}$ erg s$^{-1}$ rises sharply from 40% to 80% of the overall AGN population.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

Context. Strong gravitationally lensed quasars are among the most interesting and useful observable extragalactic phenomena. Because their study constitutes a unique tool in various fields of astronomy, they are highly sought, not without difficulty. Indeed, even in this era of all-sky surveys, their recognition remains a great challenge, with barely a few hundred currently known systems. Aims. In this work we aim to detect new strongly lensed quasar candidates in the recently published Gaia Data Release 2 (DR2), which is the highest spatial resolution astrometric and photometric all-sky survey, attaining effective resolutions from 0.4" to 2.2". Methods. We cross-matched a merged list of quasars and candidates with the Gaia DR2 and found 1,839,143 counterparts within 0.5". We then searched matches with more than two Gaia DR2 counterparts within 6". We further narrowed the resulting list using astrometry and photometry compatibility criteria between the Gaia DR2 counterparts. A supervised machine learning method, Extremely Randomized Trees, is finally adopted to assign to each remaining system a probability of being lensed. Results. We report the discovery of three quadruply-imaged quasar candidates that are fully detected in Gaia DR2. These are the most promising new quasar lens candidates from Gaia DR2 and a simple singular isothermal ellipsoid lens model is able to reproduce their image positions to within $\sim$1 mas. This letter demonstrates the gravitational lens discovery potential of Gaia.

We report the results of a systematic photometric survey of the peripheral regions of a sample of fourteen globular clusters in the outer halo of the Milky Way at distances d_GC>25 kpc from the Galactic centre. The survey is aimed at searching for the remnants of the host satellite galaxies where these clusters could originally have been formed before being accreted onto the Galactic halo. The limiting surface brightness varies within our sample, but reaches muV_lim=30-32 mag arcsec^-2. For only two globular clusters (NGC 7492 and Whiting 1; already suggested to be associated with the Sagittarius galaxy) we detect extended stellar populations that cannot be associated with either the clusters themselves or with the surrounding Galactic field population. We show that the lack of substructures around globular clusters at these Galactocentric distances is still compatible with the predictions of cosmological simulations whereby in the outer halo the Galactic globular cluster system is built up through hierarchical accretion at early epochs.

SNhunt151 was initially classified as a supernova (SN) impostor (nonterminal outburst of a massive star). It exhibited a slow increase in luminosity, lasting about 450 d, followed by a major brightening that reaches M_V ~ -18 mag. No source is detected to M_V > -13 mag in archival images at the position of SNhunt151 before the slow rise. Low-to-mid-resolution optical spectra obtained during the pronounced brightening show very little evolution, being dominated at all times by multicomponent Balmer emission lines, a signature of interaction between the material ejected in the new outburst and the pre-existing circumstellar medium. We also analyzed mid-infrared images from the Spitzer Space Telescope, detecting a source at the transient position in 2014 and 2015. Overall, SNhunt151 is spectroscopically a Type IIn SN, somewhat similar to SN2009ip. However, there are also some differences, such as a slow pre-discovery rise, a relatively broad light-curve peak showing a longer rise time (~ 50 d) and a slower decline, along with a negligible change in the temperature around the peak (T < 10^4 K). We suggest that SNhunt151 is the result of an outburst, or a SN explosion, within a dense circumstellar nebula, similar to those embedding some luminous blue variables like Eta Carinae and originating from past mass-loss events.

Astronomy light curves are sparse, gappy, and heteroscedastic. As a result standard time series methods regularly used for financial and similar datasets are of little help and astronomers are usually left to their own instruments and techniques to classify light curves. A common approach is to derive statistical features from the time series and to use machine learning methods, generally supervised, to separate objects into a few of the standard classes. In this work, we transform the time series to two-dimensional light curve representations in order to classify them using modern deep learning techniques. In particular, we show that convolutional neural networks based classifiers work well for broad characterization and classification. We use labeled datasets of periodic variables from CRTS survey and show how this opens doors for a quick classification of diverse classes with several possible exciting extensions.

We report the results of our radio, optical and infra-red studies of a peculiar radio source 4C~35.06, an extended radio-loud AGN at the center of galaxy cluster Abell 407 ($z=0.047$). The central region of this cluster hosts a remarkably tight ensemble of nine galaxies, the spectra of which resemble those of passive red ellipticals, embedded within a diffuse stellar halo of $\sim$1~arcmin size. This system (named the `Zwicky's Nonet') provides unique and compelling evidence for a multiple-nucleus cD galaxy precursor. Multifrequency radio observations of 4C~35.06 with the Giant Meterwave Radio Telescope (GMRT) at 610, 235 and 150 MHz reveal a system of 400~kpc scale helically twisted and kinked radio jets and outer diffuse lobes. The outer extremities of jets contain extremely steep spectrum (spectral index -1.7 to -2.5) relic/fossil radio plasma with a spectral age of a few$\,\times (10^7 - 10^8)$ yr. Such ultra-steep spectrum relic radio lobes without definitive hot-spots are rare, and they provide an opportunity to understand the life-cycle of relativistic jets and physics of black hole mergers in dense environments. We interpret our observations of this radio source in the context of the growth of its central black hole, triggering of its AGN activity and jet precession, all possibly caused by galaxy mergers in this dense galactic system. A slow conical precession of the jet axis due to gravitational perturbation between interacting black holes is invoked to explain the unusual jet morphology.

Supernova (SN) 2016bdu is an unusual transient resembling SN 2009ip. SN 2009ip-like events are characterized by a long-lasting phase of erratic variability which ends with two luminous outbursts a few weeks apart. The second outburst is significantly more luminous (about 3 mag) than the first. In the case of SN 2016bdu, the first outburst (Event A) reached an absolute magnitude M(r) ~ -15.3 mag, while the second one (Event B) occurred over one month later and reached M(r) ~ -18 mag. By inspecting archival data, a faint source at the position of SN 2016bdu is detectable several times in the past few years. We interpret these detections as signatures of a phase of erratic variability, similar to that experienced by SN 2009ip between 2008 and mid-2012, and resembling the currently observed variability of the luminous blue variable SN 2000ch in NGC 3432. Spectroscopic monitoring of SN 2016bdu during the second peak initially shows features typical of a SN IIn. One month after the Event B maximum, the spectra develop broad Balmer lines with P Cygni profiles and broad metal features. At these late phases, the spectra resemble those of a typical Type II SN. All members of this SN 2009ip-like group are remarkably similar to the Type IIn SN 2005gl. For this object, the claim of a terminal SN explosion is supported by the disappearance of the progenitor star. The similarity with SN 2005gl suggests that all members of this family may finally explode as genuine SNe, although the unequivocal detection of nucleosynthesised elements in their nebular spectra is still missing.

A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic variability, which is not frequently observed in the CV systems. These long-term variations are possibly caused by various mechanisms, such as the precession of the accretion disk, hierarchical triple star system, magnetic field change of the companion star, and other possible mechanisms. We discuss the possible mechanisms in this study. If the long-term period is less than several tens of days, the disk precession period scenario is favored. However, the hierarchical triple star system or the variations in magnetic field strengths are most likely the predominant mechanisms for longer periods.

There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $\Delta m = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes.

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc.

We describe here the parallels in astronomy and earth science datasets, their analyses, and the opportunities for methodology transfer from astroinformatics to geoinformatics. Using example of hydrology, we emphasize how meta-data and ontologies are crucial in such an undertaking. Using the infrastructure being designed for EarthCube - the Virtual Observatory for the earth sciences - we discuss essential steps for better transfer of tools and techniques in the future e.g. domain adaptation. Finally we point out that it is never a one-way process and there is enough for astroinformatics to learn from geoinformatics as well.

The time domain is the emerging forefront of astronomical research with new facilities and instruments providing unprecedented amounts of data on the temporal behavior of astrophysical populations. Dealing with the size and complexity of this requires new techniques and methodologies. Quasars are an ideal work set for developing and applying these: they vary in a detectable but not easily quantifiable manner whose physical origins are poorly understood. In this paper, we will review how quasars are identified by their variability and how these techniques can be improved, what physical insights into their variability can be gained from studying extreme examples of variability, and what approaches can be taken to increase the number of quasars known. These will demonstrate how astroinformatics is essential to discovering and understanding this important population.

We obtained precise line-of-sight radial velocities of 23 member stars of the remote halo globular cluster Palomar 4 (Pal 4) using the High Resolution Echelle Spectrograph (HIRES) at the Keck I telescope. We also measured the mass function of the cluster down to a limiting magnitude of V~28 mag using archival HST/WFPC2 imaging. We derived the cluster's surface brightness profile based on the WFPC2 data and on broad-band imaging with the Low-Resolution Imaging Spectrometer (LRIS) at the Keck II telescope. We find a mean cluster velocity of 72.55+/-0.22 km/s and a velocity dispersion of 0.87+/-0.18 km/s. The global mass function of the cluster, in the mass range 0.55<=M<=0.85 M_solar, is shallower than a Kroupa mass function and the cluster is significantly depleted in low-mass stars in its center compared to its outskirts. Since the relaxation time of Pal 4 is of the order of a Hubble time, this points to primordial mass segregation in this cluster. Extrapolating the measured mass function towards lower-mass stars and including the contribution of compact remnants, we derive a total cluster mass of 29800 M_solar. For this mass, the measured velocity dispersion is consistent with the expectations of Newtonian dynamics and below the prediction of Modified Newtonian Dynamics (MOND). Pal 4 adds to the growing body of evidence that the dynamics of star clusters in the outer Galactic halo can hardly be explained by MOND.

We present the spatially resolved star formation history (SFH) of the Carina dwarf spheroidal galaxy, obtained from deep, wide-field g,r imaging and a metallicity distribution from the literature. Our photometry covers $\sim2$ deg$^2$, reaching up to $\sim10$ times the half-light radius of Carina with a completeness higher than $50\%$ at $g\sim24.5$, more than one magnitude fainter than the oldest turnoff. This is the first time a combination of depth and coverage of this quality has been used to derive the SFH of Carina, enabling us to trace its different populations with unprecedented accuracy. We find that Carina's SFH consists of two episodes well separated by a star formation temporal gap. These episodes occurred at old ($>10$ Gyr) and intermediate ($2$-$8$ Gyr) ages. Our measurements show that the old episode comprises the majority of the population, accounting for $54\pm5\%$ of the stellar mass within $1.3$ times the King tidal radius, while the total stellar mass derived for Carina is $1.60\pm0.09\times 10^{6} M_{\rm{\odot}}$, and the stellar mass-to-light ratio $1.8\pm0.2$. The SFH derived is consistent with no recent star formation which hints that the observed blue plume is due to blue stragglers. We conclude that the SFH of Carina evolved independently of the tidal field of the Milky Way, since the frequency and duration of its star formation events do not correlate with its orbital parameters. This result is supported by the age/metallicity relation observed in Carina, and the gradients calculated indicating that outer regions are older and more metal poor.

We report photometric and spectroscopic observations of the optical transient LSQ13zm. Historical data reveal the presence of an eruptive episode (that we label as `2013a') followed by a much brighter outburst (`2013b') three weeks later, that we argue to be the genuine supernova explosion. This sequence of events closely resemble those observed for SN2010mc and (in 2012) SN2009ip. The absolute magnitude reached by LSQ13zm during 2013a ($M_R=-14.87\pm0.25\,\rm{mag}$) is comparable with those of supernova impostors, while that of the 2013b event ($M_R=-18.46\pm0.21\,\rm{mag}$) is consistent with those of interacting supernovae. Our spectra reveal the presence of a dense and structured circumstellar medium, probably produced through numerous pre-supernova mass-loss events. In addition, we find evidence for high-velocity ejecta, with a fraction of gas expelled at more than 20000\kms. The spectra of LSQ13zm show remarkable similarity with those of well-studied core-collapse supernovae. From the analysis of the available photometric and spectroscopic data, we conclude that we first observed the last event of an eruptive sequence from a massive star, likely a Luminous Blue Variable, which a short time later exploded as a core-collapse supernova. The detailed analysis of archival images suggest that the host galaxy is a star-forming Blue Dwarf Compact Galaxy.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

The exploitation of present and future synoptic (multi-band and multi-epoch) surveys requires an extensive use of automatic methods for data processing and data interpretation. In this work, using data extracted from the Catalina Real Time Transient Survey (CRTS), we investigate the classification performance of some well tested methods: Random Forest, MLPQNA (Multi Layer Perceptron with Quasi Newton Algorithm) and K-Nearest Neighbors, paying special attention to the feature selection phase. In order to do so, several classification experiments were performed. Namely: identification of cataclysmic variables, separation between galactic and extra-galactic objects and identification of supernovae.

Big Data are revolutionizing nearly every aspect of the modern society. One area where this can have a profound positive societal impact is the field of Health Care Informatics (HCI), which faces many challenges. The key idea behind this study is: can we use some of the experience and technical and methodological solutions from the fields that have successfully adapted to the Big Data era, namely astronomy and space science, to help accelerate the progress of HCI? We illustrate this with examples from the Virtual Observatory framework, and the NCI EDRN project. An effective sharing and reuse of tools, methods, and experiences from different fields can save a lot of effort, time, and expense. HCI can thus benefit from the proven solutions to big data challenges from other domains.

Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths.

Representing simultaneous black hole accretion during a merger, binary active galactic nuclei (AGNs) could provide valuable observational constraints to models of galaxy mergers and AGN triggering. High-resolution radio interferometer imaging offers a promising method to identify a large and uniform sample of binary AGNs, because it probes a generic feature of nuclear activity and is free from dust obscuration. Our previous search yielded 52 strong candidates of kpc-scale binaries over the 92 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 area with 2"-resolution Very Large Array (VLA) images. Here we present 0.3"-resolution VLA 6 GHz observations for six candidates that have complete optical spectroscopy. The new data confirm the binary nature of four candidates and identify the other two as line-of-sight projections of radio structures from single AGNs. The four binary AGNs at $z \sim 0.1$ reside in major mergers with projected separations of 4.2-12 kpc. Optical spectral modeling shows that their hosts have stellar masses between $10.3 < log(M_\star/M_\odot) < 11.5$ and velocity dispersions between $120 < \sigma_\star < 320$ km/s. The radio emission is compact ($<$0.4") and show steep spectrum ($-1.8 < \alpha < -0.5$) at 6 GHz. The host galaxy properties and the Eddington-scaled accretion rates broadly correlate with the excitation state, similar to the general radio-AGN population at low redshifts. Our estimated binary AGN fraction indicates that simultaneous accretion occurs $>23^{+15}_{-8}$% of the time when a kpc-scale galaxy pair is detectable as a radio-AGN. The high duty cycle of the binary phase strongly suggests that major mergers can trigger and synchronize black hole accretion.

The optical light curve of the quasar PG 1302-102 at $z = 0.278$ shows a strong, smooth 5.2 yr periodic signal, detectable over a period of $\sim 20$ yr. Although the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim 10^{5}$ yr due to gravitational wave emission alone. Here we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from {\it WISE} and {\it Akari}, we confirm that the periodic behavior reported in the optical light curve from Graham et al. (2015) is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6 \mu$m time lags of $(2219 \pm 153, 2408 \pm 148)$ days for a near face-on orientation of the torus, or $(4103\pm 153, 4292 \pm 148)$ days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the further out dusty regions. The implied distance of this dusty, assumed toroidal region is $\sim$ 1.5 pc for a near face-on geometry, or $\sim$1.1 pc for the relativistic Doppler boosted case.

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

Hierarchical assembly models predict a population of supermassive black hole (SMBH) binaries. These are not resolvable by direct imaging but may be detectable via periodic variability (or nanohertz frequency gravitational waves). Following our detection of a 5.2 year periodic signal in the quasar PG 1302-102 (Graham et al. 2015), we present a novel analysis of the optical variability of 243,500 known spectroscopically confirmed quasars using data from the Catalina Real-time Transient Survey (CRTS) to look for close (< 0.1 pc) SMBH systems. Looking for a strong Keplerian periodic signal with at least 1.5 cycles over a baseline of nine years, we find a sample of 111 candidate objects. This is in conservative agreement with theoretical predictions from models of binary SMBH populations. Simulated data sets, assuming stochastic variability, also produce no equivalent candidates implying a low likelihood of spurious detections. The periodicity seen is likely attributable to either jet precession, warped accretion disks or periodic accretion associated with a close SMBH binary system. We also consider how other SMBH binary candidates in the literature appear in CRTS data and show that none of these are equivalent to the identified objects. Finally, the distribution of objects found is consistent with that expected from a gravitational wave-driven population. This implies that circumbinary gas is present at small orbital radii and is being perturbed by the black holes. None of the sources is expected to merge within at least the next century. This study opens a new unique window to study a population of close SMBH binaries that must exist according to our current understanding of galaxy and SMBH evolution.

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model. The recent availability of large collections of astronomical time series of flux measurements (light curves) offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 $\pm$ 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of $\sim9$ years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.

We use seven year's worth of observations from the Catalina Sky Survey and the Siding Spring Survey covering most of the northern and southern hemisphere at galactic latitudes higher than 20 degrees to search for serendipitously imaged moving objects in the outer solar system. These slowly moving objects would appear as stationary transients in these fast cadence asteroids surveys, so we develop methods to discover objects in the outer solar system using individual observations spaced by months, rather than spaced by hours, as is typically done. While we independently discover 8 known bright objects in the outer solar system, the faintest having $V=19.8\pm0.1$, no new objects are discovered. We find that the survey is nearly 100% efficient at detecting objects beyond 25 AU for $V\lesssim 19.1$ ($V\lesssim18.6$ in the southern hemisphere) and that the probability that there is one or more remaining outer solar system object of this brightness left to be discovered in the unsurveyed regions of the galactic plane is approximately 32%.

Galaxy mergers play an important role in the growth of galaxies and their supermassive black holes. Simulations suggest that tidal interactions could enhance black hole accretion, which can be tested by the fraction of binary active galactic nuclei (AGNs) among galaxy mergers. But determining the fraction requires a statistical sample of binaries. We have identified kpc-scale binary AGNs directly from high-resolution radio imaging. Inside the 92 square deg covered by the high-resolution Very Large Array survey of the Sloan Digital Sky Survey (SDSS) Stripe 82 field, we identified 22 grade A and 30 grade B candidates of binary radio AGNs with angular separations less than 5" (10 kpc at z = 0.1). Eight of the candidates have optical spectra for both components from the SDSS spectroscopic surveys and our Keck program. Two grade B candidates are projected pairs, but the remaining six candidates are all compelling cases of binary AGNs based on either emission line ratios or the excess in radio power compared to the H-alpha-traced star formation rate. Only two of the six binaries were previously discovered by an optical spectroscopic search. Based on these results, we estimate that ~60% of our binary candidates would be confirmed once we obtain complete spectroscopic information. We conclude that wide-area high-resolution radio surveys offer an efficient method to identify large samples of binary AGNs. These radio-selected binary AGNs complement binaries identified at other wavelengths and are useful for understanding the triggering mechanisms of black hole accretion.

We review some aspects of the current state of data-intensive astronomy, its methods, and some outstanding data analysis challenges. Astronomy is at the forefront of "big data" science, with exponentially growing data volumes and data rates, and an ever-increasing complexity, now entering the Petascale regime. Telescopes and observatories from both ground and space, covering a full range of wavelengths, feed the data via processing pipelines into dedicated archives, where they can be accessed for scientific analysis. Most of the large archives are connected through the Virtual Observatory framework, that provides interoperability standards and services, and effectively constitutes a global data grid of astronomy. Making discoveries in this overabundance of data requires applications of novel, machine learning tools. We describe some of the recent examples of such applications.

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leq\delta\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 \sigma$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40\deg$.

Effective data visualization is a key part of the discovery process in the era of big data. It is the bridge between the quantitative content of the data and human intuition, and thus an essential component of the scientific path from data into knowledge and understanding. Visualization is also essential in the data mining process, directing the choice of the applicable algorithms, and in helping to identify and remove bad data from the analysis. However, a high complexity or a high dimensionality of modern data sets represents a critical obstacle. How do we visualize interesting structures and patterns that may exist in hyper-dimensional data spaces? A better understanding of how we can perceive and interact with multi dimensional information poses some deep questions in the field of cognition technology and human computer interaction. To this effect, we are exploring the use of immersive virtual reality platforms for scientific data visualization, both as software and inexpensive commodity hardware. These potentially powerful and innovative tools for multi dimensional data visualization can also provide an easy and natural path to a collaborative data visualization and exploration, where scientists can interact with their data and their colleagues in the same visual space. Immersion provides benefits beyond the traditional desktop visualization tools: it leads to a demonstrably better perception of a datascape geometry, more intuitive data understanding, and a better retention of the perceived relationships in the data.

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

Over six years of operation, the Catalina Real-time Transient Survey (CRTS) has identified 1043 cataclysmic variable (CV) candidates --- the largest sample of CVs from a single survey to date. Here we provide spectroscopic identification of 85 systems fainter than g<19, including three AMCVn binaries, one helium-enriched CV, one polar and one new eclipsing CV. We analyse the outburst properties of the full sample and show that it contains a large fraction of low accretion rate CVs with long outburst recurrence times. We argue that most of the high accretion rate dwarf novae in the survey footprint have already been found and that future CRTS discoveries will be mostly low accretion rate systems. We find that CVs with white dwarf dominated spectra have significantly fewer outbursts in their CRTS light curves compared to disc-dominated CVs, reflecting the difference in their accretion rates. Comparing the CRTS sample to other samples of CVs, we estimate the overall external completeness to be 23.6 per cent, but show that as much as 56 per cent of CVs have variability amplitudes that are too small to be selected using the transient selection criteria employed by current ground-based surveys. The full table of CRTS CVs, including their outburst and spectroscopic properties examined in this paper, is provided in the online materials.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We investigate the properties of 367 ultra-short period binary candidates selected from 31,000 sources recently identified from Catalina Surveys data. Based on light curve morphology, along with WISE, SDSS and GALEX multi-colour photometry, we identify two distinct groups of binaries with periods below the 0.22 day contact binary minimum. In contrast to most recent work, we spectroscopically confirm the existence of M-dwarf+M-dwarf contact binary systems. By measuring the radial velocity variations for five of the shortest-period systems, we find examples of rare cool-white dwarf+M-dwarf binaries. Only a few such systems are currently known. Unlike warmer white dwarf systems, their UV flux and their optical colours and spectra are dominated by the M-dwarf companion. We contrast our discoveries with previous photometrically-selected ultra-short period contact binary candidates, and highlight the ongoing need for confirmation using spectra and associated radial velocity measurements. Overall, our analysis increases the number of ultra-short period contact binary candidates by more than an order of magnitude.

We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.

We present 855 cataclysmic variable candidates detected by the Catalina Real-time Transient Survey (CRTS) of which at least 137 have been spectroscopically confirmed and 705 are new discoveries. The sources were identified from the analysis of five years of data, and come from an area covering three quarters of the sky. We study the amplitude distribution of the dwarf novae CVs discovered by CRTS during outburst, and find that in quiescence they are typically two magnitudes fainter compared to the spectroscopic CV sample identified by SDSS. However, almost all CRTS CVs in the SDSS footprint have ugriz photometry. We analyse the spatial distribution of the CVs and find evidence that many of the systems lie at scale heights beyond those expected for a Galactic thin disc population. We compare the outburst rates of newly discovered CRTS CVs with the previously known CV population, and find no evidence for a difference between them. However, we find that significant evidence for a systematic difference in orbital period distribution. We discuss the CVs found below the orbital period minimum and argue that many more are yet to be identified among the full CRTS CV sample. We cross-match the CVs with archival X-ray catalogs and find that most of the systems are dwarf novae rather than magnetic CVs.

We present optical photometry and spectra of the super luminous type II/IIn supernova CSS121015:004244+132827 (z=0.2868) spanning epochs from -30 days (rest frame) to more than 200 days after maximum. CSS121015 is one of the more luminous supernova ever found and one of the best observed. The photometric evolution is characterized by a relatively fast rise to maximum (~40 days in the SN rest frame), and by a linear post-maximum decline. The light curve shows no sign of a break to an exponential tail. A broad Halpha is first detected at ~ +40 days (rest-frame). Narrow, barely-resolved Balmer and [O III] 5007 A lines, with decreasing strength, are visible along the entire spectral evolution. The spectra are very similar to other super luminous supernovae (SLSNe) with hydrogen in their spectrum, and also to SN 2005gj, sometimes considered a type Ia interacting with H-rich CSM. The spectra are also similar to a subsample of H-deficient SLSNe. We propose that the properties of CSS121015 are consistent with the interaction of the ejecta with a massive, extended, opaque shell, lost by the progenitor decades before the final explosion, although a magnetar powered model cannot be excluded. Based on the similarity of CSS121015 with other SLSNe (with and without H), we suggest that the shocked-shell scenario should be seriously considered as a plausible model for both types of SLSN.

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

We compare quasar selection techniques based on their optical variability using data from the Catalina Real-time Transient Survey (CRTS). We introduce a new technique based on Slepian wavelet variance (SWV) that shows comparable or better performance to structure functions and damped random walk models but with fewer assumptions. Combining these methods with WISE mid-IR colors produces a highly efficient quasar selection technique which we have validated spectroscopically. The SWV technique also identifies characteristic timescales in a time series and we find a characteristic rest frame timescale of ~54 days, confirmed in the light curves of ~18000 quasars from CRTS, SDSS and MACHO data, and anticorrelated with absolute magnitude. This indicates a transition between a damped random walk and $P(f) \propto f^{-1/3}$ behaviours and is the first strong indication that a damped random walk model may be too simplistic to describe optical quasar variability.

We use, for the first time, the Detrend Fluctuation Analysis (DFA) to study the correlation properties of the transmitted flux fluctuations, in the Lyman-$\alpha$ (Ly$\alpha$) Forest along the lines of sight (LOS) to QSOs, at different space scales. We consider in our analysis the transmitted flux in the intergalactic medium over the redshift range 2 $\leq$ z $\leq$ 4.5 from a sample of 45 high-quality medium resolution (R $\sim$ 4300) quasar spectra obtained with Echelle Spectrograph and Imager (ESI) mounted on the Keck II 10-m telescope, and from a sample of 19 high-quality high resolution (R $\sim$ 50000) quasar spectra obtained with Ultra-Violet and Visible Echelle Spectrograh (UVES) mounted on the ESO KUEYEN 8.2 m telescope. The result of the DFA method applied to both datasets, shows that there exists a difference in the correlation properties between the short and long-range regimes: the slopes of the transmitted flux fluctuation function are different on small and large scales. The scaling exponents, $\alpha_{1}$ = 1.635$\pm$ 0.115 and $\alpha_{2}$ = 0.758$\pm$ 0.085 for the ESI/Keck sample and $\alpha_{1}$ = 1.763 $\pm 0.128$ and $\alpha_{2}$ = 0.798 $\pm 0.084$ for the UVES/VLT sample for the short and long range regime respectively. The transition between the two regims is observed at about $\sim 1.4 h^{-1} Mpc$ (comoving). The fact that $\alpha_{1}$ is always larger than $\alpha_{2}$ for each spectrum supports the common view that the Universe is smoother on large scales than on small scales. The non detection of considerable variations in the scaling exponents from LOS to LOS confirms that anisotropies cannot be ubiquitous, at least on these scales.

The amount of collected data in many scientific fields is increasing, all of them requiring a common task: extract knowledge from massive, multi parametric data sets, as rapidly and efficiently possible. This is especially true in astronomy where synoptic sky surveys are enabling new research frontiers in the time domain astronomy and posing several new object classification challenges in multi dimensional spaces; given the high number of parameters available for each object, feature selection is quickly becoming a crucial task in analyzing astronomical data sets. Using data sets extracted from the ongoing Catalina Real-Time Transient Surveys (CRTS) and the Kepler Mission we illustrate a variety of feature selection strategies used to identify the subsets that give the most information and the results achieved applying these techniques to three major astronomical problems.

We present the results of a pilot survey to find dust-reddened quasars by matching the FIRST radio catalog to the UKIDSS near-infrared survey, and using optical data from SDSS to select objects with very red colors. The deep K-band limit provided by UKIDSS allows for finding more heavily-reddened quasars at higher redshifts as compared with previous work using FIRST and 2MASS. We selected 87 candidates with K<=17.0 from the UKIDSS Large Area Survey (LAS) First Data Release (DR1) which covers 190 deg2. These candidates reach up to ~1.5 magnitudes below the 2MASS limit and obey the color criteria developed to identify dust-reddened quasars. We have obtained 61 spectroscopic observations in the optical and/or near-infrared as well as classifications in the literature and have identified 14 reddened quasars with E(B-V)>0.1, including three at z>2. We study the infrared properties of the sample using photometry from the WISE Observatory and find that infrared colors improve the efficiency of red quasar selection, removing many contaminants in an infrared-to-optical color-selected sample alone. The highest-redshift quasars (z > 2) are only moderately reddened, with E(B-V) ~ 0.2-0.3. We find that the surface density of red quasars rises sharply with faintness, comprising up to 17% of blue quasars at the same apparent K-band flux limit. We estimate that to reach more heavily reddened quasars (i.e., E(B-V) > 0.5) at z>2 and a depth of K=17 we would need to survey at least ~2.5 times more area.

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

The ominous warnings of a `data deluge' in the life sciences from high-throughput DNA sequencing data are being supplanted by a second deluge, of cliches bemoaning our collective scientific fate unless we address the genomic data `tsunami'. It is imperative that we explore the many facets of the genome, not just sequence but also transcriptional and epigenetic variability, integrating these observations in order to attain a genuine understanding of how genes function, towards a goal of genomics-based personalized medicine. Determining any individual's genomic properties requires comparison to many others, sifting out the specific from the trends, requiring access to the many in order to yield information relevant to the few. This is the central big data challenge in genomics that still requires some sort of resolution. Is there a practical, feasible way of directly connecting the scientific community to this data universe? The best answer could be in the stars overhead.

We present the discovery of a plausible super-luminous supernova (SLSN), found in the archival data of Sloan Digital Sky Survey (SDSS) Stripe 82, called PSN 000123+000504. The supernova peaked at M_g<-21.3 mag in the second half of September 2005, but was missed by the real-time supernova hunt. The observed part of the light curve (17 epochs) showed that the rise to the maximum took over 30 days, while the decline time lasted at least 70 days (observed frame), closely resembling other SLSNe of SN2007bi type. Spectrum of the host galaxy reveals a redshift of z=0.281 and the distance modulus of \mu=40.77 mag. Combining this information with the SDSS photometry, we found the host galaxy to be an LMC-like irregular dwarf galaxy with the absolute magnitude of M_B=-18.2+/-0.2 mag and the oxygen abundance of 12+log[O/H]=8.3+/-0.2. Our SLSN follows the relation for the most energetic/super-luminous SNe exploding in low-metallicity environments, but we found no clear evidence for SLSNe to explode in low-luminosity (dwarf) galaxies only. The available information on the PSN 000123+000504 light curve suggests the magnetar-powered model as a likely scenario of this event. This SLSN is a new addition to a quickly growing family of super-luminous SNe.

We present Keck/DEIMOS spectroscopy and CFHT/MegaCam photometry for the Milky Way globular cluster Palomar 13. We triple the number of spectroscopically confirmed members, including many repeat velocity measurements. Palomar 13 is the only known globular cluster with possible evidence for dark matter, based on a Keck/HIRES 21 star velocity dispersion of sigma=2.2+/-0.4 km/s. We reproduce this measurement, but demonstrate that it is inflated by unresolved binary stars. For our sample of 61 stars, the velocity dispersion is sigma=0.7(+0.6/-0.5) km/s. Combining our DEIMOS data with literature values, our final velocity dispersion is sigma=0.4(+0.4/-0.3) km/s. We determine a spectroscopic metallicity of [Fe/H]=-1.6+/-0.1 dex, placing a 1-sigma upper limit of sigma_[Fe/H]~0.2 dex on any internal metallicity spread. We determine Palomar 13's total luminosity to be M_V=-2.8+/-0.4, making it among the least luminous known globular clusters. The photometric isophotes are regular out to the half-light radius and mildly irregular outside this radius. The outer surface brightness profile slope is shallower than typical globular clusters (eta=-2.8+/-0.3). Thus at large radius, tidal debris is likely affecting the appearance of Palomar 13. Combining our luminosity with the intrinsic velocity dispersion, we find a dynamical mass of of M_1/2=1.3(+2.7/-1.3)x10^3 M_sun and a mass-to-light ratio of M/L_V=2.4(+5.0/-2.4) M_sun/L_sun. Within our measurement errors, the mass-to-light ratio agrees with the theoretical predictions for a single stellar population. We conclude that, while there is some evidence for tidal stripping at large radius, the dynamical mass of Palomar 13 is consistent with its stellar mass and neither significant dark matter, nor extreme tidal heating, is required to explain the cluster dynamics.

This paper presents a comparison of popular period finding algorithms applied to the light curves of variable stars from the Catalina Real-time Transient Survey (CRTS), MACHO and ASAS data sets. We analyze the accuracy of the methods against magnitude, sampling rates, quoted period, quality measures (signal-to-noise and number of observations), variability, and object classes. We find that measure of dispersion-based techniques - analysis-of-variance with harmonics and conditional entropy - consistently give the best results but there are clear dependencies on object class and light curve quality. Period aliasing and identifying a period harmonic also remain significant issues. We consider the performance of the algorithms and show that a new conditional entropy-based algorithm is the most optimal in terms of completeness and speed. We also consider a simple ensemble approach and find that it performs no better than individual algorithms.

We present a homogeneous study of blue straggler stars across ten outer halo globular clusters, three classical dwarf spheroidal and nine ultra-faint galaxies based on deep and wide-field photometric data taken with MegaCam on the Canada-France-Hawaii Telescope. We find blue straggler stars to be ubiquitous among these Milky Way satellites. Based on these data, we can test the importance of primordial binaries or multiple systems on blue straggler star formation in low density environments. For the outer halo globular clusters we find an anti-correlation between the specific frequency of blue straggler and absolute magnitude, similar to that previously observed for inner halo clusters. When plotted against density and encounter rate, the frequency of blue stragglers are well fitted by single trends with smooth transitions between dwarf galaxies and globular clusters, which points to a common origin for their blue stragglers. The fraction of blue stragglers stays constant and high in the low encounter rate regime spanned by our dwarf galaxies, and decreases with density and encounter rate in the range spanned by our globular clusters. We find that young stars can mimic blue stragglers in dwarf galaxies only if their ages are 2.5+/-0.5 Gyr and they represent ~1-7% of the total number of stars, which we deem highly unlikely. These results point to mass-transfer or mergers of primordial binaries or multiple systems as the dominant blue straggler formation mechanism in low density systems.

This paper presents a new period finding method based on conditional entropy that is both efficient and accurate. We demonstrate its applicability on simulated and real data. We find that it has comparable performance to other information-based techniques with simulated data but is superior with real data, both for finding periods and just identifying periodic behaviour. In particular, it is robust against common aliasing issues found with other period-finding algorithms.

High-volume feature-rich data sets are becoming the bread-and-butter of 21st century astronomy but present significant challenges to scientific discovery. In particular, identifying scientifically significant relationships between sets of parameters is non-trivial. Similar problems in biological and geosciences have led to the development of systems which can explore large parameter spaces and identify potentially interesting sets of associations. In this paper, we describe the application of automated discovery systems of relationships to astronomical data sets, focussing on an evolutionary programming technique and an information-theory technique. We demonstrate their use with classical astronomical relationships - the Hertzsprung-Russell diagram and the fundamental plane of elliptical galaxies. We also show how they work with the issue of binary classification which is relevant to the next generation of large synoptic sky surveys, such as LSST. We find that comparable results to more familiar techniques, such as decision trees, are achievable. Finally, we consider the reality of the relationships discovered and how this can be used for feature selection and extraction.

We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted.

We present the analysis of 1,207 RR Lyrae found in photometry taken by the Catalina Survey's Mount Lemmon telescope. By combining accurate distances for these stars with measurements for ~14,000 type-AB RR Lyrae from the Catalina Schmid telescope, we reveal an extended association that reaches Galactocentric distances beyond 100 kpc and overlaps the Sagittarius streams system. This result confirms earlier evidence for the existence of an outer halo tidal stream resulting from a disrupted stellar system. By comparing the RR Lyrae source density with that expected based on halo models, we find the detection has ~8 sigma significance. We investigate the distances, radial velocities, metallicities, and period-amplitude distribution of the RR Lyrae. We find that both radial velocities and distances are inconsistent with current models of the Sagittarius stream. We also find tentative evidence for a division in source metallicities for the most distant sources. Following prior analyses, we compare the locations and distances of the RR Lyrae with photometrically selected candidate horizontal branch stars and find supporting evidence that this structure spans at least 60 deg of the sky. We investigate the prospects of an association between the stream and unusual globular cluster NGC 2419.

There are many occasions when one does not have complete information in order to classify objects into different classes, and yet it is important to do the best one can since other decisions depend on that. In astronomy, especially time-domain astronomy, this situation is common when a transient is detected and one wishes to determine what it is in order to decide if one must follow it. We propose to use the Difference Boosting Neural Network (DBNN) which can boost differences between feature vectors of different objects in order to differentiate between them. We apply it to the publicly available data of the Catalina Real-Time Transient Survey (CRTS) and present preliminary results. We also describe another use with a stellar spectral library to identify spectra based on a few features. The technique itself is more general and can be applied to a varied class of problems.

We present the analysis of 12227 type-ab RR Lyrae found among the 200 million public lightcurves in the Catalina Surveys Data Release 1 (CSDR1). These stars span the largest volume of the Milky Way ever surveyed with RR Lyrae, covering ~20,000 square degrees of the sky (0 < RA < 360, -22 < Dec < 65 deg) to heliocentric distances of up to 60kpc. Each of the RR Lyrae are observed between 60 and 419 times over a six-year period. Using period finding and Fourier fitting techniques we determine periods and apparent magnitudes for each source. We find that the periods at generally accurate to sigma = 0.002% by comparison with 2842 previously known RR Lyrae and 100 RR Lyrae observed in overlapping survey fields. We photometrically calibrate the light curves using 445 Landolt standard stars and show that the resulting magnitudes are accurate to ~0.05 mags using SDSS data for ~1000 blue horizontal branch stars and 7788 of the RR Lyrae. By combining Catalina photometry with SDSS spectroscopy, we analyze the radial velocity and metallicity distributions for > 1500 of the RR Lyrae. Using the accurate distances derived for the RR Lyrae, we show the paths of the Sagittarius tidal streams crossing the sky at heliocentric distances from 20 to 60 kpc. By selecting samples of Galactic halo RR Lyrae, we compare their velocity, metallicity, and distance with predictions from a recent detailed N-body model of the Sagittarius system. We find that there are some significant differences between the distances and structures predicted and our observations.

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This presents some unusual challenges: the data are sparse, heterogeneous and incomplete; evolving in time; and most of the relevant information comes not from the data stream itself, but from a variety of archival data and contextual information (spatial, temporal, and multi-wavelength). We are exploring a variety of novel techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. The current surveys are already overwhelming our ability to effectively follow all of the potentially interesting events, and these challenges will grow by orders of magnitude over the next decade as the more ambitious sky surveys get under way. While we focus on an application in a specific domain (astrophysics), these challenges are more broadly relevant for event or anomaly detection and knowledge discovery in massive data streams.

Astronomy has been at the forefront of the development of the techniques and methodologies of data intensive science for over a decade with large sky surveys and distributed efforts such as the Virtual Observatory. However, it faces a new data deluge with the next generation of synoptic sky surveys which are opening up the time domain for discovery and exploration. This brings both new scientific opportunities and fresh challenges, in terms of data rates from robotic telescopes and exponential complexity in linked data, but also for data mining algorithms used in classification and decision making. In this paper, we describe how an informatics-based approach-part of the so-called "fourth paradigm" of scientific discovery-is emerging to deal with these. We review our experiences with the Palomar-Quest and Catalina Real-Time Transient Sky Surveys; in particular, addressing the issue of the heterogeneity of data associated with transient astronomical events (and other sensor networks) and how to manage and analyze it.

We present a sample of 120 dust-reddened quasars identified by matching radio sources detected at 1.4 GHz in the FIRST survey with the near-infrared 2MASS catalog and color-selecting red sources. Optical and/or near-infrared spectroscopy provide broad wavelength sampling of their spectral energy distributions that we use to determine their reddening, characterized by E(B-V). We demonstrate that the reddening in these quasars is best-described by SMC-like dust. This sample spans a wide range in redshift and reddening (0.1 < z < 3, 0.1 < E(B-V) < 1.5), which we use to investigate the possible correlation of luminosity with reddening. At every redshift, dust-reddened quasars are intrinsically the most luminous quasars. We interpret this result in the context of merger-driven quasar/galaxy co-evolution where these reddened quasars are revealing an emergent phase during which the heavily obscured quasar is shedding its cocoon of dust prior to becoming a "normal" blue quasar. When correcting for extinction, we find that, depending on how the parent population is defined, these red quasars make up < 15-20% of the luminous quasar population. We estimate, based on the fraction of objects in this phase, that its duration is 15-20% as long as the unobscured, blue quasar phase.

The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

The time domain has been identified as one of the most important areas of astronomical research for the next decade. The Virtual Observatory is in the vanguard with dedicated tools and services that enable and facilitate the discovery, dissemination and analysis of time domain data. These range in scope from rapid notifications of time-critical astronomical transients to annotating long-term variables with the latest modeling results. In this paper, we will review the prior art in these areas and focus on the capabilities that the VAO is bringing to bear in support of time domain science. In particular, we will focus on the issues involved with the heterogeneous collections of (ancillary) data associated with astronomical transients, and the time series characterization and classification tools required by the next generation of sky surveys, such as LSST and SKA.

Sky surveys represent a fundamental data basis for astronomy. We use them to map in a systematic way the universe and its constituents, and to discover new types of objects or phenomena. We review the subject, with an emphasis on the wide-field imaging surveys, placing them in a broader scientific and historical context. Surveys are the largest data generators in astronomy, propelled by the advances in information and computation technology, and have transformed the ways in which astronomy is done. We describe the variety and the general properties of surveys, the ways in which they may be quantified and compared, and offer some figures of merit that can be used to compare their scientific discovery potential. Surveys enable a very wide range of science; that is perhaps their key unifying characteristic. As new domains of the observable parameter space open up thanks to the advances in technology, surveys are often the initial step in their exploration. Science can be done with the survey data alone or a combination of different surveys, or with a targeted follow-up of potentially interesting selected sources. Surveys can be used to generate large, statistical samples of objects that can be studied as populations, or as tracers of larger structures. They can be also used to discover or generate samples of rare or unusual objects, and may lead to discoveries of some previously unknown types. We discuss a general framework of parameter spaces that can be used for an assessment and comparison of different surveys, and the strategies for their scientific exploration. As we move into the Petascale regime, an effective processing and scientific exploitation of such large data sets and data streams poses many challenges, some of which may be addressed in the framework of Virtual Observatory and Astroinformatics, with a broader application of data mining and knowledge discovery technologies.

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

We report the discovery of a new ultra-faint globular cluster in the constellation of Ursa Minor, based on stellar photometry from the MegaCam imager at the Canada-France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, is located at a distance of 45 +/- 5 kpc and at a projected distance of only 45 arcmin from the center of the Ursa Minor dSph galaxy. Using a Maximum Likelihood technique we measure a half-light radius of 0.5 arcmin, or equivalently 7 pc and an ellipticity consistent with being zero. We estimate its absolute magnitude to be M_V=-0.4 +/- 0.9, which corresponds to L_V=120 (+160, -65) L_sun and we measure a heliocentric radial velocity of -137 +/- 4 km/s based on Keck/DEIMOS spectroscopy. This new satellite is separate from Ursa Minor by ~30 kpc and 110 km/s suggesting the cluster is not obviously associated with the dSph, despite the very close angular separation. Based on its photometric properties and structural parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along with Segue 3 this is one of the faintest stellar clusters known to date.

Astronomical data are gathered through a very large number of heterogeneous techniques and stored in very diversified and often incompatible data repositories. Moreover in the e-science environment, it is needed to integrate services across distributed, heterogeneous, dynamic "virtual organizations" formed by different resources within a single enterprise and/or external resource sharing and service provider relationships. The DAME/VONeural project, run jointly by the University Federico II, INAF (National Institute of Astrophysics) Astronomical Observatories of Napoli and the California Institute of Technology, aims at creating a single, sustainable, distributed e-infrastructure for data mining and exploration in massive data sets, to be offered to the astronomical (but not only) community as a web application. The framework makes use of distributed computing environments (e.g. S.Co.P.E.) and matches the international IVOA standards and requirements. The integration process is technically challenging due to the need of achieving a specific quality of service when running on top of different native platforms. In these terms, the result of the DAME/VO-Neural project effort will be a service-oriented architecture, obtained by using appropriate standards and incorporating Grid paradigms and restful Web services frameworks where needed, that will have as main target the integration of interdisciplinary distributed systems within and across organizational domains.

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This problem will grow by orders of magnitude with the next generation of surveys. We are exploring a variety of novel automated classification techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. We describe briefly some of the methods used.

The Catalina Real-time Transient Survey (CRTS) currently covers 33,000 deg^2 of the sky in search of transient astrophysical events, with time baselines ranging from 10 minutes to ~7 years. Data provided by the Catalina Sky Survey provides an unequaled baseline against which >4,000 unique optical transient events have been discovered and openly published in real-time. Here we highlight some of the discoveries of CRTS.

The time domain community wants robust and reliable tools to enable production of and subscription to community-endorsed event notification packets (VOEvent). The VAO Transient Facility (VTF) is being designed to be the premier brokering service for the community, both collecting and disseminating observations about time-critical astronomical transients but also supporting annotations and the application of intelligent machine-learning to those observations. This distinguishes two types of activity associated with the facility: core infrastructure and user services. In this paper, we will review the prior art in both areas and describe the planned capabilities of the VTF. In particular, we will focus on scalability and quality-of-service issues required by the next generation of sky surveys, such as LSST and SKA.

Synoptic sky surveys are becoming the largest data generators in astronomy, and they are opening a new research frontier, that touches essentially every field of astronomy. Opening of the time domain to a systematic exploration will strengthen our understanding of a number of interesting known phenomena, and may lead to the discoveries of as yet unknown ones. We describe some lessons learned over the past decade, and offer some ideas that may guide strategic considerations in planning and execution of the future synoptic sky surveys.

Exploration of the time domain - variable and transient objects and phenomena - is rapidly becoming a vibrant research frontier, touching on essentially every field of astronomy and astrophysics, from the Solar system to cosmology. Time domain astronomy is being enabled by the advent of the new generation of synoptic sky surveys that cover large areas on the sky repeatedly, and generating massive data streams. Their scientific exploration poses many challenges, driven mainly by the need for a real-time discovery, classification, and follow-up of the interesting events. Here we describe the Catalina Real-Time Transient Survey (CRTS), that discovers and publishes transient events at optical wavelengths in real time, thus benefiting the entire community. We describe some of the scientific results to date, and then focus on the challenges of the automated classification and prioritization of transient events. CRTS represents a scientific and a technological testbed and precursor for the larger surveys in the future, including the Large Synoptic Survey Telescope (LSST) and the Square Kilometer Array (SKA).

Active galactic nuclei (AGNs) with double-peaked [O III] lines are suspected to be sub-kpc or kpc-scale binary AGNs. However, pure gas kinematics can produce the same double-peaked line profile in spatially integrated spectra. Here we combine integral-field spectroscopy and high-resolution imaging of 42 double-peaked [O III] AGNs from Sloan Digital Sky Survey to investigate the constituents of the population. We find two binary AGNs where the line-splitting is driven by the orbital motion of the merging nuclei. Such objects account for only ~2% of the double-peaked AGNs. Almost all (~98%) of the double-peaked AGNs were selected because of gas kinematics; and half of those show spatially resolved narrow-line regions that extend 4-20 kpc from the nuclei. Serendipitously, we find two spectrally unresolved binary AGNs where gas kinematics produced the double-peaked [O III] lines. The relatively frequent serendipitous discoveries indicate that only ~1% of binary AGNs would appear double-peaked in Sloan spectra and 2.2_{-0.8}^{+2.5}% of all Sloan AGNs are binary AGNs. Therefore, the double-peaked sample does not offer much advantage over any other AGN samples in finding binary AGNs. The binary AGN fraction implies an elevated AGN duty cycle (8_{-3}^{+8}%), suggesting galaxy interactions enhance nuclear accretion. We illustrate that integral-field spectroscopy is crucial for identifying binary AGNs: several objects previously classified as "binary AGNs" with long-slit spectra are most likely single AGNs with extended narrow-line regions. The formation of extended narrow-line regions driven by radiation pressure is also discussed.

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between \sigma_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

We describe the development of a system for an automated, iterative, real-time classification of transient events discovered in synoptic sky surveys. The system under development incorporates a number of Machine Learning techniques, mostly using Bayesian approaches, due to the sparse nature, heterogeneity, and variable incompleteness of the available data. The classifications are improved iteratively as the new measurements are obtained. One novel feature is the development of an automated follow-up recommendation engine, that suggest those measurements that would be the most advantageous in terms of resolving classification ambiguities and/or characterization of the astrophysically most interesting objects, given a set of available follow-up assets and their cost functions. This illustrates the symbiotic relationship of astronomy and applied computer science through the emerging discipline of AstroInformatics.

We present the discovery and follow-up observations of SN 2008jb, a core-collapse supernova in the dwarf irregular galaxy ESO 302-14 at 9.6 Mpc. This transient was missed by galaxy-targeted surveys and was only found in archival optical images obtained by CRTS and ASAS. It was detected shortly after explosion and reached a bright optical maximum, Vmax = 13.6 mag (M_Vmax ~ -16.5). The shape of the light curve shows a plateau of 100 days, followed by a drop of 1.4 mag in V-band to a decline with the approximate Co 56 decay slope, consistent with 0.04 Msun of Ni 56 synthesized in the explosion. A spectrum obtained 2 years after explosion shows a broad, boxy Halpha emission line, which is unusual for type IIP supernovae. We detect the supernova in archival Spitzer and WISE images obtained 8-14 months after explosion, which show clear signs of warm dust emission. The dwarf irregular host galaxy has a low gas-phase oxygen abundance, 12 + log(O/H) = 8.2 (~1/5 Solar), similar to those of the SMC and the hosts of long gamma-ray bursts and luminous core-collapse supernovae. We study the host environment using GALEX FUV, R-band, and Halpha images and find that the supernova occurred in a large star-formation complex. The morphology of the Halpha emission appears as a large shell (R = 350 pc) surrounding the FUV and optical emission. We estimate an age of ~9 Myr and a total mass of ~2 x 10^5 Msun for the star-formation complex. These properties are consistent with the expanding Halpha supershells observed in well-studied nearby dwarf galaxies, which are tell-tale signs of feedback from the cumulative effect of massive star winds and supernovae. The age estimated for the star-forming region suggests a relatively high-mass progenitor star with initial mass of ~20 Msun. We discuss the implications of these findings in the study of core-collapse supernova progenitors. (Abridged)

The exponential growth of astronomical data collected by both ground based and space borne instruments has fostered the growth of Astroinformatics: a new discipline laying at the intersection between astronomy, applied computer science, and information and computation (ICT) technologies. At the very heart of Astroinformatics is a complex set of methodologies usually called Data Mining (DM) or Knowledge Discovery in Data Bases (KDD). In the astronomical domain, DM/KDD are still in a very early usage stage, even though new methods and tools are being continuously deployed in order to cope with the Massive Data Sets (MDS) that can only grow in the future. In this paper, we briefly outline some general problems encountered when applying DM/KDD methods to astrophysical problems, and describe the DAME (DAta Mining & Exploration) web application. While specifically tailored to work on MDS, DAME can be effectively applied also to smaller data sets. As an illustration, we describe two application of DAME to two different problems: the identification of candidate globular clusters in external galaxies, and the classification of active galactic nuclei (AGN). We believe that tools and services of this nature will become increasingly necessary for the data-intensive astronomy (and indeed all sciences) in the 21st century.

We report the confirmation of a kpc-scale binary active galactic nucleus (AGN) with high-resolution radio images from the Expanded Very Large Array (EVLA). SDSS J150243.1+111557 is a double-peaked [O III] AGN at z = 0.39 selected from the Sloan Digital Sky Survey. Our previous near-infrared adaptive optics imaging reveals two nuclei separated by 1.4" (7.4 kpc), and our optical integral-field spectroscopy suggests that they are a type-1--type-2 AGN pair. However, these data alone cannot rule out the single AGN scenario where the narrow emission-line region associated with the secondary is photoionized by the broad-line AGN in the primary. Our new EVLA images at 1.4, 5.0, and 8.5 GHz show two steep-spectrum compact radio sources spatially coincident with the optical nuclei. The radio power of the type-2 AGN is an order-of-magnitude in excess of star-forming galaxies with similar extinction-corrected [O II] 3727 luminosities, indicating that the radio emission is powered by accretion. Therefore, SDSS J150243.1+111557 is one of the few confirmed kpc-scale binary AGN systems. Spectral-energy-distribution modeling shows that SDSS J150243.1+111557 is a merger of two ~10^{11} M_sun galaxies. With both black hole masses around 10^8 Msun, the AGNs are accreting at ~10 times below the Eddington limit.

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

In course of the search for the optical identifications associated with ROSAT X-ray sources we have found a highly variable object with the very unusual long-term behavior, color indices and high X-ray-to-optical flux ratio. We report the archival photometric light curve from the Catalina Sky Survey, optical spectroscopy from RTT150 and time-resolved photometry from Astrotel-Caucasus telescope. The object appears to be the magnetic cataclysmic variable (polar) with orbital period of P=3.20 hr.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

As a natural consequence of galaxy mergers, binary active galactic nuclei (AGNs) should be commonplace. Nevertheless, observational confirmations are rare, especially for binaries with separations less than ten kpc. Such a system may show two sets of narrow emission lines in a single spectrum owing to the orbital motion of the binary. We have obtained high-resolution near-infrared images of 50 double-peaked [O III] 5007 AGNs with the Keck II laser guide star adaptive optics system. The Sloan Digital Sky Survey sample is compiled from the literature and consists of 17 type-1 AGNs between 0.18 < z < 0.56 and 33 type-2 AGNs between 0.03 < z < 0.24. The new images reveal eight type-1 and eight type-2 sources that are apparently undergoing mergers. These are strong candidates of kpc-scale binary AGNs, because they show multiple components separated between 0.6 and 12 kpc and often disturbed morphologies. Because most of the type-1s are at higher redshifts than the type-2s, the higher merger fraction of type-1s (47+/-20%) compared to that of type-2s (24+/-10%) can be attributed to the general evolution of galaxy merger fraction with redshift. Furthermore, we show that AGN mergers are outliers of the M_BH-sigma relation because of over-estimated stellar velocity dispersions, illustrating the importance of removing mergers from the samples defining the M_BH-sigma relations. Finally, we find that the emission-line properties are indistinguishable for spatially resolved and unresolved sources, emphasizing that scenarios involving a single AGN can produce the same double-peaked line profiles and they account for at least 70% of the double-peaked [O III] AGNs.

Catalina Real-Time Transient Survey (CRTS) is a synoptic sky survey uses data streams from 3 wide-field telescopes in Arizona and Australia, covering the total area of ~30,000 deg2, down to the limiting magnitudes ~ 20 - 21 mag per exposure, with time baselines from 10 min to 6 years (and growing); there are now typically ~ 200 - 300 exposures per pointing, and coadded images reach deeper than 23 mag. The basic goal of CRTS is a systematic exploration and characterization of the faint, variable sky. The survey has detected ~ 3,000 high-amplitude transients to date, including ~ 1,000 supernovae, hundreds of CVs (the majority of them previously uncatalogued), and hundreds of blazars / OVV AGN, highly variable and flare stars, etc. CRTS has a complete open data philosophy: all transients are published immediately electronically, with no proprietary period at all, and all of the data (images, light curves) will be publicly available in the near future, thus benefiting the entire astronomical community. CRTS is a scientific and technological testbed and precursor for the grander synoptic sky surveys to come.

We present an updated determination of the z ~ 4 QSO luminosity function (QLF), improving the quality of the determination of the faint end of the QLF presented in Glikman et al. (2010). We have observed an additional 43 candidates from our survey sample, yielding one additional QSO at z = 4.23 and increasing the completeness of our spectroscopic follow-up to 48% for candidates brighter than R = 24 over our survey area of 3.76 deg2. We study the effect of using K-corrections to compute the rest-frame absolute magnitude at 1450A compared with measuring M1450 directly from the object spectra. We find a luminosity-dependent bias: template-based K-corrections overestimate the luminosity of low-luminosity QSOs, likely due to their reliance on templates derived from higher luminosity QSOs. Combining our sample with bright quasars from the Sloan Digital Sky Survey and using spectrum-based M1450 for all the quasars, we fit a double-power-law to the binned QLF. Our best fit has a bright-end slope, {\alpha} = 3.3\pm0.2, and faint-end slope, {\beta} = 1.6(+0.8/-0.6). Our new data revise the faint-end slope of the QLF down to flatter values similar to those measured at z ~ 3. The break luminosity, though poorly constrained, is at M* = -24.1(+0.7/-1.9), approximately 1 - 1.5 mag fainter than at z ~ 3. This QLF implies that QSOs account for about half the radiation needed to ionize the IGM at these redshifts.

Although white dwarfs are believed to be the end point of most stellar evolution, unlike main sequence stars, they have not yet been the subject of dedicated time-domain surveys for exoplanets. We discuss how their size and distinctive colour make them excellent targets for wide-field searches for exoplanets. In particular, we note that planets of Earth-size can give rise to multi-magnitude eclipses of massive white dwarfs. Such a large signal is almost unmistakable and would be detectable even with very low-precision photometry. For objects of smaller size, the high accuracy photometry currently being used to detect Super-Earth and smaller planets transiting Sun-sized stars, is capable of revealing minor planets down to R~100km as they transit white dwarfs. Such observations can be used to test current evidence for asteroid-size objects being the cause for dust rings which have recently been observed for a number of white dwarfs. No other current exoplanet search method is capable of detecting such exo-asteroids. As an initial test of this search strategy, we combine synoptic data from the Catalina Sky Survey with multi-colour photometry and spectra from the Sloan Digital Sky Survey to search ~12,000 white dwarf lightcurves for eclipsing events. We find 20 new eclipsing white dwarf binary systems with low-mass companions. This doubles the number of known eclipsing white dwarfs and is expected to enable the determination of accurate white dwarf radii. Three of the discoveries have radii consistent with substellar systems and show no evidence of flux from the eclipsing object in their SDSS optical spectra, or near-IR data.

Cosmic collisions on planets cause detectable optical flashes that range from terrestrial shooting stars to bright fireballs. On June 3, 2010 a bolide in Jupiter's atmosphere was simultaneously observed from the Earth by two amateur astronomers observing Jupiter in red and blue wavelengths. The bolide appeared as a flash of 2 s duration in video recording data of the planet. The analysis of the light curve of the observations results in an estimated energy of the impact of 0.9-4.0x10^{15} J which corresponds to a colliding body of 8-13 m diameter assuming a mean density of 2 g cm^{-3}. Images acquired a few days later by the Hubble Space Telescope and other large ground-based facilities did not show any signature of aerosol debris, temperature or chemical composition anomaly, confirming that the body was small and destroyed in Jupiter's upper atmosphere. Several collisions of this size may happen on Jupiter on a yearly basis. A systematic study of the impact rate and size of these bolides can enable an empirical determination of the flux of meteoroids in Jupiter with implications for the populations of small bodies in the outer Solar System and may allow a better quantification of the threat of impacting bodies to Earth. The serendipitous recording of this optical flash opens a new window in the observation of Jupiter with small telescopes.

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived from echellete spectroscopy of binary quasars. We find CIV systems at similar redshifts in paired sightlines more often than expected for sightline-independent redshifts. As the separation of the sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest separations, the CIV systems cluster similar to Lyman-break galaxies (Adelberger et al. 2005a). The CIV systems are significantly less correlated than these galaxies, however, at separations less than R_1 ~ 0.42 +/- 0.15 h-1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this length scale is significantly smaller than the break scale estimated from the line-of-sight correlation function in redshift space (Scannapieco et al. 2006a). Using a simple model, we interpret the new real-space measurement as an indication of the typical physical size of enriched regions. We adopt this size for enriched regions and fit the redshift-space distortion in the line-of-sight correlation function. The fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying mass distribution and places an upper limit on the outflow (or inflow) speed of metals. The implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies (Shapley et al. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints on the evolution of the CIV correlation function with cosmic time. This study demonstrates the potential of tomography for measuring the metal enrichment history of the circum-galactic medium.

We present the first case of strong gravitational lensing by a QSO : SDSS J0013+1523, at z = 0.120. The discovery is the result of a systematic search for emission lines redshifted behind QSOs, among 22298 spectra of the SDSS data release 7. Apart from the z = 0.120 spectral features of the foreground QSO, the spectrum of SDSS J0013+1523 also displays the OII and Hbeta emission lines and the OIII doublet, all at the same redshift, z = 0.640. Using sharp Keck adaptive optics K-band images obtained using laser guide stars, we unveil two objects within a radius of 2 arcsec from the QSO. Deep Keck optical spectroscopy clearly confirms one of these objects at z = 0.640 and shows traces of the OIII, emission line of the second object, also at z = 0.640. Lens modeling suggests that they represent two images of the same z = 0.640 emission-line galaxy. Our Keck spectra also allow us to measure the redshift of an intervening galaxy at z = 0.394, located 3.2 arcsec away from the line of sight to the QSO. If the z = 0.120 QSO host galaxy is modeled as a singular isothermal sphere, its mass within the Einstein radius is M_E(r < 1 kpc) = 2.16e10 M_Sun and its velocity dispersion is sigma_SIS = 169 km/s. This is about 1 sigma away from the velocity dispersion estimated from the width of the QSO Hbeta emission line, sigma_*(M_BH) = 124 +/- 47 km/s. Deep optical HST imaging will be necessary to constrain the total radial mass profile of the QSO host galaxy using the detailed shape of the lensed source. This first case of a QSO acting as a strong lens on a more distant object opens new directions in the study of QSO host galaxies.

We have conducted a spectroscopic survey to find faint quasars (-26.0 < M_{1450} < -22.0) at redshifts z=3.8-5.2 in order to measure the faint end of the quasar luminosity function at these early times. Using available optical imaging data from portions of the NOAO Deep Wide-Field Survey and the Deep Lens Survey, we have color-selected quasar candidates in a total area of 3.76 deg^2. Thirty candidates have R <= 23 mags. We conducted spectroscopic followup for 28 of our candidates and found 23 QSOs, 21 of which are reported here for the first time, in the 3.74 < z <5.06 redshift range. We estimate our survey completeness through detailed Monte Carlo simulations and derive the first measurement of the density of quasars in this magnitude and redshift interval. We find that the binned luminosity function is somewhat affected by the K-correction used to compute the rest-frame absolute magnitude at 1450A. Considering only our R <= 23 sample, the best-fit single power-law (Phi \propto L^beta) gives a faint-end slope beta = -1.6+/-0.2. If we consider our larger, but highly incomplete sample going one magnitude fainter, we measure a steeper faint-end slope -2 < beta < -2.5. In all cases, we consistently find faint-end slopes that are steeper than expected based on measurements at z ~ 3. We combine our sample with bright quasars from the Sloan Digital Sky Survey to derive parameters for a double-power-law luminosity function. Our best fit finds a bright-end slope, alpha = -2.4+/-0.2, and faint-end slope, beta = -2.3+/-0.2, without a well-constrained break luminosity. This is effectively a single power-law, with beta = -2.7+/-0.1. We use these results to place limits on the amount of ultraviolet radiation produced by quasars and find that quasars are able to ionize the intergalactic medium at these redshifts.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We report on the discovery and initial observations of the energetic type IIn supernova (SN), 2008fz. The optical energy emitted by SN 2008fz (based on the light curve over a 88 day period), is possibly the most ever observed for a supernova (1.4 x 10^51 erg). The event was more luminous than the type IIn SN 2006gy, but exhibited same smooth, slowly evolving light curve. As is characteristic of type IIn SN, the early spectra of 2008fz initially exhibited narrow Balmer lines which were replaced by a broader component at later times. The spectra also show a blue continuum with no signs of Ca or Na absorption, suggesting that there is little extinction due to intragalatic dust in the host or circumstellar material. No host galaxy is identified in prior coadded images reaching R ~ 22. From the supernova's redshift, z=0.133, we place an upper limit on the host of M_R=-17. The presence of the SN within such a faint host follows the majority of recently discovered highly luminous SN. A possible reason for this occurrence is the very high star formation rate occurring in low-mass galaxies in combination with the low metallicity environment, which makes the production of very massive stars possible. We determine the peak absolute magnitude of the event to be M_V = -22.3 from the initial photometry and the redshift distance, placing it among the most luminous supernovae discovered.

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables.

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an experiment in academic and scientific practices enabled by the immersive VR technologies. We describe the current and planned activities and research directions of MICA, and offer some thoughts as to what the future developments in this arena may be.

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data.

Skyalert.org is a web application to collect and disseminate observations about time-critical astronomical transients, and to add annotations and intelligent machine-learning to those observations. The information is "pushed" to subscribers, who may be either humans (email, text message etc) or they may be machines that control telescopes. Subscribers can prepare precise "trigger rules" to decide which events should reach them and their robots, rules that may be based on sky position, or on the specific vocabulary of parameters that define a particular type of observation. Our twin thrusts are automation of process, and discrimination of interesting events.

Precise radial velocity measurements from HIRES on the Keck I telescope are presented for 40 stars in the outer halo globular cluster NGC 2419. These data are used to probe the cluster's stellar mass function and search for the presence of dark matter in this cluster. NGC 2419 is one of the best Galactic globular clusters for such a study due to its long relaxation time (T_{r0} ~ 10^{10} yr) and large Galactocentric distance (R_{GC} ~ 90 kpc) -- properties that make significant evolutionary changes in the low-mass end of the cluster mass function unlikely. We find a mean cluster velocity of <v_r>=-20.3 +- 0.7 km/sec and an internal velocity dispersion of \sigma = 4.14 +- 0.48 km/sec, leading to a total mass of (9.0 +- 2.2) * 10^5 Msun and a global mass-to-light ratio of M/L_V = 2.05 +- 0.50 in solar units. This mass-to-light ratio is in good agreement with what one would expect for a pure stellar system following a standard mass function at the metallicity of NGC 2419. In addition, the mass-to-light ratio does not appear to rise towards the outer parts of the cluster. Our measurements therefore rule out the presence of a dark matter halo with mass larger than ~10^7 Msun inside the central 500 pc, which is lower than what is found for the central dark matter densities of dSph galaxies. We also discuss the relevance of our measurements for alternative gravitational theories such as MOND, and for possible formation scenarios of ultra-compact dwarf galaxies.

As part of the conceptual and preliminary design processes of the Thirty Meter Telescope (TMT), the TMT site testing team has spent the last five years measuring the atmospheric properties of five candidate mountains in North and South America with an unprecedented array of instrumentation. The site testing period was preceded by several years of analyses selecting the five candidates, Cerros Tolar, Armazones and Tolonchar in northern Chile; San Pedro Martir in Baja California, Mexico and the 13 North (13N) site on Mauna Kea, Hawaii. Site testing was concluded by the selection of two remaining sites for further consideration, Armazones and Mauna Kea 13N. It showed that all five candidates are excellent sites for an extremely large astronomical observatory and that none of the sites stands out as the obvious and only logical choice based on its combined properties. This is the first article in a series discussing the TMT site testing project.

We use the distant outer halo globular cluster Palomar 14 as a test case for classical vs. modified Newtonian dynamics (MOND). Previous theoretical calculations have shown that the line-of-sight velocity dispersion predicted by these theories can differ by up to a factor of three for such sparse, remote clusters like Pal 14. We determine the line-of-sight velocity dispersion of Palomar 14 by measuring radial velocities of 17 red giant cluster members obtained using the Very Large Telescope (VLT) and Keck telescope. The systemic velocity of Palomar 14 is (72.28+-0.12)km/s. The derived velocity dispersion of (0.38+-0.12)km/s of the 16 definite member stars is in agreement with the theoretical prediction for the classical Newtonian case according to Baumgardt et al. (2005). In order to exclude the possibility that a peculiar mass function might have influenced our measurements, we derived the cluster's main sequence mass function down to 0.53Msolar using archival images obtained with the Hubble Space Telescope. We found a mass function slope of 1.27+-0.44, which is, compared to the canonical mass function, a significantly shallower slope. The derived lower limit on the cluster's mass is higher than the theoretically predicted mass in case of MOND. Our data are consistent with a central density of 0.1 Msolar pc^-3. We need no dark matter in Palomar 14. If the cluster is on a circular orbit, our spectroscopic and photometric results argue against MOND, unless this cluster experienced significant mass loss.

We present the complete galaxy cluster catalog from the Northern Sky Optical Cluster Survey, a new, objectively defined catalog of candidate galaxy clusters at z<0.25 drawn from the Digitized Second Palomar Observatory Sky Survey (DPOSS). The data presented here cover the Southern Galactic Cap, as well as the less-well calibrated regions of the Northern Galactic Cap. In addition, due to improvements in our cluster finder and measurement methods, we provide an updated catalog for the well-calibrated Northern Galactic Cap region previously published in Paper II. The complete survey covers 11,411 square degrees, with over 15,000 candidate clusters. We discuss improved photometric redshifts, richnesses and optical luminosities which are provided for each cluster. A variety of substructure measures are computed for a subset of over 11,000 clusters. We also discuss the derivation of dynamical radii r_200 and its relation to cluster richness. A number of consistency checks between the three areas of the survey are also presented, demonstrating the homogeneity of the catalog over disjoint sky areas. We perform extensive comparisons to existing optically and X-ray selected cluster catalogs, and derive new X-ray luminosities and temperatures for a subset of our clusters. We find that the optical and X-ray luminosities are well correlated, even using relatively shallow ROSAT All Sky Survey and DPOSS data. This survey provides a good comparison sample to the MaxBCG catalog based on Sloan Digital Sky Survey Data, and complements that survey at low redshifts 0.07<z<0.1.

Physical sizes of extended radio galaxies can be employed as a cosmological "standard ruler", using a previously developed method. Eleven new radio galaxies are added to our previous sample of nineteen sources, forming a sample of thirty objects with redshifts between 0 and 1.8. This sample of radio galaxies are used to obtain the best fit cosmological parameters in a quintessence model in a spatially flat universe, a cosmological constant model that allows for non-zero space curvature, and a rolling scalar field model in a spatially flat universe. Results obtained with radio galaxies are compared with those obtained with different supernova samples, and with combined radio galaxy and supernova samples. Results obtained with different samples are consistent, suggesting that neither method is seriously affected by systematic errors. Best fit radio galaxy and supernovae model parameters determined in the different cosmological models are nearly identical, and are used to determine dimensionless coordinate distances to supernovae and radio galaxies, and distance moduli to the radio galaxies. The distance moduli to the radio galaxies can be combined with supernovae samples to increase the number of sources, particularly high-redshift sources, in the samples. The constraints obtained here with the combined radio galaxy plus supernovae data set in the rolling scalar field model are quite strong. The best fit parameter values suggest a value of omega is less than about 0.35, and the model parameter alpha is close to zero; that is, a cosmological constant provides a good description of the data. We also obtain new constraints on the physics of engines that power the large-scale radio emission.

We present a sample of $i_{775}$-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Survey QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distributions from that of GOODS at the 99% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of $i_{775}$-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 96% confidence level, depending on which selection criteria were adopted to identify $i_{775}$-dropouts; its width cannot be explained by cosmic variance alone. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest redshift QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences between lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.

Digital synoptic sky surveys pose several new object classification challenges. In surveys where real-time detection and classification of transient events is a science driver, there is a need for an effective elimination of instrument-related artifacts which can masquerade as transient sources in the detection pipeline, e.g., unremoved large cosmic rays, saturation trails, reflections, crosstalk artifacts, etc. We have implemented such an Artifact Filter, using a supervised neural network, for the real-time processing pipeline in the Palomar-Quest (PQ) survey. After the training phase, for each object it takes as input a set of measured morphological parameters and returns the probability of it being a real object. Despite the relatively low number of training cases for many kinds of artifacts, the overall artifact classification rate is around 90%, with no genuine transients misclassified during our real-time scans. Another question is how to assign an optimal star-galaxy classification in a multi-pass survey, where seeing and other conditions change between different epochs, potentially producing inconsistent classifications for the same object. We have implemented a star/galaxy multipass classifier that makes use of external and a priori knowledge to find the optimal classification from the individually derived ones. Both these techniques can be applied to other, similar surveys and data sets.

Exploration of time domain is now a vibrant area of research in astronomy, driven by the advent of digital synoptic sky surveys. While panoramic surveys can detect variable or transient events, typically some follow-up observations are needed; for short-lived phenomena, a rapid response is essential. Ability to automatically classify and prioritize transient events for follow-up studies becomes critical as the data rates increase. We have been developing such methods using the data streams from the Palomar-Quest survey, the Catalina Sky Survey and others, using the VOEventNet framework. The goal is to automatically classify transient events, using the new measurements, combined with archival data (previous and multi-wavelength measurements), and contextual information (e.g., Galactic or ecliptic latitude, presence of a possible host galaxy nearby, etc.); and to iterate them dynamically as the follow-up data come in (e.g., light curves or colors). We have been investigating Bayesian methodologies for classification, as well as discriminated follow-up to optimize the use of available resources, including Naive Bayesian approach, and the non-parametric Gaussian process regression. We will also be deploying variants of the traditional machine learning techniques such as Neural Nets and Support Vector Machines on datasets of reliably classified transients as they build up.

We report on the results from the first six months of the Catalina Real-time Transient Survey (CRTS). In order to search for optical transients with timescales of minutes to years, the CRTS analyses data from the Catalina Sky Survey which repeatedly covers twenty six thousand of square degrees on the sky. The CRTS provides a public stream of transients that are bright enough to be followed up using small telescopes. Since the beginning of the survey, all CRTS transients have been made available to astronomers around the world in real-time using HTML tables, RSS feeds and VOEvents. As part of our public outreach program the detections are now also available in KML through Google Sky. The initial discoveries include over 350 unique optical transients rising more than two magnitudes from past measurements. Sixty two of these are classified as supernovae, based on light curves, prior deep imaging and spectroscopic data. Seventy seven are due to cataclysmic variables (only 13 previously known), while an additional 100 transients were too infrequently sampled to distinguish between faint CVs and SNe. The remaining optical transients include AGN, Blazars, high proper motions stars, highly variable stars (such as UV Ceti stars) and transients of an unknown nature. Our results suggest that there is a large population of SNe missed by many current supernova surveys because of selection biases. These objects appear to be associated with faint host galaxies. We also discuss the unexpected discovery of white dwarf binary systems through dramatic eclipses.

We conducted an exploratory search for quasars at z~ 6 - 8, using the Early Data Release from United Kingdom Infrared Deep Sky survey (UKIDSS) cross-matched to panoramic optical imagery. High redshift quasar candidates are chosen using multi-color selection in i,z,Y,J,H and K bands. After removal of apparent instrumental artifacts, our candidate list consisted of 34 objects. We further refined this list with deeper imaging in the optical for ten of our candidates. Twenty-five candidates were followed up spectroscopically in the near-infrared and in the optical. We confirmed twenty-five of our spectra as very low-mass main-sequence stars or brown dwarfs, which were indeed expected as the main contaminants of this exploratory search. The lack of quasar detection is not surprising: the estimated probability of finding a single z>6 quasar down to the limit of UKIDSS in the 27.3 square degrees of the EDR is <5%. We find that the most important limiting factor in this work is the depth of the available optical data. Experience gained in this pilot project can help refine high-redshift quasar selection criteria for subsequent UKIDSS data releases.

The relationship between galaxies and supermassive black holes (SMBH) found in their cores plays a key role in the formation and evolution of both of these major constituents of the universe, as well as the evolution of the intergalactic medium. Neither can be fully understood on their own, and studies of galaxy and SMBH co-formation and co-evolution are now among the central topics of research in cosmology. Yet the very origins, and the early growth phases of the SMBH are still not firmly established. We review our current understanding of the relevant processes and their astrophysical and cosmological context, with an emphasis on the observability of the SMBH growth mechanisms at high redshifts, and their leftover progeny at low redshifts.

There is an increasing number of large, digital, synoptic sky surveys, in which repeated observations are obtained over large areas of the sky in multiple epochs. Likewise, there is a growth in the number of (often automated or robotic) follow-up facilities with varied capabilities in terms of instruments, depth, cadence, wavelengths, etc., most of which are geared toward some specific astrophysical phenomenon. As the number of detected transient events grows, an automated, probabilistic classification of the detected variables and transients becomes increasingly important, so that an optimal use can be made of follow-up facilities, without unnecessary duplication of effort. We describe a methodology now under development for a prototype event classification system; it involves Bayesian and Machine Learning classifiers, automated incorporation of feedback from follow-up observations, and discriminated or directed follow-up requests. This type of methodology may be essential for the massive synoptic sky surveys in the future.

We describe briefly the Palomar-Quest (PQ) digital synoptic sky survey, including its parameters, data processing, status, and plans. Exploration of the time domain is now the central scientific and technological focus of the survey. To this end, we have developed a real-time pipeline for detection of transient sources. We describe some of the early results, and lessons learned which may be useful for other, similar projects, and time-domain astronomy in general. Finally, we discuss some issues and challenges posed by the real-time analysis and scientific exploitation of massive data streams from modern synoptic sky surveys.

The model-independent method of using type Ia supernovae proposed and developed by Daly and Djorgovski (2003, 2004) has been applied to the Riess et al. (2007) supernovae sample. Assuming only a Robertson-Walker metric, we find that the universe is accelerating today. This result is purely kinematic, is independent of the contents of the universe, and does not require that a theory of gravity be specified. Our model-independent method allows a determination of q(z) for a particular value of space curvature. When q(z) transitions from negative to positive values, the universe transitions from an accelerating to a decelerating state. For zero space curvature, we find that the universe transitions from acceleration to deceleration at a zedshift of about = 0.35 for the Riess et al. (2007) sample. If a theory of gravity is specified, the supernovae data can be used to determine the pressure, energy density, and equation of state of the dark energy, and the potential and kinetic energy density of a dark energy scalar field as functions of redshift. The relevant equations from General Relativity are applied, and these functions are obtained. The results are consistent with predictions in the standard Lambda Cold Dark Matter model at about the two sigma level.

We extend and apply a model-independent analysis method developed earlier by Daly & Djorgovski to new samples of supernova standard candles, radio galaxy and cluster standard rulers, and use it to constrain physical properties of the dark energy as functions of redshift. Similar results are obtained for the radio galaxy and supernova data sets. The first and second derivatives of the distance are compared directly with predictions in a standard model based on General Relativity. The good agreement indicates that General Relativity provides an accurate description of the data on look-back time scales of about ten billion years. The first and second derivatives are combined to obtain the acceleration parameter, assuming only the validity of the Robertson-Walker metric, independent of a theory of gravity and of the physical nature of the dark energy. The acceleration of the universe at the current epoch is indicated by the analysis. The effect of non-zero space curvature on q(z) is explored. We solve for the pressure, energy density, equation of state, and potential and kinetic energy of the dark energy as functions of redshift assuming that General Relativity is the correct theory of gravity, and the results indicate that a cosmological constant in a spatially flat universe provides a good description of each of these quantities over the redshift range from zero to about one. We define a new function, the dark energy indicator, in terms of the first and second derivatives of the coordinate distance and show how this can be used to measure deviations of w from -1 and to obtain a new and independent measure of Omega.

In this paper we present multiband optical and UV Hubble Space Telescope photometry of the two Galactic globular clusters NGC 6388 and NGC 6441. We have obtained the first UV color-magnitude diagrams for NGC 6388 and NGC 6441. These diagrams confirm previous results, obtained in optical bands, about the presence of a sizeable stellar population of extremely hot horizontal branch stars. At least in NGC 6388, we find a clear indication that at the hot end of the horizontal branch the distribution of stars forms a hook-like feature, closely resembling those observed in NGC 2808 and Omega Centauri. We briefly review the theoretical scenarios that have been suggested for interpreting this observational feature. We also investigate the tilted horizontal branch morphology and provide further evidence that supports early suggestions that this feature cannot be interpreted as an effect of differential reddening. We show that a possible solution of the puzzle is to assume that a small fraction - ranging between 10-20% - of the stellar population in the two clusters is strongly helium enriched (Y=0.40 in NGC 6388 and Y=0.35 in NGC 6441). The occurrence of a spread in the He abundance between the canonical value (Y=0.26) and the quoted upper limits can significantly help in explaining the whole morphology of the horizontal branch and the pulsational properties of the variable stars in the target clusters.

We report the discovery of two low-luminosity quasars at z~4, both of which show prominent N IV] 1486A emission. This line is extremely rare in quasar spectra at any redshift; detecting it in two out of a sample of 23 objects (i.e., ~ 9% of the sample) is intriguing and is likely due to the low-luminosity, high-redshift quasar sample we are studying. This is still a poorly explored regime, where contributions from associated, early starbursts may be significant. One interpretation of this line posits photoionization by very massive young stars. Seeing N IV] 1486A emission in a high-redshift quasar may thus be understood in the context of co-formation and early co-evolution of galaxies and their supermassive black holes. Alternatively, we may be seeing a phenomenon related to the early evolution of quasar broad emission line regions. The non-detection (and possibly even broad absorption) of N V 1240A line in the spectrum of one of these quasars may support that interpretation. These two objects may signal a new faint quasar population or an early AGN evolutionary stage at high redshifts.

We report the discovery of the first known probable case of a physical triple quasar (not a gravitational lens). A previously known double system, QQ 1429-008 at z = 2.076, is shown to contain a third, fainter QSO component at the same redshift within the measurement errors. Deep optical and IR imaging at the Keck and VLT telescopes has failed to reveal a plausible lensing galaxy group or a cluster, and moreover, we are unable to construct any viable lensing model which could lead to the observed distribution of source positions and relative intensities of the three QSO image components. Furthermore, there are hints of differences in broad-band spectral energy distributions of different components, which are more naturally understood if they are physically distinct AGN. Therefore, we conclude that this system is most likely a physical triple quasar, the first such close QSO grouping known at any redshift. The projected component separations in the restframe are ~ 30 - 50 kpc for the standard concordance cosmology, typical of interacting galaxy systems. The existence of this highly unusual system supports the standard picture in which galaxy interactions lead to the onset of QSO activity.

In astronomy multiple images are frequently obtained at the same position of the sky for follow-up co-addition as it helps one go deeper and look for fainter objects. With large scale panchromatic synoptic surveys becoming more common, image co-addition has become even more necessary as new observations start to get compared with co-added fiducial sky in real time. The standard co-addition techniques have included straight averages, variance weighted averages, medians etc. A more sophisticated nonlinear response chi-square method is also used when it is known that the data are background noise limited and the point spread function is homogenized in all channels. A more robust object detection technique capable of detecting faint sources, even those not seen at all epochs which will normally be smoothed out in traditional methods, is described. The analysis at each pixel level is based on a formula similar to Mahalanobis distance. The method does not depend on the point spread function.

We describe briefly a novel interpretation of the physical nature of dark energy (DE), based on the vacuum fluctuations model by Gurzadyan & Xue, and describe an internally consistent solution for the behavor of DE as a function of redshift. A key choice is the nature of the upper bound used for the computation of energy density contributions by vacuum modes. We show that use of the comoving horizon radius produces a viable model, whereas use of the proper horizon radius is inconsistent with the observations. After introduction of a single phenomenological parameter, the model is consistent with all of the curently available data, and fits them as well as the standard cosmological constant model, while making testable predictions. While some substantial interpretative uncertainties remain, future developments of this model may lead to significant new insights into the physical nature of DE.

Understanding the acceleration of the universe and its cause is one of the key problems in physics and cosmology today, and is best studied using a variety of mutually complementary approaches. Daly and Djorgovski (2003, 2004) proposed a model independent approach to determine the expansion and acceleration history of the universe and a number of important physical parameters of the dark energy as functions of redshift directly from the data. Here, we apply the method to explicitly determine the first and second derivatives of the coordinate distance with respect to redshift and combine them to solve for the acceleration of the universe and the kinetic and potential energy density of the dark energy as functions of redshift. A data set of 228 supernova and 20 radio galaxy measurements with redshifts from zero to 1.79 is used for this study. The values we obtain are shown to be consistent with the values expected in a standard Lambda Cold Dark Matter model.

We have mosaiced 2MASS images to derive surface brightness profiles in JHK for 104 Galactic globular clusters. We fit these with King profiles, and show that the core radii are identical to within the errors for each of these IR colors, and are identical to the core radii at V in essentially all cases. We derive integrated light colors V-J, V-H, V-K_s, J-H and J-K_s for these globular clusters. Each color shows a reasonably tight relation between the dereddened colors and metallicity. Fits to these are given for each color. The IR--IR colors have very small errors due largely to the all-sky photometric calibration of the 2MASS survey, while the V-IR colors have substantially larger uncertainties. We find fairly good agreement with measurements of integrated light colors for a smaller sample of Galactic globular clusters by Aaronson, Malkan & Kleinmann from 1977. Our results provide a calibration for the integrated light of distant single burst old stellar populations from very low to Solar metallicities. A comparison of our dereddened measured colors with predictions from several models of the integrated light of single burst old populations shows good agreement in the low metallicity domain for V-K_s colors, but an offset at a fixed [Fe/H] of ~0.1 mag in J-K_s, which we ascribe to photometric system transformation issues. Some of the models fail to reproduce the behavior of the integrated light colors of the Galactic globular clusters near Solar metallicity.

We present optical and near-IR observations of the host galaxy of GRB 020127, for which we measure R-K=6.2 mag. This is only the second GRB host to date, which is classified as an ERO. The spectral energy distribution is typical of a dusty starburst galaxy, with a redshift, z~1.9, a luminosity, L~5L*, and an inferred stellar mass of ~10^11-10^12 Msun, two orders of magnitude more massive than typical GRB hosts. A comparison to the z~2 mass-metallicity (M-Z) relation indicates that the host metallicity is about 0.5 Zsun. This result shows that at least some GRBs occur in massive, metal-enriched galaxies, and that the proposed low metallicity bias of GRB progenitors is not as severe as previously claimed. Instead we propose that the blue colors and sub-L* luminosities of most GRB hosts reflect their young starburst populations. This, along with the locally increased fraction of starbursts at lower stellar mass, may in fact be the underlying reason for the claimed low metallicity bias of z<0.2 GRB hosts. Consequently, GRBs and their host galaxies may serve as a reliable tracer of cosmic star formation, particularly at z>1 where the M-Z relation is systematically lower, and the star formation rate is dominated by sub-L* galaxies.

We review some of the recent developments and challenges posed by the data analysis in modern digital sky surveys, which are representative of the information-rich astronomy in the context of Virtual Observatory. Illustrative examples include the problems of an automated star-galaxy classification in complex and heterogeneous panoramic imaging data sets, and an automated, iterative, dynamical classification of transient events detected in synoptic sky surveys. These problems offer good opportunities for productive collaborations between astronomers and applied computer scientists and statisticians, and are representative of the kind of challenges now present in all data-intensive fields. We discuss briefly some emergent types of scalable scientific data analysis systems with a broad applicability.

We present a comparison of optical and X-ray properties of galaxy clusters in the northern sky. We determine the recovery rate of X-ray detected clusters in the optical as a function of richness, redshift and X-ray luminosity, showing that the missed clusters are typically low contrast systems when observed optically. We employ four different statistical tests to test for the presence of substructure using optical two-dimensional data, finding that approximately 35% of the clusters show strong signs of substructure. However, the results are test-dependent, with variations also due to the magnitude range and radius utilized.We have also performed a comparison of X-ray luminosity and temperature with optical galaxy counts (richness). We find that the slope and scatter of the relations between richness and the X-ray properties are heavily dependent on the density contrast of the clusters. The selection of substructure-free systems does not improve the correlation between X-ray luminosity and richness, but this comparison also shows much larger scatter than one obtained using the X-ray temperature. In the latter case, the sample is significantly reduced because temperature measurements are available only for the most massive (and thus high contrast) systems. However, the comparison between temperature and richness is very sensitive to the exclusion of clusters showing signs of substructure. The correlation of X-ray luminosity and richness is based on the largest sample to date ($\sim$ 750 clusters), while tests involving temperature use a similar number of objects as previous works ($\lsim$100). The results presented here are in good agreement with existing literature.

We present deep Hubble Space Telescope and Spitzer Space Telescope observations of the host galaxy of GRB 050904 at z=6.295. The host is detected in the H-band and marginally at 3.6 micron. From these detections, and limits in the z'-band and 4.5 micron, we infer an extinction-corrected absolute magnitude, M(UV)=-20.7 mag, or ~L*, a substantial star formation rate of 15 solar masses per year, and a stellar mass of a few 10^9 solar masses. A comparison to the published sample of spectroscopically-confirmed galaxies at z>5.5 reveals that the host of GRB 050904 would evade detection and/or confirmation in any of the current surveys due to the lack of detectable Ly-alpha emission, which is likely the result of dust extinction (A[1200]~1.5 mag). This suggests that not all luminous starburst galaxies at z~6 are currently being accounted for. Most importantly, using the metallicity of Z~0.05 solar inferred from the afterglow absorption spectrum, our observations indicate for the first time that the observed evolution in the mass- and luminosity-metallicity relations from z=0 to z~2 continues on to z>6. The ease of measuring redshifts and metallicities from the afterglow emission suggests that in tandem with the next generation ground- and space-based telescopes, a GRB mission with dedicated near-IR follow-up can provide unique information on the evolution of stars and galaxies through the epoch of re-ionization.

Understanding the physical nature of the dark energy which appears to drive the accelerated expansion of the unvierse is one of the key problems in physics and cosmology today. This important problem is best studied using a variety of mutually complementary approaches. Daly and Djorgovski (2003, 2004) proposed a model independent approach to determine a number of important physical parameters of the dark energy as functions of redshift directly from the data. Here, we expand this method to include the determinations of its potential and kinetic energy as functions of redshift. We show that the dark energy potential and kinetic energy may be written as combinations of the first and second derivatives of the coordinate distance with respect to redshift. We expand the data set to include new supernova measurements, and now use a total of 248 coordinate distances that span the redshift range from zero to 1.79. First and second derivatives of the coordinate distance are obtained as functions of redshift, and these are combined to determine the potential and kinetic energy of the dark energy as functions of redshift. An update on the redshift behavior of the dimensionless expansion rate E(z), the acceleration rate q(z), and the dark energy pressure p(z), energy density f(z), and equation of state w(z) is also presented. We find that the standard Omega = 0.3 and Lambda = 0.7 model is in an excellent agreement with the data. We also show tentative evidence that the Cardassian and Chaplygin gas models in a spatially flat universe do not fit the data as well.

The interpretation of globular cluster horizontal branch (HB) morphology is a classical problem that can significantly blur our understanding of stellar populations. In this paper, we present a new multivariate analysis connecting the effective temperature extent of the HB with other cluster parameters. The work is based on Hubble Space Telescope photometry of 54 Galactic globular clusters. The present study reveals an important role of the total mass of the globular cluster on its HB morphology. More massive clusters tend to have HBs more extended to higher temperatures. For a set of three input variables including the temperature extension of the HB, [Fe/H] and M_V, the first two eigenvectors account for the 90% of the total sample variance. Possible effects of cluster self-pollution on HB morphology, eventually stronger in more massive clusters, could explain the results here derived.

We present observations of an optically-faint quasar, RD J114816.2+525339, discovered from deep multi-color observations of the field around the z = 6.42 quasar SDSS J1148+5251. The two quasars have a projected separation of 109 arcsec and both are outliers in r-z versus z-J color-color space. Keck spectroscopy reveals RD J114816.2+525339 to be a broad-absorption line quasar at z = 5.70. With z_AB = 23.0, RD J114816.2+525339 is 3.3 mag fainter than SDSS J1148+5251, making it the faintest quasar known at z>5.5. This object was identified in a survey of ~2.5 square degrees. The implied surface density of quasars at these redshifts and luminosities is broadly consistent with previous extrapolations of the faint end of the quasar luminosity function and supports the idea that active galaxies provide only a minor component of the reionizing ultraviolet flux at these redshifts.

Observations of QSOs at z ~ 5.7 - 6.4 show the appearance of Gunn-Peterson troughs around z ~ 6, and a change in the slope of the IGM optical depth tau(z) near z ~ 5.5. These results are interpreted as a signature of the end of the reionization era, which probably started at considerably higher redshifts. However, there also appears to be a substantial cosmic variance in the transmission of the IGM, both along some lines of sight, and among different lines of sight, in this intriguing redshift regime. We suggest that this is indicative of a spatially uneven reionization, possibly caused by the bias-driven primordial clustering of the reionization sources. There is also some independent evidence for a strong clustering of QSOs at z ~ 4 - 5 and galaxies around them, supporting the idea of the strong biasing of the first luminous sources at these redshifts. Larger samples of high-z QSOs are needed in order to provide improved, statistically significant constraints for the models of these phenomena. We expect that the Palomar-Quest (PQ) survey will soon provide a new set of QSOs to be used as cosmological probes in this redshift regime.

With a redshift of z=5.5 and an optical blue magnitude M_B ~ -24.2 mag (~4.5 10^12 L_sun), RDJ030117+002025 is the most distant optically faint (M_B > -26 mag) quasar known. MAMBO continuum observations at lambda=1.2 mm (185 micrometer rest-frame) showed that this quasar has a far-IR luminosity comparable to its optical luminosity. We present near-infrared J- and K-band photometry obtained with NIRC on the Keck I telescope, tracing the slope of the rest frame UV spectrum of this quasar. The observed spectral index is close to the value of alpha_nu ~ -0.44 measured in composite spectra of optically-bright SDSS quasars. It thus appears that the quasar does not suffer from strong dust extinction, which further implies that its low rest-frame UV luminosity is due to an intrinsically-faint AGN. The FIR to optical luminosity ratio is then much larger than that observed for the more luminous quasars, supporting the suggestion that the FIR emission is not powered by the AGN but by a massive starburst.

Color gradients in early-type galaxies contain valuable clues about their formation and evolutionary histories and mechanisms. We examine color gradients in 1,700 early-type galaxies in 159 galaxy clusters spanning a redshift range of 0.05 to 0.2. We find that color gradients strongly depend on the environment where galaxies reside, with steeper color gradients in poor rather than rich clusters. No dependence of color gradients on galaxy luminosity is found in either rich or poor clusters. The difference in color gradients can be explained by a change in the internal metallicity and/or an age gradient in these galaxies. Our results support a reasonable picture whereby young early-type galaxies form in a dissipative collapse process, and then undergo increased (either major or minor) merging activity in richer rather than in poor clusters.

All sciences, including astronomy, are now entering the era of information abundance. The exponentially increasing volume and complexity of modern data sets promises to transform the scientific practice, but also poses a number of common technological challenges. The Virtual Observatory concept is the astronomical community's response to these challenges: it aims to harness the progress in information technology in the service of astronomy, and at the same time provide a valuable testbed for information technology and applied computer science. Challenges broadly fall into two categories: data handling (or "data farming"), including issues such as archives, intelligent storage, databases, interoperability, fast networks, etc., and data mining, data understanding, and knowledge discovery, which include issues such as automated clustering and classification, multivariate correlation searches, pattern recognition, visualization in highly hyperdimensional parameter spaces, etc., as well as various applications of machine learning in these contexts. Such techniques are forming a methodological foundation for science with massive and complex data sets in general, and are likely to have a much broather impact on the modern society, commerce, information economy, security, etc. There is a powerful emerging synergy between the computationally enabled science and the science-driven computing, which will drive the progress in science, scholarship, and many other venues in the 21st century.

In this paper we present an objectively defined catalog of 459 small, high density groups of galaxies out to z~0.2 in a region of 6260 square degrees in the northern sky derived from the Digitized Second Palomar Observatory Sky Survey. Our catalog extends down to r = 19.0 and has a median redshift of z_med = 0.12, making it complementary to Hickson's catalog for the nearby universe (z_med = 0.03). The depth and angular coverage of this catalog makes it valuable for studies of the general characteristics of small groups of galaxies and how galaxies evolve in and around them. We also examine the relationship between compact groups and large scale structure.

We review the origins of the Virtual Observatory (VO) concept, and the current status of the efforts in this field. VO is the response of the astronomical community to the challenges posed by the modern massive and complex data sets. It is a framework in which information technology is harnessed to organize, maintain, and explore the rich information content of the exponentially growing data sets, and to enable a qualitatively new science to be done with them. VO will become a complete, open, distributed, web-based framework for astronomy of the early 21st century. A number of significant efforts worldwide are now striving to convert this vision into reality. The technological and methodological challenges posed by the information-rich astronomy are also common to many other fields. We see a fundamental change in the way all science is done, driven by the information technology revolution.

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

We present tentative evidence for primordial clustering, manifested as an excess of color-selected objects in the field of the QSO SDSS J1030+0524 at redshift z=6.28. We have selected objects red in i_{775}-z_{850} on the basis of Hubble Space Telescope Advanced Camera for Surveys imaging of a field centered on the QSO. Compared to data at comparable depth obtained by the GOODS survey, we find an excess of objects with (i_{775}-z_{850}) \geq 1.5 in the QSO field. The significance of the detection is estimated to be ~97% on the basis of the counts alone and increases to 99.4% if one takes into account the color distribution. If confirmed this would represent the highest redshift example of galaxy clustering and would have implications on models for the growth of structure. Bias-driven clustering of first luminous objects forming in the highest peaks of the primordial density field is expected in most models of early structure formation. The redshift of one of the candidates has been found to be z=5.970 by our spectroscopy with Keck I/LRIS, confirming the validity of our color selection.

We report new HST/WFPC2 photometry for 10 globular clusters (GC) in M31 observed in F5555W(V) and F814W(I). Additionally we have reanalyzed HST archival data of comparable quality for 2 more GCs. Extraordinary care is taken to account for the effects of blended stellar images and required field subtraction. We thus reach 1 mag fainter than the horizontal branch (HB) even in unfavorable cases. We present the color-magnitude diagrams (CMDs) and discuss their main features also in comparison with the properties of the Galactic GCs. This analysis is augmented with CMDs previously obtained and discussed by Fusi Pecci et al. (1996) on 8 other M31 clusters. We report the following significant results: 1. The locus of the red giant branches give reliable metallicity determinations which compare generally very well with ground-based integrated spectroscopic and photometric measures, as well as giving good reddening estimates. 2. The HB morphologies show the same behavior with metallicity as the Galactic GCs, with indications that the 2nd-parameter effect can be present in some GCs of our sample. However, at [Fe/H] ~ -1.7 we observe a number of GCs with red HB morphology such that the HB type versus [Fe/H] relation is offset from the MW and resembles that of the Fornax dwarf spheroidal galaxy. One explanation for the offset is that they are younger than their MW counterparts by 1-2 Gyr. 3. The Mv(HB)-[Fe/H] relationship has been determined and the slope (~0.20) is very similar to the values derived from RR Lyrae stars in the MW and the LMC. The zero-point of this relation based on the assumed distance modulus (m-M)o(M31)=24.47+/-0.03 is consistent with (m-M)o(LMC)=18.55.

The morphology of clusters of galaxies may be described with a set of parameters which contain information about the formation and evolutionary history of these systems. In this paper we present a preliminary study of the morphological parameters of a sample of 28 compact Abell clusters extracted from DPOSS data. The morphology of galaxy clusters is parameterized by their apparent ellipticity, position angle of the major axis, centre coordinates, core radius and beta-model power law index. Our procedure provides estimates of these parameters by simultaneously fitting them all, overcoming some of the difficulties induced by sparse data and low number statistics typical of this kind of analysis. The cluster parameters were fitted in a 3 x 3 h^-2 sqMpc region, measuring the background in a 2 <R< 2.5 h^-1Mpc annulus. We also explore the correlations between shape and profile parameters and other cluster properties. One third of this compact cluster sample has core radii smaller than 50 h^-1 kpc, i.e. near the limit that our data allow us to resolve, possibly consistent with cusped models. The remaining clusters span a broad range of core radii up to 750 h^-1 kpc. More than 80 per cent of this sample has ellipticity higher than 0.2. The alignment between the cluster and the major axis of the dominant galaxy is confirmed, while no correlation is observed with other bright cluster members. No significant correlation is found between cluster richness and ellipticity. Instead, cluster richness is found to correlate, albeit with large scatter, with the cluster core radius.[abridged]

One of the new frontiers of astronomical research is the exploration of time variability on the sky at different wavelengths and flux levels. We have carried out a pilot project using DPOSS data to study strong variables and transients, and are now extending it to the new Palomar-QUEST synoptic sky survey. We report on our early findings and outline the methodology to be implemented in preparation for a real-time transient detection pipeline. In addition to large numbers of known types of highly variable sources (e.g., SNe, CVs, OVV QSOs, etc.), we expect to find numerous transients whose nature may be established by a rapid follow-up. Whereas we will make all detected variables publicly available through the web, we anticipate that email alerts would be issued in the real time for a subset of events deemed to be the most interesting. This real-time process entails many challenges, in an effort to maintain a high completeness while keeping the contamination low. We will utilize distributed Grid services developed by the GRIST project, and implement a variety of advanced statistical and machine learning techniques.

The Grist project (http://grist.caltech.edu/) is developing a grid-technology based system as a research environment for astronomy with massive and complex datasets. This knowledge extraction system will consist of a library of distributed grid services controlled by a workflow system, compliant with standards emerging from the grid computing, web services, and virtual observatory communities. This new technology is being used to find high redshift quasars, study peculiar variable objects, search for transients in real time, and fit SDSS QSO spectra to measure black hole masses. Grist services are also a component of the ``hyperatlas'' project to serve high-resolution multi-wavelength imagery over the Internet. In support of these science and outreach objectives, the Grist framework will provide the enabling fabric to tie together distributed grid services in the areas of data access, federation, mining, subsetting, source extraction, image mosaicking, statistics, and visualization.

The cosmic reionization era, which includes formation of the first stars, galaxies, and AGN, is now one of the most active frontiers of cosmological research. We review briefly our current understanding of the early structure formation, and use the ideas about a joint formation of massive black holes (which power the early QSOs) and their host galaxies to employ high-redshift QSOs as probes of the early galaxy formation and primordial large-scale structure. There is a growing evidence for a strong biasing in the formation of the first luminous sources, which would lead to a clumpy reionization. Absorption spectroscopy of QSOs at z > 6 indicates the end of the reionization era at z ~ 6; yet measurements from the WMAP satellite suggest and early reionization at z ~ 10 - 20. The first generation of massive stars, perhaps aided by the early mini-quasars, may have reionized the universe at such high redshifts, but their feedback may have disrupted the subsequent star and galaxy formation, leading to an extended and perhaps multimodal reionization history ending by z ~ 6. Observations of gamma-ray bursts from the death events of these putative Population III stars may provide essential insight into the primordial structure formation, reionization, early chemical enrichment, and formation of seed black holes which may grow to become central engines of luminous quasars.

Exploration of the time variability on the sky over a broad range of flux levels and wavelengths is rapidly becoming a new frontier of astronomical research. We describe here briefly the Palomar-QUEST survey being carried out from the Samuel Oschin 48-inch Schmidt telescope at Palomar. The following features make the survey an attractive candidate for studying time variability: anticipated survey area of 12,000 - 15,000 sq. degrees in the drift scan mode, point source depth of 21st mag. in I under good conditions, near simultaneous observations in four filters, and at least four passes per year at each location covered. The survey will yield a large number of transients and highly variable sources in the near future and in that sense is a prototype of LSST and Pan-STARRS. We briefly outline our strategy for searching such objects and the proposed pipeline for detecting transients in real-time.

We report the discovery of a Lyman-alpha-emitting galaxy at redshift z=6.545 serendipitously identified in the course of spectroscopic follow-up of hard X-ray sources on behalf of the Serendipitous Extragalactic X-Ray Source Identification (SEXSI) survey. The line flux of the galaxy, 2.1e-17 erg/cm2/s, is similar to line fluxes probed by narrow-band imaging surveys; the 5.2 square-arcminutes surveyed implies a surface density of z~6.5 Lyman-alpha emitters somewhat higher than that inferred from narrow-band surveys. This source marks the sixth Lyman-alpha-emitting galaxy identified at z~6.5, a redshift putatively beyond the epoch of reionization when the damping wings of the neutral hydrogen of the intergalactic medium is capable of severely attenuating Lyman-alpha emission. By comparing the Lyman-alpha emitter luminosity functions at z~5.7 and z~6.5, we infer that the intergalactic medium may remain largely reionized from the local universe out to z~6.5.

We present an optically selected galaxy cluster catalog from ~ 2,700 square degrees of the Digitized Second Palomar Observatory Sky Survey (DPOSS), spanning the redshift range 0.1 < z < 0.5, providing an intermediate redshift supplement to the previous DPOSS cluster survey. This new catalog contains 9,956 cluster candidates and is the largest resource of rich clusters in this redshift range to date. The candidates are detected using the best DPOSS plates based on seeing and limiting magnitude. The search is further restricted to high galactic latitude (|b| > 50), where stellar contamination is modest and nearly uniform. We also present a performance comparison of two different detection methods applied to this data, the Adaptive Kernel and Voronoi Tessellation techniques. In the regime where both catalogs are expected to be complete, we find excellent agreement, as well as with the most recent surveys in the literature. Extensive simulations are performed and applied to the two different methods, indicating a contamination rate of ~ 5%. These simulations are also used to optimize the algorithms and evaluate the selection function for the final cluster catalog. Redshift and richness estimates are also provided, making possible the selection of subsamples for future studies.

We describe a method to derive the expansion and acceleration rates directly from the data, without the need for the specification of a theory of gravity, and without adopting an a priori parameterization of the form or redshift evolution of the dark energy. If one also specifies a theory of gravity we can also determine the pressure, energy density, and equation of state of the dark energy as functions of redshift. We then apply this methodology on a modern data set of distances to Supernovae and to radio galaxies. We find that the universe transitions from deceleration to acceleration at a redshift of about 0.4, and the present value of deceleration parameter is -0.35 +/- 0.15. The standard ``concordance model'' provides a reasonably good fit to the dimensionless expansion rate as a function of redshift, though it fits the dimensionless acceleration rate as a function of redshift less well. Adopting General Relativity as the theory of gravity, we obtain the redshift trends for the pressure, energy density, and equation of state of the dark energy out to a redshift of about one. They are generally consistent with the concordance model, at least out to a redshift of about 0.5, but the existing data preclude any stronger conclusions at this point.

Understanding of the nature of dark energy, which appears to drive the expansion of the universe, is one of the central problems of physical cosmology today. In an earlier paper [Daly & Djorgovski (2003)] we proposed a novel method to determine the expansion rate and the deceleration parameter in a largely model-independent way, directly from the data on coordinate distances. Here we expand this methodology to include measurements of the pressure of dark energy, its normalized energy density, and the equation of state parameter as functions of redshift. We then apply this methodology to a new, combined data set of distances to supernovae and radio galaxies.

One of the goals of current cosmological studies is the determination of the expansion and acceleration rates of the universe as functions of redshift, and the determination of the properties of the dark energy that can explain these observations. Here the expansion and acceleration rates are determined directly from the data, without the need for the specification of a theory of gravity, and without adopting an a priori parameterization of the form or redshift evolution of the dark energy. We use the latest set of distances to SN standard candles from Riess et al. (2004), supplemented by data on radio galaxy standard ruler sizes, as described by Daly and Djorgovski (2003, 2004). We find that the universe transitions from acceleration to deceleration at a redshift of about 0.4. The standard "concordance model" provides a reasonably good fit to the dimensionless expansion rate as a function of redshift, though it fits the dimensionless acceleration rate as a function of redshift less well. The expansion and acceleration rates are then combined with a theory of gravity to determine the pressure, energy density, and equation of state of the dark energy as functions of redshift. Adopting General Relativity as the correct theory of gravity, the redshift trends for the pressure, energy density, and equation of state of the dark energy out to redshifts of about one are determined, and are found to be generally consistent with the concordance model.

The advent of wide-area multicolour synoptic sky surveys is leading to data sets unprecedented in size, complexity and data throughput. VO technology offers a way to exploit these to the full but requires changes in design philosophy. The Palomar-QUEST survey is a major new survey being undertaken by Caltech, Yale, JPL and Indiana University to repeatedly observe 1/3 of the sky (~15000 sq. deg. between -27 < Dec <27 in seven passbands. Utilising the 48-inch Oschin Schmidt Telescope at the Palomar Observatory with the 112-CCD QUEST camera covering the full 4 x 4 sq. deg. field of view, it will generate \~1TB of data per month. In this paper, we review the design of QUEST as a VO resource, a federated data set and an exemplar of VO standards.

We discuss the main properties of the Galactic globular cluster (GC) blue straggler stars (BSS), as inferred from our new catalog containing nearly 3000 BSS. The catalog has been extracted from the photometrically homogeneous V vs. (B-V) color-magnitude diagrams (CMD) of 56 GCs, based on WFPC2 images of their central cores. In our analysis we used consistent relative distances based on the same photometry and calibration. The number of BSS has been normalized to obtain relative frequencies (F_{BSS}) and specific densities (N_S) using different stellar populations extracted from the CMD. The cluster F_{BSS} is significantly smaller than the relative frequency of field BSS. We find a significant anti-correlation between the BSS relative frequency in a cluster and its total absolute luminosity (mass). There is no statistically significant trend between the BSS frequency and the expected collision rate. F_{BSS} does not depend on other cluster parameters, apart from a mild dependence on the central density. PCC clusters act like normal clusters as far as the BSS frequency is concerned. We also show that the BSS luminosity function for the most luminous clusters is significantly different, with a brighter peak and extending to brighter luminosities than in the less luminous clusters. These results imply that the efficiency of BSS production mechanisms and their relative importance vary with the cluster mass.

The Atlasmaker project is using Grid technology, in combination with NVO interoperability, to create new knowledge resources in astronomy. The product is a multi-faceted, multi-dimensional, scientifically trusted image atlas of the sky, made by federating many different surveys at different wavelengths, times, resolutions, polarizations, etc. The Atlasmaker software does resampling and mosaicking of image collections, and is well-suited to operate with the Hyperatlas standard. Requests can be satisfied via on-demand computations or by accessing a data cache. Computed data is stored in a distributed virtual file system, such as the Storage Resource Broker (SRB). We expect these atlases to be a new and powerful paradigm for knowledge extraction in astronomy, as well as a magnificent way to build educational resources. The system is being incorporated into the data analysis pipeline of the Palomar-Quest synoptic survey, and is being used to generate all-sky atlases from the 2MASS, SDSS, and DPOSS surveys for joint object detection.

Hyperatlas is an open standard intended to facilitate the large-scale federation of image-based data. The subject of hyperatlas is the space of sphere-to-plane projection mappings (the FITS-WCS information), and the standard consists of coherent collections of these on which data can be resampled and thereby federated with other image data. We hope for a distributed effort that will produce a multi-faceted image atlas of the sky, made by federating many different surveys at different wavelengths and different times. We expect that hyperatlas-compliant imagery will be published and discovered through an International Virtual Observatory Alliance (IVOA) registry, and that grid-based services will emerge for the required resampling and mosaicking.

Many past studies of cosmological gamma-ray bursts (GRBs) have been limited because of the large distance to typical GRBs, resulting in faint afterglows. There has long been a recognition that a nearby GRB would shed light on the origin of these mysterious cosmic explosions, as well as the physics of their fireballs. However, GRBs nearer than z=0.2 are extremely rare, with an estimated rate of localisation of one every decade. Here, we report the discovery of bright optical afterglow emission from GRB 030329. Our prompt dissemination and the brilliance of the afterglow resulted in extensive followup (more than 65 telescopes) from radio through X-ray bands, as well as measurement of the redshift, z=0.169. The gamma-ray and afterglow properties of GRB 030329 are similar to those of cosmological GRBs (after accounting for the small distance), making this the nearest known cosmological GRB. Observations have already securely identified the progenitor as a massive star that exploded as a supernova, and we anticipate futher revelations of the GRB phenomenon from studies of this source.

Given the paucity of empirical constraints, the nature of the newly-recognized phenomena called X-Ray Flashes (XRFs) has been an open question. However, with the recent detections of radio and X-ray afterglow it is finally possible to study the large- and small-scale environments of XRFs. We present Chandra, Hubble Space Telescope (HST), and Keck observations of the fields of XRFs 011030 and 020427. Astrometric comparisons of the X-ray transient positions and the HST images reveal the XRFs to be associated with faint blue galaxies. Photometric evidence of these putative hosts suggests that these two XRFs originated from redshifts less than z ~ 3.5, and thus cannot be due to GRBs at very high redshifts. In both host-burst offsets and host properties, these XRFs could have been drawn from distributions similar to those measured of gamma-ray bursts (GRBs). We conclude with a discussion of the implications of this XRF-GRB host connection for the possible progenitors of XRFs.

Determination of the expansion and acceleration history of the universe is one of the fundamental goals of cosmology. Detailed measurements of these rates as a function of redshift can provide new physical insights into the nature and evolution of the dark energy, which apparently dominates the global dynamics of the universe at the present epoch. We present here dimensionless coordinate distances y(z) to twenty radio galaxies reaching out to z of about 1.8, the redshift range currently not covered by Supernova standard candle observations. We develop a simple numerical method for a direct determination of the expansion and acceleration rates, E(z) and q(z), from the data, which makes no assumptions about the underlying cosmological model or the equation of state parameter w. This differential method is in contrast to the traditional cosmological tests, where particular model equations are integrated and then compared with the observations. The new model-independent method is applied to currently available Supernova data and the data on radio galaxies presented here. We derive the expansion rate of the universe as a function of redshift, E(z), and for the first time obtain a direct estimate of the acceleration rate of the universe as a function of redshift, q(z), in a way that is independent of assumptions regarding the dark energy and its redshift evolution. The current observations indicate that the universe transitions from acceleration to deceleration at a redshift greater than 0.3, with a best fit estimate of about 0.45; this transition redshift and our determinations of E(z) are broadly in agreement with the currently popular Friedmann-Lemaitre cosmology with Omega = 0.3, and Lambda = 0.7, even though no model assumptions are made in deriving the fits for E(z) and q(z).

We report a discovery of a quasar at z = 4.96 +- 0.03 within a few Mpc of the quasar SDSS 0338+0021 at z = 5.02 +- 0.02. The newly found quasar has the SDSS i and z magnitudes of ~ 21.2, and an estimated absolute magnitude M_B ~ -25.2. The projected separation on the sky is 196 arcsec, and the redshift difference Delta z = 0.063 +- 0.008. The probability of finding this quasar pair by chance in the absence of clustering in this particular volume is ~ 10^-4 to 10^-3. We conclude that the two objects probably mark a large-scale structure, possibly a protocluster, at z ~ 5. This is the most distant such structure currently known. Our search in the field of 13 other QSOs at z >~ 4.8 so far has not resulted in any detections of comparable luminous QSO pairs, and it is thus not yet clear how representative is this structure at z ~ 5. However, along with the other evidence for clustering of quasars and young galaxies at somewhat lower redshifts, the observations are at least qualitatively consistent with a strong biasing of the first luminous and massive objects, in agreement with general predictions of theoretical models. More extensive searches for clustered quasars and luminous galaxies at these redshifts will provide valuable empirical constraints for our understanding of early galaxy and structure formation.

We present accurate metallicity measurements for 121 damped Lya systems at 0.5<z<5 including ~50 new measurements from our recently published Echellette Spectrograph and Imager surveys. This dataset is analysed to determine the age-metallicity relation of neutral gas in the universe. Contrary to previous datasets this sample shows statistically significant evolution in the mean metallicity. The best linear fit rate to metallicity vs. redshift is -0.26 +/- 0.07 dex corresponding to approximately a factor of 2 every Gyr at z=3. The DLA continue to maintain a floor in metallicity of ~1/700 solar independent of observational effects. This metallicity threshold limits the prevalence of primordial gas in high redshift galaxies and stresses the correspondence between damped systems and star formation (i.e. galaxy formation). This floor is significantly offset from the metallicity of the Lya forest and therefore we consider it to be more related to active star formation within these galaxies than scenarios of enrichment in the very early universe. Finally, we comment on an apparent 'missing metals problem': the mean metallicity of the damped systems is ~10x lower than the value expected from their observed star formation history. This problem is evident in current theoretical treatments of chemical evolution and galaxy formation; it may indicate a serious flaw in our understanding of the interplay between star formation and metal production.

This paper presents chemical abundance measurements for 12 new z>3 damped Lya systems discovered toward five quasars from the Palomar Sky Survey. We determine HI column densities from profile fits to the observed damped Lya profiles and measure ionic column densities and limits for all observed metal-line transitions. This dataset, acquired with the Echellette Spectrograph and Imager on the KeckII telescope, adds to the rapidly growing database of damped Lya abundances. It will impact studies of chemical evolution in the early universe and help identify candidates for detailed follow-up observations with echelle spectrographs. We report the discovery of the first quasar sightline with four cosmologically distinct damped Lya systems.

With the recent release of large (i.e., > hundred million objects), well-calibrated photometric surveys, such as DPOSS, 2MASS, and SDSS, spectroscopic identification of important targets is no longer a simple issue. In order to enhance the returns from a spectroscopic survey, candidate sources are often preferentially selected to be of interest, such as brown dwarfs or high redshift quasars. This approach, while useful for targeted projects, risks missing new or unusual species. We have, as a result, taken the alternative path of spectroscopically identifying interesting sources with the sole criterion being that they are in low density areas of the g - r and r - i color-space defined by the DPOSS survey. In this paper, we present three peculiar broad absorption line quasars that were discovered during this spectroscopic survey, demonstrating the efficacy of this approach. PSS J0052+2405 is an Iron LoBAL quasar at a redshift z = 2.4512 with very broad absorption from many species. PSS J0141+3334 is a reddened LoBAL quasar at z = 3.005 with no obvious emission lines. PSS J1537+1227 is a Iron LoBAL at a redshift of z = 1.212 with strong narrow Mgii and Feii emission. Follow-up high resolution spectroscopy of these three quasars promises to improve our understanding of BAL quasars. The sensitivity of particular parameter spaces, in this case a two-color space, to the redshift of these three sources is dramatic, raising questions about traditional techniques of defining quasar populations for statistical analysis.

Images of the CO 2-1 line emission, and the radio continuum emission, from the redshift 4.12 gravitationally lensed quasi-stellar object (QSO) PSS J2322+1944 reveal an Einstein ring with a diameter of 1.5". These observations are modeled as a star forming disk surrounding the QSO nucleus with a radius of 2 kpc. The implied massive star formation rate is 900 M_sun/year. At this rate a substantial fraction of the stars in a large elliptical galaxy could form on a dynamical time scale of 10^8 years. The observation of active star formation in the host galaxy of a high-redshift QSO supports the hypothesis of coeval formation of supermassive black holes and stars in spheroidal galaxies.

Correlations among the independently measured physical properties of globular clusters (GCs) can provide powerful tests for theoretical models and new insights into their dynamics, formation, and evolution. We review briefly some of the previous work, and present preliminary results from a comparative study of GC correlations in the Local Group galaxies. The results so far indicate that these diverse GC systems follow the same fundamental correlations, suggesting a commonality of formative and evolutionary processes which produce them.

We present the discovery of GRB 020405 made with the Inter-Planetary Network (IPN). With a duration of 60 s, the burst appears to be a typical long duration event. We observed the 75-square acrminute IPN error region with the Mount Stromlo Observatory's 50-inch robotic telescope and discovered a transient source which subsequently decayed and was also associated with a variable radio source. We identify this source as the afterglow of GRB 020405. Subsequent observations by other groups found varying polarized flux and established a redshift of 0.690 to the host galaxy. Motivated by the low redshift we triggered observations with WFPC2 on-board the Hubble Space Telescope (HST). Modeling the early ground-based data with a jet model, we find a clear red excess over the decaying optical lightcurves that is present between day 10 and day 141 (the last HST epoch). This `bump' has the spectral and temporal features expected of an underlying supernova (SN). In particular, the red color of the putative SN is similar to that of the SN associated with GRB 011121, at late time. Restricting the sample of GRBs to those with z<0.7, a total of five bursts, red bumps at late times are found in GRB 970228, GRB 011121, and GRB 020405. It is possible that the simplest idea, namely that all long duration GRBs have underlying SNe with a modest dispersion in their properties (especially peak luminosity), is sufficient to explain the non detections.

We have identified eighty-four small, high density groups of galaxies out to z ~ 0.2 in a region of ~ 2000 square degrees around the north galactic pole using DPOSS (the Digitized Second Palomar Observatory Sky Survey). The groups have at least four galaxies satisfying more stringent criteria than those used by Hickson in his pioneering work in 1982: the adopted limiting surface brightness for each group is brighter (24 mag/arcsec^2 instead of 26 mag/arcsec^2), and the spread in magnitude among the member galaxies is narrower (two magnitudes instead of three). We also adopt a slightly modified version of the isolation criterion used by Hickson, in order to avoid rejecting groups with projected nearby faint background galaxies. A 10% contamination rate due to projection effects is expected for this sample based on extensive simulations.

This is a review article for the book "Compact Stellar X-Ray Sources", Editors W. Lewin and M. van der Klis, to be published by Cambridge University Press in 2003. It covers the phenomenology of cosmic gamma-ray bursts, the theory of their afterglows, and the observations and interpretation of their host galaxies.

Studies of the cosmic gamma-ray bursts (GRBs) and their host galaxies are starting to provide interesting or even unique new insights in observational cosmology. GRBs represent a new way of identifying a population of star-forming galaxies at cosmological redshifts. GRB hosts are broadly similar to the normal field galaxy populations at comparable redshifts and magnitudes, and indicate at most a mild luminosity evolution out to z ~ 1.5 - 2. GRB optical afterglows seen in absorption provide a powerful new probe of the ISM in dense, central regions of their host galaxies, complementary to the traditional studies using QSO absorbers. Some GRB hosts are heavily obscured, and provide a new way to select a population of cosmological sub-mm sources, and a novel constraint on the total obscured fraction of star formation over the history of the universe. Finally, detection of GRB afterglows at z > 6 may provide a unique way to probe the primordial star formation, massive IMF, early IGM, and chemical enrichment at the end of the cosmic reionization era.

Studies of the cosmic gamma-ray bursts (GRBs) and their host galaxies are now starting to provide interesting or even unique new insights in observational cosmology. Observed GRB host galaxies have a median magnitude R ~ 25 mag, and show a range of luminosities, morphologies, and star formation rates, with a median redshift z ~ 1. They represent a new way of identifying a population of star-forming galaxies at cosmological redshifts, which is mostly independent of the traditional selection methods. They seem to be broadly similar to the normal field galaxy populations at comparable redshifts and magnitudes, and indicate at most a mild luminosity evolution over the redshift range they probe. Studies of GRB optical afterglows seen in absorption provide a powerful new probe of the ISM in dense, central regions of their host galaxies, which is complementary to the traditional studies using QSO absorption line systems. Some GRB hosts are heavily obscured, and provide a new way to select a population of cosmological sub-mm sources. A census of detected optical tranistents may provide an important new way to constrain the total obscured fraction of star formation over the history of the universe. Finally, detection of GRB afterglows at high redshifts (z > 6) may provide a unique way to probe the primordial star formation, massive IMF, early IGM, and chemical enrichment at the end of the cosmic reionization era.

We present a new, objectively defined catalog of candidate galaxy clusters based on the galaxy catalogs from the Digitized Second Palomar Observatory Sky Survey (DPOSS). This cluster catalog, derived from the best calibrated plates in the high latitude (|b|>30) Northern Galactic Cap region, covers 5,800 square degrees, and contains 8,155 candidate clusters. A simple adaptive kernel density mapping technique, combined with the SExtractor object detection algorithm, is used to detect galaxy overdensities, which we identify as clusters. Simulations of the background galaxy distribution and clusters of varying richnesses and redshifts allow us to optimize detection parameters, and measure the completeness and contamination rates for our catalog. Cluster richnesses and photometric redshifts are measured, using integrated colors and magnitudes for each cluster. An extensive spectroscopic survey is used to confirm the photometric results. This catalog, with well-characterized sample properties, provides a sound basis for future studies of cluster physics and large scale structure.

We report on Hubble Space Telescope WFPC2 optical and STIS near ultraviolet MAMA observations, and ground-based optical observations of GRB010222, spanning 15 hrs to 71 days. The observations are well-described by a relativistic blast-wave model with a hard electron-energy distribution, p = 1.57, and a jet transition at t_j=0.93 days. These values are slightly larger than previously found as a result of a correction for the contribution from the host galaxy to the late-time ground-based observations and the larger temporal baseline provided by the Hubble Space Telescope observations. The host galaxy is found to contain a very compact core (size <0.25 arcsec) which coincides with the position of the optical transient. The STIS near ultraviolet MAMA observations allow for an investigation of the extinction properties along the line of sight to GRB010222. We find that the far ultraviolet curvature component (c_4) is rather large. In combination with the low optical extinction A_V =0.11 mag, when compared to the Hydrogen column inferred from X-ray observations, we suggest that this is evidence for dust destruction.

Using the Echellete Spectrograph Imager (ESI) on the Keck II 10-m telescope we have measured the redshifts of the host galaxies of gamma-ray bursts GRB 990506 and GRB 000418, z=1.30658 +/- 0.00004 and 1.1181 +/- 0.0001, respectively. Thanks to the excellent spectral resolution of ESI we resolved the [O II] 3727 doublet in both cases. The measured redshift of GRB 990506 is the highest known for a dark burst GRB, though entirely consistent with the notion that dark and non-dark bursts have a common progenitor origin. The relative strengths of the [O II], He I, [Ne III], and H gamma emission lines suggest that the host of GRB 000418 is a starburst galaxy, rather than a LINER or Seyfert 2. Since the host of GRB 000418 has been detected at sub-millimeter wavelengths these spectroscopic observations suggest that the sub-millimeter emission is due to star-formation (as opposed to AGN) activity. The [O II]-derived unobscured star-formation rates are 13 and 55 M_solar/yr for the hosts of GRB 990506 and GRB 000418, respectively. In contrast, the star-formation rate of the host of GRB 000418 derived from sub-millimeter observations is twenty times larger.

A Topic Map is a structured network of hyperlinks that points into an information pool. Topic Maps have an existence independent of the information pool and hence different Topic Maps can form different layers above the same information pool and provide us with different views of it. We explore the use of Topic Maps with the Unified Column Descriptor (UCD) scheme developed in the frame of the ESO-CDS data mining project. UCD, with its multi-tier hierarchical structure, categorizes parameters reported in tables and catalogs. By using Topic Maps we show how columns from different catalogs with similar but not identical descriptions could be combined. A direct application for the Virtual Observatory community is that of merging catalogs in order to generate customized views of data.

We present revised infrared (2.2 micron) astrometry of the field containing the Galactic microquasar GRS 1758-258, using observations at the Keck I 10-m telescope. We find three candidates for the microquasar within a 3 sigma error circle, but none within 2 sigma. We show that if the 18.4 day X-ray period of GRS 1758-258 is due to a binary orbit, then only one of the three candidates, an early K-type giant, is large enough to power the microquasar via Roche lobe overflow. We therefore identify this star as the infrared counterpart of GRS 1758-258, which we classify as a low mass X-ray binary. Long term infrared monitoring of this source should provide further information about the microquasar system, including a confirmation of the X-ray period and an estimate of the compact object's mass.

We present the photometric calibration technique for the Digitized Second Palomar Observatory Sky Survey (DPOSS), used to create seamless catalogs of calibrated objects over large sky areas. After applying a correction for telescope vignetting, the extensive plate overlap regions are used to transform sets of plates onto a common instrumental photometric system. Photometric transformations to the Gunn gri system for each plate, for stars and galaxies, are derived using these contiguous stitched areas and an extensive CCD imaging library obtained for this purpose. We discuss the resulting photometric accuracy, survey depth, and possible systematic errors.

There has been an unprecedented and continuing growth in the volume, quality, and complexity of astronomical data sets over the past few years, mainly through large digital sky surveys. Virtual Observatory (VO) concept represents a scientific and technological framework needed to cope with this data flood. We review some of the applied statistics and computing challenges posed by the analysis of large and complex data sets expected in the VO-based research. The challenges are driven both by the size and the complexity of the data sets (billions of data vectors in parameter spaces of tens or hundreds of dimensions), by the heterogeneity of the data and measurement errors, the selection effects and censored data, and by the intrinsic clustering properties (functional form, topology) of the data distribution in the parameter space of observed attributes. Examples of scientific questions one may wish to address include: objective determination of the numbers of object classes present in the data, and the membership probabilities for each source; searches for unusual, rare, or even new types of objects and phenomena; discovery of physically interesting multivariate correlations which may be present in some of the clusters; etc.

The advent of the Virtual Observatory (VO) concept signals a pradigm shift in the way astronomy will be done in the era of information abundance and ubiquitous networking. Small telescopes will be playing a number of essential roles in this new research environment, probably contributing a major portion of all data taken in astronomy, both as surveying instruments, and as follow-up facilities. In this review we describe the VO concept and its background, and how small telescopes will fit in this emergent way of doing astronomy in the 21st century.

We present ground-based optical observations of GRB 020124 starting 1.6 hours after the burst, as well as subsequent Very Large Array (VLA) and Hubble Space Telescope (HST) observations. The optical afterglow of GRB 020124 is one of the faintest afterglows detected to date, and it exhibits a relatively rapid decay, $F_\nu\propto t^{-1.60\pm 0.04}$, followed by further steepening. In addition, a weak radio source was found coincident with the optical afterglow. The HST observations reveal that a positionally coincident host galaxy must be the faintest host to date, R>29.5 mag. The afterglow observations can be explained by several models requiring little or no extinction within the host galaxy, A_V~0-0.9 mag. These observations have significant implications for the statistics of the so-called dark bursts (bursts for which no optical afterglow is detected), which are usually attributed to dust extinction within the host galaxy. The faintness and relatively rapid decay of the afterglow of GRB 020124, combined with the low inferred extinction indicate that some dark bursts are intrinsically dim and not dust obscured. Thus, the diversity in the underlying properties of optical afterglows must be observationally determined before substantive inferences can be drawn from the statistics of dark bursts.

GRB 010921 was the first HETE-2 GRB to be localized via its afterglow emission. The low-redshift of the host galaxy, z=0.451, prompted us to undertake intensive multi-color observations with the Hubble Space Telescope with the goal of searching for an underlying supernova component. We do not detect any coincident supernova to a limit 1.34 mag fainter than SN 1998bw at 99.7% confidence, making this one of the most sensitive searches for an underlying SN. Analysis of the afterglow data allow us to infer that the GRB was situated behind a net extinction (Milky Way and the host galaxy) of A_V ~ 1.8 mag in the observer frame. Thus, had it not been for such heavy extinction our data would have allowed us to probe for an underlying SN with brightness approaching those of more typical Type Ib/c supernovae.

We present the complete photometric database and the color-magnitude diagrams for 74 Galactic globular clusters observed with the HST/WFPC2 camera in the F439W and F555W bands. A detailed discussion of the various reduction steps is also presented, and of the procedures to transform instrumental magnitudes into both the HST F439W and F555W flight system and the standard Johnson B and V systems. We also describe the artificial star experiments which have been performed to derive the star count completeness in all the relevant branches of the color magnitude diagram. The entire photometric database and the completeness function will be made available on the Web immediately after the publication of the present paper.

We present sequential optical spectra of the afterglow of GRB 010222 obtained one day apart using the Low Resolution Imaging Spectrometer (LRIS) and the Echellette Spectrograph and Imager (ESI) on the Keck telescopes. Three low-ionization absorption systems are spectroscopically identified at z1=1.47688, z2=1.15628, and z3=0.92747. The higher resolution ESI spectrum reveals two distinct components in the highest redshift system at z1a=1.47590 and z1b=1.47688. We interpret the z1b=1.47688 system as an absorption feature of the disk of the host galaxy of GRB 010222. The best fitted power-law optical continuum and [Zn/Cr] ratio imply low dust content or a local gray dust component near the burst site. In addition, we do not detect strong signatures of vibrationally excited states of H_2. If the GRB took place in a superbubble or young stellar cluster, there are no outstanding signatures of an ionized absorber, either. Analysis of the spectral time dependence at low resolution shows no significant evidence for absorption-line variability. This lack of variability is confronted with time-dependent photoionization simulations designed to apply the observed flux from GRB 010222 to a variety of assumed atomic gas densities and cloud radii. The absence of time dependence in the absorption lines implies that high-density environments are disfavored. In particular, if the GRB environment was dust free, its density was unlikely to exceed nH=10^2 cm^-3. If depletion of metals onto dust is similar to Galactic values or less than solar abundances are present, then nH > 2 x 10^4 cm^-3 is probably ruled out in the immediate vicinity of the burst.

We present radio observations of the afterglow of the bright gamma-ray burst GRB980329 made between one month and several years after the burst, a re-analysis of previously published submillimeter data, and late-time optical and near-infrared (NIR) observations of the host galaxy. From the absence of a spectral break in the optical/NIR colors of the host galaxy, we exclude the earlier suggestion that GRB980329 lies at a redshift of z >~5. We combine our data with the numerous multi-wavelength observations of the early afterglow, fit a comprehensive afterglow model to the entire broadband dataset, and derive fundamental physical parameters of the blast-wave and its host environment. Models for which the ejecta expand isotropically require both a high circumburst density and extreme radiative losses from the shock. No low density model (n << 10 cm^{-3}) fits the data. A burst with a total energy of ~ 10^{51} erg, with the ejecta narrowly collimated to an opening angle of a few degrees, driven into a surrounding medium with density ~ 20 cm^{-3}, provides a satisfactory fit to the lightcurves over a range of redshifts.

Using observations from an extensive monitoring campaign with the Hubble Space Telescope we present the detection of an intermediate-time flux excess which is redder in color relative to the afterglow of GRB 011121, currently distinguished as the gamma-ray burst with the lowest known redshift. The red ``bump,'' which exhibits a spectral roll-over at ~7200 Angstrom, is well described by a redshifted Type Ic supernova that occurred approximately at the same time as the gamma-ray burst event. The inferred luminosity is about half that of the bright supernova 1998bw. These results serve as compelling evidence for a massive star origin of long-duration gamma-ray bursts. Models that posit a supernova explosion weeks to months preceding the gamma-ray burst event are excluded by these observations. Finally, we discuss the relationship between spherical core-collapse supernovae and gamma-ray bursts.

Of the cosmological gamma-ray bursts, GRB 011121 has the lowest redshift, z=0.36. More importantly, the multi-color excess in the afterglow detected in the Hubble Space Telescope (HST) light curves is compelling observational evidence for an underlying supernova. Here we present near-infrared and radio observations of the afterglow. We undertake a comprehensive modeling of these observations and those reported in the literature and find good evidence favoring a wind-fed circumburst medium. In detail, we infer the progenitor had a mass loss rate of Mdot ~ 10^-7 / v_w3 Mo/yr where v_w3 is the speed of the wind from the progenitor in units of 10^3 km/s. This mass loss rate is similar to that inferred for the progenitor of SN 1998bw which has been associated with GRB 980425. Our data, taken in conjunction with the HST results of Bloom et al. (2002), provide a consistent picture: the long duration GRB 011121 had a massive star progenitor which exploded as a supernova at about the same time as the GRB event.

We have measured Keck/HIRES radial velocities for 30 candidate red giants in the direction of Palomar 13: an object traditionally cataloged as a compact, low-luminosity globular cluster. From a sample of 21 confirmed members, we find a systemic velocity of 24.1 km/s and a projected, intrinsic velocity dispersion of 2.2 km/s. Although small, this dispersion is several times larger than that expected for a globular cluster of this luminosity and central concentration. Taken at face value, this dispersion implies a mass-to-light ratio of ~ 40 (in solar units) based on the best-fit King-Michie model. The surface density profile of Palomar 13 also appears to be anomalous among Galactic globular clusters -- depending upon the details of background subtraction and model-fitting, Palomar 13 either contains a substantial population of "extra-tidal" stars, or it is far more spatially extended than previously suspected. The full surface density profile is equally well-fit by a King-Michie model having a high concentration and large tidal radius, or by a NFW model. We examine -- and tentatively reject -- a number of possible explanations for the observed characteristics of Palomar 13 (e.g., velocity "jitter" among the red giants, spectroscopic binary stars, non-standard mass functions, modified Newtonian dynamics), and conclude that the two most plausible scenarios are either catastrophic heating during a recent perigalacticon passage, or the presence of a massive dark halo. Thus, the available evidence suggests that Palomar 13 is either a globular cluster which is now in the process of dissolving into the Galactic halo, or a faint, dark-matter-dominated stellar system (ABRIDGED).

Using the IRAM interferometer we have detected J=4-3 and 5-4 CO line emission toward the radio quiet quasar PSS2322+1944. At a redshift of z_CO=4.1199 this is the fourth and strongest detection of CO at z>4. The velocity-integrated CO J=4-3 and J=5-4 line fluxes are 4.21+-0.40 and 3.74+-0.56 Jy km/s, and the linewidth is ~300 km/s. The CO J=10-9 line was searched for but not detected with an upper intensity limit of 30 mJy. The 1.35 mm (250 mum rest wavelength) continuum flux density is 7.5+-1.3 mJy, in agreement with previous bolometer measurements at 1.2 mm with the 30-m IRAM telescope. The 3mm (580 mum rest wavelength) continuum is not detected with a 3 sigma upper limit of 0.7mJy. We also report observations of the 450 mum continuum in PSS2322+1944 using the SCUBA array at the JCMT. The quasar was detected with a 450 mum flux density of 79+-19mJy. At the angular resolution of 4.8" x 2.1" at 1.3mm and 6.2" x 4.9" at 3.2 mm, the interferometer observations do not show evidence of spatial extension in the continuum or CO line emission. Assuming no gravitational magnification, we estimate a molecular gas mass of ~2.5x10^11 M_sun. The molecular gas is warm (40-100 K) and dense (10^3.5-10^4.1 cm^-3). The infrared-to-CO luminosity ratio is ~185 L_sun/(K km/s pc^2), comparable to the values found for ultraluminous infrared galaxies. The detection of CO emission in this high redshift QSO provides further evidence that the radio emission and the millimeter to submillimeter continuum emission are predominantly powered by a starburst which is coeval with the AGN activity.

We report the discovery of the optical and radio afterglow of GRB 010921, the first gamma-ray burst afterglow to be found from a localization by the High Energy Transient Explorer (HETE) satellite. We present optical spectroscopy of the host galaxy which we find to be a dusty and apparently normal star-forming galaxy at z = 0.451. The unusually steep optical spectral slope of the afterglow can be explained by heavy extinction, A_V > 0.5 mag, along the line of sight to the GRB. Dust with similar A_V for the the host galaxy as a whole appears to be required by the measurement of a Balmer decrement in the spectrum of the host galaxy. Thanks to the low redshift, continued observations of the afterglow will enable the strongest constraints, to date, on the existence of a possible underlying supernova.

We present a comprehensive study to measure the locations of gamma-ray bursts (GRBs) relative to their host galaxies. In total, we find the offsets of 20 long-duration GRBs from their apparent host galaxy centers utilizing ground-based images from Palomar and Keck and space-based images from the Hubble Space Telescope (HST). The median projected angular (physical) offset is 0.17 arcsec (1.3 kpc). The median offset normalized by the individual host half-light radii is 0.98 suggesting a strong connection of GRB locations with the UV light of their hosts. This provides strong observational evidence for the connection of GRBs to star formation. We compare the observed offset distribution with an exponential disk, a model for the location of collapsars. The statistical comparison shows good agreement, with the Kolmogorov-Smirnov probability that the observed offsets derive from the model distribution of P_KS = 0.45. We also compare the observed GRB offsets with the expected offset distribution of delayed merging remnant progenitors (black hole--neutron star and neutron star--neutron star binaries). We find that delayed merging remnant progenitors, insofar as the predicted offset distributions from population synthesis studies are representative, can be ruled out at the 2 x 10^{-3} level. This is arguably the strongest observational constraint yet against delayed merging remnants as the progenitors of long-duration GRBs. In the course of this study, we have also discovered the putative host galaxies of GRB 990510 and GRB 990308 in archival HST data. (abstract abridged)

We present early-time Keck spectroscopic observations and late-time Hubble Space Telescope (HST) imaging of GRB 000926. The HST images show a small offset between the optical transient and the compact host galaxy. Combined with the large equivalent widths measured for metallic absorption lines by the Keck Echellette Spectrograph and Imager (ESI) and the Low-Resolution Imaging Spectrometer (LRIS), this indicates that the GRB exploded near the center of its host. The ESI spectroscopy reveals two absorption systems centered at z = 2.0379 \pm 0.0008 with a velocity separation of 168 km s^{-1}, which we interpret as being due to individual clouds in the host galaxy. The ratios of chromium to zinc equivalent widths indicates the host is depleted in dust relative to local values to a similar degree as damped Lyman alpha systems at the same redshift. Further, the two clouds appear to have similar relative metal abundance and dust to gas ratio. If one cloud is associated with the GRB site, this implies the explosion did not significantly alter the surrounding environment.

Of all the well localized gamma-ray bursts, GRB 000911 has the longest duration (T_90 ~ 500 s), and ranks in the top 1% of BATSE bursts for fluence. Here, we report the discovery of the afterglow of this unique burst. In order to simultaneously fit our radio and optical observations, we are required to invoke a model involving an hard electron distribution, p ~ 1.5 and a jet-break time less than 1.5 day. A spectrum of the host galaxy taken 111 days after the burst reveals a single emission line, interpreted as [OII] at a redshift z = 1.0585, and a continuum break which we interpret as the Balmer limit at this redshift. Despite the long T_90, the afterglow of GRB 000911 is not unusual in any other way when compared to the set of afterglows studied to date. We conclude that the duration of the GRB plays little part in determining the physics of the afterglow.

One major component of the VO will be catalogs measuring gigabytes and terrabytes if not more. Some mechanism like XML will be used for structuring the information. However, such mechanisms are not good for information retrieval on their own. For retrieval we use queries. Topic Maps that have started becoming popular recently are excellent for segregating information that results from a query. A Topic Map is a structured network of hyperlinks above an information pool. Different Topic Maps can form different layers above the same information pool and provide us with different views of it. This facilitates in being able to ask exact questions, aiding us in looking for gold needles in the proverbial haystack. Here we discuss the specifics of what Topic Maps are and how they can be implemented within the VO framework. URL: http://www.astro.caltech.edu/~aam/science/topicmaps/

We present millimeter and submillimeter wavelength observations and near-infrared K-band imaging toward the bright gamma-ray burst GRB 010222. Over seven epochs the flux density of the source was constant with an average flux density 3.74+/-0.53 mJy at 350 GHz and 1.05+/-0.22 mJy at 250 GHz, giving a spectral index alpha=3.78+/-0.25. We rule out the possibility that this emission originated from the burst or its afterglow and we conclude that it is due to a dusty, high redshift starburst galaxy (SMM J14522+4301). We argue that the host galaxy of GRB 010222 is the most plausible counterpart of SMM J14522+4301, based in part on the centimeter detection of the host at the expected level. The optical/NIR properties of the host galaxy of GRB 010222 suggest that it is a blue, sub-L_*, similar to other GRB host galaxies. This contrasts with the enormous far-infrared luminosity of this galaxy based on our submillimeter detection (L_Bol~4X10^{12} L_odot). We suggest that this GRB host galaxy has a very high star formation rate, SFR 600 M_\odot/yr, most of which is unseen at optical wavelengths.

Discovery of the cosmic reionization epoch would represent a significant milestone in cosmology. We present Keck spectroscopy of the quasar SDSS 1044-0125, at z = 5.73. The spectrum shows a dramatic increase in the optical depth at observed wavelengths lambda >~7550 A, corresponding to z_abs >~ 5.2. Only a few small, narrow transmission regions are present in the spectrum beyond that point, and out to the redshifts where the quasar signal begins. We interpret this result as a signature of the trailing edge of the cosmic reionization epoch, which we estimate to occur around <z> ~ 6 (as indeed confirmed by subsequent observations by Becker et al.), and extending down to z \~ 5.2. This behavior is expected in the modern theoretical models of the reionization era, which predict a patchy and gradual onset of reionization. The remaining transmission windows we see may correspond to the individual reionization bubbles (Stromgren spheres) embedded in a still largely neutral intergalactic medium, intersected by the line of sight to the quasar. Future spectroscopic observations of quasars at comparable or larger redshifts will provide a more detailed insight into the structure and extent of the reionization era.

As a scientific discipline, Astronomy is rather unique. We only have one laboratory, the Universe, and we cannot, of course, change the initial conditions and study the resulting effects. On top of this, acquiring Astronomical data has historically been a very labor-intensive effort. As a result, data has traditionally been preserved for posterity. With recent technological advances, however, the rate at which we acquire new data has grown exponentially, which has generated a Data Tsunami, whose wave train threatens to overwhelm the field. In this conference proceedings, we present and define the concept of virtual observatories, which we feel is the only logical answer to this dilemma.

Like every other field of intellectual endeavor, astronomy is being revolutionised by the advances in information technology. There is an ongoing exponential growth in the volume, quality, and complexity of astronomical data sets, mainly through large digital sky surveys and archives. The Virtual Observatory (VO) concept represents a scientific and technological framework needed to cope with this data flood. Systematic exploration of the observable parameter spaces, covered by large digital sky surveys spanning a range of wavelengths, will be one of the primary modes of research with a VO. This is where the truly new discoveries will be made, and new insights be gained about the already known astronomical objects and phenomena. We review some of the methodological challenges posed by the analysis of large and complex data sets expected in the VO-based research. The challenges are driven both by the size and the complexity of the data sets (billions of data vectors in parameter spaces of tens or hundreds of dimensions), by the heterogeneity of the data and measurement errors, including differences in basic survey parameters for the federated data sets (e.g., in the positional accuracy and resolution, wavelength coverage, time baseline, etc.), various selection effects, as well as the intrinsic clustering properties (functional form, topology) of the data distributions in the parameter spaces of observed attributes. Answering these challenges will require substantial collaborative efforts and partnerships between astronomers, computer scientists, and statisticians.

Observations of GRB host galaxies and their environments in general can provide valuable clues about the nature of progenitors. Bursts are associated with faint, <R> ~ 25 mag, galaxies at cosmological redshifts, <z> ~ 1. The host galaxies span a range of luminosities and morphologies, and appear to be broadly typical for the normal, evolving, actively star-forming galaxy populations at comparable redshifts and magnitudes, but may have somewhat elevated SFR per unit luminosity. There are also spectroscopic hints of massive star formation, from the ratios of [Ne III] and [O II] lines. The observed, unobscured star formation rates are typically a few M_sun/yr, but a considerable fraction of the total star formation in the hosts may be obscured by dust. A census of detected optical afterglows provides a powerful new handle on the obscured fraction of star formation in the universe; the current results suggest that at most a half of the massive star formation was hidden by dust.

GRB 970828 was the first well-localized gamma-ray burst (GRB) X-ray afterglow for which no optical afterglow was found despite a prompt, deep search down to R_lim ~ 24.5 mag. We report the discovery of a short-lived radio flare within the X-ray localization error circle of this burst. Such radio flares are seen in about 25% of GRB afterglows and their origin is not well understood. The precise radio position enabled us to identify the likely host galaxy of this burst, and to measure its redshift, z = 0.9578. The host appears to be an interacting/merging system. Under the assumption that the X-ray afterglow is mainly due to synchrotron mechanism we infer the optical afterglow flux. The observed upper limits to the optical flux are easily explained by invoking an intervening dusty cloud within the host galaxy. These observations support the idea that GRBs with no detectable optical afterglows or "dark GRBs" can be due to dust extinction within the host galaxies. The census of dark GRBs can then be used to constrain the fraction of the obscured star formation in the universe. We argue that the existing data already indicate that the obscured star-formation rate is no more than 1/2 of that seen at UV and optical wavelengths.

Cosmic gamma-ray bursts are one of the great frontiers of astrophysics today. They are a playground of relativists and observers alike. They may teach us about the death of stars and the birth of black holes, the physics in extreme conditions, and help us probe star formation in the distant and obscured universe. In this review we summarise some of the remarkable progress in this field over the past few years. While the nature of the GRB progenitors is still unsettled, it now appears likely that at least some bursts originate in explosions of very massive stars, or at least occur in or near the regions of massive star formation. The physics of the burst afterglows is reasonably well understood, and has been tested and confirmed very well by the observations. Bursts are found to be beamed, but with a broad range of jet opening angles; the mean gamma-ray energies after the beaming corrections are ~ 10^51 erg. Bursts are associated with faint <R> ~ 25 mag) galaxies at cosmological redshifts, with <z> ~ 1. The host galaxies span a range of luminosities and morphologies, but appear to be broadly typical for the normal, actively star-forming galaxy populations at comparable redshifts and magnitudes. Some of the challenges for the future include: the nature of the short bursts and possibly other types of bursts and transients; use of GRBs to probe the obscured star formation in the universe, and possibly as probes of the very early universe; and their detection as sources of high-energy particles and gravitational waves.

Astronomy has a long history of acquiring, systematizing, and interpreting large quantities of data. Starting from the earliest sky atlases through the first major photographic sky surveys of the 20th century, this tradition is continuing today, and at an ever increasing rate. Like many other fields, astronomy has become a very data-rich science, driven by the advances in telescope, detector, and computer technology. Numerous large digital sky surveys and archives already exist, with information content measured in multiple Terabytes, and even larger, multi-Petabyte data sets are on the horizon. Systematic observations of the sky, over a range of wavelengths, are becoming the primary source of astronomical data. Numerical simulations are also producing comparable volumes of information. Data mining promises to both make the scientific utilization of these data sets more effective and more complete, and to open completely new avenues of astronomical research. Technological problems range from the issues of database design and federation, to data mining and advanced visualization, leading to a new toolkit for astronomical research. This is similar to challenges encountered in other data-intensive fields today. These advances are now being organized through a concept of the Virtual Observatories, federations of data archives and services representing a new information infrastructure for astronomy of the 21st century. In this article, we provide an overview of some of the major datasets in astronomy, discuss different techniques used for archiving data, and conclude with a discussion of the future of massive datasets in astronomy.

GRB 000926 has one of the best-studied afterglows to-date, with multiple X-ray observations, as well as extensive multi-frequency optical and radio coverage. Broadband afterglow observations, spanning from X-ray to radio frequencies, provide a probe of the density structure of the circumburst medium, as well as of the ejecta energetics, geometry, and the physical parameters of the relativistic blastwave resulting from the explosion. We present an analysis of {\em Chandra X-ray Observatory} observations of this event, along with {\em Hubble Space Telescope} and radio monitoring. We combine these data with ground-based optical and IR observations and fit the synthesized afterglow lightcurve using models where collimated ejecta expand into a surrounding medium. We find that we can explain the broadband lightcurve with reasonable physical parameters only if the cooling is dominated by inverse Compton scattering. Excess X-ray emission in the broadband spectrum indicates that we are directly observing a contribution from inverse Compton scattering. It is the first time this has been observed in a GRB afterglow, and it implies that the GRB exploded in a reasonably dense (n~30 cm^{-3}) medium, consistent with a diffuse interstellar cloud environment.

(Abridged version) Mayall II = G1 is one of the brightest globular clusters belonging to M31, the Andromeda galaxy. Our observations with HST/WFPC2 provide data for the (I vs. V-I) and (V vs. V-I) color-magnitude diagrams. From model fitting, we determine a rather high mean metallicity of [Fe/H] = --0.95 +- 0.09, somewhat similar to 47 Tucanae. We find a larger spread in V-I than can be explained by the measurement errors, and we attribute this to an intrinsic metallicity dispersion amongst the stars of G1. So far, only omega Centauri, the giant Galactic globular cluster, has been known to exhibit such an intrinsic metallicity dispersion. Three estimates of the total mass of this globular cluster can be obtained: King mass = 15 x 10^6 with M/Lv ~ 7.5, Virial mass = 7.3 x 10^6 with M/Lv ~ 3.6, and King-Michie mass range from 14 to 17 x 10^6. Although uncertain, all of these mass estimates make G1 more than twice as massive as omega Centauri. Such large masses relate to the metallicity spread whose origin is still unknown (either self-enrichment, an inhomogeneous proto-cluster cloud, or remaining core of a dwarf galaxy). When considering the positions of G1 in the different diagrams defined by Kormendy (1985), G1 always appears on the sequence defined by globular clusters, and definitely away from the other sequences defined by elliptical galaxies, bulges, and dwarf spheroidal galaxies. The same is true for omega Centauri and for the nucleus of the dwarf elliptical NGC 205. This does not prove that all (massive) globular clusters are the remnant cores of nucleated dwarf galaxies.

We present deep infrared ($2.2 \mu$m) imaging of the Galactic microquasars 1E1740-2942 and GRS 1758-258 using the Keck-I 10-meter telescope in June 1998. The observations were taken under excellent seeing conditions ($\sim 0.45 \arcsec$ full-width half-maximum), making them exceptionally deep for these crowded fields. We used the USNO-A2.0 catalog to astrometrically calibrate the infrared images (along with an optical CCD image in the case of GRS 1758-258), providing independent frame ties to the known radio positions of the objects. For 1E1740-2942, we confirm potential candidates for the microquasar previously identified by Marti et al., and show that none of the objects near the microquasar have varied significantly from 1998 to 1999. For GRS 1758-258, our astrometry indicates a position shifted from previous reports of candidates for the microquasar. We find no candidates inside our 90% confidence radius to a $2 \sigma$ limiting magnitude of $K_s = 20.3$ mag. We discuss the implications of these results for the nature of the microquasar binary systems.

Gamma-ray bursts (GRBs) are the most brilliant objects in the Universe but efforts to estimate the total energy released in the explosion -- a crucial physical quantity -- have been stymied by their unknown geometry: spheres or cones. We report on a comprehensive analysis of GRB afterglows and derive their conical opening angles. We find that the gamma-ray energy release, corrected for geometry, is narrowly clustered around 5x10**50 erg. We draw three conclusions. First, the central engines of GRBs release energies that are comparable to ordinary supernovae, suggesting a connection. Second, the wide variation in fluence and luminosity of GRBs is due entirely to a distribution of opening angles. Third, only a small fraction of GRBs are visible to a given observer and the true GRB rate is at least a factor of 500 times larger than the observed rate.

We present multi-color light-curves of the optical afterglow of GRB 000926. Beginning ~1.5 days after the burst, the light-curves of this GRB steepen measurably. The existence of such achromatic breaks are usually taken to be an important observational signature that the ejecta are not expanding isotropically, but rather have a collimated jet-like geometry. If we interpret the data in this context, we derive an opening angle of 5 degrees, which reduces the energy release compared to an isotropic model by a factor of ~275, to 1.7 x 10^51 erg. To fit the data with a simple jet model requires extinction along the line of sight. The derived A_V is in the range 0.11 - 0.82 mag, depending on the adopted extinction law and whether the electrons giving rise to the optical emission are undergoing synchrotron cooling or not. Since this is in excess of the expected extinction from our Galaxy, we attribute this to the GRB host. We note that this extinction is typical of a galactic disk, and therefore the event likely took place in the disk of its host.

We review some of the scientific opportunities and technical challenges posed by the exploration of the large digital sky surveys, in the context of a Virtual Observatory (VO). The VO paradigm will profoundly change the way observational astronomy is done. Clustering analysis techniques can be used to discover samples of rare, unusual, or even previously unknown types of astronomical objects and phenomena. Exploration of the previously poorly probed portions of the observable parameter space are especially promising. We illustrate some of the possible types of studies with examples drawn from DPOSS; much more complex and interesting applications are forthcoming. Development of the new tools needed for an efficient exploration of these vast data sets requires a synergy between astronomy and information sciences, with great potential returns for both fields.

Systematic exploration of the observable parameter space, covered by large digital sky surveys spanning a range of wavelengths, will be one of the primary modes of research with a Virtual Observatory (VO). This will include searches for rare, unusual, or even previously unknown types of astronomical objects and phenomena, e.g. as outliers in some parameter space of measured properties, both in the catalog and image domains. Examples from current surveys include high-redshift quasars, type-2 quasars, brown dwarfs, and a small number of objects with puzzling spectra. Opening of the time domain will be especially interesting in this regard. Data-mining tools such as unsupervised clustering techniques will be essential in this task, and should become an important part of the VO toolkit.

The composite galaxy luminosity function (hereafter LF) of 39 Abell clusters of galaxies is derived by computing the statistical excess of galaxy counts in the cluster direction with respect to control fields. Due to the wide field coverage of the digitised POSS-II plates, we can measure field counts around each cluster in a fully homogeneous way. Furthermore, the availability of virtually unlimited sky coverage allows us to directly compute the LF errors without having to rely on the estimated variance of the background. The wide field coverage also allows us to derive the LF of the whole cluster, including galaxies located in the cluster outskirts. The global composite LF has a slope alpha ~ -1.1+/-0.2 with minor variations from blue to red filters, and M* ~ -21.7,-22.2,-22.4 mag (H_0=50 km/s/Mpc) in g, r and i filters, respectively. These results are in quite good agreement with several previous determinations and in particular with the LF determined for the inner region of a largely overlapping set of clusters, but derived making use of a completely different method for background subtraction. The similarity of the two LFs suggests the existence of minor differences between the LF in the cluster outskirts and in the central region, or a negligible contribution of galaxies in the cluster outskirts to the global LF.

Astronomy is entering a new era as multiple, large area, digital sky surveys are in production. The resulting datasets are truly remarkable in their own right; however, a revolutionary step arises in the aggregation of complimentary multi-wavelength surveys (i.e., the cross-identification of a billion sources). The federation of these large datasets is already underway, and is producing a major paradigm shift as Astronomy has suddenly become an immensely data-rich field. This new paradigm will enable quantitatively and qualitatively new science, from statistical studies of our Galaxy and the large-scale structure in the universe, to discoveries of rare, unusual, or even completely new types of astronomical objects and phenomena. Federating and then exploring these large datasets, however, is an extremely challenging task. The Digital Sky project was initiated with this task in mind and is working to develop the techniques and technologies necessary to solve the problems inherent in federating these large databases, as well as the mining of the resultant aggregate data.

We examine the Butcher-Oemler effect and its cluster richness dependence in the largest sample studied to date: 295 Abell clusters. We find a strong correlation between cluster richness and the fraction of blue galaxies, f_B, at every redshift. The slope of the f_B(z) relation is similar for all richnesses, but at a given redshift, f_B is systematically higher for poor clusters. This is the chief cause of scatter in the f_B vs. z diagram: the spread caused by the richness dependence is comparable to the trend in f_B over a typical redshift baseline, so that conclusions drawn from smaller samples have varied widely. The two parameters, z, and a consistently defined projected galaxy number density, N, together account for all of the observed variation in f_B within the measurement errors. The redshift evolution of f_B is real, and occurs at approximately the same rate for clusters of all richness classes.

Our current understanding of the low surface brightness universe is quite incomplete, not only in the optical, but also in other wavelength regimes. As a demonstration of the type of science which is facilitated by a virtual observatory, we have undertaken a project utilizing both images and catalogs to explore the multi-wavelength, low surface brightness universe. Here, we present some initial results of this project. Our techniques are complimentary to normal data reduction pipeline techniques in that we focus on the diffuse emission that is ignored or removed by more traditional algorithms. This requires a spatial filtering which must account for objects of interest, in addition to observational artifacts (e.g., bright stellar halos). With this work we are exploring the intersection of the catalog and image domains in order to maximize the scientific information we can extract from the federation of large survey data.

We present the composite luminosity function (hereafter LF) of galaxies for 81 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data processed by the CRoNaRio collaboration. The derived LF is very accurate due to the use of homogeneous data both for the clusters and the control fields and to the local estimate of the background, which takes into account the presence of large-scale structures and of foreground clusters and groups. The global composite LF is quite flat down to $M^*+5$ has a slope $\alpha\sim-1.0\pm0.2$ with minor variations from blue to red filters, and $M^*\sim-21.8,-22.0,-22.3$ mag ($H_0=50$ km s$^{-1}$ Mpc$^{-1}$) in the $g, r$ and $i$ filters, respectively (errors are detailed in the text). We find a significant difference between rich and poor clusters thus arguing in favour of a dependence of the LF on the properties of the environment.

We present the composite luminosity function (hereafter LF) of galaxies for 80 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data in the framework of the CRoNaRio collaboration. Our determination of the LF has been computed with very high accuracy thanks to the use of homogeneous data both for the clusters and the control fields and to a local estimate of the background, which takes into account the presence of large-scale structures and of foreground clusters and groups. The global composite LF is quite flat down to M^+5 and it has a slope alpha=-1.0 +/- 0.2 with minor variations from blue to red filters, and M^* = -21.9,-22.0,-22.3 mag (H_0=50 km s^-1 Mpc^-1) in g, r and i filters, respectively (errors are detailed in the text). We find a significant difference between rich and poor clusters supporting the existence of an LF's dependence on the environment.

We discuss the implementation and validation of a procedure aimed to detect groups and clusters of galaxies in DPOSS catalogs.

We present optical and near-IR identification and spectroscopy of the host galaxy of GRB 980613. The burst was apparently associated with the optically (restframe UV) brightest component of a system of at least five galaxies or galaxy fragments at a redshift of z = 1.0969. The component we identify as the host galaxy shows a moderately high unobscured star formation rate, SFR = 5 M_solar/yr, but a high SFR per unit mass, indicative of a starburst. The image components show a broad range of (R-K) colors, with two of them being very red, possibly due to dust. Overall morphology of the system can be naturally interpreted as a strong tidal interaction of two or more galaxies, at a redshift where such events were much more common than now. Given the well established causal link between galaxy mergers and starbursts, we propose that this is a strong case for a GRB originating from a merger-induced starburst system. This supports the proposed link between GRBs and massive star formation.

The gamma-ray burst of 28 February 1997 (GRB 970228) ushered in the discovery of the afterglow phenomenon. Despite intense study of the host galaxy, however, the nature of the host and the distance to the burst eluded the community. Here we present the measurement of the redshift of GRB 970228, and, based on its spectroscopic and photometric properties, identify the host as an underluminous, but otherwise normal galaxy at redshift z=0.695 undergoing a modest level of star formation. At this redshift, the GRB released an isotropic equivalent energy of E = (6.8 +/- 0.5) x 10^51 erg [30-2000 keV restframe]. We find no evidence that the host is much bluer or forming stars more rigorously than the general field population. In fact, by all accounts in our analysis (color-magnitude, magnitude-radius, star-formation rate, Balmer-break amplitude) the host properties appear typical for faint blue field galaxies at comparable redshifts.

We constrain the mass-to-light ratios, gas mass fractions, baryon mass fractions and the ratios of total to luminous mass for a sample of eight nearby relaxed galaxy groups and clusters: A262, A426, A478, A1795, A2052, A2063, A2199 and MKW4s. We use ASCA spatially resolved spectroscopic X-ray observations and ROSAT PSPC images to constrain the total and gas masses of these clusters. To measure cluster luminosities we use galaxy catalogs resulting from the digitization and automated processing of the second generation Palomar Sky Survey plates calibrated with CCD images in the Gunn-Thuan g, r, and i bands. Under the assumption of hydrostatic equilibrium and spherical symmetry, we can measure the total masses of clusters from their intra-cluster gas temperature and density profiles. Spatially resolved ASCA spectra show that the gas temperature decreases with increasing distance from the center. By comparison, the assumption that the gas is isothermal results in an underestimate of the total mass at small radii, and an overestimate at large cluster radii. We have obtained luminosity functions for all clusters in our sample. After correcting for background and foreground galaxies, we estimate the total cluster luminosity using Schechter function fits to the galaxy catalogs. In the three lowest redshift clusters where we can sample to fainter absolute magnitudes, we have detected a flattening of the luminosity function at intermediate magnitudes and a rise at the faint end. These clusters were fit with a sum of two Schechter functions. The remaining clusters were well fit with a single Schechter function.

Large digital sky surveys, over a broad range of wavelengths, both from the ground and from space observatories, are becoming a major source of astronomical data. Some examples include the Sloan Digital Sky Survey (SDSS) and the Digital Palomar Observatory Sky Survey (DPOSS) in the visible, the Two-Micron All-Sky Survey (2MASS) in the near-infrared, the NRAO VLA Sky Survey (NVSS) and the Faint Images of the Radio Sky at Twenty centimeters (FIRST) in the radio. Many others surveys are planned or expected, in addition to the previously named surveys. While most surveys are exclusively imaging, large-scale spectroscopic surveys also exist. In addition, a number of experiments with specific scientific goals, e.g., microlensing surveys for MACHOs, searches for near-Earth asteroids, are generating comparable data volumes. Typical sizes of resulting data sets (as of the late 1990's) are in the range of tens of Terabytes of digital information, with detections of many millions or even billions of sources, and several tens of parameters measured for each detected source. This vast amount of new information presents both a great scientific opportunity and a great technological challenge: how to process, and calibrate the raw data; how to store, combine, and access them using modern computing hardware and networks; and how to visualize, explore and analyses these great data sets quickly and efficiently.

We report on the discovery of the radio afterglow from the gamma-ray burst (GRB) of 1999 May 6 (GRB 990506) using the Very Large Array (VLA). The radio afterglow was detected at early times (1.5 days), but began to fade rapidly sometime between 1 and 5 days after the burst. If we attribute the radio emission to the forward shock from an expanding fireball, then this rapid onset of the decay in the radio predicts that the corresponding optical transient began to decay between 1 and 5 minutes after the burst. This could explain why no optical transient for GRB 990506 was detected in spite of numerous searches. The cause of the unusually rapid onset of the decay for the afterglow is probably the result of an isotropically energetic fireball expanding into a low density circumburst environment. At the location of the radio afterglow we find a faint (R ~ 24 mag) host galaxy with a double morphology.

This paper presents photometric redshifts for 431 Abell clusters imaged as part of the Palomar Abell Cluster Optical Survey (PACOS), of which 236 are new redshi fts. We have obtained moderately deep, 3--band (Gunn gri) imaging for this sam ple at the Palomar Observatory 60'' telescope, as part of the photometric calibration of DPOSS. Our data acquisition, reduction, and photometric calibration techniques are described, and photometric accuracy and consistency is demonstrated. An empirical redshift estimator is presented, utilizing background-corrected median g-r colors and mean g magnitudes for the ensemble of galaxies in each field. We present photometric redshift estimates for the clusters in our sample with an accuracy of sigma_z=0.038. These redshift estimates provide checks on single-galaxy cluster redshifts, as well as distance information for studies of the Butcher-Oemler effect, luminosity functions, M/L ratios, and many other projects.

We estimate the empirical R parameter in 26 Galactic Globular Clusters covering a wide metallicity range, imaged by WFPC2 on board the HST. The improved spatial resolution permits a large fraction of the evolved stars to be measured and permits accurate assessment of radial populaton gradients and completeness corrections. In order to evaluate both the He abundance and the He to metal enrichment ratio, we construct a large set of evolutionary models by adopting similar metallicities and different He contents. We find an absolute He abundance which is lower than that estimated from spectroscopic measurements in HII regions and from primordial nucleosynthesis models. This discrepancy could be removed by adopting a C12O16 nuclear cross section about a factor of two smaller than the canonical value, although also different assumptions for mixing processes can introduce systematical effects. The trend in the R parameter toward solar metallicity is consistent with an upper limit to the He to metal enrichment ratio of the order of 2.5.

We exploited the large areal coverage offered by the DPOSS to analyze the outermost regions of the galactic globular cluster M 92. Two independent photometric reduction programs (SKICAT and DAOPHOT) were used to construct a color-magnitude diagram and a surface density profile for this cluster, based on J- and F-band DPOSS plates. A strong similarity has been found in the performance of the two programs in the low--crowded outermost cluster regions. After removing the background contribution, we obtained the cluster outer surface density profile down to a surface brightness magnitude of mu_V ~ 31 mag/sq.arcsec and matched it with the inner profile of Trager et al. 1995. The profile shapes match very well: since our data are uncalibrated, the shift in magnitude between the profiles has been also used to calibrate our profile. The analysis shows that the cluster has an extra tidal halo extending out to ~ 30 arcmin from the cluster center at a 3 sigma level over the background noise. This halo is revealed to be almost circular.

Gamma-ray burst astronomy has undergone a revolution in the last three years, spurred by the discovery of fading long-wavelength counterparts. We now know that at least the long duration GRBs lie at cosmological distances with estimated electromagnetic energy release of 10**51 -- 10**53 erg, making these the brightest explosions in the Universe. In this article we review the current observational state, beginning with the statistics of X-ray, optical, and radio afterglow detections. We then discuss the insights these observations have given to the progenitor population, the energetics of the GRB events, and the physics of the afterglow emission. We focus particular attention on the evidence linking GRBs to the explosion of massive stars. Throughout, we identify remaining puzzles and uncertainties, and emphasize promising observational tools for addressing them. The imminent launch of HETE-2 and the increasingly sophisticated and coordinated ground-based and space-based observations have primed this field for fantastic growth.

We report on two programs to address differential evolution between the radio-loud and radio-quiet quasar populations at high (z>4) redshift. Both programs entail studying the radio properties of optically-selected quasars. First, we have observed 32 optically-selected, high-redshift (z>4) quasars with the VLA at 6 cm (5 GHz). These sources comprise a statistically complete and well-understood sample. We detect four quasars above our 3-sigma limit of ~0.15 mJy, which is sufficiently sensitive to detect all radio-loud quasars at the probed redshift range. Second, we have correlated 134 z>4 quasars, comprising all such sources that we are aware of as of mid-1999, with FIRST and NVSS. These two recent 1.4 GHz VLA sky surveys reach 3-sigma limits of approximately 0.6 mJy and 1.4 mJy respectively. We identify a total of 15 z>4 quasars, of which six were not previously known to be radio-loud. The depth of these surveys does not reach the radio-loud/radio-quiet demarcation luminosity density (L(1.4 GHz) = 10^32.5 h(50)^(-2) ergs/s/Hz) at the redshift range considered; this correlation therefore only provides a lower limit to the radio-loud fraction of quasars at high-redshift. The two programs together identify eight new radio-loud quasars at z>4, a significant increase over the seven currently in the published literature. We find no evidence for radio-loud fraction depending on optical luminosity for -25 > M_B > -28 at z~2, or for -26>M_B>-28 at z>4. Our results also show no evolution in the radio-loud fraction between z~2 and z>4 (-26>M_B>-28).

Images of the Hubble Deep Field at 3.2 um, taken with the 10-m Keck Telescope, are presented. The images cover a total area of ~2.5 arcmin^2. To a 5-sigma limit of [3.2](total)=17.5 mag (Vega-relative), 11 sources are detected, 9 of which are extragalactic. The integrated galaxy number counts are therefore ~1.3x10^4 deg^{-2} to this depth. The galaxies detected at 3.2 um have a median redshift of z=0.56. All 3.2 um sources have 1.6 um, 1.1 um and visual counterparts, all of fairly regular morphologies; several also have 6.7 um, 15 um, 850 um, 8.5 GHz, or 1.4 GHz counterparts. No sources are found which are either anomalously red or anomalously blue in their H-[3.2] color, and there are significant detections of most of the known near-infrared HDF sources for which detection in these 3.2 um data seemed likely.

Quasars at large redshifts provide a powerful probe of structure formation in the early universe. Several arguments suggest that the formation of ellipticals and massive bulges may have involved an early quasar phase. At very large redshifts, such structures are likely to be found at the highest peaks of the density field, and would thus be highly biased tracers: the earliest (massive) galaxy formation may have occurred in the cores of future rich clusters. Preliminary results from our search for clustered protogalaxies around quasars at z > 4 support this idea. Quasars at even larger redshifts may be an important contributor to the reionisation of the universe, and signposts of the earliest galaxy and cluster formation.

The Northern Sky Optical Cluster Survey is a project to create an objective catalog of galaxy clusters over the entire high-galactic-latitude Northern sky, with well understood selection criteria. We use the object catalogs generated from the Digitized Second Palomar Sky Survey (DPOSS, Djorgovski et al. 1999) as the basis for this survey. We apply a color criterion to select against field galaxies, and use a simple adaptive kernel technique to create galaxy density maps, combined with the bootstrap technique to make significance maps, from which density peaks are selected. This survey attempts to eliminate some of the subjective criteria and assumptions of past surveys, including detection by eye (Abell 1958, ACO 1989) and assumed luminosity functions and cluster profiles (PDCS, Postman et al. 1995). We also utilize more information (especially colors) than the most similar recent survey, the APM (Dalton et al. 1992). This paper presents the details of our cluster detection technique, as well as some initial results for two small areas totaling ~60 square degrees. We find a mean surface density of ~1.5 clusters per square degree, consistent with the detection of richness class 0 and higher clusters to z~0.3. In addition, we demonstrate an effective photometric redshift estimator for our clusters.

We report the discovery of a radio transient VLA 232937.2-235553, coincident with the proposed X-ray afterglow for the gamma-ray burst GRB 981226. This GRB has the highest ratio of X-ray to gamma-ray fluence of all the GRBs detected by BeppoSAX so far and yet no corresponding optical transient was detected. The radio light curve of VLA 232937.2-235553 is qualitatively similar to that of several other radio afterglows. At the sub-arcsecond position provided by the radio detection, optical imaging reveals an extended R=24.9 mag object, which we identify as the host galaxy of GRB 981226. Afterglow models which invoke a jet-like geometry for the outflow or require an ambient medium with a radial density dependence, such as that produced by a wind from a massive star, are both consistent with the radio data. Furthermore, we show that the observed properties of the radio afterglow can explain the absence of an optical transient without the need for large extinction local to the GRB.

We present preliminary results of a simple peak finding algorithm applied to catalogues of galaxies, extracted from the Second Palomar Sky Survey in the framework of the CRoNaRio project. All previously known Abell and Zwicky clusters in a test region of 5x5 sq. deg. are recovered and new candidate clusters are also detected. This algorithm represents an alternative way of searching for galaxy clusters with respect to that implemented and tested at Caltech on the same type of data (Gal et al. 1998).

The projected multiplicity function of galaxies gives (per square degree) the density of galaxy aggregates formed by N members. We use the CRoNaRio matched catalogues to derive the low N (N<15) multiplicity function from D-POSS data. The van Albada (1982) algorithm was implemented and used to identify candidate groups in the CRoNaRio catalogues. In absence of redshift surveys complete down to the magnitude limit of the DPOSS material and covering a wide enough area, the performances of the algorithm were tested on realistically simulated catalogues. The application to a set of 13 CRoNaRio catalogues allowed us to derive a multiplicity function as accurate as those available in literature obtained from redshift surveys.

We present the composite luminosity function of galaxies for 24 Abell clusters studied in our survey of the Northern Hemisphere, using DPOSS data in the framework of the CRoNaRio collaboration. Our determination of the luminosity function has been computed with very high accuracy thanks to 1) the use of homogeneous data both in the cluster and in the control field, 2) a local estimate of the background, which takes into account the presence of large-scale structures (which tipically make the background contribution larger) and of foreground clusters and groups; 3) the inclusion in the error calculation of the variance of background counts, which is tipically larger than Poissonian fluctuation. We plan at least a tenfold increase of the number of clusters.

We have obtained near-infrared (1.6 micron) images of 11 powerful 3CR radio galaxies at redshifts 0.8 < z < 1.8 using NICMOS on board HST. The high angular resolution permits a detailed study of galaxy morphology in these systems at rest-frame optical wavelengths, where starlight dominates over the extended, aligned UV continuum. The NICMOS morphologies are mostly symmetric and are consistent with dynamically relaxed, elliptical host galaxies dominated by a red, mature stellar population. The aligned structures are sometimes faintly visible, and nuclear point sources may be present in a few cases which manifest the ``unveiled'' AGN that is obscured from view at optical wavelengths. Our observations are consistent with the hypothesis that the host galaxies of z ~ 1-2 radio galaxies are similar to modern-day gE galaxies. Their sizes are typical of gE galaxies but smaller than present-day cD and brightest cluster galaxies, and their surface brightnesses are higher, as expected given simple luminosity evolution.

We present a preliminary measurement of the angular clustering of faint ($R \le 25$) field galaxies in which we concentrate on the behavior of $w(\theta)$ on small angular scales ($\theta \ls 10''$). The galaxies are strongly clustered and $w(\theta)$ is well-characterized by a power law of the form $A_w \theta^{-\delta}$. The best-fitting value of the power law index, $\delta$, is, however, steeper than the fiducial value of $\delta = 0.8$, indicating there are more pairs of galaxies separated by $\theta \ls 10''$ in our sample than would be otherwise expected. Using the best-fitting form of $w(\theta)$, we estimate that $\sim 10%$ of the galaxies are in physically close pairs (separations $\ls 21 h^{-1}$ kpc). This is a factor of order 2 larger than local galaxy samples, but comparable to galaxy samples with $<z> \sim 0.4$. The mean redshift of our galaxies is of order 0.95, and, therefore, our result suggests that there was little or no evolution in the merger rate of galaxies between $z \sim 1$ and $z \sim 0.4$.

We report on the current status of our search for protoclusters around quasars at z > 4. While the search is still very incomplete, clustered companion galaxies are found in virtually every case examined so far. The implied comoving number densities of protogalaxies are two to four orders of magnitude higher than expected for the general field, but are comparable to the number densities in rich cluster cores. The comoving densities of star formation in these regions are also enhanced by a comparable factor. We interpret these results as an evidence for biased galaxy formation in the highest peaks of the primordial density field.

The optical transient of the faint Gamma Ray Burst 990308 was detected by the QUEST camera on the Venezuelan 1-m Schmidt telescope starting 3.28 hours after the burst. Our photometry gives $V = 18.32 \pm 0.07$, $R = 18.14 \pm 0.06$, $B = 18.65 \pm 0.23$, and $R = 18.22 \pm 0.05$ for times ranging from 3.28 to 3.47 hours after the burst. The colors correspond to a spectral slope of close to $f_{\nu} \propto \nu^{1/3}$. Within the standard synchrotron fireball model, this requires that the external medium be less dense than $10^{4} cm^{-3}$, the electrons contain $> 20%$ of the shock energy, and the magnetic field energy must be less than 24% of the energy in the electrons for normal interstellar or circumstellar densities. We also report upper limits of $V > 12.0$ at 132 s (with LOTIS), $V > 13.4$ from 132-1029s (with LOTIS), $V > 15.3$ at 28.2 min (with Super-LOTIS), and a 8.5 GHz flux of $< 114 \mu Jy$ at 110 days (with the Very Large Array). WIYN 3.5-m and Keck 10-m telescopes reveal this location to be empty of any host galaxy to $R > 25.7$ and $K > 23.3$. The lack of a host galaxy likely implies that it is either substantially subluminous or more distant than a red shift of $\sim 1.2$.

We present multi-color optical and two-frequency radio observations of the bright SAX event, GRB 990510. The well-sampled optical decay, together with the radio observations are inconsistent with simple spherical afterglow models. The achromatic optical steepening and the decay of the radio afterglow both occuring at $t \sim 1$ day are evidence for hydrodynamical evolution of the source, and can be most easily interpreted by models where the GRB ejecta are collimated in a jet. Employing a simple jet model to interpret the observations, we derive a jet opening angle of $\theta_o = 0.08$, reducing the isotropic gamma-ray emission of $2.9 \times 10^{53}$ erg by a factor $\sim 300$. If the jet interpretation is correct, we conclude that GRB observations to-date are consistent with an energy for the central source of $E \lsim 10^{52}$ erg.

We have determined the projected rotational velocities of thirteen blue horizontal-branch (BHB) stars in the globular cluster M13 via rotational broadening of metal absorption lines. Our sample spans the photometric gap observed in the horizontal branch distribution at Teff ~= 11000 K, and reveals a pronounced difference in stellar rotation on either side of this feature---bluewards of the gap, all the stars show modest rotations, vsini < 10 km/s, while to the red side of the gap, we confirm the more rapidly rotating population (vsini ~= 40 km/s) previously observed by Peterson et al (1995). Taken together with these prior results, our measurements indicate that a star's rotation is indeed related to its location along the HB, although the mechanism behind this correlation remains unknown. We explore possible connections between stellar rotation and mass loss mechanisms which influence the photometric morphology of globular cluster HBs.

We present ground-based images of the z=1.824 radio galaxy 3C 256 in the standard BVRIJHK filters and an interference filter centered at 8800A, a Hubble Space Telescope image in a filter dominated by Ly-alpha emission (F336W), and spectra covering rest-frame wavelengths from Ly-alpha to [O III] 5007. Together with published polarimetry observations, we use these to decompose the overall spectral energy distribution into nebular continuum emission, scattered quasar light, and stellar emission. The nebular continuum and scattered light together comprise half (one third) of the V-band (K-band) light within a 4-arcsec aperture, and are responsible for the strong alignment between the optical/near-infrared light and the radio emission. The stellar emission is dominated by a population estimated to be 100-200 Myr old (assuming a Salpeter IMF), and formed in a short burst with a peak star formation rate of 1-4x10^3 Msun/yr. The total stellar mass is estimated to be no more than 2x10^{11} Msun, which is far less than other luminous radio galaxies at similar redshifts, and suggests that 3C 256 will undergo further star formation or mergers.

We report on the creation of an unbiased catalog of galaxy clusters from the galaxy catalogs derived from the digitized POSS-II (DPOSS). Utilizing the g-r color information, we show that it is possible to estimate redshifts for galaxy clusters at z<0.4 with an rms accuracy of 0.01.

Cosmic gamma-ray bursts (GRBs) have been firmly established as one of the most powerful phenomena in the Universe, releasing electromagnetic energy approaching the rest-mass energy of a neutron star in a few seconds. The two currently popular models for GRB progenitors are the coalescence of two compact objects (such as neutron stars or black holes) or collapse of a massive star. An unavoidable consequence of the latter model is that a bright supernovae should accompany the GRB. The emission from this supernova competes with the much brighter afterglow produced by the relativistic shock that gives rise to the GRB itself. Here we present evidence for an unusual light curve for GRB 980326 based on new optical observations. The transient brightened ~3 weeks after the burst to a flux sixty times larger than that extrapolated from the rapid decay seen at early time. Furthermore, the spectrum changed dramatically and became extremely red. We argue that the new source is the underlying supernova. If our hypothesis is true then this would be the first evidence for a supernova connection with GRBs at cosmological distances. We suggest that GRBs with long durations are associated with death of massive stars.

High-resolution optical spectra of thirteen blue horizontal-branch (BHB) stars in the globular cluster M13 show enormous deviations in element abundances from the expected cluster metallicity. In the hotter stars (T_eff > 12000 K), helium is depleted by factors of 10 to 100 below solar, while iron is enhanced to three times the solar abundance, two orders of magnitude above the canonical metallicity [Fe/H] ~= -1.5 dex for this globular cluster. Nitrogen, phosphorus, and chromium exhibit even more pronounced enhancements, and other metals are also mildly overabundant, with the exception of magnesium, which stays very near the expected cluster metallicity. These photospheric anomalies are most likely due to diffusion --- gravitational settling of helium, and radiative levitation of the other elements --- in the stable radiative atmospheres of these hot stars. The effects of these mechanisms may have some impact on the photometric morphology of the cluster's horizontal branch and on estimates of its age and distance.

Afterglow, or long-lived emission, has now been detected from about a dozen well-positioned gamma-ray bursts. Distance determinations made by measuring optical emission lines from the host galaxy, or absorption lines in the afterglow spectrum, place the burst sources at significant cosmological distances, with redshifts ranging from ~1--3. The energy required to produce the bright gamma-ray flashes is enormous: up to ~10^{53} erg or 10 percent of the rest mass energy of a neutron star, if the emission is isotropic. Here we present the discovery of the optical afterglow and the redshift of GRB 990123, the brightest well-localized GRB to date. With our measured redshift of >1.6, the inferred isotropic energy release exceeds the rest mass of a neutron star thereby challenging current theoretical models for the origin of GRBs. We argue that the optical and IR afterglow measurements reported here may provide the first observational evidence of beaming in a GRB, thereby reducing the required energetics to a level where stellar death models are still tenable.

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.

We report on Hubble Space Telescope (HST) observations of the host galaxy of GRB 971214, about four months after the burst. The redshift of the proposed host galaxy at z=3.418, combined with optical and radio observations of the burst afterglow, implies the extremely large isotropic energy release from the burst in gamma-rays of approximately 3x10^53 ergs, some two orders of magnitude higher than the previously commonly assumed numbers. The positional offset between the optical transient observed at the Keck telescope and the centroid of the proposed host galaxy in the HST image is 0.14 -+0.07 arcsec. We find no evidence in our deep HST image for a chance foreground galaxy superposed along the line of sight to the proposed host at z=3.418. The morphology and photometric properties of this galaxy, such as the total flux, morphology, radial surface profile and scale length, are typical as compared to other, spectroscopically confirmed z>3 galaxies.

We describe DPOSS, a new digital survey of the northern sky, based on the POSS-II photographic sky atlas. The survey covers the entire sky north of delta = -3 deg in 3 bands, calibrated to the Gunn $gri$ system, reaching to equivalent limiting magnitude of B_lim ~ 22 mag. As a result of the state-of-the-art digitisation of the plates, detailed processing of the scans, and a very extensive CCD calibration program, the data quality exceedes that of the previous photographically-based efforts. The end product of the survey will be the Palomar-Norris Sky Catalog, anticipated to contain > 50 million galaxies and > 2 billion stars, down to the survey classification limit, ~ 1 mag above the flux detection limit. Numerous scientific projects utilising these data have been started, and we describe briefly some of them; they illustrate the scientific potential of the data, and serve as the scientific verification tests of the survey. Finally, we discuss some general issues posed by the advent of multi-terabyte data sets in astronomy.

We report on the discovery of the radio, infrared and optical transient coincident with an X-ray transient proposed to be the afterglow of GRB 980703. At later times when the transient has faded below detection, we see an underlying galaxy with R=22.6; this galaxy is the brightest host galaxy (by nearly 2 magnitudes) of any cosmological GRB thus far. In keeping with an established trend, the GRB is not significantly offset from the host galaxy. Interpreting the multi-wavelength data in the framework of the popular fireball model requires that the synchrotron cooling break was between the optical and X-ray bands on July 8.5 UT and that the intrinsic extinction of the transient is Av=0.9. This is somewhat higher than the extinction for the galaxy as a whole, as estimated from spectroscopy.

We present spectroscopic observations of the host galaxy of the gamma-ray burst (GRB) 980703. Several emission and absorption features are detected, making the redshift, z = 0.966, completely unambiguous. This is only the third known redshift for a GRB host. The implied isotropic gamma-ray energy release from the burst is in excess of 10^{53} erg, for a reasonable choice of cosmological parameters. The spectroscopic properties of the host galaxy are typical for a star formation powered object. Using the observed value of the Balmer decrement, we derived the extinction in the galaxy's restframe, A_V = 0.3 +- 0.3 mag. Using three different star formation rate indicators, we estimate SFR ~ 10 Msun/yr, or higher, depending on the extinction, with a lower limit of SFR > 7 Msun/yr. This is the highest value of the star formation rate measured for a GRB galaxy so far, and it gives some support to the idea that GRBs are closely related to massive star formation.

The form of the early-type galaxy scaling relation (the Fundamental Plane or FP) is a direct indicator of the underlying physical origins for the galaxy sequence. Observed properties of the FP include: (1) the slope increases with wavelength; (2) the slope deviates from the virial expectation (assuming homology and constant M/L) at all wavelengths; (3) the intercept evolves passively with redshift; and (4) the slope decreases slowly with redshift. The first property implies that stellar populations contribute to the slope of the FP, the second and fourth properties exclude metallicity effects as the sole cause of the slope, and the third implies that the stellar content of the "average" early-type galaxy formed at high redshift. A composite model--including variations in age and metallicity, as well as a wavelength-independent effect such as homology breaking--is presented which can fit all four observed properties. This model implies that the most luminous early-type galaxies contain the oldest and most metal-rich stars, while the lowest luminosity galaxies formed the bulk of their stars as recently at z_f~1.

We present late-time imaging and spectroscopic observations of the optical transient (OT) of gamma-ray burst (GRB) 970508. Imaging observations roughly 200 and 300 days after the burst provide unambiguous evidence for the flattening of the light-curve. The spectroscopic observations reveal two persistent features which we identify with [O II] 3727 Angstrom and [Ne III] 3869 Angstrom at a redshift of z= 0.835 -- the same redshift as the absorption system seen when the transient was bright. The OT was coincident with the underlying galaxy to better than 370 milliarcsec or a projected radial separation of less than 2.7 kpc. The luminosity of the [O II] line implies a minimum star-formation rate of > 1 M_solar per year. In our assumed cosmology, the implied restframe absolute magnitude is M_B = -18.55, or L_B = 0.12 L_*. This object, the likely host of GRB 970508, can thus be characterized as an actively star-forming dwarf galaxy. The close spatial connection between this dwarf galaxy and the GRB requires that at least some fraction of progenitors be not ejected in even the weakest galactic potentials.

The physical origins of the Fundamental Plane (FP) scaling relations are investigated for early-type galaxies observed at optical and near-infrared wavelengths. The slope for the FP is shown to increase systematically with wavelength from the U-band through the K-band. A distance-independent construction of the observables is described which provides an accurate measurement of the change in the FP slope between any pair of bandpasses. The variation of the FP slope with wavelength is strong evidence of systematic variations in stellar content along the elliptical galaxy sequence. The intercept of the diagnostic relationship between log(D_K/D_V) and log(sigma_0) shows no significant dependence on environment within the uncertainties of the Galactic extinction corrections, demonstrating the universality of the stellar populations contributions at the level of Delta(V-K)=0.03 mag to the zero-point of the global scaling relations. Several other constraints on the properties of early-type galaxies --- the slope of the Mg_2-sigma_0 relation, the effects of stellar populations gradients, and deviations of early-type galaxies from a dynamically homologous family --- are included to construct an empirical, self-consistent model which provides a complete picture of the underlying physical properties which are varying along the early-type galaxy sequence. This empirical approach demonstrates that there are significant systematic variations in both age and metallicity along the elliptical galaxy sequence, and that a small, but systematic, breaking of dynamical homology (or a similar, wavelength independent effect) is required. Predictions for the evolution of the slope of the FP with redshift are described. [abriged]

Near-infrared imaging data on 251 early-type galaxies in clusters and groups are used to construct the near-infrared Fundamental Plane (FP) r_eff ~ sigma_0^1.53 <Sigma_K>_eff^-0.79. The slope of the FP therefore departs from the virial expectation of r_eff ~ sigma_0^2 <Sigma>_eff^-1 at all optical and near-infrared wavelengths, which could be a result of the variation of M/L along the elliptical galaxy sequence, or a systematic breakdown of homology among the family of elliptical galaxies. The slope of the near-infrared FP excludes metallicity variations as the sole cause of the slope of the FP. Age effects, dynamical deviations from a homology, or any combination of these (with or without metallicity), however, are not excluded. The scatter of both the near-infrared and optical FP are nearly identical and substantially larger than the observational uncertainties, demonstrating small but significant intrinsic cosmological scatter for the FP at all wavelengths. The lack of a correlation of the residuals of the near-infrared FP and the residuals from the Mg_2-sigma relation indicates that the thickness of these relations cannot be ascribed only to age or metallicity effects. Due to this metallicity independence, the small scatter of the near-infrared FP excludes a model in which age and metallicity effects ``conspire'' to keep the optical FP thin. All of these results suggest that the possible physical origins of the FP relations are complicated due to combined effects of variations of stellar populations and structural parameters among elliptical galaxies.

Several arguments suggest that quasars at z > 4 may be in the cores of future giant ellipticals, and forming at the very highest peaks of the primordial density field. A strong bias-driven primordial clustering is then expected in these fields, which are naturally interpreted as the cores of future rich clusters. High-redshift quasars can thus be used as markers of some of the earliest galaxy formation sites. Recent discoveries of (proto)galaxy companions of z > 4 quasars, and hints of strong clustering among high-z quasars give some support to this idea. The frequency of companions is much higher than is expected for a general field, or the upper limits for comparable redshifts in the HDF. Clustering of quasars themselves at z > 3 or 4 may reflect the primordial large-scale structure.

An object discovered during an infrared survey of the field near the quasar B2 0149+33 has an emission line at 2.25$\mu$m that we interpret as H$\alpha$ at a redshift of 2.43. The K-band image shows two compact components 10 kpc apart surrounded by more extended emission over ~20 kpc. The H$\alpha$ emission appears to be extended over ~15 kpc (2") in a coarsely sampled (0".8/pixel) image. The star formation rate may be as high as 250 - 1000 M$_\odot$ yr$^{-1}$, depending on the extinction. Alternatively, the line may be powered by an active nucleus, although the probability of serendipitously discovering an AGN in the survey volume is only ~0.02. The increasing number of similar objects reported in the literature indicate that they may be an important, unstudied population in the high redshift universe.

New high quality Keck and ESO images of PKS 1830-211 are presented. Applying a powerful new deconvolution algorithm to these optical and infrared data, both images of the flat spectrum core of the radio source have been identified. An extended source is also detected in the optical images, consistent with the expected location of the lensing galaxy. The source counterparts are very red at I-K~7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

We discuss BIMA (Berkeley Illinois Maryland Association) data and present new high quality optical and near-IR Keck images of the bright radio ring PKS 1830-211. Applying a powerful new deconvolution algorithm we have been able to identify both images of the radio source. In addition we recover an extended source in the optical, consistent with the expected location of the lensing galaxy. The source counterparts are very red, I-K=7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.

Our understanding of the nature of DLA systems at large redshifts, ostensibly progenitors of normal disk galaxies, depends critically on their direct identifications with galaxies, and the resulting measurements of their properties. A few such objects have now been found, reaching out to z ~ 4.1. Their observed luminosities are L ~ L_*, star formation rates a few M_sun/yr, physical sizes ~ 20 kpc, and velocity fields of a few hundred km/s, implying masses > 10^11 M_sun. While their morphology remains uncertain, the observed properties are consistent with those expected of young disk galaxies in the early stages of formation. We also find a statistically significant excess of foreground galaxies near lines of sight to luminous quasars at z > 4. This suggests a systematical gravitational lensing magnification of such quasar samples, possibly with important consequences for the estimates of the quasar luminosity function at high redshifts, and the deduced Omega_b in DLA systems found in their spectra.

The analysis and an efficient scientific exploration of the Digital Palomar Observatory Sky Survey (DPOSS) represents a major technical challenge. The input data set consists of 3 Terabytes of pixel information, and contains a few billion sources. We describe some of the specific scientific problems posed by the data, including searches for distant quasars and clusters of galaxies, and the data-mining techniques we are exploring in addressing them. Machine-assisted discovery methods may become essential for the analysis of such multi-Terabyte data sets. New and future approaches involve unsupervised classification and clustering analysis in the Giga-object data space, including various Bayesian techniques. In addition to the searches for known types of objects in this data base, these techniques may also offer the possibility of discovering previously unknown, rare types of astronomical objects.

The CRONA project is a joint enterprise between Caltech and the astronomical observatory of Rome and Naples aimed to the photometric calibration of the Second Palomar Sky Survey (POSS-II) plates and to the construction of a complete photometric and astrometric catalogue of all objects present on the POSS-II. The calibrated final catalogues of objects will be made available to the community. We present here a comparison between the performances of the SKICAT and DAOPHOT packages applied to POSS-II images of the globular cluster M92.

We have used WFPC2 to construct B, V color-magnitude diagrams of four metal-rich globular clusters, NGC 104 (47 Tuc), NGC 5927, NGC 6388, and NGC 6441. All four clusters have well populated red horizontal branches (RHB), as expected for their metallicity. However, NGC 6388 and 6441 also exhibit a prominent blue HB (BHB) extension, including stars reaching as faint in V as the turnoff luminosity. This discovery demonstrates directly for the first time that a major population of hot HB stars can exist in old, metal-rich systems. This may have important implications for the interpretation of the integrated spectra of elliptical galaxies. The cause of the phenomenon remains uncertain. We examine the possibility that NGC 6388 and 6441 are older than the other clusters, but a simple difference in age may not be sufficient to produce the observed distributions along the HB. The high central densities in NGC 6388 and 6441 suggest that the existence of the blue HB (BHB) tails might be caused by stellar interactions in the dense cores of these clusters, which we calculate to have two of the highest collision rates among globular clusters in the Galaxy. Tidal collisions might act in various ways to enhance loss of envelope mass, and therefore populate the blue side of the HB. However, the relative frequency of tidal collisions does not seem large enough (compared to that of the clusters with pure RHBs) to account for such a drastic difference in HB morphology. While a combination of an age difference and dynamical interactions may help, prima facie the lack of a radial gradient in the BHB/RHB star ratio seems to argue against dynamical effects playing a role.

We present distributions of colors of stars along the horizontal branch of the globular cluster NGC 2808, from Hubble Space Telescope WFPC2 imaging in B, V, and an ultraviolet filter (F218W). This cluster's HB is already known to be strongly bimodal, with approximately equal-sized HB populations widely separated in the color-magnitude diagram. Our images reveal a long blue tail with two gaps, for a total of four nearly distinct HB groups. These gaps are very narrow, corresponding to envelope-mass differences of only \sim 0.01 Msun. This remarkable multimodality may be a signature of mass-loss processes, subtle composition variations, or dynamical effects; we briefly summarize the possibilities. The existence of narrow gaps between distinct clumps on the HB presents a challenge for models that attempt to explain HB bimodality or other peculiar HB structures.

In this paper we present preliminary results from our HST project aimed at exploring the connection between the stellar dynamics and stellar evolution in the cores of high-density globular clusters. We illustrate the discovery of horizontal branches (HB) with extended blue tails in the very metal rich globular clusters NGC 6388 ([Fe/H]=-0.6) and NGC 6441 ([Fe/H]=-0.5), and of a four-mode horizontal branch in NGC 2808. We briefly summarize the possible causes of these anomalies, which might be a signature of mass-loss processes,subtle composition variations, or dynamical effects. The existence of blue HBs in metal rich clusters and of significant narrow gaps between clumps on the HB are a challenge for the theoretical models. The composite color-magnitude diagram and luminosity function of 878 blue stragglers in our sample of 13 galactic globular clusters are also presented and briefly discussed.

In this paper we present the 13 B and V color-magnitude diagrams discussed by Piotto et al. in the same Procedings. The data have been collected with the WFPC2 camera on HST. We emphasize the discovery of blue horizontal branches in the two metal rich globular clusters NGC 6388 and NGC 6441, and the four-mode HB of NGC 2808. The significance of the gaps in the HB of NGC 2808 is also assesed.

We are conducting an effort to catalog all sources detected in the digitized POSS-II (DPOSS). This is now becoming an international collaboration including Caltech and the Observatories of Rome and Naples (project CRONA). There are also ongoing extensive CCD calibration efforts at Palomar. The resulting Palomar-Norris Sky Catalog (PNSC) is expected to contain > 5 x 10^7 galaxies, and > 2 x 10^9 stars, in 3 colors (photographic JFN bands, calibrated to CCD gri system), down to the limiting magnitude equivalent of B approx 22 mag. The star-galaxy classification is accurate to 90 - 95% down to the equivalent of B approx 21 mag. The catalog will be made available to the community via computer networks or other suitable media, probably in installments, as soon as scientific validation and quality checks are completed and the funding allows it. Analysis and manipulation software will also be freely available. A great variety of scientific projects will be possible with this vast new data base, including studies of the large-scale structure in the universe, the Galactic structure, automatic optical identifications of sources from radio through x-ray bands, generation of objectively defined catalogs of clusters and groups of galaxies, generation of statistically complete catalogs of galaxies to be used in future redshift surveys, searches for high-redshift quasars and other active objects, searches for variable or extreme-color objects, etc. We have already started work on a number of such projects, and we summarize briefly some of the results obtained to date.

We review the concept of the Fundamental Plane (FP), and present new results on the near-infrared FP. We show that the K-band FP differs both from the optical form and the virial expectation under the assumption of constant (M/L) and homology. Systematic variations of only the stellar populations parameters (age, metallicity, IMF) cannot reproduce the slopes of the FP simultaneously at both the optical and near-infrared wavelengths. There appears to be an additional effect which could be due to a systematic departure of the velocity distributions of elliptical galaxies from a homologous family. In order to distinguish between the different possible stellar populations parameters, the FP and its projections can be observed at a range of redshifts. We describe several such evolutionary tests: in the intercept of the color-magnitude relation; in the intercept of the Kormendy surface brightness-radius correlation; and in the intercept of the FP itself. All three tests are fully consistent with the population of cluster elliptical galaxies having formed at high redshift; this contradicts large, systematic age variations as a significant contributor to the slope of the FP.

We report the detection, both in line and continuum emission, of a galaxy responsible for a damped Ly-alpha absorption (DLA) system at z=3.15 projected 2.3 arcsec from QSO 2233+131. The star formation rate implied by both the Ly-alpha and restframe far-UV continuum is approximately 7 M_sun / yr. The galaxy also shows no sign of an active nucleus and has blue colors, both consistent with it being observed early in its star formation history. The Ly-alpha emission is offset approximately 200 km/s from the DLA velocity, which could be explained as the rotational signature of a disk galaxy. We also report on the progress of an ongoing near-infrared narrow-band search with the Keck telescope for line emission from a population of primeval galaxies at high redshift. We have found several promising candidates which are being followed up with optical spectroscopy. Our survey is beginning to place constraints on a possible population of slightly dusty primeval galaxies at high redshift.

Color Magnitude Diagrams (CMDs) of individual stars in 8 Globular Clusters in M31 down to about 1 mag fainter (V$\sim 26.5$) than the Horizontal Branch have been obtained with the Hubble Space Telescope. In particular, we observed G280 and G351 with the FOC (f/96+F430W, f/96+F480LP) while the WFPC2 (F555W,F814W) frames for G1,G58,G105,G108,G219+Bo468 were retrieved from the HST archive. The cluster metallicities -\fe- range from -1.8 to -0.4. Coupled with sufficiently accurate (to $\sim\pm0.1$ mag) measures of the mean brightness of the HB --\vhb--, appropriate estimates of reddening for each cluster, and the adoption of a distance modulus to M31 of \dmo = 24.43, this has allowed us to yield a direct calibration for the mean absolute magnitude of the HB at the instability strip --\mv-- with varying \fe: M_V^{HB} = (0.13 \pm 0.07)\fe + (0.95 \pm 0.09), where the associated errors result from the adopted global errors in the measure and best fitting procedures. The slope of the derived relation is fully consistent with that predicted by the standard and canonical models ($\sim0.15$) and obtained by various ground-based observations, while it is only marginally compatible with higher values ($\sim0.30$), also obtained in the past. The zero-point, which is crucial to absolute age determinations, depends on the adopted distance to M31 and is moreover affected by an additional error due to the residual uncertainties in the \hst photometric zero-points ($\sim0.05$ mag, at least). If confirmed, such a calibration of the \mv \vs \fe relationship would imply old absolute ages ($> 16$Gyr) for the oldest Galactic globulars and fairly small age spread among those having a constant magnitude difference between the Main-Sequence Turnoff and the HB.

We report on imaging and spectroscopic observations of the radio-loud, highly polarized quasar PKS 1610-771 (z = 1.71). Our long-slit spectroscopy of the companion 4.55 arcseconds NW of the quasar confirms the stellar nature of this object, so ruling out the previously suspected gravitationally lensed nature of this system. PKS 1610-771 looks fuzzy on our sub-arcsecond R and I images and appears located in a rich environment of faint galaxies. Possible magnification, without image splitting of the quasar itself, by some of these maybe foreground galaxies cannot be excluded. The continuum fuzz (made of the closest two objects, viz. A and D) is elongated in a direction orthogonal to the E vector of the optical polarization, as in high-redshift radio-galaxies. The spectrum of PKS 1610-771 appears strongly curved, in a convex way, with a maximum of intensity at ~ 7,600 A (2,800 A rest frame), possibly indicating a strong ultraviolet absorption by dust.

Galaxies believed to be responsible for damped Ly-alpha absorption (DLA) systems in the spectra of high-redshift quasars represent a viable population of progenitors of normal disk galaxies. They appear to contain a substantial fraction of the baryons known to exist in normal galaxies today. Here we report on the detection of an object, designated DLA 2233+131, responsible for a previously known DLA system at z(abs)=3.150 in the spectrum of a quasar 2233+131 [z(QSO)=3.295]. This is the first unambiguous detection of a DLA galaxy, in both emission line and stellar continuum. Its properties correspond closely to what may be expected from a young disk galaxy in the early stages of formation, with no sign of an active nucleus. This gives a strong support to the idea that DLA systems represent a population of young galaxies at high redshifts.

We use the intercept of the elliptical galaxy radius--surface brightness (SB) relation at a fixed metric radius as the standard condition for the Tolman SB test of the universal expansion. We use surface photometry in the optical and near-IR of elliptical galaxies in Abell~2390 ($z=0.23$) and Abell~851 ($z=0.41$), and compare them to the Coma cluster at $z\approx 0$. The photometric data for each cluster are well-described by the Kormendy relation $r_e \propto \Sigma_e^{A}$, where $A=-0.9$ in the optical and $A=-1.0$ in the near-IR. The scatter about this near-IR relation is only $0.076$ in $\log r_e$ at the highest redshift, which is much smaller than at low redshifts, suggesting a remarkable homogeneity of the cluster elliptical population at $z=0.41$. We use the intercept of these fixed-slope correlations at $R_e = 1$~kpc (assuming $H_0=75$~km~s$^{-1}$~Mpc$^{-1}$, $\Omega_0=0.2$, and $\Lambda_0=0$, where the results are only weakly dependent on the cosmology) to construct the Tolman SB test for these three clusters. The data are fully consistent with universal expansion if we assume simple models of passive evolution for elliptical galaxies, but are inconsistent with a non-expanding geometry (the tired light cosmology) at the $5 \, \sigma$ confidence level at $z=0.41$. These results suggest luminosity evolution in the restframe $K$-band of $0.36 \pm 0.14$~mag from $z = 0.41$ to the present, and are consistent with the ellipticals having formed at high redshift. The SB intercept in elliptical galaxy correlations is thus a powerful tool for investigating models of their evolution for significant lookback times.

We present results from a near-infrared $K$-band imaging survey of 59 elliptical galaxies in five nearby clusters. We measure photometric parameters for each galaxy using surface photometry and draw velocity dispersions from the literature. Three observables define a near-infrared Fundamental Plane (FP) of elliptical galaxies with $R_e \propto \sigma^{1.44\pm0.04} \Sigma_e^{-0.79\pm0.04}$. The scatter in the near-infrared relation is small at $16.5$\% in distance, which is equivalent to, or less than, the scatter of the optical FP. We suggest that the small deviation of the near-infrared FP relation from the optical FP is due to the reduction of metallicity effects in the near-infrared bandpass. While the small scatter of the optical FP could be consistent with compensating effects of age and metallicity, the similarly small scatter of the near-infrared FP is nearly independent of metallicity and hence places a strong constraint on possible age spreads among elliptical galaxies at every point along the FP. We suggest that the departure of the near-infrared FP from the pure virial form $R_e \propto \sigma^2 \Sigma_e^{-1}$, and the corresponding observed relation $\left( M/L \right) \propto M^{0.16\pm 0.01}$, may be explained by slight systematic departures of the structure and dynamics of elliptical galaxies from a homology.

We present the first results from a near-infrared narrow-band search for primeval galaxies (PGs) using the Keck 10~m telescope. We have targeted our fields for three QSOs and one radio galaxy at redshifts $z = 2.28$ to $4.70$. We selected narrow-band filters in the $K$-band centered on strong emision lines at the redshift of the targets. We reach limiting fluxes between $1.6\times 10^{-17}$~erg~s$^{-1}$~cm$^{-2}$ and $1.0\times 10^{-16}$~erg~s$^{-1}$~cm$^{-2}$, and cover a total area of $\sim 4$~arcmin$^2$. At the redshifts of interest, these flux limits correspond to typical restframe line luminosities of $\sim 10^{42} - 10^{43}$~erg~s$^{-1}$, unobscured star formation rates of $\sim 1 - 100$ $M_\odot$~yr$^{-1}$, and a sampled comoving volume of several tens of Mpc$^3$. We have found no candidate PGs at a $2\sigma$ confidence level. We demonstrate that for moderate dust absorption ($A_V \gsim 2^m$) in a simple dust-screen model, our preliminary survey puts a strong constraint on intrinsic PG luminosities during an intial burst of star formation. In the case of the QSOs, we have used PSF-subtraction to search for the presence of faint, extended line emission surrounding these radio-quiet AGNs. We find no evidence for extended emission down to $3.7 \times 10^{-17}$~erg~s$^{-1}$~cm$^{-2}$. This imposes limits on the reprocessed fraction of the QSO continuum emission ranging between $0.02$\% and $0.29$\%; if beaming effects are not important, then this implies similarly small covering factors for clouds optically thick to the QSO continuum emission.

We report the discovery with the Keck telescope of two new multiply-imaged arcs in the luminous X-ray cluster A2219 ($z=0.225$). The brighter arc in the field is red and we use spectroscopic and photometric information to identify it as a $z \sim 1$ moderately star-forming system. The brightness of this arc implies that it is formed from two merging images of the background source, and we identify possible candidates for the third image of this source. The second giant arc in this cluster is blue, and while fainter than the red arc it has a similarly large angular extent (32 arcsec). This arc comprises three images of a single nucleated source -- the relative parities of the three images are discernible in our best resolution images. The presence of several bright multiply imaged arcs in a single cluster allows detailed modelling of the cluster mass distribution, especially when redshift information is available. We present a lensing model of the cluster which explains the properties of the various arcs, and we contrast this model with the optical and X-ray information available on the cluster. We uncover significant differences between the distributions of mass and X-ray gas in the cluster. We suggest that such discrepancies may indicate an on-going merger event in the cluster core, possibly associated with a group around the second brightest cluster member. The preponderance of similar merger signatures in a large fraction of the moderate redshift clusters would indicate their dynamical immaturity.

We present the initial results of a study of globular clusters in M31, using the Faint Object Camera (FOC) on the Hubble Space Telescope (HST). The sample of objects consists of 13 clusters spanning a range of properties. Three inde- pendent image deconvolution techniques were used in order to compensate for the optical problems of the HST, leading to mutually fully consistent results. We present detailed tests and comparisons to determine the reliability and limits of these deconvolution methods, and conclude that high-quality surface photome- try of M31 globulars is possible with the HST data. Surface brightness profiles have been extracted, and core radii, half-light radii, and central surface brightness values have been measured for all of the clusters in the sample. Their comparison with the values derived from ground-based observations indi- cates the later to be systematically and strongly biased by the seeing effects, as it may be expected. A comparison of the structural parameters with those of the Galactic globulars shows that the structural properties of the M31 globu- lars are very similar to those of their Galactic counterparts. A candidate for a post-core-collapse cluster, Bo 343 = G 105, has been already identified from these data; this is the first such detection in the M31 g. c. system.