We analyzed multiwavelength observations of the previously identified Galactic center X-ray binary CXO 174528.79-290942.8 (XID 6592) and determine that the near-infrared counterpart is a red supergiant based on its spectrum and luminosity. Scutum X-1 is the only previously known X-ray binary with a red supergiant donor star and closely resembles XID 6592 in terms of X-ray luminosity (L$_{\mathrm{X}}$), absolute magnitude, and IR variability (L$_{\mathrm{IR,var}}$), supporting the conclusion that XID 6592 contains a red supergiant donor star. The XID 6592 infrared counterpart shows variability of $\sim$0.5 mag in the Wide-field Infrared Survey Explorer-1 band (3.4 $\mu$m) on timescales of a few hours. Other infrared data sets also show large-amplitude variability from this source at earlier epochs but do not show significant variability in recent data. We do not expect red supergiants to vary by $\sim50\%$ in luminosity over these short timescales, indicating that the variability should be powered by the compact object. However, the X-ray luminosity of this system is typically $\sim1000\times$ less than the variable luminosity in the infrared and falls below the Chandra detection limit. While X-ray reprocessing can produce large-amplitude fast infrared variability, it typically requires L$_{{\mathrm{X}}} >> $ L$_{\mathrm{IR,var}}$ to do so, indicating that another process must be at work. We suggest that this system may be a supergiant fast X-ray transient (SFXT), and that a large ($\sim10^{38}$ ergs s$^{-1}$), fast ($10^{2-4}$ s) X-ray flare could explain the rapid IR variability and lack of a long-lasting X-ray outburst detection. SFXTs are normally associated with blue supergiant companions, so if confirmed, XID 6592 would be the first red supergiant SFXT, as well as the second X-ray red supergiant binary.

We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $\nu_1$ and $\nu_3$ vibrational bands at 2.7 $\mu$m, from the $\nu_2$ vibrational band at 6.1 $\mu$m, and from pure rotational transitions near 10-13 $\mu$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

In 2019 October Betelgeuse began a decline in V-band brightness that went beyond the minimum expected from its quasi-periodic ~420 day cycle, becoming the faintest in recorded photometric history. Observations obtained in 2019 December with VLT/SPHERE (Montarges 2020) have shown that the southern half of the star has become markedly fainter than in 2019 January indicating that a major change has occurred in, or near, the photosphere. We present SOFIA-EXES high spectral-resolution observations of [Fe II] 25.99 mic and [S I] 25.25 mic emission lines from Betelgeuse obtained during the unprecedented 2020 February V-band brightness minimum to investigate potential changes in the circumstellar flow. These spectra are compared to observations obtained in 2015 and 2017 when the V magnitude was typical of brighter phases. We find only very small changes in the gas velocities reflected by either of the line profiles, no significant changes in the flux to continuum ratios, and hence no significant changes in the [Fe II]/[S I] flux ratios. There is evidence that absorption features have appeared in the 2020 continuum. The Alfven wave-crossing time from the upper-photosphere is sufficiently long that one would not expect a change in the large scale magnetic field to reach the circumstellar [Fe II] and [S I] line forming regions, 3 < R(R*) < 20. However, the light-crossing time is of order a few hours and a reduction in luminosity may reduce the dust-gas heating rate and [O I] 63.19 mic emission which has contributions from R > 20R*, where significant circumstellar oxygen-rich dust is observed.

