We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic labels from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise -- and by construction accurate -- prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

For the past 150 years, the prevailing view of the local Interstellar Medium (ISM) was based on a peculiarity known as the Gould's Belt, an expanding ring of young stars, gas, and dust, tilted about 20$^\circ$ to the Galactic plane. Still, the physical relation between local gas clouds has remained practically unknown because the distance accuracy to clouds is of the same order or larger than their sizes. With the advent of large photometric surveys and the Gaia satellite astrometric survey this situation has changed. Here we report the 3-D structure of all local cloud complexes. We find a narrow and coherent 2.7 kpc arrangement of dense gas in the Solar neighborhood that contains many of the clouds thought to be associated with the Gould Belt. This finding is inconsistent with the notion that these clouds are part of a ring, disputing the Gould Belt model. The new structure comprises the majority of nearby star-forming regions, has an aspect ratio of about 1:20, and contains about 3 million solar masses of gas. Remarkably, the new structure appears to be undulating and its 3-D distribution is well described by a damped sinusoidal wave on the plane of the Milky Way, with an average period of about 2 kpc and a maximum amplitude of about 160 pc. Our results represent a first step in the revision of the local gas distribution and Galactic structure and offer a new, broader context to studies on the transformation of molecular gas into stars.

Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method which combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of $\sim5\%$. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, vols I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with $\lesssim10\%$ scatter for clouds across our entire distance spectrum (150 pc $-$ 2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.

Using a population of large-scale filaments extracted from an AREPO simulation of a Milky Way-like galaxy, we seek to understand the extent to which observed large-scale filament properties (with lengths $\gtrsim 100$ pc) can be explained by galactic dynamics alone. From an observer's perspective in the disk of the galaxy, we identify filaments forming purely due to galactic dynamics, without the effects of feedback or local self-gravity. We find that large-scale Galactic filaments are intrinsically rare, and we estimate that at maximum approximately one filament per $\rm kpc^{2}$ should be identified in projection, when viewed from the direction of our Sun in the Milky Way. In this idealized scenario, we find filaments in both the arm and interarm regions, and hypothesize that the former may be due to gas compression in the spiral-potential wells, with the latter due to differential rotation. Using the same analysis pipeline applied previously to observations, we analyze the physical properties of large-scale Galactic filaments, and quantify their sensitivity to projection effects and galactic environment (i.e. whether they lie in the arm or interarm regions). We find that observed "Giant Molecular Filaments" are consistent with being non-self-gravitating structures dominated by galactic dynamics. Straighter, narrower, and denser "Bone-like" filaments, like the paradigmatic Nessie filament, have similar column densities, velocity gradients, and Galactic plane heights ($z\approx$ 0 pc) to those in our simple model, but additional physical effects (such as feedback and self-gravity) must be invoked to explain their lengths and widths.

We present a new three-dimensional map of dust reddening, based on Gaia parallaxes and stellar photometry from Pan-STARRS 1 and 2MASS. This map covers the sky north of a declination of -30 degrees, out to a distance of several kiloparsecs. This new map contains three major improvements over our previous work. First, the inclusion of Gaia parallaxes dramatically improves distance estimates to nearby stars. Second, we incorporate a spatial prior that correlates the dust density across nearby sightlines. This produces a smoother map, with more isotropic clouds and smaller distance uncertainties, particularly to clouds within the nearest kiloparsec. Third, we infer the dust density with a distance resolution that is four times finer than in our previous work, to accommodate the improvements in signal-to-noise enabled by the other improvements. As part of this work, we infer the distances, reddenings and types of 799 million stars. We obtain typical reddening uncertainties that are ~30% smaller than those reported in the Gaia DR2 catalog, reflecting the greater number of photometric passbands that enter into our analysis. Our 3D dust map can be accessed at https://doi.org/10.7910/DVN/2EJ9TX or through the Python package "dustmaps," and can be queried interactively at http://argonaut.skymaps.info. Our catalog of stellar parameters can be accessed at https://doi.org/10.7910/DVN/AV9GXO.

We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. (2014). First, we infer the distance and extinction to stars along sightlines towards the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance-extinction estimates as being caused by a dust screen with a 2-D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. (2014) sightlines towards the Dame et al. (2001) and Magnani et al. (1985) clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is approximately 5-6%, though this can be higher for clouds at farther distances, due to the limitations of our single-cloud model.

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and $\rm ^{12} CO$ spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and 2MASS photometry, based on a technique presented in Green et al. 2014; Green et al. 2015 and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine the Green et al. technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight towards these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices towards each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC348, B1, NGC1333, L1448, and L1451 star-forming regions and find that individual clouds are located between $\approx 275-300$ pc, with typical combined uncertainties of $\approx 5\%$. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of $294\pm 17$ pc, about 60 pc higher than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond.

We present RadFil, a publicly available Python package that gives users full control over how to build and fit radial density profiles for interstellar filaments. RadFil builds filament profiles by taking radial cuts across the spine of a filament, thereby preserving the radial structure of the filament across its entire length. Pre-existing spines can be inputted directly into RadFil, or can be computed using the FilFinder package as part of the RadFil workflow. We provide Gaussian and Plummer built-in fitting functions, in addition to a background subtraction estimator, which can be fit to the entire ensemble of radial cuts or an average radial profile for the filament. Users can tweak parameters like the radial cut sampling interval, the background subtraction estimation radii, and the Gaussian/Plummer fitting radii. As a result, RadFil can provide treatment of how the resulting filament properties rely on systematics in the building and fitting process. We walk through the typical RadFil workflow and compare our results to those from an independent radial density profile code obtained using the same data; we find that our results are entirely consistent. RadFil is open source and available on GitHub. We also provide a complete working tutorial of the code available as a Jupyter notebook which users can download and run themselves.

The characterization of our Galaxy's longest filamentary gas features has been the subject of several studies in recent years, producing not only a sizeable sample of large-scale filaments, but also confusion as to whether all these features (e.g. "Bones", "Giant Molecular Filaments") are the same. They are not. We undertake the first standardized analysis of the physical properties ($\rm{H_2}\;$column densities, dust temperatures, morphologies, radial column density profiles) and kinematics of large-scale filaments in the literature. We expand and improve upon prior analyses by using the same data sets, techniques, and spiral arm models to disentangle the filaments' inherent properties from selection criteria and methodology. Our results suggest that the myriad filament finding techniques are uncovering different physical structures, with length ($11-269\,\rm\,pc$), width ($1-40\,\rm\,pc$), mass ($\rm3\times10^3\; M_\odot-1.1\times10^{6}\;M_\odot$), aspect ratio (3:1-117:1), and high column density fraction ($0.2-100\%$) varying by over an order of magnitude across the sample of 45 filaments. We develop a radial profile fitting code, $\texttt{RadFil}$, which is publicly available. We also perform a $\textit{position-position-velocity}$ ($\textit{p-p-v}$) analysis on a subsample and find that while $60-70\%$ lie spatially in the plane of the Galaxy, only $30-45\%$ concurrently exhibit spatial $\textit{and}$ kinematic proximity to spiral arms. In a parameter space defined by aspect ratio, dust temperature, and column density, we broadly distinguish three filament categories, which could indicate different formation mechanisms or histories. Highly elongated "Bone-like" filaments show the most potential for tracing gross spiral structure (e.g. arms, spurs), while other categories could be large concentrations of molecular gas (GMCs, core complexes).

Compact groups (CGs) provide an environment in which interactions between galaxies and with the intra-group medium enable and accelerate galaxy transitions from actively star forming to quiescent. Galaxies in transition from active to quiescent can be selected, by their infrared (IR) colors, as canyon or infrared transition zone (IRTZ) galaxies. We used a sample of CG galaxies with IR data from the Wide Field Infrared Survey Explorer (WISE) allowing us to calculate the stellar mass and star formation rate (SFR) for each galaxy. Furthermore, we present new CO(1-0) data for 27 galaxies and collect data from the literature to calculate the molecular gas mass for a total sample of 130 galaxies. This data set allows us to study the difference in the molecular gas fraction (Mmol/Mstar) and star formation efficiency (SFE=SFR/Mmol) between active, quiescent, and transitioning (i.e., canyon and IRTZ) galaxies. We find that transitioning galaxies have a mean molecular gas fraction and a mean SFE that are significantly lower than those of actively star-forming galaxies. The molecular gas fraction is higher than that of quiescent galaxies, whereas the SFE is similar. These results indicate that the transition from actively star-forming to quiescent in CG galaxies goes along with a loss of molecular gas, possibly due to tidal forces exerted from the neighboring galaxies or a decrease in the gas density. In addition, the remaining molecular gas loses its ability to form stars efficiently, possibly owing to turbulence perturbing the gas, as seen in other, well-studied examples such as Stephan's Quintet and HCG~57. Thus, the amount and properties of molecular gas play a crucial role in the environmentally driven transition of galaxies from actively star forming to quiescent.

Recently, Goodman et al. (2014) argued that the very long, very thin infrared dark cloud "Nessie" lies directly in the Galactic mid-plane and runs along the Scutum-Centaurus arm in position-position-velocity ($p-p-v$) space as traced by lower density $\rm {CO}$ and higher density ${\rm NH}_3$ gas. Nessie was presented as the first "bone" of the Milky Way, an extraordinarily long, thin, high-contrast filament that can be used to map our Galaxy's "skeleton." Here, we present evidence for additional bones in the Milky Way Galaxy, arguing that Nessie is not a curiosity but one of several filaments that could potentially trace Galactic structure. Our ten bone candidates are all long, filamentary, mid-infrared extinction features which lie parallel to, and no more than 20 pc from, the physical Galactic mid-plane. We use $\rm {CO}$, ${\rm N}_2{\rm H}^+$, $\rm {HCO}^+$, and ${\rm NH}_3$ radial velocity data to establish the three-dimensional location of the candidates in ${\it p-p-v}$ space. Of the ten candidates, six also: have a projected aspect ratio of $\geqq$50:1; run along, or extremely close to, the Scutum-Centaurus arm in ${\it p-p-v}$ space; and exhibit no abrupt shifts in velocity. The evidence presented here suggests that these candidates are marking the locations of significant spiral features, with the bone called filament 5 ("BC_18.88-0.09") being a close analog to Nessie in the Northern Sky. As molecular spectral-line and extinction maps cover more of the sky at increasing resolution and sensitivity, it should be possible to find more bones in future studies.

Compact groups provide an environment to study the growth of galaxies amid multiple prolonged interactions. With their dense galaxy concentrations and relatively low velocity dispersions, compact groups mimic the conditions of hierarchical galaxy assembly. Compact group galaxies are known to show a bimodality in $Spitzer$ IRAC infrared colorspace: galaxies are preferentially either quiescent with low specific star formation rates, or are prolifically forming stars---galaxies with moderate levels of specific star formation are rare. Previous $Spitzer$ IRAC studies identifying this "canyon" have been limited by small number statistics. We utilize whole-sky WISE data to study 163 compact groups, thereby tripling our previous sample and including more galaxies with intermediate mid-IR colors indicative of moderate specific star formation rates (SSFRs). We define a distinct WISE mid-IR color-space ($\log[{\frac{\rm f_{12}}{\rm f_{4.6}}}]$ vs. $\log[{\frac{\rm f_{22}}{\rm f_{3.4}}}]$) that we use to identify canyon galaxies from the larger sample. We confirm that compact group galaxies show a bimodal distribution in the mid-infrared and identify 37 canyon galaxies with reliable photometry and intermediate mid-IR colors. Morphologically, we find that the canyon harbors a large population of both Sa-Sbc and E/S0 type galaxies, and that they fall on the optical red sequence rather than the green valley. Finally, we provide a catalog of WISE photometry for 567 of 652 galaxies selected from the sample of 163 compact groups.

Compact groups of galaxies provide conditions similar to those experienced by galaxies in the earlier universe. Recent work on compact groups has led to the discovery of a dearth of mid-infrared transition galaxies (MIRTGs) in IRAC (3.6 - 8.0 micron) color space (Johnson et al. 2007; Walker et al. 2012) as well as at intermediate specific star formation rates (Tzanavaris et al. 2010). However, we find that in compact groups these mid-infrared (mid-IR) transition galaxies in the mid-infrared dearth have already transitioned to the optical ([g-r]) red sequence. We investigate the optical color-magnitude diagram (CMD) of 99 compact groups containing 348 galaxies and compare the optical CMD with mid-IR color space for compact group galaxies. Utilizing redshifts available from SDSS, we identified new galaxy members for 6 groups. By combining optical and mid-IR data, we obtain information on both the dust and the stellar populations in compact group galaxies. We also compare with more isolated galaxies and galaxies in the Coma cluster, which reveals that, similar to clusters, compact groups are dominated by optically red galaxies. While we find that compact group transition galaxies lie on the optical red sequence, LVL+SINGS mid-IR transition galaxies span the range of optical colors. The dearth of mid-IR transition galaxies in compact groups may be due to a lack of moderately star forming low mass galaxies; the relative lack of these galaxies could be due to their relatively small gravitational potential wells. This makes them more susceptible to this dynamic environment, thus causing them to more easily lose gas or be accreted by larger members.