Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

PG 1553+113 is the first blazar showing an approximately two-year quasi-periodic pattern in its gamma-ray light curve. Such quasi-periodicity might have a geometrical origin, possibly related to the precessing nature of the jet, or could be intrinsic to the source and related to pulsational accretion flow instabilities. By means of a ~2yr very long baseline array (VLBA) monitoring at 15, 24, and 43 GHz we investigate the source pc-scale properties during an entire cycle of gamma-ray activity in the period 2015-2017. In contrast to the well-defined periodicity in the gamma-ray emission, at radio frequencies no clear periodic pattern can be recognized. The jet position angle, constrained by means of the total intensity ridge line, varies across the different observing epochs in the range 40-60 deg. We also investigate the time evolution of the source polarization properties, including the rotation measure. The brightness temperature is found to decrease as the frequency increases with an intrinsic value of ~1.5 x 10^10 K and the estimated Doppler factor is ~1.4.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

Investigating the magnetic field structure in the innermost regions of relativistic jets is fundamental to shed light on the crucial physical processes giving rise to the jet formation, as well as to its extraordinary radiation output up to gamma-ray energies. We study the magnetic field structure of the quasar CTA 102 with 3 and 7 mm-VLBI polarimetric observations, reaching an unprecedented resolution (~50 microarcsec). We also investigate the variability and physical process occurring in the source during the observing period which coincides with a very active state of the source till high-energies. The Faraday rotation analysis between 3 and 7mm shows a gradient in rotation measure with a maximum value of ~6X10^4 rad/m^2 and intrinsic electric vector position angles (EVPAs) oriented around the centroid of the core, suggesting the presence of large-scale helical magnetic fields. Such a magnetic field structure is also visible in 7 mm images when a new superluminal component is crossing the core region. The 7mm EVPAs orientation is different when the component is exiting the core or crossing a stationary feature at ~0.1 mas. The interaction between the superluminal component and a recollimation shock at ~0.1 mas could have triggered the multi-wavelengths flares. The variability Doppler factor associated with such interaction is large enough to explain the high energy emission, as we infer from the analysis of gamma-ray and X-ray data, and it is in agreement with the Doppler factor obtained to explain the extraordinary optical flare by Raiteri et al.(2017).

We describe the design and performance of the near-infrared (1.51--1.70 micron), fiber-fed, multi-object (300 fibers), high resolution (R = lambda/delta lambda ~ 22,500) spectrograph built for the Apache Point Observatory Galactic Evolution Experiment (APOGEE). APOGEE is a survey of ~ 10^5 red giant stars that systematically sampled all Milky Way populations (bulge, disk, and halo) to study the Galaxy's chemical and kinematical history. It was part of the Sloan Digital Sky Survey III (SDSS-III) from 2011 -- 2014 using the 2.5 m Sloan Foundation Telescope at Apache Point Observatory, New Mexico. The APOGEE-2 survey is now using the spectrograph as part of SDSS-IV, as well as a second spectrograph, a close copy of the first, operating at the 2.5 m du Pont Telescope at Las Campanas Observatory in Chile. Although several fiber-fed, multi-object, high resolution spectrographs have been built for visual wavelength spectroscopy, the APOGEE spectrograph is one of the first such instruments built for observations in the near-infrared. The instrument's successful development was enabled by several key innovations, including a "gang connector" to allow simultaneous connections of 300 fibers; hermetically sealed feedthroughs to allow fibers to pass through the cryostat wall continuously; the first cryogenically deployed mosaic volume phase holographic grating; and a large refractive camera that includes mono-crystalline silicon and fused silica elements with diameters as large as ~ 400 mm. This paper contains a comprehensive description of all aspects of the instrument including the fiber system, optics and opto-mechanics, detector arrays, mechanics and cryogenics, instrument control, calibration system, optical performance and stability, lessons learned, and design changes for the second instrument.

