We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the Initial-Final Mass Relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color-magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub and ReadtheDocs, respectively.

The new NASA Astrophysics Data System (ADS) is designed with a serviceoriented architecture (SOA) that consists of multiple customized Apache Solr search engine instances plus a collection of microservices, containerized using Docker, and deployed in Amazon Web Services (AWS). For complex systems, like the ADS, this loosely coupled architecture can lead to a more scalable, reliable and resilient system if some fundamental questions are addressed. After having experimented with different AWS environments and deployment methods, we decided in December 2017 to go with Kubernetes as our container orchestration. Defining the best strategy to properly setup Kubernetes has shown to be challenging: automatic scaling services and load balancing traffic can lead to errors whose origin is difficult to identify, monitoring and logging the activity that happens across multiple layers for a single request needs to be carefully addressed, and the best workflow for a Continuous Integration and Delivery (CI/CD) system is not self-evident. We present here how we tackle these challenges and our plans for the future.

astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.

OSIRIS is a near-infrared (1.0--2.4 $\mu$m) integral field spectrograph operating behind the adaptive optics system at Keck Observatory, and is one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics (LGS AO) papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP including characterization of the flux assignment algorithm; exploration of spatial rippling in the reduced data cubes; and improvements to several calibration files, including the rectification matrix, the bad pixel mask, and the wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

We present a kinematic study of the nuclear stellar disk in M31 at infrared wavelengths using high spatial resolution integral field spectroscopy. The spatial resolution achieved, FWHM = 0."12 (0.45 pc at the distance of M31), has only previously been equaled in spectroscopic studies by space-based long-slit observations. Using adaptive optics-corrected integral field spectroscopy from the OSIRIS instrument at the W. M. Keck Observatory, we map the line-of-sight kinematics over the entire old stellar eccentric disk orbiting the supermassive black hole (SMBH) at a distance of r<4 pc. The peak velocity dispersion is 381+/-55 km/s , offset by 0.13 +/- 0.03 from the SMBH, consistent with previous high-resolution long-slit observations. There is a lack of near-infrared (NIR) emission at the position of the SMBH and young nuclear cluster, suggesting a spatial separation between the young and old stellar populations within the nucleus. We compare the observed kinematics with dynamical models from Peiris & Tremaine (2003). The best-fit disk orientation to the NIR flux is [$\theta_l$, $\theta_i$, $\theta_a$] = [-33 +/- 4$^{\circ}$, 44 +/- 2$^{\circ}$, -15 +/- 15$^{\circ}$], which is tilted with respect to both the larger-scale galactic disk and the best-fit orientation derived from optical observations. The precession rate of the old disk is $\Omega_P$ = 0.0 +/- 3.9 km/s/pc, lower than the majority of previous observations. This slow precession rate suggests that stellar winds from the disk will collide and shock, driving rapid gas inflows and fueling an episodic central starburst as suggested in Chang et al. (2007).

We present a high spatial resolution optical and infrared study of the circumnuclear region in Arp 220, a late-stage galaxy merger. Narrowband imaging using HST/WFC3 has resolved the previously observed peak in H$\alpha$+[NII] emission into a bubble-shaped feature. This feature measures 1.6" in diameter, or 600 pc, and is only 1" northwest of the western nucleus. The bubble is aligned with the western nucleus and the large-scale outflow axis seen in X-rays. We explore several possibilities for the bubble origin, including a jet or outflow from a hidden active galactic nucleus (AGN), outflows from high levels of star formation within the few hundred pc nuclear gas disk, or an ultraluminous X-ray source. An obscured AGN or high levels of star formation within the inner $\sim$100 pc of the nuclei are favored based on the alignment of the bubble and energetics arguments.

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

We test analytic predictions from different models of magnetospheric accretion, which invoke disk-locking, using stellar and accretion parameters derived from models of low resolution optical spectra of 36 T Tauri stars (TTSs) in NGC 2264 (age~3 Myrs). Little evidence is found for models that assume purely dipolar field geometries; however, strong support is found in the data for a modified version of the X-wind model (Shu et al. 1994) which allows for non-dipolar field geometries. The trapped flux concept in the X-wind model is key to making the analytic predictions which appear supported in the data. By extension, our analysis provides support for the outflows predicted by the X-wind as these also originate in the trapped flux region. In addition, we find no support in the data for accretion powered stellar winds from young stars. By comparing the analysis presented here of NGC 2264 with a similar analysis of stars in Taurus (age~1-2 Myr), we find evidence that the equilibrium interaction between the magnetic field and accretion disk in TTS systems evolves as the stars grow older, perhaps as the result of evolution of the stellar magnetic field geometry. We compare the accretion rates we derive with accretion rates based on U-band excess, finding good agreement. In addition, we use our accretion parameters to determine the relationship between accretion and H-beta luminosity, again finding good agreement with previously published results; however, we also find that care must be used when applying this relationship due to strong chromospheric emission in young stars which can lead to erroneous results in some cases.

We present a study of the stellar and circumstellar properties of Class I sources using low-resolution (R~1000) near-infrared K- and L-band spectroscopy. We measure prominent spectral lines and features in 8 objects and use fits to standard star spectra to determine spectral types, visual extinctions, K-band excesses, and water ice optical depths. Four of the seven systems studied are close binary pairs; only one of these systems, Haro 6-10, was angularly resolvab le. For certain stars some properties found in our analysis differ substantially from published values; we analyze the origin of these differences. We determine extinction to each source using three different methods and compare and discuss the resulting values. One hypothesis that we were testing, that extinction dominates over the K-band excess in obscuration of the stellar photospheric absorption lines, appears not to be true. Accretion luminosities and mass accretion rates calculated for our targets are highly uncertain, in part the result of our inexact knowledge of extinction. For the six targets we were able to place on an H-R diagram, our age estimates, <2 Myr, are somewhat younger than those from comparable studies. Our results underscore the value of low-resolution spectroscopy in the study of protostars and their environments; however, the optimal approach to the study of Class I sources likely involves a combination of high- and low-resolution near-infrared, mid-infrared, and millimeter wavelength observations. Accurate and precise measurements of extinction in Class I protostars will be key to improving our understanding of these objects.