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 hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.5 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*.

We report new infrared measurements of the supermassive black hole at the Galactic Center, Sgr A*, over a decade that was previously inaccessible at these wavelengths. This enables a variability study that addresses variability timescales that are ten times longer than earlier published studies. Sgr A* was initially detected in the near-infrared with adaptive optics observations in 2002. While earlier data exists in form of speckle imaging (1995 - 2005), Sgr A* was not detected in the initial analysis. Here, we improved our speckle holography analysis techniques. This has improved the sensitivity of the resulting speckle images by up to a factor of three. Sgr A* is now detectable in the majority of epochs covering 7 years. The brightness of Sgr A* in the speckle data has an average observed K magnitude of 16.0, which corresponds to a dereddened flux density of $3.4$ mJy. Furthermore, the flat power spectral density (PSD) of Sgr A* between $\sim$80 days and 7 years shows its uncorrelation in time beyond the proposed single power-law break of $\sim$245 minutes. We report that the brightness and its variability is consistent over 22 years. This analysis is based on simulations using Witzel et al. (2018) model to characterize infrared variability from 2006 to 2016. Finally, we note that the 2001 periapse of the extended, dusty object G1 had no apparent effect on the near-infrared emission from accretion flow onto Sgr A*. The result is consistent with G1 being a self-gravitating object rather than a disrupting gas cloud.

The electromagnetic counterpart to the Galactic center supermassive black hole, Sgr A*, has been observed in the near-infrared for over 20 years and is known to be highly variable. We report new Keck Telescope observations showing that Sgr A* reached much brighter flux levels in 2019 than ever measured at near-infrared wavelengths. In the K$^\prime$ band, Sgr A* reached flux levels of $\sim6$ mJy, twice the level of the previously observed peak flux from $>13,000$ measurements over 130 nights with the VLT and Keck Telescopes. We also observe a factor of 75 change in flux over a 2-hour time span with no obvious color changes between 1.6 $\mu$m and 2.1 $\mu$m. The distribution of flux variations observed this year is also significantly different than the historical distribution. Using the most comprehensive statistical model published, the probability of a single night exhibiting peak flux levels observed this year, given historical Keck observations, is less than $0.3\%$. The probability to observe the flux levels similar to all 4 nights of data in 2019 is less than $0.05\%$. This increase in brightness and variability may indicate a period of heightened activity from Sgr A* or a change in its accretion state. It may also indicate that the current model is not sufficient to model Sgr A* at high flux levels and should be updated. Potential physical origins of Sgr A*'s unprecedented brightness may be from changes in the accretion-flow as a result of the star S0-2's closest passage to the black hole in 2018 or from a delayed reaction to the approach of the dusty object G2 in 2014. Additional multi-wavelength observations will be necessary to both monitor Sgr A* for potential state changes and to constrain the physical processes responsible for its current variability.

The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3\sigma$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale.

Nearby dwarf galaxies are local analogues of high-redshift and metal-poor stellar populations. Most of these systems ceased star formation long ago, but they retain signatures of their past that can be unraveled by detailed study of their resolved stars. Archaeological examination of dwarf galaxies with resolved stellar spectroscopy provides key insights into the first stars and galaxies, galaxy formation in the smallest dark matter halos, stellar populations in the metal-free and metal-poor universe, the nature of the first stellar explosions, and the origin of the elements. Extremely large telescopes with multi-object R=5,000-30,000 spectroscopy are needed to enable such studies for galaxies of different luminosities throughout the Local Group.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

The stellar initial mass function (IMF) is a fundamental property of star formation, offering key insight into the physics driving the process as well as informing our understanding of stellar populations, their by-products, and their impact on the surrounding medium. While the IMF appears to be fairly uniform in the Milky Way disk, it is not yet known how the IMF might behave across a wide range of environments, such as those with extreme gas temperatures and densities, high pressures, and low metallicities. We discuss new opportunities for measuring the IMF in such environments in the coming decade with JWST, WFIRST, and thirty-meter class telescopes. For the first time, we will be able to measure the high-mass slope and peak of the IMF via direct star counts for massive star clusters across the Milky Way and Local Group, providing stringent constraints for star formation theory and laying the groundwork for understanding distant and unresolved stellar systems.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

We present improved relative astrometry for stars within the central half parsec of our Galactic Center based on data obtained with the 10 m W. M. Keck Observatory from 1995 to 2017. The new methods used to improve the astrometric precision and accuracy include correcting for local astrometric distortions, applying a magnitude dependent additive error, and more carefully removing instances of stellar confusion. Additionally, we adopt jackknife methods to calculate velocity and acceleration uncertainties. The resulting median proper motion uncertainty is 0.05 mas/yr for our complete sample of 1184 stars in the central 10'' (0.4 pc). We have detected 24 accelerating sources, 2.6 times more than the number of previously published accelerating sources, which extend out to 4'' (0.16 pc) from the black hole. Based on S0-2's orbit, our new astrometric analysis has reduced the systematic error of the supermassive black hole (SMBH) by a factor of 2. The linear drift in our astrometric reference frame is also reduced in the North-South direction by a factor of 4. We also find the first potential astrometric binary candidate S0-27 in the Galactic center. These astrometric improvements provide a foundation for future studies of the origin and dynamics of the young stars around the SMBH, the structure and dynamics of the old nuclear star cluster, the SMBH's properties derived from orbits, and tests of General Relativity (GR) in a strong gravitational field.