We report results from a deep polarization imaging of the nearby radio galaxy 3C$\,$84 (NGC$\,$1275). The source was observed with the Global Millimeter VLBI Array (GMVA) at 86$\,$GHz at an ultra-high angular resolution of $50\mu$as (corresponding to 250$R_{s}$). We also add complementary multi-wavelength data from the Very Long Baseline Array (VLBA; 15 & 43$\,$GHz) and from the Atacama Large Millimeter/submillimeter Array (ALMA; 97.5, 233.0, and 343.5$\,$GHz). At 86$\,$GHz, we measure a fractional linear polarization of $\sim2$% in the VLBI core region. The polarization morphology suggests that the emission is associated with an underlying limb-brightened jet. The fractional linear polarization is lower at 43 and 15$\,$GHz ($\sim0.3-0.7$% and $<0.1$%, respectively). This suggests an increasing linear polarization degree towards shorter wavelengths on VLBI scales. We also obtain a large rotation measure (RM) of $\sim10^{5-6}~{\rm rad/m^{2}}$ in the core at $\gtrsim$43$\,$GHz. Moreover, the VLBA 43$\,$GHz observations show a variable RM in the VLBI core region during a small flare in 2015. Faraday depolarization and Faraday conversion in an inhomogeneous and mildly relativistic plasma could explain the observed linear polarization characteristics and the previously measured frequency dependence of the circular polarization. Our Faraday depolarization modeling suggests that the RM most likely originates from an external screen with a highly uniform RM distribution. To explain the large RM value, the uniform RM distribution, and the RM variability, we suggest that the Faraday rotation is caused by a boundary layer in a transversely stratified jet. Based on the RM and the synchrotron spectrum of the core, we provide an estimate for the magnetic field strength and the electron density of the jet plasma.

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ~287 as revealed by multi-epoch Very Long Baseline Array (VLBA) images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle (EVPA) of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave "core" or, at other epochs, coincident with a bright, quasi-stationary emission feature $\sim0.2$~milliarcsec ($\sim$0.9~pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and $\gamma$-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies $\geq10$~pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies.

Low ($\lesssim 1\%$) levels of circular polarization (CP) detected at radio frequencies in the relativistic jets of some blazars can provide insight into the underlying nature of the jet plasma. CP can be produced through linear birefringence, in which initially linearly polarized emission produced in one region of the jet is altered by Faraday rotation as it propagates through other regions of the jet with varying magnetic field orientation. Marscher has begun a study of jets with such magnetic geometries using the Turbulent Extreme Multi-Zone (TEMZ) model, in which turbulent plasma crossing a standing shock in the jet is represented by a collection of thousands of individual plasma cells, each with distinct magnetic field orientations. Here we develop a radiative transfer scheme that allows the numerical TEMZ code to produce simulated images of the time-dependent linearly and circularly polarized intensity at different radio frequencies. In this initial study, we produce synthetic polarized emission maps that highlight the linear and circular polarization expected within the model.

Blazars exhibit flares across the entire electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at longer wavelengths; however, a small subset appears to occur in isolation, with little or no correlated variability at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. Macdonald et al. have developed the Ring of Fire model to explain the origin of orphan $\gamma$-ray flares from within blazar jets. In this model, electrons contained within a blob of plasma moving relativistically along the spine of the jet inverse-Compton scatter synchrotron photons emanating off of a ring of shocked sheath plasma that enshrouds the jet spine. As the blob propagates through the ring, the scattering of the ring photons by the blob electrons creates an orphan $\gamma$-ray flare. This model was successfully applied to modeling a prominent orphan $\gamma$-ray flare observed in the blazar PKS 1510$-$089. To further support the plausibility of this model, Macdonald et al. presented a stacked radio map of PKS 1510$-$089 containing the polarimetric signature of a sheath of plasma surrounding the spine of the jet. In this paper, we extend our modeling and stacking techniques to a larger sample of blazars: 3C 273, 4C 71$.$01, 3C 279, 1055$+$018, CTA 102, and 3C 345, the majority of which have exhibited orphan $\gamma$-ray flares. We find that the model can successfully reproduce these flares, while our stacked maps reveal the existence of jet sheaths within these blazars.

We analyze the parsec-scale jet kinematics from 2007 June to 2013 January of a sample of $\gamma$-ray bright blazars monitored roughly monthly with the Very Long Baseline Array at 43~GHz. In a total of 1929 images, we measure apparent speeds of 252 emission knots in 21 quasars, 12 BL~Lacertae objects (BLLacs), and 3 radio galaxies, ranging from 0.02$c$ to 78$c$; 21\% of the knots are quasi-stationary. Approximately 1/3 of the moving knots execute non-ballistic motions, with the quasars exhibiting acceleration along the jet within 5~pc (projected) of the core, and knots in the BLLacs tending to decelerate near the core. Using apparent speeds of components and timescales of variability from their light curves, we derive physical parameters of 120 superluminal knots, including variability Doppler factors, Lorentz factors, and viewing angles. We estimate the half-opening angle of each jet based on the projected opening angle and scatter of intrinsic viewing angles of knots. We determine characteristic values of physical parameters for each jet and AGN class based on the range of values obtained for individual features. We calculate intrinsic brightness temperatures of the cores, $T_{\rm b,int}^{\rm core}$, at all epochs, finding that the radio galaxies usually maintain equipartition conditions in the cores, while $\sim$30\% of $T_{\rm b,int}^{\rm core}$ measurements in the quasars and BLLacs deviate from equipartition values by a factor $>$10. This probably occurs during transient events connected with active states. In the Appendix we briefly describe the behavior of each blazar during the period analyzed.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We present a sample of local red giant stars observed using the New Mexico State University 1 m telescope with the APOGEE spectrograph, for which we estimate stellar ages and the age distribution from the high-resolution spectroscopic stellar parameters and accurate distance measurements from Hipparcos. The high-resolution (R ~ 23,000), near infrared (H-band, 1.5-1.7 micron) APOGEE spectra provide measurements of the stellar atmospheric parameters (temperature, surface gravity, [M/H], and [alpha/M]). Due to the smaller uncertainties in surface gravity possible with high-resolution spectra and accurate Hipparcos distance measurements, we are able to calculate the stellar masses to within 40%. For red giants, the relatively rapid evolution of stars up the red giant branch allows the age to be constrained based on the mass. We examine methods of estimating age using both the mass-age relation directly and a Bayesian isochrone matching of measured parameters, assuming a constant star formation history (SFH). To improve the prior on the SFH, we use a hierarchical modeling approach to constrain the parameters of a model SFH from the age probability distribution functions of the data. The results of an alpha dependent Gaussian SFH model shows a clear relation between age and [alpha/M] at all ages. Using this SFH model as the prior for an empirical Bayesian analysis, we construct a full age probability distribution function and determine ages for individual stars. The age-metallicity relation is flat, with a slight decrease in [M/H] at the oldest ages and a ~ 0.5 dex spread in metallicity. For stars with ages < 1 Gyr we find a smaller spread, consistent with radial migration having a smaller effect on these young stars than on the older stars.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

We present the design and test results of a compact optical fiber double-scrambler for high-resolution Doppler radial velocity instruments. This device consists of a single optic: a high-index $n$$\sim$2 ball lens that exchanges the near and far fields between two fibers. When used in conjunction with octagonal fibers, this device yields very high scrambling gains and greatly desensitizes the fiber output from any input illumination variations, thereby stabilizing the instrument profile of the spectrograph and improving the Doppler measurement precision. The system is also highly insensitive to input pupil variations, isolating the spectrograph from telescope illumination variations and seeing changes. By selecting the appropriate glass and lens diameter the highest efficiency is achieved when the fibers are practically in contact with the lens surface, greatly simplifying the alignment process when compared to classical double-scrambler systems. This prototype double-scrambler has demonstrated significant performance gains over previous systems, achieving scrambling gains in excess of 10,000 with a throughput of $\sim$87% using uncoated Polymicro octagonal fibers. Adding a circular fiber to the fiber train further increases the scrambling gain to $>$20,000, limited by laboratory measurement error. While this fiber system is designed for the Habitable-zone Planet Finder spectrograph, it is more generally applicable to other instruments in the visible and near-infrared. Given the simplicity and low cost, this fiber scrambler could also easily be multiplexed for large multi-object instruments.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is an integral-field spectroscopic survey of 10,000 nearby galaxies that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). MaNGA's 17 pluggable optical fiber-bundle integral field units (IFUs) are deployed across a 3 deg field, they yield spectral coverage 3600-10,300 Ang at a typical resolution R ~ 2000, and sample the sky with 2" diameter fiber apertures with a total bundle fill factor of 56%. Observing over such a large field and range of wavelengths is particularly challenging for obtaining uniform and integral spatial coverage and resolution at all wavelengths and across each entire fiber array. Data quality is affected by the IFU construction technique, chromatic and field differential refraction, the adopted dithering strategy, and many other effects. We use numerical simulations to constrain the hardware design and observing strategy for the survey with the aim of ensuring consistent data quality that meets the survey science requirements while permitting maximum observational flexibility. We find that MaNGA science goals are best achieved with IFUs composed of a regular hexagonal grid of optical fibers with rms displacement of 5 microns or less from their nominal packing position, this goal is met by the MaNGA hardware, which achieves 3 microns rms fiber placement. We further show that MaNGA observations are best obtained in sets of three 15-minute exposures dithered along the vertices of a 1.44 arcsec equilateral triangle, these sets form the minimum observational unit, and are repeated as needed to achieve a combined signal-to-noise ratio of 5 per Angstrom per fiber in the r-band continuum at a surface brightness of 23 AB/arcsec^2. (abbrev.)

