We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present high-sensitivity ($\sigma_I\sim0.2-0.5$ mJy, $\sigma_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $\mu$m dust emission from the young multiple system VLA 1623 in $\rho$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W, and L1527 strongly in the polarized emission, at the $\sim1-3\%$ level. We spatially resolve VLA 1623A into two sources, VLA 1623Aa and VLA 1623Ab, separated by $\sim30$ au and located within a cavity of radius $\sim50$ au within the circumbinary Keplerian disk, as well as the edge-on disk of VLA 1623W. The polarization angle of the emission is uniform across each protostellar source and nearly coincides with each disk's minor axis. The offsets between the minor axis and the polarization angle are not uniformly distributed at the $P\lesssim2\times10^{-4}$ level. The circumbinary disk surrounding VLA 1623Aab is azimuthally symmetrically polarized. Each compact source's emission is partially optically thick ($\tau\gtrsim1$) at 872 $\mu$m, complicating interpretations of polarization involving aligned grains. We find evidence against alignment by radiative flux in each source, particularly in the edge-on VLA 1623W and L1527. We detect astrometric offsets between the polarized emission and the total intensity in VLA 1623Aa, VLA 1623Ab, and VLA 1623B, as predicted if self-scattering in the optically thick limit operates. We conclude that self-scattering is likely responsible for disk-scale polarization at 872 $\mu$m in these systems.

We analyze the dust morphology of 29 transition disks (TDs) observed with ALMA at (sub-) millimeter-emission. We perform the analysis in the visibility plane to characterize the total flux, cavity size, and shape of the ring-like structure. First, we found that the $M_{\rm{dust}}-M_\star$ relation is much flatter for TDs than the observed trends from samples of class II sources in different star forming regions. This relation demonstrates that cavities open in high (dust) mass disks, independent of the stellar mass. The flatness of this relation contradicts the idea that TDs are a more evolved set of disks. Two potential reasons (not mutually exclusive) may explain this flat relation: the emission is optically thick or/and millimeter-sized particles are trapped in a pressure bump. Second, we discuss our results of the cavity size and ring width in the context of different physical processes for cavity formation. Photoevaporation is an unlikely leading mechanism for the origin of the cavity of any of the targets in the sample. Embedded giant planets or dead zones remain as potential explanations. Although both models predict correlations between the cavity size and the ring shape for different stellar and disk properties, we demonstrate that with the current resolution of the observations, it is difficult to obtain these correlations. Future observations with higher angular resolution observations of TDs with ALMA will help to discern between different potential origins of cavities in TDs.

We present 870 $\mu$m ALMA dust polarization observations of 10 young Class 0/I protostars in the Perseus Molecular Cloud. At $\sim$ 0.35$"$ (80 au) resolution, all of our sources show some degree of polarization, with most (9/10) showing significantly extended emission in the polarized continuum. Each source has incredibly intricate polarization signatures. In particular, all three disk-candidates have polarization vectors roughly along the minor axis, which is indicative of polarization produced by dust scattering. On $\sim$ 100 au scales, the polarization is at a relatively low level ($\lesssim 1\%$) and is quite ordered. In sources with significant envelope emission, the envelope is typically polarized at a much higher ($\gtrsim 5\%$) level and has a far more disordered morphology. We compute the cumulative probability distributions for both the small (disk-scale) and large (envelope-scale) polarization percentage. We find that the two are intrinsically different, even after accounting for the different detection thresholds in the high/low surface brightness regions. We perform Kolmogorov-Smirnov and Anderson-Darling tests on the distributions of angle offsets of the polarization from the outflow axis. We find disk-candidate sources are different from the non-disk-candidate sources. We conclude that the polarization on the 100 au scale is consistent with the signature of dust scattering for disk-candidates and that the polarization on the envelope-scale in all sources may come from another mechanism, most likely magnetically aligned grains.

We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

Dust polarization in millimeter (and centimeter) has been mapped in disks around an increasing number of young stellar objects. It is usually thought to come from emission by magnetically aligned (non-spherical) grains, but can also be produced by dust scattering. We present a semi-analytic theory of disk polarization that includes both the direction emission and scattering, with an emphasis on their relative importance and how they are affected by the disk inclination. For face-on disks, both emission and scattering tend to produce polarization in the radial direction, making them difficult to distinguish, although the scattering-induced polarization can switch to the azimuthal direction if the incident radiation is beamed strongly enough in the radial direction in the disk plane. Disk inclination affects the polarizations from emission and scattering differently, especially on the major axis where, in the edge-on limit, the former vanishes while the latter reaches a polarization fraction as large as 1/3. The polarizations from the two competing mechanisms tend to cancel each other on the major axis, producing two low polarization "holes" (one on each side of the center) under certain conditions. We find tantalizing evidence for at least one such "hole" in NGC1333 IRAS4A1, whose polarization observed at 8 mm on the 100 AU scale is indicative of a pattern dominated by scattering close to the center and by direction emission in the outer region. If true, it would imply not only that a magnetic field exists on the disk scale, but that it is strong enough to align large, possibly mm-sized, grains.

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at {\lambda} = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with {\beta} ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

We present multiple epochs of near-IR spectroscopy for a sample of 25 young stars, including T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra of the BrGamma transition of hydrogen, with a resolution of ~3500. Epochs were taken over a span of >1 year, sampling time-spacings of roughly one day, one month, and one year. The majority of our targets show BrGamma emission, and in some cases these are the first published detections. Time-variability is seen in approximately half of the targets showing BrGamma emission. We compare the observed variability with expectations for rotationally-modulated accretion onto the central stars and time-variable continuum emission or extinction from matter in the inner disk. Our observations are not entirely consistent with models of rotationally-modulated magnetospheric accretion. Further monitoring, over a larger number of epochs, will facilitate more quantitative constraints on variability timescales and amplitudes, and a more conclusive comparison with theoretical models.

We present near-IR spectra of a sample of T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra with a resolution of ~3500. Here we present spectra of the v=2->0 and v=3->1 bandheads of ro-vibrational transitions of carbon monoxide. We observed these spectra over multiple epochs spaced by a few days and approximately one month. Several of our targets show CO emission or absorption features. However we see little evidence of variability in these features across multiple epochs. We compare our results with previous observations, and discuss the physical implications of non-variable CO emission across the sampled timescales.

Oscillations were observed across the whole solar disk using the Doppler shift and line depth of spectral lines from the CO molecule near 4666~nm with the National Solar Observatory's McMath/Pierce solar telescope. Power, coherence, and phase spectra were examined, and diagnostic diagrams reveal power ridges at the solar global mode frequencies to show that these oscillations are solar p-modes. The phase was used to determine the height of formation of the CO lines by comparison with the IR continuum intensity phase shifts as measured in Kopp et al., 1992; we find the CO line formation height varies from 425 < z < 560 km as we move from disk center towards the solar limb 1.0 > mu > 0.5. The velocity power spectra show that while the sum of the background and p-mode power increases with height in the solar atmosphere as seen in previous work, the power in the p-modes only (background subtracted) decreases with height, consistent with evanescent waves. The CO line depth weakens in regions of stronger magnetic fields, as does the p-mode oscillation power. Across most of the solar surface the phase shift is larger than the expected value of 90 degrees for an adiabatic atmosphere. We fit the phase spectra at different disk positions with a simple atmospheric model to determine that the acoustic cutoff frequency is about 4.5 mHz with only small variations, but that the thermal relaxation frequency drops significantly from 2.7 to 0 mHz at these heights in the solar atmosphere.