We report studies on mitigation of optical effects of bright Low-Earth-Orbit (LEO) satellites on Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST). These include options for pointing the telescope to avoid satellites, laboratory investigations of bright trails on the Rubin Observatory LSST Camera sensors, algorithms for correcting image artifacts caused by bright trails, experiments on darkening SpaceX Starlink satellites, and ground-based follow-up observations. Starlink satellites with no brightness mitigation are presently g ~ 5.1 mag, and an initial experiment "DarkSat" is g ~ 6.1 mag. Future Starlink darkening plans may reach g ~ 7 mag, a brightness level that enables non-linear image artifact correction to the same level as background noise. However, the satellite trails will still exist at S/N ~ 100, generating systematics that may impact data analysis and limiting some science. LEO satellite trails are wider than a point-spread function because satellites are slightly out of focus due to their finite distance; for Rubin Observatory's 8.4-m mirror and a satellite at 550 km, this widening effect is about 3 arcsec, which helps avoid saturation by decreasing the trail's peak surface brightness.

The LSST survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive weak lensing shear systematics on a baseline operational simulation, using the $\rho-$statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to $\sim 3-4$ better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to weak lensing) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the weak lensing statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

Accurate photometric redshift (photo-$z$) estimates are essential to the cosmological science goals of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this work we use simulated photometry for mock galaxy catalogs to explore how LSST photo-$z$ estimates can be improved by the addition of near-infrared (NIR) and/or ultraviolet (UV) photometry from the Euclid, WFIRST, and/or CASTOR space telescopes. Generally, we find that deeper optical photometry can reduce the standard deviation of the photo-$z$ estimates more than adding NIR or UV filters, but that additional filters are the only way to significantly lower the fraction of galaxies with catastrophically under- or over-estimated photo-$z$. For Euclid, we find that the addition of ${JH}$ $5{\sigma}$ photometric detections can reduce the standard deviation for galaxies with $z>1$ ($z>0.3$) by ${\sim}20\%$ (${\sim}10\%$), and the fraction of outliers by ${\sim}40\%$ (${\sim}25\%$). For WFIRST, we show how the addition of deep ${YJHK}$ photometry could reduce the standard deviation by ${\gtrsim}50\%$ at $z>1.5$ and drastically reduce the fraction of outliers to just ${\sim}2\%$ overall. For CASTOR, we find that the addition of its ${UV}$ and $u$-band photometry could reduce the standard deviation by ${\sim}30\%$ and the fraction of outliers by ${\sim}50\%$ for galaxies with $z<0.5$. We also evaluate the photo-$z$ results within sky areas that overlap with both the NIR and UV surveys, and when spectroscopic training sets built from the surveys' small-area deep fields are used.

In this white paper, we present the scientific cases for adding narrowband optical filters to the Large Synoptic Survey Telescope (LSST). LSST is currently planning to observe the southern sky in 6 broadband optical filters. Three of the four LSST science themes would benefit from adding narrowband filter observations. We discuss the technical considerations of using narrowband filters with the LSST and lay out the scientific impact that would result on the study of AGB stars, emission line nebula (e.g., supernova remnants and planetary nebulae), photometric redshifts of galaxies, and the determination of stellar parameters.

The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis.

Telescope scheduling is the task of determining the best sequence of observations (pointings and filter choices) for a survey system. Because it is computationally intractable to optimize over all possible multi-year sequences of observations, schedulers use heuristics to pick the best observation at a given time. A greedy scheduler selects the next observation by choosing whichever one maximizes a scalar merit function, which serves as a proxy for the scientific goals of the telescope. This sort of bottom-up approach for scheduling is not guaranteed to produce a schedule for which the sum of merit over all observations is maximized. As an alternative to greedy schedulers, we introduce ALTSched, which takes a top-down approach to scheduling. Instead of considering only the next observation, ALTSched makes global decisions about which area of sky and which filter to observe in, and then refines these decisions into a sequence of observations taken along the meridian to maximize SNR. We implement ALTSched for the Large Synoptic Survey Telescope (LSST), and show that it equals or outperforms the baseline greedy scheduler in essentially all quantitative performance metrics. Due to its simplicity, our implementation is considerably faster than OpSim, the simulated greedy scheduler currently used by the LSST Project: a full ten year survey can be simulated in 4 minutes, as opposed to tens of hours for OpSim. LSST's hardware is fixed, so improving the scheduling algorithm is one of the only remaining ways to optimize LSST's performance. We see ALTSched as a prototype scheduler that gives a lower bound on the performance achievable by LSST.