Blazars exhibit flares across the electromagnetic spectrum. Many $\gamma$-ray flares are highly correlated with flares detected at optical wavelengths; however, a small subset appears to occur in isolation, with little or no variability detected at longer wavelengths. These "orphan" $\gamma$-ray flares challenge current models of blazar variability, most of which are unable to reproduce this type of behavior. We present numerical calculations of the time-variable emission of a blazar based on a proposal by Marscher et al. (2010) to explain such events. In this model, a plasmoid ("blob") propagates relativistically along the spine of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a spatially varying magnetic field, and acceleration of the blob's bulk velocity. Synthetic light curves produced by our model are compared to the observed light curves from an orphan flare that was coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510$-$089. In addition, we present Very Long Baseline Array polarimetric observations that point to the existence of a jet sheath in PKS 1510$-$089, thus providing further observational support for the plausibility of our model. An estimate of the bolometric luminosity of the sheath within PKS 1510$-$089 is made, yielding $L_{\rm sh} \sim 3 \times 10^{45} ~ \rm erg ~ \rm s^{-1}$. This indicates that the sheath within PKS 1510$-$089 is potentially a very important source of seed photons.

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 describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet.

We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent \gamma-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each \gamma-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the \gamma-ray outbursts, although the \gamma-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the \gamma-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the \gamma-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that \gamma-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells.

