We present high resolution near-infrared spectro-polarimetric observations using the SPIRou instrument at CFHT during a transit of the recently detected young planet AU Mic b, with supporting spectroscopic data from iSHELL at IRTF. We detect Zeeman signatures in the Stokes V profiles, and measure a mean longitudinal magnetic field of $\overline{B}_\ell=46.3\pm0.7$~G. Rotationally modulated magnetic spots cause long-term variations of the field with a slope of $d{B_\ell}/dt=-108.7\pm7.7$~G/d. We apply the cross-correlation technique to measure line profiles and obtain radial velocities through CCF template matching. We find an empirical linear relationship between radial velocity and $B_\ell$, which allows us to calibrate out the radial velocity variations which stellar activity induces through rotational modulation of spots. We model the corrected radial velocities for the classical Rossiter-McLaughlin effect, using MCMC to sample the posterior distribution of the model parameters. This analysis shows that the orbit of AU Mic b is prograde and aligned with the stellar rotation axis with a sky-projected spin-orbit obliquity of $\lambda=0^{\circ}\pm5^{\circ}$. The aligned orbit of AU Mic b indicates that it formed in the protoplanetary disk that evolved to the current debris disk around AU Mic.

We describe a flexible data reduction package for high resolution cross-dispersed echelle data. This open-source package is developed in Python and includes optional GUIs for most of the steps. It does not require any pre-knowledge about the form or position of the echelle-orders. It has been tested on cross-dispersed echelle spectrographs between 13k and 115k resolution (bifurcated fiber-fed spectrogaph ESO-HARPS and single fiber-fed spectrograph TNT-MRES). HiFLEx can be used to determine radial velocities and is designed to use the TERRA package but can also control the radial velocity packages such as CERES and SERVAL to perform the radial velocity analysis. Tests on HARPS data indicates radial velocities results within 3m/s of the literature pipelines without any fine tuning of extraction parameters.

Our understanding of the brown dwarf population in star forming regions is dependent on knowing distances and proper motions, and therefore will be improved through the Gaia space mission. In this paper, we select new samples of very low mass objects (VLMOs) in Upper Scorpius using UKIDSS colors and optimised proper motions calculated using Gaia DR1. The scatter in proper motions from VLMOs in Upper Scorpius is now (for the first time) dominated by the kinematic spread of the region itself, not by the positional uncertainties. With age and mass estimates updated using Gaia parallaxes for early type stars in the same region, we determine masses for all VLMOs. Our final most complete sample includes 453 VLMOs of which $\sim$125 are expected to be brown dwarfs. The cleanest sample is comprised of 131 VLMOs, with $\sim$105 brown dwarfs. We also compile a joint sample from the literature which includes 415 VLMOs, out of which 152 are likely brown dwarfs. The disc fraction among low-mass brown dwarfs ($M<0.05 M_{\odot}$) is substantially higher than in more massive objects, indicating that discs around low-mass brown dwarfs survive longer than in low-mass stars overall. The mass function for $0.01<M<0.1$ $M_{\odot}$ is consistent with the Kroupa IMF. We investigate the possibility that some 'proper motion outliers' have undergone a dynamical ejection early in their evolution. Our analysis shows that the color-magnitude cuts used when selecting samples introduce strong bias into the population statistics due to varying level of contamination and completeness.

We explore the potential of \emph{Gaia} for the field of benchmark ultracool/brown dwarf companions, and present the results of an initial search for metal-rich/metal-poor systems. A simulated population of resolved ultracool dwarf companions to \emph{Gaia} primary stars is generated and assessed. Of order $\sim$24,000 companions should be identifiable outside of the Galactic plane ($|b| > 10\,$deg) with large-scale ground- and space-based surveys including late M, L, T, and Y types. Our simulated companion parameter space covers $0.02 \le M/M_{\odot} \le 0.1$, $0.1 \le {\rm age/Gyr} \le 14$, and $-2.5 \le {\rm [Fe/H]} \le 0.5$, with systems required to have a false alarm probability $<10^{-4}$, based on projected separation and expected constraints on common-distance, common-proper motion, and/or common-radial velocity. Within this bulk population we identify smaller target subsets of rarer systems whose collective properties still span the full parameter space of the population, as well as systems containing primary stars that are good age calibrators. Our simulation analysis leads to a series of recommendations for candidate selection and observational follow-up that could identify $\sim$500 diverse \emph{Gaia} benchmarks. As a test of the veracity of our methodology and simulations, our initial search uses UKIDSS and SDSS to select secondaries, with the parameters of primaries taken from Tycho-2, RAVE, LAMOST and TGAS. We identify and follow-up 13 new benchmarks. These include M8-L2 companions, with metallicity constraints ranging in quality, but robust in the range $-0.39 \le {\rm [Fe/H]} \le +0.36$, and with projected physical separation in the range $0.6\,<\,s/{\rm kau}\,<76$. Going forward, \emph{Gaia} offers a very high yield of benchmark systems, from which diverse sub-samples may be able to calibrate a range of foundational ultracool/sub-stellar theory and observation.