The Large Synoptic Survey Telescope (LSST) project will conduct a ten year multi-band survey starting in 2022. Observing strategies for this survey are being actively investigated, and the science capabilities can be best forecasted on the basis of simulated strategies from the LSST Operations Simulator (OpSim). OpSim simulates a stochastic realization of the sequence of LSST pointings over the survey duration, and is based on a model of the observatory (including telescope) and historical data of observational conditions. OpSim outputs contain a record of each simulated pointing of the survey along with a complete characterization of the pointing in terms of observing conditions, and some useful quantities derived from the characteristics of the pointing. Thus, each record can be efficiently used to derive the properties of observations of all astrophysical sources found in that pointing. However, in order to obtain the time series of observations (light curves) of a set of sources, it is often more convenient to compute all observations of an astrophysical source, and iterate over sources. In this document, we describe the open source python package OpSimSummary which allows for a convenient reordering. The objectives of this package are to provide users with an Application Programming Interface (API) for accessing all such observations and summarizing this information in the form intermediate data products usable by third party software such as SNANA, thereby also bridging the gap between official LSST products and pre-existing simulation codes.

We report the discovery of an SN1988Z-like type IIn supernova KISS15s found in a low-mass star-forming galaxy at redshift z=0.038 during the course of the Kiso Supernova Survey (KISS). KISS15s shows long-duration optical continuum and emission line light curves, indicating that KISS15s is powered by a continuous interaction between the expanding ejecta and dense circumstellar medium (CSM). The H$\alpha$ emission line profile can be decomposed into four Gaussians of narrow, intermediate, blue-shifted intermediate, and broad velocity width components, with a full width at half maximum of $\lesssim 100$, $\sim 2,000$, and $\sim 14,000$ km s${}^{-1}$ for the narrow, intermediate, and broad components, respectively. The presence of the blue-shifted intermediate component, of which the line-of-sight velocity relative to the systemic velocity is about $-5,000$ km s${}^{-1}$, suggests that the ejecta-CSM interaction region has an inhomogeneous morphology and anisotropic expansion velocity. We found that KISS15s shows increasing infrared continuum emission, which can be interpreted as hot dust thermal emission of $T \sim 1,200$ K from newly formed dust in a cool, dense shell in the ejecta-CSM interaction region. The progenitor mass-loss rate, inferred from bolometric luminosity, is $\dot{M} \sim 0.4 M_{\odot} \text{yr}^{-1} (v_{w}/40 \text{km}~\text{s}^{-1})$, where $v_{w}$ is the progenitor's stellar wind velocity. This implies that the progenitor of KISS15s was a red supergiant star or a luminous blue variable that had experienced a large mass-loss in the centuries before the explosion.

Cosmology is one of the four science pillars of LSST, which promises to be transformative for our understanding of dark energy and dark matter. The LSST Dark Energy Science Collaboration (DESC) has been tasked with deriving constraints on cosmological parameters from LSST data. Each of the cosmological probes for LSST is heavily impacted by the choice of observing strategy. This white paper is written by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy that will benefit all cosmological analyses with LSST. It is accompanied by the DESC DDF (Deep Drilling Fields) white paper (Scolnic et al.). We use a variety of metrics to understand the effects of the observing strategy on measurements of weak lensing, large-scale structure, clusters, photometric redshifts, supernovae, strong lensing and kilonovae. In order to reduce systematic uncertainties, we conclude that the current baseline observing strategy needs to be significantly modified to result in the best possible cosmological constraints. We provide some key recommendations: moving the WFD (Wide-Fast-Deep) footprint to avoid regions of high extinction, taking visit pairs in different filters, changing the 2x15s snaps to a single exposure to improve efficiency, focusing on strategies that reduce long gaps (>15 days) between observations, and prioritizing spatial uniformity at several intervals during the 10-year survey.