We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results.

We present multiwavelength data of the blazar 3C 454.3 obtained during an extremely bright outburst from November 2010 through January 2011. These include flux density measurements with the Herschel Space Observatory at five submillimeter-wave and far-infrared bands, the Fermi Large Area Telescope at gamma-ray energies, Swift at X-ray, ultraviolet (UV), and optical frequencies, and the Submillimeter Array at 1.3 mm. From this dataset, we form a series of 52 spectral energy distributions (SEDs) spanning nearly two months that are unprecedented in time coverage and breadth of frequency. Discrete correlation anlaysis of the millimeter, far-infrared, and gamma-ray light curves show that the variations were essentially simultaneous, indicative of co-spatiality of the emission, at these wavebands. In contrast, differences in short-term fluctuations at various wavelengths imply the presence of inhomegeneities in physical conditions across the source. We locate the site of the outburst in the parsec-scale core, whose flux density as measured on 7 mm Very Long Baseline Array images increased by 70 percent during the first five weeks of the outburst. Based on these considerations and guided by the SEDs, we propose a model in which turbulent plasma crosses a conical standing shock in the parsec-scale region of the jet. Here, the high-energy emission in the model is produced by inverse Compton scattering of seed photons supplied by either nonthermal radiation from a Mach disk, thermal emission from hot dust, or (for X-rays) synchrotron radiation from plasma that crosses the standing shock. For the two dates on which we fitted the model SED to the data, the model corresponds very well to the observations at all bands except at X-ray energies, where the spectrum is flatter than observed.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We analyze total and polarized intensity images of the quasar 3C 454.3 obtained monthly with the VLBA at 43 GHz within the ongoing Boston U. monitoring program of gamma-ray blazars started in June 2007. The data are supplemented by VLBA observations performed during intense campaigns of 2 week duration when the quasar was observed 3 times per campaign. We find a strong increase of activity in the parsec-scale jet of the quasar during high gamma-ray states in December 2009, April 2010, and November 2010. We detect new superluminal knots, K09 and K10, associated with the autumn 2009 and 2010 outbursts, respectively, and compare their kinematic parameters. We analyze optical polarimetric behavior along with polarization parameters of the parsec-scale jet and outline similarities and differences in polarization properties across wavelengths. The results of the analysis support the conclusions that the optical polarized emission is produced in a region located in the vicinity of the mm-wave core of the jet of the quasar, and that the gamma-ray outbursts occur when a superluminal disturbance passes through the core.

Analysis of comprehensive monitoring of 34 gamma-ray bright quasars, BL Lac objects, and radio galaxies reveals a close connection between events in the millimeter-wave emission imaged with the VLBA at 43 GHz and flares at gamma-ray and lower frequencies. Roughly 2/3 of the flares are coincident with the appearance of a new superluminal knot and/or a flare in the millimeter-wave "core'" located parsecs from the central engine. This presents a theoretical challenge to explain how the gamma-ray flux can often be variable on intra-day time-scales. Possible answers to this include very narrow opening angles of the jet, small volume filling factors of the highest energy electrons, chaotic magnetic fields, and turbulent velocity fields relative to the mean jet flow.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.