Precision measurements of the stars in short-period orbits around the supermassive black hole at the Galactic Center are now being used to constrain general relativistic effects, such as the gravitational redshift and periapse precession. One of the largest systematic uncertainties in the measured orbits has been errors in the astrometric reference frame, which is derived from seven infrared-bright stars associated with SiO masers that have extremely accurate radio positions, measured in the Sgr A*-rest frame. We have improved the astrometric reference frame within 14'' of the Galactic Center by a factor of 2.5 in position and a factor of 5 in proper motion. In the new reference frame, Sgr A* is localized to within a position of 0.645 mas and proper motion of 0.03 mas/yr. We have removed a substantial rotation (2.25 degrees per decade), that was present in the previous less-accurate reference frame used to measure stellar orbits in the field. With our improved methods and continued monitoring of the masers, we predict that orbital precession predicted by General Relativity will become detectable in the next ~5 years.

We present a $\approx 11.5$ year adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central $\sim 0.4$ pc ($\sim 10''$) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager ($K'$-band, $\lambda_0 = 2.124 \text{ } \mu \text{m}$). We achieve a photometric uncertainty floor of $\Delta m_{K'} \sim 0.03$ ($\approx 3\%$), comparable to the highest precision achieved in other AO studies. Approximately half of our sample ($50 \pm 2 \%$) shows variability. $52 \pm 5\%$ of known early-type young stars and $43 \pm 4 \%$ of known late-type giants are variable. These variability fractions are higher than those of other young, massive star populations or late-type giants in globular clusters, and can be largely explained by two factors. First, our experiment time baseline is sensitive to long-term intrinsic stellar variability. Second, the proper motion of stars behind spatial inhomogeneities in the foreground extinction screen can lead to variability. We recover the two known Galactic center eclipsing binary systems: IRS 16SW and S4-258 (E60). We constrain the Galactic center eclipsing binary fraction of known early-type stars to be at least $2.4 \pm 1.7\%$. We find no evidence of an eclipsing binary among the young S-stars nor among the young stellar disk members. These results are consistent with the local OB eclipsing binary fraction. We identify a new periodic variable, S2-36, with a 39.43 day period. Further observations are necessary to determine the nature of this source.

As a young massive cluster in the Central Molecular Zone, the Arches cluster is a valuable probe of the stellar Initial Mass Function (IMF) in the extreme Galactic Center environment. We use multi-epoch Hubble Space Telescope observations to obtain high-precision proper motion and photometric measurements of the cluster, calculating cluster membership probabilities for stars down to 1.8 M$_{\odot}$ between cluster radii of 0.25 pc -- 3.0 pc. We achieve a cluster sample with just ~8% field contamination, a significant improvement over photometrically-selected samples which are severely compromised by the differential extinction across the field. Combining this sample with K-band spectroscopy of 5 cluster members, we forward model the Arches cluster to simultaneously constrain its IMF and other properties (such as age and total mass) while accounting for observational uncertainties, completeness, mass segregation, and stellar multiplicity. We find that the Arches IMF is best described by a 1-segment power law that is significantly top-heavy: $\alpha$ = 1.80 $\pm$ 0.05 (stat) $\pm$ 0.06 (sys), where dN/dm $\propto$ m$^{-\alpha}$, though we cannot discount a 2-segment power law model with a high-mass slope only slightly shallower than local star forming regions ($\alpha$ = 2.04$^{+0.14}_{-0.19}$ $\pm$ 0.04) but with a break at 5.8$^{+3.2}_{-1.2}$ $\pm$ 0.02 M$_{\odot}$. In either case, the Arches IMF is significantly different than the standard IMF. Comparing the Arches to other young massive clusters in the Milky Way, we find tentative evidence for a systematically top-heavy IMF at the Galactic Center.

A precise extinction law is a critical input when interpreting observations of highly reddened sources such as young star clusters and the Galactic Center (GC). We use Hubble Space Telescope observations of a region of moderate extinction and a region of high extinction to measure the optical and near-infrared extinction law (0.8 $\mu$m -- 2.2 $\mu$m). The moderate extinction region is the young massive cluster Westerlund 1 (Wd1; A$_{Ks} \sim$ 0.6 mag), where 453 proper motion-selected main-sequence stars are used to measure the shape of the extinction law. To quantify the shape we define the parameter $\mathcal{S}_{1/\lambda}$, which behaves similarly to a color excess ratio but is continuous as a function of wavelength. The high extinction region is the GC (A$_{Ks} \sim$ 2.5 mag), where 819 red clump stars are used to determine the normalization of the law. The best-fit extinction law is able to reproduce the Wd1 main sequence colors, which previous laws misestimate by 10%-30%. The law is inconsistent with a single power law, even when only the near-infrared filters are considered, and has A$_{F125W}$/A$_{Ks}$ and A$_{F814W}$/A$_{Ks}$ values that are 18% and 24% larger than the commonly used \citet{Nishiyama:2009fc} law, respectively. Using the law we recalculate the Wd1 distance to be 3896 $\pm$ 328 pc from published observations of eclipsing binary W13. This new extinction law should be used for highly reddened populations in the Milky Way, such as the Quintuplet cluster and Young Nuclear Cluster. A python code is provided to generate the law for future use.