We review the measurements of dark energy enabled by observations of the Deep Drilling Fields and the optimization of survey design for cosmological measurements. This white paper is the result of efforts by the LSST DESC Observing Strategy Task Force (OSTF), which represents the entire collaboration, and aims to make recommendations on observing strategy for the DDFs that will benefit all cosmological analyses with LSST. It is accompanied by the DESC-WFD white paper (Lochner et al.). We argue for altering the nominal deep drilling plan to have $>6$ month seasons, interweaving $gri$ and $zy$ observations every 3 days with 2, 4, 8, 25, 4 visits in $grizy$, respectively. These recommendations are guided by metrics optimizing constraints on dark energy and mitigation of systematic uncertainties, including specific requirements on total number of visits after Y1 and Y10 for photometric redshifts (photo-$z$) and weak lensing systematics. We specify the precise locations for the previously-chosen LSST deep fields (ELAIS-S1, XMM-LSS, CDF-S, and COSMOS) and recommend Akari Deep Field South as the planned fifth deep field in order to synergize with Euclid and WFIRST. Our recommended DDF strategy uses $6.2\%$ of the LSST survey time. We briefly discuss synergy with white papers from other collaborations, as well as additional mini-surveys and Target-of-Opportunity programs that lead to better measurements of dark energy.

How ground-based telescopes schedule their observations in response to competing science priorities and constraints, variations in the weather, and the visibility of a particular part of the sky can significantly impact their efficiency. In this paper we introduce the Feature-Based telescope scheduler that is an automated, proposal-free decision making algorithm that offers \textit{controllability} of the behavior, \textit{adjustability} of the mission, and quick \textit{recoverability} from interruptions for large ground-based telescopes. By framing this scheduler in the context of a coherent mathematical model the functionality and performance of the algorithm is simple to interpret and adapt to a broad range of astronomical applications. This paper presents a generic version of the Feature-Based scheduler, with minimal manual tailoring, to demonstrate its potential and flexibility as a foundation for large ground-based telescope schedulers which can later be adjusted for other instruments. In addition, a modified version of the Feature-Based scheduler for the Large Synoptic Survey Telescope (LSST) is introduced and compared to previous LSST scheduler simulations.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We use SDSS imaging data in $ugriz$ passbands to study the shape of the point spread function (PSF) profile and the variation of its width with wavelength and time. We find that the PSF profile is well described by theoretical predictions based on von K\'{a}rm\'{a}n's turbulence theory. The observed PSF radial profile can be parametrized by only two parameters, the profile's full width at half maximum (FWHM) and a normalization of the contribution of an empirically determined "instrumental" PSF. The profile shape is very similar to the "double gaussian plus power-law wing" decomposition used by SDSS image processing pipeline, but here it is successfully modeled with two free model parameters, rather than six as in SDSS pipeline. The FWHM variation with wavelength follows the $\lambda^{\alpha}$ power law, where $\alpha \approx-0.3$ and is correlated with the FWHM itself. The observed behavior is much better described by von K\'{a}rm\'{a}n's turbulence theory with the outer scale parameter in the range 5$-$100 m, than by the Kolmogorov's turbulence theory. We also measure the temporal and angular structure functions for FWHM and compare them to simulations and results from literature. The angular structure function saturates at scales beyond 0.5$-$1.0 degree. The power spectrum of the temporal behavior is found to be broadly consistent with a damped random walk model with characteristic timescale in the range $\sim5-30$ minutes, though data show a shallower high-frequency behavior. The latter is well fit by a single power law with index in the range $-1.5$ to $-1.0$. A hybrid model is likely needed to fully capture both the low-frequency and high-frequency behavior of the temporal variations of atmospheric seeing.

