Astrophysics lies at the crossroads of big datasets (such as the Large Synoptic Survey Telescope and Gaia), open source software to visualize and interpret high dimensional datasets (such as Glue, WorldWide Telescope, and OpenSpace), and uniquely skilled software engineers who bridge data science and research fields. At the same time, more than 4,000 planetariums across the globe immerse millions of visitors in scientific data. We have identified the potential for critical synergy across data, software, hardware, locations, and content that -- if prioritized over the next decade -- will drive discovery in astronomical research. Planetariums can and should be used for the advancement of scientific research. Current facilities such as the Hayden Planetarium in New York City, Adler Planetarium in Chicago, Morrison Planetarium in San Francisco, the Iziko Planetarium and Digital Dome Research Consortium in Cape Town, and Visualization Center C in Norrkoping are already developing software which ingests catalogs of astronomical and multi-disciplinary data critical for exploration research primarily for the purpose of creating scientific storylines for the general public. We propose a transformative model whereby scientists become the audience and explorers in planetariums, utilizing software for their own investigative purposes. In this manner, research benefits from the authentic and unique experience of data immersion contained in an environment bathed in context and equipped for collaboration. Consequently, in this white paper we argue that over the next decade the research astronomy community should partner with planetariums to create visualization-based research opportunities for the field. Realizing this vision will require new investments in software and human capital.

The physical properties of brown dwarf companions found to orbit nearby, solar-type stars can be benchmarked against independent measures of their mass, age, chemical composition, and other parameters, offering insights into the evolution of substellar objects. The TRENDS high-contrast imaging survey has recently discovered a (mass/age/metallicity) benchmark brown dwarf orbiting the nearby (d=18.69+/-0.19 pc), G8V/K0V star HD 4747. We have acquired follow-up spectroscopic measurements of HD 4747 B using the Gemini Planet Imager to study its spectral type, effective temperature, surface gravity, and cloud properties. Observations obtained in the H-band and K1-band recover the companion and reveal that it is near the L/T transition (T1+/-2). Fitting atmospheric models to the companion spectrum, we find strong evidence for the presence of clouds. However, spectral models cannot satisfactorily fit the complete data set: while the shape of the spectrum can be well-matched in individual filters, a joint fit across the full passband results in discrepancies that are a consequence of the inherent color of the brown dwarf. We also find a $2\sigma$ tension in the companion mass, age, and surface gravity when comparing to evolutionary models. These results highlight the importance of using benchmark objects to study "secondary effects" such as metallicity, non-equilibrium chemistry, cloud parameters, electron conduction, non-adiabatic cooling, and other subtleties affecting emergent spectra. As a new L/T transition benchmark, HD 4747 B warrants further investigation into the modeling of cloud physics using higher resolution spectroscopy across a broader range of wavelengths, polarimetric observations, and continued Doppler radial velocity and astrometric monitoring.

The nearby Sun-like star GJ 758 hosts a cold substellar companion, GJ 758 B, at a projected separation of $\lesssim$30 AU, previously detected in high-contrast multi-band photometric observations. In order to better constrain the companion's physical characteristics, we acquired the first low-resolution ($R \sim 50$) near-infrared spectrum of it using the high-contrast hyperspectral imaging instrument Project 1640 on Palomar Observatory's 5-m Hale telescope. We obtained simultaneous images in 32 wavelength channels covering the $Y$, $J$, and $H$ bands ($\sim$952-1770 nm), and used data processing techniques based on principal component analysis to efficiently subtract chromatic background speckle-noise. GJ 758 B was detected in four epochs during 2013 and 2014. Basic astrometric measurements confirm its apparent northwest trajectory relative to the primary star, with no clear signs of orbital curvature. Spectra of SpeX/IRTF observed T dwarfs were compared to the combined spectrum of GJ 758 B, with ${\chi}^2$ minimization suggesting a best fit for spectral type T7.0$\pm$1.0, but with a shallow minimum over T5-T8. Fitting of synthetic spectra from the BT-Settl13 model atmospheres gives an effective temperature $T_{\text{eff}}=741 \pm 25$ K and surface gravity $\log g = 4.3 \pm 0.5$ dex (cgs). Our derived best-fit spectral type and effective temperature from modeling of the low-resolution spectrum suggest a slightly earlier and hotter companion than previous findings from photometric data, but do not rule out current results, and confirm GJ 758 B as one of the coolest sub-stellar companions to a Sun-like star to date.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

We have used the Swedish ESO Submillimetre Telescope (SEST) to search for new class II methanol maser transitions towards the southern source G345.01+1.79. Over a period of 5 days we observed 11 known or predicted class II methanol maser transitions. Emission with the narrow line width and characteristic velocity of class II methanol masers (in this source) was detected in 8 of these transitions, two of which have not previously been reported as masers. The new class II methanol maser transitions are the 13(-3)-12(-4)E transition at 104.1 GHz and the 5(1)-4(2)E transition at 216.9 GHz. Both of these are from transition series for which there are no previous known class II methanol maser transitions. This takes the total number of known class II methanol maser series to 10, and the total number of transitions (or transition groups) to 18. The observed 104.1 GHz maser suggests the presence of two or more regions of masing gas with similar line of sight velocities, but quite different physical conditions. Although these newly discovered transitions are likely to be relatively rare, where they are observed combined studies using the Australia Telescope Compact Array and the Atacama Large Millimeter Array offer the prospect to be able to undertake multi-transition methanol maser studies with unprecedented detail.