Supernovae from core-collapse of massive stars drive shocks into the molecular clouds from which the stars formed. Such shocks affect future star formation from the molecular clouds, and the fast-moving, dense gas with compressed magnetic fields is associated with enhanced cosmic rays. This paper presents new theoretical modeling, using the Paris-Durham shock model, and new observations, using the Stratospheric Observatory for Infrared Astronomy (SOFIA), of the H$_2$ S(5) pure rotational line from molecular shocks in the supernova remnant IC443. We generate MHD models for non-steady-state shocks driven by the pressure of the IC443 blast wave into gas of densities $10^3$ to $10^5$ cm$^{-3}$. We present the first detailed derivation of the shape of the velocity profile for emission from H$_2$ lines behind such shocks, taking into account the shock age, preshock density, and magnetic field. For preshock densities $10^3$-$10^5$ cm$^{-3}$, the the predicted shifts of line centers, and the line widths, of the H$_2$ lines range from 20-2, and 30-4 km/s, respectively. The a priori models are compared to the observed line profiles, showing that clumps C and G can be explained by shocks into gas with density 10$^3$ to $2\times 10^4$ cm$^{-3}$ and strong magnetic fields. For clump B2 (a fainter region near clump B), the H$_2$ spectrum requires a J-type shock into moderate density (~100 cm$^{-3}$) with the gas accelerated to 100 km/s from its pre-shock location. Clump B1 requires both a magnetic-dominated C-type shock (like for clumps C and G) and a J-type shock (like for clump B1) to explain the highest observed velocities. The J-type shocks that produce high-velocity molecules may be locations where the magnetic field is nearly parallel to the shock velocity, which makes it impossible for a C-type shock (with ions and neutrals separated) to form.

Spectrally-resolved observations of three pure rotational lines of H$_2$, conducted with the EXES instrument on SOFIA toward the classic bow shock HH7, reveal systematic velocity shifts between the S(5) line of ortho-H$_2$ and the two para-H$_2$ lines [S(4) and S(6)] lying immediately above and below it on the rotational ladder. These shifts, reported here for the first time, imply that we are witnessing the conversion of para-H$_2$ to ortho-H$_2$ within a shock wave driven by an outflow from a young stellar object. The observations are in good agreement with the predictions of models for non-dissociative, C-type molecular shocks. They provide a clear demonstration of the chemical changes wrought by interstellar shock waves, in this case the conversion of para-H$_2$ to ortho-H$_2$ in reactive collisions with atomic hydrogen, and provide among the most compelling evidence yet obtained for C-type shocks in which the flow velocity changes continuously.

We present mid-infrared SOFIA/EXES spectroscopy of Europa, seeking direct evidence of the presence of water vapor arising from plumes venting from the surface of Europa. We place quantitatively useful upper limits on the strength of water vibrational-rotational emission lines. Conversion to water mass limits is dependent on the rotational temperature of the vapor. For low rotational temperature, the limits lie below the inferred water mass from previous HST plume observations. For higher temperatures, the limits are comparable. We also present coordinated HST transit observations obtained close in time to the SOFIA observations. There is evidence for a feature close to the location of the previously seen feature north of the crater Pwyll in one of the HST images, although it was not observable by EXES given its location. We conclude that if a water plume had been active at the time of the SOFIA observation, with the strength implied by previous HST observations, then under the right Earth atmospheric and geometric conditions, the plume could have been detected by EXES, however no IR water vibrational-rotational emission was detected.

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $\nu_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

We have performed a 5-8 $\mu$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $\mu$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $\mu$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance.

Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in molecular gas. If there is dense gas that has been retained in close proximity to BN, it may be evidence that the latter scenario is favored. We observed BN at high spectral resolution in three windows near 6 $\mu$m using SOFIA/EXES targeting the $\nu_2$ vibrational band of H$_2$O. Absorption from only three transitions of H$_2$O is detected, and through kinematic analysis is associated with cool, dense foreground gas, not BN itself. We find no evidence for H$_2$O absorption or emission at the systemic velocity of BN.

We present high resolution spectroscopy of mid-infrared molecular emission from two very active T Tauri stars, AS 205 N and DR Tau. In addition to measuring high signal-to-noise line profiles of water, we report the first spectrally resolved mid-infrared line profiles of HCN emission from protoplanetary disks. The similar line profiles and temperatures of the HCN and water emission indicate that they arise in the same volume of the disk atmosphere, within 1-2AU of the star. The results support the earlier suggestion that the observed trend of increasing HCN/water emission with disk mass is a chemical fingerprint of planetesimal formation and core accretion in action. In addition to directly constraining the emitting radii of the molecules, the high resolution spectra also help to break degeneracies between temperature and column density in deriving molecular abundances from low resolution mid-infrared spectra. As a result, they can improve our understanding of the extent to which inner disks are chemically active. Contrary to predictions from HCN excitation studies carried out for AS 205 N, the mid-infrared and near-infrared line profiles of HCN are remarkably similar. The discrepancy may indicate that HCN is not abundant beyond a couple of AU or that infrared pumping of HCN does not dominate at these distances.