We develop a method to identify the spectroscopic signature of unresolved L-dwarf ultracool companions, which compares the spectra of candidates and their associated control stars using spectral ratio differences and residual spectra. We present SpeX prism-mode spectra (0.7-2.5 micron) for a pilot sample of 111 mid M dwarfs, including 28 that were previously identified as candidates for unresolved ultracool companionship (a sub-sample from Cook et al. 2016; paper 1) and 83 single M dwarfs that were optically colour-similar to these candidates (which we use as `control stars'). We identify four candidates with evidence for near-infrared excess. One of these (WISE J100202.50+074136.3) shows strong evidence for an unresolved L dwarf companion in both its spectral ratio difference and its residual spectra, two most likely have a different source for the near-infrared excess, and the other may be due to spectral noise. We also establish expectations for a null result (i.e. by searching for companionship signatures around the M dwarf control stars), as well as determining the expected outcome for ubiquitous companionship (as a means of comparison with our actual results), using artificially generated unresolved M+L dwarf spectra. The results of these analyses are compared to those for the candidate sample, and reasonable consistency is found. With a full follow-up programme of our candidates sample from Cook et al., we might expect to confirm up to 40 such companions in the future, adding extensively to the known desert population of M3-M5 dwarfs.

The task of flux calibration for LAMOST (Large sky Area Multi-Object Spectroscopic Telescope) spectra is difficult due to many factors. For example, the lack of standard stars, flat fielding for large field of view, and variation of reddening between different stars especially at low galactic latitudes etc. Poor selection, bad spectral quality, or extinction uncertainty of standard stars not only might induce errors to the calculated spectral response curve (SRC), but also might lead to failures in producing final 1D spectra. In this paper, we inspected spectra with Galactic latitude |b|>=60 degree and reliable stellar parameters, determined through the LAMOST Stellar Parameter Pipeline (LASP), to study the stability of the spectrograph. To guarantee the selected stars had been observed by each fiber, we selected 37,931 high quality exposures of 29,000 stars from LAMOST DR2, and more than 7 exposures for each fiber. We calculated the SRCs for each fiber for each exposure, and calculated the statistics of SRCs for spectrographs with both the fiber variations and time variations. The result shows that the average response curve of each spectrograph (henceforth ASPSRC) is relatively stable with statistical errors <= 10%. From the comparison between each ASPSRC and the SRCs for the same spectrograph obtained by 2D pipeline, we find that the ASPSRCs are good enough to use for the calibration. The ASPSRCs have been applied to spectra which were abandoned by LAMOST 2D pipeline due to the lack of standard stars, increasing the number of LAMOST spectra by 52,181 in DR2. Comparing those same targets with SDSS, the relative flux differences between SDSS spectra and that of LAMOST spectra with the ASPSRC method are less than 10%, which underlines that the ASPSRC method is feasible for LAMOST flux calibration.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

During the planet formation process, billions of comets are created and ejected into interstellar space. The detection and characterization of such interstellar comets (also known as extra-solar planetesimals or extra-solar comets) would give us in situ information about the efficiency and properties of planet formation throughout the galaxy. However, no interstellar comets have ever been detected, despite the fact that their hyperbolic orbits would make them readily identifiable as unrelated to the solar system. Moro-Mart\'in et al. 2009 have made a detailed and reasonable estimate of the properties of the interstellar comet population. We extend their estimates of detectability with a numerical model that allows us to consider "close" interstellar comets, e.g., those that come within the orbit of Jupiter. We include several constraints on a "detectable" object that allow for realistic estimates of the frequency of detections expected from the Large Synoptic Survey Telescope (LSST) and other surveys. The influence of several of the assumed model parameters on the frequency of detections is explored in detail. Based on the expectation from Moro-Mart\'in et al. 2009, we expect that LSST will detect 0.001-10 interstellar comets during its nominal 10-year lifetime, with most of the uncertainty from the unknown number density of small (nuclei of $\sim$0.1-1 km) interstellar comets. Using simulated LSST-like astrometric data, we study the problem of orbit determination for these bodies, finding that LSST could identify their orbits as hyperbolic and determine an ephemeris sufficiently accurate for follow-up in about 4-7 days. We give the hyperbolic orbital parameters of the most detectable interstellar comets. Taking the results into consideration, we give recommendations to future searches for interstellar comets.

Locating ultra-cool companions to M dwarfs is important for constraining low-mass formation models, the measurement of sub-stellar dynamical masses and radii, and for testing ultra-cool evolutionary models. We present an optimised method for identifying M dwarfs which may have unresolved ultra-cool companions. We construct a catalogue of 440,694 candidates, from WISE, 2MASS and SDSS, based on optical and near-infrared colours and reduced proper motion. With strict reddening, photometric and quality constraints we isolate a sub-sample of 36,898 M dwarfs and search for possible mid-infrared M dwarf + ultra-cool dwarf candidates by comparing M dwarfs which have similar optical/near-infrared colours (chosen for their sensitivity to effective temperature and metallicity). We present 1,082 M dwarf + ultra-cool dwarf candidates for follow-up. Using simulated ultra-cool dwarf companions to M dwarfs, we estimate that the occurrence of unresolved ultra-cool companions amongst our M dwarf + ultra-cool dwarf candidates should be at least four times the average for our full M dwarf catalogue. We discuss possible contamination and bias and predict yields of candidates based on our simulations.

We present a new Y dwarf, WISE J030449.03-270508.3, confirmed from a candidate sample designed to pick out low temperature objects from the WISE database. The new object is typed Y0pec following a visual comparison with spectral standards, and lies at a likely distance of 10-17 pc. Its tangential velocity suggests thin disk membership, but it shows some spectral characteristics that suggest it may be metal-poor and/or older than previously identified Y0 dwarfs. Based on trends seen for warmer late type T dwarfs, the Y-band flux peak morphology is indicative of sub-solar metallicity, and the enhanced red wing of the J-band flux peak offers evidence for high gravity and/or low metallicity (with associated model trends suggesting an age closer to ~10 Gyr and mass in the range 0.02-0.03 Mo). This object may thus be extending the population parameter-space of the known Y0 dwarfs.