At a projected distance of ~26 pc from Sgr A*, the Arches cluster provides insight to star formation in the extreme Galactic Center (GC) environment. Despite its importance, many key properties such as the cluster's internal structure and orbital history are not well known. We present an astrometric and photometric study of the outer region of the Arches cluster (R > 6.25") using HST WFC3IR. Using proper motions we calculate membership probabilities for stars down to F153M = 20 mag (~2.5 M_sun) over a 120" x 120" field of view, an area 144 times larger than previous astrometric studies of the cluster. We construct the radial profile of the Arches to a radius of 75" (~3 pc at 8 kpc), which can be well described by a single power law. From this profile we place a 3-sigma lower limit of 2.8 pc on the observed tidal radius, which is larger than the predicted tidal radius (1 - 2.5 pc). Evidence of mass segregation is observed throughout the cluster and no tidal tail structures are apparent along the orbital path. The absence of breaks in the profile suggests that the Arches has not likely experienced its closest approach to the GC between ~0.2 - 1 Myr ago. If accurate, this constraint indicates that the cluster is on a prograde orbit and is located front of the sky plane that intersects Sgr A*. However, further simulations of clusters in the GC potential are required to interpret the observed profile with more confidence.

Low resolution ESO VLT/FORS spectra of blue supergiant stars are analyzed to determine stellar metallicities (based on elements such as Fe, Ti, Mg) in the extended disk of the spiral galaxy NGC3621. Mildly subsolar metallicity (-0.30 dex) is found for the outer objects beyond 7 kpc independent of galactocentric radius and compatible with the absence of a metallicity gradient confirming the results of a recent investigation of interstellar medium HII region gas oxygen abundances. The stellar metallicities are slightly higher than those from the HII regions when based on measurements of the weak forbidden auroral oxygen line at 4363 AE but lower than the ones obtained with the R23 strong line method. It is shown that the present level of metallicity in the extended disk cannot be the result of chemical evolution over the age of the disk with the present rate of in situ star formation. Additional mechanisms must be involved. In addition to metallicity, stellar effective temperatures, gravities, interstellar reddening, and bolometric magnitudes are determined. After application of individual reddening corrections for each target the flux-weighted gravity-luminosity relationship of blue supergiant stars is used to obtain a distance modulus of 29.07+/-0.09 mag (distance D=6.52+/-0.28 Mpc). This new distance is discussed in relation to Cepheid and tip of the red giant branch distances.

We present a quantitative analysis of the low-resolution (4.5 A) spectra of 12 late-B and early-A blue supergiants (BSGs) in the metal-poor dwarf galaxy NGC 3109. A modified method of analysis is presented which does not require use of the Balmer jump as an independent temperature indicator, as used in previous studies. We determine stellar effective temperatures, gravities, metallicities, reddening, and luminosities, and combine our sample with the early-B type BSGs analyzed by Evans et al. (2007) to derive the distance to NGC 3109 using the Flux-weighted Gravity-Luminosity Relation (FGLR). Using primarily Fe-group elements, we find an average metallicity of [Z] = -0.67 +/- 0.13, and no evidence of a metallicity gradient in the galaxy. Our metallicities are higher than those found by Evans et al. (2007) based on the oxygen abundances of early-B supergiants ([Z] = -0.93 +/- 0.07), suggesting a low alpha/Fe ratio for the galaxy. We adjust the position of NGC 3109 on the BSG-determined galaxy mass-metallicity relation accordingly and compare it to metallicity studies of HII regions in star-forming galaxies. We derive an FGLR distance modulus of 25.55 +/- 0.09 (1.27 Mpc) that compares well with Cepheid and tip of the red giant branch (TRGB) distances. The FGLR itself is consistent with those found in other galaxies, demonstrating the reliability of this method as a measure of extragalactic distances.

We present the first direct determination of a stellar metallicity in the spiral galaxy NGC4258 (D=7.6 Mpc) based on the quantitative analysis of a low-resolution (~5 AE) Keck LRIS spectrum of a blue supergiant star located in its disk. A determination of stellar metallicity in this galaxy is important for the absolute calibration of the Cepheid Period-Luminosity relation as an anchor for the extragalactic distance scale and for a better characterization of its dependence as a function of abundance. We find a value 0.2 dex lower than solar metallicity at a galactocentric distance of 8.7 kpc, in agreement with recent HII region studies using the weak forbidden auroral oxygen line at 4363 AE. We determine the effective stellar temperature, gravity, luminosity and line-of-sight extinction of the blue supergiant being studied. We show that it fits well on the flux-weighted gravity--luminosity relation (FGLR), strengthening the potential of this method as a new extragalactic distance indicator.