In this paper we present and characterize a nearest-neighbors color-matching photometric redshift estimator that features a direct relationship between the precision and accuracy of the input magnitudes and the output photometric redshifts. This aspect makes our estimator an ideal tool for evaluating the impact of changes to LSST survey parameters that affect the measurement errors of the photometry, which is the main motivation of our work (i.e., it is not intended to provide the "best" photometric redshifts for LSST data). We show how the photometric redshifts will improve with time over the 10-year LSST survey and confirm that the nominal distribution of visits per filter provides the most accurate photo-$z$ results. The LSST survey strategy naturally produces observations over a range of airmass, which offers the opportunity of using an SED- and $z$-dependent atmospheric affect on the observed photometry as a color-independent redshift indicator. We show that measuring this airmass effect and including it as a prior has the potential to improve the photometric redshifts and can ameliorate extreme outliers, but that it will only be adequately measured for the brightest galaxies, which limits its overall impact on LSST photometric redshifts. We furthermore demonstrate how this airmass effect can induce a bias in the photo-$z$ results, and caution against survey strategies that prioritize high-airmass observations for the purpose of improving this prior. Ultimately, we intend for this work to serve as a guide for the expectations and preparations of the LSST science community with regards to the minimum quality of photo-$z$ as the survey progresses.

Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

We present a study of the kinematics of the extraplanar ionized gas around several dozen galaxies observed by the Mapping of Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We considered a sample of 67 edge-on galaxies out of more than 1400 extragalactic targets observed by MaNGA, in which we found 25 galaxies (or 37%) with regular lagging of the rotation curve at large distances from the galactic midplane. We model the observed $H\alpha$ emission velocity fields in the galaxies, taking projection effects and a simple model for the dust extinction into the account. We show that the vertical lag of the rotation curve is necessary in the modeling, and estimate the lag amplitude in the galaxies. We find no correlation between the lag and the star formation rate in the galaxies. At the same time, we report a correlation between the lag and the galactic stellar mass, central stellar velocity dispersion, and axial ratio of the light distribution. These correlations suggest a possible higher ratio of infalling-to-local gas in early-type disk galaxies or a connection between lags and the possible presence of hot gaseous halos, which may be more prevalent in more massive galaxies. These results again demonstrate that observations of extraplanar gas can serve as a potential probe for accretion of gas.

In this work we quantify the performance of $LSST$ on the detection of eclipsing binaries. We use $Kepler$ observed binaries to create a large sample of simulated pseudo-$LSST$ binary light curves. From these light curves, we attempt to recover the known binary signal. The success rate of period recovery from the pseudo-$LSST$ light curves is indicative of $LSST$'s expected performance. Using an off-the-shelf Analysis of Variance (AoV) routine, we successfully recover 71% of the targets in our sample. We examine how the binary period impacts recovery success and see that for periods longer than 10~days the chance of successful binary recovery drops below 50%.

The Large Synoptic Survey Telescope (LSST) will survey the southern sky from 2022--2032 with unprecedented detail. Since the observing strategy can lead to artifacts in the data, we investigate the effects of telescope-pointing offsets (called dithers) on the $r$-band coadded 5$\sigma$ depth yielded after the 10-year survey. We analyze this survey depth for several geometric patterns of dithers (e.g., random, hexagonal lattice, spiral) with amplitude as large as the radius of the LSST field-of-view, implemented on different timescales (per season, per night, per visit). Our results illustrate that per night and per visit dither assignments are more effective than per season. Also, we find that some dither geometries (e.g., hexagonal lattice) are particularly sensitive to the timescale on which the dithers are implemented, while others like random dithers perform well on all timescales. We then model the propagation of depth variations to artificial fluctuations in galaxy counts, which are a systematic for large-scale structure studies. We calculate the bias in galaxy counts caused by the observing strategy, accounting for photometric calibration uncertainties, dust extinction, and magnitude cuts; uncertainties in this bias limit our ability to account for structure induced by the observing strategy. We find that after 10 years of the LSST survey, the best dither strategies lead to uncertainties in this bias smaller than the minimum statistical floor for a galaxy catalog as deep as $r$$<$27.5. A few of these strategies bring the uncertainties close to the statistical floor for $r$$<$25.7 after only one year of survey.

