Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars, with $L_X > 10^{44}\,$erg/s and $M_B < -23$. We define a new population of quasars which are unobscured, have X-ray luminosities in excess of $10^{44}\,$erg/s, have broad emission lines, and yet are also bright in the far-infrared, with a 250$\mu$m flux density of $S_{\rm 250}>30$mJy. We refer to these Herschel-detected, unobscured quasars as "Cold Quasars". A mere 4% (21) of the X-ray- and optically-selected unobscured quasars in Stripe 82X are detected at 250$\mu$m. These Cold Quasars lie at $z\sim1-3$, have $L_{\rm IR}>10^{12}\,L_\odot$, and have star formation rates of $\sim200-1400\,M_\odot$/yr. Cold Quasars are bluer in the mid-IR than the full quasar population, and 72% of our Cold Quasars have WISE W3 $<$ 11.5 [Vega], while only 19% of the full quasar sample meets this criteria. Crucially, Cold Quasars have on average $\sim9\times$ as much star formation as the main sequence of star forming galaxies at similar redshifts. Although dust-rich, unobscured quasars have occasionally been noted in the literature before, we argue that they should be considered as a separate class of quasars due to their high star formation rates. This phase is likely short-lived, as the central engine and immense star formation consume the gas reservoir. Cold Quasars are type-1 blue quasars that reside in starburst galaxies.

There is no unique way to describe the dark energy-dark matter interaction, as we have little information about the nature and dynamics of the dark sector. Hence, in many of the phenomenological dark matter fluid interaction models in the literature, the interaction strength $Q_{\nu}$ in the dark sector is introduced by hand. Demanding that the interaction strength $Q_{\nu}$ in the dark sector must have a field theory description, we obtain a unique form of interaction strength. We show the equivalence between the fields and fluids for the $f(R,\chi)$ model where $f$ is an arbitrary, smooth function of $R$ and the scalar field $\chi$, which represents dark matter. Up to first order in perturbations, we show that the one-to-one mapping between the field theory description and the phenomenological fluid description of interacting dark energy and dark matter exists \emph{only} for this unique form of interaction. We then classify the interacting dark energy models considered in the literature into two categories based on the field-theoretic description.

Gamma-ray burst (GRB) data suggest that the jets from GRBs in the high redshift universe are more narrowly collimated than those at lower redshifts. This implies that we detect relatively fewer long GRB progenitor systems (i.e. massive stars) at high redshifts, because a greater fraction of GRBs have their jets pointed away from us. As a result, estimates of the star formation rate (from the GRB rate) at high redshifts may be diminished if this effect is not taken into account. In this paper, we estimate the star formation rate (SFR) using the observed GRB rate, accounting for an evolving jet opening angle. We find that the SFR in the early universe (z > 3) can be up to an order of magnitude higher than the canonical estimates, depending on the severity of beaming angle evolution and the fraction of stars that make long gamma-ray bursts. Additionally, we find an excess in the SFR at low redshifts, although this lessens when accounting for evolution of the beaming angle. Finally, under the assumption that GRBs do in fact trace canonical forms of the cosmic SFR, we constrain the resulting fraction of stars that must produce GRBs, again accounting for jet beaming-angle evolution. We find this assumption suggests a high fraction of stars in the early universe producing GRBs - a result that may, in fact, support our initial assertion that GRBs do not trace canonical estimates of the SFR.

van den Heuvel & Tauris argue that if the red giant star in the system 2MASS J05215658+4359220 has a mass of 1 solar mass (M$_\odot$), then its unseen companion could be a binary composed of two 0.9 M$_\odot$ stars, making a triple system. We contend that the existing data are most consistent with a giant of mass $3.2^{+1.0}_{-1.0}$ M$_\odot$, implying a black hole companion of $3.3^{+2.8}_{-0.7}$ M$_\odot$.

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We investigate the Fe, C, N, O, Mg, Al, Si, K, Ca, Ce and Nd abundances of 2283 red giant stars in 31 globular clusters from high-resolution spectra observed in both the northern and southern hemisphere by the SDSS-IV APOGEE-2 survey. This unprecedented homogeneous dataset, largest to date, allows us to discuss the intrinsic Fe spread, the shape and statistics of Al-Mg and N-C anticorrelations as a function of cluster mass, luminosity, age and metallicity for all 31 clusters. We find that the Fe spread does not depend on these parameters within our uncertainties including cluster metallicity, contradicting earlier observations. We do not confirm the metallicity variations previously observed in M22 and NGC 1851. Some clusters show a bimodal Al distribution, while others exhibit a continuous distribution as has been previously reported in the literature. We confirm more than 2 populations in $\omega$ Cen and NGC 6752, and find new ones in M79. We discuss the scatter of Al by implementing a correction to the standard chemical evolution of Al in the Milky Way. After correction, its dependence on cluster mass is increased suggesting that the extent of Al enrichment as a function of mass was suppressed before the correction. We observe a turnover in the Mg-Al anticorrelation at very low Mg in $\omega$ Cen, similar to the pattern previously reported in M15 and M92. $\omega$ Cen may also have a weak K-Mg anticorrelation, and if confirmed, it would be only the third cluster known to show such a pattern.

Jet opening angles of long gamma-ray bursts (lGRBs) appear to evolve in cosmic time, with lGRBs at higher redshifts being on average more narrowly beamed than those at lower redshifts. We examine the nature of this anti-correlation in the context of collimation by the progenitor stellar envelope. First, we show that the data indicate a strong correlation between gamma-ray luminosity and jet opening angle, and suggest this is a natural selection effect - only the most luminous GRBs are able to successfully launch jets with large opening angles. Then, by considering progenitor properties expected to evolve through cosmic time, we show that denser stars lead to more collimated jets; we argue that the apparent anti-correlation between opening angle and redshift can be accounted for if lGRB massive star progenitors at high redshifts have higher average density compared to those at lower redshifts. This may be viable for an evolving IMF - under the assumption that average density scales directly with mass, this relationship is consistent with the form of the IMF mass evolution suggested in the literature. The jet angle-redshift anti-correlation may also be explained if the lGRB progenitor population is dominated by massive stars at high redshift, while lower redshift lGRBs allow for a greater diversity of progenitor systems (that may fail to collimate the jet as acutely). Overall, however, we find both the jet angle-redshift anti-correlation and jet angle-luminosity correlation are consistent with the conditions of jet launch through, and collimation by, the envelope of a massive star progenitor.

The clustering of active galactic nuclei (AGN) sheds light on their typical large (Mpc-scale) environments, which can constrain the growth and evolution of supermassive black holes. Here we measure the clustering of luminous X-ray-selected AGN in the Stripe 82X and XMM-XXL-North surveys around the peak epoch of black hole growth, in order to investigate the dependence of luminosity on large-scale AGN environment. We compute the auto-correlation function of AGN in two luminosity bins, $10^{43}\leq L_X<10^{44.5}$ erg s$^{-1}$ at $z\sim 0.8$ and $L_X\geq 10^{44.5}$ erg s$^{-1}$ at $z\sim 1.8$, and calculate the AGN bias taking into account the redshift distribution of the sources using three different methods. Our results show that while the less luminous sample has an inferred typical halo mass that is smaller than for the more luminous AGN, the host halo mass may be less dependent on luminosity than suggested in previous work. Focusing on the luminous sample, we calculate a typical host halo mass of $\sim 10^{13}$ M$_{\odot}~h^{-1}$, which is similar to previous measurements of moderate-luminosity X-ray AGN and significantly larger than the values found for optical quasars of similar luminosities and redshifts. We suggest that the clustering differences between different AGN selection techniques are dominated by selection biases, and not due to a dependence on AGN luminosity. We discuss the limitations of inferring AGN triggering mechanisms from halo masses derived by large-scale bias.

Lithium abundances are presented for 91 dwarf and subgiant stars in the Galactic bulge. The analysis is based on line synthesis of the 7Li line at 6707 {\AA} in high-resolution spectra obtained during gravitational microlensing events, when the brightnesses of the targets were highly magnified. Our main finding is that the bulge stars at sub-solar metallicities, and that are older than about eight billion years, does not show any sign of Li production, that is, the Li trend with metallicity is flat (or even slightly declining). This indicates that no lithium was produced during the first few billion years in the history of the bulge. This finding is essentially identical to what is seen for the (old) thick disk stars in the Solar neighbourhood, and adds another piece of evidence for a tight connection between the metal-poor bulge and the Galactic thick disk. For the bulge stars younger than about eight billion years, the sample contains a group of stars at very high metallicities at [Fe/H]~+0.4 that have lithium abundances in the range A(Li)=2.6-2.8. In the Solar neighbourhood the lithium abundances have been found to peak at a A(Li)~3.3 at [Fe/H]~ +0.1 and then decrease by 0.4-0.5 dex when reaching [Fe/H]~+0.4. The few bulge stars that we have at these metallicities, seem to support this declining A(Li) trend. This could indeed support the recent claim that the low A(Li) abundances at the highest metallicities seen in the Solar neighbourhood could be due to stars from the inner disk, or the bulge region, that have migrated to the Solar neighbourhood.

We present spectroscopic determinations of the effective temperatures, surface gravities and metallicities for 21 M-dwarfs observed at high-resolution (R $\sim$ 22,500) in the \textit{H}-band as part of the SDSS-IV APOGEE survey. The atmospheric parameters and metallicities are derived from spectral syntheses with 1-D LTE plane parallel MARCS models and the APOGEE atomic/molecular line list, together with up-to-date H$_{2}$O and FeH molecular line lists. Our sample range in $T_{\rm eff}$ from $\sim$ 3200 to 3800K, where eleven stars are in binary systems with a warmer (FGK) primary, while the other 10 M-dwarfs have interferometric radii in the literature. We define an $M_{K_{S}}$--Radius calibration based on our M-dwarf radii derived from the detailed analysis of APOGEE spectra and Gaia DR2 distances, as well as a mass-radius relation using the spectroscopically-derived surface gravities. A comparison of the derived radii with interferometric values from the literature finds that the spectroscopic radii are slightly offset towards smaller values, with $\Delta$ = -0.01 $\pm$ 0.02 $R{\star}$/$R_{\odot}$. In addition, the derived M-dwarf masses based upon the radii and surface gravities tend to be slightly smaller (by $\sim$5-10\%) than masses derived for M-dwarf members of eclipsing binary systems for a given stellar radius. The metallicities derived for the 11 M-dwarfs in binary systems, compared to metallicities obtained for their hotter FGK main-sequence primary stars from the literature, shows excellent agreement, with a mean difference of [Fe/H](M-dwarf - FGK primary) = +0.04 $\pm$ 0.18 dex, confirming the APOGEE metallicity scale derived here for M-dwarfs.

Our knowledge of the populations and occurrence rates of planets orbiting evolved intermediate-mass stars lags behind that for solar-type stars by at least a decade. Some radial velocity surveys have targeted these low-luminosity giant stars, providing some insights into the properties of their planetary systems. Here we present the final data release of the Pan-Pacific Planet Search, a 5-year radial velocity survey using the 3.9m Anglo-Australian Telescope. We present 1293 precise radial velocity measurements for 129 stars, and highlight six potential substellar-mass companions which require additional observations to confirm. Correcting for the substantial incompleteness in the sample, we estimate the occurrence rate of giant planets orbiting low-luminosity giant stars to be approximately 7.8$^{+9.1}_{-3.3}$\%. This result is consistent with the frequency of such planets found to orbit main-sequence A-type stars, from which the PPPS stars have evolved.

We present a computationally tractable implementation of spectro-perfectionism, a method which minimizes error imparted by spectral extraction. We develop our method in conjunction with a full raw reduction pipeline for the MINiature Exoplanet Radial Velocity Array (MINERVA), capable of performing both optimal extraction and spectro-perfectionism. Although spectro-perfectionism remains computationally expensive, our implementation can extract a MINERVA exposure in approximately $30\,\text{min}$. We describe our localized extraction procedure and our approach to point spread function fitting. We compare the performance of both extraction methods on a set of 119 exposures on HD122064, an RV standard star. Both the optimal extraction and spectro-perfectionism pipelines achieve nearly identical RV precision under a six-exposure chronological binning. We discuss the importance of reliable calibration data for point spread function fitting and the potential of spectro-perfectionism for future precise radial velocity exoplanet studies.