Interstellar masers are unique probes of the environments in which they arise. In studies of high-mass star formation their primary function has been as signposts of these regions and they have been used as probes of the kinematics and physical conditions in only a few sources. With a few notable exceptions, we know relatively little about the evolutionary phase the different maser species trace, nor their location with respect to other star formation tracers. While detailed studies of a small number of maser regions can reveal much about them, other information can only be obtained through large, systematic searches. In particular, such surveys are vital in efforts to determine an evolutionary sequence for the common maser species, and there is growing evidence that methanol masers may trace an earlier phase than the other common maser species of OH and water.

We have used the Mopra Telescope to search for glycine and the simple chiral molecule propylene oxide in the Sgr B2 (LMH) and Orion KL, in the 3-mm band. We have not detected either species, but have been able to put sensitive upper limits on the abundances of both molecules. The 3-sigma upper limits derived for glycine conformer I are 3.7 x 10^{14} cm^{-2} in both Orion-KL and Sgr B2 (LMH), comparable to the reported detections of conformer I by Kuan et al. However, as our values are 3-sigma upper limits rather than detections we conclude that this weighs against confirming the detection of Kuan et al. We find upper limits for the glycine II column density of 7.7 x 10^{12} cm^{-2} in both Orion-KL and Sgr B2 (LMH), in agreement with the results of Combes et al. The results presented here show that glycine conformer II is not present in the extended gas at the levels detected by Kuan et al. for conformer I. Our ATCA results (Jones et al.) have ruled out the detection of glycine (both conformers I and II) in the compact hot core of the LMH at the levels reported, so we conclude that it is unlikely that Kuan et al. have detected glycine in either Sgr B2 or Orion-KL. We find upper limits for propylene oxide abundance of 3.0 x 10^{14} cm^{-2} in Orion-KL and 6.7 x 10^{14} cm^{-2} in Sgr B2 (LMH). We have detected fourteen features in Sgr B2 and four features in Orion-KL which have not previously been reported in the ISM, but have not be able to plausibly assign these transitions to any carrier.

We have used the Australia Telescope Compact Array to conduct a search for the simplest amino acid, glycine (conformers I and II), and the simple chiral molecule propylene oxide at 3-mm in the Sgr B2 LMH. We searched 15 portions of spectrum between 85 and 91 GHz, each of 64 MHz bandwidth, and detected 58 emission features and 21 absorption features, giving a line density of 75 emission lines and 25 absorption lines per GHz stronger than the 5 sigma level of 110 mJy. Of these, 19 are transitions previously detected in the interstellar medium, and we have made tentative assignments of a further 23 features to molecular transitions. However, as many of these involve molecules not previously detected in the ISM, these assignments cannot be regarded with confidence. Given the median line width of 6.5 km/s in Sgr B2 LMH, we find that the spectra have reached a level where there is line confusion, with about 1/5 of the band being covered with lines. Although we did not confidently detect either glycine or propylene oxide, we can set 3 sigma upper limits for most transitions searched. We also show that if glycine is present in the Sgr B2 LMH at the level of N = 4 x 10^{14} cm^{-2} found by Kuan et al. (2003) in their reported detection of glycine, it should have been easily detected with the ATCA synthesized beam size of 17.0 x 3.4 arcsec^{2}, if it were confined to the scale of the LMH continuum source (< 5 arcsec). This thus puts a strong upper limit on any small-scale glycine emission in Sgr B2, for both of conformers I and II.

The class II masers of methanol are associated with the early stages of formation of high-mass stars. Modelling of these dense, dusty environments has demonstrated that pumping by infrared radiation can account for the observed masers. Collisions with other molecules in the ambient gas also play a significant role, but have not been well modelled in the past. Here we examine the effects on the maser models of newly available collision rate coefficients for methanol. The new collision data does not alter which transitions become masers in the models, but does influence their brightness and the conditions under which they switch on and off. At gas temperatures above 100 K the effects are broadly consistent with a reduction in the overall collision cross-section. This means, for example, that a slightly higher gas density than identified previously can account for most of the observed masers in W3(OH). We have also examined the effects of including more excited state energy levels in the models, and find that these play a significant role only at dust temperatures above 300 K. An updated list of class II methanol maser candidates is presented.

In studies of the environment of massive young stellar objects, recent progress in both observations and theory allows a unified treatment of data for maser and 'quasi-thermal' lines. Interferometric maser images provide information on the distribution and kinematics of masing gas on small spatial scales. Observations of multiple masing transitions provide constraints on the physical parameters. Interferometric data on 'quasi-thermal' molecular lines permits an investigation of the overall distribution and kinematics of the molecular gas in the vicinity of young stellar objects, including those which are deeply embedded. Using multiple transitions of different molecules, one can obtain good constraints on the physical and chemical parameters. Combining these data enables the construction of unified models, which take into account spatial scales differing by orders of magnitude. Here we present such a combined analysis of the environment around the ultracompact HII region in W3(OH). This includes the structure of the methanol masing region, physical structure of the near vicinity of W3(OH), detection of new masers in the large-scale shock front and embedded sources in the vicinity of the TW young stellar object.