We present the detection of Cepheids in the barred spiral galaxy NGC 1313, using the Wide Field and Planetary Camera 2 on the Hubble Space Telescpoe. Twenty B (F450W) and V (F555W) epochs of observations spanning over three weeks were obtained, on which the profile-fitting photometry of all stars in the monitored field was performed using the package HSTphot. A sample of 26 variable stars have been identified to be Cepheids, with periods between 3 and 14 days. Based on the derived period-luminosity relations in B and V bands, we obtain an extinction-corrected distance modulus of mu = 28.32 +- 0.08 (random) +- 0.06 (systematic), employing the Large Magellanic Cloud as the distance zero point calibrator. The above moduli correspond to a distance of 4.61 +- 0.17 (random) +- 0.13 (systematic) Mpc, consistent with previous measurements reported in the literature with uncertainties. In addition, the reddening to NGC 1313 is found to be small.

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.

The formation mechanisms of thick discs are under discussion. Thick discs might have formed either at high redshift on a short time-scale or might have been built slowly over time. They may have an internal or an external origin. Here we study in detail the kinematics and the stellar populations of the thick disc of ESO533-4. ESO533-4 is a nearby bulgeless galaxy. We present the first ever IFU study of an edge-on galaxy with enough depth to study the thick disc. We exposed ESO533-4 with VIMOS@VLT for 6.5hours. The FOV covered an axial extent 0.1-0.7r_25 (1-7kpc). We used pPXF and the MILES library to obtain velocity and stellar population maps. We compared our kinematic data with simple GADGET-2 models. The apparent rotational lag of the thick disc of ESO533-4 is compatible with that expected from the combinations of two effects: differential asymmetric drift and the projection effects arising from studying a disc a few degrees (2-3) away from edge-on. Thus, ESO533-4 contains little or no retrograde material. This is compatible with three formation scenarii: the secular heating of an initially thin disc, the formation of the thick disc at high redshift in a turbulent disc phase, and its creation in a major merger event. If happening in all galaxies, this last mechanism would cause retrograde thick discs in half of them. Retrograde discs have not been observed in the five massive disc galaxies (v_c>120km s^-1) for which thick disc kinematics are known. The populations of the thin and the thick discs are separated in the Age-log(Z/Z_Sun) plane. Thus, the thin and thick discs are made of two distinct stellar populations. Although the stellar population results are not conclusive due to the high dust extinction in ESO533-4, they do not favour a secular evolution origin for the thick disc. Hence, we suggest that the thick disc of ESO533-4 formed in a relatively short event (Abridged).

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.

We conduct a survey of satellite galaxies around the nearby spiral NGC 4258 by combining spectroscopic observations from the Apache Point Observatory 3.5-meter telescope with SDSS spectra. New spectroscopy is obtained for 15 galaxies. Of the 47 observed objects, we categorize 8 of them as probable satellites, 8 as possible satellites, and 17 as unlikely to be satellites. We do not speculate on the membership of the remaining 14 galaxies due to a lack of velocity and distance information. Radially integrating our best fit NFW profile for NGC 4258 yields a total mass of 1.8e12 Msun within 200 kpc. We find that the angular distribution of the satellites appears to be random, and not preferentially aligned with the disk of NGC 4258. In addition, many of the probable satellite galaxies have blue u-r colors and appear to be star-forming irregulars in SDSS images; this stands in contrast to the low number of blue satellites in the Milky Way and M31 systems at comparable distances.

Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data-cubes have 5.6'' FWHM spatial resolution, ~5A FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling galaxies out to ~0.7 R_25. These data-cubes can be used to produce 2D maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data-cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

We use the Mitchell Spectrograph (formerly VIRUS-P) to observe 12 nearby disk galaxies. We successfully measure ages in the outer disk in six systems. In three cases (NGC 2684, NGC 6155, and NGC 7437), we find that a downward break in the disk surface brightness profile corresponds with a change in the dominant stellar population with the interior being dominated by active star formation and the exterior having older stellar populations that are best-fit with star formation histories that decline with time. The observed increase in average stellar ages beyond a profile break is similar to theoretical models that predict surface brightness breaks are caused by stellar migration, with the outer disk being populated from scattered old interior stars. In three more cases (IC 1132, NGC 4904, and NGC 6691), we find no significant change in the stellar population as one crosses the break radius. In these galaxies, both the inner and outer disks are dominated by active star formation and younger stellar populations. While radial migration can contribute to the stellar populations beyond the break, it appears more than one mechanism is required to explain all of our observed stellar profile breaks.

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format.