We investigate the impact of bursts in star formation on the predictions of one-zone chemical evolution models, adopting oxygen (O), iron (Fe), and strontium (Sr), as representative $\alpha$, iron-peak, and s-process elements, respectively. To this end, we develop the Versatile Integrator for Chemical Evolution (VICE), a python package. Starbursts driven by a temporary boost of gas accretion rate create loops in [O/Fe]-[Fe/H] evolutionary tracks and a peak in the stellar [O/Fe] distribution at intermediate values. Bursts driven by a temporary boost of star formation efficiency have a similar effect, and they also produce a population of $\alpha$-deficient stars during the depressed star formation phase that follows the burst. This $\alpha$-deficient population is more prominent if the outflow rate is tied to a time-averaged star formation rate (SFR) instead of the instantaneous SFR. Theoretical models of Sr production predict a strong metallicity dependence of supernova and asymptotic giant branch (AGB) star yields, though comparison to data suggests an additional source that is nearly metallicity-independent. Evolution of [Sr/Fe] and [Sr/O] during a starburst is complex because of the yield metallicity dependence and the multiple timescales in play. Moderate amplitude (10-20\%) sinusoidal oscillations in SFR produce loops in [O/Fe]-[Fe/H] tracks and multiple peaks in [O/Fe] distributions, which could be one source of intrinsic scatter in observed sequences. We investigate models that have a factor of ~2 enhancement of SFR at t = 12 Gyr, as suggest by some recent Milky Way observations. A late episode of enhanced star formation could help explain the existence of young stars with moderate $\alpha$-enhancements and the surprisingly young median age found for solar metallicity stars in the solar neighborhood, while also raising the possibility that this starburst has not fully decayed.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($\mu_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer.

The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at $R_o$. In this Letter, we show that a non-uniform density at $R_o$ causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient $|\partial{\ln{\rho}}/\partial{\theta}|$ above which this change of the solution occurs. For realistic $R_o$, this critical gradient is less than 0.01 and tends to 0 as $R_o \rightarrow \infty$. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogenous density boundary condition remains constant at nearly the Bondi rate $\dot{M}_B$, while the outflow rate can greatly exceed $\dot{M}_B$.

The direct formation of a massive black hole is a potential seeding mechanism of the earliest observed supermassive black holes. We investigate how the existence of a massive black hole seed impacts the ionization and thermal state of its pre-galactic host halo and subsequent star formation. We show that its X-ray radiation ionizes and heats the medium, enhancing $\rm{H}_2$ formation in shielded regions, within the nuclear region in the span of a million years. The enhanced molecular cooling triggers the formation of a $\sim 10^4~{\rm M}_\odot$ metal-free stellar cluster at a star formation efficiency of $\sim 0.1\%$ in a single event. Star formation occurs near the edges of the H II region that is partially ionized by X-rays, thus the initial size depends on the black hole properties and surrounding environment. The simulated metal-free galaxy has an initial half-light radius of $\sim 10$ pc but expands to $\sim 50$ pc after 10 million years because of the outward velocities of their birth clouds. Supernova feedback then quenches any further star formation for tens of millions of years, allowing the massive black hole to dominate the spectrum once the massive metal-free stars die.

Using a sample of 70 924 stars from the second data release of the GALAH optical spectroscopic survey, we construct median sequences of [X/Mg] vs. [Mg/H] for 21 elements, separating the high-$\alpha$/``low-Ia'' and low-$\alpha$/``high-Ia'' stellar populations through cuts in [Mg/Fe]. Previous work with the near-IR APOGEE survey has shown that such sequences are nearly independent of location in the Galactic disk, implying that they are determined by stellar nucleosynthesis yields with little sensitivity to other chemical evolution aspects. The separation between the two [X/Mg] sequences indicates the relative importance of prompt and delayed enrichment mechanisms, while the sequences' slopes indicate metallicity dependence of the yields. GALAH and APOGEE measurements agree for some of their common elements, but differ in sequence separation or metallicity trends for others. GALAH offers access to nine new elements. We infer that about $75\%$ of solar C comes from core collapse supernovae and $25\%$ from delayed mechanisms. We find core collapse fractions of $60-80\%$ for the Fe-peak elements Sc, Ti, Cu, and Zn, with strong metallicity dependence of the core collapse Cu yield. For the neutron capture elements Y, Ba, and La, we infer large delayed contributions with non-monotonic metallicity dependence. The separation of the [Eu/Mg] sequences implies that at least $\sim30\%$ of Eu enrichment is delayed with respect to star formation. We compare our results to predictions of several supernova and AGB yield models; C, Na, K, Mn, and Ca all show discrepancies with models that could make them useful diagnostics of nucleosynthesis physics.

A large number of dark energy and modified gravity models lead to the same expansion history of the Universe, hence, making it difficult to distinguish them from observations. To make the calculations transparent, we consider $f(R)$ gravity with a pressureless matter without making any assumption about the form of $f(R)$. Using the late-time expansion history realizations constructed by Shafieloo et al~\cite{2018-Shafieloo.etal-PRD}, we explicitly show for any $f(R)$ model that the Bardeen potentials $\Psi$ and $\Phi$ evolve differently. For an arbitrary $f(R)$ model that leads to late-time accelerated expansion, we explicitly show that $|\Psi + \Phi|$ and its time-derivative evolves differently than the $\Lambda$CDM model at lower redshifts. We show that the $\Psi/\Phi$ has a significant deviation from unity for larger wave-numbers. We discuss the implications of the results for the cosmological observations.

Analyses of data from spectroscopic and astrometric surveys have led to conflicting results concerning the vertical characteristics of the Milky Way. Ages are often used to provide clarity, but typical uncertainties of $>$ 40\,\% restrict the validity of the inferences made. Using the \textit{Kepler} APOKASC sample for context, we explore the global population trends of two K2 campaign fields (3 and 6), which extend further vertically out of the Galactic plane than APOKASC. We analyse the properties of red giant stars utilising three asteroseismic data analysis methods to cross-check and validate detections. The Bayesian inference tool PARAM is used to determine the stellar masses, radii and ages. Evidence of a pronounced red giant branch bump and an [$\alpha$/Fe] dependence on the position of the red clump is observed from the radii distribution of the K2 fields. Two peaks in the age distribution centred at $\sim$5 and and $\sim$12 Gyr are found using a sample with $\sigma_{\rm{age}}$ $<$ 35\,\%. In a comparison with \textit{Kepler}, we find the older peak to be more prominent for K2. This age bimodality is also observed based on a chemical selection of low- ($\leq$ 0.1) and high- ($>$ 0.1) [$\alpha$/Fe] stars. As a function of vertical distance from the Galactic mid-plane ($|Z|$), the age distribution shows a transition from a young to old stellar population with increasing $|Z|$ for the K2 fields. Further coverage of campaign targets with high resolution spectroscopy is required to increase the yield of precise ages achievable with asteroseismology.

The internal working of low-mass stars is of great significance to both the study of stellar structure and the history of the Milky Way. Asteroseismology has the power to directly sense the internal structure of stars and allows for the determination of the evolutionary state -- i.e. has helium burning commenced or is the energy generated only by the fusion in the hydrogen-burning shell? We use observational data from red-giant stars in a combination (known as APOKASC) of asteroseismology (from the \textit{Kepler} mission) and spectroscopy (from SDSS/APOGEE). The new feature of the analysis is that the APOKASC evolutionary state determination is based on the comparison of diverse approaches to the investigation of the frequency-power spectrum. The high level of agreement between the methods is a strong validation of the approaches. Stars for which there is not a consensus view are readily identified. The comparison also facilitates the identification of unusual stars including those that show evidence for very strong coupling between p and g cavities. The comparison between the classification based on the spectroscopic data and asteroseismic data have led to a new value for the statistical uncertainty in APOGEE temperatures. These consensus evolutionary states will be used as an input for methods that derive masses and ages for these stars based on comparison of observables with stellar evolutionary models (`grid-based modeling') and as a training set for machine-learning and other data-driven methods of evolutionary state determination

The MINiature Exoplanet Radial Velocity Array (MINERVA) is a dedicated observatory of four 0.7m robotic telescopes fiber-fed to a KiwiSpec spectrograph. The MINERVA mission is to discover super-Earths in the habitable zones of nearby stars. This can be accomplished with MINERVA's unique combination of high precision and high cadence over long time periods. In this work, we detail changes to the MINERVA facility that have occurred since our previous paper. We then describe MINERVA's robotic control software, the process by which we perform 1D spectral extraction, and our forward modeling Doppler pipeline. In the process of improving our forward modeling procedure, we found that our spectrograph's intrinsic instrumental profile is stable for at least nine months. Because of that, we characterized our instrumental profile with a time-independent, cubic spline function based on the profile in the cross dispersion direction, with which we achieved a radial velocity precision similar to using a conventional "sum-of-Gaussians" instrumental profile: 1.8 m s$^{-1}$ over 1.5 months on the RV standard star HD 122064. Therefore, we conclude that the instrumental profile need not be perfectly accurate as long as it is stable. In addition, we observed 51 Peg and our results are consistent with the literature, confirming our spectrograph and Doppler pipeline are producing accurate and precise radial velocities.

We present the discovery of KELT-24 b, a massive hot Jupiter orbiting a bright (V=8.3 mag, K=7.2 mag) young F-star with a period of 5.6 days. The host star, KELT-24 (HD 93148), has a $T_{\rm eff}$ =$6509^{+50}_{-49}$ K, a mass of $M_{*}$ = $1.460^{+0.055}_{-0.059}$ $M_{\odot}$, radius of $R_{*}$ = $1.506\pm0.022$ $R_{\odot}$, and an age of $0.78^{+0.61}_{-0.42}$ Gyr. Its planetary companion (KELT-24 b) has a radius of $R_{\rm P}$ = $1.272\pm0.021$ $R_{\rm J}$, a mass of $M_{\rm P}$ = $5.18^{+0.21}_{-0.22}$ $M_{\rm J}$, and from Doppler tomographic observations, we find that the planet's orbit is well-aligned to its host star's projected spin axis ($\lambda$ = $2.6^{+5.1}_{-3.6}$). The young age estimated for KELT-24 suggests that it only recently started to evolve from the zero-age main sequence. KELT-24 is the brightest star known to host a transiting giant planet with a period between 5 and 10 days. Although the circularization timescale is much longer than the age of the system, we do not detect a large eccentricity or significant misalignment that is expected from dynamical migration. The brightness of its host star and its moderate surface gravity make KELT-24b an intriguing target for detailed atmospheric characterization through spectroscopic emission measurements since it would bridge the current literature results that have primarily focused on lower mass hot Jupiters and a few brown dwarfs.

We examine the relationship between a number of long gamma-ray burst (lGRB) properties (isotropic emitted energy, luminosity, intrinsic duration, jet opening angle) and redshift. We find that even when accounting for conservative detector flux limits, there appears to be a significant correlation between isotropic equivalent energy and redshift, suggesting cosmological evolution of the lGRB progenitor. Analyzing a sub-sample of lGRBs with jet opening angle estimates, we find the beaming-corrected lGRB emitted energy does not correlate with redshift, but jet opening angle does. Additionally, we find a statistically significant anti-correlation between the intrinsic prompt duration and redshift, even when accounting for potential selection effects. We also find that - for a given redshift - isotropic energy is positively correlated with intrinsic prompt duration. None of these GRB properties appear to be correlated with galactic offset. From our selection-effect-corrected redshift distribution, we estimate a co-moving rate density for lGRBs, and compare this to the global cosmic star formation rate (SFR). We find the lGRB rate mildly exceeds the global star formation rate between a redshift of 3 and 5, and declines rapidly at redshifts above this (although we cannot constrain the lGRB rate above a redshift of about 6 due to sample incompleteness). We find the lGRB rate diverges significantly from the SFR at lower redshifts. We discuss both the correlations and lGRB rate density in terms of various lGRB progenitor models and their apparent preference for low-metallicity environments.

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

The problem of the origin of the elements is a fundamental one in astronomy and one that has many open questions. Prominent examples include (1) the nature of Type Ia supernovae and the timescale of their contributions; (2) the observational identification of elements such as titanium and potassium with the $\alpha$-elements in conflict with core-collapse supernova predictions; (3) the number and relative importance of r-process sites; (4) the origin of carbon and nitrogen and the influence of mixing and mass loss in winds; and (5) the origin of the intermediate elements, such as Cu, Ge, As, and Se, that bridge the region between charged-particle and neutron-capture reactions. The next decade will bring to maturity many of the new tools that have recently made their mark, such as large-scale chemical cartography of the Milky Way and its satellites, the addition of astrometric and asteroseismic information, the detection and characterization of gravitational wave events, 3-D simulations of convection and model atmospheres, and improved laboratory measurements for transition probabilities and nuclear masses. All of these areas are key for continued improvement, and such improvement will benefit many areas of astrophysics.

A small set of chemically old stars that appear young by their independently derived masses has been detected. These are so-called $\alpha$-rich young stars. For a sample of 51 red-giant stars, for which spectra are available from SDSS/APOGEE and masses are available from asteroseismic measures based on \textit{Kepler} lightcurves, we derive the C, N and O abundances through an independent analysis. These stars span a wide range of N/C surface number density ratios. We interpret the high-mass stars with low N/C as being products of mergers or mass transfer during or after first dredge up, because the dredge-up features are the same as for low-mass stars. The $\alpha$-rich young stars with high N/C follow the expected trend of N/C for their mass, and could be either genuine young stars (leaving their high [$\alpha$/Fe] unexplained) or the results of mergers on the main sequence.

Exoplanet surveys of evolved stars have provided increasing evidence that the formation of giant planets depends not only on stellar metallicity ([Fe/H]), but also the mass ($M_\star$). However, measuring accurate masses for subgiants and giants is far more challenging than it is for their main-sequence counterparts, which has led to recent concerns regarding the veracity of the correlation between stellar mass and planet occurrence. In order to address these concerns we use HIRES spectra to perform a spectroscopic analysis on an sample of 245 subgiants and derive new atmospheric and physical parameters. We also calculate the space velocities of this sample in a homogeneous manner for the first time. When reddening corrections are considered in the calculations of stellar masses and a -0.12 M$_{\odot}$ offset is applied to the results, the masses of the subgiants are consistent with their space velocity distributions, contrary to claims in the literature. Similarly, our measurements of their rotational velocities provide additional confirmation that the masses of subgiants with $M_\star \geq 1.6$ M$_{\odot}$ (the "Retired A Stars") have not been overestimated in previous analyses. Using these new results for our sample of evolved stars, together with an updated sample of FGKM dwarfs, we confirm that giant planet occurrence increases with both stellar mass and metallicity up to 2.0 M$_{\odot}$. We show that the probability of formation of a giant planet is approximately a one-to-one function of the total amount of metals in the protoplanetary disk $M_\star 10^{[Fe/H]}$. This correlation provides additional support for the core accretion mechanism of planet formation.

High redshift quasars emit copious X-ray photons which heat the intergalactic medium to temperatures up to $\sim$ 10$^6$ K. At such high temperatures the primordial gas will not form stars until it is assembled into dark matter haloes with masses of up to $\sim$ 10$^{11}$ M$_{\odot}$, at which point the hot gas collapses and cools under the influence of gravity. Once this occurs, there is a massive reservoir of primordial gas from which stars can form, potentially setting the stage for the brightest Population (Pop) III starbursts in the early Universe. Supporting this scenario, recent observations of quasars at z $\sim$ 6 have revealed a lack of accompanying Lyman $\alpha$ emitting galaxies, consistent with suppression of primordial star formation in haloes with masses below $\sim$ 10$^{10}$ M$_{\odot}$. Here we model the chemical and thermal evolution of the primordial gas as it collapses into such a massive halo irradiated by a nearby quasar in the run-up to a massive Pop III starburst. We find that within $\sim$ 100 kpc of the highest redshift quasars discovered to date the Lyman-Werner flux produced in the quasar host galaxy may be high enough to stimulate the formation of a direct collapse black hole (DCBH). A survey with single pointings of the NIRCam instrument at individual known high-z quasars may be a promising strategy for finding Pop III stars and DCBHs with the James Webb Space Telescope.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

We report the discovery of the first likely black hole in a non-interacting binary system with a field red giant. By combining radial velocity measurements from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) with photometric variability data from the All-Sky Automated Survey for Supernovae (ASAS-SN), we identified the bright rapidly-rotating giant 2MASS J05215658+4359220 as a binary system with a massive unseen companion. Subsequent radial velocity measurements reveal a system with an orbital period of 83 days and near-zero eccentricity. The photometric variability period of the giant is consistent with the orbital period, indicative of star spots and tidal synchronization. Constraints on the giant's mass and radius from its luminosity, surface gravity, and temperature imply an unseen companion with mass of $3.3^{+2.8}_{-0.7}$ M$_\odot$, indicating a low-mass black hole or an exceedingly massive neutron star. Measurement of the astrometric binary motion by {\it Gaia} will further characterize the system. This discovery demonstrates the potential of massive spectroscopic surveys like APOGEE and all-sky, high-cadence photometric surveys like ASAS-SN to revolutionize our understanding of the compact object mass function, and to test theories of binary star evolution and the supernova mechanism.

The evolution of a star is driven by the physical processes in its interior making the theory of stellar structure and evolution the most crucial ingredient for not only stellar evolution studies, but any field of astronomy which relies on the yields along stellar evolution. High-precision time-series photometric data assembled by recent space missions revealed that current models of stellar structure and evolution show major shortcomings already in the two earliest nuclear burning phases, impacting all subsequent phases prior to the formation of the end-of-life remnant. This white paper focuses specifically on the transport of chemical elements and of angular momentum in the stellar structure and evolution models of stars born with convective core, as revealed by their gravity-mode oscillations.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

Interacting binaries containing white dwarfs can lead to a variety of outcomes that range from powerful thermonuclear explosions, which are important in the chemical evolution of galaxies and as cosmological distance estimators, to strong sources of low frequency gravitational wave radiation, which makes them ideal calibrators for the gravitational low-frequency wave detector LISA mission. However, current theoretical evolution models still fail to explain the observed properties of the known populations of white dwarfs in both interacting and detached binaries. Major limitations are that the existing population models have generally been developed to explain the properties of sub-samples of these systems, occupying small volumes of the vast parameter space, and that the observed samples are severely biased. The overarching goal for the next decade is to assemble a large and homogeneous sample of white dwarf binaries that spans the entire range of evolutionary states, to obtain precise measurements of their physical properties, and to further develop the theory to satisfactorily reproduce the properties of the entire population. While ongoing and future all-sky high- and low-resolution optical spectroscopic surveys allow us to enlarge the sample of these systems, high-resolution ultraviolet spectroscopy is absolutely essential for the characterization of the white dwarfs in these binaries. The Hubble Space Telescope is currently the only facility that provides ultraviolet spectroscopy, and with its foreseeable demise, planning the next ultraviolet mission is of utmost urgency.

Chemical abundances for 15 elements (C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are presented for 83 stellar members of the 4 Gyr old solar-metallicity open cluster M67. The sample contains stars spanning a wide range of evolutionary phases, from G dwarfs to red clump stars. The abundances were derived from near-IR ($\lambda$1.5 -- 1.7$\mu$m) high-resolution spectra ($R$ = 22,500) from the SDSS-IV/APOGEE survey. A 1-D LTE abundance analysis was carried out using the APOGEE synthetic spectral libraries, via chi-square minimization of the synthetic and observed spectra with the qASPCAP code. We found significant abundance differences ($\sim$0.05 -- 0.30 dex) between the M67 member stars as a function of the stellar mass (or position on the HR diagram), where the abundance patterns exhibit a general depletion (in [X/H]) in stars at the main-sequence turnoff. The amount of the depletion is different for different elements. We find that atomic diffusion models provide, in general, good agreement with the abundance trends for most chemical species, supporting recent studies indicating that measurable atomic diffusion operates in M67 stars.

Internal mixing on the giant branch is an important process which affects the evolution of stars and the chemical evolution of the galaxy. While several mechanisms have been proposed to explain this mixing, better empirical constraints are necessary. Here, we use [C/N] abundances in 26097 evolved stars from the SDSS-IV/APOGEE-2 Data Release 14 to trace mixing and extra mixing in old field giants with -1.7< [Fe/H] < 0.1. We show that the APOGEE [C/N] ratios before any dredge-up occurs are metallicity dependent, but that the change in [C/N] at the first dredge-up is metallicity independent for stars above [Fe/H] ~ -1. We identify the position of the red giant branch (RGB) bump as a function of metallicity, note that a metallicity-dependent extra mixing episode takes place for low-metallicity stars ([Fe/H] <-0.4) 0.14 dex in log g above the bump, and confirm that this extra mixing is stronger at low metallicity, reaching $\Delta$ [C/N] = 0.58 dex at [Fe/H] = -1.4. We show evidence for further extra mixing on the upper giant branch, well above the bump, among the stars with [Fe/H] < -1.0. This upper giant branch mixing is stronger in the more metal-poor stars, reaching 0.38 dex in [C/N] for each 1.0 dex in log g. The APOGEE [C/N] ratios for red clump (RC) stars are significantly higher than for stars at the tip of the RGB, suggesting additional mixing processes occur during the helium flash or that unknown abundance zero points for C and N may exist among the red clump RC sample. Finally, because of extra mixing, we note that current empirical calibrations between [C/N] ratios and ages cannot be naively extrapolated for use in low-metallicity stars specifically for those above the bump in the luminosity function.

Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c.

Super-Earths belong to a class of planet not found in the Solar System, but which appear common in the Galaxy. Given that some super-Earths are rocky, while others retain substantial atmospheres, their study can provide clues as to the formation of both rocky planets and gaseous planets, and - in particular - they can help to constrain the role of photo-evaporation in sculpting the exoplanet population. GJ 9827 is a system already known to host 3 super-Earths with orbital periods of 1.2, 3.6 and 6.2 days. Here we use new HARPS-N radial velocity measurements, together with previously published radial velocities, to better constrain the properties of the GJ 9827 planets. Our analysis can't place a strong constraint on the mass of GJ 9827 c, but does indicate that GJ 9827 b is rocky with a composition that is probably similar to that of the Earth, while GJ 9827 d almost certainly retains a volatile envelope. Therefore, GJ 9827 hosts planets on either side of the radius gap that appears to divide super-Earths into pre-dominantly rocky ones that have radii below $\sim 1.5 R_\oplus$, and ones that still retain a substantial atmosphere and/or volatile components, and have radii above $\sim 2 R_\oplus$. That the less heavily irradiated of the 3 planets still retains an atmosphere, may indicate that photoevaporation has played a key role in the evolution of the planets in this system.

We present [C/N]-[Fe/H] abundance trends from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey, Data Release 14 (DR14), for red giant branch stars across the Milky Way Galaxy (MW, 3 kpc $<$ R $<$ 15 kpc). The carbon-to-nitrogen ratio (often expressed as [C/N]) can indicate the mass of a red giant star, from which an age can be inferred. Using masses and ages derived by Martig et al., we demonstrate that we are able to interpret the DR14 [C/N]-[Fe/H] abundance distributions as trends in age-[Fe/H] space. Our results show that an anti-correlation between age and metallicity, which is predicted by simple chemical evolution models, is not present at any Galactic zone. Stars far from the plane ($|$Z$|$ $>$ 1 kpc) exhibit a radial gradient in [C/N] ($\sim$ $-$0.04 dex/kpc). The [C/N] dispersion increases toward the plane ($\sigma_{[C/N]}$ = 0.13 at $|$Z$|$ $>$ 1 kpc to $\sigma_{[C/N]}$ = 0.18 dex at $|$Z$|$ $<$ 0.5 kpc). We measure a disk metallicity gradient for the youngest stars (age $<$ 2.5 Gyr) of $-$0.060 dex/kpc from 6 kpc to 12 kpc, which is in agreement with the gradient found using young CoRoGEE stars by Anders et al. Older stars exhibit a flatter gradient ($-$0.016 dex/kpc), which is predicted by simulations in which stars migrate from their birth radii. We also find that radial migration is a plausible explanation for the observed upturn of the [C/N]-[Fe/H] abundance trends in the outer Galaxy, where the metal-rich stars are relatively enhanced in [C/N].

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We investigate the role of angular momentum in the clustering of dark matter haloes. We make use of data from two high-resolution N-body simulations spanning over four orders of magnitude in halo mass, from $10^{9.8}$ to $10^{14}\ h^{-1}\ \text{M}_\odot$. We explore the hypothesis that mass accretion in filamentary environments alters the angular momentum of a halo, thereby driving a correlation between the spin parameter $\lambda$ and the strength of clustering. However, we do not find evidence that the distribution of matter on large scales is related to the spin of haloes. We find that a halo's spin is correlated with its age, concentration, sphericity, and mass accretion rate. Removing these correlations strongly affects the strength of secondary spin bias at low halo masses. We also find that high spin haloes are slightly more likely to be found near another halo of comparable mass. These haloes that are found near a comparable mass neighbour - a \textit{twin} - are strongly spatially biased. We demonstrate that this \textit{twin bias}, along with the relationship between spin and mass accretion rates, statistically accounts for halo spin secondary bias.

We present the discovery of 7 new planets and 8 planet candidates around subgiant stars, as additions to the known sample of planets around "retired A stars" (Johnson et al. 2006). Among these are the possible first 3-planet systems around subgiant stars, HD 163607 and HD 4917. Additionally, we present calculations of possible transit times, durations, depths, and probabilities for all known planets around subgiant (3 < log g < 4) stars, focused on possible transits during the TESS mission. While most have transit probabilities of 1-2%, we find that there are 3 planets with transit probabilities > 9%.

Jupiter-mass planets with large semi-major axes ($a > 1.0$ AU) occur at a higher rate around evolved intermediate mass stars. There is a pronounced paucity of close-in ($a < 0.6$ AU), intermediate period ($5 < P < 100$ days), low-mass ($M_{\rm planet} < 0.7M_{\rm Jup} $) planets, known as the `Planet Desert'. Current radial velocity methods have yet to yield close-in, low-mass planets around these stars because the planetary signals could be hidden by the (5-10) m s$^{-1}$ radial velocity variations caused by acoustic oscillations. We find that by implementing an observing strategy of taking three observations per night separated by an optimal $\Delta t$, which is a function of the oscillation periods and amplitudes, we can average over the stellar jitter and improve our sensitivity to low-mass planets. We find $\Delta t$ can be approximated using the stellar mass and radius given by the relationship $\Delta t = $1.79 $(M/M_{\odot})^{-0.82} ~(R/R_{\odot})^{1.92}$. We test our proposed method by injecting planets into very well sampled data of a subgiant star, $\gamma$ Cep. We compare the fraction of planets recovered by our method to the fraction of planets recovered using current radial velocity observational strategies. We find that our method decreases the RMS of the stellar jitter due to acoustic oscillations by a factor of three over current single epoch observing strategies used for subgiant stars. Our observing strategy provides a means to test whether the Planet Desert extends to lower mass planets.

We present an overview of the distributions of 11 elemental abundances in the Milky Way's inner regions, as traced by APOGEE stars released as part of SDSS Data Release 14/15 (DR14/DR15), including O, Mg, Si, Ca, Cr, Mn, Co, Ni, Na, Al, and K. This sample spans ~4000 stars with R_GC<4 kpc, enabling the most comprehensive study to date of these abundances and their variations within the innermost few kiloparsecs of the Milky Way. We describe the observed abundance patterns ([X/Fe]-[Fe/H]), compare to previous literature results and to patterns in stars at the solar Galactic radius, and discuss possible trends with DR14/DR15 effective temperatures. We find that the position of the [Mg/Fe]-[Fe/H] "knee" is nearly constant with R_GC, indicating a well-mixed star-forming medium or high levels of radial migration in the early inner Galaxy. We quantify the linear correlation between pairs of elements in different subsamples of stars and find that these relationships vary; some abundance correlations are very similar between the alpha-rich and alpha-poor stars, but others differ significantly, suggesting variations in the metallicity dependencies of certain supernova yields. These empirical trends will form the basis for more detailed future explorations and for the refinement of model comparison metrics. That the inner Milky Way abundances appear dominated by a single chemical evolutionary track and that they extend to such high metallicities underscore the unique importance of this part of the Galaxy for constraining the ingredients of chemical evolution modeling and for improving our understanding of the evolution of the Galaxy as a whole.

Recent work indicates that the nearby Galactic halo is dominated by the debris from a major accretion event. We confirm that result from an analysis of APOGEE-DR14 element abundances and $\textit{Gaia}$-DR2 kinematics of halo stars. We show that $\sim$2/3 of nearby halo stars have high orbital eccentricities ($e \gtrsim 0.8$), and abundance patterns typical of massive Milky Way dwarf galaxy satellites today, characterised by relatively low [Fe/H], [Mg/Fe], [Al/Fe], and [Ni/Fe]. The trend followed by high $e$ stars in the [Mg/Fe]-[Fe/H] plane shows a change of slope at [Fe/H]$\sim-1.3$, which is also typical of stellar populations from relatively massive dwarf galaxies. Low $e$ stars exhibit no such change of slope within the observed [Fe/H] range and show slightly higher abundances of Mg, Al and Ni. Unlike their low $e$ counterparts, high $e$ stars show slightly retrograde motion, make higher vertical excursions and reach larger apocentre radii. By comparing the position in [Mg/Fe]-[Fe/H] space of high $e$ stars with those of accreted galaxies from the EAGLE suite of cosmological simulations we constrain the mass of the accreted satellite to be in the range $10^{8.5}\lesssim M_*\lesssim 10^{9}\mathrm{M_\odot}$. We show that the median orbital eccentricities of debris are largely unchanged since merger time, implying that this accretion event likely happened at $z\lesssim1.5$. The exact nature of the low $e$ population is unclear, but we hypothesise that it is a combination of $\textit{in situ}$ star formation, high $|z|$ disc stars, lower mass accretion events, and contamination by the low $e$ tail of the high $e$ population. Finally, our results imply that the accretion history of the Milky Way was quite unusual.

We map the trends of elemental abundance ratios across the Galactic disk, spanning R = 3-15 kpc and midplane distance |Z|= 0-2 kpc, for 15 elements in a sample of 20,485 stars measured by the SDSS/APOGEE survey (O, Na, Mg, Al, Si, P, S, K, Ca, V, Cr, Mn, Fe, Co, Ni). Adopting Mg rather than Fe as our reference element, and separating stars into two populations based on [Fe/Mg], we find that the median trends of [X/Mg] vs. [Mg/H] in each population are nearly independent of location in the Galaxy. The full multi-element cartography can be summarized by combining these nearly universal median sequences with our measured metallicity distribution functions and the relative proportions of the low-[Fe/Mg] (high-alpha) and high-[Fe/Mg] (low-alpha) populations, which depend strongly on R and |Z|. We interpret the median sequences with a semi-empirical "2-process" model that describes both the ratio of core collapse and Type Ia supernova contributions to each element and the metallicity dependence of the supernova yields. These observationally inferred trends can provide strong tests of supernova nucleosynthesis calculations. Our results lead to a relatively simple picture of abundance ratio variations in the Milky Way, in which the trends at any location can be described as the sum of two components with relative contributions that change systematically and smoothly across the Galaxy. Deviations from this picture and future extensions to other elements can provide further insights into the physics of stellar nucleosynthesis and unusual events in the Galaxy's history.

This paper reports on the validation and mass measurement of K2-263b, a sub-Neptune orbiting a quiet G9V star. Using K2 data from campaigns C5 and C16, we find this planet to have a period of $50.818947\pm 0.000094$ days and a radius of $2.41\pm0.12$ R$_{\oplus}$. We followed this system with HARPS-N to obtain 67 precise radial velocities. A combined fit of the transit and radial velocity data reveals that K2-263b has a mass of $14.8\pm3.1$ M$_{\oplus}$. Its bulk density ($5.7_{-1.4}^{+1.6}$ g cm$^{-3}$) implies that this planet has a significant envelope of water or other volatiles around a rocky core. EPIC211682544b likely formed in a similar way as the cores of the four giant planets in our own Solar System, but for some reason, did not accrete much gas. The planetary mass was confirmed by an independent Gaussian process-based fit to both the radial velocities and the spectroscopic activity indicators. K2-263b belongs to only a handful of confirmed K2 exoplanets with periods longer than 40 days. It is among the longest periods for a small planet with a precisely determined mass using radial velocities.

Supermassive black holes (SMBHs) are the central engines of luminous quasars and are found in most massive galaxies today. But the recent discoveries of ULAS J1120+0641, a $2 \times 10^9$ M$_{\odot}$ BH at $z =$ 7.1, and ULAS J1342+0928, a $8.0 \times 10^{8}$ M$_{\odot}$ BH at $z =$ 7.5, now push the era of quasar formation up to just 690 Myr after the Big Bang. Here we report new cosmological simulations of SMBHs with X-rays fully coupled to primordial chemistry and hydrodynamics that show that J1120 and J1342 can form from direct collapse black holes (DCBHs) if their growth is fed by cold, dense accretion streams, like those thought to fuel rapid star formation in some galaxies at later epochs. Our models reproduce all of the observed properties of J1120: its mass, luminosity, and H II region as well as star formation rates and metallicities in its host galaxy. They also reproduce the dynamical mass of the innermost 1.5 kpc of its emission region recently measured by ALMA and J-band magnitudes that are in good agreement with those found by the VISTA Hemisphere Survey.

We describe the detection of water vapor in the atmosphere of the transiting hot Jupiter KELT-2Ab by treating the star-planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet's motion with deep combined flux observations of the star and the planet. In total, six epochs of Keck NIRSPEC $L$-band observations were obtained, and the full data set was subjected to a cross correlation analysis with a grid of self-consistent atmospheric models. We measure a radial projection of the Keplerian velocity, $K_P$, of 148 $\pm$ 7 km s$^{-1}$, consistent with transit measurements, and detect water vapor at 3.8$\sigma$. We combine NIRSPEC $L$-band data with $Spitzer$ IRAC secondary eclipse data to further probe the metallicity and carbon-to-oxygen ratio of KELT-2Ab's atmosphere. While the NIRSPEC analysis provides few extra constraints on the $Spitzer$ data, it does provide roughly the same constraints on metallicity and carbon-to-oxygen ratio. This bodes well for future investigations of the atmospheres of non-transiting hot Jupiters.

The InfraRed Imaging Spectrograph (IRIS) is the first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope (TMT). IRIS includes three natural guide star (NGS) On-Instrument Wavefront Sensors (OIWFS) to measure tip/tilt and focus errors in the instrument focal plane. NFIRAOS also has an internal natural guide star wavefront sensor, and IRIS and NFIRAOS must precisely coordinate the motions of their wavefront sensor positioners to track the locations of NGSs while the telescope is dithering (offsetting the telescope to cover more area), to avoid a costly re-acquisition time penalty. First, we present an overview of the sequencing strategy for all of the involved subsystems. We then predict the motion of the telescope during dithers based on finite-element models provided by TMT, and finally analyze latency and jitter issues affecting the propagation of position demands from the telescope control system to individual motor controllers.

Data and analysis methodology used for the SDSS/APOGEE Data Releases 13 and 14 are described, highlighting differences from the DR12 analysis presented in Holtzman (2015). Some improvement in the handling of telluric absorption and persistence is demonstrated. The derivation and calibration of stellar parameters, chemical abundances, and respective uncertainties are described, along with the ranges over which calibration was performed. Some known issues with the public data related to the calibration of the effective temperatures (DR13), surface gravity (DR13 and DR14), and C and N abundances for dwarfs (DR13 and DR14) are highlighted. We discuss how results from a data-driven technique, The Cannon (Casey 2016), are included in DR14, and compare those with results from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). We describe how using The Cannon in a mode that restricts the abundance analysis of each element to regions of the spectrum with known features from that element leads to Cannon abundances can lead to significantly different results for some elements than when all regions of the spectrum are used to derive abundances.

We consider an $f(R)$ action that is non-minimally coupled to a massive scalar field. The model closely resembles scalar-tensor theory and by conformal transformation can be transformed to Einstein frame. To avoid the ambiguity of the frame dependence, we obtain an exact analytical solution in Jordan frame and show that the model leads to a period of accelerated expansion with an exit. Further, we compute the scalar and tensor power spectrum for the model and compare them with observations.

The first detailed chemical abundance analysis of the M dwarf (M4.0) exoplanet-hosting star Ross 128 is presented here, based upon near-infrared (1.5--1.7 \micron) high-resolution ($R$$\sim$22,500) spectra from the SDSS-APOGEE survey. We determined precise atmospheric parameters $T_{\rm eff}$=3231$\pm$100K, log$g$=4.96$\pm$0.11 dex and chemical abundances of eight elements (C, O, Mg, Al, K, Ca, Ti, and Fe), finding Ross 128 to have near solar metallicity ([Fe/H] = +0.03$\pm$0.09 dex). The derived results were obtained via spectral synthesis (1-D LTE) adopting both MARCS and PHOENIX model atmospheres; stellar parameters and chemical abundances derived from the different adopted models do not show significant offsets. Mass-radius modeling of Ross 128b indicate that it lies below the pure rock composition curve, suggesting that it contains a mixture of rock and iron, with the relative amounts of each set by the ratio of Fe/Mg. If Ross 128b formed with a sub-solar Si/Mg ratio, and assuming the planet's composition matches that of the host-star, it likely has a larger core size relative to the Earth. The derived planetary parameters -- insolation flux (S$_{\rm Earth}$=1.79$\pm$0.26) and equilibrium temperature ($T_{\rm eq}$=294$\pm$10K) -- support previous findings that Ross 128b is a temperate exoplanet in the inner edge of the habitable zone.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present a catalog of stellar properties for a large sample of 6676 evolved stars with APOGEE spectroscopic parameters and \textit{Kepler} asteroseismic data analyzed using five independent techniques. Our data includes evolutionary state, surface gravity, mean density, mass, radius, age, and the spectroscopic and asteroseismic measurements used to derive them. We employ a new empirical approach for combining asteroseismic measurements from different methods, calibrating the inferred stellar parameters, and estimating uncertainties. With high statistical significance, we find that asteroseismic parameters inferred from the different pipelines have systematic offsets that are not removed by accounting for differences in their solar reference values. We include theoretically motivated corrections to the large frequency spacing ($\Delta \nu$) scaling relation, and we calibrate the zero point of the frequency of maximum power ($\nu_{\rm max}$) relation to be consistent with masses and radii for members of star clusters. For most targets, the parameters returned by different pipelines are in much better agreement than would be expected from the pipeline-predicted random errors, but 22\% of them had at least one method not return a result and a much larger measurement dispersion. This supports the usage of multiple analysis techniques for asteroseismic stellar population studies. The measured dispersion in mass estimates for fundamental calibrators is consistent with our error model, which yields median random and systematic mass uncertainties for RGB stars of order 4\%. Median random and systematic mass uncertainties are at the 9\% and 8\% level respectively for RC stars.

Many field stars reside in binaries, and the analysis and interpretation of photometric and spectroscopic surveys must take this into account. We have developed a model to predict how binaries influence the scientific results inferred from large spectroscopic surveys. Based on the rapid binary evolution code BSE, it allows us to model a representative population of binaries and generate synthetic survey observations. We describe this model in detail, and apply it to the radial velocity variation of subgiant and giant stars in the Galactic disc, as observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III. APOGEE provides an excellent data set for testing our binary models since a large fraction of the stars have been observed multiple times. By comparing our model to the APOGEE observations we constrain the initial binary fraction of solar-metallicity stars in the sample to be $f_{\rm b,0}=0.35\pm0.01$, in line with the solar neighbourhood. We find that the binary fraction is higher at lower metallicities, consistent with other observational studies. Our model matches the shape of the high-velocity scatter in APOGEE, which suggests that most velocity variability above 0.5 km/s comes from binaries. Our exploration of binary initial properties shows that APOGEE is mostly sensitive to binaries with periods between 3 and 3000 years, and is largely insensitive to the detailed properties of the population. We can, however, rule out a population where the mass of the lower-mass star is drawn from the IMF independently of its companion.

We report on a program to monitor classical novae (CNe) to determine if they produced dust in the ejecta created by their outbursts. Of the ten systems we followed, five produced dust. We also present limited infrared and optical spectroscopy of these objects. We present a complete $JHKLM$ spectrum for V992 Sco. V992 Sco was one of the brightest CNe in the infrared of all time, and our $M$-band spectrum of this object shows strong emission from the CO fundamental. We believe this to be the first, and only, spectroscopic observation of this feature in a CNe.

We present the confirmation of a small, moderately-irradiated (F = 155 +/- 7 Fearth) Neptune with a substantial gas envelope in a P=11.8728787+/-0.0000085-day orbit about a quiet, Sun-like G0V star Kepler-1655. Based on our analysis of the Kepler light curve, we determined Kepler-1655b's radius to be 2.213+/-0.082 Rearth. We acquired 95 high-resolution spectra with TNG/HARPS-N, enabling us to characterize the host star and determine an accurate mass for Kepler-1655b of 5.0+3.1/-2.8 Mearth via Gaussian-process regression. Our mass determination excludes an Earth-like composition with 98\% confidence. Kepler-1655b falls on the upper edge of the evaporation valley, in the relatively sparsely occupied transition region between rocky and gas-rich planets. It is therefore part of a population of planets that we should actively seek to characterize further.

We present the discovery of KELT-22Ab, a hot Jupiter from the KELT-South survey. KELT-22Ab transits the moderately bright ($V\sim 11.1$) Sun-like G2V star TYC 7518-468-1. The planet has an orbital period of $P = 1.3866529 \pm 0.0000027 $ days, a radius of $R_{P} = 1.285_{-0.071}^{+0.12}~R_{J}$, and a relatively large mass of $M_{P} = 3.47_{-0.14}^{+0.15}~ M_{J}$. The star has $R_{\star} = 1.099_{-0.046}^{+0.079}~ R_{\odot}$, $M_{\star} = 1.092_{-0.041}^{+0.045}~ M_{\odot}$, ${T_{\rm eff}\,} = 5767_{-49}^{+50}~$ K, ${\log{g_\star}} = 4.393_{-0.060}^{+0.039}~$ (cgs), and [m/H] = $+0.259_{-0.083}^{+0.085}~$, and thus, other than its slightly super-solar metallicity, appears to be a near solar twin. Surprisingly, KELT-22A exhibits kinematics and a Galactic orbit that are somewhat atypical for thin disk stars. Nevertheless, the star is rotating quite rapidly for its estimated age, shows evidence of chromospheric activity, and is somewhat metal rich. Imaging reveals a slightly fainter companion to KELT-22A that is likely bound, with a projected separation of 6\arcsec ($\sim$1400 AU). In addition to the orbital motion caused by the transiting planet, we detect a possible linear trend in the radial velocity of KELT-22A suggesting the presence of another relatively nearby body that is perhaps non-stellar. KELT-22Ab is highly irradiated (as a consequence of the small semi-major axis of $a/R_{\star} = 4.97$), and is mildly inflated. At such small separations, tidal forces become significant. The configuration of this system is optimal for measuring the rate of tidal dissipation within the host star. Our models predict that, due to tidal forces, the semi-major axis of KELT-22Ab is decreasing rapidly, and is thus predicted to spiral into the star within the next Gyr.

Detailed chemical abundance distributions for fourteen elements are derived for eight high-probability stellar members of the solar metallicity old open cluster M67 with an age of $\sim$4 Gyr. The eight stars consist of four pairs, with each pair occupying a distinct phase of stellar evolution: two G-dwarfs, two turnoff stars, two G-subgiants, and two red clump K-giants. The abundance analysis uses near-IR high-resolution spectra ($\lambda$1.5 -- 1.7$\mu$m) from the APOGEE survey and derives abundances for C, N, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe. Our derived stellar parameters and metallicity for 2M08510076+113115 suggest that this star is a solar-twin, exhibiting abundance differences relative to the Sun of $\leq$ 0.04 dex for all elements. Chemical homogeneity is found within each class of stars ($\sim$0.02 dex), while significant abundance variations ($\sim$0.05 -- 0.20 dex) are found across the different evolutionary phases; the turnoff stars typically have the lowest abundances, while the red clump tend to have the largest. Non-LTE corrections to the LTE-derived abundances are unlikely to explain the differences. A detailed comparison of the derived Fe, Mg, Si, and Ca abundances with recently published surface abundances from stellar models that include chemical diffusion, provides a good match between the observed and predicted abundances as a function of stellar mass. Such agreement would indicate the detection of chemical diffusion processes in the stellar members of M67.

M-dwarf stars are promising targets for identifying and characterizing potentially habitable planets. K2-3 is a nearby (45 pc), early-type M dwarf hosting three small transiting planets, the outermost of which orbits close to the inner edge of the stellar (optimistic) habitable zone. The K2-3 system is well suited for follow-up characterization studies aimed at determining accurate masses and bulk densities of the three planets. Using a total of 329 radial velocity measurements collected over 2.5 years with the HARPS-N and HARPS spectrographs and a proper treatment of the stellar activity signal, we aim to improve measurements of the masses and bulk densities of the K2-3 planets. We use our results to investigate the physical structure of the planets. We analyse radial velocity time series extracted with two independent pipelines by using Gaussian process regression. We adopt a quasi-periodic kernel to model the stellar magnetic activity jointly with the planetary signals. We use Monte Carlo simulations to investigate the robustness of our mass measurements of K2-3\,c and K2-3\,d, and to explore how additional high-cadence radial velocity observations might improve them. Despite the stellar activity component being the strongest signal present in the radial velocity time series, we are able to derive masses for both planet b ($M_{\rm b}=6.6\pm1.1$ $M_{\rm \oplus}$) and planet c ($M_{\rm c}=3.1^{+1.3}_{-1.2}$ $M_{\rm \oplus}$). The Doppler signal due to K2-3\,d remains undetected, likely because of its low amplitude compared to the radial velocity signal induced by the stellar activity. The closeness of the orbital period of K2-3\,d to the stellar rotation period could also make the detection of the planetary signal complicated. [...]

Probing the connection between a star's metallicity and the presence and properties of any associated planets offers an observational link between conditions during the epoch of planet formation and mature planetary systems. We explore this connection by analyzing the metallicities of Kepler target stars and the subset of stars found to host transiting planets. After correcting for survey incompleteness, we measure planet occurrence: the number of planets per 100 stars with a given metallicity $M$. Planet occurrence correlates with metallicity for some, but not all, planet sizes and orbital periods. For warm super-Earths having $P = 10-100$ days and $R_P = 1.0-1.7~R_E$, planet occurrence is nearly constant over metallicities spanning $-$0.4 dex to +0.4 dex. We find 20 warm super-Earths per 100 stars, regardless of metallicity. In contrast, the occurrence of warm sub-Neptunes ($R_P = 1.7-4.0~R_E$) doubles over that same metallicity interval, from 20 to 40 planets per 100 stars. We model the distribution of planets as $d f \propto 10^{\beta M} d M$, where $\beta$ characterizes the strength of any metallicity correlation. This correlation steepens with decreasing orbital period and increasing planet size. For warm super-Earths $\beta = -0.3^{+0.2}_{-0.2}$, while for hot Jupiters $\beta = +3.4^{+0.9}_{-0.8}$. High metallicities in protoplanetary disks may increase the mass of the largest rocky cores or the speed at which they are assembled, enhancing the production of planets larger than 1.7 $R_E$. The association between high metallicity and short-period planets may reflect disk density profiles that facilitate the inward migration of solids or higher rates of planet-planet scattering.

Hot Jupiters were the first exoplanets to be discovered around main sequence stars and astonished us with their close-in orbits. They are a prime example of how exoplanets have challenged our textbook, solar-system inspired story of how planetary systems form and evolve. More than twenty years after the discovery of the first hot Jupiter, there is no consensus on their predominant origin channel. Three classes of hot Jupiter creation hypotheses have been proposed: in situ formation, disk migration, and high eccentricity tidal migration. Although no origin channel alone satisfactorily explains all the evidence, two major origins channels together plausibly account for properties of hot Jupiters themselves and their connections to other exoplanet populations.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a chronology of the evolution of the disk throughout the age of the Galaxy. We chemically dissect the Milky Way disk population using a sample of red giant stars spanning out to 2~kpc in the solar annulus observed by the {\it Kepler} satellite, with the added dimension of asteroseismic ages. Our results reveal a clear difference in age between the low- and high-$\alpha$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$\alpha$/Fe] for the high-$\alpha$ disk stars. Our results indicate that this component formed over a period of more than 2~Gyr with a wide range of [M/H] and [$\alpha$/Fe] independent of time. Our findings show that the kinematic properties of young $\alpha$-rich stars are consistent with the rest of the high-$\alpha$ population and different from the low-$\alpha$ stars of similar age, rendering support to their origin being old stars that went through a mass transfer or stellar merger event, making them appear younger, instead of migration of truly young stars formed close to the Galactic bar.

