We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

We present the highest resolution (1") $^{12}$CO(1-0) observations of molecular gas in the dwarf starburst galaxy NGC 625, to date, obtained with ALMA. Molecular gas is distributed in discrete clouds within an area of $0.4$ kpc$^2$ and does not have well-ordered large-scale motions. We measure a total molecular mass in NGC 625 of $5.3\times 10^6 M_\odot$, assuming a Milky Way CO-to-H$_2$ conversion factor. We use the CPROPS package to identify molecular clouds and measure their properties. The 19 resolved CO clouds have a median radius of 20 pc, a median line width 2.5 km s$^{-1}$, and a median surface density of 169 $M_\odot$pc$^{-2}$. Larson scaling relations suggest that molecular clouds in NGC 625 are mostly in virial equilibrium. Comparison of our high-resolution CO observations with a star formation rate map, inferred from ancillary optical observations observations, suggests that about 40% of the molecular clouds coincide with the brightest HII regions. These bright HII regions have a range of molecular gas depletion timescales, all within a factor of $\sim 3$ of the global depletion time in NGC 625 of 106-134 Myr. The highest surface density molecular clouds towards the southwest of the galaxy, in a region we call the Butterfly, do not show strong star formation activity and suggest a depletion time scale longer than 5 Gyr.

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from ground-based data at the distance of 2 Mpc. For the two star-forming complexes, in which no O-type stars had previously been identified, we find a total of 13 newly identified O-type stars in the mass range 15-50 M$_{\odot}$, as well as 4 Wolf-Rayet stars. We use the derived massive stellar content to analyze the impact of stellar feedback on the HII regions. As already found for HII regions in the Magellanic Clouds, the dynamics of the analyzed NGC 300 HII regions are dominated by a combination of the pressure of the ionized gas and stellar winds. By comparing the derived ionized gas mass loading factors to the total gas mass loading factor across the NGC 300 disk, we find that the latter is an order of magnitude higher, either indicating very early evolutionary stages for these HII regions, or being a direct result of the multi-phase nature of feedback-driven bubbles. Moreover, we analyze the relation between the star formation rate and the pressure of the ionized gas as derived from small (<100 pc) scales, as both quantities are systematically overestimated when derived on galactic scales. With the wealth of ongoing and upcoming IFU instruments and programs, this study serves as a pathfinder for the systematic investigation of resolved stellar feedback in nearby galaxies, and it delivers the necessary analysis tools to enable massive stellar content and feedback studies sampling an unprecedented range of HII region properties across entire galaxies in the nearby Universe.

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$\alpha$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$\alpha$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.

It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular ring, northern and southern spiral arm segments). While our measurements roughly fall on the well-known correlation between total infrared and HCN luminosity, bridging the gap between Galactic and extragalactic observations, we find systematic offsets from that relation for different dynamical environments probed in M51, e.g. the southern arm segment is more quiescent due to low star formation efficiency (SFE) of the dense gas, despite having a high dense gas fraction. Combining our results with measurements from the literature at 100pc scales, we find that the SFE of the dense gas and the dense gas fraction anti-correlate and correlate, respectively, with the local stellar mass surface density. This is consistent with previous kpc-scale studies. In addition, we find a significant anti-correlation between the SFE and velocity dispersion of the dense gas. Finally, we confirm that a correlation also holds between star formation rate surface density and the dense gas fraction, but it is not stronger than the correlation with dense gas surface density. Our results are hard to reconcile with models relying on a universal gas density threshold for star formation and suggest that turbulence and galactic dynamics play a major role in setting how efficiently dense gas converts into stars.