We present high-spectral resolution observations from 12.96 - 13.33 microns towards Orion IRc2 using the mid-infrared spectrograph, EXES, on SOFIA. These observations probe the physical and chemical conditions of the Orion hot core, which is sampled by a bright, compact, mid-infrared background continuum source in the region, IRc2. All ten of the rovibrational C2H2 transitions expected in our spectral coverage, are detected with high S/N, yielding continuous coverage of the R-branch lines from J=9-8 to J=18-17, including both ortho and para species. Eight of these rovibrational transitions are newly reported detections. The isotopologue, 13CCH2, is clearly detected with high signal-to-noise. This enabled a direct measurement of the 12C/13C isotopic ratio for the Orion hot core of 14 +/- 1 and an estimated maximum value of 21. We also detected several HCN rovibrational lines. The ortho and para C2H2 ladders are clearly separate and tracing two different temperatures, 226 K and 164 K, respectively, with a non-equilibrium ortho to para ratio (OPR) of 1.7 +\- 0.1. Additionally, the ortho and para V_LSR values differ by about 1.8 +/- 0.2 km/s, while, the mean line widths differ by 0.7 +/- 0.2 km/s, suggesting that these species are not uniformly mixed along the line of sight to IRc2. We propose that the abnormally low C2H2 OPR could be a remnant from an earlier, colder phase, before the density enhancement (now the hot core) was impacted by shocks generated from an explosive event 500 yrs ago.

We present high spectral resolution (~3 km/s) observations of the nu_2 ro-vibrational band of H2O in the 6.086--6.135 micron range toward the massive protostar AFGL 2591 using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten absorption features are detected in total, with seven caused by transitions in the nu_2 band of H2O, two by transitions in the first vibrationally excited nu_2 band of H2O, and one by a transition in the nu_2 band of H2{18}O. Among the detected transitions is the nu_2 1(1,1)--0(0,0) line which probes the lowest lying rotational level of para-H2O. The stronger transitions appear to be optically thick, but reach maximum absorption at a depth of about 25%, suggesting that the background source is only partially covered by the absorbing gas, or that the absorption arises within the 6 micron emitting photosphere. Assuming a covering fraction of 25%, the H2O column density and rotational temperature that best fit the observed absorption lines are N(H2O)=(1.3+-0.3)*10^{19} cm^{-2} and T=640+-80 K.

The physical structures of the outer atmospheres of red giants are not known. They are certainly complex and a range of recent observations are showing that we need to embrace to non-classical atmosphere models to interpret these regions. This region's properties is of importance, not the least, for the understanding of the mass-loss mechanism for these stars, which is not still understood. Here, we present observational constraints of the outer regions of red giants, based on mid-IR, high spectral resolution spectra. We also discuss possible non-LTE effects and highlight a new non-LTE code that will be used to analyse the spectra of these atmospheric layers. We conclude by mentioning our new SOFIA/EXES observations of red giants at 6 microns, where the vibration-rotation lines of water vapour can be detected and spectrally resolved for the first time.

We have conducted a near-infrared spectroscopic survey of 47 candidate counterparts to X-ray sources discovered by the Chandra X-ray Observatory near the Galactic Center (GC). Though a significant number of these astrometric matches are likely to be spurious, we sought out spectral characteristics of active stars and interacting binaries, such as hot, massive spectral types or emission lines in order to corroborate the X-ray activity and certify the authenticity of the match. We present three new spectroscopic identifications, including a Be high mass X-ray binary (HMXB) or a gamma Cassiopeiae (Cas) system, a symbiotic X-ray binary and an O-type star of unknown luminosity class. The Be HMXB/gamma Cas system and the symbiotic X-ray binary are the first of their classes to be spectroscopically identified in the GC region.