Like the majority of spiral galaxies, NGC 6155 exhibits an exponential surface brightness profile that steepens significantly at large radii. Using the VIRUS-P IFU spectrograph, we have gathered spatially resolved spectra of the system. Modifying the GANDALF spectral fitting routine for use on the complex stellar populations found in spirals, we find that the average stellar ages increase significantly beyond the profile break radius. This result is in good agreement with recent simulations that predict the outskirts of disk galaxies are populated through stellar migration. With the ability to bin multiple fibers, we are able to measure stellar population ages down to mu_V~24 mag/sq arcsec.

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.

The ACS Nearby Galaxy Survey Treasury (ANGST) has acquired deep ACS imaging of a field in the outer disk of the large spiral galaxy M81. These data were obtained over a total of 20 HST orbits, providing a baseline long enough to reliably identify Cepheid variable stars in the field. Fundamental mode and first overtone types have been distinguished through comparative fits with corresponding Cepheid light curve templates derived from principal component analysis of confirmed Cepheids in the LMC, SMC, and Milky Way. A distance modulus of 27.78 pm 0.05_random pm 0.14_systematic with a corresponding distance of 3.60 pm 0.23 Mpc has been calculated from a sample of 11 fundamental mode and 2 first overtone Cepheids (assuming an LMC distance modulus of mu_LMC=18.41 pm 0.10_r pm 0.13_s).

We have generated accurate V and I template light curves using a combination of Fourier decomposition and principal component analysis for a large sample of Cepheid light curves. Unlike previous studies, we include short period Cepheids and stars pulsating in the first overtone mode in our analysis. Extensive Monte Carlo simulations show that our templates can be used to precisely measure Cepheid magnitudes and periods, even in cases where there are few observational epochs. These templates are ideal for characterizing serendipitously discovered Cepheids and can be used in conjunction with surveys such as Pan-Starrs and LSST where the observational sampling may not be optimized for Cepheids.

We have measured Lick index equivalent widths to derive luminosity weighted stellar ages and metallicities for thin and thick disk dominated regions of 9 edge-on disk galaxies with the ARC 3.5 meter telescope at Apache Point Observatory. In all cases, the thick disks are confirmed to be old stellar populations, with typical ages between 4 and 10 Gyr. The thin disks are uniformly younger than the thick disks, and show strong radial age gradients, with the outer regions of the disks being younger than 1 Gyr. We do not detect any significant metallicity differences or alpha-element enhancement in the thick disk stars compared to the thin disk, due to the insensitivity of the Lick indices to these differences at low metallicity. We compare these results to thick disks measured in other systems and to predictions from thick disk formation models.

We present kinematic measurements of thin and thick disk components in a sample of nine edge-on galaxies. We extract stellar and ionized gas rotation curves at and above the galaxies' midplanes using the Ca II triplet absorption features and H-alpha emission lines measured with the GMOS spectrographs on Gemini North and South. For the higher mass galaxies in the sample, we fail to detect differences between the thin and thick disk kinematics. In the lower mass galaxies, there is a wide range of thick disk behavior including thick disks with substantial lag and one counter-rotating thick disk. We compare our rotation curves with expectations from thick disk formation models and conclude that the wide variety of thick disk kinematics favors a formation scenario where thick disk stars are accreted or formed during merger events as opposed to models that form thick disks through gradual thin disk heating.

I briefly review the growing evidence that thick stellar disks surround most edge-on disk galaxies. Recent studies show that these extragalactic thick disks have old ages, low metallicities, long scale lengths, and moderately flattened axial ratios, much like the thick disk of the Milky Way. However, the properties of thick disks change systematically with the mass of the galaxy. The thick disks of low mass galaxies are more prominent and somewhat more metal-poor than those surrounding massive disk galaxies. Given the strong evidence that thick disks are fossils from an early epoch of merging, these trends place tight constraints on the early assembly of disk galaxies. The data suggest that the vast majority of the accretion onto the early disk was gaseous, but that the majority of thick disk stars were directly accreted.