We present new ${ }^{12}$CO(J=1-0) observations of Henize 2-10, a blue compact dwarf galaxy about 8.7 Mpc away, taken with the Atacama Large Millimeter Array. These are the highest spatial and spectral resolution observations, to date, of the molecular gas in this starburst galaxy. We measure a molecular mass of $1.2\times10^8 M_\odot$ in Henize 2-10, and most of the molecular gas is contained within a region having a size of about 310 pc. We use the CPROPS algorithm to identify 119 resolved giant molecular clouds distributed throughout the galaxy, and the molecular gas contained within these clouds make up between 45 to 70% of the total molecular mass. The molecular clouds in Henize 2-10 have similar median sizes (~26 pc), luminous masses (~$ 4\times 10^5$ $M_\odot$), and surface densities (~$180$ $M_\odot$ pc$^{-2}$) to Milky Way clouds. We provide evidence that Henize 2-10 clouds tend to be in virial equilibrium, with the virial and luminous masses scaling according to $M_{vir}\propto M_{lum}^{1.2\pm0.1}$, similar to clouds in the Milky Way. However, we measure a scaling relationship between luminous mass and size, $M_{vir}\propto R^{3.0\pm0.3}$, that is steeper than what is observed in Milky Way clouds. Assuming Henize 2-10 molecular clouds are virialized, we infer values of the CO-to-H$_2$ conversion factor ranging from 0.5 to 13 times the standard value in the Solar Neighborhood. Given star formation efficiencies as low as 5%, the most massive molecular clouds in Henize 2-10 currently have enough mass to form the next generation of super-star clusters in the galaxy.

We present the first spatially resolved (~0.5 kpc) measurements of the molecular gas depletion time $\tau_{depl}$ across the disk of the interacting spiral galaxy NGC\,2276, a system with an asymmetric morphology in various SFR tracers. To estimate $\tau_{depl}$, we use new NOEMA observations of the $^{12}$CO(1-0) emission tracing the bulk molecular gas reservoir in NGC 2276, and extinction-corrected H$\alpha$ measurements obtained with the PMAS/PPaK integral field unit for robust estimates of the SFR. We find a systematic decrease in $\tau_{depl}$ of 1-1.5 dex across the disk of NGC 2276, with a further, abrupt drop in $\tau_{depl}$ of ~1 dex along the galaxy's western edge. The global $\tau_{depl}$ in NGC 2776 is $\tau_{depl}=0.55$ Gyr, insistent with literature measurements for the nearby galaxy population. Such a large range in $\tau_{depl}$ on sub-kpc scales has never previously been observed within an individual isolated or pre-merger system. When using a metallicity-dependent molecular gas conversion factor X$\rm_{CO}$ the variation decreases by 0.5 dex. We attribute the variation in $\tau_{depl}$ to the influence of galactic-scale tidal forces and ram pressure on NGC 2276's molecular interstellar medium (ISM). Our observations add to the growing body of numerical and observational evidence that galaxy-galaxy interactions significantly modify the molecular gas properties and star-forming activity within galactic disks throughout the interaction, and not just during the final merger phase.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local (<kpc) and global galaxy environment. We present first results from two large observing programs on ALMA and VLT/MUSE, mapping cloud scales (1arcsec = 47pc) in both molecular gas and star forming tracers across 90 kpc^2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and HII regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find the increase of scatter in the star formation relation on small scales is consistent with gas and HII regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and HII regions, which we test by directly matching our catalog of 1502 HII regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar mass surface density, molecular gas boundedness) on 50 pc scales.

We have conducted a $^{12}$CO(2-1) survey of several molecular gas complexes in the vicinity of H II regions within the spiral galaxy NGC 300 using the Atacama Large Millimeter Array. Our observations attain a resolution of 10 pc and 1 km s$^{-1}$, sufficient to fully resolve Giant Molecular Clouds (GMCs), and are the highest to date obtained beyond the Local Group. We use the CPROPS algorithm to identify and characterize 250 GMCs across the observed regions. GMCs in NGC 300 appear qualitatively and quantitatively similar to those in the Milky Way disk: they show an identical scaling relationship between size $R$ and linewidth $\Delta V$ ($\Delta V \propto R^{0.48\pm0.05}$), appear to be mostly in virial equilibrium, and are consistent with having a constant surface density of $60$ M$_{\odot}$ pc$^{-2}$. The GMC mass spectrum is similar to those in the inner disks of spiral galaxies (including the Milky Way). Our results suggest that global galactic properties such as total stellar mass, morphology, and average metallicity may not play a major role in setting GMC properties, at least within the disks of galaxies on the star-forming main sequence. Instead, GMC properties may be more strongly influenced by local environmental factors such as the mid-plane disk pressure. In particular, in the inner disk of NGC 300 we find this pressure to be similar to that in the local Milky Way but markedly lower than that in the disk of M51 where GMCs are characterized by systematically higher surface densities and a higher coefficient for the size-linewidth relation.