The Chandra X-ray Observatory has now discovered nearly 10,000 X-ray point sources in the 2 x 0.8 degree region around the Galactic Center (Muno 2009). The sources are likely to be a population of accreting binaries in the Galactic Center, but little else is known of their nature. We obtained JHKs imaging of the 17'x 17' region around Sgr A*, an area containing 4339 of these X-ray sources, with the ISPI camera on the CTIO 4-m telescope. We cross-correlate the Chandra and ISPI catalogs to find potential IR counterparts to the X-ray sources. The extreme IR source crowding in the field means that it is not possible to establish the authenticity of the matches with astrometry and photometry alone. We find 2137 IR/X-ray astrometrically matched sources: statistically we estimate that our catalog contains 289 +/- 13 true matches to soft X-ray sources and 154 +/- 39 matches to hard X-ray sources. However, the fraction of true counterparts to candidate counterparts for hard sources is just 11 %, compared to 60 % for soft sources, making hard source NIR matches particularly challenging for spectroscopic follow-up. We calculate a color-magnitude diagram (CMD) for the matches to hard X-ray sources, and find regions where significant numbers of the IR matches are real. We use their CMD positions to place limits on the absolute Ks band magnitudes of the potential NIR counterparts to hard X-ray sources. We find regions of the counterpart CMD with 9 +/- 3 likely Wolf-Rayet/supergiant binaries (with 4 spectroscopically confirmed in the literature) as well as 44 +/- 13 candidates that could consist of either main sequence high mass X-ray binaries or red giants with an accreting compact companion. (abridged)

We have constructed a thermally compensated field-widened monolithic Michelson interferometer that can be used with a medium-resolution spectrograph to measure precise Doppler radial velocities of stars. Our prototype monolithic fixed-delay interferometer is constructed with off-the-shelf components and assembled using a hydrolysis bonding technique. We installed and tested this interferometer in the Exoplanet Tracker (ET) instrument at the Kitt Peak 2.1m telescope, an instrument built to demonstrate the principles of dispersed fixed delay interferometry. An iodine cell allows the interferometer drift to be accurately calibrated, relaxing the stability requirements on the interferometer itself. When using our monolithic interferometer, the ET instrument has no moving parts (except the iodine cell), greatly simplifying its operation. We demonstrate differential radial velocity precision of a few m s$^{-1}$ on well known radial velocity standards and planet bearing stars when using this interferometer. Such monolithic interferometers will make it possible to build relatively inexpensive instruments that are easy to operate and capable of precision radial velocity measurements. A larger multi-object version of the Exoplanet Tracker will be used to conduct a large scale survey for planetary systems as part of the Sloan Digital Sky Survey III (SDSS III). Variants of the techniques and principles discussed in this paper can be directly applied to build large monolithic interferometers for such applications, enabling the construction of instruments capable of efficiently observing many stars simultaneously at high velocity-precision.

I describe the IR and X-ray observational campaign we have undertaken for the purpose of determining the nature of the faint discrete X-ray source population discovered by Chandra in the Galactic Center (GC). Data obtained for this project includes a deep Chandra survey of the Galactic Bulge; deep, high resolution IR imaging from VLT/ISAAC, CTIO/ISPI, and the UKIDSS Galactic Plane Survey (GPS); and IR spectroscopy from VLT/ISAAC and IRTF/SpeX. By cross-correlating the GC X-ray imaging from Chandra with our IR surveys, we identify candidate counterparts to the X-ray sources via astrometry. Using a detailed IR extinction map, we are deriving magnitudes and colors for all the candidates. Having thus established a target list, we will use the multi-object IR spectrograph FLAMINGOS-2 on Gemini-South to carry out a spectroscopic survey of the candidate counterparts, to search for emission line signatures which are a hallmark of accreting binaries. By determining the nature of these X-ray sources, this FLAMINGOS-2 Galactic Center Survey will have a dramatic impact on our knowledge of the Galactic accreting binary population.

I describe the IR and X-ray campaign we have undertaken to determine the nature of the faint discrete X-ray source population discovered by Chandra in the Galactic Center. These results will provide the input to the FLAMINGOS-2 Galactic Center Survey (F2GCS). With FLAMINGOS-2's multi-object IR spectrograph we will obtain 1000s of IR spectra of candidate X-ray source counterparts, allowing us to efficiently identify the nature of these sources, and thus dramatically increase the number of known X-ray binaries and CVs in the Milky Way.

