CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

CASA, the Common Astronomy Software Applications, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), and is often used also for other radio telescopes. CASA has always been distributed as a single, integrated application, including a Python interpreter and all the libraries, packages and modules. As part of the ongoing development of CASA 6, and the switch from Python 2 to 3, CASA will provide greater flexibility for users to integrate CASA into existing Python workflows by using a modular architecture and standard pip wheel installation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference will give an overview of the CASA 6 project.

We present the results of a survey for damped (DLA, log N(H I) > 20.3) and sub-damped Lyman-? systems (19.5 < log N(H I) < 20.3) at z > 2.55 along the lines-of-sight to 77 quasars with emission redshifts in the range 4 < zem < 6.3. Intermediate resolution (R ? 4300) spectra have been obtained with the Echellette Spectrograph and Imager (ESI) mounted on the Keck telescope. A total of 100 systems with log N(H I) > 19.5 are detected of which 40 systems are damped Lyman-? systems for an absorption length of ?X = 378. About half of the lines of sight of this homogeneous survey have never been investigated for DLAs. We study the evolution with redshift of the cosmological density of the neutral gas and find, consis- tently with previous studies at similar resolution, that ?DLA,H I decreases at z > 3.5. The overall cosmological evolution of ?HI shows a peak around this redshift. The H I column density distribution for log N(H I) ? 20.3 is ?tted, consistently with previous surveys, with a single power-law of index ? ? -1.8$\pm$0.25. This power-law overpredicts data at the high-end and a second, much steeper, power-law (or a gamma function) is needed. There is a flattening of the function at lower H I column densities with an index of ? ? ?1.4 for the column density range log N(H I) = 19.5?21. The fraction of H I mass in sub-DLAs is of the order of 30%. The H column density distribution does not evolve strongly from z ? 2.5 to z ? 4.5.

We study the density field around z_em > 4 quasars using high quality medium spectral resolution ESI-Keck spectra (R~4300, SNR > 25) of 45 high-redshift quasars selected from a total of 95 spectra. This large sample considerably increases the statistics compared to previous studies. The redshift evolution of the mean photo-ionization rate and the median optical depth of the intergalactic medium (IGM) are derived statistically from the observed transmitted flux and the pixel optical depth probability distribution function respectively. This is used to study the so-called proximity effect, that is, the observed decrease of the median optical depth of the IGM in the vicinity of the quasar caused by enhanced photo-ionization rate due to photons emitted by the quasar. We show that the proximity effect is correlated with the luminosity of the quasars, as expected. By comparing the observed decrease of the median optical depth with the theoretical expectation we find that the optical depth does not decrease as rapidly as expected when approaching the quasar if the gas in its vicinity is part of the standard IGM. We interpret this effect as revealing gaseous overdensities on scales as large as ~15 Mpc/h. The mean overdensity is of the order of two and five within, respectively, 10 and 3 Mpc/h. If true, this would indicate that high redshift quasars are located in the center of overdense regions that could evolve with time into massive clusters of galaxies. The overdensity is correlated with luminosity: brighter quasars show higher overdensities.

We report the abundance analysis of new high S/N spectra of the most metal-poor ([Fe/H] $= -2.95$) star presently known to be a member of a dwarf galaxy, the Draco dSph red giant, D119. No absorption lines for elements heavier than Ni are detected in two Keck HIRES spectra covering the $\lambda\lambda$ 3850--6655 \AA{} wavelength range, phenomenon not previously noted in any other metal-poor star. We present upper limits for several heavy element abundances. The most stringent limits, based on the non-detection of \ion{Sr}{2} and \ion{Ba}{2} lines, indicate that the total s- and r-process enrichment of D119 is at least 100 times smaller than Galactic stars of similar metallicity. The light element abundances are consistent with the star having formed out of material enciched primarily by massive Type II supernovae (M $> 20$--25 M$_{\odot}$). If this is the case, we are forced to conclude that massive, metal-poor Type II supernovae did not contribute to the r-process in the proto-Draco environment. We compare the abundance pattern observed in D119 to current predictions of prompt enrichement and pair-instability supernovae and find that the model predictions fail by an order or maginitude or more for many elements.

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.

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.

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.

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.

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 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 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.

Using the NIRSPEC spectrograph at Keck II, we have obtained echelle spectra covering the range 1.5-1.8um for 2 of the brightest giants in Liller 1 and NGC 6553, old metal rich globular clusters in the Galactic bulge. We use spectrum synthesis for the abundance analysis, and find [Fe/H]=-0.3 +/- 0.2 and [O/H]=+0.3 +/- 0.2 dex. The composition of the clusters is similar to that of field stars in the bulge and is consistent with a sceanrio in which the clusters formed early, with rapid enrichment. We have dificulty achieveing a good fit to the spectrum of NGC 6553 using either the low or the high values recently reported in the literature, unless unusually large, or no alpha-element enhancements are adopted, respectively.

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.

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.

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.

We present two Keck HIRES spectra (R=40,000) of a bulge K3 giant taken on successive nights during the second caustic crossing of the binary microlensing event EROS BLG-2000-5. This caustic crossing served to effectively resolve the surface of the source star: the spectrum from the second night is dominated by the limb, while the spectrum from the first night is comprised of light from a broader range of radii. To demonstrate that the spectra are adequate to resolve the differences between them, we analyze the H-alpha line. The equivalent width is ~8% smaller on the second night, and the signal-to-noise ratio per resolution element (165 and 75 on the two nights respectively) is sufficient to show that the difference is approximately constant over the ~2A (~15 resolution element) extent of the line. The sign of the difference is in the expected direction since the limb is the coolest part of the star and therefore should have the weakest H-alpha. We invite atmosphere modelers to predict the difference spectrum from the entire spectral range 5500<lambda<7900A so that these predictions can be compared to our observations.

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.

We present Cu and Ba abundances for 7 G-K dwarf stars, members of the solar-metallicity, 0.3 Gyr old Ursa Major Moving Group. All analyzed member stars show [Ba/Fe] excesses of +0.3-plus, associated with [Cu/Fe] deficiencies of up to -0.23 dex. The present results suggest that there is an anti-correlation between the abundances of Cu and the heavy elements produced by the main component of the neutron capture s-process. Other possible anomalies are Na and C deficiencies with respect to normal solar-metallicity stars. The new data do not confirm the recent claim that the group member HR6094 is a Ba dwarf star.

We report detailed analysis of high-resolution spectra of nine high velocity metal-rich dwarfs in the solar neighborhood. The stars are super metal-rich and 5 of them have [Fe/H]>=+0.4, making them the most metal-rich stars currently known. We find that alpha-elements decrease with increasing metallicity; s-elements are underabundant by about [s-elements/Fe]=-0.3. While exceeding the [Fe/H] of current bulge samples, the chemistry of these stars has important similarities and differences. The near-solar abundances of the alpha-capture elements places these stars on the metal-rich extension of McWilliam & Rich (1994) [ApJS, 91, 749], but their s-process abundances are much lower than those of the bulge giants. These low s-process values have been interpreted as the hallmark of an ancient stellar population.