We used multi-band observations by the Herschel Space Observatory to study the dust emission properties of the nearby spiral galaxy NGC 300. We compiled a first catalogue of the population of giant dust clouds (GDCs) in NGC 300 and give an estimate of the total dust mass of the galaxy. We carried out source detection with the multiwavelength source extraction algorithm getsources and calculated physical properties of the GDCs, including mass and temperature, from five-band Herschel PACS and SPIRE observations from 100-500 $\mu$m; the final size and mass estimates are based on the observations at 250 $\mu$m that have an effective spatial resolution of $\sim$170 pc. We correlated our final catalogue of GDCs to pre-existing catalogues of HII regions to infer the number of GDCs associated with high-mass star formation and determined the H$\alpha$ emission of the GDCs. Our final catalogue of GDCs includes 146 sources, 90 of which are associated with known HII regions. We find that the dust masses of the GDCs are completely dominated by the cold dust component and range from $\sim$1.1$\cdot$10$^{3}$ - 1.4$\cdot$10$^{4}$ M$_{\odot}$. The GDCs have effective temperatures of $\sim$13-23 K and show a distinct cold dust effective temperature gradient from the centre towards the outer parts of the stellar disk. We find that the population of GDCs in our catalogue constitutes $\sim$16% of the total dust mass of NGC 300, which we estimate to be about 5.4$\cdot$10$^{6}$ M$_{\odot}$. At least about 87% of our GDCs have a high enough average dust mass surface density to provide sufficient shielding to harbour molecular clouds. We compare our results to previous pointed molecular gas observations in NGC 300 and results from other nearby galaxies and also conclude that it is very likely that most of our GDCs are associated with complexes of giant molecular clouds.

A primary goal of many undergraduate- and graduate-level courses in the physical sciences is to prepare students to engage in scientific research, or to prepare students for careers that leverage skillsets similar to those used by research scientists. Even for students who may not intend to pursue a career with these characteristics, exposure to the context of applications in modern research can be a valuable tool for teaching and learning. However, a persistent barrier to student participation in research is familiarity with the technical language, format, and context that academic researchers use to communicate research methods and findings with each other: the literature of the field. Astrobites, an online web resource authored by graduate students, has published brief and accessible summaries of more than 1,300 articles from the astrophysical literature since its founding in 2010. This article presents three methods for introducing students at all levels within the formal higher education setting to approaches and results from modern research. For each method, we provide a sample lesson plan that integrates content and principles from Astrobites, including step-by-step instructions for instructors, suggestions for adapting the lesson to different class levels across the undergraduate and graduate spectrum, sample student handouts, and a grading rubric.

We present the first high angular resolution study of giant molecular clouds (GMCs) in the nearby spiral galaxy NGC 300, based on observations from the Submillimeter Array (SMA). We target eleven 500 pc-sized regions of active star formation within the galaxy in the CO(J=2-1) line at 40 pc spatial and 1 km/s spectral resolution and identify 45 individual GMCs. We characterize the physical properties of these GMCs, and find that they are similar to GMCs in the disks of the Milky Way and other nearby spiral galaxies. For example, the GMC mass spectrum in our sample has a slope of 1.80+/-0.07. Twelve clouds are spatially resolved by our observations, of which ten have virial mass estimates that agree to within a factor of two with mass estimates derived directly from CO integrated intensity, suggesting that the majority of these GMCs are bound. The resolved clouds show consistency with Larson's fundamental relations between size, linewidth, and mass observed in the Milky Way. We find that the linewidth scales with the size as DeltaV ~ R^(0.52+/-0.20), and the median surface density in the subsample is 54 Msun/pc^(-2). We detect 13CO in four GMCs and find a mean 12CO/13CO flux ratio of 6.2. Our interferometric observations recover between 30% and 100% of the integrated intensity from the APEX single dish CO observations of Faesi et al. 2014, suggesting the presence of low-mass GMCs and/or diffuse gas below our sensitivity limit. The fraction of APEX emission recovered increases with the SMA total intensity as well as with the star formation rate.