Ultra-short period (USP) planets are a class of low mass planets with periods shorter than one day. Their origin is still unknown, with photo-evaporation of mini-Neptunes and in-situ formation being the most credited hypotheses. Formation scenarios differ radically in the predicted composition of USP planets, it is therefore extremely important to increase the still limited sample of USP planets with precise and accurate mass and density measurements. We report here the characterization of an USP planet with a period of 0.28 days around K2-141 (EPIC 246393474), and the validation of an outer planet with a period of 7.7 days in a grazing transit configuration. We derived the radii of the planets from the K2 light curve and used high-precision radial velocities gathered with the HARPS-N spectrograph for mass measurements. For K2-141b we thus inferred a radius of $1.51\pm0.05~R_\oplus$ and a mass of $5.08\pm0.41~M_\oplus$, consistent with a rocky composition and lack of a thick atmosphere. K2-141c is likely a Neptune-like planet, although due to the grazing transits and the non-detection in the RV dataset, we were not able to put a strong constraint on its density. We also report the detection of secondary eclipses and phase curve variations for K2-141b. The phase variation can be modeled either by a planet with a geometric albedo of $0.30 \pm 0.06$ in the Kepler bandpass, or by thermal emission from the surface of the planet at $\sim$3000K. Only follow-up observations at longer wavelengths will allow us to distinguish between these two scenarios.