We analyze the global structure of 34 late-type, edge-on, undisturbed, disk galaxies spanning a wide range of mass. We measure structural parameters for the galaxies using two-dimensional least-squares fitting to our $R$-band photometry. The fits require both a thick and a thin disk to adequately fit the data. The thick disks have larger scale heights and longer scale lengths than the embedded thin disks, by factors of ~2 and ~1.25, respectively. The observed structural parameters agree well with the properties of thick and thin disks derived from star counts in the Milky Way and from resolved stellar populations in nearby galaxies. We find that massive galaxies' luminosities are dominated by the thin disk. However, in low mass galaxies (Vc < 120 km/s), thick disk stars contribute nearly half of the luminosity and dominate the stellar mass. Thus, although low mass dwarf galaxies appear blue, the majority of their stars are probably quite old. Our data are most easily explained by a formation scenario where the thick disk is assembled through direct accretion of stellar material from merging satellites while the thin disk is formed from accreted gas. The baryonic fraction in the thin disk therefore constrains the gas-richness of the merging pre-galactic fragments. If we include the mass in HI as part of the thin disk, the thick disk contains <10% of the baryons in high mass galaxies, and ~25-30% of the baryons in low-mass galaxies. We discuss how our trends can be explained by supernova-driven outflow at early times as well as the possibilities for predicting abundance trends in thick disks, and for removing discrepancies between semi-analytic galaxy formation models and the observed colors of low mass galaxies. (abstract abridged)

We present kinematic measurements of the thick and thin disks in two edge-on galaxies. We have derived stellar rotation curves at and above the galaxies' midplanes using Ca II triplet features measured with the GMOS spectrograph on Gemini North. In one galaxy, FGC 1415, the kinematics above the plane show clear rotation that lags that of the midplane by ~20-50%, similar to the behavior seen in the Milky Way. However, the kinematics of the second galaxy, FGC 227, are quite different. The rotation above the plane is extremely slow, showing <25% of the rotation speed of the stars at the midplane. We decompose the observed rotation curves into a superposition of thick and thin disk kinematics, using 2-dimensional fits to the galaxy images to determine the fraction of thick disk stars at each position. We find that the thick disk of FGC 1415 rotates at 30-40% of the rotation speed of the thin disk. In contrast, the thick disk of FGC 227 is very likely counter-rotating, if it is rotating at all. These observations are consistent with the velocity dispersion profiles we measure for each galaxy. The detection of counter-rotating thick disks conclusively rules out models where the thick disk forms either during monolithic collapse or from vertical heating of a previous thin disk. Instead, the data strongly support models where the thick disk forms from direct accretion of stars from infalling satellites.

We find that disk galaxies show a sharp, mass-dependent transition in the structure of their dusty ISM. Dust lanes are a generic feature of massive disks with V_rot>120km/s, but are completely absent in galaxies with V_rot<120km/s. The transition reflects an increase in the scale height of the cold ISM in low mass galaxies, driven by larger turbulent velocities supporting the gas layer, rather than sharp drops in the gas surface density. We identify the V_rot=120km/s transition with the onset of gravitational instabilities in high mass galaxies. The instabilities lead to fragmentation and gravitational collapse along spiral arms, smaller gas scale heights, lower turbulent velocities, and thus to narrow dust lanes. The drop in velocity dispersion may be due either to a switch in the driving mechanism for turbulence or to a change in the response of the ISM to supernovae after the ISM has collapsed to a dense layer. The resulting smaller gas scale height can lead to significant increases in the star formation rate when disk instabilities are present, and may explain the Kennicutt surface density threshold for star formation. Our data suggest that star formation will be systematically less efficient in low mass disks with V_c<120km/s, leading to star formation timescales longer than the gas accretion timescale. This effect can suppress the metallicity and nucleosynthetic yields of low mass disks, and thus explain the disk mass-metallicity relationship without invoking galactic SN-driven outflows. The transitions in disk stability, dust structure, and/or star formation efficiency may also be responsible for observed changes in the slope of the Tully-Fisher relation, in the sharp increase in the thickness of dwarf galaxy disks, and in the onset of bulges in galaxies with V_rot>120km/s. (Abridged)