We present the results of a galaxy-wide study of molecular gas and star formation in a sample of 76 HII regions in the nearby spiral galaxy NGC 300. We have measured the molecular gas at 250 pc scales using pointed CO(J=2-1) observations with the APEX telescope. We detect CO in 42 of our targets, deriving molecular gas masses ranging from our sensitivity limit of ~10^5 Msun to 7x10^5 Msun. We find a clear decline in the CO detection rate with galactocentric distance, which we attribute primarily to the decreasing radial metallicity gradient in NGC 300. We combine GALEX FUV, Spitzer 24 micron, and H-alpha narrowband imaging to measure the star formation activity in our sample. We have developed a new direct modeling approach for computing star formation rates that utilizes these data and population synthesis models to derive the masses and ages of the young stellar clusters associated with each of our HII region targets. We find a characteristic gas depletion time of 230 Myr at 250 pc scales in NGC 300, more similar to the results obtained for Milky Way Giant Molecular Clouds than the longer (>2 Gyr) global depletion times derived for entire galaxies and kpc-sized regions within them. This difference is partially due to the fact that our study accounts for only the gas and stars within the youngest star forming regions. We also note a large scatter in the NGC 300 SFR-molecular gas mass scaling relation that is furthermore consistent with the Milky Way cloud results. This scatter likely represents real differences in giant molecular cloud physical properties such as the dense gas fraction.

Educational technology has attained significant importance as a mechanism for supporting experiential learning of science concepts. However, the growth of this mechanism is limited by the significant time and technical expertise needed to develop such products, particularly in specialized fields of science. We sought to test whether interactive, educational, online software modules can be developed effectively by students as a curriculum component of an advanced science course. We discuss a set of fifteen such modules developed by Harvard University graduate students to demonstrate various concepts related to astronomy and physics. Their successful development of these modules demonstrates that online software tools for education and outreach on specialized topics can be produced while simultaneously fulfilling project-based learning objectives. We describe a set of technologies suitable for module development and present in detail four examples of modules developed by the students. We offer recommendations for incorporating educational software development within a graduate curriculum and conclude by discussing the relevance of this novel approach to new online learning environments like edX

Recent studies have identified several young stellar objects (YSOs) which exhibit significant mid-infrared (mid-IR) variability. A wide range of physical mechanisms may be responsible for these variations, including changes in a YSO's accretion rate or in the extinction or emission from the inner disk. We have obtained and analyzed multi-epoch near-infrared (NIR) spectra for five actively accreting YSOs in the $\rho$ Oph star-forming region along with contemporaneous mid-IR light curves obtained as part of the YSOVAR Spitzer/IRAC survey. Four of the five YSOs exhibit mid-IR light curves with modest ($\sim 0.2$--0.4 mag) but statistically significant variations over our 40-day observation window. Measuring the strengths of prominent photospheric absorption lines and accretion sensitive \ion{H}{1} and \ion{He}{1} lines in each NIR spectrum, we derive estimates of each YSO's spectral type, effective temperature ({\Teff}), and $H$ band extinction ($A_H$), and analyze the time evolution of their NIR veiling ($r_H$ and $r_K$) and mass accretion rates ({\Macc}). Defining a YSO's evolutionary stage such that heavily veiled, high accretion rate objects are less evolved than those with lower levels of veiling and ongoing accretion, we infer that GY 314 is the most evolved YSO in our sample, with GY 308 and GY 292 at progressively earlier evolutionary stages. Leveraging our multi-epoch, multi-wavelength dataset, we detect significant variations in mass accretion rates over timescales of days to weeks, but find that extinction levels in these YSOs remain relatively constant. We find no correlation between these YSO mid-IR light curves and time-resolved veiling or mass accretion rates, such that we are unable to link their mid-IR variability with physical processes localized near the inner edge of the circumstellar disk or within regions which are directly responsive to mass accretion...