We have used the NASA Tidbinbilla 70-m antenna to search for emission from the 2(1)-3(0)E 19.9-GHz transition of methanol. The search was targeted towards 22 star formation regions which exhibit maser emission in the 107.0-GHz 3(1)-4(0)A+ methanol transition, as well as in the 6.6-GHz 5(1)-6(0)A+ transition characteristic of class II methanol maser sources. A total of 7 sources were detected in the 19.9-GHz transition, 6 of these being new detections. Many of the new detections are weak (< 0.5 Jy), however, they appear to be weak masers rather than thermal or quasi-thermal emission. We find a strong correlation between sources which exhibit 19.9-GHz methanol masers and those which both have the class II methanol masers projected against radio continuum emission and have associated 6035-MHz OH masers. This suggests that the 19.9-GHz methanol masers arise in very specific physical conditions, probably associated with a particular evolutionary phase. In the model of (Cragg, Sobolev & Godfrey 2002) these observations are consistent with gas temperatures of 50 K, dust temperatures of 150-200 K and gas densities of 10^6.5-10^7.5 cm^-3.

Methanol masers and molecular shock tracers were observed in the W3(OH)/W3(H$_2$O) region with the BIMA array and the Onsala 20m radiotelescope. Characteristics of the outflows in the region are discussed. A model of the W3(OH) methanol maser formation region is constructed.

In the early days of methanol maser discoveries the 9(2)-10(1) A+ transition at 23.1 GHz was found to exhibit maser characteristics in the northern star-forming region W3(OH), and probable maser emission in two other sources. Attention subsequently turned to the 6.6-GHz 5(1)-6(0) A+ methanol maser transition, which has proved a valuable tracer of early high-mass star formation. We have undertaken a new search for 23.1-GHz methanol masers in 50 southern star formation regions using the Parkes radiotelescope. The target sources all exhibit class II methanol maser emission at 6.6 GHz, with 20 sources also displaying maser features in the 107.0-GHz 3(1)-4(0) A+ methanol line. Strong emission at 23.1 GHz in NGC 6334F was confirmed, but no emission was detected in the remaining sources. Thus the 23.1-GHz methanol masers are rare. A maser model in which methanol molecules are pumped to the second torsionally excited state by radiation from warm dust can account for class II maser activity in all the transitions in which it is observed. According to this model the 23.1-GHz maser is favoured by conditions representing low gas temperature, high external dust temperature, low gas density, and high column density of methanol; the scarcity of this maser indicates that such combinations of conditions are uncommon. We have undertaken new model calculations to examine the range of parameters compatible with the upper limits on 23.1-GHz emission from our survey.

We have used the Australia Telescope National Facility Mopra 22m millimetre telescope to search for emission from the 85.5-GHz and 86.6-GHz transitions of methanol. The search was targeted towards 22 star formation regions which exhibit maser emission in the 107.0-GHz methanol transition, as well as in the 6.6-GHz transition characteristic of class II methanol maser sources. A total of 22 regions were searched at 85.5 GHz resulting in 5 detections, of which 1 appears to be a newly discovered maser. For the 86.6-GHz transition observations were made of 18 regions which yielded 2 detections, but no new maser sources. This search demonstrates that emission from the 85.5- and 86.6-GHz transitions is rare. Detection of maser emission from either of these transitions therefore indicates the presence of special conditions, different from those in the majority of methanol maser sources. We have observed temporal variability in the 86.6-GHz emission towards 345.010+1.792, which along with the very narrow line width, confirms that the emission is a maser in this source. We have combined our current observations with published data for the 6.6-, 12.1-, 85.5-, 86.6-, 107.0-, 108.8- and 156.6-GHz transitions for comparison with the maser model of Sobolev & Deguchi (1994). This has allowed us to estimate the likely ranges of dust temperature, gas density, and methanol column density, both for typical methanol maser sources and for those sources which also show 107.0-GHz emission.

We have used the ATNF Mopra antenna and the SEST antenna to search in the directions of several class II methanol maser sources for emission from six methanol transitions in the frequency range 85-115 GHz. The transitions were selected from excitation studies as potential maser candidates. Methanol emission at one or more frequencies was detected from five of the maser sources, as well as from Orion KL. Although the lines are weak, we find evidence of maser origin for three new lines in G345.01+1.79, and possibly one new line in G9.62+0.20. The observations, together with published maser observations at other frequencies, are compared with methanol maser modelling for G345.01+1.79 and NGC6334F. We find that the majority of observations in both sources are consistent with a warm dust (175 K) pumping model at hydrogen density ~10^6 cm^-3 and methanol column density ~5 x 10^17 cm^-2. The substantial differences between the maser spectra in the two sources can be attributed to the geometry of the maser region.