We present the discovery of the giant planet KELT-19Ab, which transits the moderately bright $(\mathrm{V} \sim 9.9)$ A8V star TYC 764-1494-1 with an orbital period of 4.61 days. We confirm the planetary nature of the companion via a combination of radial velocities, which limit the mass to $< 4.1\,\mathrm{M_J}$ $(3\sigma)$, and a clear Doppler tomography signal, which indicates a retrograde projected spin-orbit misalignment of $\lambda = -179.7^{+3.7}_{-3.8}$ degrees. Global modeling indicates that the $\rm{T_{eff}} =7500 \pm 110\,\mathrm{K}$ host star has $\mathrm{M_*} = 1.62^{+0.25}_{-0.20}\,\mathrm{M_\odot}$ and $\mathrm{R_*} = 1.83 \pm 0.10\,\mathrm{R_\odot}$. The planet has a radius of $\mathrm{R_P}=1.91 \pm 0.11\,\mathrm{R_J}$ and receives a stellar insolation flux of $\sim 3.2\times 10^{9}\,\mathrm{erg\,s^{-1}\,cm^{-2}}$, leading to an inferred equilibrium temperature of $\rm{T_{EQ}} = \sim 1935\,\rm{K}$ assuming zero albedo and complete heat redistribution. With a $v\sin{I_*}=84.8\pm 2.0\,\mathrm{km\,s^{-1}}$, the host is relatively slowly rotating compared to other stars with similar effective temperatures, and it appears to be enhanced in metallic elements but deficient in calcium, suggesting that it is likely an Am star. KELT-19A would be the first detection of an Am host of a transiting planet of which we are aware. Adaptive optics observations of the system reveal the existence of a companion with late G9V/early K1V spectral type at a projected separation of $\approx 160\,\mathrm{AU}$. Radial velocity measurements indicate that this companion is bound. Most Am stars are known to have stellar companions, which are often invoked to explain the relatively slow rotation of the primary. In this case, the stellar companion is unlikely to have caused the tidal braking of the primary. However, it may have emplaced the transiting planetary companion via the Kozai-Lidov mechanism.

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.

(Abridged) We present the first APOKASC catalog of spectroscopic and asteroseismic data for 415 dwarfs and subgiants. Asteroseismic data have been obtained by Kepler in short cadence. The spectroscopic parameters are based on spectra taken as part of APOGEE and correspond to DR13 of SDSS. We analyze our data using two Teff scales, the spectroscopic values from DR13 and those derived from SDSS griz photometry. We use the differences in our results arising from these choices as a test of systematic Teff, and find that they can lead to significant differences in the derived stellar properties. Determinations of surface gravity ($\log{g}$), mean density ($\rho$), radius ($R$), mass ($M$), and age ($\tau$) for the whole sample have been carried out with stellar grid-based modeling. We have assessed random and systematic error sources in the spectroscopic and seismic data, as well as in the grid-based modeling determination of the stellar quantities in the catalog. We provide stellar properties for both Teff scales. The median combined (random and systematic) uncertainties are 2% (0.01 dex; $\log{g}$), 3.4% ($\rho$), 2.6% ($R$), 5.1% ($M$), and 19% ($\tau$) for the photometric Teff scale and 2% ($\log{g}$), 3.5% ($\rho$), 2.7% ($R$), 6.3% ($M$), and 23% ($\tau$) for the spectroscopic scale. Comparisons with stellar quantities in the catalog by Chaplin et al.(2014) highlight the importance of metallicity measurements for determining stellar parameters accurately. We compare our results with those from other sources, including stellar radii determined from TGAS parallaxes and asteroseismic analyses based on individual frequencies. We find a very good agreement in all cases. Comparisons give strong support to the determination of stellar quantities based on global seismology, a relevant result for future missions such as TESS and PLATO. Table 5 corrected (wrongly listed SDSS Teff before).

We derive water vapor column abundances and aerosol properties from Mars Science Laboratory (MSL) ChemCam passive mode observations of scattered sky light. Each ChemCam passive sky observation acquires spectra at two different elevation angles. We fit these spectra with a discrete-ordinates multiple scattering radiative transfer model, using the correlated-k approximation for gas absorption bands. The retrieval proceeds by first fitting the continuum of the ratio of the two elevation angles to solve for aerosol properties, and then fitting the continuum-removed ratio to solve for gas abundances. The final step of the retrieval makes use of the observed CO2 absorptions and the known CO2 abundance to correct the retrieved water vapor abundance for the effects of the vertical distribution of scattering aerosols and to derive an aerosol scale height parameter. The ChemCam-retrieved water abundances show, with only a few exceptions, the same seasonal behavior and the same timing of seasonal minima and maxima as the TES, CRISM, and REMS-H data sets that we compare them to. However ChemCam-retrieved water abundances are generally lower than zonal and regional scale from-orbit water vapor data, while at the same time being significantly larger than pre-dawn REMS-H abundances. Pending further analysis of REMS-H volume mixing ratio uncertainties, the differences between ChemCam and REMS-H pre-dawn mixing ratios appear to be much too large to be explained by large scale circulations and thus they tend to support the hypothesis of substantial diurnal interactions of water vapor with the surface. Our preliminary aerosol results, meanwhile, show the expected seasonal pattern in dust particle size but also indicate a surprising inter-annual increase in water-ice cloud opacities.

We report the discovery of a new ultra-short-period planet and summarize the properties of all such planets for which the mass and radius have been measured. The new planet, EPIC~228732031b, was discovered in {\it K2} Campaign 10. It has a radius of 1.81$^{+0.16}_{-0.12}~R_{\oplus}$ and orbits a G dwarf with a period of 8.9 hours. Radial velocities obtained with Magellan/PFS and TNG/HARPS-N show evidence for stellar activity along with orbital motion. We determined the planetary mass using two different methods: (1) the "floating chunk offset" method, based only on changes in velocity observed on the same night, and (2) a Gaussian process regression based on both the radial-velocity and photometric time series. The results are consistent and lead to a mass measurement of $6.5 \pm 1.6~M_{\oplus}$, and a mean density of $6.0^{+3.0}_{-2.7}$~g~cm$^{-3}$.

We present precise radial velocity observations of WASP-47, a star known to host a hot Jupiter, a distant Jovian companion, and, uniquely, two additional transiting planets in short-period orbits: a super-Earth in a ~19 hour orbit, and a Neptune in a ~9 day orbit. We analyze our observations from the HARPS-N spectrograph along with previously published data to measure the most precise planet masses yet for this system. When combined with new stellar parameters and reanalyzed transit photometry, our mass measurements place strong constraints on the compositions of the two small planets. We find unlike most other ultra-short-period planets, the inner planet, WASP-47 e, has a mass (6.83 +/- 0.66 Me) and radius (1.810 +/- 0.027 Re) inconsistent with an Earth-like composition. Instead, WASP-47 e likely has a volatile-rich envelope surrounding an Earth-like core and mantle. We also perform a dynamical analysis to constrain the orbital inclination of WASP-47 c, the outer Jovian planet. This planet likely orbits close to the plane of the inner three planets, suggesting a quiet dynamical history for the system. Our dynamical constraints also imply that WASP-47 c is much more likely to transit than a geometric calculation would suggest. We calculate a transit probability for WASP-47 c of about 10%, more than an order of magnitude larger than the geometric transit probability of 0.6%.

Planets in the mass range from 2 to 15 MEarth are very diverse. Some of them have low densities, while others are very dense. By measuring the masses and radii, the mean densities, structure, and composition of the planets are constrained. These parameters also give us important information about their formation and evolution, and about possible processes for atmospheric loss.We determined the masses, radii, and mean densities for the two transiting planets orbiting K2-106. The inner planet has an ultra-short period of 0.57 days. The period of the outer planet is 13.3 days.Although the two planets have similar masses, their densities are very different. For K2-106b we derive Mb=8.36-0.94+0.96 MEarh, Rb=1.52+/-0.16 REarth, and a high density of 13.1-3.6+5.4 gcm-3. For K2-106c, we find Mc=5.8-3.0+3.3 MEarth, Rc=2.50-0.26+0.27 REarth and a relatively low density of 2.0-1.1+1.6 gcm-3.Since the system contains two planets of almost the same mass, but different distances from the host star, it is an excellent laboratory to study atmospheric escape. In agreement with the theory of atmospheric-loss processes, it is likely that the outer planet has a hydrogen-dominated atmosphere. The mass and radius of the inner planet is in agreement with theoretical models predicting an iron core containing 80+20-30% of its mass. Such a high metal content is surprising, particularly given that the star has an ordinary (solar) metal abundance. We discuss various possible formation scenarios for this unusual planet.

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.

APOGEE-2 is a high-resolution, near-infrared spectroscopic survey observing roughly 300,000 stars across the entire sky. It is the successor to APOGEE and is part of the Sloan Digital Sky Survey IV (SDSS-IV). APOGEE-2 is expanding upon APOGEE's goals of addressing critical questions of stellar astrophysics, stellar populations, and Galactic chemodynamical evolution using (1) an enhanced set of target types and (2) a second spectrograph at Las Campanas Observatory in Chile. APOGEE-2 is targeting red giant branch (RGB) and red clump (RC) stars, RR Lyrae, low-mass dwarf stars, young stellar objects, and numerous other Milky Way and Local Group sources across the entire sky from both hemispheres. In this paper, we describe the APOGEE-2 observational design, target selection catalogs and algorithms, and the targeting-related documentation included in the SDSS data releases.

The Kepler Space Telescope is currently searching for planets transiting stars along the ecliptic plane as part of its extended K2 mission. We processed the publicly released data from the first year of K2 observations (Campaigns 0, 1, 2, and 3) and searched for periodic eclipse signals consistent with planetary transits. Out of 59,174 targets we searched, we detect 234 planetary candidates around 208 stars. These candidates range in size from gas giants to smaller than the Earth, and range in orbital periods from hours to over a month. We conducted initial reconnaissance spectroscopy of 68 of the brighter candidate host stars, and present high resolution optical spectra for these stars. We make all of our data products, including light curves, spectra, and vetting diagnostics available to users online.

We present our models of the effect of binaries on high-resolution spectroscopic surveys. We want to determine how many binary stars will be observed, whether unresolved binaries will contaminate measurements of chemical abundances, and how we can use spectroscopic surveys to better constrain the population of binary stars in the Galaxy. Using a rapid binary-evolution algorithm that enables modelling of the most complex binary systems we generate a series of large binary populations in the Galactic disc and evaluate the results. As a first application we use our model to study the binary fraction in APOGEE giants. We find tentative evidence for a change in binary fraction with metallicity.