We demonstrate the ability to measure precise stellar barycentric radial velocities with the dispersed fixed-delay interferometer technique using the Exoplanet Tracker (ET), an instrument primarily designed for precision differential Doppler velocity measurements using this technique. Our barycentric radial velocities, derived from observations taken at the KPNO 2.1 meter telescope, differ from those of Nidever et al. by 0.047 km/s (rms) when simultaneous iodine calibration is used, and by 0.120 km/s (rms) without simultaneous iodine calibration. Our results effectively show that a Michelson interferometer coupled to a spectrograph allows precise measurements of barycentric radial velocities even at a modest spectral resolution of R ~ 5100. A multi-object version of the ET instrument capable of observing ~500 stars per night is being used at the Sloan 2.5 m telescope at Apache Point Observatory for the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS), a wide-field radial velocity survey for extrasolar planets around TYCHO-2 stars in the magnitude range 7.6<V<12. In addition to precise differential velocities, this survey will also yield precise barycentric radial velocities for many thousands of stars using the data analysis techniques reported here. Such a large kinematic survey at high velocity precision will be useful in identifying the signature of accretion events in the Milky Way and understanding local stellar kinematics in addition to discovering exoplanets, brown dwarfs and spectroscopic binaries.

We report the detection of the first extrasolar planet, ET-1 (HD 102195b), using the Exoplanet Tracker (ET), a new generation Doppler instrument. The planet orbits HD 102195, a young star with solar metallicity that may be part of the local association. The planet imparts radial velocity variability to the star with a semiamplitude of $63.4\pm2.0$ m s$^{-1}$ and a period of 4.11 days. The planetary minimum mass ($m \sin i$) is $0.488\pm0.015$ $M_J$.

The Exoplanet Tracker is a prototype of a new type of fibre-fed instrument for performing high precision relative Doppler measurements to detect extra-solar planets. A combination of Michelson interferometer and medium resolution spectrograph, this low-cost instrument facilitates radial velocity measurements with high throughput over a small bandwidth (~ 300 Angstroms), and has the potential to be designed for multi-object operation with moderate bandwidths (~1000 Angstroms). We present the first planet detection with this new type of instrument, a successful confirmation of the well established planetary companion to 51 Peg, showing an rms precision of 11.5m/s over five days. We also show comparison measurements of the radial velocity stable star, Eta Cas, showing an rms precision of 7.9m/s over seven days. These new results are starting to approach the precision levels obtained with traditional radial velocity techniques based on cross-dispersed echelles. We anticipate that this new technique could have an important impact in the search for extra-solar planets.

We propose to use a high throughput and high precision multi-object dispersed fixed-delay interferometer for all sky Doppler surveys for extrasolar planets. This instrument, a combination of a fixed-delay interferometer with a moderate resolution spectrometer,is completely different from current echelle spectrometers. Doppler RV is measured through monitoring interference fringe shifts of stellar absorption lines over a broad band. Coupling this multi-object instrument with a wide field telescope (a few degree, such as Sloan and WIYN) and UV, visible and near-IR detectors will allow to simultaneously obtain hundreds of stellar fringing spectra for searching for planets. The RV survey speed can be increased by more than 2 orders of magnitude over that for the echelles. A prototype dispersed fixed-delay interferometer has been observed at the Hobby-Eberly 9m and Palomar 5m telescopes in 2001 and demonstrated photo noise limited Doppler precision with Aldebaran. Our recent observations at the KPNO 2.1m telescope in 2002 demonstrate a short term Doppler precision of ~ 3 m/s with eta Cas (V = 3.5), a RV stable star and also obtained a RV curve for 51 Peg. (V = 5.5), confirming previous planet detection with an independent RV technique. The total measured detection efficiency including the sky, telescope and fiber transmission losses, the instrument and iodine transmission losses and detector quantum efficiency is 3.4% under 1.5 arcsec seeing conditions, which is comparable to all of the echelle spectrometers for planet detection.