Time domain astronomy and the increasing number of exoplanet candidates call for reliable, robust, and automatic wavelength calibration. We present an algorithm for wavelength calibrating \'echelle spectrographs that uses order-by-order extracted spectra and a list of laboratory wavelengths. Our approach is fully automatic and does not need the pixel locations of certain spectral features with which to anchor the wavelength solution, nor the true order number of each diffraction order. We use spectral features that are duplicated in adjacent orders to establish the scale-invariant component of the wavelength solution. We then match the central wavelengths of spectral features to laboratory wavelengths to establish the scale and higher order components of the wavelength solution. We demonstrate our method on the four spectrographs of Las Cumbres Observatory's Network of Robotic \'Echelle Spectrographs (NRES), on the High Accuracy Radial Velocity Planet Searcher (HARPS) spectrograph, and on synthetic data. We obtain a velocity-equivalent precision of 10 m/s on NRES. We achieve 1 m/s on HARPS, which agrees with the precision reported by the HARPS team. On synthetic data, we achieve the velocity precision set by Gaussian centroiding errors. Our algorithm likely holds for a wide range of spectrographs beyond the five presented here. We provide an open-source Python package, xwavecal (https://github.com/gmbrandt/xwavecal/), which outputs wavelength calibrated spectra as well as the wavelengths of spectral features.

We combine Keck/HIRES radial velocities, imaging with HiCIAO/Subaru and the Hubble Space Telescope, and absolute astrometry from Hipparcos and Gaia to measure a dynamical mass of $70 \pm 5$ $M_{\rm Jup}$ for the brown dwarf companion to Gl 229. Gl 229B was the first imaged brown dwarf to show clear signs of methane in its atmosphere. Cooling models have been used to estimate a mass in the range of 20-55 $M_{\rm Jup}$, much lower than our measured value. We argue that our high dynamical mass is unlikely to be due to perturbations from additional unseen companions or to Gl 229B being itself a binary, and we find no evidence of a previously claimed radial velocity planet around Gl 229A. Future Gaia data releases will confirm the reliability of the absolute astrometry, though the data pass all quality checks in both Hipparcos and Gaia. Our dynamical mass implies a very old age for Gl 229, in some tension with kinematic and activity age indicators, and/or shortcomings in brown dwarf cooling models. Gl 229B joins a small but growing list of T dwarfs with masses approaching the minimum mass for core hydrogen ignition.

We present nebular-phase spectra of the Type Ia supernova (SN Ia) 2016brx, a member of the 1991bg-like subclass that lies at the faint end of the SN Ia luminosity function. Nebular spectra are available for only three other 1991bg-like SNe, and their Co line centers are all within <~ 500 km/s of each other. In contrast, the nebular Co line center of SN 2016brx is blue-shifted by >1500 km/s compared to them and by ~1200 km/s compared to the rest frame. This is a significant shift relative to the narrow nebular line velocity dispersion of <~ 2000 km/s of these SNe. The large range of nebular line shifts implies that the Ni56 in the ejecta of SN 1991bg-like events is off-center by ~1000 km/s rather than universally centrally confined as previously suggested. With the addition of SN 2016brx, the Co nebular line shapes of 1991bg-like objects appear to connect with the brighter SNe Ia that show double-peak profiles, hinting at a continuous distribution of line profiles among SNe Ia. One class of models to produce both off-center and bi-modal Ni56 distributions is collisions of white dwarfs with unequal and equal masses.