The four main findings about the age and abundance structure of the Milky Way bulge based on microlensed dwarf and subgiant stars are: (1) a wide metallicity distribution with distinct peaks at [Fe/H]=-1.09, -0.63, -0.20, +0.12, +0.41; (2) a high fraction of intermediate-age to young stars where at [Fe/H]>0 more than 35 % are younger than 8 Gyr, (3) several episodes of significant star formation in the bulge 3, 6, 8, and 11 Gyr ago; (4) the `knee' in the alpha-element abundance trends of the sub-solar metallicity bulge appears to be located at a slightly higher [Fe/H] (about 0.05 to 0.1 dex) than in the local thick disk.

We present a detailed elemental abundance study of 90 F and G dwarf, turn-off and subgiant stars in the Galactic bulge. Based on high-resolution spectra acquired during gravitational microlensing events, stellar ages and abundances for 11 elements (Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Zn, Y and Ba) have been determined. We find that the Galactic bulge has a wide metallicity distribution with significant peaks at [Fe/H]=-1.09, -0.63, -0.20, +0.12, +0.41. We also find a high fraction of intermediate-age to young stars: at [Fe/H]>0 more than 35 % are younger than 8 Gyr. For [Fe/H]<-0.5 most stars are 10 Gyr or older. We have also identified several episodes when significant star formation in the bulge happened: 3, 6, 8, and 12 Gyr ago. We further find that the "knee" in the alpha-element abundance trends of the sub-solar metallicity bulge is located at about 0.1 dex higher [Fe/H] than in the local thick disk. The Galactic bulge has complex age and abundance properties that appear to be tightly connected to the main Galactic stellar populations. In particular, the peaks in the metallicity distribution, the star formation episodes, and the abundance trends, show similarities with the properties of the Galactic thin and thick disks. At the same time there are additional components not seen outside the bulge region, and that most likely can be associated with the Galactic bar. For instance, the star formation rate appears to have been slightly faster in the bulge than in the local thick disk, which most likely is an indication of the denser stellar environment closer to the Galactic centre. Our results strengthen the observational evidence that support the idea of a secular origin for the Galactic bulge, formed out of the other main Galactic stellar populations present in the central regions of our Galaxy.

The upsilon Andromedae system was the first multi-planet system discovered orbiting a main sequence star. We describe the detection of water vapor in the atmosphere of the innermost non-transiting gas giant ups~And~b by treating the star-planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet's motion and break the mass-inclination degeneracy for this non-transiting planet via deep combined flux observations of the star and the planet. In total, seven epochs of Keck NIRSPEC $L$ band observations, three epochs of Keck NIRSPEC short wavelength $K$ band observations, and three epochs of Keck NIRSPEC long wavelength $K$ band observations of the ups~And~system were obtained. We perform a multi-epoch cross correlation of the full data set with an atmospheric model. We measure the radial projection of the Keplerian velocity ($K_P$ = 55 $\pm$ 9 km/s), true mass ($M_b$ = 1.7 $^{+0.33}_{-0.24}$ $M_J$), and orbital inclination \big($i_b$ = 24 $\pm$ 4$^{\circ}$\big), and determine that the planet's opacity structure is dominated by water vapor at the probed wavelengths. Dynamical simulations of the planets in the ups~And~system with these orbital elements for ups~And~b show that stable, long-term (100 Myr) orbital configurations exist. These measurements will inform future studies of the stability and evolution of the ups~And~system, as well as the atmospheric structure and composition of the hot Jupiter.

Earth's tectonic processes regulate the formation of continental crust, control its unique deep water and carbon cycles, and are vital to its surface habitability. A major driver of steady-state plate tectonics on Earth is the sinking of the cold subducting plate into the underlying mantle. This sinking is the result of the combined effects of the thermal contraction of the lithosphere and of metamorphic transitions within the basaltic oceanic crust and lithospheric mantle. The latter of these effects is dependent on the bulk composition of the planet, e.g., the major, terrestrial planet-building elements Mg, Si, Fe, Ca, Al, and Na, which vary in abundance across the Galaxy. We present thermodynamic phase-equilibria calculations of planetary differentiation to calculate both melt composition and mantle mineralogy, and show that a planet's refractory and moderately-volatile elemental abundances control a terrestrial planet's likelihood to produce mantle-derived, melt-extracted crusts that sink. Those planets forming with a higher concentration of Si and Na abundances are less likely to undergo sustained tectonics compared to the Earth. We find only 1/3 of the range of stellar compositions observed in the Galaxy is likely to host planets able to sustain density-driven tectonics compared to the Sun/Earth. Systems outside of this compositional range are less likely to produce planets able to tectonically regulate their climate and may be inhospitable to life as we know it.

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.

The size of a planet is an observable property directly connected to the physics of its formation and evolution. We used precise radius measurements from the California-Kepler Survey (CKS) to study the size distribution of 2025 $\textit{Kepler}$ planets in fine detail. We detect a factor of $\geq$2 deficit in the occurrence rate distribution at 1.5-2.0 R$_{\oplus}$. This gap splits the population of close-in ($P$ < 100 d) small planets into two size regimes: R$_P$ < 1.5 R$_{\oplus}$ and R$_P$ = 2.0-3.0 R$_{\oplus}$, with few planets in between. Planets in these two regimes have nearly the same intrinsic frequency based on occurrence measurements that account for planet detection efficiencies. The paucity of planets between 1.5 and 2.0 R$_{\oplus}$ supports the emerging picture that close-in planets smaller than Neptune are composed of rocky cores measuring 1.5 R$_{\oplus}$ or smaller with varying amounts of low-density gas that determine their total sizes.

The California-Kepler Survey (CKS) is an observational program to improve our knowledge of the properties of stars found to host transiting planets by NASA's Kepler Mission. The improvement stems from new high-resolution optical spectra obtained using HIRES at the W. M. Keck Observatory. The CKS stellar sample comprises 1305 stars classified as Kepler Objects of Interest, hosting a total of 2075 transiting planets. The primary sample is magnitude-limited (Kp < 14.2) and contains 960 stars with 1385 planets. The sample was extended to include some fainter stars that host multiple planets, ultra short period planets, or habitable zone planets. The spectroscopic parameters were determined with two different codes, one based on template matching and the other on direct spectral synthesis using radiative transfer. We demonstrate a precision of 60 K in effective temperature, 0.10 dex in surface gravity, 0.04 dex in [Fe/H], and 1.0 km/s in projected rotational velocity. In this paper we describe the CKS project and present a uniform catalog of spectroscopic parameters. Subsequent papers in this series present catalogs of derived stellar properties such as mass, radius and age; revised planet properties; and statistical explorations of the ensemble. CKS is the largest survey to determine the properties of Kepler stars using a uniform set of high-resolution, high signal-to-noise ratio spectra. The HIRES spectra are available to the community for independent analyses.

We present stellar and planetary properties for 1305 Kepler Objects of Interest (KOIs) hosting 2025 planet candidates observed as part of the California-Kepler Survey. We combine spectroscopic constraints, presented in Paper I, with stellar interior modeling to estimate stellar masses, radii, and ages. Stellar radii are typically constrained to 11%, compared to 40% when only photometric constraints are used. Stellar masses are constrained to 4%, and ages are constrained to 30%. We verify the integrity of the stellar parameters through comparisons with asteroseismic studies and Gaia parallaxes. We also recompute planetary radii for 2025 planet candidates. Because knowledge of planetary radii is often limited by uncertainties in stellar size, we improve the uncertainties in planet radii from typically 42% to 12%. We also leverage improved knowledge of stellar effective temperature to recompute incident stellar fluxes for the planets, now precise to 21%, compared to a factor of two when derived from photometry.

HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.

The critical Lyman--Werner flux required for direct collapse blackholes (DCBH) formation, or J$_{crit}$, depends on the shape of the irradiating spectral energy distribution (SED). The SEDs employed thus far have been representative of {{realistic}} single stellar populations. We study the effect of binary stellar populations on the formation of DCBH, as a result of their contribution to the Lyman--Werner radiation field. Although binary populations with ages $>$ 10 Myr yield a larger LW photon output, we find that the corresponding values of J$_{crit}$ can be up to 100 times higher than single stellar populations. We attribute this to the shape of the binary SEDs as they produce a sub--critical rate of H$^-$ photodetaching 0.76 eV photons as compared to single stellar populations, reaffirming the role that H$^-$ plays in DCBH formation. This further corroborates the idea that DCBH formation is better understood in terms of a critical region in the H$_2$--H$^-$ photo--destruction rate parameter space, rather than a single value of LW flux.

In this study we focus on the $z\sim6.6$ Lyman-$\alpha$ CR7 consisting of clump A that is host to a potential direct collapse black hole (DCBH), and two metal enriched star forming clumps B and C. In contrast to claims that signatures of metals rule out the existence of DCBHs, we show that metal pollution of A from star forming clumps clumps B and C is inevitable, and that A can form a DCBH well before its metallicity exceeds the critical threshold of $10^{-5}-10^{-6}\ \rm Z_{\odot}$. Assuming metal mixing happens instantaneously, we derive the metallicity of A based on the star formation history of B and C. We find that treating a final accreting black hole of $10^6-10^7\ \rm M_{\odot}$ in A for nebular emission already pushes its $H_{160}$ - [3.6] and [3.6]-[4.5] colours into the 3$\sigma$ limit of observations. Hence, we show that the presence of metals in DCBH hosts is inevitable, and that it is the coevolution of the LW radiation field and metals originating from neighbouring galaxies that governs DCBH formation in a neighbouring {initially} pristine atomic cooling haloes.

We investigate the environment in which direct-collapse black holes may form by analysing a cosmological, hydrodynamical simulation that is part of the First Billion Years project. This simulation includes the most relevant physical processes leading to direct collapse of haloes, most importantly, molecular hydrogen depletion by dissociation of $H_2$ and $H^-$ from the evolving Lyman-Werner radiation field. We selected a sample of pristine atomic cooling haloes that have never formed stars in their past, have not been polluted with heavy elements and are cooling predominantly via atomic hydrogen lines. Amongst them we identified six haloes that could potentially harbour massive seed black holes formed via direct collapse (with masses in the range of $10^{4-6} M_{sun}$). These potential hosts of direct-collapse black holes form as satellites and are found within 15 physical kpc of proto-galaxies, with stellar masses in the range $10^{5-7} M_{sun}$ and maximal star formation rates of 0.1 Msun/yr over the past 5 Myr, and are exposed to the highest flux of Lyman-Werner radiation emitted from the neighbouring galaxies. It is the proximity to these proto-galaxies that differentiates these haloes from rest of the sample.

Doppler-based planet surveys point to an increasing occurrence rate of giant planets with stellar mass. Such surveys rely on evolved stars for a sample of intermediate-mass stars (so-called retired A stars), which are more amenable to Doppler observations than their main-sequence progenitors. However, it has been hypothesised that the masses of subgiant and low-luminosity red-giant stars targeted by these surveys --- typically derived from a combination of spectroscopy and isochrone fitting --- may be systematically overestimated. Here, we test this hypothesis for the particular case of the exoplanet-host star HD 212771 using K2 asteroseismology. The benchmark asteroseismic mass ($1.45^{+0.10}_{-0.09}\:\text{M}_{\odot}$) is significantly higher than the value reported in the discovery paper ($1.15\pm0.08\:\text{M}_{\odot}$), which has been used to inform the stellar mass-planet occurrence relation. This result, therefore, does not lend support to the above hypothesis. Implications for the fates of planetary systems are sensitively dependent on stellar mass. Based on the derived asteroseismic mass, we predict the post-main-sequence evolution of the Jovian planet orbiting HD 212771 under the effects of tidal forces and stellar mass loss.

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity dependent convective mixing length is used, with $ \Delta\alpha_{\rm ML, YREC} \sim 0.2$ per dex in metallicity, a trend inconsistent with the predictions of three dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of 2 even at modest deviations from solar metallicity ([Fe/H]=$-$0.5).

We report a detailed characterization of the Kepler-19 system. This star was previously known to host a transiting planet with a period of 9.29 days, a radius of 2.2 R$_\oplus$ and an upper limit on the mass of 20 M$_\oplus$. The presence of a second, non-transiting planet was inferred from the transit time variations (TTVs) of Kepler-19b, over 8 quarters of Kepler photometry, although neither mass nor period could be determined. By combining new TTVs measurements from all the Kepler quarters and 91 high-precision radial velocities obtained with the HARPS-N spectrograph, we measured through dynamical simulations a mass of $8.4 \pm 1.6$ M$_\oplus$ for Kepler-19b. From the same data, assuming system coplanarity, we determined an orbital period of 28.7 days and a mass of $13.1 \pm 2.7$ M$_\oplus$ for Kepler-19c and discovered a Neptune-like planet with a mass of $20.3 \pm 3.4$ M$_\oplus$ on a 63 days orbit. By comparing dynamical simulations with non-interacting Keplerian orbits, we concluded that neglecting interactions between planets may lead to systematic errors that could hamper the precision in the orbital parameters when the dataset spans several years. With a density of $4.32 \pm 0.87$ g cm$^{-3}$ ($0.78 \pm 0.16$ $\rho_\oplus$) Kepler-19b belongs to the group of planets with a rocky core and a significant fraction of volatiles, in opposition to low-density planets characterized by transit-time variations only and the increasing number of rocky planets with Earth-like density. Kepler-19 joins the small number of systems that reconcile transit timing variation and radial velocity measurements.

We assess the impact of trapped Lyman alpha cooling radiation on the formation of direct collapse black holes (DCBHs). We apply a one-zone chemical and thermal evolution model, accounting for the photodetachment of H- ions, precursors to the key coolant H2, by Lyman alpha photons produced during the collapse of a cloud of primordial gas in an atomic cooling halo at high redshift. We find that photodetachment of H- by trapped Lyman alpha photons may lower the level of the H2-dissociating background radiation field required for DCBH formation substantially, dropping the critical flux by up to a factor of a few. This translates into a potentially large increase in the expected number density of DCBHs in the early Universe, and supports the view that DCBHs may be the seeds for the BHs residing in the centers of a significant fraction of galaxies today. We find that detachment of H- by Lyman alpha has the strongest impact on the critical flux for the relatively high background radiation temperatures expected to characterize the emission from young, hot stars in the early Universe. This lends support to the DCBH origin of the highest redshift quasars.

The prospects for the habitability of M-dwarf planets have long been debated, due to key differences between the unique stellar and planetary environments around these low-mass stars, as compared to hotter, more luminous Sun-like stars. Over the past decade, significant progress has been made by both space- and ground-based observatories to measure the likelihood of small planets to orbit in the habitable zones of M-dwarf stars. We now know that most M dwarfs are hosts to closely-packed planetary systems characterized by a paucity of Jupiter-mass planets and the presence of multiple rocky planets, with roughly a third of these rocky M-dwarf planets orbiting within the habitable zone, where they have the potential to support liquid water on their surfaces. Theoretical studies have also quantified the effect on climate and habitability of the interaction between the spectral energy distribution of M-dwarf stars and the atmospheres and surfaces of their planets. These and other recent results fill in knowledge gaps that existed at the time of the previous overview papers published nearly a decade ago by Tarter et al. (2007) and Scalo et al. (2007). In this review we provide a comprehensive picture of the current knowledge of M-dwarf planet occurrence and habitability based on work done in this area over the past decade, and summarize future directions planned in this quickly evolving field.

We report the first detailed chemical abundance analysis of the exoplanet-hosting M-dwarf stars Kepler-138 and Kepler-186 from the analysis of high-resolution ($R$ $\sim$ 22,500) $H$-band spectra from the SDSS IV - APOGEE survey. Chemical abundances of thirteen elements - C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, and Fe - are extracted from the APOGEE spectra of these early M-dwarfs via spectrum syntheses computed with an improved line list that takes into account H$_{2}$O and FeH lines. This paper demonstrates that APOGEE spectra can be analyzed to determine detailed chemical compositions of M-dwarfs. Both exoplanet-hosting M-dwarfs display modest sub-solar metallicities: [Fe/H]$_{Kepler-138}$ = -0.09 $\pm$ 0.09 dex and [Fe/H]$_{Kepler-186}$ = -0.08 $\pm$ 0.10 dex. The measured metallicities resulting from this high-resolution analysis are found to be higher by $\sim$0.1-0.2 dex than previous estimates from lower-resolution spectra. The C/O ratios obtained for the two planet-hosting stars are near-solar, with values of 0.55 $\pm$ 0.10 for Kepler-138 and 0.52 $\pm$ 0.12 for Kepler-186. Kepler-186 exhibits a marginally enhanced [Si/Fe] ratio.

We present the discovery of two extended $\sim$0.12 mag dimming events of the weak-lined T-Tauri star V1334. The start of the first event was missed but came to an end in late 2003, and the second began in February 2009, and continues as of November 2016. Since the egress of the current event has not yet been observed, it suggests a period of $>$13 years if this event is periodic. Spectroscopic observations suggest the presence of a small inner disk, although the spectral energy distribution shows no infrared excess. We explore the possibility that the dimming events are caused by an orbiting body (e.g. a disk warp or dust trap), enhanced disk winds, hydrodynamical fluctuations of the inner disk, or a significant increase in the magnetic field flux at the surface of the star. We also find a $\sim$0.32 day periodic photometric signal that persists throughout the 2009 dimming which appears to not be due to ellipsoidal variations from a close stellar companion. High precision photometric observations of V1334 Tau during K2 campaign 13, combined with simultaneous photometric and spectroscopic observations from the ground, will provide crucial information about the photometric variability and its origin.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Finding the first cosmic explosions in the universe will only be possible with the James Webb Space Telescope, the Wide-Field Infrared Survey Telescope and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Cluster Lensing and Supernova Survey with Hubble (CLASH), could reveal such events now by magnifying their flux by factors of 10 or more. We find that CLASH will likely discover at least 2 - 3 core-collapse supernovae at 5 < z < 12 and perhaps as many as ten. Future surveys of cluster lenses similar in scope to CLASH by the James Webb Space Telescope might find hundreds of these events out to z ~ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae will also constrain cosmic star formation rates in the era of first galaxy formation.

Chemical evolution models are powerful tools for interpreting stellar abundance surveys and understanding galaxy evolution. However, their predictions depend heavily on the treatment of inflow, outflow, star formation efficiency (SFE), the stellar initial mass function, the Type Ia supernova delay time distribution, stellar yields, and stellar population mixing. Using flexCE, a flexible one-zone chemical evolution code, we investigate the effects of and trade-offs between parameters. Two critical parameters are SFE and the outflow mass-loading parameter, which shift the knee in [O/Fe]-[Fe/H] and the equilibrium abundances that the simulations asymptotically approach, respectively. One-zone models with simple star formation histories follow narrow tracks in [O/Fe]-[Fe/H] unlike the observed bimodality (separate high-alpha and low-alpha sequences) in this plane. A mix of one-zone models with inflow timescale and outflow mass-loading parameter variations, motivated by the inside-out galaxy formation scenario with radial mixing, reproduces the two sequences better than a one-zone model with two infall epochs. We present [X/Fe]-[Fe/H] tracks for 20 elements assuming three different supernova yield models and find some significant discrepancies with solar neighborhood observations, especially for elements with strongly metallicity-dependent yields. We apply principal component abundance analysis (PCAA) to the simulations and existing data to reveal the main correlations amongst abundances and quantify their contributions to variation in abundance space. For the stellar population mixing scenario, the abundances of alpha-elements and elements with metallicity-dependent yields dominate the first and second principal components, respectively, and collectively explain 99% of the variance in the model. flexCE is a python package available at https://github.com/bretthandrews/flexCE.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

The occurrence rate of hot Jupiters from the Kepler transit survey is roughly half that of radial velocity surveys targeting solar neighborhood stars. One hypothesis to explain this difference is that the two surveys target stars with different stellar metallicity distributions. To test this hypothesis, we measure the metallicity distribution of the Kepler targets using the Hectochelle multi-fiber, high-resolution spectrograph. Limiting our spectroscopic analysis to 610 dwarf stars in our sample with log(g)>3.5, we measure a metallicity distribution characterized by a mean of [M/H]_{mean} = -0.045 +/- 0.00, in agreement with previous studies of the Kepler field target stars. In comparison, the metallicity distribution of the California Planet Search radial velocity sample has a mean of [M/H]_{CPS, mean} = -0.005 +/- 0.006, and the samples come from different parent populations according to a Kolmogorov-Smirnov test. We refit the exponential relation between the fraction of stars hosting a close-in giant planet and the host star metallicity using a sample of dwarf stars from the California Planet Search with updated metallicities. The best-fit relation tells us that the difference in metallicity between the two samples is insufficient to explain the discrepant Hot Jupiter occurrence rates; the metallicity difference would need to be $\simeq$0.2-0.3 dex for perfect agreement. We also show that (sub)giant contamination in the Kepler sample cannot reconcile the two occurrence calculations. We conclude that other factors, such as binary contamination and imperfect stellar properties, must also be at play.

We report the detection of two new planets orbiting the K giants HD 86950 and HD 222076, based on precise radial velocities obtained with three instruments: AAT/UCLES, FEROS, and CHIRON. HD 86950b has a period of 1270$\pm$57 days at $a=2.72\pm$0.08 AU, and m sin $i=3.6\pm$0.7 Mjup. HD 222076b has $P=871\pm$19 days at $a=1.83\pm$0.03 AU, and m sin $i=1.56\pm$0.11 Mjup. These two giant planets are typical of the population of planets known to orbit evolved stars. In addition, we find a high-amplitude periodic velocity signal ($K\sim$50 m/s) in HD 29399, and show that it is due to stellar variability rather than Keplerian reflex motion. We also investigate the relation between planet occurrence and host-star metallicity for the 164-star Pan-Pacific Planet Search sample of evolved stars. In spite of the small sample of PPPS detections, we confirm the trend of increasing planet occurrence as a function of metallicity found by other studies of planets orbiting evolved stars.

We report chemical abundances obtained by SDSS-III/APOGEE for giant stars in five globular clusters located within 2.2 kpc of the Galactic centre. We detect the presence of multiple stellar populations in four of those clusters (NGC 6553, NGC 6528, Terzan 5, and Palomar 6) and find strong evidence for their presence in NGC 6522. All clusters present a significant spread in the abundances of N, C, Na, and Al, with the usual correlations and anti-correlations between various abundances seen in other globular clusters. Our results provide important quantitative constraints on theoretical models for self-enrichment of globular clusters, by testing their predictions for the dependence of yields of elements such as Na, N, C, and Al on metallicity. They also confirm that, under the assumption that field N-rich stars originate from globular cluster destruction, they can be used as tracers of their parental systems in the high- metallicity regime.

We report on detection, from observations obtained with the APOGEE spectroscopic survey, of a metal-poor ([Fe/H] $= -1.3$ dex) field giant star with an extreme Mg-Al abundance ratio ([Mg/Fe] $= -0.31$ dex; [Al/Fe] $= 1.49$ dex). Such low Mg/Al ratios are seen only among the second-generation population of globular clusters, and are not present among Galactic disk field stars. The light element abundances of this star, 2M16011638-1201525, suggest that it could have been born in a globular cluster. We explore several origin scenarios, in particular studying the orbit of the star to check the probability of it being kinematically related to known globular clusters. We performed simple orbital integrations assuming the estimated distance of 2M16011638-1201525 and the available six-dimensional phase-space coordinates of 63 globular clusters, looking for close encounters in the past with a minimum distance approach within the tidal radius of each cluster. We found a very low probability that 2M16011638-1201525 was ejected from most globular clusters; however, we note that the best progenitor candidate to host this star is globular cluster $\omega$ Centauri (NGC 5139). Our dynamical investigation demonstrates that 2M16011638-1201525 reaches a distance $|Z_{max}| < 3 $ kpc from the Galactic plane and a minimum and maximum approach to the Galactic center of $R_{min}<0.62$ kpc and $R_{max}<7.26$ kpc in an eccentric ($e\sim0.53$) and retrograde orbit. Since the extreme chemical anomaly of 2M16011638-1201525 has also been observed in halo field stars, this object could also be considered a halo contaminant, likely been ejected into the Milky Way disk from the halo. We conclude that, 2M16011638-20152 is also kinematically consistent with the disk but chemically consistent with halo field stars.

We target the thermal emission spectrum of the non-transiting gas giant HD 88133 b with high-resolution near-infrared spectroscopy, by treating the planet and its host star as a spectroscopic binary. For sufficiently deep summed flux observations of the star and planet across multiple epochs, it is possible to resolve the signal of the hot gas giant's atmosphere compared to the brighter stellar spectrum, at a level consistent with the aggregate shot noise of the full data set. To do this, we first perform a principal component analysis to remove the contribution of the Earth's atmosphere to the observed spectra. Then, we use a cross-correlation analysis to tease out the spectra of the host star and HD 88133 b to determine its orbit and identify key sources of atmospheric opacity. In total, six epochs of Keck NIRSPEC L band observations and three epochs of Keck NIRSPEC K band observations of the HD 88133 system were obtained. Based on an analysis of the maximum likelihood curves calculated from the multi-epoch cross correlation of the full data set with two atmospheric models, we report the direct detection of the emission spectrum of the non-transiting exoplanet HD 88133 b and measure a radial projection of the Keplerian orbital velocity of 40 $\pm$ 15 km/s, a true mass of 1.02$^{+0.61}_{-0.28}M_J$, a nearly face-on orbital inclination of 15${^{+6}_{-5}}^{\circ}$, and an atmosphere opacity structure at high dispersion dominated by water vapor. This, combined with eleven years of radial velocity measurements of the system, provides the most up-to-date ephemeris for HD 88133.

HD 179070, aka Kepler-21, is a V = 8.25 F6IV star and the brightest exoplanet host discovered by Kepler. An early detailed analysis by Howell et al. (2012) of the first thirteen months (Q0 - Q5) of Kepler light curves revealed transits of a planetary companion, Kepler-21b, with a radius of about 1.60 +/- 0.04 R_earth and an orbital period of about 2.7857 days. However, they could not determine the mass of the planet from the initial radial velocity observations with Keck-HIRES, and were only able to impose a 2-sigma upper limit of 10 M_earth. Here we present results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points. We detect the Doppler signal of Kepler-21b with a radial velocity semi-amplitude K = 2.00 +/- 0.65 m/s, which corresponds to a planetary mass of 5.1 +/- 1.7 M_earth. We also measure an improved radius for the planet of 1.639 (+0.019, -0.015) R_earth, in agreement with the radius reported by Howell et al. (2012). We conclude that Kepler-21b, with a density of 6.4 +/- 2.1 g/cm^3, belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape. The radial velocity analysis presented in this paper serves as example of the type of analysis that will be necessary to confirm the masses of TESS small planet candidates.

Much of the inner Milky Way's (MW) global rotation and velocity dispersion patterns can be reproduced by models of secularly-evolved, bar-dominated bulges. More sophisticated constraints, including the higher moments of the line-of-sight velocity distributions (LOSVDs) and limits on the chemodynamical substructure, are critical for interpreting observations of the unresolved inner regions of extragalactic systems and for placing the MW in context with other galaxies. Here, we use SDSS-APOGEE data to develop these constraints, by presenting the first maps of the LOSVD skewness and kurtosis of metal-rich and metal-poor inner MW stars (divided at [Fe/H] = -0.4), and comparing the observed patterns to those that are seen both in N-body models and in extragalactic bars. Despite closely matching the mean velocity and dispersion, the models do not reproduce the observed LOSVD skewness patterns in different ways, which demonstrates that our understanding of the detailed orbital structure of the inner MW remains an important regime for improvement. We find evidence in the MW of the skewness-velocity correlation that is used as a diagnostic of extragalactic bar/bulges. This correlation appears in metal-rich stars only, providing further evidence for different evolutionary histories of chemically differentiated populations. We connect these skewness measurements to previous work on high-velocity "peaks" in the inner Galaxy, confirming the presence of that phenomenon, and we quantify the cylindrical rotation of the inner Galaxy, finding that the latitude-independent rotation vanishes outside of lon ~ 7 deg. Finally, we evaluate the MW data in light of select extragalactic bar diagnostics and discuss progress and challenges of using the MW as a resolved analog of unresolved stellar populations.

While the vast majority of multiple-planet systems have their orbital angular momentum axes aligned with the spin axis of their host star, Kepler-56 is an exception: its two transiting planets are coplanar yet misaligned by at least 40 degrees with respect to their host star. Additional follow-up observations of Kepler-56 suggest the presence of a massive, non-transiting companion that may help explain this misalignment. We model the transit data along with Keck/HIRES and HARPS-N radial velocity data to update the masses of the two transiting planets and infer the physical properties of the third, non-transiting planet. We employ a Markov Chain Monte Carlo sampler to calculate the best-fitting orbital parameters and their uncertainties for each planet. We find the outer planet has a period of 1002 $\pm$ 5 days and minimum mass of 5.61 $\pm$ 0.38 Jupiter masses. We also place a 95% upper limit of 0.80 m/s/yr on long-term trends caused by additional, more distant companions.

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $\Lambda$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

Kepler-20 is a solar-type star (V = 12.5) hosting a compact system of five transiting planets, all packed within the orbital distance of Mercury in our own Solar System. A transition from rocky to gaseous planets with a planetary transition radius of ~1.6 REarth has recently been proposed by several publications in the literature (Rogers 2015;Weiss & Marcy 2014). Kepler-20b (Rp ~ 1.9 REarth) has a size beyond this transition radius, however previous mass measurements were not sufficiently precise to allow definite conclusions to be drawn regarding its composition. We present new mass measurements of three of the planets in the Kepler-20 system facilitated by 104 radial velocity measurements from the HARPS-N spectrograph and 30 archival Keck/HIRES observations, as well as an updated photometric analysis of the Kepler data and an asteroseismic analysis of the host star (MStar = 0.948+-0.051 Msun and Rstar = 0.964+-0.018 Rsun). Kepler-20b is a 1.868+0.066-0.034 REarth planet in a 3.7 day period with a mass of 9.70+1.41-1.44 MEarth resulting in a mean density of 8.2+1.5-1.3 g/cc indicating a rocky composition with an iron to silicate ratio consistent with that of the Earth. This makes Kepler-20b the most massive planet with a rocky composition found to date. Furthermore, we report the discovery of an additional non-transiting planet with a minimum mass of 19.96+3.08-3.61 MEarth and an orbital period of ~34 days in the gap between Kepler-20f (P ~ 11 days) and Kepler-20d (P ~ 78 days).

With the advent of the space missions CoRoT and Kepler, it has become feasible to determine precise asteroseismic masses and ages for large samples of red-giant stars. In this paper, we present the CoRoGEE dataset -- obtained from CoRoT lightcurves for 606 red giant stars in two fields of the Galactic disc which have been co-observed for an ancillary project of APOGEE. We have used the Bayesian parameter estimation code PARAM to calculate distances, extinctions, masses, and ages for these stars in a homogeneous analysis, resulting in relative statistical uncertainties of $\sim2\%$ in distance, $\sim4\%$ in radius, $\sim9\%$ in mass and $\sim25\%$ in age. We also assess systematic age uncertainties due to different input physics and mass loss. We discuss the correlation between ages and chemical abundance patterns of field stars over a large radial range of the Milky Way's disc (5 kpc $<R_{\rm Gal}<$ 14 kpc), focussing on the [$\alpha$/Fe]-[Fe/H]-age plane in five radial bins of the Galactic disc. We find an overall agreement with the expectations of chemical-evolution models computed before the present data were available, especially for the outer regions. However, our data also indicate that a significant fraction of stars now observed near and beyond the Solar Neighbourhood migrated from inner regions. Mock CoRoGEE observations of a chemo-dynamical Milky Way disc model show that the number of high-metallicity stars in the outer disc is too high to be accounted for even by the strong radial mixing present in the model. The mock observations also reveal that the age distribution of the [$\alpha$/Fe]-enhanced sequence in the CoRoGEE inner-disc field is much broader than expected from a combination of radial mixing and observational errors. We suggest that a thick disc/bulge component that formed stars for more than 3 Gyr may account for these discrepancies.

The mass and age of substellar objects are degenerate parameters leaving the evolutionary state of brown dwarfs ambiguous without additional information. Theoretical models are normally used to help distinguish between old, massive brown dwarfs and young, low mass brown dwarfs but these models have yet to be properly calibrated. We have carried out an infrared high-contrast imaging program with the goal of detecting substellar objects as companions to nearby stars to help break degeneracies in inferred physical properties such as mass, age, and composition. Rather than using imaging observations alone, our targets are pre-selected based on the existence of dynamical accelerations informed from years of stellar radial velocity (RV) measurements. In this paper, we present the discovery of a rare benchmark brown dwarf orbiting the nearby ($d=18.69\pm0.19$ pc), solar-type (G9V) star HD 4747 ([Fe/H]=$-0.22\pm0.04$) with a projected separation of only $\rho=11.3\pm0.2$ AU ($\theta \approx$ 0.6"). Precise Doppler measurements taken over 18 years reveal the companion's orbit and allow us to place strong constraints on its mass using dynamics ($m \sin(i) = 55.3\pm1.9M_J$). Relative photometry ($\Delta K_s=9.05\pm0.14$, $M_{K_s}=13.00\pm0.14$, $K_s - L' = 1.34\pm0.46$) indicates that HD 4747 B is most-likely a late-type L-dwarf and, if near the L/T transition, an intriguing source for studying cloud physics, variability, and polarization. We estimate a model-dependent mass of $m=72^{+3}_{-13}M_J$ for an age of $3.3^{+2.3}_{-1.9}$ Gyr based on gyrochronology. Combining astrometric measurements with RV data, we calculate the companion dynamical mass ($m=60.2\pm3.3M_J$) and orbit ($e=0.740\pm0.002$) directly. As a new mass, age, and metallicity benchmark, HD 4747 B will serve as a laboratory for precision astrophysics to test theoretical models that describe the emergent radiation of brown dwarfs.

We present the first ab initio cosmological simulations of a CR7-like object which approximately reproduce the observed line widths and strengths. In our model, CR7 is powered by a massive ($3.23 \times 10^7$ $M_\odot$) black hole (BH) the accretion rate of which varies between $\simeq$ 0.25 and $\simeq$ 0.9 times the Eddington rate on timescales as short as 10$^3$ yr. Our model takes into account multi-dimensional effects, X-ray feedback, secondary ionizations and primordial chemistry. We estimate Ly-$\alpha$ line widths by post-processing simulation output with Monte Carlo radiative transfer and calculate emissivity contributions from radiative recombination and collisional excitation. We find the luminosities in the Lyman-$\alpha$ and He II 1640 angstrom lines to be $5.0\times10^{44}$ and $2.4\times10^{43}$ erg s$^{-1}$, respectively, in agreement with the observed values of $>$ $8.3\times10^{43}$ and $2.0\times10^{43}$ erg s$^{-1}$. We also find that the black hole heats the halo and renders it unable to produce stars as required to keep the halo metal free. These results demonstrate the viability of the BH hypothesis for CR7 in a cosmological context. Assuming the BH mass and accretion rate that we find, we estimate the synchrotron luminosity of CR7 to be $P \simeq 10^{40} - 10^{41}$ erg s$^{-1}$, which is sufficiently luminous to be observed in $\mu$Jy observations and would discriminate this scenario from one where the luminosity is driven by Population III stars.

We report the discovery of a transiting exoplanet, KELT-11b, orbiting the bright ($V=8.0$) subgiant HD 93396. A global analysis of the system shows that the host star is an evolved subgiant star with $T_{\rm eff} = 5370\pm51$ K, $M_{*} = 1.438_{-0.052}^{+0.061} M_{\odot}$, $R_{*} = 2.72_{-0.17}^{+0.21} R_{\odot}$, log $g_*= 3.727_{-0.046}^{+0.040}$, and [Fe/H]$ = 0.180\pm0.075$. The planet is a low-mass gas giant in a $P = 4.736529\pm0.00006$ day orbit, with $M_{P} = 0.195\pm0.018 M_J$, $R_{P}= 1.37_{-0.12}^{+0.15} R_J$, $\rho_{P} = 0.093_{-0.024}^{+0.028}$ g cm$^{-3}$, surface gravity log ${g_{P}} = 2.407_{-0.086}^{+0.080}$, and equilibrium temperature $T_{eq} = 1712_{-46}^{+51}$ K. KELT-11 is the brightest known transiting exoplanet host in the southern hemisphere by more than a magnitude, and is the 6th brightest transit host to date. The planet is one of the most inflated planets known, with an exceptionally large atmospheric scale height (2763 km), and an associated size of the expected atmospheric transmission signal of 5.6%. These attributes make the KELT-11 system a valuable target for follow-up and atmospheric characterization, and it promises to become one of the benchmark systems for the study of inflated exoplanets.

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4\sigma$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7\sigma$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

Formation of globular clusters (GCs), the Galactic bulge, or galaxy bulges in general, are important unsolved problems in Galactic astronomy. Homogeneous infrared observations of large samples of stars belonging to GCs and the Galactic bulge field are one of the best ways to study these problems. We report the discovery by APOGEE of a population of field stars in the inner Galaxy with abundances of N, C, and Al that are typically found in GC stars. The newly discovered stars have high [N/Fe], which is correlated with [Al/Fe] and anti-correlated with [C/Fe]. They are homogeneously distributed across, and kinematically indistinguishable from, other field stars in the same volume. Their metallicity distribution is seemingly unimodal, peaking at [Fe/H]~-1, thus being in disagreement with that of the Galactic GC system. Our results can be understood in terms of different scenarios. N-rich stars could be former members of dissolved GCs, in which case the mass in destroyed GCs exceeds that of the surviving GC system by a factor of ~8. In that scenario, the total mass contained in so-called "first-generation" stars cannot be larger than that in "second-generation" stars by more than a factor of ~9 and was certainly smaller. Conversely, our results may imply the absence of a mandatory genetic link between "second generation" stars and GCs. Last, but not least, N-rich stars could be the oldest stars in the Galaxy, the by-products of chemical enrichment by the first stellar generations formed in the heart of the Galaxy.

We report on the discovery and characterization of the transiting planet K2-39b (EPIC 206247743b). With an orbital period of 4.6 days, it is the shortest-period planet orbiting a subgiant star known to date. Such planets are rare, with only a handful of known cases. The reason for this is poorly understood, but may reflect differences in planet occurrence around the relatively high-mass stars that have been surveyed, or may be the result of tidal destruction of such planets. K2-39 is an evolved star with a spectroscopically derived stellar radius and mass of $3.88^{+0.48}_{-0.42}~\mathrm{R_\odot}$ and $1.53^{+0.13}_{-0.12}~\mathrm{M_\odot}$, respectively, and a very close-in transiting planet, with $a/R_\star = 3.4$. Radial velocity (RV) follow-up using the HARPS, FIES and PFS instruments leads to a planetary mass of $50.3^{+9.7}_{-9.4}~\mathrm{M_\oplus}$. In combination with a radius measurement of $8.3 \pm 1.1~\mathrm{R_\oplus}$, this results in a mean planetary density of $0.50^{+0.29}_{-0.17}$ g~cm$^{-3}$. We furthermore discover a long-term RV trend, which may be caused by a long-period planet or stellar companion. Because K2-39b has a short orbital period, its existence makes it seem unlikely that tidal destruction is wholly responsible for the differences in planet populations around subgiant and main-sequence stars. Future monitoring of the transits of this system may enable the detection of period decay and constrain the tidal dissipation rates of subgiant stars.

We present the effective temperatures ($T_{\rm eff}$), metallicities, and colours in SDSS, 2MASS, and WISE filters, of a sample of 3834 late-K and early-M dwarfs selected from the Sloan Digital Sky Survey APOGEE spectroscopic survey ASPCAP catalog. We confirm that ASPCAP $T_{\rm eff}$ values between 3550 K$<T_{\rm eff}<$4200 K are accurate to $\sim$100 K compared to interferometric $T_{\rm eff}$ values. In that same $T_{\rm eff}$ range, ASPCAP metallicities are accurate to 0.18 dex between $-1.0<$[M/H]$<0.2$. For these cool dwarfs, nearly every colour is sensitive to both $T_{\rm eff}$ and metallicity. Notably, we find that $g-r$ is not a good indicator of metallicity for near-solar metallicity early-M dwarfs. We confirm that $J-K_S$ colour is strongly dependent on metallicity, and find that $W1-W2$ colour is a promising metallicity indicator. Comparison of the late-K and early-M dwarf colours, metallicities, and $T_{\rm eff}$ to those from three different model grids shows reasonable agreement in $r-z$ and $J-K_S$ colours, but poor agreement in $u-g$, $g-r$, and $W1-W2$. Comparison of the metallicities of the KM dwarf sample to those from previous colour-metallicity relations reveals a lack of consensus in photometric metallicity indicators for late-K and early-M dwarfs. We also present empirical relations for $T_{\rm eff}$ as a function of $r-z$ colour combined with either [M/H] or $W1-W2$ colour, and for [M/H] as a function of $r-z$ and $W1-W2$ colour. These relations yield $T_{\rm eff}$ to $\sim$100 K and [M/H] to $\sim$0.18 dex precision with colours alone, for $T_{\rm eff}$ in the range of 3550-4200 K and [M/H] in the range of $-$0.5-0.2.

In an effort to measure the masses of planets discovered by the NASA {\it K2} mission, we have conducted precise Doppler observations of five stars with transiting planets. We present the results of a joint analysis of these new data and previously published Doppler data. The first star, an M dwarf known as K2-3 or EPIC~201367065, has three transiting planets ("b", with radius $2.1~R_{\oplus}$; "c", $1.7~R_{\oplus}$; and "d", $1.5~R_{\oplus}$). Our analysis leads to the mass constraints: $M_{b}=8.1^{+2.0}_{-1.9}~M_{\oplus}$ and $M_{c}$ < $ 4.2~M_{\oplus}$~(95\%~conf.). The mass of planet d is poorly constrained because its orbital period is close to the stellar rotation period, making it difficult to disentangle the planetary signal from spurious Doppler shifts due to stellar activity. The second star, a G dwarf known as K2-19 or EPIC~201505350, has two planets ("b", $7.7~R_{\oplus}$; and "c", $4.9~R_{\oplus}$) in a 3:2 mean-motion resonance, as well as a shorter-period planet ("d", $1.1~R_{\oplus}$). We find $M_{b}$= $28.5^{+5.4}_{-5.0} ~M_{\oplus}$, $M_{c}$= $25.6^{+7.1}_{-7.1} ~M_{\oplus}$ and $M_{d}$ < $14.0~M_{\oplus} $~(95\%~conf.). The third star, a G dwarf known as K2-24 or EPIC~203771098, hosts two transiting planets ("b", $5.7~R_{\oplus}$; and "c", $7.8~R_{\oplus}$) with orbital periods in a nearly 2:1 ratio. We find $M_{b}$= $19.8^{+4.5}_{-4.4} ~M_{\oplus}$ and $M_{c}$ = $26.0^{+5.8}_{-6.1}~M_{\oplus}$.....

We present spectroscopic stellar parameters for the complete target list of 164 evolved stars from the Pan-Pacific Planet Search, a five-year radial velocity campaign using the 3.9m Anglo-Australian Telescope. For 87 of these bright giants, our work represents the first determination of their fundamental parameters. Our results carry typical uncertainties of 100 K, 0.15 dex, and 0.1 dex in $T_{\rm eff}$, $\log g$, and [Fe/H] and are consistent with literature values where available. The derived stellar masses have a mean of $1.31^{+0.28}_{-0.25}$ Msun, with a tail extending to $\sim$2 Msun, consistent with the interpretation of these targets as "retired" A-F type stars.

Future generations of precise radial velocity (RV) surveys aim to achieve sensitivity sufficient to detect Earth mass planets orbiting in their stars' habitable zones. A major obstacle to this goal is astrophysical radial velocity noise caused by active areas moving across the stellar limb as a star rotates. In this paper, we quantify how stellar activity impacts exoplanet detection with radial velocities as a function of orbital and stellar rotational periods. We perform data-driven simulations of how stellar rotation affects planet detectability and compile and present relations for the typical timescale and amplitude of stellar radial velocity noise as a function of stellar mass. We show that the characteristic timescales of quasi-periodic radial velocity jitter from stellar rotational modulations coincides with the orbital period of habitable zone exoplanets around early M-dwarfs. These coincident periods underscore the importance of monitoring the targets of RV habitable zone planet surveys through simultaneous photometric measurements for determining rotation periods and activity signals, and mitigating activity signals using spectroscopic indicators and/or RV measurements at different wavelengths.

There are no field brown dwarf analogs with measured masses, radii, and luminosities, precluding our ability to connect the population of transiting brown dwarfs with measurable masses and radii and field brown dwarfs with measurable luminosities and atmospheric properties. LHS 6343 C, a weakly-irradiated brown dwarf transiting one member of an M+M binary in the Kepler field, provides the first opportunity to probe the atmosphere of a non-inflated brown dwarf with a measured mass and radius. Here, we analyze four Spitzer observations of secondary eclipses of LHS 6343 C behind LHS 6343 A. Jointly fitting the eclipses with a Gaussian process noise model of the instrumental systematics, we measure eclipse depths of 1.06 \pm 0.21 ppt at 3.6 microns and 2.09 \pm 0.08 ppt at 4.5 microns, corresponding to brightness temperatures of 1026 \pm 57 K and 1249 \pm 36 K, respectively. We then apply brown dwarf evolutionary models to infer a bolometric luminosity log(L_star / L_sun) = -5.16 \pm 0.04. Given the known physical properties of the brown dwarf and the two M dwarfs in the LHS 6343 C system, these depths are consistent with models of a 1100 K T dwarf at an age of 5 Gyr and empirical observations of field T5-6 dwarfs with temperatures of 1070 \pm 130 K. We investigate the possibility that the orbit of LHS 6343 C has been altered by the Kozai-Lidov mechanism and propose additional astrometric or Rossiter-McLaughlin measurements of the system to probe the dynamical history of the system.

Given that low-mass stars have intrinsically low luminosities at optical wavelengths and a propensity for stellar activity, it is advantageous for radial velocity (RV) surveys of these objects to use near-infrared (NIR) wavelengths. In this work we describe and test a novel RV extraction pipeline dedicated to retrieving RVs from low mass stars using NIR spectra taken by the CSHELL spectrograph at the NASA Infrared Telescope Facility, where a methane isotopologue gas cell is used for wavelength calibration. The pipeline minimizes the residuals between the observations and a spectral model composed of templates for the target star, the gas cell, and atmospheric telluric absorption; models of the line spread function, continuum curvature, and sinusoidal fringing; and a parameterization of the wavelength solution. The stellar template is derived iteratively from the science observations themselves without a need for separate observations dedicated to retrieving it. Despite limitations from CSHELL's narrow wavelength range and instrumental systematics, we are able to (1) obtain an RV precision of 35 m/s for the RV standard star GJ 15 A over a time baseline of 817 days, reaching the photon noise limit for our attained SNR, (2) achieve ~3 m/s RV precision for the M giant SV Peg over a baseline of several days and confirm its long-term RV trend due to stellar pulsations, as well as obtain nightly noise floors of ~2 - 6 m/s, and (3) show that our data are consistent with the known masses, periods, and orbital eccentricities of the two most massive planets orbiting GJ 876. Future applications of our pipeline to RV surveys using the next generation of NIR spectrographs, such as iSHELL, will enable the potential detection of Super-Earths and Mini-Neptunes in the habitable zones of M dwarfs.

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

We present the results of a precise near-infrared (NIR) radial velocity (RV) survey of 32 low-mass stars with spectral types K2-M4 using CSHELL at the NASA IRTF in the $K$-band with an isotopologue methane gas cell to achieve wavelength calibration and a novel iterative RV extraction method. We surveyed 14 members of young ($\approx$ 25-150 Myr) moving groups, the young field star $\varepsilon$ Eridani as well as 18 nearby ($<$ 25 pc) low-mass stars and achieved typical single-measurement precisions of 8-15 m s$^{-1}$ with a long-term stability of 15-50 m s$^{-1}$. We obtain the best NIR RV constraints to date on 27 targets in our sample, 19 of which were never followed by high-precision RV surveys. Our results indicate that very active stars can display long-term RV variations as low as $\sim$ 25-50 m s$^{-1}$ at $\approx$ 2.3125 $\mu$m, thus constraining the effect of jitter at these wavelengths. We provide the first multi-wavelength confirmation of GJ 876 bc and independently retrieve orbital parameters consistent with previous studies. We recovered RV variability for HD 160934 AB and GJ 725 AB that are consistent with their known binary orbits, and nine other targets are candidate RV variables with a statistical significance of 3-5$\sigma$. Our method combined with the new iSHELL spectrograph will yield long-term RV precisions of $\lesssim$ 5 m s$^{-1}$ in the NIR, which will allow the detection of Super-Earths near the habitable zone of mid-M dwarfs.

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

With detections of quasars powered by increasingly massive black holes (BHs) at increasingly early times in cosmic history over the past decade, there has been correspondingly rapid progress made on the theory of early BH formation and growth. Here we review the emerging picture of how the first massive BHs formed from the primordial gas and then grew to supermassive scales. We discuss the initial conditions for the formation of the progenitors of these seed BHs, the factors dictating the initial masses with which they form, and their initial stages of growth via accretion, which may occur at super-Eddington rates. Finally, we briefly discuss how these results connect to large-scale simulations of the growth of supermassive BHs over the course of the first billion years following the Big Bang.

We conducted a Doppler survey at Keck combined with NIRC2 K-band AO imaging to search for massive, long-period companions to 123 known exoplanet systems with one or two planets detected using the radial velocity (RV) method. Our survey is sensitive to Jupiter mass planets out to 20 AU for a majority of stars in our sample, and we report the discovery of eight new long-period planets, in addition to 20 systems with statistically significant RV trends indicating the presence of an outer companion beyond 5 AU. We combine our RV observations with AO imaging to determine the range of allowed masses and orbital separations for these companions, and account for variations in our sensitivity to companions among stars in our sample. We estimate the total occurrence rate of companions in our sample to be 52 +/- 5% over the range 1 - 20 M_Jup and 5 - 20 AU. Our data also suggest a declining frequency for gas giant planets in these systems beyond 3-10 AU, in contrast to earlier studies that found a rising frequency for giant planets in the range 0.01-3 AU. This suggests either that the frequency of gas giant planets peaks between 3-10 AU, or that outer companions in these systems have a different semi-major axis distribution than the overall gas giant planet population. Our results also suggest that hot gas giants may be more likely to have an outer companion than cold gas giants. We find that planets with an outer companion have higher average eccentricities than their single counterparts, suggesting that dynamical interactions between planets may play an important role in these systems.

Metallicity gradients provide strong constraints for understanding the chemical evolution of the Galaxy. We report on radial abundance gradients of Fe, Ni, Ca, Si, and Mg obtained from a sample of 304 red-giant members of 29 disk open clusters, mostly concentrated at galactocentric distances between ~8 - 15 kpc, but including two open clusters in the outer disk. The observations are from the APOGEE survey. The chemical abundances were derived automatically by the ASPCAP pipeline and these are part of the SDSS III Data Release 12. The gradients, obtained from least squares fits to the data, are relatively flat, with slopes ranging from -0.026 to -0.033 dex/kpc for the alpha-elements [O/H], [Ca/H], [Si/H] and [Mg/H] and -0.035 dex/kpc and -0.040 dex/kpc for [Fe/H] and [Ni/H], respectively. Our results are not at odds with the possibility that metallicity ([Fe/H]) gradients are steeper in the inner disk (R_GC ~7 - 12 kpc) and flatter towards the outer disk. The open cluster sample studied spans a significant range in age. When breaking the sample into age bins, there is some indication that the younger open cluster population in our sample (log age < 8.7) has a flatter metallicity gradient when compared with the gradients obtained from older open clusters.

The tension between the Hipparcos parallax of the Pleiades and other independent distance estimates continues even after the new reduction of the Hipparcos astrometric data and the development of a new geometric distance measurement for the cluster. A short Pleiades distance from the Hipparcos parallax predicts a number of stars in the solar neighborhood that are sub-luminous at a given photospheric abundance. We test this hypothesis using spectroscopic abundances for a subset of stars in the Hipparcos catalog, which occupy the same region as the Pleiades in the color-magnitude diagram. We derive stellar parameters for 170 nearby G and K type field dwarfs in the Hipparcos catalog based on high-resolution spectra obtained using KPNO 4-m echelle spectrograph. Our analysis shows that, when the Hipparcos parallaxes are adopted, most of our sample stars follow empirical color-magnitude relations. A small fraction of stars are too faint compared to main-sequence fitting relations by $\Delta M_V \geq 0.3$ mag, but the differences are marginal at a $2\sigma$ level partly due to relatively large parallax errors. On the other hand, we find that photometric distances of stars showing signatures of youth as determined from lithium absorption line strengths and $R'_{\rm HK}$ chromospheric activity indices are consistent with the Hipparcos parallaxes. Our result is contradictory to a suggestion that the Pleiades distance from main-sequence fitting is significantly altered by stellar activity and/or the young age of its stars, and provides an additional supporting evidence for the long distance scale of the Pleiades.

We report the discovery of two giant planets orbiting the K giant HD 33844 based on radial velocity data from three independent campaigns. The planets move on nearly circular orbits with semimajor axes $a_b=1.60\pm$0.02 AU and $a_c=2.24\pm$0.05 AU, and have minimum masses (m sin $i$) of $M_b=1.96\pm$0.12 Mjup and $M_c=1.76\pm$0.18 Mjup. Detailed N-body dynamical simulations show that the two planets remain on stable orbits for more than $10^6$ years for low eccentricities, and are most likely trapped in a mutual 3:5 mean-motion resonance.

We present the stellar kinematics across the Galactic bulge and into the disk at positive longitudes from the SDSS-III APOGEE spectroscopic survey of the Milky Way. APOGEE includes extensive coverage of the stellar populations of the bulge along the mid-plane and near-plane regions. From these data, we have produced kinematic maps of 10,000 stars across longitudes 0 deg < l < 65 deg, and primarily across latitudes of |b| < 5 deg in the bulge region. The APOGEE data reveal that the bulge is cylindrically rotating across all latitudes and is kinematically hottest at the very centre of the bulge, with the smallest gradients in both kinematic and chemical space inside the inner-most region (l,|b|) < (5,5) deg. The results from APOGEE show good agreement with data from other surveys at higher latitudes and a remarkable similarity to the rotation and dispersion maps of barred galaxies viewed edge on. The thin bar that is reported to be present in the inner disk within a narrow latitude range of |b| < 2 deg appears to have a corresponding signature in [Fe/H] and [alpha/Fe]. Stars with [Fe/H] > -0.5 have dispersion and rotation profiles that are similar to that of N-body models of boxy/peanut bulges. There is a smooth kinematic transition from the thin bar and boxy bulge (l,|b|) < (15,12) deg out into the disk for stars with [Fe/H] > -1.0, and the chemodynamics across (l,b) suggests the stars in the inner Galaxy with [Fe/H] > -1.0 have an origin in the disk.

Kepler-454 (KOI-273) is a relatively bright (V = 11.69 mag), Sun-like starthat hosts a transiting planet candidate in a 10.6 d orbit. From spectroscopy, we estimate the stellar temperature to be 5687 +/- 50 K, its metallicity to be [m/H] = 0.32 +/- 0.08, and the projected rotational velocity to be v sin i <2.4 km s-1. We combine these values with a study of the asteroseismic frequencies from short cadence Kepler data to estimate the stellar mass to be 1.028+0:04-0:03 M_Sun, the radius to be 1.066 +/- 0.012 R_Sun and the age to be 5.25+1:41-1:39 Gyr. We estimate the radius of the 10.6 d planet as 2.37 +/- 0.13 R_Earth. Using 63 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 36 observations made with the HIRES spectrograph at Keck Observatory, we measure the mass of this planet to be 6.8 +/- 1.4M_Earth. We also detect two additional non-transiting companions, a planet with a minimum mass of 4.46 +/- 0.12 M_J in a nearly circular 524 d orbit and a massive companion with a period >10 years and mass >12.1M_J . The twelve exoplanets with radii <2.7 R_Earth and precise mass measurements appear to fall into two populations, with those <1.6 R_Earth following an Earth-like composition curve and larger planets requiring a significant fraction of volatiles. With a density of 2.76 +/- 0.73 g cm-3, Kepler-454b lies near the mass transition between these two populations and requires the presence of volatiles and/or H/He gas.

We present the $H$-band spectral line lists adopted by the Apache Point Observatory Galactic Evolution Experiment (APOGEE). The APOGEE line lists comprise astrophysical, theoretical, and laboratory sources from the literature, as well as newly evaluated astrophysical oscillator strengths and damping parameters. We discuss the construction of the APOGEE line list, which is one of the critical inputs for the APOGEE Stellar Parameters and Chemical Abundances Pipeline, and present three different versions that have been used at various stages of the project. The methodology for the newly calculated astrophysical line lists is reviewed. The largest of these three line lists contains 134,457 molecular and atomic transitions. In addition to the format adopted to store the data, the line lists are available in MOOG, Synspec and Turbospectrum formats. We also present a list of $H$-band spectral features that are either poorly represented or completely missing in our line list. This list is based on the average of a large number of spectral fit residuals for APOGEE observations spanning a wide range of stellar parameters.

We present griz observations for the clusters M92, M13 and NGC 6791 and gr photometry for M71, Be 29 and NGC 7789. In addition we present new membership identifications for all these clusters, which have been observed spectroscopically as calibrators for the SDSS/SEGUE survey; this paper focuses in particular on the red giant branch stars in the clusters. In a number of cases, these giants were too bright to be observed in the normal SDSS survey operations, and we describe the procedure used to obtain spectra for these stars. For M71, also present a new variable reddening map and a new fiducial for the gr giant branch. For NGC 7789, we derived a transformation from Teff to g-r for giants of near solar abundance, using IRFM Teff measures of stars with good ugriz and 2MASS photometry and SEGUE spectra. The result of our analysis is a robust list of known cluster members with correctly dereddened and (if needed) transformed gr photometry for crucial calibration efforts for SDSS and SEGUE.

We statistically quantify the amount of substructure in the Milky Way stellar halo using a sample of 4568 halo K giant stars at Galactocentric distances ranging over 5-125 kpc. These stars have been selected photometrically and confirmed spectroscopically as K giants from the Sloan Digital Sky Survey's SEGUE project. Using a position-velocity clustering estimator (the 4distance) and a model of a smooth stellar halo, we quantify the amount of substructure in the halo, divided by distance and metallicity. Overall, we find that the halo as a whole is highly structured. We also confirm earlier work using BHB stars which showed that there is an increasing amount of substructure with increasing Galactocentric radius, and additionally find that the amount of substructure in the halo increases with increasing metallicity. Comparing to resampled BHB stars, we find that K giants and BHBs have similar amounts of substructure over equivalent ranges of Galactocentric radius. Using a friends-of-friends algorithm to identify members of individual groups, we find that a large fraction (~33%) of grouped stars are associated with Sgr, and identify stars belonging to other halo star streams: the Orphan Stream, the Cetus Polar Stream, and others, including previously unknown substructures. A large fraction of sample K giants (more than 50%) are not grouped into any substructure. We find also that the Sgr stream strongly dominates groups in the outer halo for all except the most metal-poor stars, and suggest that this is the source of the increase of substructure with Galactocentric radius and metallicity.

We present the discovery of HAT-P-57b, a P = 2.4653 day transiting planet around a V = 10.465 +- 0.029 mag, Teff = 7500 +- 250 K main sequence A8V star with a projected rotation velocity of v sin i = 102.1 +- 1.3 km s^-1. We measure the radius of the planet to be R = 1.413 +- 0.054 R_J and, based on RV observations, place a 95% confidence upper limit on its mass of M < 1.85 M_J . Based on theoretical stellar evolution models, the host star has a mass and radius of 1.47 +- 0.12 M_sun, and 1.500 +- 0.050 R_sun, respectively. Spectroscopic observations made with Keck-I/HIRES during a partial transit event show the Doppler shadow of HAT-P-57b moving across the average spectral line profile of HAT-P- 57, confirming the object as a planetary system. We use these observations, together with analytic formulae that we derive for the line profile distortions, to determine the projected angle between the spin axis of HAT-P-57 and the orbital axis of HAT-P-57b. The data permit two possible solutions, with -16.7 deg < lambda < 3.3 deg or 27.6 deg < lambda < 57.4 deg at 95% confidence, and with relative probabilities for the two modes of 26% and 74%, respectively. Adaptive optics imaging with MMT/Clio2 reveals an object located 2.7" from HAT-P-57 consisting of two point sources separated in turn from each other by 0.22". The H and L -band magnitudes of the companion stars are consistent with their being physically associated with HAT-P-57, in which case they are stars of mass 0.61 +- 0.10 M_sun and 0.53 +- 0.08 M_sun. HAT-P-57 is the most rapidly rotating star, and only the fourth main sequence A star, known to host a transiting planet.

We present precise Doppler observations of WASP-47, a transiting planetary system featuring a hot Jupiter with both inner and outer planetary companions. This system has an unusual architecture and also provides a rare opportunity to measure planet masses in two different ways: the Doppler method, and the analysis of transit-timing variations (TTV). Based on the new Doppler data, obtained with the Planet Finder Spectrograph on the Magellan/Clay 6.5m telescope, the mass of the hot Jupiter is $370 \pm 29~M_{\oplus}$. This is consistent with the previous Doppler determination as well as the TTV determination. For the inner planet WASP-47e, the Doppler data lead to a mass of $12.2\pm 3.7~ M_{\oplus}$, in agreement with the TTV-based upper limit of $<$22~$M_{\oplus}$ ($95\%$ confidence). For the outer planet WASP-47d, the Doppler mass constraint of $10.4\pm 8.4~M_{\oplus}$ is consistent with the TTV-based measurement of $15.2^{+6.7}_{-7.6}~ M_{\oplus}$.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) has built the largest moderately high-resolution (R=22, 500) spectroscopic map of the stars across the Milky Way, and including dust-obscured areas. The APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) is the software developed for the automated analysis of these spectra. ASPCAP determines atmospheric parameters and chemical abundances from observed spectra by comparing observed spectra to libraries of theoretical spectra, using chi-2 minimization in a multidimensional parameter space. The package consists of a fortran90 code that does the actual minimization, and a wrapper IDL code for book-keeping and data handling. This paper explains in detail the ASPCAP components and functionality, and presents results from a number of tests designed to check its performance. ASPCAP provides stellar effective temperatures, surface gravities, and metallicities precise to 2%, 0.1 dex, and 0.05 dex, respectively, for most APOGEE stars, which are predominantly giants. It also provides abundances for up to 15 chemical elements with various levels of precision, typically under 0.1 dex. The final data release (DR12) of the Sloan Digital Sky Survey III contains an APOGEE database of more than 150,000 stars. ASPCAP development continues in the SDSS-IV APOGEE-2 survey.

We present multiwavelength, multi-telescope, ground-based follow-up photometry of the white dwarf WD 1145+017, that has recently been suggested to be orbited by up to six or more, short-period, low-mass, disintegrating planetesimals. We detect 9 significant dips in flux of between 10% and 30% of the stellar flux from our ground-based photometry. We observe transits deeper than 10% on average every ~3.6 hr in our photometry. This suggests that WD 1145+017 is indeed being orbited by multiple, short-period objects. Through fits to the multiple asymmetric transits that we observe, we confirm that the transit egress timescale is usually longer than the ingress timescale, and that the transit duration is longer than expected for a solid body at these short periods, all suggesting that these objects have cometary tails streaming behind them. The precise orbital periods of the planetesimals in this system are unclear from the transit-times, but at least one object, and likely more, have orbital periods of ~4.5 hours. We are otherwise unable to confirm the specific periods that have been reported, bringing into question the long-term stability of these periods. Our high precision photometry also displays low amplitude variations suggesting that dusty material is consistently passing in front of the white dwarf, either from discarded material from these disintegrating planetesimals or from the detected dusty debris disk. For the significant transits we observe, we compare the transit depths in the V- and R-bands of our multiwavelength photometry, and find no significant difference; therefore, for likely compositions the radius of single-size particles in the cometary tails streaming behind the planetesimals in this system must be ~0.15 microns or larger, or ~0.06 microns or smaller, with 2-sigma confidence.

White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.

We explore the radial velocity performance benefits of coupling starlight to a fast-scanning interferometer and a fast-readout spectrometer with zero readout noise. By rapidly scanning an interferometer we can decouple wavelength calibration errors from precise radial velocity measurements, exploiting the advantages of lock-in amplification. In a Bayesian framework, we investigate the correlation between wavelength calibration errors and resulting radial velocity errors. We construct an end-to-end simulation of this approach to address the feasibility of achieving 10 cm/s radial velocity precision on a typical Sun-like star using existing, 5-meter-class telescopes. We find that such a precision can be reached in a single night, opening up possibilities for ground-based detections of Earth-Sun analog systems.

The first supernovae will soon be visible at the edge of the observable universe, revealing the birthplaces of Population III stars. With upcoming near-infrared missions, a broad analysis of the detectability of high-$z$ supernovae is paramount. We combine cosmological and radiation transport simulations, instrument specifications, and survey strategies to create synthetic observations of primeval core-collapse, Type IIn and pair-instability supernovae with the James Webb Space Telescope ($JWST$). We show that a dedicated observational campaign with the $JWST$ can detect up to $\sim 15$ pair-instability explosions, $\sim 300$ core-collapse supernovae, but less than one Type IIn explosion per year, depending on the Population III star formation history. Our synthetic survey also shows that $\approx 1-2 \times10^2$ supernovae detections, depending on the accuracy of the classification, are sufficient to discriminate between a Salpeter and flat mass distribution for high redshift stars with a confidence level greater than 99.5 per cent. We discuss how the purity of the sample affects our results and how supervised learning methods may help to discriminate between CC and PI SNe.

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons at the highest redshifts. In this, the latest of a series of studies on Population III supernovae, we present numerical simulations of 25 - 50 M$_{\odot}$ hypernovae and their light curves done with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible at z = 10 - 15 to the James Webb Space Telescope (JWST) and z = 4 - 5 to the Wide-Field Infrared Survey Telescope (WFIRST), tracing star formation rates in the first galaxies and at the end of cosmological reionization. If, however, the hypernova crashes into a dense shell ejected by its progenitor, it is expected that a superluminous event will occur that may be seen at z ~ 20, in the first generation of stars.

We report a new giant planet orbiting the K giant HD 155233, as well as four stellar-mass companions from the Pan-Pacific Planet Search, a southern hemisphere radial velocity survey for planets orbiting nearby giants and subgiants. We also present updated velocities and a refined orbit for HD 47205b (7 CMa b), the first planet discovered by this survey. HD 155233b has a period of 885$\pm$63 days, eccentricity e=0.03$\pm$0.20, and m sin i=2.0$\pm$0.5 M_jup. The stellar-mass companions range in m sin i from 0.066 M_sun to 0.33 M_sun. Whilst HD 104358B falls slightly below the traditional 0.08 M_sun hydrogen-burning mass limit, and is hence a brown dwarf candidate, we estimate only a 50% a priori probability of a truly substellar mass.

We report the discovery of KELT-4Ab, an inflated, transiting Hot Jupiter orbiting the brightest component of a hierarchical triple stellar system. The host star is an F star with $T_{\rm eff}=6206\pm75$ K, $\log g=4.108\pm0.014$, $\left[{\rm Fe}/{\rm H}\right]=-0.116_{-0.069}^{+0.065}$, ${\rm M_*}=1.201_{-0.061}^{+0.067} \ {\rm M}_{\odot}$, and ${\rm R_*}=1.610_{-0.068}^{+0.078} \ {\rm R}_{\odot}$. The best-fit linear ephemeris is $\rm {BJD_{TDB}} = 2456193.29157 \pm 0.00021 + E\left(2.9895936 \pm 0.0000048\right)$. With a magnitude of $V\sim10$, a planetary radius of $1.699_{-0.045}^{+0.046} \ {\rm R_J}$, and a mass of $0.902_{-0.059}^{+0.060} \ {\rm M_J}$, it is the brightest host among the population of inflated Hot Jupiters ($R_P > 1.5R_J$), making it a valuable discovery for probing the nature of inflated planets. In addition, its existence within a hierarchical triple and its proximity to Earth ($210$ pc) provides a unique opportunity for dynamical studies with continued monitoring with high resolution imaging and precision radial velocities. In particular, the motion of the binary stars around each other and of both stars around the primary star relative to the measured epoch in this work should be detectable when it rises in October 2015.

The Apache Point Observatory Galactic Evolution Experiment (APOGEE), one of the programs in the Sloan Digital Sky Survey III (SDSS-III), has now completed its systematic, homogeneous spectroscopic survey sampling all major populations of the Milky Way. After a three year observing campaign on the Sloan 2.5-m Telescope, APOGEE has collected a half million high resolution (R~22,500), high S/N (>100), infrared (1.51-1.70 microns) spectra for 146,000 stars, with time series information via repeat visits to most of these stars. This paper describes the motivations for the survey and its overall design---hardware, field placement, target selection, operations---and gives an overview of these aspects as well as the data reduction, analysis and products. An index is also given to the complement of technical papers that describe various critical survey components in detail. Finally, we discuss the achieved survey performance and illustrate the variety of potential uses of the data products by way of a number of science demonstrations, which span from time series analysis of stellar spectral variations and radial velocity variations from stellar companions, to spatial maps of kinematics, metallicity and abundance patterns across the Galaxy and as a function of age, to new views of the interstellar medium, the chemistry of star clusters, and the discovery of rare stellar species. As part of SDSS-III Data Release 12, all of the APOGEE data products are now publicly available.

The extended Kepler mission, K2, is now providing photometry of new fields every three months in a search for transiting planets. In a recent study, Foreman-Mackey and collaborators presented a list of 36 planet candidates orbiting 31 stars in K2 Campaign 1. In this contribution, we present stellar and planetary properties for all systems. We combine ground-based seeing-limited survey data and adaptive optics imaging with an automated transit analysis scheme to validate 21 candidates as planets, 17 for the first time, and identify 6 candidates as likely false positives. Of particular interest is K2-18 (EPIC 201912552), a bright (K=8.9) M2.8 dwarf hosting a 2.23 \pm 0.25 R_Earth planet with T_eq = 272 \pm 15 K and an orbital period of 33 days. We also present two new open-source software packages which enable this analysis. The first, isochrones, is a flexible tool for fitting theoretical stellar models to observational data to determine stellar properties using a nested sampling scheme to capture the multimodal nature of the posterior distributions of the physical parameters of stars that may plausibly be evolved. The second is vespa, a new general-purpose procedure to calculate false positive probabilities and statistically validate transiting exoplanets.

The properties of the first galaxies are shaped in large part by the first generations of stars, which emit high energy radiation and unleash both large amounts of mechanical energy and the first heavy elements when they explode as supernovae. We survey the theory of the formation of the first galaxies in this context, focusing on the results of cosmological simulations to illustrate a number of the key processes that define their properties. We first discuss the evolution of the primordial gas as it is incorporated into the earliest galaxies under the influence of the high energy radiation emitted by the earliest stars; we then turn to consider how the injection of heavy elements by the first supernovae transforms the evolution of the primordial gas and alters the character of the first galaxies. Finally, we discuss the prospects for the detection of the first galaxies by future observational missions, in particular focusing on the possibility that primordial star-forming galaxies may be uncovered.

We measure the binarity of detached M-dwarfs in the Kepler field with orbital periods in the range of 1-90 days. Kepler's photometric precision and nearly continuous monitoring of stellar targets over time baselines ranging from 3 months to 4 years make its detection efficiency for eclipsing binaries nearly complete over this period range and for all radius ratios. Our investigation employs a statistical framework akin to that used for inferring planetary occurrence rates from planetary transits. The obvious simplification is that eclipsing binaries have a vastly improved detection efficiency that is limited chiefly by their geometric probabilities to eclipse. For the M-dwarf sample observed by the Kepler Mission, the fractional incidence of eclipsing binaries implies that there are $0.11 ^{+0.02} _{-0.04}$ close stellar companions per apparently single M-dwarf. Our measured binarity is higher than previous inferences of the occurrence rate of close binaries via radial velocity techniques, at roughly the 2$\sigma$ level. This study represents the first use of eclipsing binary detections from a high quality transiting planet mission to infer binary statistics. Application of this statistical framework to the eclipsing binaries discovered by future transit surveys will establish better constraints on short-period M$+$M binary rate, as well as binarity measurements for stars of other spectral types.

Characterizing the physical properties of exoplanets, and understanding their formation and orbital evolution requires precise and accurate knowledge of their host stars. Accurately measuring stellar masses is particularly important because they likely influence planet occurrence and the architectures of planetary systems. Single main-sequence stars typically have masses estimated from evolutionary tracks, which generally provide accurate results due to their extensive empirical calibration. However, the validity of this method for subgiants and giants has been called into question by recent studies, with suggestions that the masses of these evolved stars could have been overestimated. We investigate these concerns using a sample of 59 benchmark evolved stars with model-independent masses (from binary systems or asteroseismology) obtained from the literature. We find very good agreement between these benchmark masses and the ones estimated using evolutionary tracks. The average fractional difference in the mass interval $\sim$0.7 - 4.5 M$_{\odot}$, is consistent with zero (-1.30 $\pm$ 2.42%), with no significant trends in the residuals relative to the input parameters. A good agreement between model-dependent and -independent radii (-4.81 $\pm$ 1.32%) and surface gravities (0.71 $\pm$ 0.51%) is also found. The consistency between independently determined ages for members of binary systems adds further support for the accuracy of the method employed to derive the stellar masses. Taken together, our results indicate that determination of masses of evolved stars using grids of evolutionary tracks is not significantly affected by systematic errors, and is thus valid for estimating the masses of isolated stars beyond the main sequence.

Cosmological simulations of Population (Pop) III star formation suggest that the primordial initial mass function may have extended to sub-solar masses. If Pop III stars with masses < 0.8 M_Sun did form, then they should still be present in the Galaxy today as either main sequence or red giant stars. Despite broad searches, however, no primordial stars have yet been identified. It has long been recognized that the initial metal-free nature of primordial stars could be masked due to accretion of metal-enriched material from the interstellar medium (ISM). Here we point out that while gas accretion from the ISM may readily occur, the accretion of dust from the ISM can be prevented due to the pressure of the radiation emitted from low-mass stars. This implies a possible unique chemical signature for stars polluted only via accretion, namely an enhancement in gas phase elements relative to those in the dust phase. Using Pop III stellar models, we outline the conditions in which this signature could be exhibited, and we derive the expected signature for the case of accretion from the local ISM. Intriguingly, due to the large fraction of iron depleted into dust relative to that of carbon and other elements, this signature is similar to that observed in many of the so-called carbon-enhanced metal-poor (CEMP) stars. We therefore suggest that some fraction of the observed CEMP stars may, in fact, be accretion-polluted Pop III stars.

We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD219134. This is the first result of the Rocky Planet Search program with HARPS-N on the TNG in La Palma. The inner planet orbits the star in 3.0937 +/-0.0004 days, on a quasi-circular orbit with a semi-major axis of 0.0382 +/- 0.0003 AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 +/- 0.24 m/s) and observed depth of the transit (359 +/- 38 ppm), the planet mass and radius are estimated to be 4.46 +/- 0.47 M_{\oplus} and 1.606 +/- 0.086 R_{\oplus} leading to a mean density of 5.89 +/- 1.17 g/cc, suggesting a rocky composition. One additional planet with minimum mass of 2.67 +/- 0.59 M_{\oplus} moves on a close-in, quasi-circular orbit with a period of 6.765 +/- 0.005 days. The third planet in the system has a period of 46.78 +/- 0.16 days and a minimum mass of 8.7 +/- 1.1 M{\oplus}, at 0.234 +/- 0.002 AU from the star. Its eccentricity is 0.32 +/- 0.14. The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 +/- 0.1 days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 +/- 6 M_{\oplus} orbits the star in 1190 days, on an eccentric orbit (e=0.27 +/- 0.11) at a distance of 2.14 +/- 0.27 AU.

NASA's re-purposed Kepler mission -- dubbed K2 -- has brought new scientific opportunities that were not anticipated for the original Kepler mission. One science goal that makes optimal use of K2's capabilities, in particular its 360-degree ecliptic field of view, is galactic archaeology -- the study of the evolution of the Galaxy from the fossil stellar record. The thrust of this research is to exploit high-precision, time-resolved photometry from K2 in order to detect oscillations in red giant stars. This asteroseismic information can provide estimates of stellar radius (hence distance), mass and age of vast numbers of stars across the Galaxy. Here we present the initial analysis of a subset of red giants, observed towards the North Galactic Gap, during the mission's first full science campaign. We investigate the feasibility of using K2 data for detecting oscillations in red giants that span a range in apparent magnitude and evolutionary state (hence intrinsic luminosity). We demonstrate that oscillations are detectable for essentially all cool giants within the $\log g$ range $\sim 1.9-3.2$. Our detection is complete down to $\mathit{Kp}\sim 14.5$, which results in a seismic sample with little or no detection bias. This sample is ideally suited to stellar population studies that seek to investigate potential shortcomings of contemporary Galaxy models.

We measure the two-point correlation function of G-dwarf stars within 1-3 kpc of the Sun in multiple lines-of-sight using the Schlesinger et al. G-dwarf sample from the SDSS SEGUE survey. The shapes of the correlation functions along individual SEGUE lines-of-sight depend sensitively on both the stellar-density gradients and the survey geometry. We fit smooth disk galaxy models to our SEGUE clustering measurements, and obtain strong constraints on the thin- and thick-disk components of the Milky Way. Specifically, we constrain the values of the thin- and thick-disk scale heights with 3% and 2% precision, respectively, and the values of the thin- and thick-disk scale lengths with 20% and 8% precision, respectively. Moreover, we find that a two-disk model is unable to fully explain our clustering measurements, which exhibit an excess of clustering at small scales (< 50 pc). This suggests the presence of small-scale substructure in the disk system of the Milky Way.

The MINiature Exoplanet Radial Velocity Array (MINERVA) is a US-based observational facility dedicated to the discovery and characterization of exoplanets around a nearby sample of bright stars. MINERVA employs a robotic array of four 0.7 m telescopes outfitted for both high-resolution spectroscopy and photometry, and is designed for completely autonomous operation. The primary science program is a dedicated radial velocity survey and the secondary science objective is to obtain high precision transit light curves. The modular design of the facility and the flexibility of our hardware allows for both science programs to be pursued simultaneously, while the robotic control software provides a robust and efficient means to carry out nightly observations. In this article, we describe the design of MINERVA including major hardware components, software, and science goals. The telescopes and photometry cameras are characterized at our test facility on the Caltech campus in Pasadena, CA, and their on-sky performance is validated. New observations from our test facility demonstrate sub-mmag photometric precision of one of our radial velocity survey targets, and we present new transit observations and fits of WASP-52b -- a known hot-Jupiter with an inflated radius and misaligned orbit. The facility is now in the process of being relocated to its final destination at the Fred Lawrence Whipple Observatory in southern Arizona, and science operations will begin in 2015.

We derive age constraints for 1639 red giants in the APOKASC sample for which seismic parameters from Kepler, as well as effective temperatures, metallicities and [alpha/Fe] values from APOGEE DR12 are available. We investigate the relation between age and chemical abundances for these stars, using a simple and robust approach to obtain ages. We first derive stellar masses using standard seismic scaling relations, then determine the maximum possible age for each star as function of its mass and metallicity, independently of its evolutionary stage. While the overall trend between maximum age and chemical abundances is a declining fraction of young stars with increasing [alpha/Fe], at least 14 out of 241 stars with [alpha/Fe]>0.13 are younger than 6 Gyr. Five stars with [alpha/Fe]>0.2 have ages below 4 Gyr. We examine the effect of modifications in the standard seismic scaling relations, as well as the effect of very low helium fractions, but these changes are not enough to make these stars as old as usually expected for alpha-rich stars (i.e., ages greater than 8-9 Gyr). Such unusual alpha-rich young stars have also been detected by other surveys, but defy simple explanations in a galaxy evolution context.

Light curves of microlensing events involving stellar binaries and planetary systems can provide information about the orbital elements of the system due to orbital modulations of the caustic structure. Accurately measuring the orbit in either the stellar or planetary case requires detailed modeling of subtle deviations in the light curve. At the same time, the natural, Cartesian parameterization of a microlensing binary is partially degenerate with the microlens parallax. Hence, it is desirable to perform independent tests of the predictions of microlens orbit models using radial velocity time series of the lens binary system. To this end, we present 3.5 years of RV monitoring of the binary lens system OGLE-2009-BLG-020L, for which Skowron et al. (2011) constrained all internal parameters of the 200--700 day orbit. Our RV measurements reveal an orbit that is consistent with the predictions of the microlens light curve analysis, thereby providing the first confirmation of orbital elements inferred from microlensing events.

From pulsating stars to transiting exoplanets, the search for periodic signals in K2 data, Kepler's 2-wheeled extension, is relevant to a long list of scientific goals. Systematics affecting K2 light curves due to the decreased spacecraft pointing precision inhibit the easy extraction of periodic signals from the data. We here develop a method for producing periodograms of K2 light curves that are insensitive to pointing-induced systematics; the Systematics-Insensitive Periodogram (SIP). Traditional sine-fitting periodograms use a generative model to find the frequency of a sinusoid that best describes the data. We extend this principle by including systematic trends, based on a set of 'Eigen light curves', following Foreman-Mackey et al. (2015), in our generative model as well as a sum of sine and cosine functions over a grid of frequencies. Using this method we are able to produce periodograms with vastly reduced systematic features. The quality of the resulting periodograms are such that we can recover acoustic oscillations in giant stars and measure stellar rotation periods without the need for any detrending. The algorithm is also applicable to the detection of other periodic phenomena such as variable stars, eclipsing binaries and short-period exoplanet candidates. The SIP code is available at https://github.com/RuthAngus/SIPK2.